Universal edge bands induced by linearly polarized irradiation on phosphorene

Phosphorene (a monolayer of black phosphorus) is a large gap semiconductor with high mobility and has great application potential. Numerical calculations reveal that phosphorene is a topologically trivial material and can only host edge bands on specified edges such as the zigzag edge. A linearly polarized irradiation on the phosphorene lattice results in the dynamic gaps in the quasi-energy spectrum. We found that the irradiation polarized in the zigzag direction induces new edge bands within the dynamic gaps on any type of edge (zigzag, armchair, or other bearded edge). We proposed a new gauge independent quantity, δ + g , to account for the appearence of universal edge bands, where δ is the detune and g is the light induced valence-conduction band transition element. The number of edge bands in the dynamic gaps is reflected by the winding number of it.


Introduction
Periodically driving can turn a normal material into the so-called Floquet topological insulator [1][2][3][4][5]. Light irradiation is an important periodic driving source and its impact on materials draws much attention. Recently, the Floquet topologic phase of condensed matter materials and photonic crystals were observed experimentally [6,7], and verified the validity of the theoretical prediction and showed the potential application of irradiation for topologic phase engineering. Light irradiation affects the bulk band structure of materials by two mechanisms. First, the electron near the neutral point of gapless materials, such as graphene, emits a photon and re-absorbs it to renormalize the band structure, and so a mini-gap is generated at the neutral point [8][9][10][11], which is a second order effect. Secondly, light induces resonant transition between conduction band and valence band states, and results in dynamic gaps in the quasi-energy spectrum at w =  E 2, where ÿ is the reduced Planck constant and ω is the angular frequency of light [12][13][14][15]. The latter is more attractive because it is a first order process and will be considered in this paper. If the light is circularly polarized, it can change the material's topological property, which is reflected by the arising of new helical edge bands in the dynamic gaps [14][15][16][17]. These photon-exited spectrums can be detected by the Angle-resolved photoemission spectra measurement [7]. Because the illuminated system is not an equilibrium one, the occupation rule of the quasi-energy spectrum is complicated, and the edge band induced by irradiation does not carry exactly one conductance quanta in electronic transport [18][19][20][21][22].
Phosphorene, a monolayer of black phosphorus, has great potential for electronic application because its field-effect transistors have a large on-off current ratio at room temperatures [23][24][25]. Its energy gap is moderately large and its anisotropy is quite strong. By applying static strain, phosphorene can be tuned into a band-inversion material [26,27], which could host topological node-line states near the Fermi energy [28]. Vertical electric field can make phosphorene become a Dirac semimetal [29] and induces a topological phase transition in few-layer phosphorene [30]. The linearly polarized laser illumination results in dynamic gaps in the quasi-spectrum of phosphorene, but cannot induce topological phase transition [31]. In this paper, we investigated the band structure of phosphorene ribbons under linearly polarized irradiation. We found that in the dynamic gaps, there exist edge bands on both zigzag and armchair edges if the polarization is along the zigzag direction, but no edge band appears when the polarization is along the armchair direction. Whether or not edge bands exist in the dynamic gaps is governed by a gauge independent quantity, d + g , with δ being the detune and g being the light induced valence-conduction band transition element. For the polarization along other directions, the zigzag component and the armchair component of the polarization vector compete with each other. When the polarization is deviated from the zigzag direction, the k-interval of edge band existence shrinks and it vanishes for a critical polarization orientation.

Hamiltonian
The phosphorene lattice is composed of puckered-up and puckered-down phosphorus atoms as demonstrated in figure 1. The in-plane structure of the phosphorene is a deformed honeycomb lattice with the in-plane parameters a=0.8014, b x =1.515, and b y =1.674 (in units of Å) [32]. The x-axis and y-axis are set to be parallel to the armchair and zigzag directions. Because there is no difference between the on-site potentials on the puckered-up and puckered-down atoms, the phosphorene lattice can be regarded as a bipartite system consisting of A and B sublattices, as shown in the inset. The five types of nonzero hoppings are labeled by t t t t , , , where the summation runs over all nonzero hoppings. In basis of A-B sublattices, the k-space Hamiltonian is and the corresponding eigen states where ñ |c and ñ |v are the conduction and valence band sates respectively and f is the complex angle of h. When a polarized light irradiates normally on the phosphorene sheet, the light can be described by a time-dependent vector potential w A t cos , where A is the amplitude vector and t is the time. The vector potential can be involved into the Hamiltonian by the Peierls substitution (Form here on, ÿ and e in equation expressions are set to be unit to simplify the derivation) This leads to the additional time-dependent Hamiltonian In basis of ñ |c and ñ |v , we apply the rotating wave approximation and have the Hamiltonian The time-dependent Shrödinger equation of the Hamiltonian H t can be numerically solved in the Floquet space in basis of with n being an integer. The quasi-spectrum can be obtained by averaging the density of states over a period of driving. In this paper, we instead use the method described in [31] to calculate the quasispectrum of phosphorene ribbons. The studied armchair and zigzag ribbons lie along x and y directions, respectively, according to the axis definition in figure 1.
show the quasi-spectrum for phosphorene ribbons under the irradiation polarized along xand y-directions 3 . In the band gap of the zigzag ribbon, the nearly flat edge bands survive but are disturbed a little by the irradiation, which is not our target in this paper. (The origination of the edge bands of un-irradiated phosphorene ribbons was discussed in [34] and [35], which can be explained by the winding number analysis of h by choosing proper gauges for the armchair and zigzag edges respectively. The winding number reflects the number of edge bands and was proven to be 0 for armchair ribbons and 1 for zigzag ribbons). Within the conduction or valence band, the irradiation induces the so-called dynamic gap. Interestingly, one can find that new edge bands appear in the dynamic gaps if the polarization is along y-direction (the zigzag direction) regardless of the edge type of the ribbons. This means there must be a new gauge invariance to govern the existence of those edge bands, which will be explored in the following.

Spectrum and edge bands
To find out the physics within the dynamic gaps, we use the time-dependent unitary transformation is the detune, to reduce the time-dependent Hamiltonian H t into a static one, saying we transform the reduced Hamiltonian H rdc into The eigen values of the effective Hamiltonian H eff are . Because H eff is deduced from the original time-dependent Hamiltonian H t without approximation, all the information of the irradiated system is encoded in the effective Hamiltonian. Figures 2(e)-(h) show the dispersion of H eff for the irradiated phosphorene ribbons with zigzag and armchair edges. The features of the quasi-energy bands of H t can be found their reflection in the band structure of H eff : the two dynamic gaps for H t merge in the spectrum of H eff and so do the edge bands.
The winding number of the effective Hamiltonian can be used to handle the edge bands of the ribbons. The complex quantity η has the real part δ and the imaginary part · A v. For a ribbon lying along x direction (armchair ribbon), k x is a good quantum number. At a given quantum number k x , η is a function of parameter k y . When the parameter k y changes though out the Brillouin zone, the endpoint of the complex vector of η on the complex plane moves and evolves as a closed curve. The turn of times of the curve around the origin is called the winding number. Similarly, if the ribbon lies in y-direction (zigzag ribbon), k y is a good quantum number, and trajectory of the endpoint of η is a closed curve when the parameter k x changes through out Brillouin zone. To reflect the number of edge bands, the winding number analysis should use different gauges for different types of edge. We consider the gauge transformation U g , which transforms the operator  H k k to  + ( ) U H U g k k g and translates the state mñ | (m = c v , ) to mñ | U g . According to equation (10), g keeps unchanged under the gauge transformation, and so η is a gauge invariable. This means that the existence of edge bands in the dynamic gaps is independent of the edge type.
Firstly we consider the irradiation is polarized in x-direction (the armchair direction). For this case we have h d = + v i x (A x is set to be a unit which has no influence on the winding number analysis). The trajectory of the endpoint of η for the armchair ribbon at the wavevector denoted by the bold dot in figure 2(e) is shown in figure 2(m). Because v x is even with respect to k x and k y (see the appendix), the endpoint moves forth and back in a round trip and the closed trajectory curve is reduced into a single line. The curve contains nothing, the winding number is 0, and so no edge state can be found in the dynamic gap. The trajectory for the ribbon of zigzag type at the given wavevector labeled by the bold dot in figure 2(g) is illustrated as figure 2(o). The trajectory curve for this case is reduced to a line again, the winding number is zero, and no edge state exists in the dynamic gap. For both armchair and zigzag ribbons, the intercept of the trajectory of η on the imaginary axis means the value of v x at resonance d = 0. Because the trajectory curve modifies when the wavevector (k x or k y , depending on the ribbon orientation) varies, the resonant value of v x is a function of k x (or k y ) [31]. We define the resonant response quantities to understand the property of dynamic gap. Figures 2(i) and (k) show the function of v x 0 versus k x for the armchair ribbon and that versus k y for the zigzag ribbon, respectively. One can see that consistency between the curve shape and the profile of the dynamic gap in the spectrum.
Secondly, we assume that the irradiation polarization is in the y-direction (the zigzag direction), and we have h d = + v i y (A y is set to be a unit). Figures 2(n) and (p) show the winding number pictures at the wavevectors represented respectively by the bold dots in figures 2(f) and (h). The trajectories of η are up-down symmetric closed curves, because v y is a odd function of k x and k y (see the appendix). The origin is inclosed by the trajectory for armchair and zigzag ribbons, and correspondingly, edge states can be found in the band structure at the selected wavevectors. The intercept of the trajectory on the imaginary axis is v y 0 , and a finite v y 0 means the inclusion of the origin. The vertical dashed lines in the figure indicate the interval for v y 0 , also the interval for edge band existence. Figures 2(j) and (l) show the curves of v y 0 as function of wavevector k x for armchair ribbon and k y for zigzag ribbon. The curve for either case hits the axis three times, and correspondingly there are three degenerate points in the spectrum (the degeneracy is not exact because of the lateral confinement effect).
We also check the spectrum of the effective Hamiltonian for phosphorene ribbons with other edges (not shown), such as bearded zigzag and bearded armchair, under linearly polarized illumination, and verify our conclusion: if the polarization is along the zigzag direction, edge bands arise in the dynamic gap, and if along the armchair direction, no edge band turns up.
When the polarization is along neither x-nor y-direction, we have h d y y . The complex number η can be decomposed to d + A v i y y and A v i x x . The trajectory of the former for the armchair or zigzag ribbon is just the curve in figure 2(n) or (p) amplified along the imaginary axis by the factor A y . The distance between the intersect points on the imaginary axis is A v 2 y y 0 . If we pick the latter part up, it lift the amplified curve vertically by the amount A v x x , but the intersect point distance A v 2 y y 0 is not affected by the lift. Figures 3(e) and (f) show the trajectory of η for armchair and zigzag ribbons at the wavevectors denoted by the bold dots. The shape of the trajectory curve of η is deformed and asymmetric, because the lift amount A v x x is not a constant when the wavevector varies. In figure 3(e), the lift amount is large enough so that the origin is exposed outside out the trajectory curve. Therefore, there is no edge state at the selected wavevector. For the case of figure 3(f), the lift amount is not such large, the origin is included by the trajectory, and one can find the edge band states in the spectrum at the corresponding wavevector.
The inclusion of the origin is determined by the competition between A v

Conclusion
Linearly polarized irradiation induces edge bands in phosphorene ribbons on any type of edge in the dynamic gaps if the polarization is along the zigzag direction. The appearing of the universal edge bands is determined by the winding number of d + g . It is gauge independent, and thus is unrelated with the edge type.

Appendix. Symmetry of v
The winding number of η has tight relation with the symmetry of v. According to the amiltonian in equation (2) and the eigen states in equation (6), the quantity á  ñ is calculated as x y x y and correspondingly, we have a r g , A 6 x y x y x y x y Now we consider the symmetry of the x component of v, We proved that v x is even about k x and k y . Turn to the y component, we have So we verified that v y is an odd function of k x and k y .