arXiv:1705.06421v4 [cond-mat.str-el] 14 Jan 2019 a nihrn ednytwrstefraino spin- of formation surface), the towards TI tendency the inherent liq- (e.g., an Dirac has locking to spin-momentum helical contrast with 2D sharp uid the exhibits in lead Furthermore semiconductors, that conventional insulator can anomaly. excitonic the excitons parity topological the of on novel condensation a study to and that a show formation We devising fermions. the Dirac by of instability states particle-hole topological (2D) surface [16–19]. TI field the the Zeeman in and a effect with 15], (QAH) Hall [14, anomalous effect quantum proximity (SC) superconductor eettpso astrsaecniee,icuigthe including dif- considered, when are terms phenomena p mass exotic of parity of types the number ferent exhibit a to (TI) and predicted insulator anomaly been Indeed, topological have of states effects. fermions surface these Dirac test gapless to the provide platform they quintessential and have [4–13], a physics attention matter broad drawn condensed recently in Quasiparticles anomaly. dispersion parity linear called gauge with is external effect the This or of field. Abelian an theory in Chern-Simons through resulting nonAbelian [1–3], symmetries correction these radiative break the inevitability will mass infinitesimal an However, (TRS). symmetry time-reversal aietesraeeetosaddietesystem the proposed drive the and towards electrons surface renor- the can topologi- malize pumping optical the interband spin-momentum the realize the locking, With to insulator. platform excitonic cal promising a state surface is insulator topological mag- The in-plane spontaneous netization. a and defects at charge eato mn h lcrn n h oe in holes the to and leads electrons usually in- semiconductors the attractive Isotropic among the since teraction considered condensates states. commonly excitonic are BCS-like for functions into wave pairs electron spatial of hole condensation coherent and the by formed ymtyadi hrceie yacia Chern chiral a by number characterized is and symmetry ee epooeanwtp fectncinsu- with excitonic parameter of order type exhibits new that a lator propose we Here, a rpris rcinlzdecttoswt e with topologi- excitations novel fractionalized two properties: in cal results which anomaly, ity + nti ok eepn h el ftwo-dimensional of realm the expand we work, this In asesDrcprilspeev ohprt and parity both preserve particles Dirac Massless xioi nuaosaeisltn states insulating are insulators Excitonic ip yetplgclsprodco nue ythe by induced superconductor topological type rdcino Topological a of Prediction C c 2 eateto hsc n srnm,Saga ioTn Uni Tong Jiao Shanghai Astronomy, and Physics of Department 1 = / 2 hssaedsly h par- the displays state This . 3 eateto hsc,AioaSaeUiest,Tme AZ Tempe, University, State Arizona Physics, of Department n eateto hsc,NnigUiest,Nnig210 Nanjing University, Nanjing Physics, of Department and u Wang Rui p + olbrtv noainCne fAvne Microstructu Advanced of Center Innovation Collaborative ip instability. 1 ainlLbrtr fSldSaeMicrostructures, State Solid of Laboratory National 1 , 2 s nrErten Onur , wv excitons. -wave p + p ip 3 + , xioi nuao ihPrt Anomaly Parity with Insulator Excitonic ∗ / ip 2 agn Wang Baigeng ltsteTS ocaatrz uheoi tts we number states, Chern exotic vi- chiral such only characterize a but To propose symmetry TRS. gauge the U(1) olates global the break not iaeteectncisaiiyo togycreae TI strongly-correlated a inves- of we instability Then, excitonic the state. tigate ground of magnetization in-plane hsi ute epnil o i rcinlzdexcita- fractionalized e (i) charge for with responsible tions further defects. is topological Chern- with This effective associated an theory gauge deriving Simons by anomaly parity derlying the unlike However, basis. chiral the ufc,a ilb hw eo,cnb sdfrgenera- for unique used a TI be of of can stabilization below, property and shown optical a tion be therefore unique will and as This surface, polarization, spin state. same spin-triplet the and with electron an hole of the a consists of exciton spin-flip to resultant in- the the Due direct by electron; accompanied the is Fig.1(b). surface, however, excitation TI by is, terband of shown situation locking as The spin-momentum surfaces, the pair. TI a in into different exciton form spin-singlet they a in after resulting (Fig.1(a)), spins posite otayt h ovninlectncisltr 2]de- In [20] a by insulators model. scribed excitonic conventional theoretical the minimal to a contrary studying by state tonic aec oewt h aemomentum same the with direct hole the and valence electron coupling, a conduction spin-orbit a generates the excitation of interband semiconduc- absence conventional For with Fig.1. tors pro- by absorption shown optical as interband it- cesses, the manifests by clearly property most remarkable self This excitons. triplet has efisl xlr h ain etr fthe of feature salient the explore firstly We ahmtclyrsmlstespinless the resembles mathematically ao.Tentrltopological natural The lator. etoa eiodcos b h eia pntxueo T of texture spin favors helical state The surface (b) semiconductors. ventional I.1 ceai lto h ietitradotclab- optical interband direct the (a) of plot process. Schematic sorption 1: FIG. 澻濴澼 p + ip 1 ymtyihrn oterodrparameters. order their to inherent symmetry , † s n .Y Xing Y. D. and wv C-ieter 2] the [21], theory BCS-like -wave est,Saga 020 China 200240, Shanghai versity, p / wv ahrthan rather -wave s tvrie n i)aspontaneous a (ii) and vortices at 2 wv xiosaefre ntecon- the in formed are excitons -wave 9,China. 093, 55 USA. 85257 res, 澻濵澼 1 p wv xioi insulator excitonic -wave C p c + s n nelteun- the unveil and wv excitons. -wave p ip p + + xioi insu- excitonic ip k ip n h op- the and p C tdoes it SC, p C[2 in [22] SC + + ip ip state exci- I 2 surface model where we show that a p-wave type inter- (a) (b) d action component emerges intrinsically, which stabilizes d y y the p + ip excitonic phase in a significant parameter re- gion of the phase diagram. Moreover, we also consider the effects of interband optical processes and we demon- d d strate that the electron-photon coupling can effectively z z drive the system into the proposed p + ip excitonic insu- lator as a result of the spin-flip mechanism demonstrated in Fig.1(b). d x d Results x Model. We start from a spontaneous symmetry FIG. 2: The mapping of dˆ(k). (a) The mapping in the breaking of 2D Dirac fermions with spin-momentum original spin basis. (b) The mapping in the chiral basis. The locking and a momentum-dependent mass term, H = difference between the two situations shows that the two topo- HDirac + Mˆ (k), where logical invariants are needed to completely characterize the p + ip excitonic insulator, i.e., (C,Cc)=(0, 1/2). H = Φ† v σ kΦk. (1) Dirac k F · k X obtained after a Bogoliubov transformation as, T k is the 2D wave vector, Φk = [ψk, , ψk, ] is the † ↑ ↓ Pk gkck ck,− Ψ = uke ,+ 0 , (3) fermionic spinor and ψk†,σ (ψk,σ) denotes the cre- | i | i k ation (annihilation) operator of electrons with spin σ. Y ˆ We firstly focus on a concrete mass term M(k) = with gk = vk/uk, uk = ∆c(kx iky)/kξk, vk = vFk/ξk, z 2 x y 2 − k Φk† (∆c cos θσ + ∆c sin θσ ∆c sin θ cos θσ )Φk, where ξ = ∆2 + v k2 is the energy spectrum of − k c F while under which conditions Mˆ (k) can be spontaneously quasi-particles. Ψ clearly describes a coherent state of P p| i generated will be investigated below. θ is the polar angle electron-hole pairs with the distribution factor gk of the of k, and ∆c is the overall gap scale. It is straightforward p + ip symmetry. 1 to verify that Mˆ (k) breaks TRS, TˆMˆ (k)Tˆ− = Mˆ ( k), Topological classification and chiral Chern num- 6 − where Tˆ = iσyK being the time-reversal operator, with ber. To examine the topological classification of the K the complex conjugation operator. In contrast with p + ip excitonic insulator, we write Eq.(2) into the com- the TI surface with a Zeeman field, where the diago- 3 i i pact form as, = i=1 d(k)τ , with τ being the Pauli nal term mσz breaks TRS [23], here the diagonal term matrices definedH in the chiral basis. For the continuum z ˆ P ∆c cos θσ in M(k) respects the TRS. The TRS is how- model, where k is envisioned as a single point, the ˆ → ∞ ever violated by the two off-diagonal terms in M(k). three-vector d(k) = ( ∆c cos θ, ∆c sin θ, vFk) defines a One can unveil the physical meaning of Mˆ (k) by mak- mapping from the k-space− (S2)− to the d-vector space. ing a unitary transformation onto the chiral basis with Since the z-component vFk 0, the normalized d-vector ˆ ˆ ≥ HDirac being diagonalized: Ck = R(k)Φk, where R(k) is dˆ(k) = d(k)/ d(k) resides in a hemisphere surface. For ˆ ˆ iθ | | a 2 by 2 rotation matrix with R11 = R21 = e /√2, dˆ(k) covering the hemisphere n times, the winding num- − T Rˆ12 = Rˆ22 =1/√2, and Ck = [ck,+,ck, ] is the spinor − ber becomes n/2. This winding number is a topological T composed of the conduction and valence band fermions. invariant defined in the chiral basis Ck = [ck,+,ck, ] , − In the chiral basis, H becomes which we term as the chiral Chern number Cc. Recalling that for the TI surface with a Zeeman field, a iθ ck,+ † vFk ∆ce− ck,+ conventional Chern number C is defined in the same way H = iθ − . ck, ∆ce vFk ck, T k − − [16, 23] but in the spin basis, Φk = [ψk, , ψk, ] . The X    − −    ↑ ↓ (2) difference between C and Cc becomes clear from the per- iθ The off-diagonal ∆ce− describes the process where an spective of Berry phase as following. In the original spin electron in the valence band is scattered up to the con- basis, following a closed loop in momentum space, the duction band, leaving a particle-hole excitation on top of accumulation of Berry phase is φ = i dk uk k uk , − · h |∇ | i the many-body ground state. For ∆c = 0, the particle- where uk is the eigenstate of the valence band in hole pairs condense and gap out the Dirac6 cone, gener- the spin| basis.i For a Bloch state definedH in the chi- ating an excitonic insulator. Another interesting feature ral basis, a similar geometric phase can be defined as iθ c c c c of the off-diagonal term is the factor e± = (kx iky)/k φ = i dk uk k uk ,where uk is the eigenstate [24], similar to the spinless p + ip-wave SC [25],± but we of the− valence· h band|∇ | ini the chiral| basis.i The Bloch note that the order parameter resides in the particle-hole states inH the above two formalisms are related to each c ˆ rather than the particle-particle basis. The model Hamil- other by the unitary transformation uk = R(k) uk . tonian can be analyzed similarly as the BCS supercon- It then becomes obvious that the two| Berryi phases| sat-i c ductor theory. The ground state wave function Ψ is isfy φ φ = i dk uk Rˆ†(k) kRˆ(k) uk , leading to | i − − · h | ∇ | i H 3

C C = (i/2π) dk uk Rˆ†(k) kRˆ(k) uk =1/2, after . A unique normalized zero energy bound state [32, 33] c − − ·h | ∇ | i inserting the Rˆ(k) matrix we used to diagonalize H . can be obtained as, u = Af(r)ei(α+π/2)/2 and v = u⋆, H Dirac To explicitly demonstrate the topological invariants, where A is the normalization constant and the function we plot the normalized vector dˆ(k) with traversing the f(r) = exp[ ∆c r/ vF ] decays exponentially. −| | | | whole k-space. Fig.2(a) and (b) show the mapping of We note that the single particle Hamiltonian h in 2 2 the dˆ(k) vector in the spin and chiral basis, respectively. Eq.(2) satisfies (iτ )h(iτ ) = h. This allows one to de- 2 It is found that dˆ(k) in latter basis completely covers the fine a time-reversal-type operator Θ = iτ K, with which, it is clear that for any eigenstate ψ with energy E, hemisphere while dˆ(k) in the former does not, leading to E C = 1/2 and C = 0. Therefore, two topological invari- ΘψE is also a eigenstate with energy E. Besides, the c total number of fermion states are conserved,− namely, ants are needed to completely characterize the p + ip ex- ∞ dEδρ(E) = 0, where δρ(E) is the change of den- citonic state, namely, (C, C )=(0, 1/2). The zero Chern c sity−∞ of states. Applying this formula to the case where number C = 0 means that the p+ip excitonic insulator is Ra topological defect is in presence, and then taking into topologically distinct from the TI surface with a Zeeman field that has C = 1/2 [16]. This is clear from the Lan- account the fact that ρ(E) is a symmetric function with respect to E = 0, one comes to the conclusion that both dau quantization of the two states. For the TI surface the conduction and valence band lose half a state. This with a Zeeman field, the Landau level (LL) for n 1 reads as E = 2ev2 B n + m2, where m is the| | Zee- ≥ state is compensated by the zero mode that is bound to n F the topological defect. Consequently, at half filling, the man gap. In addition,± there| | is an unconventional zeroth p electron charge associated to the defect is either e/2 or LL E = sgn(m)m which is located exactly at the gap 0 e/2, depending on whether the zero mode is occupied− or edge [26–28]. This zeroth LL is a manifestation of the empty. nontrivial Zeeman mass as a result of C = 1/2 [27, 28]. For the p + ip excitonic insulator phase, however, one We have shown that the p + ip excitonic insulator is cannot find the unconventional zeroth LL located at the characterized by Cc =1/2 and displays fractional charge gap edge, consistent with C = 0. e/2 at vortices. However, it is still unclear how the two observations are related to each other. This can be an- Chern-Simons theory and fractionalization. The swered by the parity anomaly [34] of the p + ip excitonic nonzero chiral Chern number C = 1/2 indicates that c insulator and the effective gauge theory associated with the p + ip excitonic insulator must be topologically dis- the topological defects. First, we couple an external U(1) tinct from normal insulators. Recalling the p + ip SCs gauge field A (with µ = 0, 1, 2) to the p + ip excitonic where Majorana fermions are encoded in vortices [22], µ state [35], and we introduce the Dirac matrices defined we shift our attention to topological defects in the p + ip as γ0 = τ 3, γ1,2 = iτ 1,2. In order to extract the effec- excitonic state. To this end, we can solve the Bogoliubov tive theory for the vortices, we include the vortex de- de Gennes (BdG)-type equations looking for zero energy grees of freedom χ in the order parameter. Second, we solutions. In real space, the equations can be derived as perform a singular gauge transformation, c e iχ/2c [See methods for details]: ± [36], which results in the following Lagrangian± → [See meth-± ods for details], iα 1 1 e ( i∂r + r− ∂α)u(r) ∆c(r)v(r) =0, (4) − − vF µ µ = k(iγ ∂µ + γ aµ + m) k, (6) iα 1 1 ⋆ L C C e− ( i∂r r− ∂α)v(r)+ ∆c (r)u(r) =0, (5) − − vF µ where k ( k) is thee Grassmanne field. γ are defined as γ0 =C γ0C, γ1 = γ2, γ2 = γ1. In Eq.(6), the vortex where α is the polar angle in real space and u, v are − the two components of the spinor wave function. Similar degrees of freedom have been absorbed intoe an effective equations have been studied in Ref.[22, 29, 30]. How- gaugee field aeµ after thee singular transformation. For a static defect with ∂ χ = 0, a = A + 1 ∂ χ. Finally, ever, note that the physics in the present model is differ- 0 µ µ 2 µ ent from previous works. Ref.[22] derives the equations after one integrates out the Grassmann field, a Chern- for SCs, where charge conservation is broken and Ma- Simons term with respect to the effective gauge field aµ jorana fermions are the emergent excitations. For our can be obtained in the second order expansion: system, the U(1) gauge symmetry is preserved, therefore I µνρ no Majorana modes are allowed. Ref.[29] investigates CS = ǫ aµ∂ν aρ, (7) L 2π graphene with a chirality-mixing texture which induces the mass term, Ref.[30] discusses the exciton condensa- where the coefficient I is found to be 1/2 for the proposed tion between two opposite TI surfaces. In contrary, here p + ip excitonic insulator. we focus on the excitonic order in a Dirac liquid state We note that for the TI surface with a Zeeman field and the mass term is generated by the condensation of [23], a similar Chern-Simons Lagrangian arises as a re- p + ip excitons. To solve the equations in the presence sult of the parity anomaly, where the coefficient in front iα of a topological defect [31], we assume ∆c(r) = ∆c(r)e of the Chern-Simons term is exactly the Chern number 4

C. Here, the coefficient I = 1/2 is identified to be the chiral Chern number Cc defined above. The similarity between the two gauge field theories is due to the un- derlying parity anomaly, however, for the p + ip exci- tonic insulator, it is the effective gauge field aµ rather than the external field Aµ that obeys the Chern-Simons 1 form. Further inserting aµ = Aµ + 2 ∂µχ, the full La- grangian = I ǫµνρA ∂ A A jµ + consists of LCS 2π µ ν ρ − µ Lv a Lagrangian describing the vortex v and a current op- erator associated with the vortex, L I jµ = ǫµνρ∂ ∂ χ. (8) 2π ν ρ For a static vortex (with vortex number n = 1) that cen- ij tered at r0, we have ǫ ∂i∂j χ = 2πδ(r r0), inserting which into Eq.(8), the electron charge− of the defect is then obtained as Q = d2xj0 = I C = 1/2. There- ≡ c fore, we have proved that the chiral Chern number Cc equals to the fractionalR charge (in unit of e) of a topolog- ical defect, justifying Cc to be the essential topological invariant characterizing the p + ip excitonic insulator. Spontaneous in-plane magnetization. The order pa- iθ FIG. 3: The in-plane spin texture of the p + ip excitonic rameter ∆ e− in Eq.(2) breaks TRS. This suggests − c insulator. (a),(b),(c),(d) are the results for different global that the ground state carries a magnetization associated phase of the exciton order parameter, β = 0, β = π/2,β = with it. Before proceeding, we note that the order param- π,β = 3π/2, respectively. The total net magnetization M is iθ eter ∆ce− still enjoys a redundant global phase. In shown by the arrow for each cases. the following,− we take this into account and write the off- iθ′ diagonal term in Eq.(2) as ∆ce− , where θ′ = θ β with β being a constant angle− independent of k. The− 2π when traversing a closed path around k = 0 point. global phase β is analogous to U(1) phase of a supercon- Furthermore, we calculate the total in-plane magnetiza- ductor. tion M by integrating ( σx , σy ) over the momentum h i h i To clearly show the feature of magnetization, we calcu- space. It is found that M always points towards the op- late the spin expectation values σx , σy , σz from the posite direction of the global phase β, as shown by the single particle eigenstates and studyh i h thei spinh i texture of thick red arrows in Fig.3. We also note that the emergent the p+ip excitonic ground state Ψ (Eq.(3)). We obtain magnetization of ground state can lead to anisotropic in- | i plane transport of excitons which serves as an evident

2 experimental signature of the p + ip excitonic insulator. x (∆c cos θ′ + ξk)(vFk cos θ′ + ∆c sin θ′) Furthermore, since β determines the direction of mag- σ = 2 2 2 , h i − vFk + ∆c + ∆c cos θ′ξk netization of ground state, one expects to control β by y (∆c cos θ′ + ξk)(vFk sin θ′ ∆c sin θ′ cos θ′) applying a weak in-plane magnetic field during the pro- σ = 2 2 2 − , h i − vFk + ∆c + ∆c cos θ′ξk cess of TRS symmetry breaking. 2 2 2 2 2 The excitonic phases on interacting TI surfaces. z ∆c sin θ′ ∆c cos θ′ ∆c 2∆c cos θ′ξk σ = − − − . After careful characterization of the topological proper- h i − v2 k2 + ∆2 + ∆ cos θ ξ F c c ′ k ties of the p + ip excitonic insulator, we now study its As has been discussed before, σz does not break TRS. possible realization on interacting TI surface states. Al- The integral of σz over the wholeh ik-space vanishes. The though various instabilities of interacting topological in- in-plane spin textureh i ( σx , σy ) is however nonzero, and sulators and topological superconductors have been ex- we plot it in Fig.3 forh differenti h i values of β, i.e., β = 0, tensively studied [37–48], we focus on the excitonic in- β = π/2, β = π, β =3π/2. Interestingly, it is found that stability of the surface states and propose to realize the the in-plane spin configuration is tilted from the stan- minimal model of the topological excitonic insulators via dard helical spin texture of the original TI surface, and spontaneous symmetry breaking that enjoys a fractional the spins at two TRS-related momentum points, k and chiral Chern number Cc = 1/2, which is different from k, are no longer always opposite to each other, violat- earlier proposals based on semiconductor heterostruc- ing− the TRS. Though being tilted, the chiral nature of ture with an integer Thouless-Kohmoto-Nightingale-Nijs electron’s spin is preserved as we can see that the direc- (TKNN) number [49]. tion of the vector ( σx , σy ) still undergoes a change of The effective model Eq.(1) describing the surface state h i h i 5

where we assume a momentum cutoff Λ implicit in the sum of k. n is the band index. The interaction term (1) V ˆ in the second line, H = k k′ [kck† ck, n] eff − 4 , ,n ,n − · ˆ [k′ck† ′, nck′,n], is of p-wave-type, andP it arises from the spin-momentum-locking− of the surface state. After intro- ducing a complex bosonic vector field ∆ and perform- ing the Hubbard-Stratonovich transformation to decou- (1) ple Heff , we arrive at a mean-field theory with the or- 1 ˆ der parameter given by ∆ = 2 V k kck† , ck,+ . When ∆ = 0, the condensation of excitonsh gaps− out thei Dirac point,6 resulting in the excitonic insulatorP phase. Similar to 3He superfluid theory[25], symmetry analysis implies two possible configurations of the order parameter that are the extrema of the free energy: (1) (chiral) p + ip excitonic insulator with ∆ = 1 ∆ (eˆ ieˆ ) and (2) the √2 c x − y polar excitonic insulator with ∆ = ∆peˆ, where eˆ is a unit vector in the xy plane. The polar excitonic phase is topologically trivial with zero chiral Chern number. By self-consistently solving the mean-field gap equa- tion, we obtain the condensation energy and the or- FIG. 4: The results from mean-field calculations. (a) Con- densation energy versus interaction V for p + ip excitonic in- der parameter (∆c, ∆p), which are plot in Fig.4(a) and sulator and polar excitonic insulator. (b) Calculated mean- Fig.4(b), respectively. Fig.4(b) shows the emergence of field excitonic order parameter for p + ip excitonic insulator a threshold Vc, beyond which the system develops the (∆c) and polar excitonic insulator (∆p) versus interaction V . excitonic instabilities with ∆c/p = 0. Fig.4(a) indicates 2 2 2 | | 6 |∆c/p| /vFΛ ≪ 1 is assumed. (c) The phase diagram based that even though both the p + ip excitonic insulator and on unbiased mean-field with considering both the s- and p+ip- polar excitonic insulator are stable compared to the gap- wave interaction channel. The dash-dotted red line denotes less Dirac cone for V > Vc, the p+ip excitonic insulator is the Hubbard interaction case with Vs = Vp. more energetically favored due to the lower condensation energy. Moveover, we also performed addition mean-field is usually modified by a higher order wrapping effect, calculations by treating the p + ip and polar excitonic in- reducing the rotational symmetry of the Dirac cone down stability on equal footing. It is found that the p + ip to a discrete C3 subgroup. Our further studies on the channel always suppresses the polar order down to zero. wrapping effect shows that the exciton gap remains very Therefore, even though the ground state energies of the robust and therefore one can neglect the wrapping effect two states are close, the mass term of the polar exciton in- for the following analysis.We firstly consider the particle- sulator is negligible in the mean-field theory. Hence, the hole symmetric case, while the generalization to finite resultant mean-field Hamiltonian exactly matches with chemical potential µ is very straightforward and will be Eq.(2), and the proposed mass term Mˆ (k) is generated discussed below. by spontaneous TRS breaking as a result of the p-wave First we start with demonstrating that a p-wave-type (1) interaction Heff . Moreover, we note that for particle- interaction that spontaneously generates the p + ip- hole asymmetric case with a finite chemical potential µ, wave exciton insulator naturally arises from the Hubbard the critical Vc increases with increasing µ. Although the model, which reads in k-space as: exciton insulator is always stable for strong enough in- V teraction V > Vc, it is more experimentally favorable to HI = ψk†+q,σψk† ′ q,σψk′,σψk,σ. (9) tune µ close to the Dirac point [50, 51]. 2 − k k′ q σ X, , X In the above calculations, we temporarily focused on (1) To facilitate the study of excitonic instability, we trans- the p-wave interaction Heff without considering the s- form HI onto the chiral basis. Then, considering only wave component. Nevertheless, the results clearly show interband forward scattering channel, the low-energy ef- that the TI surface state indeed has a strong tendency fective vectorial interaction component can be derived as towards the p + ip excitonic instability via spontaneous (see Methods for details), TRS breaking. We now take into account the interaction in the s-wave channel, which can also be derived from the H = nvFkck† ,nck,n local Hubbard interaction as (see Methods for details), k,n (2) V ′ X (10) Heff = 4 k,k′,n ck† ,nck, nck† ′, nck ,n. In order to be V − − − ˆ ˆ ′ more general, we assume in the following two indepen- [kck† ,nck, n] [k′ck† ′, nck ,n], − 4 − · − P k k′ dent couplings Vp and Vs for the p-wave and s-wave inter- X, ,n 6 action respectively, and treat both terms on equal foot- of energy and momentum requires k q in H , | | ≫ | | el-ph ing. An s-wave bosonic scalar field ∆s = Vs k ck† ck+ and therefore ensures the direct interband excitation with (2) h − i k q k. Furthermore, the photon fields are Gaussian is introduced to decouple Heff in addition toP the p+ip ex- − ≃ with the Hamiltonian H = ωq(a† aq +1/2), citonic order parameter ∆c. The coupled self-consistent ph qλ q,λ ,λ one can therefore exactly integrate them out. This pro- equations with respect to ∆s and ∆c are obtained in the P standard approach which leads to the calculated phase cedure leads to the effective action Sint describing the diagram shown by Fig.4(c). It is found that both the renormalization from photons: s-wave and p + ip-wave excitonic insulator can be spon- 2 2 2πe vF σ σ taneously generated from the Dirac liquid phase with Sint = dt ( αβ eλ)( ργ eλ) − ωq · · kk′q strong enough interactions. The former and latter are Z X αβργλX stabilized for large Vs and Vp respectively. We note that 1 ψk†,αψk,βψk† ′,ρψk′,γ. although the p + ip excitonic phase, which is absent in × i∂ ωq t − conventional semiconductors, becomes stable in a large (13) parameter region of the phase diagram, it does not oc- where (σ e )(σ e )= σx σx + σy σy , insert- cur as the ground state for the discussed local Hubbard λ αβ · λ ργ · λ αβ ργ αβ ργ ing which, it is found that Sint consists of an equal-time model, which lies in the parameter line Vs = Vp that only P enters into the s-wave excitonic condensation. renormalization term that produces the effective interac- The stabilization of the p + ip excitonic insulator. tion between electrons, Now we study the stabilization of the p + ip excitonic in- H = V ψ† ψk ψ† ′ ψk′ , (14) sulator on TI surfaces by going beyond the natural Dirac int − eff k,σ ,σ k ,σ ,σ kk′σ semimetals and taking into account external perturba- X 2 2 2 tions. Recalling the interband optical absorption pro- where Veff = 4πe vFnph/ω¯ and nph is the number den- cess shown by Fig.1(b), one expects that additional pho- sity of photons that interact with the electrons, which is tons should prefer the p-wave pairing due to the spin-flip proportional to the light intensity I. We have utilized the mechanism that gives rise to spin-triplet excitons. Since fact that all the wave vectors q satisfy q = qeˆz =ω ¯eˆz/c excitons are usually generated by absorption of photons for photons with the frequencyω ¯. The first two operators in a material, we consider a linearly-polarized perpendic- in Hint, ψk†,σψk,σ, are generated by the optical interband ular incident laser with frequencyω ¯ and light intensity absorption process with spin-flip depicted by Fig.1(b).

I. The corresponding microscopic theory can be con- Similarly, the last two operators, ψk† ′,σ¯ψk′,σ, are origi- structed by quantization of the vector potential A(r), nated from the interband relaxation process with spin- i.e., flip. After a unitary transformation to the chiral basis, Eq.(14) is exactly cast into: 2π iq r A(r)= eλ(aq,λ + a† q,λ)e− · , (11) Veff ωq − ˆ † ˆ † ′ ′ q s Hint = [kck,nck, n] [k′ck , nck ,n]. (15) λ − 2 − · − X k k′ X, ,n where a (a† ) is the annihilation (creation) operator q,λ q,λ Interestingly, H enjoys the p-wave factor kˆ kˆ , exactly of photons, e (λ = 1, 2) are the polarization vectors, int ′ λ the same as the p-wave interaction component· reduced ωq = c q is the light frequency, with c being the light | | (1) velocity and q the wave vector of photons. from the Hubbard model, Heff . The interband optical The coupling of the external field to the surface elec- processes automatically select the p-wave channel due to trons is described within the formalism of the minimal the spin-flip scatterings (Fig.1(b)). More importantly, coupling, which leads to the electron-photon interaction, Hint has a strength Veff tunable by external field param- eters,ω ¯ and I. According to the mean-field study above and the calculated phase diagram in Fig.4, it is known 2π σ Hel-ph = evF ( σσ′ eλ)(ψk†,σaq,λψk q,σ′ that the tunable external field can renormalize the Dirac − ωq · − kq ′ s Xλσσ cone state and drive it into the proposed p + ip-wave ′ excitonic insulator. The resultant state displays the pre- + ψk†,σa† q,λψk q,σ ), − − (12) dicted fractional charges associated with vortex excita- tions, the spontaneous in-plane magnetization and thus where the term (σσσ′ eλ) arises due to the helical spin the anisotropic transport of excitons on the TI surface. texture and the spin-momentum· locking of TI surface Discussion states. For linearly-polarized perpendicular incident pho- Excitonic insulators require strong interactions since tons, e1 = eˆx and e2 = eˆy. Thus, the spin-flip matrices the topological surface states become unstable only when x y σ and σ take place in Hel-ph. This verifies the spin- the interaction is as large as the ultraviolet cutoff (band flip mechanism illustrated by Fig.1(b) in a microscopic gap) as shown in Fig.4(b). Therefore we argue that topo- perspective. Moreover, since c v , the conservation logical Kondo insulators (TKIs) are good candidates for ≫ F 7 excitonic instability on TI surface. For instance, in the matter system, and may stimulate different directions to case of SmB6, which is the archetypal TKI, the renormal- search for more exotic topological state of matters. We ized Coulomb interactions (U 1eV) is much larger than also note that a laser-induced non-equilibrium version of the direct band gap (Λv 20meV).∼ For these parame- our proposal has been experimentally observed recently F ∼ ters, our mean field theory predicts an excitonic gap of in Bi2Se3 [63], where a millimetre-long transport distance order ∆ 2 10meV. Indeed, magnetism and hystere- of excitons up to 40 K is reported, supporting an exci- s/c ≃ − sis in magnetotransport have been observed in SmB6[52]. tonic condensation state formed on the TI surface. Our Whether these experiments indeed observe an excitonic study on the p + ip excitonic insulator stabilized by the insulator require future work. We would like to note optical pumping further predicts an anisotropic transport that our theory is not related with the recent works on of excitons, which we believe would be of great interest excitonic contributions to quantum oscillations[53] or the for future experimental investigations. intervalley excitonic instabilities [54] in SmB6. Methods In addition to the short-range interaction, we also con- The zero-energy BdG-type equations. From Eq.(2) sidered a long-range Coulomb interaction on the TI sur- in the main text, the zero-energy eigen-equations can be face and make generalization to finite temperature [63]. written down as:

It is found that the critical interaction Vc for BCS exciton iθ condensation is significantly reduced for long-range in- vFkuk ∆(r)e− vk = 0, iθ ⋆ − teraction. At finite temperature, taking into account the e ∆ (r)uk v kvk = 0. − − F thermal fluctuation, the strict long-range order is absent, Here, we assume a real-space dependence of the order pa- the Kosterliz-Thouless transition temperature TKT can iα be calculated by evaluating the phase stiffness of the ex- rameter ∆(r). With a vortex ,we have ∆(r) = ∆c(r)e , citon order based on the mean-field calculations [55]. Ex- where α is the polar angle of r in real-space. To simplify iθ perimentally, the driven TI surface offers a rare platform the eigen-equations, one can multiply a factor e from iθ for realization of new optical devices and facilitates the the right in the first equation, and a factor e− from the study of possible collective instabilities of Dirac electrons. left in the second equation. Then, a Fourier transforma- Experiments measuring bulk metallic TI materials have tion yields, obtained the lifetime of the excited states to be around 1 a few picoseconds [56–58], which is significantly larger ( i∂x + ∂y)ur ∆(r)vr = 0, − − v than than that in graphene. Moreover, a remarkably F 1 ⋆ long lifetime exceeding 4µs has been observed by exper- ∆ (r)ur + ( i∂x ∂y)vr = 0. vF − − iment for the surface states in bulk insulating Bi2Te2Se [59]. These indicate promising observation of a transient After transformation to the polar coordinate, we arrive or quasi-equilibrium excitonic insulator [58, 60]. Besides, at the zero-energy BdG-like equations in real-space. the surface has a reduced dielectric constant which can be tuned even smaller by doping or gating [61], this possibly iα 1 1 e ( i∂r + r− ∂α)u(r) ∆c(r)v(r) = 0, makes the long-range interaction more effective. − − vF To summarize, we proposed a new type of topological iα 1 1 ⋆ e− ( i∂r r− ∂α)v(r)+ ∆c (r)u(r) = 0. p+ip excitonic insulator that achieves the parity anomaly − − vF and is characterized by the nontrivial chiral Chern num- Effective gauge field theory. In the second quantized ber. We provided a concrete example and discussed form, the effective Hamiltonian of the p + ip excitonic its possible realization in strongly-correlated TI surface. insulator reads as Interestingly enough, similar to the TIs due to band z iθ + ⋆ iθ inversion, the predicted topological behaviors, i.e., the H = C† [v kτ ∆ e− τ ∆ e τ −]Ck, k F − c − c charge fractionalization and the emergent magnetization k of ground state, are robust and immune to continuous X Formally, one can make transformation Ck k = deformation of the Hamiltonian, as long as the exciton → C gap remains unclosed. This suggests that realizations of Ck/√k, while keeping the total Hamiltonian unchanged. topological p-wave excitonic insulators are likely in many The action in the functional representation is then cast other physical systems. For example, the spectrum under into: study can be continuously deformed from the Dirac cone z 2 ⋆ ⋆ + S = dτdk k[i∂ + τ v k + ∆ k′γ− + ∆ k′ τ ] k, to a quadratic band crossing point (QBCP) [62]. This C 0 F c c C indicates that the instability is also possible and may be Z more feasible in interacting 2D electron gas with spin- where we have defined k′ = kx iky and introduced the momentum-locking and QBCP. In this sense, the present Dirac matrices in 2D, γ0 = τ 3,−γ1 = iτ 1, γ2 = iτ 2, and 2 1 results are general and may provide new physical settings γ± = (γ iγ )/2. The above action formally describes ± 2 for the realization of quantum anomalies in condensed a 2D Dirac fermion with a mass vFk . Since the sign of 8

(1) (2) mass determines the topology of 2D Dirac systems, and Both Heff and Heff , p-wave and s-wave type respectively, v k2 0 for all k, one can continuously deform the mass are reduced from the original Hubbard interaction, which F ≥ into a constant m while keeping the topological prop- are the least irrelevant channels that can possibly gap out erty unchanged. With a U(1) gauge transformation of the Dirac cone. the Grassmann field, the gauge invariance requires the Acknowledgments: Authors are grateful to A. Burkov, minimal coupling k = kx Ax i(ky Ay). With fur- T. Sedrakyan, C. S. Ting, Q. H. Wang, L. Hao, L. B. ′ − − − ther including the vortex degrees of freedom in the order Shao for fruitful discussions. This work was supported parameter, the Lagrangiane can be written as, by 973 Program under Grant No. 2017YFA0303200, and by NSFC (Grants No. 60825402, No. 11574217). OE 0 + iχ ⋆ iχ = k(γ i∂ + γ e k + γ−k e− + m) k, acknowledges the support from ASU startup grant. L C 0 ′ ′ C Author Contributions: All authors performed the where χ is the vortexf phase whiche goese from 0 to 2π when calculations, discussed the results and prepared the one completes a closed loop in real space. Without los- manuscript. ing generality, we set ∆c = 1. The behavior of vortex Competing Financial Interests: The authors declare | | degrees of freedom can be extracted from the above La- no competing financial interests. grangian after integrating out the Grassmann fields. To Data Availability: The data that support the findings the second order expansion, the Chern-Simons action in of this study are available from the corresponding author Eq.(7) is obtained. upon reasonable request. Mean-field analysis of the Hubbard interaction. Let us denote the four k vectors in the Hubbard in- teraction Eq.(9) as k = k, k = k′, k = k′ q 1 2 3 − and k4 = k + q, respectively. For zero temperature, ∗ Electronic address: [email protected] only those degrees of freedom near µ become dominant. † Due to conservation of momentum in 2D, only the for- Electronic address: [email protected] [1] Redlich, A. N. Gauge noninvariance and parity noncon- ward scattering is least irrelevant for repulsive inter- servation of three-dimensional fermions, Phys. Rev. Lett. actions in the low-energy window in the renormaliza- 52, 18 (1984). tion group sense, which consists two scattering channels: [2] Redlich, A. N. Parity violation and gauge noninvariance k1 = k3 (k2 = k4) and k1 = k4 (k2 = k3). Then of the effective gauge field action in three dimensions, we make a unitary transformation into the chiral basis. Phys. Rev. D 29, 2366 (1984). For the scattering channel k = k (k = k ), after ex- [3] Semenoff, G. W. Condensed-matter simulation of a three- 1 3 2 4 53 pansion, there are sixteen different combinations in to- dimensional anomaly, Phys. Rev. Lett. , 2449 (1984). [4] Castro Neto, A. H., Guinea, F., Peres, N. M. R., tal. The requirement of conservation of energy greatly Novoselov, K. S. & Geim, A. K. The electronic properties simplifies the expression, leaving us six terms that read of graphene, Rev. Mod. Phys. 81, 109 (2009). ˆ ˆ ˆ ˆ as: (k′ k)c† ′ c† ck′,+ck,+, (k′ k)c† ′ c† ck′,+ck, , [5] Novoselov, K. S. et al. Two-dimensional gas of massless · k ,+ k,+ − · k ,+ k, − ˆ ˆ ˆ ˆ − ˆ Dirac fermions in graphene Nature, 438, 197-200 (2005). (k′ k)ck† ′,+ck,† ck′, ck,+, (k′ k)ck† ′, ck,† +ck′,+ck, , (k′ · − − · − − − · [6] Liu, Z. K. et al. Discovery of a three-dimensional topo- ˆ † † ′ ˆ ˆ † † ′ k)ck′, ck,+ck , ck,+, (k′ k)ck′, ck, ck , ck, . Only those logical Dirac semimetal, Na3Bi, Science, 343, 864-867 − − · − − − − scattering terms that involves two bands can open a gap (2014). in the Dirac state, while the terms that within one band [7] Wang, Z., Sun, Y., Chen, X.-Q., Franchini, C. et al. Dirac only renormalizes the Fermi velocity but do not con- semimetal and topological phase transitions in A3Bi (A = 85 tribute to any instabilities. Moveover, the translation Na, K, Rb), Phys. Rev. B , 195320 (2012). invariant mean-field treatment to the second and fifth [8] Burkov, A. A. & Balents L. Weyl semimetal in a topo- logical insulator multilayer, Phys. Rev. Lett. 107, 127205 term contributes to a self-energy that only renormalizes (2011). the chemical potential. With these consideration, the [9] Xu, S.-Y. et al. Discovery of a Weyl fermion semimetal forward scattering channel k1 = k3 (k2 = k4) produces and topological Fermi arcs, Science, 349, 613-617 (2015). the following reduced interaction, [10] Lv, B. Q. et al. Experimental discovery of Weyl semimetal TaAs, Phys. Rev. X 5, 031013 (2015). (1) V [11] Zhang, H., Liu, C.-X., Qi, X.-L., Dai, X. et al. Topologi- ˆ ˆ † † ′ Heff = (k k′)ck,nck, nck′, nck ,n. (16) cal insulators in Bi2Se3, Bi2Te3 and Sb2Te3 with a single − 4 · − − k k′ 5 X, ,n Dirac cone on the surface, Nat. Phys. , 438-442 (2009). [12] Chen, Y. L. et al. Experimental Realization of a Three- For the scattering channel k1 = k4 (k2 = k3), follow- Dimensional Topological Insulator, Bi2Te3, Science, 325, ing the same approach above, one obtains the following 178-181 (2009). reduced interaction, [13] Xia, Y. et al. Observation of a large-gap topological- insulator class with a single Dirac cone on the surface, Nat. Phys., 5, 398-402 (2009). (2) V Heff = ck† ,nck, nck† ′, nck′,n. (17) [14] Fu, L. & Kane, C. L. Superconducting proximity effect − 4 − − k k′ and Majorana fermions at the surface of a topological X, ,n 9

insulator, Phys. Rev. Lett. 100, 096407 (2008). [38] Wang, J. et al. Evidence for electron-electron interac- [15] Sun, H.-H. et al. Majorana zero mode detected with spin tion in topological insulator thin films, Phys. Rev. B 83, selective Andreev reflection in the vortex of a topological 245438 (2011). superconductor, Phys. Rev. Lett. 116, 257003 (2016). [39] Wray, L. A. et al. A topological insulator surface under [16] Qi, X.-L. & Zhang, S.-C. Topological insulators and su- strong Coulomb, magnetic and disorder perturbations, perconductors, Rev. Mod. Phys. 83, 1057 (2011). Nat. Phys. 7, 32-37 (2011). [17] Chen, Y. L. et al. Massive Dirac fermion on the surface [40] Raghu, S., Qi, X.-L., Honerkamp, C. & Zhang, S.- of a magnetically doped topological insulator, Science, C. Topological Mott Insulators, Phys. Rev. Lett. 100, 329, 659-662 (2010). 156401. [18] Chang, C.-Z. et al. Experimental observation of the quan- [41] Rachel, S. & Hur, K. L. Topological insulators and Mott tum anomalous Hall effect in a magnetic topological in- physics from the Hubbard interaction, Phys. Rev. B 82, sulator, Science, 340, 167-170 (2013). 075106 (2010). [19] Yu, R., Zhang, W., Zhang, H.-J., Zhang, S.-C., Dai, X., [42] Wang, Z., Qi, X.-L., & Zhang, S.-C. Topological order & Fang, Z., Quantized anomalous Hall effect in magnetic parameters for interacting topological insulators, Phys. topological insulators, Science, 329, 61-64 (2010). Rev. Lett. 105, 256803 (2010). [20] Bucher B., Steiner P., & Wachter P. Excitonic insula- [43] Slagle, K., You, Y.-Z., & Xu, C. Exotic quantum phase tor phase in TmSe0.45Te0.55, Phys. Rev. Lett. 67, 2717 transitions of strongly interacting topological insulators, (1991). Phys. Rev. B, 91, 115121 (2015). [21] J´erome D., Rice T. M., & Kohn W. Excitonic Insulator, [44] Gurarie, V. Single-particle Green’s functions and inter- Phys. Rev. 158, 462 (1967). acting topological insulators, Phys. Rev. B 83, 085426 [22] Read, N. & Green, D. Paired states of fermions in two di- (2011). mensions with breaking of parity and time-reversal sym- [45] Stoudenmire, E. M., Alicea, J., Starykh, O. A. & Fisher, metries and the fractional quantum Hall effect, Phys. M. P. A. Interaction effects in topological superconduct- Rev. B 61, 10267 (2000). ing wires supporting Majorana fermions, Phys. Rev. B [23] Qi, X.-L., Hughes, T., & Zhang, S.-C. Topological field 84, 014503 (2011). theory of time-reversal invariant insulators, Phys. Rev. B [46] Wang, C. & Senthil, T. Interacting fermionic topological 78, 195424 (2008). insulators/superconductors in three dimensions, Phys. [24] Wei, H., Chao, S.-P., & Aji, V. Excitonic Phases from Rev. B 89, 195124 (2015). Weyl Semimetals, Phys. Rev. Lett. 109, 196403 (2012). [47] Yao, H. & Ryu, S. Interaction effect on topological classi- [25] Volovik, G. E. Universe in a Helium droplet (Clarendon, fication of superconductors in two dimensions, Phys. Rev. Oxford, 2003). B 88, 064507 (2013). [26] Jiang, Y. et al. Mass acquisition of Dirac fermions in [48] Baum, Y. & Stern, A. Magnetic instability on the surface magnetically doped topological insulator Sb2Te3 films, of topological insulators, Phys. Rev. B 85, 121105(R) Phys. Rev. B 92, 195418 (2015). (2012). [27] Yoshimi, R. et al. Quantum Hall states stabilized in semi- [49] Hao, N., Zhang, P., Li, J., Wang, Z., Zhang, W., & Wang, magnetic bilayers of topological insulators, Nat. Com- Y. Chiral topological excitonic insulator in semiconduc- mun., 6, 8530 (2015). tor quantum wells, Phys. Rev. B 82, 195324 (2010). [28] Zhang, S. et al. Anomalous quantization trajectory and [50] Chen, J. et al. Gate-Voltage Control of Chemical Poten- parity anomaly in Co cluster decorated BiSbTeSe2 nan- tial and Weak Antilocalization in Bi2Se3, Phys. Rev. Lett. odevices, Nat. Commun., 8, 977 (2017). 105, 176602 (2010). [29] Hou, C. Y., Chamon, C., & Mudry, C. Electron frac- [51] Checkelsky J. G., Hor Y. S., Cava R. J. & Ong N. P. tionalization in two-dimensional graphenelike structures, Bulk Band Gap and Surface State Conduction Observed Phys. Rev. Lett. 98, 186809 (2007). in Voltage-Tuned Crystals of the Topological Insulator [30] Seradjeh, B., Moore, J. E., & Franz, M. Exciton con- Bi2Se3, Phys. Rev. Lett. 106, 196801 (2011). densation and charge fractionalization in a topological [52] Nakajima, Y., Syers, P., Wang, X., Wang, R. & Paglione, insulator film, Phys. Rev. Lett. 103, 066402 (2009). J. One-dimensional edge state transport in a topological [31] Nayak, C., Simon, S. H., Stern, A., Freedman, M. & Kondo insulator, Nat. Phys. 12, 213-217 (2016). Sarma, S. Das Non-Abelian anyons and topological quan- [53] Knolle, J. & Cooper, N. Excitons in topological Kondo tum computation, Rev. Mod. Phys. 80, 1083 (2008). insulators: Theory of thermodynamic and transport [32] Su, W. P., Schrieffer, J. R., & Heeger, A. J. Solitons in anomalies in SmB6, Phys. Rev. Lett. 118, 096604 (2017). polyacetylene, Phys. Rev. Lett. 42, 1698 (1979). [54] Roy, B., Hofmann, J., Stanev, V., Sau, J. & Galitski, [33] Su, W. P., Schrieffer, J. R. & Heeger, A. J. Soliton exci- V. Excitonic and nematic instabilities on the surface of tations in polyacetylene, Phys. Rev. B 22, 2099 (1980). topological Kondo insulators, Phys. Rev. B. 92, 245431 [34] Chamon, C. et al. Electron fractionalization for two- (2015). dimensional Dirac fermions, Phys. Rev. B 77, 235431 [55] Min, H., Bistritzer, R., Su, J.-J. & MacDonald, A.H. (2008). Room-Temperature Superfluidity in Graphene Bilayers, [35] Kurkov, M. & Vassilevich, D. Parity anomaly in four di- Phys. Rev. B 78, 121401(R). mensions, Phys. Rev. D 96, 025011 (2017). [56] Sobota, J. A., et al. Ultrafast Optical Excitation of a Per- [36] Franz, M. & Teˇsanovi´cZ. Quasiparticles in the vortex sistent Surface-State Population in the Topological Insu- lattice of unconventional superconductors: Bloch waves lator Bi2Se3, Phys. Rev. Lett. 108, 117403 (2012). or Landau levels, Phys. Rev. Lett. 84, 554 (2000). [57] Wang, Y. H., et al. Measurement of Intrinsic Dirac [37] Wang, C., Potter, A. C., & Senthil T. Classification of Fermion Cooling on the Surface of the Topological In- interacting electronic topological insulators in three di- sulator Bi2Se3 Using Time-Resolved and Angle-Resolved mensions, Science 343, 629-631 (2014). Photoemission Spectroscopy, Phys. Rev. Lett. 109, 10

127401 (2012). fermions on the surface of topological insulators, Nat. [58] Hajlaoui, M., et al. Tuning a Schottky barrier in a pho- Phys. 7, 939-943 (2011). toexcited topological insulator with transient Dirac cone [62] Sun, K., Yao, H., Fradkin, E. & Kivelson, S. A. Topo- electron-hole asymmetry, Nat. Commun. 5, 3003 (2014). logical insulators and nematic phases from spontaneous [59] Neupane M., et al., Gigantic Surface Lifetime of an symmetry breaking in 2D Fermi systems with a quadratic Intrinsic Topological Insulator, Phys. Rev. Lett., 115, band crossing, Phys. Rev. Lett. 103, 046811 (2009) 116801 (2015). [63] Hou, Y., et al. Emergence of excitonic super- [60] Triola, C., Pertsova, A., Markiewicz, Robert S., & Bal- fluid at topological-insulator surfaces. Preprint at atsky, Alexander V. Excitonic gap formation in pumped https://arxiv.org/abs/1810.10653 (2018). Dirac materials, Phys. Rev. B 95, 205410 (2017) [61] Beidenkopf, H. et al. Spatial fluctuations of helical Dirac