a topological insulating ( ti ) phase is a new state of quantum matter featuring massless dirac - like boundary states whose robustness is guaranteed by the time - reversal symmetry ( trs)@xcite . in this rapidly growing field ,
search for new topological phases has recently turned to extending the consideration of symmetry - protected states to include non - spatial symmetries@xcite such as , the particle - hole and chiral symmetries as well as the spatial symmetries@xcite . in particular , the nontrivial band topology can be shown to be protected by certain crystal symmetries , leading to the new class of tis called topological crystalline insulators ( tcis)@xcite . like the tis , a three dimensional ( 3d ) tci is also predicted to host metallic surface states on surfaces which preserve appropriate crystal symmetries@xcite .    to date , the semiconducting pb@xmath0sn@xmath1(te , se ) alloys are the only experimentally realized tci materials , whose topological nature has been verified through direct observation of dirac - like surface states via angle - resolved photoemission spectroscopy ( arpes ) experiments@xcite .
these gapless surface states are protected by the mirror symmetry with respect to the ( 110 ) or ( @xmath510 ) lattice plane .
the characteristic features of the associated dirac states such as linear dispersion , lifshitz transition , spin / orbital texture , among others , have been examined via scanning tunneling microscopy / spectroscopy ( stm / sts ) and transport measurements@xcite . moreover , it has been found that , below a critical temperature , which depends on chemical composition , the cubic lattice structure can be distorted by strains , resulting in an orthorhombic or a rhombohedral structure@xcite .
remarkably , in a recent stm experiment by okada _
et al._@xcite , two unexpected landau levels have been observed for the ( 001 ) surface of pb@xmath0sn@xmath1se under a perpendicular magnetic field .
the presence of these two extra levels is believed to be associated with the gap opening of the two surface dirac cones induced by a ferroelectric - like lattice distortion , which breaks the corresponding mirror symmetry .
it is clear that it is important to understand effects of strains on the electronic structure of tcis in order to gain a deeper handle on the nature and origin of their topological states , and how these states could be manipulated for practical applications of the tcis .    despite some strain - related studies in the literature@xcite , we are not aware of a systematic investigation of how spin - textures and other properties of dirac states in the tcis evolve under various symmetry breaking strains .
so motivated , here we examine strain effects on the ( 001 ) surface of the pb@xmath0sn@xmath1(te , se ) tcis .
our analysis is based on an effective four - band model , which is shown to capture all essential features of the topological surface states .
the form of strain related perturbations on the hamiltonian is clarified through general symmetry considerations , allowing us to delineate how characteristic features of the topological states in tcis evolve under strains .
significantly , our analysis shows that certain mirror - symmetry - breaking strains induce hedgehog spin texture with out - of - the - plane spin - polarization at the dirac point .
such a spin texture of dirac states has been reported previously in a 3d ti , but only in the presence of an exchange field@xcite .
our findings thus identify a possible new pathway for realizing tci - based spintronics devices without requiring a trs breaking field@xcite .
moreover , we show that the interplay of strains and applied zeeman field can be used to tune the chern number of the surface ground state , which is a topological invariant characterizing a quantum anomalous hall ( qah ) or a quantum hall insulator .
electron - electron interactions effects are ubiquitous in condensed matter systems , especially when charge screening is relatively poor as is the case generally for dirac - like states in graphene@xcite or tis and tcis.@xcite in this connection , we discuss possible symmetry breaking orders generated under short - range repulsion @xmath6 as a first step toward understanding correlation driven effects on dirac states in the presence of applied strains .
we numerically obtain the zero - temperature phase diagram as a function of @xmath6 and the particle density @xmath7 and delineate its evolution with increasing strength of strain .
our results not only give insight into the nature of competing orders , but also build the foundation for what may be called `` straintronics '' applications driven by electron correlations .    the paper is organized as follows . in sec .
ii , we briefly discuss the effective four - band model for the ( 001 ) surface states in the low - energy regime , and its extensions for addressing strain effects in tcis .
iii turns to consider strain effects on properties of the dirac states , including the interplay between applied strains and perpendicular zeeman fields .
these results then allow us to investigate in section iv the possible interaction - driven competing orders in dirac states in the absence as well as presence of strains .
finally , section v comments on potential applications and implications of our study , and concludes with a summary of our results .
we start by reviewing the four - band @xmath2 model for the ( 001 ) surface states in a tci developed previously in ref . , and discuss its generalization to account for effects of strains . on the ( 001 ) surface
, the low - energy surface states can be viewed as two sets of interacting coaxial dirac cones , originating from the interface between the inverted bands of the tci and the vacuum : one set is centered at @xmath3 and the other at @xmath4 in the surface brillouin zone ( sbz ) , see fig .
[ fig : lattice ] ( b ) . defining [ 110 ] , [ @xmath510 ] , and [ 001 ] as @xmath8 , @xmath9 , and @xmath10 directions
[ see fig .
[ fig : lattice ] ( a ) ] , respectively , the effective model @xmath11 around @xmath3 must obey the following three essential symmetries which leave @xmath3 invariant : the mirror reflection about the @xmath12-plane , the mirror reflection about the @xmath13-plane , and time - reversal symmetry with the corresponding symmetry operations , i.e. m_xzh_|x(k_x , k_y)m_xz^-1=h_|x(k_x ,-
k_y ) , + m_yzh_|x(k_x , k_y)m_yz^-1=h_|x(-k_x , k_y ) , + th_|x(k_x , k_y)t^-1=h_|x(-k_x ,-
k_y ) , where @xmath14 is measured from @xmath3 .    informed by the orbital characters of surface states as revealed by first - principles calculations,@xcite one may choose , for instance in the case of snte , @xmath15 , @xmath16 , latexmath:[$    @xmath11 because one of the interacting dirac cones is @xmath18-orbital ( sn ) rich , while the other is @xmath19-orbital ( te ) rich . as a result , the symmetry operations can be represented by the matrices , @xmath20 , @xmath21 , and @xmath22 , where @xmath23 with the pauli matrices @xmath24 and @xmath25 acting on orbital and spin spaces , respectively , and @xmath26 denotes complex conjugation .
these considerations restrict @xmath11 to the form h_|x&= & m_30+m^_22+k_x(v_1x_02+v_2x _
20+v_3x_32 ) + & + & k_y(-v_1y_01+v_2y_13-v_3y_31),[h ] [ eq : surfaceh ] up to first order in @xmath27 .
( @xmath9 ) denotes the direction along [ 110 ] ( [ @xmath510 ] ) .
( b ) surface brillouin zone of snte . without strain
, there are four gapless dirac cones , protected by the mirror symmetry either with respect to @xmath28 or @xmath29.,width=302 ]    for @xmath30=0 and @xmath31=@xmath32=@xmath33=@xmath34=0 , @xmath11 becomes block diagonal , reflecting the structure of the two underlying parent dirac cones with two associated dirac points , @xmath35 ( @xmath18-rich ) and @xmath36 ( @xmath19-rich ) at @xmath3 [ see fig .  [
fig : b1 ] ( a ) ] . a non - vanishing value of any of the parameters @xmath37 , @xmath31 , @xmath32 , @xmath33 and @xmath34 induces interaction between the two cones , leading to two immediate consequences : ( 1 ) the original dirac points are shifted to @xmath38 ; and ( 2 ) all degeneracies along the intersection of the two cones are lifted except for the two points of a time - reversal symmetric pair on the @xmath28 line , indicating the emergence of two child dirac cones protected by the mirror reflection about the @xmath12-plane [ see fig .  [
fig : lattice ] ( b ) ] ; hereafter , referred to as the low - energy dirac points associated with @xmath39=0 .
our effective four - band @xmath2 model correctly capture two key features of surface states of tcis .
firstly , for @xmath40 topology of the constant energy contours changes from two separate fermi circles at low energies to two concentric ellipses at high energies .
this lifshitz transition@xcite indicates the presence of van hove singularities ( vhss ) in the underlying electronic spectrum , which are located along the @xmath41 line at e_vhs+ = |m^|,= ( 0 , ) , if @xmath31 , @xmath32 , @xmath33 , and @xmath34 are neglected .
a similar situation arises for @xmath42 with @xmath43 . in all cases
, we have a logarithmically diverging density of states proportional to @xmath44 with @xmath45 and @xmath46 representing an energy scale away from the vhs and a cutoff energy for surface states , respectively .
secondly , there is no out - of - the - plane spin component [ fig .
[ fig : b1 ] ( b)-(c ) ] by using symmetry arguments in that the net out - of - the - plane spin polarization , @xmath47 , yielding @xmath48 .
we turn next to generalize our four - band effective model to include low - energy states throughout the sbz in the presence of inter - cone interactions .
this can be done by noticing that the surface states in the vicinity of @xmath3 and @xmath4 are related by @xmath49 symmetry , so that the effective model around @xmath4 can be obtained explicitly by a @xmath49 rotation , h_|y(k_x , k_y)=_4 h_|x(k_y ,-
k_x)_4 ^ -1 , where @xmath50 and the @xmath19 orbital for te atoms in the basis states is now replaced by @xmath51 orbital in @xmath52 .
the total surface hamiltonian then is h_(001)=h_|x(k_x , k_y)h_|y(k_x , k_y ) .
although @xmath53 describes an even number of dirac cones at low energies like the case of a weak ti@xcite , note that the four dirac points in a tci do not locate at time - reversal invariant momenta .
this key difference can lead to a rather different phase diagram in the presence of the electron - electron interactions@xcite in a tci , as discussed in sec .
we now discuss effects on the surface states due to strains , which could be either intrinsic or extrinsic .
since the gapless surface states in the pb@xmath0sn@xmath1(te , se ) tcis are mainly protected by mirror symmetries , a perturbation , which breaks one of these symmetries can be expected to at least partly gap the surface spectrum , and modify the spin textures around the associated ( massive ) dirac points .
although this is generally true , we will see below that this is not always the case due to other symmetry considerations . for this purpose
, we will carry out a systematic analysis of strain induced effects along the lines of sec .
ii above .        a general strain can be described in terms of a symmetric strain tensor @xmath54 , where @xmath55 denotes the displacement field and @xmath56 are chosen to coincide with the principal crystal axes as defined in sec .
any effective , strain - induced perturbation term corresponding to @xmath57 will transform in the same way as the strain tensor under time - reversal and mirror symmetries . by examining the transformation properties of @xmath58 up to first order in @xmath27 ,
the possible resulting perturbations around @xmath3 are listed in table [ t1 ] .    [ cols="^,^,^,^,^,^,^,^,^,^,^,^,^ " , ]     as expected , the strain induced perturbations in the first four columns of table i can open gaps at both low - energy dirac points along @xmath28 due to @xmath59 broken symmetry . more significantly , the resulting spin texture becomes hedgehog - like at low energies , similar to the spin reorientation phenomenon found on the surfaces of manganese - doped bi@xmath60se@xmath61 thin films@xcite .    as a concrete example , let us add the perturbation @xmath62 in eq.([eq : surfaceh ] ) , where @xmath63 denotes the electron - phonon coupling strength .
clearly , in fig .
[ fig : f1](a ) , a gap opens at the dirac point with its magnitude proportional to @xmath64 .
the spin - texture of the conduction band of the massive dirac cone is also depicted in fig .
[ fig : f1](b ) . as
the absolute value of the eigenenergy @xmath65 becomes smaller ( _ i.e. _ , closer to the bottom of the upper cone ) , the induced out - of - the - plane spin component , _ 03= , [ eq : sz ] becomes larger , where we have set @xmath66 for simplicity .
several points should be noted here as follows .
firstly , since the strain perturbation does not break trs , the spin texture for the two massive cones around @xmath3 must form time - reversed partners .
this can be inferred from eq .
( [ eq : sz ] ) by changing the sign of @xmath67 .
secondly , the out - of - the - plane spin component of the lower cone is basically opposite to that of the upper cone at the same @xmath27 .
finally , the whole spin texture resulting from the strain differs from the case where the hedgehog - like texture is induced by the perpendicular zeeman field when both massive dirac cones around @xmath3 are considered , see fig .
[ fig : f1](c ) .    other strain perturbations in table [ t1 ] , which break either the @xmath13-mirror symmetry or both the mirror symmetries , neither open a gap at the dirac point nor induce an out - of - the - plane spin component .
the robustness of this gapless dirac point originates from an underlying symmetry : the former type of perturbation is due to the presence of @xmath59 symmetry , while the latter type ( shear deformation ) is protected by a `` space - time '' symmetry , @xmath68 , with a rotation @xmath69 followed by a time - reversal operation .
the space - time protection now allows persistence of gapless dirac points sitting at generic @xmath70 points away from the mirror line @xmath28 .
this special feature has been noted previously in refs . , and it leads to a new type of tci@xcite .
vanishing out - of - the - plane spin polarization can be proven via symmetry arguments .
for instance , consider the surface states around @xmath3 under the strain perturbation @xmath71 .
the corresponding hamiltonian now reads : @xmath72 . assuming @xmath73 is an eigenstate of @xmath74 , the @xmath12-mirror symmetry guarantees that 0=|[_02,_|x]|=2iv_1yk_y |_03| .
thus , the out - of - the - plane spin component @xmath75 vanishes everywhere around @xmath3 .
other perturbations can be analyzed in a similar manner .
there other cases deserve some comment .
the first two cases involve a uniform expansion and a uniaxial stretch ( @xmath49 breaking ) , which correspond to @xmath76 and @xmath77 , respectively@xcite . because both these cases respect trs as well as the mirror symmetries ,
their net effect is only to renormalize the parameters in the original @xmath11 .
consequently , these perturbations do not open a gap at dirac points , although positions of the dirac points could shift in opposite directions along the @xmath28 line .
the third case is an experimentally observed ferroelectric - like distortion@xcite , in which two kinds of atoms are displaced along a certain direction in an opposite manner . denoting the displacement vector @xmath78 , the non - vanishing component of this distortion
, @xmath79 ( @xmath80 ) , preserves the mirror symmetry with respect to the principal axis @xmath8 ( @xmath9 ) , but breaks the rotation and the mirror symmetries along a perpendicular direction . to the zeroth order in @xmath70 around @xmath3 , the perturbation due to such distortion
can be straightforwardly shown to have the form : h_f^= _ fxd_x_10+_fyd_y_23 , where @xmath81 and @xmath82 denote the coupling strengths .
in fact , the former term is similar to the effect of strain @xmath83 , preserving the gapless dirac points , while the latter term is similar to the effect of strain @xmath84 , resulting in opening a gap instead .     with a fixed zeeman field , where @xmath85 .
the strength of the strain perturbation is given by @xmath86 ( a ) 0 , ( b ) 0.03 , ( c ) 0.05 , and ( d ) 0.08 . from ( b ) to ( d ) , one of the dirac cones around @xmath4 closes and reopens again due to a band inversion .
the model parameters are the same as those used in fig .
[ fig : b1].,width=264 ]      like a perpendicular zeeman field , the strain field on the surface of a tci may not only give rise to spin reorientation , but it may also induce a charge gap in the dirac cones . however , there are key differences in the effects of strain and zeeman fields . an applied zeeman field , which respects @xmath69 rotational symmetry around the @xmath10-axis , opens gaps with the _ same _ sign in the pair of dirac cones around @xmath3 , while an applied strain field , which breaks @xmath69 rotational symmetry , would induce gaps with _ opposite _ sign .
this observation provides the foundation for tenability of the chern number via the interplay of applied strain and zeeman fields , suggested in ref . using a simplified two - band model for each dirac point .    as a concrete example based on our four - band model , consider an applied zeeman field along the @xmath10 direction along with a strain field which breaks @xmath87 but preserves @xmath59 around @xmath4 on the ( 001 ) surface . from symmetry considerations ,
the strain and zeeman fields around @xmath3 and @xmath4 become coupled as h_|x & + & h_z_03 , + h_|y & + & h_z_03+v_s_23 , where @xmath88 ( @xmath89 ) is the field strength of zeeman ( strain ) field .
note that for a given @xmath88 , all four low - energy dirac points open gaps in the absence of an applied strain [ see fig .
[ fig : biv](a ) ] .
when a non - vanishing @xmath90 is introduced , the gaps of the two massive dirac cones around @xmath4 evolve in an opposite manner in the sense that one cone increases while the other decreases .
this is consistent with a picture in which the strain field induces out - of - the - plane spin components around the two massive dirac cones which form time - reversed partners with respect to their spin textures [ see fig .  [
fig : f1](b ) ] , and therefore , respond oppositely to the existing @xmath88 field .
with increasing strength of @xmath89 , the gap in the dirac cone with a decreasing gap continues to decrease further , becomes gapless at @xmath91 , and reopens inverted again , as shown in figs .
[ fig : biv](b)-(d ) .
this band inversion indicates that the system undergoes a topological phase transition .
the topological nature of a ground state with broken time - reversal symmetry can be characterized by the chern number .
a non - zero value of the chern number , @xmath92 , given by @xmath93}{(e_n - e_m)^2},\end{aligned}\ ] ] where the summation over momentum space is implicit , indicates a finite hall conductance @xmath94 , which can be obtained by integrating the berry curvature of the wave functions of the occupied bands over the momentum space@xcite .
[ fig : chern ] presents the chern number of the system as a function of the zeeman field strength @xmath88 and the strain strength @xmath89 within the framework of our four - band model . in the presence of a zeeman field , by varying the strength of strain from small to large values is seen to drive the system from @xmath95 to @xmath96 , demonstrating the tunability of the system to the strain field . note that in our consideration here we have implicitly assumed that the sample is thick enough so that the hybridization between the top and bottom surfaces of the sample can be neglected ; with @xmath97 the bottom surface will contribute another chern number @xmath95 , which is not shown in fig .
[ fig : chern ] for simplicity .     and
model parameters are fixed as : @xmath98=-2 , @xmath37=-1 , @xmath99=-1 , @xmath100=-1.5 , with other parameters taken to be zero .
the phase boundary ( blue line ) is given by @xmath101 . ]
in discussing electron - electron interactions , we consider short - range repulsive interactions between electrons on the ( 001 ) surface of a tci , and focus on delineating how the phase diagram evolves in the presence of a strain , which breaks one of the mirror symmetries , say , @xmath87 , based on the the interacting hamiltonian : @xmath102 where ^()&=&[c^_|x , p_z(),c^_|x , p_z ( ) , c^_|x , p_x(),c^_|x , p_x ( ) , + & & c^_|y , p_z(),c^_|y , p_z ( ) , c^_|y , p_y(),c^_|y , p_y ( ) ] , [ eq : basis ] with the subscript @xmath103 ( @xmath104 ) denoting the momentum point involved in the @xmath70 expansion . for @xmath53
, we have added the strain perturbation and hence , @xmath105 .
the interaction term , @xmath106 ; @xmath107 and @xmath108 denote repulsive contact interactions between like and unlike orbitals , respectively , with _ u & = & ud^2 r_=p_z , p_x , p_y _ ( ) _ ( ) , + _ v & = & d^2 r_^ _ ( ) _ ^ ( ) , [ hint ] where @xmath109 ( @xmath110 , @xmath111 ) [ @xmath112 are defined in eq .
( [ eq : ftc ] ) below ] .
since we are treating only the low - energy , long - wavelength physics , it is reasonable to consider field operators @xmath113 , which vary slowly on the scale of the lattice constant .
this can be done by fourier transforming the operators : _ , s ( ) = _ e^i c_,s ( ) = e^i c_,s ( ) , [ eq : ftc ] where @xmath114 , expanding from the origin ( 0,0 ) . in terms of these field and density operators , and the identities given in the appendix ,
the full hamiltonian @xmath115 can be rewritten as @xmath116 where @xmath117 has the same form as eq .
( [ eq : basis ] ) except that the @xmath70-dependence is now replaced by @xmath118 .
in order to determine the ground state of the system as a function of the interaction strength and fermion density , we employ the self - consistent mean - field ( mf ) approach , which is expected to be reliable as long as the interaction strength is much smaller than the bulk band gap .
taking @xmath119 as a representative case , we decouple the @xmath120 into bilinear fermion terms . after some straightforward but tedious algebra ,
we obtain 52 order parameters as well as 12 renormalized band parameters as detailed in the appendix , which can be generally expressed as o_=^ _ [ eq : order ] where a summation over momentum and occupied states is implicit , and @xmath121 , with the pauli matrices , @xmath122 , acting on @xmath123 pseudo - spin ( valley ) space .
note that , here we consider only the symmetry breaking orders in the particle - hole channel : superconductivity in the particle - particle channel due to weak onsite repulsion can only be achieved via beyond mf treatment , usually with exponentially small transition temperature@xcite .    in this connection
, we define the `` particle density '' at a given chemical potential @xmath124 as the density deviation from the total particle density where @xmath125 , namely , n=_k
< k_^e()= -_k < k_^e_0()=0 , [ eq : density ] where @xmath126 ( @xmath127 ) is the eigenenergy of the non - interacting ( mf - decoupled interacting ) system ; @xmath128 is a cutoff momentum , which is chosen such that our results are insensitive to its magnitude . note that by any induced valley `` polarization '' , @xmath129 , we mean @xmath130 that is , the `` valley polarization '' due to external strain in the absence of interactions is subtracted in defining @xmath129 .     versus @xmath7 phase diagram without strain .
@xmath131 denotes the particle density at @xmath132 .
solid and dashed lines mark phase boundaries of continuous and abrupt phase transitions , respectively . ]    for a given @xmath6 and @xmath7 , the ground state can then be determined by minimizing the mf free energy with respect to various order parameters using the form detailed in the appendix .
we numerically solve the resulting set of coupled equations self - consistently to obtain the zero - temperature @xmath6-@xmath7 phase diagrams both with and without the applied strain .
note that , in order to make the underlying physics more transparent , we use the following model parameters in computations : @xmath133 , @xmath134 , @xmath135 , @xmath136 , and a cutoff energy , @xmath137 , which determines @xmath128 . however , we expect our results to be generic and relevant more generally , detailed effects of material - specific parameters notwithstanding .
we begin by considering a system without strain , _
i.e. _ , @xmath138 . in the order parameter space we explored in this case ,
the leading orders ( aside from the metal phase ) are basically the ferromagnetic state ( fm ) , @xmath139 , and the @xmath18-orbital rich @xmath140 charge density wave ( @xmath18-cdw ) : = _
( -1)^x+y[n_p_z,()+n_p_z , ( ) ] , where the length scale is in units of the lattice constant .
the zero - temperature @xmath6-@xmath7 phase diagram is shown in fig .
[ fig : pd1 ] .
the @xmath18-cdw phase corresponds to a non - vanishing @xmath141 ( the leading piece ) , and the fm phase represents a non - vanishing value of @xmath139 ( the leading piece ) .
finally , in the metal phase , all symmetries are preserved with vanishing values of all 52 order parameters .    at @xmath142 ,
the `` @xmath18-cdw '' phase appears when the interaction strength @xmath6 is larger than the critical value @xmath143 .
this phase involves a finite value of @xmath141 ; it is associated with broken translation symmetry , while the trs and both the mirror symmetries remain intact [ see schematic fig .  [
fig : realspace](b ) ] , so that each dirac cone remains gapless .
the existence of a critical value of @xmath6 can be understood by noticing that in the non - interacting limit , there are only four dirac points with zero density of states at the fermi level .
the necessity of a threshold value of @xmath6 has also been predicted theoretically in 2d systems with linear energy dispersion , such as graphene and the surfaces of 3d weak tis@xcite .    upon electron doping ,
the system assumes either the fm phase [ see fig .  [
fig : realspace](a ) ] at relatively high particle densities , or the @xmath18-cdw phase [ see fig .  [
fig : realspace](b ) ] at low particle densities .
the fm phase is associated with broken trs as well as both the mirror symmetries simultaneously , resulting in gaps to open up at all the dirac points , although the corresponding ground state is still metallic with finite @xmath124 .
note that the critical @xmath6 , above which the system enters into the gapless @xmath18-cdw phase , is less than the critical value for the fm phase , indicating a favorable free energy gain compared to the competing fm phase at low doping .
-orbital rich , @xmath140 charge density wave state ( @xmath18-cdw ) .
different colors label the two distinct atoms involved ; arrows depict the net out - of - the - plane spin polarization.,width=302 ]    we further note several points in connection with the phase diagram of fig .
[ fig : pd1 ] as follows : ( 1 ) the spin density wave ( sdw ) phase , as anticipated usually in dirac systems like graphene , does not occur here at @xmath142 for large @xmath6 .
this could be attributed to the non - trivial orbital textures and strong spin - orbit coupling in our case ; ( 2 ) we find that the dynamically generated strain - like order , which could gap the spectrum , is relatively disfavored as the system tries to preserve the combined @xmath144 symmetry , although such a phase becomes favored in the presence of external strain ; ( 3 ) if we suppress inter - valley scattering , our preliminary results suggest that the cdw phase is also suppressed@xcite ; and finally , ( 4 ) as we pointed out already in sec .
ii above , the density of states at an energy scale @xmath45 near a van hove singularity diverges like @xmath145 , indicating propensity for the occurrence of more symmetry breaking orders ( including those in the particle - particle channel ) .
the present mean - field treatment would then be inappropriate , even in the weak - coupling limit .
here we consider an applied strain with a moderate coupling strength , @xmath146 [ of order @xmath147 ) ] , which breaks the @xmath13-mirror symmetry .
the resulting zero - temperature phase diagram is shown in fig .  [
fig : pd2 ] .
in contrast to the case without strain , a coexisting phase lying between the @xmath18-cdw and fm phases is now seen to emerge .
suppression of the @xmath18-cdw phase with strain can be anticipated on physical grounds .
the reason is that any @xmath140 order requires good nesting between the valleys , but the strain weakens the nesting tendency by opening gaps at dirac points along the @xmath13-mirror line .
notably , in the strained system , spin polarization in the fm phase is mainly contributed by the @xmath19 orbital due to the broken mirror symmetry , which lifts the @xmath51 orbital to higher energies .     versus @xmath7 phase diagram in the presence of a moderate strain .
lines mark boundaries between various continuous phases . ]
when the coupling strength of the applied strain @xmath89 becomes the largest energy scale in the system , all other degrees of freedom can be ignored at the @xmath104 valley as states around @xmath104 are all gapped out , leading to the phase diagram of fig .
[ fig : pd3 ] .
all @xmath140 orders are seen to disappear , and only the fm phase survives for @xmath148 at a given density @xmath7 .
the presence of the concave boundary region separating the metallic and fm phases in fig .
8 can be understood from energetics : at any given density in this region , the chemical potential lies within the gap induced by the fm order at @xmath149 , yielding a gain in free energy and thus lowers @xmath150 .
as the particle density increases , the chemical potential is eventually unable to locate within the gap and this argument fails .
interestingly , at @xmath142 , although the ground state is insulating with a small gap due to trs breaking , the system does not exhibit quantum anomalous hall effect .
this should be contrasted sharply with the case of a zeeman field applied perpendicular to the ( 001 ) surface of a tci in the non - interacting limit , where such a broken trs leads to non - vanishing chern number as shown in sec .
iiib above .     versus @xmath7 phase diagram in the case of large strain .
solid line marks the boundary between the two continuous phases . ]
some of the strains discussed in this study could be realized experimentally through chemical doping .
since pbse and snse assume cubic and orthorhombic structures@xcite , respectively , strains or lattice distortions could likely be manipulated by adjusting chemical compositions .
in fact , a ferroelectric - like lattice distortion has been observed in the tcis via stm / sts spectroscopy@xcite .
mechanical tuning could provide another approach .
for example , a piezoelectric layer could be grown on the surface of a tci to allow manipulation of strains via electric field as proposed in ref ..    as to practical applications , we have shown clearly that the chern number of the surface states can be tuned through the interplay of applied strain and zeeman fields .
this provides a potential new handle for controlling the topological conducting edge channels for low power - consuming , next generation electronic technologies .
furthermore , our analysis shows that in the presence of electron - electron interactions , ferromagnetism could be turned on / off via strain fields , which , offering a novel pathway toward spintronics applications .    in summary
, we have systematically investigated the effects of strains , which break either one or both mirror symmetries without breaking time - reversal symmetry , on the ( 001 ) surface of the tcis . under mirror - symmetry breaking strains , the low energy dirac cones not only become massive , but also develop hedgehog - like spin textures .
we show that the chern number of the surface states can be tuned via applied strain fields in the presence of a zeeman field .
finally , we delineate the competing orders that can result from varying strengths of short - range ( repulsive ) electron - electron interactions at the mean - field level , and show that the resulting correlated phases are amenable to control through strain fields .
the tunability of the interacting as well as the non - interacting electronic structures and topological states of the tcis via strain fields and their interplay with zeeman fields as revealed by our study suggests new pathways for developing spintronics and other applications platforms based on the tcis .
we thank chen fang and hong yao for useful discussions .
w.f.t . and c.y.h .
acknowledge the support from most in taiwan under grant no.103 - 2112-m-110 - 008-my3 .
the work at northeastern university was supported by the us department of energy ( doe ) , office of science , basic energy sciences grant number de - fg02 - 07er46352 , and benefited from northeastern university s advanced scientific computation center ( ascc ) and the nersc supercomputing center through doe grant number de - ac02 - 05ch11231 . h.l .
acknowledges the singapore national research foundation for the support under nrf award no .
nrf - nrff2013 - 03 .
99 m. z. hasan and c. l. kane , rev .
phys . * 82 * , 3045 ( 2010 ) .
x. l. qi and s. c. zhang , rev .
* 83 * , 1057 ( 2011 ) .
j. e. moore , nature ( london ) * 464 * , 194 ( 2010 ) .
y. ando , j. phys .
. jpn . * 82 * , 102001 ( 2013 ) .
a. p. schnyder , s. ryu , a. furusaki , and a. w. w. ludwig , phys . rev .
b * 78 * , 195125 ( 2008 ) .
a. kitaev , aip conf . proc . * 1134 * , 22 ( 2009 ) . t. l. hughes , e. prodan , and b. a. bernevig , phys . rev .
b * 83 * , 245132 ( 2010 ) .
a. m. turner , y. zhang , r. s. k. mong , and a. vishwanath , phys .
b * 85 * , 165120 ( 2012 ) . c. fang , m. j. gilbert , and b. a. bernevig , phys . rev .
b * 86 * , 115112 ( 2012 ) .
r. slager , a. mesaros , v. juricic , and j. zannen , nature physics * 9 * , 98 ( 2013 ) . k. shiozaki and m. sato , phys .
b * 90 * , 165114 ( 2014 ) .
l. fu , phys .
* 106 * , 106802 ( 2011 ) .
j. c. y. teo , l. fu , and c. l. kane , phys . rev .
b * 78 * , 045426 ( 2008 ) .
k . chiu , h. yao , and s. ryu , phys .
b * 88 * , 075142 ( 2013 ) .
t. h. hsieh , h. lin , j. liu , w. duan , a. bansil , and l. fu , nat .
commun . * 3 * , 982 ( 2012 ) .
xu , c. liu , n.  alidoust , m.  neupane , d.  qian , i.  belopolski , j.  d. denlinger , y.  j. wang , h.  lin , l.  a. wray , g.  landolt , b.  slomski , j.  h. dil , a.  marcinkova , e.  morosan , q.  gibson , r.  sankar , f.  c. chou , r.  j. cava , a.  bansil and m.  z. hasan , nat .
commun . * 3 * , 1192 ( 2012 ) .
p.  dziawa , b.  j. kowalski , k.  dybko , r.  buczko , a.  szczerbakow , m.  szot , e.  usakowska , t.  balasubramanian , b.  m. wojek , m.  h. berntsen , o.  tjernberg and t.  story , nat . mater .
* 11*,1023 ( 2012 )
. y. tanaka , z. ren , t. sato , k. nakayama , s. souma , t. takahashi , k. segawa , and y. ando , nat .
* 8 * , 800 ( 2012 ) .
a. gyenis , i. k. drozdov , s. nadj - perge , o. b. jeong , j. seo , i. pletikosic , t. valla , g. d. gu , and a. yazdani , phys .
b * 88 * , 125414 ( 2013 ) .
a. a. taskin , s. sasaki , k. segawa , and y. ando , phys .
b * 89 * , 121302(r ) ( 2014 ) .
y. okada , m. serbyn , h. lin , d. walkup , w. zhou , c. dhital , m. neupane , s. xu , y. j. wang , r.  sankar , f. chou , a. bansil , m.  z. hasan , s. d.  wilson , l. fu and v. madhavan , science * 341 * , 1496 ( 2013 ) .
i. zeljkovic , y. okada , m. serbyn , r. sankar , d. walkup , w. zhou , j. liu , g. chang , y. j. wang , m. z. hasan , f. chou , h. lin , a. bansil , l. fu , and v. madhavan , nature materials * 14 * , 318 ( 2015 ) .
i. zeljkovic , y. okada , c .- y .
huang , r. sankar , d. walkup , w. zhou , m. serbyn , f. chou , w .- f .
tsai , h. lin , a. bansil , l. fu , m. z. hasan , and v. madhavan , nature physics * 10 * , 572 ( 2014 ) .
m. iizumi , y. hamaguchi , k. komatsubara , and y. kato , j. phys .
jpn . * 38 * , 443 ( 1975 ) . c. fang , m. j. gilbert and b. a. bernevig , phys .
. lett . * 112 * , 046801 ( 2014 ) .
m. serbyn and l. fu , phys .
b * 90 * , 035402 ( 2014 ) . s .-
xu , m. neupane1 , c. liu _
* 8 * , 616 ( 2012 ) .
j. liu , t. h. hsieh , p. wei , w. duan , j. moodera and l. fu , nat .
* 13*,178 ( 2014 ) .
a. h. castro neto , f. guinea , n. m. r. peres , k. s. novoselov , and a. k. geim .
* 81 * , 109 ( 2009 ) .
t. o. wehling , e. sasioglu , c. friedrich , a. i. lichtenstein , m. i. katsnelson , and s. blugel , phys .
lett . * 106 * , 236805 ( 2011 ) .
importance of the coulomb interaction for relativistic quasiparticles may be assessed approximately by the dimensionless interaction strength , @xmath151 , where @xmath152 denotes the dielectric constant , and @xmath153 is the fermi velocity . for graphene , @xmath154 , while for snte @xmath155 is even higher due to a smaller value of @xmath153 .
wang , w .- f .
tsai , h. lin , s .- y .
xu , m.  neupane , m.  z. hasan and a.  bansil , phys .
b * 87 * , 235317 ( 2013 ) .
l. fu , c. l. kane , and e. j. mele , phys .
* 98 * , 106803 ( 2007 ) .
l. fu and c. l. kane , phys .
b * 76 * , 045302 ( 2007 ) .
g. yang , j. liu , l. fu , w. duan , and c. liu , phys .
b * 89 * , 085312 ( 2014 ) .
binghai yan , lukas mu chler , and claudia felser , phys .
* 109 * , 116406 ( 2012 ) .
qi and s .- c .
zhang , physica e : low- dimensional systems and nanostructures * 44 * , 906 ( 2012 ) . c. fang and l. fu , phys .
b * 91 * , 161105 ( 2015 ) .
i. zeljkovic , d. walkup , b. assaf , k. l scipioni , r. sankar , f. chou , v. madhavan , arxiv:1501.01299 . b. andrei bernevig and taylor l. hughes , topological insulators and topological superconductors ( princeton university press , 2013 ) . w. kohn and j.m .
luttinger , phys .
. lett . * 15 * , 524 ( 1965 ) .
s. raghu , x .-
qi , c. honerkamp , and s .- c .
zhang , phys .
lett . * 100 * , 156401 ( 2008 ) .
effects of long - range coulomb interaction should be examined further .
a. j. strauss , phys .
* 157 * , 608 ( 1967 ) .
here we present further details of the mf theory used in sec .
iv . in terms of the slowly varying field operators defined in eq .  ( [ eq : ftc ] ) , the relevant number operators are given by _
p_z(*r*)&=&c_|x , p_z^c_|x , p_z+c_|y , p_z^c_|y , p_z + & + & ( -1)^x+y(c_|x , p_z^c_|y , p_z+c_|y , p_z^c_|x , p_z ) , + _
p_z(*r*)&=&c_|x , p_z^c_|x , p_z+c_|y , p_z^c_|y , p_z + & + & ( -1)^x+y(c_|x , p_z^c_|y , p_z+c_|y , p_z^c_|x , p_z ) , + _ p_x(*r*)&=&c_|x , p_x^(*r*)c_|x , p_x(*r * ) , + _ p_x(*r*)&=&c_|x , p_x^(*r*)c_|x , p_x(*r * ) , + _
p_y(*r*)&=&c_|y , p_y^(*r*)c_|y , p_y(*r * ) , + _
p_y(*r*)&=&c_|y , p_y^(*r*)c_|y , p_y(*r * ) , + _ p_z(*r*)&=&_p_z(*r*)+_p_z(*r * ) , + _ p_x(*r*)&=&_p_x(*r*)+_p_x(*r * ) , + _ p_y(*r*)&=&_p_y(*r*)+_p_y(*r * ) .
we illustrate our mf treatment with the example of @xmath156 . for this purpose
, we decouple bilinear fermion operators composed of @xmath157 or @xmath158 in the particle - hole channel into @xmath159 , where @xmath160 represents fluctuations away from the mean - field values , which are neglected .
the condensation energy , @xmath161 , for @xmath156 is thus given by @xmath162 where the @xmath163-dependence is implicit .
in fact , @xmath164 can be further rewritten as @xmath165 where @xmath166 was defined in eq .
( [ eq : order ] ) and @xmath167 is a unit area . notice that @xmath168 is the particle density defined in eq .
( [ eq : density ] ) . when the similar procedure is applied to other terms in eq .
( [ demo1 ] ) , it turns out that @xmath161 term becomes : @xmath169\notag\\ + s\frac{u}{64}&[(o_{001}-o_{031}+o_{301}-o_{331})^2\notag\\ & + ( o_{002}-o_{032}+o_{302}-o_{332})^2].\end{aligned}\ ] ] this allows us to obtain the mf decoupled @xmath156 as @xmath170\notag\\ -\frac{u}{32}&[(o_{001}-o_{031}+o_{301}-o_{331})\notag\\ & \times\hat{\psi}^{\dag}(\sigma_{001}-\sigma_{031}+\sigma_{301}-\sigma_{331})\hat{\psi}\notag\\ & + ( o_{002}-o_{032}+o_{302}-o_{332})\notag\\ & \times\hat{\psi}^{\dag}(\sigma_{002}-\sigma_{032}+\sigma_{302}-\sigma_{332})\hat{\psi}]\bigg\}\notag\\ -s\frac{u}{64}&[(o_{000}-o_{030}+o_{300}-o_{330})^2\notag\\ & -(o_{003}-o_{033}+o_{303}-o_{333})^2]\notag\\ + s\frac{u}{64}&[(o_{001}-o_{031}+o_{301}-o_{331})^2\notag\\ & + ( o_{002}-o_{032}+o_{302}-o_{332})^2].\end{aligned}\ ] ] by applying the same trick for the other interaction terms in eq .
( [ hint ] ) , we finally find 64 mf undetermined parameters ( summarized in table .
[ t2 ] ) , and obtain the mf free energy per unit area at zero temperature as : f = d^2r^(h_(001)+h_int^mf)-n+e_con , where @xmath171 is the mf decoupled interaction hamiltonian for @xmath172 and @xmath173 is the mf condensation energy .
the ground state can then be determined by solving the set of coupled mean - field equations resulting from minimizing @xmath174 with respect to all possible order parameters .
m2.5cm|m5.5 cm type & mf parameters + renormalized band parameters&@xmath175 , @xmath176,@xmath177 , , this is an order parameter for spontaneous mirror symmetry breaking . ]
@xmath178,@xmath179 , , this is an order parameter for spontaneous mirror symmetry breaking . ]
@xmath180,@xmath181,@xmath182,@xmath183 , , @xmath181 , @xmath182 and @xmath183 is nonzero , this would be an order parameter for spontaneous mirror symmetry breaking . ]
@xmath129,@xmath184 , or @xmath184 is nonzero , this would be an order parameter for spontaneous @xmath49 symmetry breaking . ] @xmath168 + fm & @xmath185 , @xmath186 , @xmath187 , @xmath188 , @xmath139 , @xmath189 , @xmath190 , @xmath191 , @xmath192 , @xmath193 , @xmath194 , @xmath195 , @xmath196 , @xmath197 , @xmath198 , @xmath199 , @xmath200 , @xmath201 , @xmath202 , @xmath203 + cdw & @xmath204 , @xmath205 , @xmath206 , @xmath207 , @xmath208 , @xmath209 , @xmath210 , @xmath211 , @xmath212 , @xmath213 , @xmath214 , @xmath215 , @xmath216 , @xmath217 , @xmath218 , @xmath219 + sdw & @xmath220 , @xmath221 , @xmath222 , @xmath223 , @xmath224 , @xmath225 , @xmath226 , @xmath227 , @xmath228 , @xmath229 , @xmath230 , @xmath231 , @xmath232 , @xmath233 , @xmath234 , @xmath235