Hall conductance for open two-band system beyond rotating-wave approximation.

The response of the open two-band system to external fields would in general be different from that of a strictly isolated one. In this paper, we systematically study the Hall conductance of a two-band model under the influence of its environment by treating the system and its environment on equal footing. In order to clarify some well-established conclusions about the Hall conductance, we do not use the rotating wave approximation (RWA) in obtaining an effective Hamiltonian. Specifically, we first derive the ground state of the whole system (the system plus the environment) beyond the RWA, then calculate an analytical expression for Hall conductance of this open system in the ground state. We apply the expression to two examples, including a magnetic semiconductor with Rashba-type spin-orbit coupling and an electron gas on a square two-dimensional lattice. The calculations show that the transition points of topological phase are robust against the environment. Our results suggest a way to the controlling of the whole system response, which has potential applications for condensed matter physics and quantum statistical mechanics.

Introduction

Edwin. H. Hall discovered a famous phenomenon, Hall Effect, that when a conductor carrying longitudinal current was placed in a vertical magnetic field, the carriers would be pressed towards the transverse side of the conductor, which led to observed transverse voltage[1]. Its quantum counterpart–the quantum Hall effect (QHE)[2] and the later studies on the geometric phase and the topological properties of quantum Hall states promoted the study of anomalous Hall effect, which together with the Hall and spin Hall effects completed the quantum Hall trio. The integer quantum Hall effect can be understood in the single-particle framework[3,4] and can be characterized by the topological invariant called Chern number[5,6] that was introduced in the QHE to explain the jump of the Hall conductance in magnetic fields. Niu, Thouless, and Wu found, for the 2D quantum Hall, a topological invariant (the first Chern number) expressed in terms of the ground-state wave function, which is valid in the presence of an arbitrary interaction and disorder[6]. Hall conductance can be represented in terms of the topological invariant (or Chern number) in the linear response theory.

The response of topological insulators (TIs) to an external weakly classical field can be expressed in terms of Kubo formula[7], which predicts quantized Hall conductivity of the quantum Hall family. TIs were theoretically predicted and have been experimentally discovered[8-12]. In contrast to ordinary band insulators, TIs are a broad class of unconventional materials that are insulating in the interior but conductive along edges. Although time-reversal (TR) invariance is essential in the quantum spin Hall (QSH) insulator, there is a TR symmetry-breaking state of matter which is closely related to the QSH insulator: the quantum anomalous Hall (QAH) insulator[12,13]. The quantum spin Hall effect can be understood as two copies of the quantum anomalous Hall effect; each breaks TR symmetry while the whole system remains invariant under TR. In the past few years, scientific community has paid much attention to these topological materials due to their specific features. Most recently, efforts have made to investigate topological materials for open systems. For example, zero-temperature Hall conductance subjected to decoherence[14], topological order by dissipation[15-17], density-matrix Chern insulator subjected to thermal noise[18,19], topological phases induced photocurrent[20,21], and optical Hall conductivity[22-25]. Taking the quantum nature of the external field into account, a Hall conductance to characterize the linear response of a quantum system to a single-mode quantized electronmagnetic field was also defined and explored[26]. However, most of these studies are based on the reduced density matrix of the system, hence the environment is traced out before proceeding to analyze the response of the system. In many works on the linear response theory of open quantum systems, it is assumed that the effective Hamiltonian is obtained within the RWA[19,27]. This motivates us to develop a quantum response theory[28,29] for open system including the influence of the counter-rotating terms on Hall conductance, and apply it to study two-band model subjected to environments.

For this purpose, we will study the response of the open two-band system to the external field by calculating the Hall conductance when the whole system is in its ground state. The total system consists of a two-band model and its environment described by multi-mode quantized fields. We will consider the multi-mode fields entering the system via a vector potential[30-34].

In this paper, we will derive the ground state of the whole system without the RWA, and then calculate the Hall conductance when the whole system is in its ground state. This approach is not limited to any specific systems, as long as they can be described by the two-band model. We take the tight-binding model and an electron gas on a square two-dimensional lattice as examples to illustrate the theory. The calculations show that while the transition points of topological phase are robust against the environment, the Hall conductance is smaller than that without the environment.

Methods

System and Effective Hamiltonian

The model we study is a generic two-band Hamiltonian,where I is the identity matrix, are Pauli matrices and stands for the Bloch wave vector of the electron, and depend on the materials under study and determine its band structure. The term is just a shift of zero energy level, and then can be neglected[13] for simplicity. The two band may represent different physical degrees of freedom, for example, if they stand for the components of electron with spin 1/2, are the momentum-dependent coefficients which describe the spin-orbit interactions and exchange interaction of magnetic impurities; and if they are equal to the orbital degree of freedoms, then describes the hybridization between bands[35].

The system in a field can still be described by the two-band model by changing the crystal momentum, , where the electromagnetic field is represented by the vector potential . Considering an electromagnetic environment, the vector potential can be written as [36], stands for the unit vector. Substituting into the Hamiltonian and assuming the environment (electromagnetic background) is very weak, we may expend the Hamiltonian up to the first order in ,

Eq. (2) is the model Hamiltonian in this study. Straightforward calculations show that,where is the annihilation (creation) operator of the boson mode with frequency ω and g is the coupling between the system and environmental mode n. The total Hamiltonian H contains the counter-rotating terms, and is not exactly solvable even for the simple cases of single mode or single excitation.

The eigenstates of H 0 take

H can be diagonalized by a unitary matrix U: ,where we have omitted the term because this term represents a energy shift to the system much smaller than the energy difference in H 0 [26]. Here, , and , with , , , and

In the following, we will use the generalized version of the Fröhlich-Nakajima transformation[37-40] to eliminate the high-frequency terms in the effective Hamiltonian: , with the operator S where λ will be determined later. The transformation can be done order by order, and the terms of order and higher will be ignored[38]. Up to the second order in the system-environment couplings, the effective Hamiltonian is given by,where and respectively take

By choosing the following form for λ ,terms with and can be eliminated from . Noticing that contains only the rotating wave term, because the unitary transformation eliminate the counter-rotating terms, we write down the effective Hamiltonian H as,with . Here we have omitted the multiboson nondiagonal transition like and , because the contributions of these nondiagonal terms to the physical quantities are [37,38].

The ground state of the whole system

Obviously, the effective Hamiltonian can be solved exactly, and the ground state is

The ground state of the original Hamiltonian is then given by , leading towhere

In the continuum limit, the sum can be replaced by an integral. In the following we would consider the spectral density of the environment and a particular formHere α is the dimensionless coupling strength, the index s accounts for various physical situations. For example, the spectrum is sub-Ohmic when , and its is Ohmic when or super-Ohmic when . With these knowledge, we have

So far, we have presented a derivation for the ground state of a two-band model coupled to a multi-mode quantized field(or an environment) beyond the RWA.

Results

Hall conductance

In order to derive the Hall conductance as a response to a constant electric field , we consider the case that the electric field can be represented by a time-dependent vector potential, i.e., , and we take the field to be along the x-direction. This gives rise to a Hall current proportional and perpendicular to the electric field. The Hall conductivity, defined as the radio of the current density and the electric field, is therefore given by[41] where we use to denote an eigenstate of the Hamiltonian with energy in the v-th magnetic Bloch band. The general relationship between the momentum-space topology and the quantization of physical responses has been discussed extensively by Volovik[42]. The Hall conductance is related to the Chern number C, and the Chern number C is a topological invariant, which is robust against a local deformation of the Hamiltonian and can only take an integer value. The Chern number could describe the topological property of the ground-state wave function, and then it is a measurable physical quantity[43,44].

Substituting the ground state into Eq. (15), we obtainwhere the zeroth- and first-order of respectively takewith

Here, defines the Berry curvature of the lower bare band . Hall conductance in Eq. (16) for the two-band model subject to the environment has two different terms. The first term represents the linear response which returns to Eq. (15) when the system is closed, while the second one stands for a correction of the environment to the Hall conductance, and depends on k , k , the components of the Bloch vector.

Noticing that in Eq. (16) can be considered as the Hall conductance under the influence of environment, which is derived through solving the ground state of the whole system including the effect of the counter-rotating terms on it. As quantifies the response of the two-band model to a driving electric field, and the response is calculated by treating the system and environment on equal footing, we can claim that we have developed the other approach to study the response of open systems. This is different from the earlier studies base on the master equation, see for example[14].

In the next subsection, we will present two examples that together exemplify the response of this quantum open system to the classical electric field.

Examples

Consider a tight-binding model describing a magnetic semiconductor with Rashba-type spin-orbit coupling, spin-dependent effective mass, and a uniform magnetization on z direction[35]. This model can be described by Hamiltonian Eq. (1) with

In a 2D band structure, the integral over the Brillouin zone (BZ) of the Berry curvature of the lower bare band is a topological invariant that is the well-known Chern number[5,43]. For the tight-binding model, when , the Chern number is 1, while for , the Chern number is −1, and for and , the Chern number is 0.

The Hall conductance (in units of ) defined in Eq. (16) as a function of χ is shown in Fig. 1. The features observed from Fig. 1 demonstrate that the phase transition points, i.e., χ = 0, 2, 4, remain unchanged. The Hall conductance may not be an integer due to the presence of environment, the topological phase transition, however, survives in the environment. In contrast with the well-known Hall conductance shown in Fig. 1 (red sold line), is not a constant in regions, , , and , which results from the correction in Eq. (16). Physically, the correction comes from the counter-rotating terms, which would mix the two bands and then leads to a different response to the field. We also find that at the critical point χ = 0, the Hall conductance increases compared to the traditional Hall conductance , while it decreases at χ = 4. Except at these critical points, the Hall conductance almost always decreases due to the environment.

Figure 1

σ (in units of ) as a function of χ in the tight-binding model with . For comparison, the conventional Hall conductance (the red solid line) is also shown, corresponding to the coupling strength . The other parameters chosen are and .

In Fig. 2, we numerically show the Hall conductance against s for various spectral density characterized by the parameter χ, and the environment can change the topology of the system. The features found from Fig. 2 also support the conclusion that the topological phase transition survives for the two-mode system in environments.

Figure 2

(in units of ) for various index s of the spectral density is plotted as a function of s. In this plot, (a) When we choose and respectively, the Hall conductance at different coupling strength or of the quantized field compares with the integer Hall conductance (red solid line). The other parameters chosen are and for (b) and (c) .

The two-band model used in this paper describes a TR symmetry-breaking system[35]. The Hall conduction of the system is quantized and it is determined by the first Chern number of the Berry phase gauge field in the BZ[5,12]. Theories of topological materials are often formulated using tight-binding models as we used in this paper, which can be simulated by ultra-cold atomic gases[45] in optical lattices. Besides, scanning tunnelling microscopy (STM) can also be used to fabricate and characterize lattice structures with atomic precision in the solid state. Very recently, it has been shown that the tight-binding model can be realized in a vacancy lattice on the substrate surface[46].

In the second example, we consider a two-dimensional model with a tight-binding Hamiltonian for electron gas on a square lattice[47-49],where and are respectively the fermion annihilation and creation operators at the lattice site , respectively. is an on-site chemical potential. describes the hopping amplitude of electron from site l to its nearest neighbor site j. , and . When the model is for the case where the sites l and j are next-nearest neighbors, and , see Fig. 3.

Figure 3

The hopping amplitudes of the Hamiltonian: (1) black solid links t 1, (2) red dashed links −t 1, (3) diagonal links it 2 in the direction of the arrow.

Here, we will choose a two-site unit cell and rewrite the Hamiltonian as , with . then takes the form of Eq. (1) with , , , and . Recall that the Chern number[5,43] of this system is 1(−1) for [50].

In Fig. 4, we plot as a function of m for the different values of the coupling strength α. Note that both the Hall conductance σ in the closed system and in the open system change the sign when m crosses zero, see Fig. 4, which shows explicitly that the phase transition point for the closed and open system is at . This suggests that the topological phase transition survives in the open system considered here. We find that the environment suppresses the Hall conductance , but it does not change the phase transition points, see Fig. 4. With the increasing of α, the effect of environment on the Hall conductance becomes evident.

Figure 4

(in units of ) as a function of m at various coupling strength α in Eq. (13). In this model, , , . Here, the range of the parameter m is in . The other parameters, and have been used in this plot.

The topological invariant was introduced to explain phase transitions that are beyond the conventional framework of symmetry breaking. Clearly, taking the environment into account, the Hall conductance (equivalently the topological invariant–Chern number) may not be an integer. This will not be an obstacle to use the Hall conductance to detect topological properties of insulator states for quantum open system, since the phase transition points are still there. Vatsal Dwivedi and Victor Chua constructed a generalized transfer matrix corresponding to noninteracting tight-binding lattice models, which can subsequently be used to compute the bulk bands as well as the edge states[51]. The Hall conductance as a function of the gap parameter in the post-quench Hamiltonian displays a universal nonanalytic behavior in a generic two-band Chern insulator[13], such as the Dirac model, the Haldane model, or the Kitaev honeycomb model in the fermionic basis[44]. Though it is difficult to observe quantum transport of the surface states, which are usually covered by bulk carriers caused by material defects. Recently, several explicit experimental observations of QHE based on surface states have been obtained in topological materials[52-55]. These experimental schemes may provide us with a platform for realizing the prediction in this paper.

Discussion

We have studied the Hall conductance of two-band system in the presence of multi-mode quantized field. We can recover the traditional Hall conductance in Eq. (16) for the two-band model in closed system, whereas it can not return to the usual Hall conductance when the environment is taken into account. We treat the first-order of Hall conductance as the correction. This correction comes from the counter-rotating terms, which would mix the two bands and then leads to a different response to the external field. The Hall conductance for open system can not be written as a multiple of a Chern number and a constant, or as a weighted sum of Chern numbers, in this paper, there is no topological invariant for open systems.

It is well known that the change of the topology in the ground state of a system is accompanied by a topological phase transition. In the celebrated paper by Thouless[5], the Hall conductance (equivalently the Chern number) can be probed, and the Chern number usually keeps invariant as the Hamiltonian changes. The Chern number then has a jump, which indicates a topological phase transition in the ground state[43].

In summary, we find that the jump remains unchanged, but the Hall conductance is suppressed, when a two-band system is in contact with an environment modeled by a multi-mode quantized field. The calculation is based on the Kubo formula by treating the system and environment on equal footing. We do not use the rotating wave approximation in obtaining an effective Hamiltonian and the ground state of the whole system. This is different from the earlier study, where the starting point is the master equation for the open system within the RWA. This result suggests that the transition points of topological phase are robust against the environment. Although the analysis has been restricted to the two-band model, we believe that the observation makes the topological materials immune to the influence of the quantum field, and then supports its application into quantum optics and condensed matter physics.