Arxiv:2104.02428V2 [Cond-Mat.Quant-Gas] 11 Apr 2021 Physics [13–19] Has Become Quite Active, Mostly for Its General Theory
Total Page:16
File Type:pdf, Size:1020Kb
Kubo formula for non-Hermitian systems and tachyon optical conductivity Doru Sticlet,1, ∗ Bal´azsD´ora,2, 3 and C˘at˘alinPa¸scuMoca4, 5 1National Institute for R&D of Isotopic and Molecular Technologies, 67-103 Donat, 400293 Cluj-Napoca, Romania 2MTA-BME Lend¨ulet Topology and Correlation Research Group, Budapest University of Technology and Economics, 1521 Budapest, Hungary 3Department of Theoretical Physics, Budapest University of Technology and Economics, 1521 Budapest, Hungary 4MTA-BME Quantum Dynamics and Correlations Research Group, Institute of Physics, Budapest University of Technology and Economics, 1521 Budapest, Hungary 5Department of Physics, University of Oradea, 410087, Oradea, Romania Linear response theory plays a prominent role in various fields of physics and provides us with extensive information about the thermodynamics and dynamics of quantum and classical systems. Here we develop a general theory for the linear response in non-Hermitian systems with non-unitary dynamics and derive a modified Kubo formula for the generalized susceptibility for arbitrary (Her- mitian and non-Hermitian) system and perturbation. As an application, we evaluate the dynamical response of a non-Hermitian, one-dimensional Dirac model with imaginary and real masses, per- turbed by a time-dependent electric field. The model has a rich phase diagram, and in particular, features a tachyon phase, where excitations travel faster than an effective speed of light. Surprisingly, we find that the dc conductivity of tachyons is finite, and the optical sum rule is exactly satisfied for all masses. Our results highlight the peculiar properties of the Kubo formula for non-Hermitian systems and are applicable for a large variety of settings. Introduction. Linear response theory is a fundamen- periments, e.g. in a quantum gas microscopy experi- tal framework in statistical physics that describes the way ment [35, 36], the measurement backaction can quali- in which a physical system responds to small external tatively alter the low-energy physics, so it is crucial to perturbations [1]. The response is expressed in term of understand and describe the effect of small perturba- a generalized susceptibility of the unperturbed system. tions in such non-Hermitian systems. In this regard, The change in the expectation value of an operator A only recently, a non-Hermitian perturbation of a Her- due to a time-dependent perturbation V (t) = Bf(t) is mitian Hamiltonian has been considered, and a modi- assumed to be linear in the perturbative force f(t), fied Kubo formula for the generalized susceptibility has Z emerged [37]. Still, at present, a complete non-Hermitian 0 0 0 δhA(t)i = dt χAB(t; t )f(t ) ; (1) linear response theory, in which the Hamiltonian gov- erning the system and the perturbation are both non- 0 with χAB(t; t ) known as the response function. In gen- Hermitian is missing. eral, χ (t; t0) is real, and for systems invariant under AB In the present work we address this problem, and fill time translations χ (t; t0) = χ (t − t0). Traditionally, AB AB this gap, by developing a general formalism and pro- the explicit expression for the response function in Eq. (1) viding a generalized Kubo formula for the susceptibility is given by the Kubo formula [2], and has been for decades χ (t; t0) introduced in Eq. (1) which is suitable for non- a standard tool in the study of transport phenomena in AB Hermitian systems under very general conditions. We il- typical many-body [3{5], closed [6,7] or open systems [8{ lustrate our approach by computing the correlation func- 10]. Through the fluctuation-dissipation theorem, the tion for the current operator, in response to an exter- linear response provides information about the equilib- nal electric field, and the associated optical conductiv- rium fluctuations of a system, and can reveal e.g. topo- ity for a non-Hermitian tachyon model. Although non- logical properties through the quantum Hall response [11] Hermiticity is in general associated with dissipation, we or demonstrate the presence of fractional charge carriers find remarkable features such as a finite dc conductivity for the fractional quantum Hall effect via noise measure- in the absence of scatterers and the conservation of the ments [12]. optical conductivity sum rule. Recently, a newly emerging field of non-Hermitian arXiv:2104.02428v2 [cond-mat.quant-gas] 11 Apr 2021 physics [13{19] has become quite active, mostly for its General theory. We consider a general setup in which promising applications [17] and for exhibiting unique fea- the system is modeled by a time-independent non- tures, which are absent in the Hermitian realm, such Hermitian Hamiltonian H0 with non-degenerate spec- as exceptional points [20{24], non-Hermitian skin ef- trum subjected to a time-dependent perturbation V (t) = fect [25{29], non-Bloch phase transitions [30], and uni- Bf(t), which starts to act at t = 0, where B is an op- directional invisibility [31] to mention a few. By now, erator which can be non-Hermitian and f(t), the pertur- quantum optics and cold atoms setups provide fertile bation force, assumed to be small and real-valued. We grounds where such non-Hermitian concepts are probed are interested in the linear response of a system operator experimentally [21, 32{34]. In such state-of-the-art ex- δhA(t)i = hA(t)i− hAi0. The expectation value hA(t)i at 2 time t is given by bation as tr[δρI (t)AI (t)] tr[δρ(t)] δhA(t)i = − hA(t)i0 ; (3) tr[ρ(t)A] tr[ρ (t)] tr[ρ (t)] hA(t)i = ; (2) 0 0 tr[ρ(t)] y iH t= −iH0t= where AI (t) = e 0 ~Ae ~ is written in the interac- tion representation in terms of the unperturbed system, y iH0t= −iH t= where ρ(t) is the density matrix of the system in the and δρI (t) = e ~δρ(t)e 0 ~. The second, norm Schr¨odingerrepresentation obeying the non-Hermitian correction term, is specific to non-unitary dynamics and y von Neumann equation i~@tρ(t) = Hρ(t) − ρ(t)H , with is absent in the regular unitary evolution, but it is impor- H = H0 + V (t), the full Hamiltonian. Its solution tant to reproduce correctly the dynamics of the expecta- is decomposed as ρ(t) = ρ0(t) + δρ(t) with ρ0(t) = tion values [39]. Keeping in mind that AI (t) is identical y −iH0t= iH t= e ~ρ(0)e 0 ~, the density matrix of the unper- to the Heisenberg picture time evolution in terms of H0, turbed system, and δρ(t) arising from the external per- i.e. A(t) = AI (t), we finally obtain a compact expression turbation. Rotating to the interaction picture, the varia- for the generalized susceptibility [38] defined in Eq. (1) tion δhA(t)i is evaluated [38] to first order in the pertur- as 0 i y ρ (t ) 0 iH τ= −iH0τ= 0 χAB(t; t ) = − θ(τ)tr [A(τ);B]∼ − hA(t)i0[e 0 ~e ~;B]∼ · ; (4) ~ tr[ρ0(t)] in terms of a modified commutator defined as [X; Y ]∼ = system Hamiltonian in momentum space is XY − Y yX, where τ = t − t0. The first term is the X 2 natural generalization of the Hermitian Kubo formula H0 = p c σx + ∆σy − imc σz; (5) for the non-Hermitian setting. The second term arises p solely from the non-unitary dynamics, and hA(t)i0 = which has two momentum-dependent energy bands p 2 2 2 2 4 tr[ρ0(t)A]=tr[ρ0(t)]. In contrast to the Hermitian Kubo E±(p) = ±"(p), with "(p) = p c + ∆ − m c [see formula, the non-Hermitian counterpart is not time- Fig.1(a)], with the Fermi velocity as the effective light 0 translation invariant and depends separately on t and t speed c. The Hamiltonian is invariant with respect to a due to the appearance of the system density matrix ρ0(t) PT -symmetry operator σxK, with K performing com- at various times, which may have a non-unitary evolution plex conjugation, [σxK;H0] = 0. The gapped non- in the absence of the perturbation. Hermitian Dirac Hamiltonian at j∆j > mc2 is in a PT - symmetric phase since the eigenvectors are invariant to When the system Hamiltonian H0 is Hermitian, the the PT -symmetry operator and the spectrum consists expression (4) simplifies considerably since ρ0(t) = ρ(0) of two real energy bands [13, 42, 43]. In this regime, is a time-independent equilibrium density matrix. For P the density of states, g(!) = p δ(! − E±(p)) van- example, if additionally, the perturbation proves Her- ishes below the gapp edge and features a sharp thresh- mitian, the norm corrections drop out, and assum- old singularity at ∆2 − m2c4. Exactly at j∆j = mc2, ing that the system is invariant under time transla- the spectrum becomes gapless with a linear dispersion tion, one readily recovers the regular Kubo formula [1], "(p) = cjpj. The p = 0 point marks an exceptional i χ(τ) = − θ(τ)h[A(τ);B]i0. If the perturbation is anti- point (EP), where not only the eigenvalues, but also the ~ Hermitian By = −B, the correlation function recov- eigenvectors coalesce [20, 22]. Due to the linear disper- ers the central result of Ref. [37], where the modi- sion, the density of states is constant for all energies. fied commutator becomes an anticommutator, χ(τ) = For j∆j < mc2, a tachyon-like spectrum with a hyper- i − θ(τ)(hfA(τ);Bgi0 − 2hAi0hBi0). If H0 is non- bolic gapless dispersion develops. Depending on p, the ~ Hermitian, the norm correction term is crucial and it system is in a PT -symmetric phase with real eigenval- cannot be discarded as it may encompass the dominant ues for p2c2 > m2c4 − ∆2 or in a PT -symmetry broken contribution to linear response.