
2
Figure 2. The entanglement spectra of anomalous surface
states exhibit a large number of low-lying states (λi≈0.5),
that increases with the cutoff. Panels (a) and (b) show the
spectra for a gapless and massive Dirac cone, respectively,
with cutoff Λ = 10. In both cases, the number of states
within the shaded low-energy window scales linearly with the
cutoff and is thus non-universal (c).
logical order32–36.
Entanglement spectra of anomalous surfaces.
Any surface wave function of a free-fermion SPT may
be expressed as
|ΦiΛ=ˆ
ΦΛ(c†
E,i)|0i.(1)
Here ˆ
ΦΛis an operator-valued function, and c†
E,i creates
an electron in a single-particle eigenstate with energy E,
and additional quantum numbers i. The ‘empty state’
|0idenotes a Fermi sea filled up to the top of the va-
lence band at energy Λ. The symmetric, non-interacting
surface is thus represented by |Φi0
Λ=QiQµ
E=Λ c†
E,i |0i
with chemical potential µ. In general, |ΦiΛmust reduce
to |Φi0
Λfor cE→Λto describe a surface state that does
not hybridize with the bulk. The cutoff Λ is inevitable
due to the anomalous nature of the surface state, but its
precise value does not affect local observables.
To test our expectations for the ES of anomalous
states, we study the surface of a 3D TI, which hosts a
single Dirac cone governed by
HM=X
k
φ†
k[k·σ+Mσz]φk.(2)
Here, φ†
kcreates electrons with momentum k= (kx, ky),
σx, σy, σzare Pauli matrices, and Mis a time-reversal
symmetry-breaking mass. The un-normalized single-
particle eigenstates are
vk=√+M
√−Me−iϕk,uk=√−Meiϕk
−√+M,(3)
where =√k2+M2and k, ϕkare the magnitude and
polar angle of the two-dimensional momentum. The cor-
responding single-particle energies are Ev=−Eu=,
and the ground state of HMis |ΦMiΛ=Qkukφ†
k|0i. In
Appendix A, we show the spectrum of Eq. (2) for spher-
ical geometry.
We obtain the ES by decomposing the Hilbert space of
the surface into two parts, H=Hα⊗ Hβ. The Schmidt
decomposition
|Φi=Xne−λn/2|Φα,ni⊗|Φβ,ni,(4)
with |Φα(β),ni∈Hα(β)yields the ‘entanglement energies’
λi. We numerically compute these numbers for |ΦMiΛ
with α, β, the two hemispheres of a sphere with radius
Rand (Fig. 2). For local observables, we would expect
universal results when R−1MΛ. For the ES, we
instead find that the number of low-lying λigrows lin-
early with the cutoff Λ. In Fig. 2panel (c), we show the
number of pseudo energy states as a function of the cutoff
for massless and massive Dirac fermions. We attribute
this cutoff dependence and the large number of low-lying
states to the exposure of the underlying bulk (cf. Fig. 1).
As anticipated in the introduction, straightforward com-
putation of a surface state’s ES is not suitable for its
identification.
Relative entanglement spectra. We adapt entan-
glement spectroscopy to SPT surfaces by drawing on in-
sights into physical surface spectra (cf. Appendix B).
Recall that edge-energy spectra are only meaningful for
two surface states with the same anomaly. However, their
boundary state can equivalently arise at the physical edge
of a specific non-anomalous phase, which is subject to the
standard ES-edge correspondence. Our strategy is thus
to construct non-anomalous wave functions that encode
the boundary between two surface states Aand B.
To obtain the desired wave functions, we begin with
a gapped free-fermion surface state A, whose particle-
and hole-like excitations are created by cA,†
+,i , cA,†
−,i . By
construction, Acorresponds to a non-anomalous wave
function of particles and holes, i.e., their trivial vacuum.
Next, we consider a puddle of any other surface state
embedded within A. It can be created locally from par-
ticle and hole excitations on top of the uniform Astate,
which does not introduce any anomaly. Likewise, any
surface state Bcan be encoded in a non-anomalous wave
function of the excitations cA
±,i. Such wave functions de-
pend on both Aand B; we thus refer to them as relative
wave functions from which we obtain the relative ES. We
now proceed by numerically calculating the relative ES of
various 3D TI surface states and comparing them with
the physical edge spectra. Subsequently, we elaborate
on the relative wave functions and provide an analytical
perspective on their ES.
Numerical results I: Free fermions. We compute
the relative ES described above for various surface states
of a spherical 3D TI. For |ΦAiΛ, we take states with a
magnetic or superconducting gap, i.e., the ground states
of Eq. (2) or of
HSC
∆s=H0+ ∆sX
k
[φk,↑φ−k,↓+ H.c.] .(5)