| J. Phys. A: Math. Theor. 41 No 39 (3 October 2008) 392003 (9pp) |
| doi:10.1088/1751-8113/41/39/392003 |
Quantum hydrodynamics and nonlinear differential equations for degenerate Fermi gas
Eldad Bettelheim1, Alexander G Abanov2 and Paul B Wiegmann3
1 Racah Institute of Physics, The Hebrew University of Jerusalem, Safra Campus, Givat Ram, Jerusalem 91904, Israel
2 Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA
3 James Franck Institute of the University of Chicago, 5640 S Ellis Avenue, Chicago, IL 60637, USA
E-mail: eldadb@phys.huji.ac.il, alexandre.abanov@sunysb.edu and wiegmann@uchicago.edu
Received 26 June 2008, in final form 7 August 2008
Published 2 September 2008
| Abstract. We present new nonlinear differential equations for spacetime correlation functions of Fermi gas in one spatial dimension. The correlation functions we consider describe non-stationary processes out of equilibrium. The equations we obtain are integrable equations. They generalize known nonlinear differential equations for correlation functions at equilibrium [1–4] and provide vital tools for studying non-equilibrium dynamics of electronic systems. The method we developed is based only on Wick's theorem and the hydrodynamic description of the Fermi gas. Differential equations appear directly in bilinear form. PACS numbers: 02.30.Ik, 71.10.Pm, 05.45.Yv |
1. Introduction
Various spacetime correlation functions of free Fermi gas in one spatial dimension (and models related to it—impenetrable Bose and Fermi gases, XY-spin chain, 2D Ising model, etc) are Fredholm determinants or their minors and, therefore, generally speaking, obey integrable nonlinear differential equations with respect to spacetime and other parameters such as temperature, chemical potential, etc. In practice, a derivation of these equations is a complicated task. Few equations are known. However, once obtained, they are indispensable in studies of off-shell properties of Fermi gas.
Historically nonlinear differential equations for correlation functions first appeared in studies of the 2D Ising model [1]. In [2] these equations have been extended, applied to the XY-spin chain, and have been derived in explicitly bilinear form with the help of the Wick theorem generalized for these purposes in [5]. In [3] it has been shown that at zero temperature the two-point, equal-time correlation function of impenetrable bosons and of the XY-spin chain can be expressed through Painlevé transcendents. These equations have been extended to time- and temperature-dependent correlation functions at equilibrium for impenetrable bosons and the XY-chain in [4, 6].
The derivation of the known equations was essentially based on equilibrium and translational invariance properties of correlation functions. Although they have been recognized as integrable, their relations to integrable hierarchies remained unclear.
Recently, interest focused on electronic systems out of equilibrium [7] and especially on non-stationary properties of propagating localized states. There the physical properties change drastically and the dynamics acquires essentially nonlinear features such as hydrodynamic instabilities [9]. Apart from being of a fundamental interest, research in non-equilibrium degenerate Fermi gas is driven by a quest to transmit quantum information in electronic nanodevices and to control and manipulate entangled quantum many-body states. Integrable equations describing the non-equilibrium—non-stationary states of electronic gas, once available, would be a valuable tool in studies of a complex dynamics of electronic gas. They are derived in this communication.
We will obtain nonlinear differential equations with respect to spacetime for the correlation functions of vertex operators
where
R,L are the right/left handed chiral Bose fields of fermionic currents and aL,R are arbitrary (not necessarily real) parameters. The definition of these fields in terms of free fermions is given below. Equations with respect to parameters of the density matrix (e.g., temperature or chemical potential) can also be obtained, but are not discussed here.
The correlation functions of vertex operators (1) are important in many applications. In the case of equilibrium they also satisfy some previously known equations. A few particular cases are worth mentioning. The case aL = aR = 1/2 describes a one-point function of impenetrable bosons [3, 4]. The imaginary aL = aR and aL = –aR yield a generating function of moments of the distribution function of the density and velocity of electronic gas. The chiral correlation functions aL = 0 or aR = 0 describe a tunneling into an electron gas (fermi edge singularity) [9], etc.
In this communication, we give a simple physical derivation of equations on these correlation functions. Our derivation is based solely on a general form of the Wick theoremNote1, the representation of Fermi-operators through Bose fields (see, e.g., [10]), and the hydrodynamic description of the Fermi gas [11]. Our equations appear directly in Hirota's bilinear formNote2. They are known as modified Kadomtsev–Petviashvili (or mKP) equations and are the first and the second equation of mKP hierarchy. In a particular case, in free fermion systems these equations have previously appeared in [8, 12]. The nonlinear equations being integrable, are proven to be an effective tool in computing asymptotes of correlation functions in different regimes. These computations are specific to each particular problem and will not be discussed here (see [8] for some practical applications).
In the following after necessary preliminaries on the Fermi–Bose correspondence, we describe the correlation functions (5), present the nonlinear differential equations they obey (16), (18) and (19), and then give a short derivation of those equations.
2. Fermions and a Bose field
We consider free fermions on a circle of a circumference L,
where ψp is a mode of a fermion field
.
A central point of our approach is the hydrodynamic description of the Fermi gas. We briefly review it. The properties of the Fermi gas are fully described in terms of canonical hydrodynamic variables, density and velocity,
where the chiral (right, left) currents are
Tomonaga's equal time commutation relations of currents (see, e.g., [10])
lead to a canonical relation between density and velocity,
and to the canonical Bose field:
We will be interested in the spacetime dependence of correlations of two vertex operators:
Here
is the density matrix and ξ = (x, t). The colon denotes normal ordering with respect to the vacuumNote3.
Below we drop the chiral subscripts R, L in all formulae and assume that the upper (lower) sign corresponds to the right (left) sectors.
3. Coherent states
We will consider correlation functions with respect to coherent states. This means that the density matrix is an element of gl(∞), i.e., it is an exponent bilinear in fermionic modes:
This choice of the density matrix allows us to use Wick's theorem.
The set of coherent states is rather general. It exhausts most interesting applications. For example, the density matrix can be a Boltzmann distribution, or any other non-equilibrium distribution as in [7]. Another choice would be
= |V![]()
V|, where |V
= exp(∑k tkJk)|0
is a coherent state (i.e., an eigenstate of current operators). Such states appear as a result of shake up of the electronic gas by electromagnetic potential with harmonics tk.
4. Quantum hydrodynamics
The Hamiltonian of free fermions (2) in the sector of coherent states can be expressed solely in terms of density and velocityNote4:
The evolution of density and velocity follows
The hydrodynamic equations are decoupled into two independent chiral quantum Riemann equations
5. Vertex operator and correlation functions
The Heisenberg evolution equations for the chiral vertex operators follow from the quantum Riemann equation (7) (see appendix A):
where T = :(∇
)2: – i∇2
and
are holomorphic (antiholomorphic) components of the stress–energy tensor of a chiral Bose field (with the central charge 1/2).
A particular consequence of these equations is that the vertex operators
obey the Schroedinger equation. This yields to a familiar formula representing fermions as exponents of a Bose field:
A useful representation for current J = ±:e–i
(–i∇) ei
: and the stress–energy tensor reads
Current and stress–energy tensor can be further cast in fermionic form using (10) and OPE for vertex operators as follows from the current algebra (4):Note5
With the help of (12), (11) and (10) we write
Finally, we cast the equation of motion of the vertex operator in a suggestive bilinear form (quantum Hirota equation)
Here D is the Hirota derivative Df · g = ∂fg – f∂g and D2f · g = ∂2fg – 2∂f∂g + f∂2g. This equation follows from the quantum Riemann equation (7).
6. Nonlinear differential equations
The nonlinear equations themselves do not depend on the parameters aL,R, bL,R in (5). For brevity we put L = 2π in formulae below and consider only the chiral and neutral (a = b) correlation function, which we denote as
Being presented in a bilinear form the differential equations for correlation functions look identical to the quantum equation (14):
The equations must be viewed as equations for τa + n(x, t; x ', t ') in the continuous variables x, t and the discrete variable n, where |a| ≤ 1/2. As such they are a closed set of equations known as the mKP equations. The correlation function can be found given initial conditions at t = 0, for all integer n and any x. Certain asymptotes or analytical conditions may also be helpful in solving the problem for particular cases. This program is taken up in practice in [8], where in [8] the density matrix had the form |A![]()
0| where |0
is the ground state of electron gas and |A
is a chiral coherent state of the form |A
=
|0
. We also comment that τa where |a| ≤ 1/2 is a Fredholm determinant, while τa + n are its minors.
The mKP equations (16) have many different forms discussed in the literature. We will mention one of them. Setting
equation (16) becomes
In this form the classical equation (17) resembles the quantum equation (7). The last term in (17) can be considered as a quantum correction to the semiclassical Riemann equations.
In addition to the evolution (16) there is another equation which does not involve time evolution. This equation connects correlation functions with different density matrices:
and
, where apq is an arbitrary non-degenerate matrix. Then τa and
are related by
This equation is a minor generalization of the familiar 2D Toda lattice equation.
An important particular case occurs when the modes p, q in
lie very far from the Fermi surface. Then
and equation (18) becomes the 1D Toda chain equation
This equation, together with (16) being written for pairs τa, τa + 1 and τa–1, τa give a closed set which involves only τa–1, τa, τa + 1.
7. Translation invariant states
Further closure can be achieved for special cases. For example, if the density matrix commutes with momentum, the correlation function depends on the difference x – x ' (e.g., the density matrix depends only on occupation numbers np = ψ†pψp,
= exp ∑pλp np. Then it also commutes with the Hamiltonian, is a function of t – t ', and describes stationary processes). In this case, we denote the tau function as
As follows from (16), (19) three functions Fa–1, Fa, Fa + 1 obey three equations with respect to (x, t):
Here we used
.
These equations are the bilinear form of the nonlinear Schroedinger equation (without a complex involution). Introducing Ψ = Fa + 1/Fa and
we have
These equations were obtained for time- and temperature-dependent two-point correlation functions for impenetrable bosons at equilibrium (
= e–H/T) [4]. The impenetrable bosons are equivalent to free fermions. The creation operator of the boson is the vertex operator (1) with aL = aR = 1/2.
8. Derivation of nonlinear equations
The fact that the quantum equation for vertex operators (14) and the classical equations for the correlation functions (16) have the same form may not be an accident but may reflect a general property of integrable equations written in the bilinear form. Here we limit our discussion by demonstrating this phenomena for equation (16).
We now prove the main formula (16). First, we adopt the notation
where only the spatial dependence has been specified on the rhs.
Let us multiply (8) and (9) by
and
respectively, take the trace and substitute into equation (16). We find that (16) holds due to the identity
where 1 is the identity operator, and we further suppressed the arguments x, x '.
To prove this identity we first write
in terms of fermions. Using (11) we write
, where we split spatial points. We now use formula (12) for vertex operators and equation (11) to write
This formula allows us to write the first term in (22) in terms of a correlator with four fermions insertion:
Now we use a general form of the Wick theoremNote6 to express the four fermion correlator in terms of correlators containing two fermions:
After this, we trace back the path which led us to (24). With the help of equations (10) and (12), we obtain the formulae
We use them in order to write the rhs of (25) in terms of bosons. After that, one applies the operation
to the rhs of (25), takes the derivative and merges the points. This leads to (22). The calculations are somewhat cumbersome, but straightforwardNote7.
We briefly comment on the proof of (19), the proof of the more general (18) goes along the same lines. First, we write (19) in the form
The proof of this equation is easier: after replacing bosonic exponentials and the current operator by fermionic operators according to (10), (13) one recognizes Wick's theorem.
9. Integrable hierarchy of nonlinear equations
As a final remark we mention that higher equations of the mKP hierarchy describe the evolution of Fermi gas with an arbitrary spectrum H = ∑pεpψ†pψp. They are generated by the Hirota equation described in [12]. Equations for magnetic chains equivalent to the Fermi gas can be obtained in a similar manner.
Acknowledgments
PW thanks V Korepin for discussion of the results of this paper. We thank J H H Perk for pointing out to us and discussing the results of [2, 5]. EB was supported by ISF grant number 206/07. The work of AGA was supported by the NSF under the grant DMR-0348358. PW was supported by NSF under the grants DMR-0540811 and MRSEC DMR-0213745.
Appendix A. Equation of motion for vertex operators
In this appendix we sketch the derivation of equations (8), (9) from the Riemann equation (7). First of all we show how to calculate the time derivative of the vertex operator. A subtlety here is that
does not commute with
.
We recall a simple consequence of the Hadamard lemma. If H and A are two operators such that [[H, A], A] commutes with A, then
If H is a Hamiltonian ∂t A ~ [H, A] and we obtain for the time derivative of an exponent
Let us now compute an evolution of eia
. We assume that
is a right chiral field. The left chiral field obeys the identical equation.
We note that the commutator
vanishes, and then apply (A.1). We obtain
Let us now compute the commutator
. Using (7) and (3) we obtain for Fourier components
Then we proceed as
Repeating the same calculation for
– we obtain
Finally, using (A.4) in (A.2) we obtain (8), (9)
References
Notes
, where
and
.
, where fermionic operators are placed in the middle between the positive and negative bosonic modes.
Eldad Bettelheim et al 2008 J. Phys. A: Math. Theor. 41 392003
J McKenna et al 2009 J. Phys. B: At. Mol. Opt. Phys. 42 121003
G. Bryden et al. 2009 ApJ 705 1226
A. V. Rybin et al 2008 EPL 81 40009
Emilio Elizalde and Jaume Haro 2009 J. Phys. A: Math. Theor. 42 472001
The HERMES collaboration JHEP11(2009)083
J. Hwang et al 2008 EPL 82 27002
Gary T. Horowitz and Juan Maldacena JHEP02(2004)008
Roberto Iengo et al JHEP11(2009)020
Pin-Gao Gu and Takeru K. Suzuki 2009 ApJ 705 1189