Prediction of a topological p + ip excitonic insulator with parity anomaly

Excitonic insulators are insulating states formed by the coherent condensation of electron and hole pairs into BCS-like states. Isotropic spatial wave functions are commonly considered for excitonic condensates since the attractive interaction among the electrons and the holes in semiconductors usually leads to s-wave excitons. Here, we propose a new type of excitonic insulator that exhibits order parameter with p + ip symmetry and is characterized by a chiral Chern number Cc = 1/2. This state displays the parity anomaly, which results in two novel topological properties: fractionalized excitations with e/2 charge at defects and a spontaneous in-plane magnetization. The topological insulator surface state is a promising platform to realize the topological excitonic insulator. With the spin-momentum locking, the interband optical pumping can renormalize the surface electrons and drive the system towards the proposed p + ip instability.

M assless Dirac particles preserve both parity and timereversal symmetry (TRS). However, an infinitesimal mass will inevitability break these symmetries through the radiative correction [1][2][3] , resulting in an Abelian or nonAbelian Chern-Simons theory of the external gauge field. This effect is called parity anomaly. Quasiparticles with linear dispersion in condensed matter physics have recently drawn broad attention [4][5][6][7][8][9][10][11][12][13] , and they provide a quintessential platform to test these effects. Indeed, the gapless Dirac fermions of topological insulator (TI) surface states have been predicted to exhibit the parity anomaly and a number of exotic phenomena when different types of mass terms are considered, including the p + ip type topological superconductor induced by the superconductor (SC) proximity effect 14,15 , and the quantum anomalous Hall effect in the TI surface with a Zeeman field [16][17][18][19] .
In this work, we expand the realm of two-dimensional (2D) topological states by devising a study on the particle--hole instability of Dirac fermions. We show that the formation and condensation of excitons can lead to a novel topological excitonic insulator that exhibits the parity anomaly. Furthermore, in sharp contrast to conventional semiconductors, the 2D helical Dirac liquid with spin-momentum locking (e.g., the TI surface), has an inherent tendency towards the formation of spin-triplet excitons. This remarkable property manifests itself most clearly by the interband optical absorption processes, as shown by Fig. 1.
For conventional semiconductors with absence of the spin-orbit coupling, the direct interband excitation generates a conduction electron and a valence hole with the same momentum k and the opposite spins (Fig. 1a), resulting in a spin-singlet exciton after they form into a pair. The situation is, however, different in TI surfaces, as shown by Fig. 1b. Due to the spinmomentum locking of TI surface, the direct interband excitation is accompanied by the spin-flip of the electron; the resultant exciton consists of an electron and a hole with the same spin polarization, and therefore a spin-triplet state. This unique optical property of TI surface, as will be shown below, can be used for generation and stabilization of a unique p + ip excitonic insulator. The natural topological p-wave excitonic insulator has p + ip symmetry inherent to their order parameters. We firstly explore the salient feature of the p + ip excitonic state by studying a minimal theoretical model. In contrary to the conventional excitonic insulators 20 described by a s-wave Bardeen-Cooper-Schrieffe (BCS)-like theory 21 , the p + ip state mathematically resembles the spinless p + ip SC 22 in the chiral basis. However, unlike the p + ip SC, it does not break the global U(1) gauge symmetry but only violates the TRS. To characterize such exotic states, we propose a chiral Chern number C c and unveil the underlying parity anomaly by deriving an effective Chern-Simons gauge theory associated with topological defects. This is further responsible for (i) fractionalized excitations with charge e/2 at vortices and (ii) a spontaneous inplane magnetization of ground state. Then, we investigate the excitonic instability of a strongly correlated TI surface model where we show that a p-wave-type interaction component emerges intrinsically, which stabilizes the p + ip excitonic phase in a significant parameter region of the phase diagram. Moreover, we also consider the effects of interband optical processes and we demonstrate that the electron-photon coupling can effectively drive the system into the proposed p + ip excitonic insulator as a result of the spin-flip mechanism demonstrated in Fig. 1b.

Results
Model. We start from a spontaneous symmetry breaking of 2D Dirac fermions with spin-momentum locking and a momentumdependent mass term, H ¼ H Dirac þMðkÞ, where k is the 2D wave vector, Φ k = [ψ k,↑ , ψ k,↓ ] T is the fermionic spinor and ψ y k;σ (ψ k,σ ) denotes the creation (annihilation) operator of electrons with spin σ. We firstly focus on a concrete mass term while under which conditionsMðkÞ can be spontaneously generated will be investigated below. θ is the polar angle of k, and Δ c is the overall gap scale. It is straightforward to verify thatMðkÞ breaks TRS,TMðkÞT À1 ≠MðÀkÞ, whereT ¼ iσ y K being the time-reversal operator, with K the complex conjugation operator. In contrast with the TI surface with a Zeeman field, where the diagonal term mσ z breaks TRS 23 , here the diagonal term Δ c cos θσ z inMðkÞ respects the TRS. The TRS is however violated by the two off-diagonal terms inMðkÞ.
One can unveil the physical meaning ofMðkÞ by making a unitary transformation onto the chiral basis with H Dirac being diagonalized: C k ¼RðkÞΦ k , whereRðkÞ is a 2 by 2 rotation matrix withR 11 , and C k = [c k,+ , c k,− ] T is the spinor composed of the conduction and valence band fermions. In the chiral basis, H becomes The off-diagonal Δ c e −iθ describes the process where an electron in the valence band is scattered up to the conduction band, leaving a particle-hole excitation on top of the many-body ground state. For Δ c ≠ 0, the particle-hole pairs condense and gap out the Dirac cone, generating an excitonic insulator. Another interesting feature of the off-diagonal term is the factor e ±iθ = (k x ± ik y )/k 24 , similar to the spinless p + ip-wave SC 25 , but we note that the order parameter resides in the particle-hole rather than the particle-particle basis. The model Hamiltonian can be analyzed similarly as the BCS superconductor theory. The ground state wave function |Ψ〉 is obtained after a Bogoliubov transformation as, with is the energy spectrum of quasiparticles. |Ψ〉 clearly describes a coherent state of electron-hole pairs with the distribution factor g k of the p + ip symmetry.
Topological classification and chiral Chern number. To examine the topological classification of the p + ip excitonic insulator, we write Eq. (2) into the compact form as,  Recalling that for the TI surface with a Zeeman field, a conventional Chern number C is defined in the same way 16,23 but in the spin basis, Φ k = [ψ k,↑ ,ψ k,↓ ] T . The difference between C and C c becomes clear from the perspective of Berry phase as following. In the original spin basis, following a closed loop in momentum space, the accumulation of Berry phase is is the eigenstate of the valence band in the spin basis. For a Bloch state defined in the chiral basis, a similar geometric phase can be defined as is the eigenstate of the valence band in the chiral basis. The Bloch states in the above two formalisms are related to each other by the unitary transformation u c k ¼RðkÞ u k j i. It then becomes obvious that the two Berry phases satisfy ϕ c − ϕ = Ài H dk Á u k h jR y ðkÞ∇ kR ðkÞ u k j i, leading to C c − C = Àði=2πÞ H dk Á u k h jR y ðkÞ∇ kR ðkÞ u k j i= 1/2, after inserting theRðkÞ matrix we used to diagonalize H Dirac .
To explicitly demonstrate the topological invariants, we plot the normalized vectordðkÞ with traversing the whole k-space. Figure 2a, b shows the mapping of thedðkÞ vector in the spin and chiral basis, respectively. It is found thatdðkÞ in latter basis completely covers the hemisphere whiledðkÞ in the former does not, leading to C c = 1/2 and C = 0. Therefore, two topological invariants are needed to completely characterize the p + ip excitonic state, namely, (C, C c ) = (0, 1/2). The zero Chern number C = 0 means that the p + ip excitonic insulator is topologically distinct from the TI surface with a Zeeman field that has C = 1/2 16 . This is clear from the Landau quantization of the two states. For the TI surface with a Zeeman field, the Landau level (LL) for |n| ≥ 1 reads as , where m is the Zeeman gap. In addition, there is an unconventional zeroth LL E 0 = sgn(m)m which is located exactly at the gap edge [26][27][28] . This zeroth LL is a manifestation of the nontrivial Zeeman mass as a result of C = 1/2 27,28 . For the p + ip excitonic insulator phase, however, one cannot find the unconventional zeroth LL located at the gap edge, consistent with C = 0.
Chern-Simons theory and fractionalization. The nonzero chiral Chern number C c = 1/2 indicates that the p + ip excitonic insulator must be topologically distinct from normal insulators.
Recalling the p + ip SCs where Majorana fermions are encoded in vortices 22 , we shift our attention to topological defects in the p + ip excitonic state. To this end, we can solve the Bogoliubov de Gennes (BdG)-type equations looking for zero-energy solutions. In real space, the equations can be derived as (see Methods for details): where α is the polar angle in real space and u, v are the two components of the spinor wave function. Similar equations have been studied in refs. 22,29,30 . However, note that the physics in the present model is different from previous works. Read and Green 22 derive the equations for SCs, where charge conservation is broken and Majorana fermions are the emergent excitations. For our system, the U(1) gauge symmetry is preserved, therefore no Majorana modes are allowed. Hou et al. 29 investigate graphene with a chirality mixing texture which induces the mass term, and Seradjeh et al. 30 discuss the exciton condensation between two opposite TI surfaces. In contrary, here we focus on the excitonic order in a Dirac liquid state and the mass term is generated by the condensation of p + ip excitons. To solve the equations in the presence of a topological defect 31 , we assume Δ c (r) = Δ c (r)e iα . A unique normalized zero-energy bound state 32,33 can be obtained as, u = Af(r)e i(α+π/2)/2 and v = u * , where A is the normalization constant and the function f(r) = exp[−|Δ c |r/|v F |] decays exponentially. We note that the single particle Hamiltonian h in Eq. (2) satisfies (iτ 2 )h(iτ 2 ) = h. This allows one to define a time-reversaltype operator Θ = iτ 2 K, with which it is clear that for any eigenstate ψ E with energy E, Θψ E is also a eigenstate with energy −E. Besides, the total number of fermion states are conserved, namely, is the change of density of states. Applying this formula to the case where a topological defect is in presence, and then taking into account the fact that ρ (E) is a symmetric function with respect to E = 0, one comes to the conclusion that both the conduction and valence band lose half a state. This state is compensated by the zero mode that is bound to the topological defect. Consequently, at half filling, the electron charge associated to the defect is either −e/2 or e/2, depending on whether the zero mode is occupied or empty.
We have shown that the p + ip excitonic insulator is characterized by C c = 1/2 and displays fractional charge e/2 at vortices. However, it is still unclear how the two observations are related to each other. This can be answered by the parity anomaly 34 of the p + ip excitonic insulator and the effective gauge theory associated with the topological defects. First, we couple an external U(1) gauge field A μ (with μ = 0, 1, 2) to the p + ip excitonic state 35 , and we introduce the Dirac matrices defined as γ 0 = τ 3 , γ 1,2 = iτ 1,2 . In order to extract the effective theory for the vortices, we include the vortex degrees of freedom χ in the order parameter. Second, we perform a singular gauge transformation, c ± → e ±iχ/2 c ± 36 , which results in the following Lagrangian (see Methods for details), where C k C k À Á is the Grassmann field.γ μ are defined asγ 0 ¼ γ 0 , γ 1 ¼ γ 2 ,γ 2 ¼ Àγ 1 . In Eq. (6), the vortex degrees of freedom have been absorbed into an effective gauge field a μ after the singular transformation. For a static defect with ∂ 0 χ = 0, a μ ¼ A μ þ 1 2 ∂ μ χ. Finally, after one integrates out the Grassmann field, a Chern-Simons term with respect to the effective gauge field a μ can be obtained in the second order expansion: where the coefficient I is found to be 1/2 for the proposed p + ip excitonic insulator. We note that for the TI surface with a Zeeman field 23 , a similar Chern-Simons Lagrangian arises as a result of the parity anomaly, where the coefficient in front of the Chern-Simons term is exactly the Chern number C. Here, the coefficient I = 1/2 is identified to be the chiral Chern number C c defined above. The similarity between the two gauge field theories is due to the underlying parity anomaly; however, for the p + ip excitonic insulator, it is the effective gauge field a μ rather than the external field A μ that obeys the Chern-Simons form. Further inserting consists of a Lagrangian describing the vortex L v and a current operator associated with the vortex, For a static vortex (with vortex number n = 1) that centered at r 0 , we have ϵ ij ∂ i ∂ j χ ¼ 2πδ r À r 0 ð Þ, inserting which into Eq. (8), the electron charge of the defect is then obtained as . Therefore, we have proved that the chiral Chern number C c equals to the fractional charge (in unit of e) of a topological defect, justifying C c to be the essential topological invariant characterizing the p + ip excitonic insulator.
Spontaneous in-plane magnetization. The order parameter −Δ c e −iθ in Eq. (2) breaks TRS. This suggests that the ground state carries a magnetization associated with it. Before proceeding, we note that the order parameter −Δ c e −iθ still enjoys a redundant global phase. In the following, we take this into account and write the off-diagonal term in Eq. (2) as −Δ c e −iθ′ , where θ′ = θ − β with β being a constant angle independent of k. The global phase β is analogous to U(1) phase of a superconductor.
To clearly show the feature of magnetization, we calculate the spin expectation values 〈σ x 〉, 〈σ y 〉, 〈σ z 〉 from the single particle eigenstates and study the spin texture of the p + ip excitonic ground state |Ψ〉 (Eq. (3)). We obtain As has been discussed before, 〈σ z 〉 does not break TRS. The integral of 〈σ z 〉 over the whole k-space vanishes. The in-plane spin texture (〈σ x 〉, 〈σ y 〉) is however nonzero, and we plot it in Fig. 3 for different values of β, i.e., β = 0, β = π/2, β = π, β = 3π/2. Interestingly, it is found that the in-plane spin configuration is tilted from the standard helical spin texture of the original TI surface, and the spins at two TRS-related momentum points, k and −k, are no longer always opposite to each other, violating the TRS. Though being tilted, the chiral nature of electron's spin is preserved as we can see that the direction of the vector (〈σ x 〉, 〈σ y 〉) still undergoes a change of 2π when traversing a closed path around k = 0 point. Furthermore, we calculate the total in-plane magnetization M by integrating (〈σ x 〉, 〈σ y 〉) over the momentum space. It is found that M always points towards the opposite direction of the global phase β, as shown by the thick red arrows in Fig. 3. We also note that the emergent magnetization of ground state can lead to anisotropic in-plane transport of excitons which serves as an evident experimental signature of the p + ip excitonic insulator. Furthermore, since β determines the direction of magnetization of ground state, one expects to control β by applying a weak in-plane magnetic field during the process of TRS symmetry breaking.
The excitonic phases on interacting TI surfaces. After careful characterization of the topological properties of the p + ip excitonic insulator, we now study its possible realization on interacting TI surface states. Although various instabilities of interacting topological insulators and topological superconductors have been extensively studied [37][38][39][40][41][42][43][44][45][46][47][48] , we focus on the excitonic instability of the surface states and propose to realize the minimal model of the topological excitonic insulators via spontaneous symmetry breaking that enjoys a fractional chiral Chern number C c = 1/2, which is different from earlier proposals based on semiconductor heterostructure with an integer Thouless-Kohmoto-Nightingale-Nijs (TKNN) number 49 .
The effective model Eq. (1) describing the surface state is usually modified by a higher order wrapping effect, reducing the rotational symmetry of the Dirac cone down to a discrete C 3 subgroup. Our further studies on the wrapping effect shows that the exciton gap remains very robust and therefore one can neglect the wrapping effect for the following analysis. We firstly consider the particle-hole symmetric case, while the generalization to finite chemical potential μ is very straightforward and will be discussed below.
First we start with demonstrating that a p-wave-type interaction that spontaneously generates the p + ip-wave exciton insulator naturally arises from the Hubbard model, which reads in k-space as: where we assume a momentum cutoff Λ implicit in the sum of k. n is the band index. The interaction term in the second line, H and it arises from the spin-momentum locking of the surface state. After introducing a complex bosonic vector field Δ and performing the Hubbard-Stratonovich transformation to decouple H ð1Þ eff , we arrive at a mean-field theory with the order parameter given by Δ ¼ 1 2 V P kk c y k;À c k;þ D E . When Δ ≠ 0, the condensation of excitons gaps out the Dirac point, resulting in the excitonic insulator phase. Similar to 3 He 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 ffiffi 2 p Δ c b e x À ib e y and (2) the polar excitonic insulator with Δ ¼ Δ p b e, where b 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 equation, we obtain the condensation energy and the order parameter (Δ c , Δ p ), which are plot in Fig. 4a, b, respectively. Figure 4b shows the emergence of a threshold V c , beyond which the system develops the excitonic instabilities with |Δ c/p | ≠ 0. Figure 4a indicates that even though both the p + ip excitonic insulator and polar excitonic insulator are stable compared to the gapless Dirac cone for V > V c , the p + ip excitonic insulator is more energetically favored due to the lower condensation energy. Moreover, we also performed addition mean-field calculations by treating the p + ip and polar excitonic instability on equal footing. It is found that the p + ip channel always suppresses the polar order down to zero. Therefore, even though the ground state energies of the two states are close, the mass term of the polar exciton insulator is negligible in the mean-field theory. Hence, the resultant meanfield Hamiltonian exactly matches with Eq. (2), and the proposed mass termMðkÞ is generated by spontaneous TRS breaking as a result of the p-wave interaction H ð1Þ eff . Moreover, we note that for particle-hole asymmetric case with a finite chemical potential μ, the critical V c increases with increasing μ. Although the exciton insulator is always stable for strong enough interaction V > V c , it is more experimentally favorable to tune μ close to the Dirac point 50,51 . In the above calculations, we temporarily focused on the pwave interaction H ð1Þ eff without considering the s-wave component. Nevertheless, the results clearly show that the TI surface state indeed has a strong tendency towards the p + ip excitonic instability via spontaneous TRS breaking. We now take into account the interaction in the s-wave channel, which can also be derived from the local Hubbard interaction as (see Methods for details), H ð2Þ eff ¼ À V 4 P k;k′;n c y k;n c k;Àn c y k′;Àn c k′;n . In order to be more general, we assume in the following two independent couplings V p and V s for the p-wave and s-wave interaction respectively, and treat both terms on equal footing. An s-wave bosonic scalar field eff in addition to the p + ip excitonic order parameter Δ c . The coupled self-consistent equations with respect to Δ s and Δ c are obtained in the standard approach which leads to the calculated phase diagram shown by Fig. 4c. It is found that both the s-wave and p + ip-wave excitonic insulator can be spontaneously generated from the Dirac liquid phase with strong enough interactions. The former and latter are stabilized for large V s and V p respectively. We note that although the p + ip excitonic phase, which is absent in conventional semiconductors, becomes stable in a large parameter region of the phase diagram, it does not occur as the ground state for the discussed local Hubbard model, which lies in the parameter line V s = V p that only enters into the s-wave excitonic condensation. The stabilization of the p + ip excitonic insulator. Now we study the stabilization of the p + ip excitonic insulator on TI surfaces by going beyond the natural Dirac semimetals and taking into account external perturbations. Recalling the interband optical absorption process shown by Fig. 1b, one expects that additional photons should prefer the p-wave pairing due to the spin-flip mechanism that gives rise to spin-triplet excitons. Since excitons are usually generated by absorption of photons in a material, we consider a linearly polarized perpendicular incident laser with frequency ω and light intensity I. The corresponding microscopic theory can be constructed by quantization of the vector potential A(r), i.e., where a q,λ a y q;λ is the annihilation (creation) operator of photons, e λ (λ = 1, 2) are the polarization vectors, ω q = c|q| is the light frequency, with c being the light velocity and q the wave vector of photons.
The coupling of the external field to the surface electrons is described within the formalism of the minimal coupling, which leads to the electron-photon interaction, Þ ψ y k;σ a q;λ ψ kÀq;σ′ þψ y k;σ a y Àq;λ ψ kÀq;σ′ ; where the term (σ σσ′ · e λ ) arises due to the helical spin texture and the spin-momentum locking of TI surface states. For linearly polarized perpendicular incident photons, e 1 ¼ b e x and e 2 ¼ b e y . Thus, the spin-flip matrices σ x and σ y take place in H el-ph . This verifies the spin-flip mechanism illustrated by Fig. 1b in a microscopic perspective. Moreover, since c ) v F , the conservation of energy and momentum requires k j j ) q j j in H el-ph , and therefore ensures the direct interband excitation with k − q ≃ k. Furthermore, the photon fields are Gaussian with the Hamiltonian H ph ¼ P qλ ω q a y q;λ a q;λ þ 1=2 , one can therefore exactly integrate them out. This procedure leads to the effective action S int describing the renormalization from photons: Where P λ σ αβ Á e λ σ ργ Á e λ ¼ σ x αβ σ x ργ þ σ y αβ σ y ργ , inserting which it is found that S int consists of an equal-time renormalization term that produces the effective interaction between electrons, where V eff ¼ 4πe 2 v 2 F n ph = ω 2 and n ph is the number density of photons that interact with the electrons, which is proportional to the light intensity I. We have utilized the fact that all the wave vectors q satisfy q ¼ qb e z ¼ ωb e z =c for photons with the frequency ω. The first  Fig. 4 The results from mean-field calculations. a Condensation energy versus interaction V for p + ip excitonic insulator and polar excitonic insulator. b Calculated mean-field excitonic order parameter for p + ip excitonic insulator (Δ c ) and polar excitonic insulator (Δ p ) versus interaction V. Δ c=p 2 =v 2 F Λ 2 ( 1 is assumed. c The phase diagram based on unbiased mean-field with considering both the sand p + ip-wave interaction channel. The dash-dotted red line denotes the Hubbard interaction case with V s = V p two operators in H int , ψ y k;σ ψ k;σ , are generated by the optical interband absorption process with spin-flip depicted by Fig. 1b. Similarly, the last two operators, ψ y k′; σ ψ k′;σ , are originated from the interband relaxation process with spin-flip. After a unitary transformation to the chiral basis, Eq. (14) is exactly cast into: Interestingly, H int enjoys the p-wave factork Ák′, exactly the same as the p-wave interaction component reduced from the Hubbard model, H eff . The interband optical processes automatically select the p-wave channel due to the spin-flip scatterings (Fig. 1b). More importantly, H int has a strength V eff tunable by external field parameters, ω and I. According to the mean-field study above and the calculated phase diagram in Fig. 4, it is known that the tunable external field can renormalize the Dirac cone state and drive it into the proposed p + ip-wave excitonic insulator. The resultant state displays the predicted fractional charges associated with vortex excitations, the spontaneous in-plane magnetization and thus the anisotropic transport of excitons on the TI surface.

Discussion
Excitonic insulators require strong interactions since the topological surface states become unstable only when the interaction is as large as the ultraviolet cutoff (band gap) as shown in Fig. 4b. Therefore, we argue that topological Kondo insulators (TKIs) are good candidates for excitonic instability on TI surface. For instance, in the case of SmB 6 , which is the archetypal TKI, the renormalized Coulomb interaction (U~1 eV) is much larger than the direct band gap (Λv F~2 0 meV). For these parameters, our mean-field theory predicts an excitonic gap of order Δ s/c ≃ 2-10 meV. Indeed, magnetism and hysteresis in magnetotransport have been observed in SmB 6 52 . Whether these experiments indeed observe an excitonic insulator require future work. We would like to note that our theory is not related with the recent works on excitonic contributions to quantum oscillations 53 or the intervalley excitonic instabilities 54 in SmB 6 .
In addition to the short-range interaction, we also considered a long-range Coulomb interaction on the TI surface and make generalization to finite temperature 55 . It is found that the critical interaction V c for BCS exciton condensation is significantly reduced for long-range interaction. At finite temperature, taking into account the thermal fluctuation, the strict long-range order is absent, the Kosterliz-Thouless transition temperature T KT can be calculated by evaluating the phase stiffness of the exciton order based on the mean-field calculations 56 . Experimentally, the driven TI surface offers a rare platform for realization of new optical devices and facilitates the study of possible collective instabilities of Dirac electrons. Experiments measuring bulk metallic TI materials have obtained the lifetime of the excited states to be around a few picoseconds [57][58][59] , which is significantly larger than than that in graphene. Moreover, a remarkably long lifetime exceeding 4 μs has been observed by experiment for the surface states in bulk insulating Bi 2 Te 2 Se 60 . These indicate promising observation of a transient or quasi-equilibrium excitonic insulator 59,61 . Besides, the surface has a reduced dielectric constant which can be tuned even smaller by doping or gating 62 , and this possibly makes the long-range interaction more effective.
To summarize, we proposed a new type of topological p + ip excitonic insulator that achieves the parity anomaly and is characterized by the nontrivial chiral Chern number. We provided a concrete example and discussed its possible realization in strongly correlated TI surface. Interestingly enough, similar to the TIs due to band inversion, the predicted topological behaviors, i.e., the charge fractionalization and the emergent magnetization of ground state, are robust and immune to continuous deformation of the Hamiltonian, as long as the exciton gap remains unclosed. This suggests that realizations of topological pwave excitonic insulators are likely in many other physical systems. For example, the spectrum under study can be continuously deformed from the Dirac cone to a quadratic band crossing point (QBCP) 63 . This indicates that the instability is also possible and may be more feasible in interacting 2D electron gas with spinmomentum locking and QBCP. In this sense, the present results are general and may provide new physical settings for the realization of quantum anomalies in condensed matter system, and may stimulate different directions to search for more exotic topological state of matters. We also note that a laser-induced non-equilibrium version of our proposal has been experimentally observed recently in Bi 2 Se 3 55 , where a millimeter-long transport distance of excitons up to 40 K is reported, supporting an excitonic condensation state formed on the TI surface. Our study on the p + ip excitonic insulator stabilized by the optical pumping further predicts an anisotropic transport of excitons, which we believe would be of great interest for future experimental investigations.

Methods
The zero-energy BdG-type equations. From Eq. (2) in the main text, the zeroenergy eigen-equations can be written down as: Here, we assume a real-space dependence of the order parameter Δ(r). With a vortex,we have Δ(r) = Δ c (r)e iα , where α is the polar angle of r in real space. To simplify the eigen-equations, one can multiply a factor e iθ from the right in the first equation, and a factor e −iθ from the left in the second equation. Then, a Fourier transformation yields, After transformation to the polar coordinate, we arrive at the zero-energy BdG-like equations in real space.
Effective gauge field theory. In the second quantized form, the effective Hamiltonian of the p + ip excitonic insulator reads as Formally, one can make transformation C k ! C k ¼ C k = ffiffi ffi k p , while keeping the total Hamiltonian unchanged. The action in the functional representation is then cast into: where we have defined k′ = k x − ik y and introduced the Dirac matrices in 2D, γ 0 = τ 3 , γ 1 = iτ 1 , γ 2 = iτ 2 , and γ ± = (γ 2 ± iγ 1 )/2. The above action formally describes a 2D Dirac fermion with a mass v F k 2 . Since the sign of mass determines the topology of 2D Dirac systems, and v F k 2 ≥ 0 for all k, one can continuously deform the mass into a constant m while keeping the topological property unchanged. With a U(1) gauge transformation of the Grassmann field, the gauge invariance requires the minimal couplingk′ ¼ k x À A x À i k y À A y . With further including the vortex degrees of freedom in the order parameter, the Lagrangian can be written as, where χ is the vortex phase which goes from 0 to 2π when one completes a closed loop in real space. Without losing generality, we set |Δ c | = 1. The behavior of vortex degrees of freedom can be extracted from the above Lagrangian after integrating out the Grassmann fields. To the second order expansion, the Chern-Simons action in Eq. (7) is obtained.
Mean-field analysis of the Hubbard interaction. Let us denote the four k vectors in the Hubbard interaction Eq. (9) as k 1 = k, k 2 = k′, k 3 = k′ − q and k 4 = k + q, respectively. For zero temperature, only those degrees of freedom near μ become dominant. Due to conservation of momentum in 2D, only the forward scattering is least irrelevant for repulsive interactions in the low-energy window in the renormalization group sense, which consists of two scattering channels: k 1 = k 3 (k 2 = k 4 ) and k 1 = k 4 (k 2 = k 3 ). Then we make a unitary transformation into the chiral basis. For the scattering channel k 1 = k 3 (k 2 = k 4 ), after expansion, there are 16 different combinations in total. The requirement of conservation of energy greatly simplifies the expression, leaving us six terms that read as: ðk′ ÁkÞc y k′;þ c y k;þ c k′;þ c k;þ , Àðk′ ÁkÞc y k′;þ c y k;À c k′;þ c k;À , ðk′ ÁkÞc y k′;þ c y k;À c k′;À c k;þ , ðk′ ÁkÞc y k′;À c y k;þ c k′;þ c k;À , Àðk′ ÁkÞc y k′;À c y k;þ c k′;À c k;þ , ðk′ ÁkÞc y k′;À c y k;À c k′;À c k;À . Only those scattering terms that involve two bands can open a gap in the Dirac state, while the scatterings within one band only renormalize the Fermi velocity but do not contribute to any instabilities. Moreover, the translation invariant mean-field treatment to the second and fifth term contributes to a self-energy that only renormalizes the chemical potential. With these consideration, the forward scattering channel k 1 = k 3 (k 2 = k 4 ) produces the following reduced interaction, For the scattering channel k 1 = k 4 (k 2 = k 3 ), following the same approach above, one obtains the following reduced interaction,

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.