Tracing the nonequilibrium topological state of Chern insulators
Abstract
Chern insulators exhibit fascinating properties which originate from the topologically nontrivial state characterized by the Chern number. How these properties change if the system is quenched between topologically distinct phases has however not been systematically explored. In this work, we investigate the quench dynamics of the prototypical massive Dirac model for topological insulators in two dimensions. We consider both dissipation-less dynamics and the effect of electron-phonon interactions, and ask how the transient dynamics and nonequilibrium steady states affect simple observables. Specifically, we discuss a time-dependent generalization of the Hall effect and the dichroism of the photoexcitation probability between left and right circularly polarized light. We present optimized schemes based on these observables, which can reveal the evolution of the topological state of the quenched system.
pacs:
I Introduction
Topologically nontrivial phases of matter are a subject of intense current research Hasan and Kane (2010); Moore (2010). They exhibit a range of intriguing and potentially useful properties, such as the quantum anomalous Hall (QAH) effect. The usual notion is that the intrinsic topology of a system cannot be altered by local perturbations, which results in the protection of certain properties due to time-reversal symmetry. In particular, this effect leads to stable surface or edge states and their extraordinary transport properties.
The relation between the topological phase and corresponding observables is still under active investigation. Originally, such a correspondance has been established in non-interacting systems in equilibrium via the Thouless-Kohmoto-Nightingale-Tijs formula Thouless et al. (1982). In this case, the Chern number of each band (labelled by ) fully characterizes the QAH effect: the Hall conductivity amounts to . The bulk-insulating system (characterized by and ) is then called a Chern or quantum Hall insulator (QHI). The Chern number, on the other hand, is obtained via the Berry curvature from the wave-function directly – a quantity that is, strictly speaking, available for non-interacting systems (or within mean-field treatments) only. A possible extension in the context of many-body perturbation theory can be obtained by constructing an effective Hamiltonian involving the self-energy at zero frequency Wang and Zhang (2012), which allows to study the interplay of topological properties and correlation effects in QHIs Budich et al. (2013); Amaricci et al. (2015); Kumar et al. (2016). This approach is based on the connection of the Chern number to the winding number Wang et al. (2010). Alternatively, topological states can be classified by studying the response of a system to external gauge fields Fröhlich and Werner (2013). While these approaches work for noninteracting as well as for interacting electrons, the underlying assumption is that the system is in its ground state. Hence the established concepts are not necessarily applicable to finite temperature or nonequilibrium scenarios, which involve excited states.
This is particularly true for global perturbations such as quenches of the Hamiltonian parameters. For instance, a straightforward definition of a time-dependent Chern number from the time-evolving wave-functions of a non-interacting system will remain constant under unitary evolution. On the other hand, the Hall conductivity can change, e. g. after a quench, which means that it is generally not identical to the Chern number (up to the prefactor ), in contrast to the equilibrium case Caio et al. (2015); Wang et al. (2016); Ünal et al. (2016). Similarly, the bulk-boundary-correspondance might be lost after a quench Schmitt and Wang (2017). It is thus a relevant task to identify experimentally accessible quantities which allow to trace the nonequilibrium evolution of topologically nontrivial systems.
In this work, we investigate different schemes that enable us to study the nonequilibrium and transient dynamics of QHIs. We focus on (i) the time-resolved Hall effect using appropriately shaped electromagnetic pulses, and (ii) the photoabsorption asymmetry with respect to left/right circularly polarized light – a novel approach which has been suggested recently Tran et al. (2017). Both methods are based on directly observable quantities and are thus well suited for the study of (effectively) noninteracting as well as correlated and/or dissipative systems. We demonstrate their applicability by considering the well-known massive Dirac model (MDM) on a square lattice Qi et al. (2008), which captures Liu et al. (2008) the topoligical phase transition in HgTe quantum wells Bernevig et al. (2006) as a generic example. We focus on quench dynamics and demonstrate how a transition between phases of distinct topological character manifests itself in these observables. Furthermore, we study the influence of dissipation due to electron-phonon (el-ph) coupling on the transient dynamics to demonstrate the robustness of the proposed schemes.
Ii Model
As a paradigm model for two-dimensional systems we consider the massive Dirac model on a square lattice. The electronic Hamiltonian reads
(1) |
where the -dependent single-particle Hamiltonian has the generic form . Here, the denote the pseudo-spin operators with respect to the underlying bands. The coefficients are defined by
(2) |
The eigenstates of the single-particle Hamiltonian are denoted by .
Note that we limit ourselves to a spin-restricted model here, as the Hamiltonian (1) is spin-independent. As shown in Ref. Liu et al., 2008, in HgTe quantum wells – which are well modelled by Eqs. (1) and (II) – doping with manganese allows to shift the spin up/down bands in such a way that only one spin channel remains important. Here we assume such a situation and thus focus on the charge QAH effect instead of the usual quantum spin Hall effect (QSH).
It is straightforward to see that diagonalizing the Hamiltonian (1) gives rise to a trivial band insulator (BI) for , while the system corresponds to a topological insulator (TI) ^{1}^{1}1Even though the term topological insulator (TI) refers to a more general concept than the QAH insulator, we use the abbreviation TI throughout the text. for . As usual, the topological character of the bands can be determined by the Chern number , which is defined (for each band , respectively) by the integral over the Berry curvature
(3) |
Here, is the Berry connection corresponding to the upper () or lower () band.
ii.1 Quench dynamics
A quench from the BI into the TI phase (or vice-versa) offers insight into the underlying topological properties. For instance, the insulating state with nontrivial Chern number cannot be altered under unitary time evolution (which preserves time-reversal symmetry) – hence, the question arises to which state the system is driven. Furthermore, quenches may induce dynamical phase transitions, which have recently been studied for the MDM Heyl and Budich (2017). Quenches can be realized, for instance, by photodoping pulses in interacting electronic systems, leading to transient band shifts Eckstein and Werner (2013); Golež et al. (2017), or by changing the strength of the periodic driving in Floquet topological insulators Lindner et al. (2011); D’Alessio and Rigol (2015).
Here we study quenches of the gap parameter between the values and . The band hybridization is fixed at in what follows. The transition from phase A () to B () is triggered by a softened ramp of the form
(4) |
with , defined by the quench time and the duration ().
In the absence of any further interactions, the dynamics can be captured by solving the time-dependent Schrödinger equation
(5) |
where we have assumed half filling. The time-dependent single-particle Hamiltonian is defined in analogy to Eq. (4). As mentioned above, a straightforward generalization of the definition of the Chern number (3) as
(6) |
with the time-dependent Berry connection will be invariant under any unitary time evolution Caio et al. (2015); Schmitt and Wang (2017). As discussed by Wang et al. (ref. Wang et al., 2016), the Hall conductance may however undergo a change. Assuming that decoherence effects have suppressed the off-diagonal elements of the density matrix after sufficiently long time, the steady-state Hall conductance can be defined as
(7) |
Here, stands for the Berry connection of the post-quench Hamiltonian , whereas is the occupation with respect to the post-quench band structure. For the two-band MDM, is easily expressed in terms of the overlap of the pre- and post-quench basis and does not depend on the quench details.
ii.2 Electron-phonon coupling
The nonequilibrium Hall conductance (7) assumes the system to have lost the coherences due to environmental coupling; otherwise, the coherent oscillations after excitations induced by the quench hamper a unambiguous definition of . The most important intrinsic source for such dephasing effects in real systems is (besides structural defects) el-ph coupling. Since the coupling to the lattice vibrations entails – besides the decoherence effects – dissipative population dynamics, the steady state reached by the system is – in general – different from the above quench scenario Wolff et al. (2016). Therefore, we treat the el-ph coupling explicitly by extending the Hamiltonian to
(8) |
For the interaction term, we use the Fröhlich coupling van Leeuwen (2004) in two dimensions:
(9) |
Here, is a constant determining the overall coupling strength and denotes the number of points. The phonon modes with momentum are represented by the corresponding mode vector and the coordinate (momentum) operators (). The latter define the phonon Hamiltonian by
(10) |
To compute the phonon modes in the square lattice, we assume a diatomic basis characterized by masses and three force constants describing the nearest neighbour, next-nearest neighbour and diagonal interaction, respectively. Diagonalizing the dynamical matrix yields two acoustic and two optical (longitudinal and transverse) phonon modes. The parameters are chosen to produce a similar phonon dispersion as is known for HgTe Ouyang and Hu (2015) – up to an overall scaling constant. We have slightly increased the phonon energies to make the effects due to the el-ph coupling more visible. However, it is important to note that we stay in the realistic parameter regime where only intraband transitions can be induced by scattering from phonons. The phonon dispersion and the corresponding density of states (DOS) is shown in Fig. 1. In what follows, we assume low temperatures by fixing the inverse temperature at , such that only the acoustic phonons around the point are thermally activated.
It should be mentioned that the role of the el-ph coupling in topological insulators, in particular at the surface, is a topic of recent discussions. While some works point out the strong influence of el-ph coupling for inelastic scattering Giraud and Egger (2011); Giraud et al. (2012); Costache et al. (2014) and for intraband relaxation of photoexcited systems Hatch et al. (2011), other measurements suggest weak el-ph interaction effects Pan et al. (2012). Here, we take a different angle and treat the coupling strength as a parameter.
ii.3 Equations of motion in the presence of Electron-phonon interactions
The full Hamiltonian (8) constitutes an interacting electron-boson model. The numerically exact solution is out of reach for the typical number of points in reciprocal space. Furthermore, since the phonon energies are much smaller than the electronic energy scale, a weak-coupling treatment is suitable. In this context, the nonequilibrium Green’s function (NEGF) approach in its time-dependent formulation has become an important tool recently Sentef et al. (2013); Kemper et al. (2013); Murakami et al. (2015); Säkkinen et al. (2015); Murakami et al. (2016); Schüler et al. (2016); Tuovinen et al. (2016). However, most of the approaches resort to a local approximation to the el-ph interaction. This approximation excludes intraband transitions, which are, as discussed above, the only available relaxation channel in our case. Therefore, a momentum-dependent treatment of the phonons is required, which increases the computational costs of the NEGF method considerably. Since neither renormalization effects of the electronic structure due to the coupling to the phonons nor the backaction on the phononic degrees of freedom is of particular significance for the relaxation in the MDM, one can employ a simplified theory based on a master-equation approach within the Markovian approximation Breuer and Petruccione (2002). The usual Lindblad equation, however, is formulated in terms of the many-body density matrix Usenko et al. (2016). For calculating the time evolution of the single-particle density matrix, , additional conditions are required, reflecting not only the overall particle conservation, but the fermionic statistics, as well. This can be achieved by extending the linear Lindblad equation to a nonlinear master equation Rosati et al. (2014, 2015). Alternatively, equivalent master equations can also be obtained from with NEGF formulation within the generalized Kadanoff-Baym ansatz Schlünzen and Bonitz (2016) and applying the Markovian approximation Langreth and Nordlander (1991).
Within the master-equation approach, the equation of motion (EOM) for the single-particle density matrix reads
(11) |
where we have employed a more compact matrix notation. Besides the unitary time evolution captured by the first term, the el-ph coupling is described by the scattering term Rosati et al. (2014)
(12) |
Here, is the density matrix of the hole states, while denotes the scattering rate which determines the transition probability. As the scattering events are limited to intraband transitions governed by momentum (and energy) conservation, the scattering rates simplify to Rosati et al. (2015)
(13) |
with
(14) |
Here, the spectral function of the unoccupied phonon modes (greater Green’s function in the NEGF context) is defined by
(15) |
where denotes the Bose distribution. In practice, the Dirac- functions in Eq. (15) are replaced by Gaussians with a small broadening parameter .
Iii Observables in equilibrium
Let us now proceed by defining the observables which will be used to trace the nonequilibrium dynamics. We first consider the equilibrium case, where the system is either in the band insulating or QHI phase and then discuss how to extend the schemes to the time-dependent case.
iii.1 Circular asymmetry
As suggested by Tran et al. (Ref. Tran et al., 2017), the depletion rate () upon irradiation of left (right) circularly polarized light with frequency yields a direct measure of the Chern number of band . In particular, the frequency-integrated asymmetry signal is proportional to . This property was derived for a noninteracting system; however, it is generic and can also be exploited in interacting systems. A closely related effect is the pronounced polarization dependence observed in angle-resolved photoemission (ARPES) from graphene, mapping out the Berry phase of the individual bands Hwang et al. (2011); Liu et al. (2011). Alternative ways of extracting the topological character of the system in experiments is, besides the aforementioned Hall effect, are Aharonov-Bohm-type interferometry Duca et al. (2015) in optical traps and spin-polarized measurements Wu et al. (2016). However, here we focus on observables that can be easily extended to the transient regime.
For the two-band MDM, Fermi’s golden rule yields
(16) |
Here, the transition operator, derived from the standard Peierls substitution, reads
(17) |
Provided the transition is permitted by energy conservation, the excitation probability is determined by the matrix elements . The asymmetry is presented for the BI and the TI in Fig. 2 (a) and (b), respectively. We neglect the el-ph coupling at this point ().
As can be seen in Fig. 2, the main difference between BI and TI is the strong negative asymmetry in the vicinity of the point in the latter case (panel (b)), which is most pronounced at the -points where the avoided crossing between the two bands occurs. This can be understood from the fact that is proportional to the Berry curvature in the two-band case Tran et al. (2017).
In view of transient measurements, it would be most useful to extract the topological character of the system by applying short (circularly polarized) pulses rather than the continuous-wave (plus integrating over frequencies) approach presented in Ref. Tran et al., 2017. Figure 2 suggests to consider transitions close to the point, where the difference in the Berry curvature between the BI and TI is most pronounced. The thus required spectral resolution has to be balanced against time resolution determined by the pulse duration. We have chosen electric field pulses of the form
(18) |
for . The left or right polarization vector (superscript or , respectively) is given by . In order to determine the corresponding asymmetry signal, we propagate the time-dependent Schrödinger equation in the presence of the electric field (18) by the Peierls substitution ( is the vector potential corresponding to the field) in the weak-field limit (). The expected left-right asymmetry of the depletion rate translates into an asymmetry of the photoexcitation probability. After some tests we found that pulses as short as offer a good comprise between sharpness in frequency space and pulse duration. The excitation probability for both the BI and the TI case is depicted in Fig. 3 as a function of the central frequency .
In line with Ref. Tran et al., 2017, integrating over all frequencies yields zero for the BI (since ), while a nonzero value is obtained for the TI. As expected from the behavior of the matrix elements (see Fig. 2), the region in reciprocal space where the Berry curvature is the strongest in the TI case is particularly suited for mapping out the topological character of the system. This leads to the optimal frequency where the BI predominantly absorbs left circularly polarized light (corresponding to the red region around the origin in Fig. 2(a)), but where the Berry curvature leads to a strongly enhanced absorption of right circularly polarized radiation in the TI case. Choosing we obtain a field pulse which is ideally suited for tracing the transient dynamics of the system upon photoexcitation or after a quench. Note that the absorbance of the BI is reduced compared to the TI due to the larger band gap.
The next important question to address is if a similar behavior can be expected in the presence of el-ph coupling. To address this issue, we solved the master equation (11) including the electromagnetic field and computed, as for the non-interacting case, the photoexcitation probability. The result for several moderate coupling strengths is presented in Fig. 4. Besides an overall suppression of the absorbance and the less pronounced difference between the the excitation probabilities for left/right circularly polarized light, the qualitative behavior is consistent with the dissipation-less case (Fig. 3). The visible difference of the asymmetry for the relatively strong el-ph couplings demonstrated in Fig. 4 implies that one can obtain valuable information on the topological character even for weak to moderate strength of dissipative effects.
iii.2 Time-resolved Hall effect
The emergence of the integer Hall effect provides direct access to the topological character in equilibrium. A possible extension of this concept to a time-dependent scenario is – similarly as in subsection III.1 – to apply suitably shaped pulses with parameters optimized in the equilibrium case. To this end, we computed the optical conductivity for both the BI and the TI phase. The real part is shown in Fig. 5.
The plateau in at small frequencies gives us some guidance in choosing the spectral features of a suitable probe pulse which allows to map out the topological character: (i) the pulse needs to be short to enable us to trace the transient dynamics, (ii) the pulse in frequency space needs to have a maximal overlap with the region , and (iii) is required within the dipole approximation. Electromagnetic field pulses which optimally fulfil the criteria (i)–(iii) are half-cycle pulses (HCPs) Moskalenko et al. (2017). HCPs are pulses with a dominant, short peak and a weak and long tail (which hardly influences the dynamics). The dominant peak makes the field effectively unipolar, as the spectral weight is maximal in the vicinity of . Here we employ the parameterization with
(19) |
The parameters are the effective pulse duration and the shape parameter , which we fix at in accordance Moskalenko et al. (2006) with typical pulses generated in experiments You et al. (1993); Jones et al. (1993). The analytical expression in Eq. (19) fulfils exactly. The corresponding power spectrum in frequency space is shown in Fig. 5. We have found to be a good compromise between a short pulse duration and maximal overlap with the plateau region of the TI, as can be seen in Fig. 5.
The properties of the HCPs translate – in combination with the frequency dependence of the optical conductivity – into a distinct behavior of the time-dependent current
(20) |
Here, is the velocity matrix. Within the weak-field regime, linear response theory applies and relates the current to the driving field and the optical conductivity:
(21) |
Assuming a linearly polarized HCP in the -direction, the linear response relation (21) yields for the current in the -direction in frequency space
(22) |
which implies for the time-dependent current
(23) |
Thus, the current orthogonal to the field polarization is expected to closely resemble the driving pulse in the TI case, while the current will almost vanish for the BI. This behavior is confirmed by the numerical solution of the Schrödinger equation in the presence of the HCP (19) and the resulting time-dependent current , displayed in Fig. 6 for the TI case. The Hall current has a strong peak for times where the HCP has its maximum, while the current shows an oscillatory behavior. As expected, the simulations show that the Hall current is negligible for the BI phase (not shown).
Since detecting a time-dependent current on the typical time scale of the pulse (which is in the femto- to picosecond range) is difficult in experiments, we propose to analyze the behavior of the pulse-current correlation function
(24) |
This signal could be detected similarly to the total induced charge, but weighted with the known driving pulse. The behavior of the Hall current observed in Fig. 6 can thus be characterized by a peak at , while a transient current response originating from a pronounced variation with respect to the frequency will not possess this feature. This is confirmed in Fig. 6.
Iv Tracing the quench dynamics – unitary time evolution
Let us now investigate the dynamics of the system after a quench across the phase boundary and the manifestation of this transition in the observables discussed in subsection III. We focus on the noninteracting case first.
Assuming that the system is initially in equilibrium () with the lower band completely filled and the upper band empty, the gap parameter is ramped up or down in a time (see subsection II.1). This short ramp time corresponds to an almost ideal quench, i. e. the occupations of the post-quench bands is given by
(25) |
The post-quench occupation (25) is shown along the standard path in the Brillouin zone of the square lattice in Fig. 7, marked by points where . In fact, is close to one for most points.
iv.1 Time-dependent Hall effect
Examining the structure of the density matrix after the quench, one realizes that coherent superpositions of the two bands play a major role, such that the system is far away from a steady state for which the Hall conductance (7) can be defined. However, to have a measure of the post-quench state the system is driven to, we can investigate the response to HCPs as discussed in subsection III.2. The current induced by a HCP polarized along the direction after performing the quench is shown in Fig. 8(a). In comparison to the current response in equilibrium (Fig. 6), the coherent oscillations of the current dominate the hump at times when the electric field of the pulse is strong. The magnitude of these oscillations is considerably larger than in the equilibrium case, showing that the superposition state after the quench is quite different from the equilibrium TI state.
At first glance, is seems that the current in the direction is not displaying the behavior discussed in subsection III.2, as the magnitude of the current is quite small even when the electric field reaches its maximum. Nevertheless, the pulse-current correlation function defined in Eq. (24) and presented in Fig. 8(b), exhibits the distinct feature of the QHI: possesses a clear maximum at , which indicates a plateau behavior of the optical conductivity at and thus the presence of a nonzero Hall effect. The magnitude of is, however, significantly (approximately by a factor of four) reduced with respect to the equilibrium case (Fig. 6(b)). Hence one would expect a static Hall conductance of about .
We also performed analogous simulations for the quench and found the emergence of a very small, but finite Hall effect after the quench.
The coherent superposition present in the post-quench state results in an internal dynamics which might interfere with transient measurements. However, one can expect decoherence effects to diminish these coherences, leading to a mixed steady state. In this case, the static Hall effect yields valuable information on the post-quench state, as discussed in the next subsection.
iv.2 Steady-state conductance
As can be inferred from Fig. 7, the occupation after the quench reflects the pre-quench situation. For instance, the complete filling of the lower BI band is preserved when switching to the TI, apart from the avoided crossing points. This illustrates the conservation of the topological character of the system as the occupation of bands with the same orbital character (which is interchanged at the crossing points in the TI band structure, see Fig. 7) remains constant. Nevertheless, the Hall effect deviates from the equilibrium behavior, which can be seen by evaluating the steady-state optical conductivity in accordance with Ref. Wang et al., 2016 by
(26) |
Here, denotes the momentum-resolved current operator. The resulting conductivity is shown in Fig. 9. Note that computing the optical conductivity analogously to Eq. (7) assumes that all off-diagonal elements of the density matrix, which capture the coherent oscillations of the system after the excitation, are zero.
As can be inferred from Fig. 9, the Hall conductance significantly deviates from the equilibrium value. In the quench, the system acquires a finite Hall conductance , while the quench leads to . These values are consistent with the time-dependent response discussed in subsection IV.1. Without the unit factor , the latter can be regarded as a nonequilibrium generalization of the Chern number Wang et al. (2016). Note that the concrete numbers depend on both the pre- and the post-quench gap parameter .
iv.3 Transient circular asymmetry
Let us now proceed to transient properties. As discussed in subsection III.1, the circular asymmetry is a very promising candidate for tracing the dynamics in real time. To find a suitable analogue to the equilibrium scenario, we performed test calculations of the population dynamics driven by circularly polarized pulses (with the same parameters as in subsection III.1 and with central frequency ) after the system has been quenched. For the case one finds that the system is preferably excited by right circular pulses, while a left circular pulse results in a weaker depletion of the post-quench lower band. In contrast, circularly polarized pulses applied to the BI after the quench from the TI result in a depletion of the upper band instead (irrespective of the polarization). Furthermore, one observes an oscillatory dependence of the photoexcitation probability after a quench, which originates – analogously to the current discussed in subsection IV.1 – from the coherent superposition of the states belonging to the upper and lower band, respectively.
For these reasons, we propose to utilize the absorbed energy of the left or right circular probe pulses as a footprint of the topological character in nonequilibrium instead. Importantly, the energy absorption can be measured in experiments directly by placing photon detectors behind the sample. It should be noted that can also be negative when the system is excited, which corresponds to stimulated emission rather than absorption. This difference to the equilibrium case needs to be taken into account. We thus define a transient generalization of the circular asymmetry by
(27) |
where is the energy of the probe pulse which is absorbed (or emitted) by the system. The time-resolved circular asymmetry signal (27) for both quench scenarios is presented in Fig. 10 as a function of the delay between the starting time of the pulse () and the time when the system is quenched (). The absorbed energy is computed by taking the energy difference , where denotes the total energy of the quenched system in absence of the probe pulse, while is the total energy of the quenched and probed target. The energy is measured at a sufficiently large reference time after the quench and the pulse.
Figure 10 demonstrates that the quench dynamics can be traced in the time domain by the circular asymmetry. For the system with before switching, the asymmetry is positive (see Fig. 3). Hence, for one observes in Fig. 10 for the case (purple curve). For , the asymmetry assumes negative values (with larger modulus, as well) which – in accordance to Fig. 3 – indicates the transition to the TI. The opposite behavior can be observed when switching (blue curve). Further features are weakly pronounced coherent oscillations of the asymmetry after the quench to the TI which originate from the off-diagonal elements of the time-dependent density matrix. The corresponding time scale is determined by the energy difference between the states whose occupation is changed by the quench.
As our time-dependent simulations demonstrate, the time-resolved quench-probe asymmetry signal based on the absorbed energy provides a robust tool to trace the transient dynamics of the system after a quench. It primarily maps out the circular asymmetry of the underlying bands and is less sensitive to the nonequilibrium occupation. These features clearly distinguish the time-resolved asymmetry from the nonequilibrium Hall effect discussed above and render it a powerful complementary tool.
V Tracing the quench dynamics – dissipative time evolution
After having analyzed the unitary dynamics of the system after a quench, we now investigate how the picture changes if el-ph interactions, as discussed in subsection II.2, are present. Generally, the effect of coupling to the phonon modes is expected to give rise to dissipative dynamics, lowering the energy after the quench excitation. Revisiting Fig. 7 one can expect a qualitatively different behavior for the quench as compared to the case . If the system is quenched from the TI to the BI, the occupation in the upper band (see Fig. 7(a)) is located around the energy minimum at the point. Hence, no energy can be extracted from the system after the quench. The effect of the el-ph coupling is in this case primarily the dephasing of the coherences induced by the quench.
v.1 Transient dynamics probed by time-resolved photoemission
We performed numerical simulations of the quench dynamics by solving the master EOM (11) as described in subsection II.3. The parameters are – apart from the el-ph interaction – the same as in section IV. To understand the time evolution of the band structure and the nonequilibrium occupation, the most convenient quantity to look at is the time-dependent occupation with respect to the post-quench basis and – as complementary information – the transient photoelectron spectrum. Time-resolved ARPES (tr-ARPES) has recently become a standard tool for tracing the time evolution in correlated systems Schmitt et al. (2008); Graf et al. (2011); Bovensiepen and Kirchmann (2012); Smallwood et al. (2012), in parallel with a rapid development of state-of-the-art theoretical descriptions within the NEGF framework Sentef et al. (2013); Kemper et al. (2013); Murakami et al. (2016). Modeling the photoemission process by a pump pulse (frequency , pulse envelope ) yields the photocurrent Kemper et al. (2013)
(28) |
Here, denotes the lesser Green’s function, which we express within the generalized Kadanoff-Baym ansatz as with the time-evolution operator corresponding to the single-particle Hamiltonian . Varying the delay between the excitation (the quench at in our case) and the time when the probe pulse is applied, Eq. (V.1) provides a - and energy-resolved pump-probe spectrum.
In Fig. 11 we present the time-dependent occupation (upper row) along with the corresponding tr-ARPES spectra (lower row) as a function time delay for the quench . In the calculations of the tr-ARPES signals, we used the same parameters for the pulse as in section IV. The short pulse length gives rise to the broadening of the spectra in Fig. 11. At , the system is still in equilibrium and the ARPES signal follows the pre-quench band structure (orange dashed line). The post-quench basis is the BI – therefore, the occupation with respect to the BI bands exhibits a hole around the point up to the region where the avoided crossing occurs. This demonstrates a different aspect of the topological state: the occupation of the first Brillouin zone with respect to the dominant orbital character of the lower band is not singly-connected. Following the time evolution right after the quench () we see a similar picture as in Fig. 7: the ARPES spectrum now follows the post-quench band structure (solid purple lines), while the upper band is populated around the point, whereas the lower band is empty in this region. The population assumes values between zero and one around the crossing region, such that a slight lowering of the total energy by el-ph relaxation is possible. This effect can be observed for later times (). The steady state () however shows only a slight blurring of the occupation compared to directly after the quench.
The nonequilibrium dynamics is considerably more pronounced in the quench scenario , analyzed again in terms of the time-dependent population of the post-quench lower band in the tr-ARPES spectra in Fig. 12.
Right after the quench (), the ARPES spectrum closely resembles the BI band structure, apart from a shift to larger energies. However, as the occupation of the lower band shows, the post-quench TI band is empty between the point and the crossing region, while the upper TI band is populated in this region in the Brillouin zone. It is clear from Fig. 12 that this nonequilibrium occupation does not correspond to an energy minimum as filling the occupation hole in the lower band and a relaxation towards the energy minimum at the crossing points in the upper band result in a lowering of the total electronic energy. These dissipation processes are efficiently mediated by the el-ph interaction as the time-dependent occupation and the ARPES spectra for later times () demonstrate. The steady state reached at has a peculiar configuration: the occupation hole around the point has been filled, while the population of the upper band has relaxed down to the crossing point. Interestingly, from the ARPES spectrum alone one could suspect that the system has fully relaxed to a TI. However, the nonequilibrium occupation which involves both bands gives rise to a steady-state optical conductivity (Fig. 13) which deviates considerably from the equilibrium behavior (Fig. 5) in its strongly suppressed conductance and Hall conductance. The direct current Hall conductance is reduced to for the quench to the BI, while we find in the scenario. Note that the value of the Hall conductance in the TI final state is considerably smaller than in the quench scenario without el-ph interaction (Fig. 9). This can be explained by the distinct occupation in the post-quench steady-state in the presence of the el-ph coupling. The integer Hall conductance originates from interband transitions in the crossing region of the upper and lower band. Exactly those transitions are strongly suppressed in the steady-state of Fig. 12, as the occupation in the lower band in the crossing region is depleted, while the available states in the upper band are mostly occupied.
v.2 Transient circular asymmetry
Let us now investigate if the transition from the topologically trivial BI to the QHI or vice-versa can be traced in the time domain by circularly polarized pulses in an analogous fashion as for the dissipation-less case discussed in subsection IV.3. We employ the same scheme: a left or right circularly polarized pulse photoexcites the system before, during or after the quench. The energy absorbed by the pulse and, in particular, the left-right asymmetry should then reveal the character of the steady state or transient state. We have applied this recipe, using the definition (27), and present the respective energy absorption asymmetry in Fig. 14 for two representative cases of the el-ph interaction: (i) weak coupling (Fig. 14(a)–(b)) and moderate coupling strength (Fig. 14(c)–(d)). We represent the absorption asymmetry as a function of the quench-pulse delay as in Fig. 10.
For , Fig. 14(a) shows a transition similarly to the unitary case. The asymmetry changes from positive values (with small magnitude) to negative, indicating the switch from BI to TI. On the other hand, the quench (Fig. 14(b)) is accompanied by a switch of from negative to positive values (with strongly reduced magnitude). Generally, the scheme for tracing the quench dynamics works very well for the dissipative case with weak el-ph interaction.
Turning to stronger el-ph couplings, the picture slightly changes. For the quench (Fig. 14(c)), one can observe the asymmetry switching from small positive values to the negative region across the quench. The almost linear decrease of towards larger is attributed to the el-ph induced relaxation after the pulse (see Fig. 12). One can readily check that the photoexcitation probability due to the specific pulse is larger in the post-quench relaxed steady state than right after the quench, because more states are occupied in the energy window defined by the pulse frequency. Therefore, the absorption of the left-circular probe pulse becomes more efficients as increases, until the steady state is reached (at ). Generally, the dependence of on is much more pronounced in the TI phase, since the influence of the el-ph coupling (as discussed in subsection V.1) is stronger.
Vi Conclusions
We have studied the quench dynamics of the MDM as a generic model for two-dimensional topological insulators with a special emphasis on how the nonequilibrium and transient properties are reflected in experimentally accessible quantities. We have focused on two promising observables which reveal the topological character of the system: the steady-state and time-dependent Hall effect and the asymmetry of photoexcitation with respect to left or right circularly polarized pulses. Based on a realistic model for two-dimensional topological insulators we have defined suitable probe-pulse shapes by considering the equilibrium model, both for the dissipation-less case and including electron-phonon interactions. We then applied these optimized pulses to trace the nonequilibrium dynamics after a quench. Both the time-dependent Hall effect and the circular dichroism of the absorbed energy provide valuable information on the system. While the Hall current can deliver important insights into the steady state, it turns out to be less suitable for the analysis of transient states. This is due to the coherent superpositions in the system, which give rise to an intrinsic transient dynamics. The circular asymmetry of the absorbed energy, on the other hand, is much less sensitive to these effects, since it is based on the occupation of the bands only. The latter approach is thus particularly well suited for the study of the switching process between different phases.
In the presence of electron-phonon coupling, the quench dynamics can significantly differ from the dissipation-less case, provided the energy of the system is reduced by scattering from phonons. We analyzed these effects in terms of the time-resolved ARPES spectra. We investigated the steady-state properties by the nonequilibrium Hall conductance and, finally, analyzed the quench dynamics in terms of the circular asymmetry. As an important point for the potential experimental realization, we have demonstrated the robustness of our prosed transient measurement in the presence of weak to moderate electron-phonon coupling. The transient circular dichroism has thus been shown as a very promising tool to obtain insights in the up to now little explored field of nontrivial topological phases in nonequilibrium.
Appendix A Numerical solution of the master equation
For a stable numerical solution of a master equation of the type of Eq. (11) we adapt a method for propagating NEGFs in time Stan et al. (2009). The interval is discretized into an equidistant grid . In order to perform the step , we separate the unitary time evolution from the scattering term by the ansatz
(29) |
Here, denotes the time-evolution operator, which we approximate by . Inserting Eq. (A) into the EOM (11) then yields
(30) |
Apart from the approximation to the time-evolution operator , Eq. (A) is still exact. A simple and numerically stable propagation scheme is obtained by approximating . Using the Baker-Hausdorff formula, the time step (A) can be expressed as
(31) |
with and . In practice, we truncate Eq. (A) after the fourth order (). The half-step scattering term is obtained by fourth-order polynomial interpolation using , . This requires knowing the scattering term at the next time step, , which is a function of the yet unknown density matrix . Therefore, we employ a predictor-corrector scheme where is first estimated by third-order polynomial extrapolation, allowing to compute and thus obtain . The latter two steps are then iterated at each time step until is converged.
Acknowledgements.
The calculations have been performed on the Beo04 cluster at the University of Fribourg. This work has been supported by the Swiss National Science Foundation through NCCR MARVEL and ERC Consolidator Grant No. 724103. We thank Markus Schmitt and Stefan Kehrein for fruitful discussions.References
- Hasan and Kane (2010) M. Z. Hasan and C. L. Kane, “Colloquium: Topological insulators,” Rev. Mod. Phys. 82, 3045–3067 (2010).
- Moore (2010) J. E. Moore, “The birth of topological insulators,” Nature 464, 194–198 (2010).
- Thouless et al. (1982) D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, “Quantized Hall Conductance in a Two-Dimensional Periodic Potential,” Phys. Rev. Lett. 49, 405–408 (1982).
- Wang and Zhang (2012) Z. Wang and S.-C. Zhang, ‘‘Simplified Topological Invariants for Interacting Insulators,” Physical Review X 2, 031008 (2012).
- Budich et al. (2013) J. C. Budich, B. Trauzettel, and G. Sangiovanni, “Fluctuation-driven topological Hund insulators,” Phys. Rev. B 87, 235104 (2013).
- Amaricci et al. (2015) A. Amaricci, J. C. Budich, M. Capone, B. Trauzettel, and G. Sangiovanni, “First-Order Character and Observable Signatures of Topological Quantum Phase Transitions,” Phys. Rev. Lett. 114, 185701 (2015).
- Kumar et al. (2016) P. Kumar, T. Mertz, and W. Hofstetter, “Interaction-induced topological and magnetic phases in the Hofstadter-Hubbard model,” Phys. Rev. B 94, 115161 (2016).
- Wang et al. (2010) Z. Wang, X.-L. Qi, and S.-C. Zhang, “Topological Order Parameters for Interacting Topological Insulators,” Phys. Rev. Lett. 105, 256803 (2010).
- Fröhlich and Werner (2013) J. Fröhlich and P. Werner, “Gauge theory of topological phases of matter,” EPL (Europhysics Letters) 101, 47007 (2013).
- Caio et al. (2015) M. D. Caio, N. R. Cooper, and M. J. Bhaseen, “Quantum Quenches in Chern Insulators,” Phys. Rev. Lett. 115, 236403 (2015).
- Wang et al. (2016) P. Wang, M. Schmitt, and S. Kehrein, “Universal nonanalytic behavior of the Hall conductance in a Chern insulator at the topologically driven nonequilibrium phase transition,” Phys. Rev. B 93, 085134 (2016).
- Ünal et al. (2016) F. Nur Ünal, E. J. Mueller, and M. Ö. Oktel, “Nonequilibrium fractional Hall response after a topological quench,” Phys. Rev. A 94, 053604 (2016).
- Schmitt and Wang (2017) M. Schmitt and P. Wang, “Universal non-analytic behavior of the non-equilibrium Hall conductance in Floquet topological insulators,” arXiv:1703.01113 [cond-mat] (2017).
- Tran et al. (2017) D. T. Tran, A. Dauphin, A. G. Grushin, P. Zoller, and N. Goldman, “Probing topology by ”heating”,” arXiv:1704.01990 [cond-mat] (2017).
- Qi et al. (2008) X.-L. Qi, T. L. Hughes, and S.-C. Zhang, “Topological field theory of time-reversal invariant insulators,” Phys. Rev. B 78, 195424 (2008).
- Liu et al. (2008) C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, “Quantum Anomalous Hall Effect in Quantum Wells,” Phys. Rev. Lett. 101, 146802 (2008).
- Bernevig et al. (2006) B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, “Quantum Spin Hall Effect and Topological Phase Transition in HgTe Quantum Wells,” Science 314, 1757–1761 (2006).
- (18) Even though the term topological insulator (TI) refers to a more general concept than the QAH insulator, we use the abbreviation TI throughout the text.
- Heyl and Budich (2017) M. Heyl and J. C. Budich, “Dynamical Topological Quantum Phase Transitions for Mixed States,” arXiv:1705.08980 [cond-mat, physics:quant-ph] (2017).
- Eckstein and Werner (2013) M. Eckstein and P. Werner, “Photoinduced States in a Mott Insulator,” Phys. Rev. Lett. 110, 126401 (2013).
- Golež et al. (2017) D. Golež, L. Boehnke, H. U. R. Strand, M. Eckstein, and P. Werner, “Nonequilibrium : Antiscreening and Inverted Populations from Nonlocal Correlations,” Phys. Rev. Lett. 118, 246402 (2017).
- Lindner et al. (2011) N. H. Lindner, G. Refael, and V. Galitski, “Floquet topological insulator in semiconductor quantum wells,” Nature Phys. 7, 490–495 (2011).
- D’Alessio and Rigol (2015) L. D’Alessio and M. Rigol, “Dynamical preparation of Floquet Chern insulators,” Nature Communications 6, 9336 (2015).
- Wolff et al. (2016) S. Wolff, A. Sheikhan, and C. Kollath, “Dissipative time evolution of a chiral state after a quantum quench,” Phys. Rev. A 94, 043609 (2016).
- van Leeuwen (2004) R. van Leeuwen, “First-principles approach to the electron-phonon interaction,” Phys. Rev. B 69, 115110 (2004).
- Ouyang and Hu (2015) T. Ouyang and M. Hu, “First-principles study on lattice thermal conductivity of thermoelectrics HgTe in different phases,” Journal of Applied Physics 117, 245101 (2015).
- Giraud and Egger (2011) S. Giraud and R. Egger, “Electron-phonon scattering in topological insulators,” Phys. Rev. B 83, 245322 (2011).
- Giraud et al. (2012) S. Giraud, A. Kundu, and R. Egger, “Electron-phonon scattering in topological insulator thin films,” Phys. Rev. B 85, 035441 (2012).
- Costache et al. (2014) M. V. Costache, I. Neumann, J. F. Sierra, V. Marinova, M. M. Gospodinov, S. Roche, and S. O. Valenzuela, “Fingerprints of Inelastic Transport at the Surface of the Topological Insulator BiSe: Role of Electron-Phonon Coupling,” Phys. Rev. Lett. 112, 086601 (2014).
- Hatch et al. (2011) R. C. Hatch, M. Bianchi, D. Guan, S. Bao, J. Mi, B. B. Iversen, L. Nilsson, L. Hornekær, and P. Hofmann, “Stability of the BiSe topological state: Electron-phonon and electron-defect scattering,” Phys. Rev. B 83, 241303 (2011).
- Pan et al. (2012) Z.-H. Pan, A. V. Fedorov, D. Gardner, Y. S. Lee, S. Chu, and T. Valla, “Measurement of an Exceptionally Weak Electron-Phonon Coupling on the Surface of the Topological Insulator BiSe Using Angle-Resolved Photoemission Spectroscopy,” Phys. Rev. Lett. 108, 187001 (2012).
- Sentef et al. (2013) M. A. Sentef, A. F. Kemper, B. Moritz, J. K. Freericks, Z.-X. Shen, and T. P. Devereaux, “Examining Electron-Boson Coupling Using Time-Resolved Spectroscopy,” Physical Review X 3, 041033 (2013).
- Kemper et al. (2013) A. F. Kemper, M. Sentef, B. Moritz, C. C. Kao, Z. X. Shen, J. K. Freericks, and T. P. Devereaux, “Mapping of unoccupied states and relevant bosonic modes via the time-dependent momentum distribution,” Phys. Rev. B 87, 235139 (2013).
- Murakami et al. (2015) Y. Murakami, P. Werner, N. Tsuji, and H. Aoki, “Interaction quench in the Holstein model: Thermalization crossover from electron- to phonon-dominated relaxation,” Phys. Rev. B 91, 045128 (2015).
- Säkkinen et al. (2015) N. Säkkinen, Y. Peng, H. Appel, and R. van Leeuwen, “Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer,” J. Chem. Phys. 143, 234102 (2015).
- Murakami et al. (2016) Y. Murakami, P. Werner, N. Tsuji, and H. Aoki, “Damping of the collective amplitude mode in superconductors with strong electron-phonon coupling,” Phys. Rev. B 94, 115126 (2016).
- Schüler et al. (2016) M. Schüler, J. Berakdar, and Y. Pavlyukh, “Time-dependent many-body treatment of electron-boson dynamics: Application to plasmon-accompanied photoemission,” Phys. Rev. B 93, 054303 (2016).
- Tuovinen et al. (2016) R. Tuovinen, N. Säkkinen, D. Karlsson, G. Stefanucci, and R. van Leeuwen, ‘‘Phononic heat transport in the transient regime: An analytic solution,” Phys. Rev. B 93, 214301 (2016).
- Breuer and Petruccione (2002) H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 2002).
- Usenko et al. (2016) S. Usenko, M. Schüler, A. Azima, M. Jakob, L. L. Lazzarino, Y. Pavlyukh, A. Przystawik, M. Drescher, T. Laarmann, and J. Berakdar, “Femtosecond dynamics of correlated many-body states in C 60 fullerenes,” New J. Phys. 18, 113055 (2016).
- Rosati et al. (2014) R. Rosati, R. C. Iotti, Fabrizio Dolcini, and F. Rossi, “Derivation of nonlinear single-particle equations via many-body Lindblad superoperators: A density-matrix approach,” Phys. Rev. B 90, 125140 (2014).
- Rosati et al. (2015) R. Rosati, F. Dolcini, and F. Rossi, “Electron-phonon coupling in metallic carbon nanotubes: Dispersionless electron propagation despite dissipation,” Phys. Rev. B 92, 235423 (2015).
- Schlünzen and Bonitz (2016) N. Schlünzen and M. Bonitz, “Nonequilibrium Green Functions Approach to Strongly Correlated Fermions in Lattice Systems,” Contributions to Plasma Physics 56, 5–91 (2016).
- Langreth and Nordlander (1991) D. C. Langreth and P. Nordlander, “Derivation of a master equation for charge-transfer processes in atom-surface collisions,” Phys. Rev. B 43, 2541 (1991).
- Hwang et al. (2011) C. Hwang, C.-H. Park, D. A. Siegel, A. V. Fedorov, S. G. Louie, and A. Lanzara, “Direct measurement of quantum phases in graphene via photoemission spectroscopy,” Phys. Rev. B 84, 125422 (2011).
- Liu et al. (2011) Y. Liu, G. Bian, T. Miller, and T.-C. Chiang, ‘‘Visualizing Electronic Chirality and Berry Phases in Graphene Systems Using Photoemission with Circularly Polarized Light,” Phys. Rev. Lett. 107, 166803 (2011).
- Duca et al. (2015) L. Duca, T. Li, M. Reitter, I. Bloch, M. Schleier-Smith, and U. Schneider, “An Aharonov-Bohm interferometer for determining Bloch band topology,” Science 347, 288–292 (2015).
- Wu et al. (2016) Zhan Wu, L. Zhang, W. Sun, X.-T, Xu, B.-Z. Wang, S.-C. Ji, Y. Deng, S. Chen, X.-J. Liu, and J.-W. Pan, “Realization of two-dimensional spin-orbit coupling for Bose-Einstein condensates,” Science 354, 83–88 (2016).
- Moskalenko et al. (2017) A. S. Moskalenko, Z.-G. Zhu, and J. Berakdar, “Charge and spin dynamics driven by ultrashort extreme broadband pulses: A theory perspective,” Phys. Rep. Charge and spin dynamics driven by ultrashort extreme broadband pulses: a theory perspective, 672, 1–82 (2017).
- Moskalenko et al. (2006) A. S. Moskalenko, A. Matos-Abiague, and J. Berakdar, “Revivals, collapses, and magnetic-pulse generation in quantum rings,” Phys. Rev. B 74, 161303 (2006).
- You et al. (1993) D. You, D. R. Dykaar, R. R. Jones, and P. H. Bucksbaum, ‘‘Generation of high-power sub-single-cycle 500-fs electromagnetic pulses,” Optics Lett. 18, 290–292 (1993).
- Jones et al. (1993) R. R. Jones, D. You, and P. H. Bucksbaum, “Ionization of Rydberg atoms by subpicosecond half-cycle electromagnetic pulses,” Phys. Rev. Lett. 70, 1236–1239 (1993).
- Schmitt et al. (2008) F. Schmitt, P. S. Kirchmann, U. Bovensiepen, R. G. Moore, L. Rettig, M. Krenz, J.-H. Chu, N. Ru, L. Perfetti, D. H. Lu, M. Wolf, I. R. Fisher, and Z.-X. Shen, “Transient Electronic Structure and Melting of a Charge Density Wave in TbTe3,” Science 321, 1649–1652 (2008).
- Graf et al. (2011) J. Graf, C. Jozwiak, C. L. Smallwood, H. Eisaki, R. A. Kaindl, D.-H. Lee, and A. Lanzara, “Nodal quasiparticle meltdown in ultrahigh-resolution pump-probe angle-resolved photoemission,” Nature Phys. 7, 805–809 (2011).
- Bovensiepen and Kirchmann (2012) U. Bovensiepen and P.s. Kirchmann, “Elementary relaxation processes investigated by femtosecond photoelectron spectroscopy of two-dimensional materials,” Laser & Photonics Reviews 6, 589–606 (2012).
- Smallwood et al. (2012) C. L. Smallwood, J. P. Hinton, C. Jozwiak, W. Zhang, J. D. Koralek, H. Eisaki, D.-H. Lee, J. Orenstein, and A. Lanzara, “Tracking Cooper Pairs in a Cuprate Superconductor by Ultrafast Angle-Resolved Photoemission,” Science 336, 1137–1139 (2012).
- Stan et al. (2009) A. Stan, N. E. Dahlen, and R. van Leeuwen, “Time propagation of the Kadanoff-Baym equations for inhomogeneous systems,” J. Chem. Phys. 130, 224101 (2009).