Fermi Liquid Notes
Fermi Liquid Notes
Erkki Thuneberg
Department of physics
University of Oulu
1. Introduction in more length. Various generalizations of the theory are
The principal problem of physics is to determine how briefly mentioned.
bodies behave when they interact. Most basic courses of
classical and quantum mechanics treat the problem of one
or two particles or bodies. (An external potential can be
considered as one very heavy body.) The problem gets mo-
re difficult when the number of bodies involved is larger.
In particular, in condensed matter we are dealing with a
macroscopic number N ∼ 1023 of particles, and typically
hundreds of them directly interact with each other. This
problem is commonly known as the many-body problem.
There is no general solution to the many-body problem.
Instead there is a great number of approximations that
successfully explain various limiting cases. Here we discuss
one of them, the Fermi-liquid theory. This type of approxi-
mation for a fermion many-body problem was invented by
Landau (1957). It was originally proposed for liquid 3 He
at very low temperatures. Soon it was realized that a simi-
lar approach could be used to other fermion systems, most
notably to the conduction electrons of metals. The Fermi-
liquid theory allows to understand very many properties
of metals. A generalization of the Fermi-liquid theory al-
so allows to understand the superconducting state, which
occurs in many metals at low temperatures. Even when
Landau’s theory is not valid, it forms the standard against
which to compare more sophisticated theories. Thus Fermi-
liquid theory is a paradigm of many-body theories, and it
is presented in detail in many books and articles discussing
the many-body problem.
1
2. Preliminary topics gas, where we assume no interactions, V ≡ 0. Below we
concentrate on ideal spin-half (s = 1/2) Fermi gas.
Many-body problem In the absence of interactions we can assume a facto-
The many-body problem for identical particles can be rizable form
formulated as follows. Consider particles of mass m labeled Ψ0 (r1 , σ1 , r2 , σ2 , . . . ; t) = φa (r1 , σ1 )φb (r2 , σ2 ) . . . . (7)
by index i = 1, 2, . . . , N . Their locations and momenta are
written as rk and pk . The Hamiltonian is consisting of a product of single-particle wave functions
N
φα (r, σ). This does not yet satisfy the antisymmetry requi-
X p2k rement (5) but permuting the arguments in Ψ0 and sum-
H= + V (r1 , r2 , . . .). (1)
2m ming them all together multiplied by (−1)nP , where nP is
k=1
the number of pairwise permutations in a permutation P ,
Here V describes interactions between any particles, and one can generate a proper wave function
it could in many cases be written as a sum of pairwise
1 X
interactions V = V12 + V13 + . . . + V23 + . . .. The classical Ψ(r1 , σ1 , r2 , σ2 , . . .) = √ (−1)nP
many-body problem is to solve the Newton’s equations. N P
2
where npσ = 1 for an occupied state and is zero otherwise. Next we calculate the specific heat of the ideal Fermi gas
at low temperatures. The average energy is given by
The ground state of a system with N particles has N
lowest energy single-particle states occupied and others
Z
XX 2
empty. Here the maximal kinetic energy of an occupied E = p n( p ) = 3
p n(p )d3 p
σ p
(2πh̄/L)
state is called Fermi energy F . We also define the Fermi Z ∞
wave vector kF and Fermi momentum pF = h̄kF so that 8π
= 3
p2 p n(p )dp. (19)
2 2
(2πh̄/L) 0
2
p h̄ kF
F = F = . (14) Changing = p2 /2m as the integration variable we get
2m 2m
√
In momentum space this defines the Fermi surface (p = 2m3 ∞ 3/2
Z
E
pF ). All states inside the Fermi surface (p < pF ) are occu- = n()d
V π 2 h̄3 0
pied, and the ones outside are empty. Z ∞
= g() n()d. (20)
py
× × × × × × × × × × × × × × ×
0
× × × × × × × × × × × × × × ×
× × × × × × × × × × × × × × × In the second line we have expressed √ the same result by
× × × × × × × × × × × × × × ×
defining a density of states g() = m 2m/π 2 h̄3 .
pF
× × × × × × × × × × × × × × ×
× × × × × × × × × × × × × × ×
× × × × × × × × × × × × × × × px The specific heat is now obtained as the derivative of
× × × × × × × × × × × × × × ×
× × × × × × × × × × × × × × ×
energy
× × × × × × × × × × × × × × ×
× × × × × × × × × × × × × × ×
- ∂E(T, V, N )
× × × × × × × × × × × × × × × 2πh C= . (21)
× × × × × × × × × × × × × × ×
L ∂T
× × × × × × × × × × × × × × ×
× × × × × × × × × × × × × × × In order eliminate µ appearing in the distribution function
(17) one has to simultaneously satisfy
The number of particles can be calculated as Z ∞
4 3
N
X
3 πpF
= g()n()d. (22)
N =2 1=2 , (15) V 0
p<pF
(2πh̄/L)3
The result calculated with Mathematica is shown below.
where the factor 2 comes from spin. From this we get a 2C
relation between the Fermi wave vector and the particle
3 N kB
density, 1.0
N p3F
= . (16)
V 3π 2 h̄3 0.8
3
(For detailed derivation see Ashcroft-Mermin, Solid state 3. Construction of the theory
physics.)
Landau’s idea
Liquid 3 He
The experiment above raises the following idea. Could it
Helium has two stable isotopes, 4 He and 3 He. The former be possible that low temperature liquid 3 He would effec-
is by far more common in naturally occurring helium. It is tively be like an ideal gas? This was the problem Landau
a boson since in the ground state both the two electrons started thinking. He had to answer the following questions
have total spin zero, and also the nuclear spin is zero. It has
a lot of interesting properties that could be discussed, but • How could dense helium atoms behave like an ideal
here we concentrate on the other isotope. 3 He is a fermion gas?
because the nuclear spin is one half, s = 1/2. Studies of 3 He
were started as it became available in larger quantities in • If there is explanation to the first question, how one
the nuclear age after world war II as the decay product of can understand the difference by 2.7 in the density of
tritium. The two isotopes of helium are the only substances states?
that remain liquid even at the absolute zero of temperature.
• If previous questions have positive answers, are any
other modifications needed compared to the ideal gas?
Weak interactions
As a first attempt to answer the questions, consider
point-like particles (instead of real 3 He atoms). In an ideal
gas the particles fly straight trajectories without ever col-
liding. If we now allow some small size for the particles,
they will collide with each other.
Figure: the specific heat of liquid 3 He at two different
a
densities [D.S. Greywall, Phys. Rev. B 27, 2747 (1983)]. U(r) c
We see that at low temperatures, the specific heat is li-
near in temperature. This resembles the ideal gas discussed
above, but is quite puzzling since the atoms in a liquid are d
more like hard balls continuously touching each other! 0 b
0 r
Figure: illustration of various particle-particle potentials
U (r): (a) the potential between two 3 He atoms, (b) ideal
gas potential U ≡ 0, (c) potential used in scattering ap-
proach, (d) potential used in perturbation theory.
We have to consider two particles with momenta p1 and
p2 colliding and leaving with momenta p01 and p02 . In such
a process the momentum and energy has to be conserved,
2π X
Γ= |hf |Hint |ii|2 δ(Ef − Ei ). (27)
h̄
f
4
We see that the rate is proportional to the number of avai- (29). Such excitations are called quasiparticles: they deve-
lable final states f . Consider specifically the case of filled lop continuously from single-particle excitations when the
Fermi sphere plus one particle at energy 1 > F . We wish interactions are ”turned on”, but they consist of correlated
to estimate the allowed final states when particle 1 collides motion of the whole liquid. ii) The quasiparticles have long
with any particle inside the Fermi sphere, 2 < F . The life time at low energy, like in the scattering approximation
final state has to have two particles outside the Fermi sp- above.
here (01 > F , 02 > F ) since the Pauli principle forbids all
It should be noticed that Landau’s theory is phenome-
states inside. We see that the number final states gets very
nological. At this stage it has one parameter, m∗ , whose
small when the initial particle is close to the Fermi energy, value is unknown theoretically, but can be obtained from
namely both 2 and 01 have to be chosen in an energy shell experiments.
of thickness ∝ 1 − F . This means that the final states are
limited by factor ∝ (1 − F )2 . Thus the scattering of low Although the detailed structure of the quasiparticle re-
energy particles is indeed suppressed and thus resembles mained undetermined, we can develop a qualitative pic-
the one in an ideal gas. ture with a model. Consider a spherical object moving in
otherwise stationary liquid. The details of this model are
But 3 He atoms are not point particles, rather they touch discussed in the appendix. The main result is that associa-
each other continuously. Thus for one particle to move,
ted with the moving object, there is momentum in the fluid
the others must give the way. This is one of the hardest
in the same direction. In the literature this is sometimes
problems in many body theory even today, but one can called ”back flow”, but I find this name misleading. Rat-
get some idea of what happens with a model: instead of her it should be called ”forward flow”or that the moving
true 3 He-3 He interaction potential, one assumes a weak object drags with itself part of the surrounding fluid. Thin-
potential, whose effect can be calculated using quantum- king now that the total momentum of the quasiparticle is
mechanical perturbation theory. We will skip this calcula-
fixed, this means that switching on the interactions slows
tion here (see Landau-Lifshitz). The result is that the exci-
the original fermion down, since part of the momentum
tation spectrum remains qualitatively similar as in free Fer- goes into the surrounding fluid and less is left for the ori-
mi gas but there is a shift in energies. Consider specifically ginal fermion.
the case, already mentioned above, of a filled Fermi sphe-
re plus one particle at momentum p, with p > pF . This
excited state of the ideal gas corresponds to the excitation
energy
p2 pF
p − F = − F ≈ (p − pF ), (28)
2m m
where the approximation is good if p is not far from the
Fermi surface (p − pF pF ). The effect of the weak inte-
ractions is now that the excitation energy still is linear in
p−pF , but the coefficient is no more pF /m. It is customary
to write the new excitation energy in the form The same picture is obtained by quantum mechanical
analysis. In order to get the velocity of the quasiparticle,
pF we have to form a localized wave packet. This travels with
p − F = ∗ (p − pF ), (29)
m the group velocity. Based on the dispersion relation (29)
the group velocity is
where we have defined the effective mass m∗ . Note that
the Fermi momentum pF is not changed, equation (16) dEp pF
vgroup = = ∗. (31)
still remains valid. With the new dispersion relation (29) dp m
we get a new density of states
This means that the momentum of the original fermion in
m∗ pF the interacting system is mvgroup p̂ = (m/m∗ )p, i.e. the
g(F ) = 2 3 . (30)
π h̄ original fermion contributes fraction m/m∗ of the momen-
tum p and the fraction 1 − m/m∗ is contributed by other
This is determined by the effective mass m∗ , not the ba- fermions surrounding the original one. The velocity of the
re particle mass m as for ideal gas (24). We now see that quasiparticle (31) is known as the Fermi velocity
weak interactions can explain that the specific heat coef-
ficient (23) differs from its ideal gas value. However, the pF
vF = ∗ . (32)
theory is valid for small perturbations, say 10%, and thus m
is insufficient to explain the factor 2.7.
Quasiparticle interactions
Quasiparticles
Thus far we have arrived at the picture that the low
Landau now made the following assumptions. i) Even for energy properties of a Fermi liquid can be understood as
strong interactions, the excitation spectrum remains as in an ideal gas with the difference that the effective mass m∗
5
appears instead of the particle mass m. In the following conventional to define F (p̂ · p̂0 ) = g(F )f (p̂ · p̂0 ). Such a
we show that this cannot be the whole story, and one more function can be expanded in Legendre polynomials
ingredient has to be added in order to arrive at a consistent ∞
theory. X
F (p̂ · p̂0 ) = Fls Pl (p̂ · p̂0 ), (36)
A general requirement of any physical theory is that the l=0
predictions of the theory should be independent of the coor- 2
dinate system chosen. In the present case, one has to pay where P0 (x) = 1, P1 (x) = x, P2 (x) = (3x − 1)/2, etc.
attention to Galilean invariance. That means that the phy- Using these we can now reduce the requirement (35) to
sics should be the same in two coordinate frames that mo-
m∗ Fs
ve at constant velocity with respect to each other. To be =1+ 1. (37)
m 3
specific consider a coordinate system O, and a second coor-
dinate system O0 that moves with velocity u as seen in the
We have arrived at the result that in order to have
frame O. We assume to study a system of N particles (in-
m∗ 6= m, we also should include an interaction between
teracting or not) of mass m. If the total momentum of this
the quasiparticles of the form of the F1s term in (36). The
system in O0 is P 0 , then the momentum seen in frame O
other interaction terms with coefficients Fls are not requi-
has to be P = P 0 + N mu. Now the Galilean invariance
red for internal consistency of the theory but some of them
requires that if one determines the state of the system in
appear as a result of perturbation theory. In order to make
O0 at fixed total momentum P 0 , it is the same as one would
general phenomenology, they all should be retained.
do in O with momentum P .
Until now we have not considered the fermion spin except
The ideal gas obviously satisfies Galilean invariance.
∗ that they produced factors of two. In magnetic field the sys-
However, when we replace the particle mass m by m in
tem becomes spin polarized, and this is no more sufficient.
(29), the Galilean invariance is broken. The cure for this
In the general case we have to consider the spin as a fully
problem is that we have to allow interactions between the
quantum object. This means that we have to replace np
quasiparticles. Thus we rewrite (29) into the form
by a 2 × 2 density matrix in spin space. All the previous
pF
p − F = ∗ (p − pF ) + δp , (33) analysis can be generalized to include the spin dependence.
m This means that also the quasiparticle energy p becomes
1 XX 0 (0) a 2 × 2 matrix. Equation (33) has to be generalized to
δp = f (p, p )(np0 − np0 ). (34)
V σ 0
p pF
ην (p) − F δην = ∗ (p − pF )δην
(0) m
Here np is the distribution of the quasiparticles, np = 1 XX
Θ(pF − p) is the distribution function in the ground sta- + fηα,νβ (p, p0 )[nαβ (p0 ) − n(0) (p0 )δαβ ], (38)
V 0
te, where Θ(x) is the step function (Θ(x) = 0 for x < 0 p α,β
6
4. Equation of motion where Π is the momentum flux tensor. The expressions for
such fluxes can be obtained as follows. One takes the ti-
The central quantity in the theory is the quasiparticle
distribution function np (r, t). Its equation of motion is de- me derivative of the corresponding density, for example ρ
rived similarly as the Boltzmann equation, by calculating (44), one applies the kinetic equation (43) to the partial ti-
the total time derivative me derivative. Because of conservation (42), the contribu-
tion from collision integral vanishes. Thus one is left with
dnp ∂np ∂np the divergence term in (43), which allows to identify the
= + ṙ · ∇np + ṗ ·
dt ∂t ∂p expression for the flux, J in case of ρ. In the case of ρ, we
∂np ∂p ∂np can verify the consistency of equations (44)-(46).
= + · ∇np − ∇p · , (40)
∂t ∂p ∂p Above we have considered the case of no spin-
dependence. The general case is more complicated since
where ṙ ≡ dr/dt = ∂p /∂p and ṗ ≡ dp/dt = −∇p . one has to use quantum equation of motion (von Neumann
Equating this with the rate of change caused by collisions equation) for the spin density matrix. This allows to obtain
Ip gives the Landau-Boltzmann equation a kinetic equation and spin dependent observables. We give
∂np ∂p ∂np as examples the spin density and the spin current density
+ · ∇np − ∇p · = Ip . (41) tensor and the conservation law for spin,
∂t ∂p ∂p
h̄ 1 X
This differs from the ordinary Boltzmann equation that S= Tr[σ(np − n(0)
p )], (48)
2V p
the energy p (r, t) contains, in addition to external poten-
tials, the interaction energy δp (34). The collisions have spin h̄ 1 a 1 X
J = 1 + F p Tr[σ(np − n(0)
p )], (49)
to conserve the fermion number, momentum and energy. 2m∗ 3 1 V p
This implies on the collision term the conditions
∂S
X X X + ∇ · J spin = 0. (50)
Ip = 0, pIp = 0, p Ip = 0. (42) ∂t
p p p
Here Tr(. . .) denotes the trace of the 2 × 2 matrix and
σ = x̂σx + ŷσy + ẑσz is the vector of Pauli matrices. In the
The kinetic equation equation (41) is a nonlinear equa- spin-dependent case we have to generalize (44) and (45) to
tion. For most, if not all, applications a linearized form of
it is sufficient. In equilibrium both ∇np and ∇p vanish, 1 X
ρ = ρ0 + m Tr(np − n(0)
p ), (51)
and thus in a linearized theory their multipliers in (41) can V p
be evaluated in the equilibrium state. Thus linearization of 1 X
the left hand side gives J= p Tr(np − n(0)
p ). (52)
V p
!
(0)
∂np dnp
+ vF p̂ · ∇ np − δp̂ = Ip , (43)
∂t dp
Properties of a quasiparticle
(0)
where dnp /dp denotes the derivative of the equilibrium In order to develop intuition into the formulas above,
distribution function (17) evaluated at the unperturbed let us consider a single localized quasiparticle. It could be
energy (29). defined by distribution np (r, t) that differs from the ground
(0)
state distribution np only for momentum values p around
Once we have solved for the distribution function
some p0 , where p0 ≈ pF . We normalize the distribution by
np (r, t), we can calculate measurable quantities. For
example, the mass density and momentum density of the
Z
1 X
fluid are given by d3 r Tr(np − n(0)
p ) = 1, (53)
V p
Z
2 X 1 X
ρ = ρ0 + m (np − n(0)
p ), (44) d3 r Tr[σ(np − n(0)
p )] = ẑ, (54)
V p V p
2 X
J= p (np − n(0)
p ). (45) where we have chosen the spin polarization in the z di-
V p
rection. Applying these to equations (51) and (48) gives
that the quasiparticle has mass m and spin 21 h̄ẑ. This is
Here the factors of 2 come from spin. The momentum den- the same as for an additional particle in a noninteracting
sity is the same as the mass current density. Because of system. The quasiparticle propagates at the Fermi veloci-
conservation of mass and momentum, these have to obey ty (32), v = v p̂ . From (52) we get that the momentum
F 0
conservation laws of the excitation is p0 . As discussed above, the additional
∂ρ fermion now contributes only mvF to the momentum and
+ ∇ · J = 0, (46) the rest of the momentum (m∗ − m)vF = 1 F1s mvF comes
∂t 3
∂J from forward flow of other fermions. Now we can state also
+ ∇ · Π = 0, (47) the spin-content of this cloud. Namely, equations (52) and
∂t
7
(49) together imply that the effective number of spin up 5. Applications
and spin down particles in the quasiparticle are
The Fermi-liquid theory can be applied to calculate se-
1 1 veral measurable properties of a Fermi liquid. Below we list
n↑ = 1 + (F1s + F1a ), n↓ = (F1s − F1a ). (55) some of them. The standard applications are only briefly
6 6
mentioned as they can be found in most reviews of the
That is, the total mass current in the quasiparticle carried Fermi-liquid theory [see references mentioned in the intro-
by spin-up particles is n↑ mvF and correspondingly n↓ mvF duction].
is carried by spin-down particles.
Specific heat. For a Fermi gas this was calculated in Eq.
(23):
π2 2
C= g(F )kB T + O(T 2 ). (56)
3
g(F )
χ = µ0 µ2m , (57)
1 + F0a
r
1 1
c = vF (1 + F0s )(1 + F1s ). (58)
3 3
8
In the temperature range immediately above T
VGLUMK 17, NUMBER 2 PHYSICAL REVIEW LETTERS
frequencies are nearly constant,
11 JULY 1966 while the peak hei
I
'
I
'I I I ltTI I
'
I
'I I I IIII Table I. Values of the coefficients &0 and A~ of Eqs.
(1) and (2) for severali runs. i i i
400—
300— Cd /27l Ao A(
Run (MHz) {106 (K 2
cm) ] [].0 T / T cK =/c m]
sec
200
E
15.4 1.44 2. 74
I-
C3
I
15.4 1.57 2. 65
à so
IOO
45. 5 l. 58 2. 66
gb 15.4 1.62 2. 65
60
O
O
z0 40
3O
LLI
shows that, in the case of attenuation measure-
20
I—
ments, in which Tmax* is well determined,
tLI
C3 the frequency dependence is quantitatively veri-
lo
fied. It is difficult to estimate Tm~ from
8
the velocity measurements, so the scatter is
greater. However, it seems clear from these
I I I I II '
I
data that the transition observed in Ref. 3 is
T —MILLIDEGREES KELVIN indeed to be attributed to a transition from first
196—
to zero sound.
C3
I94- ~0Cn t9~
~ go
00 Our results can also be compared with the
l92—
0 0 theory of Khalatnikov and Abrikosov' in which
0 oo Fig. 3. Temperature depende
0
the velocity of zero sound is found to be given
I
I—
O 190— 0 the normal phase at zero pres
O 0
by the implicit equation
I I I
0
o
o o 0
of 2 MHz. Temperatures are
I88— 1 + ,'E, (1 kHz/div)
Frequency —
0
08 0DCO 0 O OD
transition temperature Tc ~ 1
I86 I I iI I I IIII I I I I I IIII I Eo(l + —,'E, ) + so'E, '
FIG. 2. Amplitude attenuation coefficient and sound -l)]f-l,
[W.R. Abel,
propagation A.C.as Anderson
velocity a function ofand Wheatley, where
J.C. tempera-
magnetic Phys. zo(s, )Kinetic
=1(s,/2) coefficients:
ln[(s, +1)/(s, viscosity,spthermal conductivity, spin
=cp/vF, and vF is The
diffusion. the Fermi method Recent-
generalvelocity. to calculate the thermal con-
ture Lett.
Rev.
is
in pure17,
line its
from
74 He3
liquid
waves. ofConsider
the average
drawn
at 0.32 atm and for frequencies
(1966).]
of 15.4 and 45. 5 MHz. Each point shown on the graph
Spin several raw a data
localpoints.
low-temperature
through the value.
equilibrium
deviation
In a Fermi
ly determined'
of spin density
The straight
attenuation
gas this m/sec
da- relaxes by forcollision
values
ductivity,
eters areequation 10 77
forofexample,
(41). For
the Fermi-liquid
6 25 and v
is to solve
that one has
a pressure of 0. 28 atm (close to the
param-
53 to
the Landau-Boltzmann
8 take a closer look at the Eo: +, : F:
term Ip . A general form for it is
ta represents Eq. (1), while the straight lines drawn
diffusion of the particles, called spin diffusion. In a present one).
Fermi At this pressure c, =187.2 m/sec.
through the high-temperature data represent Eq. (2), Solving Eq. (3) for s, one finds s, = 3.597 [and 0 0
X
liquid the interactions
with u/27r equal to 15.4 and intervene
45. 5 MHz. so that
With instead
the presentof diffusion, Ip = −
M (s, ) = 0. 027 033]. This leads to [(cp-cl )/
w(p, p1 ; p , p1 )
gap it was notwaves
propagating possibleof tospin
measure the 45.
density can5-MHz atten-
be generated. One 0 ,p0
uation coefficient above e =200 cm . The smooth cl]p 28 atm 0. 034, in remarkable
—— p1 ,p agreement
1
particularly interesting case is the following.
curve just above the attenuation data for 15.5 MHz is a
Consider sta-
with the measured value at
×[n 0. 32
n
p p1 atm
(1 −of n0.
p0035
)(1 − np01 )
ticplot
magnetic
of Eq. (6)field H ==15.
with cd/27I H 4ẑ MHz
andand theeomagnetization
and 0. given + 0.003. Using a theory for energy transfer
vector
−(1 − n )(1 − n )n 0 n 0 ]. (60)&
9
where nF = ρ/m is the number density of fermions. qu
pb
as
ipa
Acoustic impedance. The analytical methods to solve the r tic
le quasiparticle trajectory
Landau-Boltzmann equation were developed to high sop-
histication in 1960’s. Besides the exact calculation of the be s
p am
kinetic coefficients, another such problem was the acoustic
impedance first solved by I. L. Bekarevich and I. M. Kha-
latnikov [Sov. Phys. JETP 12, 1187 (1961)] and further
developed by E. G. Flowers and R. W. Richardson [Phys. In a Fermi gas, a quasiparticle on the trajectory would
Rev. B 17, 1238 (1978)]. Consider a semi-infinite Fermi not react to the beam. In a Fermi liquid it experiences the
liquid bounded by a planar wall. Suppose small oscilla- potential change δ (34) caused by the beam. In the simple
tions of the wall (either transverse or perpendicular to the case that F0s is the only relevant Fermi-liquid parameter,
wall). The acoustic impedance Z is defined as the ratio of the potential δ = [F0 /2N (0)]δn is determined by the par-
the force on the liquid F divided by the velocity of the wall ticle density in the beam.
u,
ε ε ε
µ pF
F = Zu, (63)
εF
δε
With harmonic time dependence ∝ exp(−iωt), the impe-
dance is complex valued, Z = Z 0 + iZ 00 . Here Z 0 gives the p p p s
dissipation and Z 00 the mass of the fluid coupled to oscil-
lation. A basic assumption of the Fermi-liquid theory is that the
potential is small compared to the Fermi energy F , i.e.
δ F . This means that a quasiparticle is slightly dece-
lerated when it enters the beam and it is accelerated back
when it leaves the crossing region. In the energy point of
view (figure above), the quasiparticle flies at constant ener-
gy ≈ µ and the potential is effectively compensated by
depletion of fermions with the same momentum direction
in the crossing region.
Let us now consider the case that the intensity of the
quasiparticle beam is varying in time. This means that
the potential seen on the crossing trajectory changes. This
changes the number of particles stored in the crossing re-
gion, and thus leads to emission of particle or hole like qua-
siparticles from the crossing region. This takes place on all
Experimental data from P.R. Roach and J.B. Ketterson crossing trajectories.
[Phys. Rev. Lett. 36, 736 (1976)]. Figure: E. G. Flowers, The effect is similar as sailing in an ocean. Usually one
R. W. Richardson and S.J. Willamson [Phys. Rev. Lett 37, is not interested how deep the water is, as long as it is
309 (1976)]. sufficiently deep. The situation is different if the ocean floor
is changing in time: one creates a tsunami.
10
generates a beam of quasiparticles. This beam interacts 6. Generalizations
with the quasiparticles that are coming to the wire. In case
Landau’s theory was originally proposed of a pure, non-
of negative F0 , which is the case in 3 He-4 He mixture, this superfluid Fermi liquid. Various generalizations are pos-
increases the restoring force. Thus the oscillation frequency sible.
is increased by the effect of Fermi-liquid interactions. In
order to determine this quantitatively, one has to solve the The theory can be generalized to the presence of conden-
Landau-Boltzmann equation (43). sed bosons [I. M. Khalatnikov, Sov. Phys. JETP 28, 1014
(1969); E.V. Thuneberg and T.H. Virtanen, Phys. Rev. B
dissipation 83, 245137 (2011)]. The main application of this theory is
Z' 2 4
1 3 to liquid mixtures of 4 He and 3 He. However, it can also
3 an p
F F clarify some conceptual points of Landau’s theory. In par-
quid
1 decoupling of fermions ticular, the relation (37) between m∗ and F1s makes some
te l i
infin 2 decoupling of bosons formulas of the theory to appear as trivial. In fermion-
2 ll t bound to quasiparticles
c e e n
3 effect of interactions, boson mixtures (37) has to be generalized and thus may
a im
in p er "Landau force" lead to better appreciation of the theory.
u ide x
liq g 4 reflected quasiparticles
in re
1 as tu from cell wallsThe Fermi-liquid theory can be generalized to the su-
cre era
de mp perfluid phase by taking into account the pairing of fer-
hydr te
odyn mions to form Cooper pairs. The theoretical basis of the re-
amic
limit Z''-Z''ideal
0 sulting Fermi-liquid theory of superfluidity has been discus-
-0.3 -0.2 -0.1 0 0.1 0.2 anF pF
sed by J.W. Serene and D. Rainer [Phys. Rep. 101, 221
frequency shift wire radius (1983)]. This applies in particular to the superfluid pha-
ses of liquid 3 He. The theory and its results form a vast
Experimental data: J. Martikainen, J. Tuoriniemi, T. field, too large to be discussed here. We simply state that
Knuuttila, and G. Pickett, J. Low Temp. Phys. 126, 139 essentially all quantities of the superfluid state are affected
(2002), theory: T. Virtanen and E.T., Phys. Rev. Lett. 106, by the Fermi-liquid interaction parameters. Let us mention
055301 (2011). here our recent work, where the Fermi-liquid interactions
Effect of Fermi-liquid interactions in semi-infinite fluid: lead to a qualitative change of quasiparticle states in a vor-
see acoustic impedance (above). In Fermi liquid film: J. tex core [M. A. Silaev, E. V. Thuneberg and M. Fogelström,
Kuorelahti and E.T., poster in QFS2015. arXiv:1505.02136].
The Fermi-liquid theory can also be generalized other
fermion systems than 3 He. The most notable case is the
conduction electrons of metals. It applies to both the nor-
mal conducting and the superconducting phases. In metals
the presence of the crystal lattice leads to complications
and additional features compared to a translationally in-
variant Fermi system. Another important difference to unc-
harged system is the electric charge of the electron.
All the above considered three-dimensional Fermi liquid.
Landau’s theory can also be applied to interacting fermions
limited to move in two dimensions. I one dimensions a dif-
ferent behavior is expected.
11
7. Derivation and limitations 8. Conclusion
Landau developed the Fermi-liquid theory intuitively as In this introduction to the Landau’s Fermi-liquid theory
the low energy expansion of a Fermi system. What it means we have concentrated on its starting point. Especially we
that the Fermi-liquid state will appear if nothing else hap- have tried to clarify the nature of the quasiparticle. Out
pens. The possible alternative is that if the interactions are of the many applications of the theory to 3 He and elsew-
strong enough, the particles will be localized and the sys- here we have only briefly mentioned a few. We hope that
tem becomes solid instead of a liquid state. An interesting understanding the elements of the theory gives good star-
approach to this is presented by R. Shankar [Rev. Mod. ting point to study the known applications, or finding new
Phys. 66, 129 (1994)]. He uses the renormalization-group ones.
method and finds that zooming into the neighborhood of
the Fermi surface, the Fermi-liquid model appears as one
alternative possibility.
The Fermi-liquid theory should appear as the low energy
limit of a microscopic many-body theory. This means that
one should be able to calculate the Fermi-liquid parame-
ters starting from a microscopic theory. At the moment it
seems that the quantum-fluids community is using experi-
mental values of the parameters. My hope is that some day
the many-body community would calculate more accurate
values of the parameters.
12
Appendix: hydrodynamic model of a quasipar-
ticle
The starting point is Euler’s equation and the equation
of continuity
∂v 1
+ v · ∇v = − ∇p
∂t ρ
∂ρ
+ ∇ · (ρv) = 0, (64)
∂t
where p is now the pressure. We assume a small veloci-
ty so that the nonlinear term can be dropped and assume
incompressible fluid, ρ = constant. The boundary condi-
tion on the surface of a sphere of radius a is n · v = n · u,
where u is the velocity of the sphere and n the surface
normal. We assume that far from the sphere v → 0. The
velocity can be represented using the potential v = ∇χ
where χ = −a3 u · r/2r3 . By the Euler equation the pres-
sure p = −ρχ̇, and the force exerted by the sphere on the
fluid
2πa3 ρ
Z
F = da p = u̇ (65)
3
is proportional to the acceleration u̇. Therefore, associated
with a moving sphere, there is momentum in the fluid
2πa3 ρ
p= u (66)
3
corresponding to half of the fluid displaced by the sphere
and moving in the same direction and at the same velocity
as the sphere. Note that this may differ from the total
momentum of the fluid, which is not essential here, and is
undetermined in the present limit of unlimited fluid.
13