Spin degrees of freedom and flattening of the spectra of singleparticle
excitations in strongly correlated Fermi systems
Date
20030815
Description
Comment: 7 pages, 15 figures
The impact of longrange spinspin correlations on the structure of a flat
portion in singleparticle spectra $\xi(p)$, which emerges beyond the point,
where the Landau state loses its stability, is studied. We supplement the
wellknown Nozieres model of a Fermi system with limited scalar longrange
forces by a similar longrange spindependent term and calculate the spectra
versus its strength $g$. It is found that Nozieres results hold as long as
$g>0$. However, with $g$ changing its sign, the spontaneous magnetization is
shown to arise at any nonzero $g$. The increase of the strength $g$ is
demonstrated to result in shrinkage of the domain in momentum space, occupied
by the flat portion of $\xi(p)$, and, eventually, in its vanishing.
portion in singleparticle spectra $\xi(p)$, which emerges beyond the point,
where the Landau state loses its stability, is studied. We supplement the
wellknown Nozieres model of a Fermi system with limited scalar longrange
forces by a similar longrange spindependent term and calculate the spectra
versus its strength $g$. It is found that Nozieres results hold as long as
$g>0$. However, with $g$ changing its sign, the spontaneous magnetization is
shown to arise at any nonzero $g$. The increase of the strength $g$ is
demonstrated to result in shrinkage of the domain in momentum space, occupied
by the flat portion of $\xi(p)$, and, eventually, in its vanishing.
Type
Identifier
doi:10.1134/1.1619498
Phys. At. Nucl. 66, 1871 (2003)
Database
Link to record
Show preview
Hide preview
ar
X
iv :c
on d
m at
/0 30
83 09
v1 [
co nd
m at.
str el
] 15
A ug
20 03
Spin degrees of freedom and flattening of the spectra of singleparticle
excitations in strongly correlated Fermi systems
V. A. Khodel1, P. Schuck2, and M. V. Zverev1 1 Kurchatov Institute, Russian Research Center,
Moscow, 123182 Russia 2 Groupe de Physique The´orique,
Institut de Physique Nucle´aire,
F01406 Orsay Cedex, France
(Dated: August 13, 2013)
The impact of longrange spinspin correlations on the structure of a flat portion in singleparticle spectra ξ(p), which emerges beyond the point, where the Landau state loses its stability, is studied. We supplement the wellknown Nozieres model of a Fermi system with limited scalar longrange forces by a similar longrange spindependent term and calculate the spectra versus its strength g. It is found that Nozieres results hold as long as g > 0. However, with g changing its sign, the spontaneous magnetization is shown to arise at any nonzero g. The increase of the strength g is demonstrated to result in shrinkage of the domain in momentum space, occupied by the flat portion of ξ(p), and, eventually, in its vanishing.
The investigation of flattening of singleparticle (sp) spectra of Fermi liquids is dated back to Ref. [1], where longrange correlations, enhanced in the vicinity of an impending ferromagnetic phase transition, were shown to result in the divergence of the effective mass M∗ at the transition point. Later in [2] an idea of the socalled fermion condensation, i.e. a rearrangement of the Lan dau state, occurring beyond a critical point in strongly correlated Fermi systems with longrange effective inter actions, was suggested. A striking feature of this rear rangement is ”swelling” of the Fermi surface (FS), i.e. the occurrence of a completely flat portion ξ(p) = 0, called the fermion condensate (FC), in the spectrum ξ(p), mea sured from the FS. To gain insight into the problem of fermion condensa
tion, let us turn to the Dyson equation, rewriting it in the form
ξ(p) = ξ0 p +Σ(p, ξ(p)) , (1)
appropriate for finding the FC solutions ξ(p) = 0. The sp mass operator Σ is usually determined by the formula
Σ(p, ε) =
∫ Wp, ε,p1, ε1)G(p1, ε1)
d4p1 (2pi)4i
, (2)
whereW is an effective interaction between particles, and G is the sp Green function. It is worth noting that the imaginary part of Σ(p, ε) vanishes at ε = 0. Therefore, in searching for the FC solutions, only the real part of Σ is relevant. For a long time, it was reckoned that in homogeneous
Fermi systems, there exists a onetoone correspondence between the momentum p and the sp energy ξ, at least, close to the FS, i.e. the derivative (dξ/dp)F is always positive, a postulate, virtually being a cornerstone of the Landau theory of Fermi liquid [3]. However, in sys tems with longrange forces it might be incorrect. Com pelling evidence for that is provided by a phenomenologi cal model, suggested by Nozieres [4], whereW is taken as
0(ξξp)/f
0
1
n(ξ)
ξ
FIG. 1: Graphical illustration of the solution of Eq. (3).
a constant in the coordinate space. Thus, in this model Eq. (1) takes the form
ξ − ξ0 p = fn(ξ) , (3)
where f is an effective coupling constant, the value of which is taken to be positive, while n(ξ) is the Landau quasiparticle momentum distribution, being 0 at positive ξ, and 1, at negative ones. This equation is easily solved, but we concentrate on its graphical solution (see Fig. 1). Let us draw both the r.h.s. and the l.h.s. of Eq. (3) as functions of ξ, taking p as an input parameter. The l.h.s. of (3), depicted at different p, provides a set of par allel straight lines, while the r.h.s. forms a kink, the ver tical segment of which is located at ξ = 0, no matter what is the input. Crossing points yield the sp spectrum ξ(p). If they lie on any of the two horizontal pieces of the kink, ξ(p) does coincide with the sp spectrum of ideal Fermi gas. An unconventional situation occurs when straight lines cross the steep section of the kink. In this case, the intersection point remains 0 in a finite momentum inter val [pi, pf ]. This plateau ξ(p) = 0, lying exactly at the FS, does form the fermion condensate. Since ξ(p) = 0 in the finite volume of momentum space, the density of states ρ(ε) acquires a singular term ρ(ε) ∼ δ(ε). Another salient feature of the fermion condensation phenomenon
2 is that in the FC domain, the quasiparticle momentum distribution n0(p) differs from the Landau one. Indeed, upon setting ξ = 0 in Eq.(3) one finds
n0(p) = − ξ0 p
f , pi < p < pf , (4)
in contrast to the conventional step function nF (p) = θ(pF − p). Thus, the FS does swell, whereas the basic assertion of Landau theory fails. In dealing with the FC problem, attention is usually
paid to the spinindependent part of the effective inter action W (see e.g. [2, 4, 5, 6, 7, 8]). Here we investigate effects, associated with longrange spinspin components of W , which involve the spinup quasiparticle distribu tion n+(p) and spindown one n−(p). It is instructive to start such an analysis with a generalization of the Nozieres model [4], supplementing it by longrange spin spin terms. As a result, one obtains two equations
ξ+(p) = ξ 0 p +
1
2 f (n+(p) + n−(p))
+ 1
2 g (n+(p)− n−(p)) ,
ξ−(p) = ξ 0 p +
1
2 f (n+(p) + n−(p))
− 1
2 g (n+(p)− n−(p)) , (5)
with a new constant g, specifying the longrange spin spin component of the model effective interaction W . It is worth noting that in the case of spin fluctuations with the nonzero critical momentum qc ≪ pF , presumably relevant to twodimensional liquid He3, the constants in Eq. (5) are related with each other: g = −f/3. In conventional Fermi liquids, the spontaneous spin
S arises only if strength of the spinspin interaction obeys the Pomeranchuk condition Gρ(0) > 1, where ρ(0) = pFM
∗/pi2 is the density of states at the FS in the Landau theory. In the case at issue, the density of states is infinite, so that one can expect the emergence of the spontaneous magnetization at any g < 0. To facilitate the solution of the problem, we recast the
system (5) to the form
ξ+ = ( 1−
a
b
) ξ0 p + a
b ξ− +
( b−
a2
b
) n− ,
ξ− = ( 1−
a
b
) ξ0 p + a
b ξ+ +
( b−
a2
b
) n+ , (6)
where a = (f + g)/2 and b = (f − g)/2. Let us now draw the plot ξ+(ξ−), proceeding from the first of Eqs.(6) and treating p as an input parameter. At ξ− > 0 the function ξ+(ξ−) is given by the straight line ξ+ = (1 − a/b)ξ
0 p +
aξ−/b, while at ξ− < 0 this straight line is slightly shifted: ξ+ = (1−a/b)ξ
0 p +(b2−a2)/b+aξ−/b. The drawing of its
counterpart ξ− from the second of Eqs.(6) yields another couple of straight lines: ξ− = (1 − a/b)ξ
0 p + aξ+/b at
ξ+ > 0, and ξ− = (1 − a/b)ξ 0 p − (b2 − a2)/b + aξ−/b at
ξ 
ξ 
ξ +
ξ +
ξ +
ξ 
(a)
A
ξ 
0
0
ξ +
(b)
0
B0
B1
'B1
ξ 
0
ξ +
0
C1
'C1 (c)
C0
0
0
D1
(d)
'D1
D0
0
E 0
(e) 0
FIG. 2: Graphical illustration of the dynamics of solutions of the set (6) with increasing ξ0p in the case of g < f .
ξ+ < 0. As we shall see, both these curves cross each other either at one or at three points, providing several solutions for the spectra ξ±(p).
Some elucidation of the procedure is ensured by Fig. 2, where the evolution of solutions of the set (6) versus the input parameter ξ0
p is shown. The solid line corresponds
to the dependence ξ+(ξ−), given by the first equation of the set (6), while the dashed one shows the same depen dence resulting from the second equation of the set. The intersection points are indicated by letters. Five panels of Fig. 2 show five different cases referring to different ξ0 p . In the first case (panel (a)), two zigzag lines have
one intersection point A, corresponding to the single so lution where both ξ+(p) and ξ−(p) are negative. With increasing ξ0
p , the intersection point moves to the origin.
At a certain value of this variable, depending on the pa rameter f only, bifurcation emerges, and there appear three intersection points (B0, B1 and B
′ 1 on the panel
(b)). One of them, B0, lies at the origin, while the other two, B1 and B
′ 1, on the ξ+ and ξ− axes, respectively. A
similar situation occurs in the panels (c) and (d). The intersection at the origin persists, while the other two crossing points lie inside the quadrants on the panel (c) (C1 and C
′ 1) and on the ξ− and ξ+ axes on the panel (d)
(D1 and D ′ 1). The last (e) panel shows the fifth case with
one intersection point at positive values of ξ+ and ξ−.
The chain A−B0−C0−D0−E is associated with the solution (denoted as Φ0), for which the sp spectra coin cide: ξ+(p) = ξ−(p) = 0 within the interval [pi, pf ], the length of which is determined by the parameter f only, as if there was no spinspin interaction at all. The analysis
3 0,1
0,0
0,1
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.1 f=0.3
(a) 
+
ξ +,
(p) p/pF
(b)
+n+,(p)
(c)
p/pF
s(p)
FIG. 3: Singleparticle spectra ξ±(p) in units of ε0F = p 2 F/2M
(panel (a)), occupation numbers n±(p) (panel (b)) and spon taneous spin s(p) = n−(p)− n+(p) (panel (c)), calculated for f = 0.3, g = −0.1 in units of ε0F .
0,1
0,0
0,1
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.01 f=0.3
(a)

+
ξ +,
(p) p/pF
(b)
 +
n +,
(p)
(c)
p/pF
s(p)
FIG. 4: The same as in Fig. 3 for f = 0.3, g = −0.01.
shows that, as long as the constant g remains positive, it is the solution Φ0 that has the lowest energy compared with the others. As for solutions, related to the chains
A−B1−C1−D1−E and A−B ′ 1−C
′ 1−D
′ 1−E (Φ1 and
Φ′1, respectively), the situation is different. The solution Φ1, corresponding to the first chain, is shown in Fig. 3. We see that the FC region, the boundaries of which were insensitive to g at any g > 0, is destroyed: the spectra
0,1
0,0
0,1
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.033 f=0.1
(a)

+
ξ +,
(p) p/pF
(b)
+ n
+, (p)
(c)
p/pF
s(p)
FIG. 5: The same as in Fig. 3 for f = 0.1, g = −0.033.
����
���
���
��� ��� ���
�
�
��� ��� ���
�
�
J ���� I ���
�D� �
�
ξ ���
�S�
S�S )
�
�
�E�
��
�� Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 6: The same as in Fig. 3 for f = 0.3, g = −0.3.
ξ+(p) and ξ−(p) repell each other. As the momentum p moves from the lower point, where these spectra simultaneously attain the FS, only the FC plateau at the spectrum ξ+(p) survives, with the energy splitting between the spectra ξ+(p) and ξ−(p) growing linearly with the p increase. Attaining the maximum at the point, where the plateau at ξ+(p) vanishes, this splitting ceases to increase and remains constant until the point, where a new FC plateau emerges at the spectrum ξ−(p), and, finally, both the spectra once again merge at the point, where the FC disappears forever.
4 ����
���
���
��� ��� ���
�
�
��� ��� ���
�
�
J ���� I ���
�D� �
�
ξ ���
�S�
S�S )
�
�
�E�
��
�� Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 7: The same as in Fig. 3 for f = 0.1, g = −0.1.
ξ 
ξ 
ξ +ξ+
ξ +
ξ 
(a)
A
ξ 
0
0
ξ +
(b)
0
B0 B2
'B2
B1
'B1
ξ 
0
ξ +
0
C2
'C2 (c)
C0
0
0
D2 D1
(d)
'D2 'D1
D0
0
E
0
(e) 0
FIG. 8: The same as in Fig. 2 for the case of g > f .
The second chain refers to the solution Φ′1, symmetric to the previous one, with changing± by ∓ and, thus, cor responds to the opposite sign of the projection of sponta neous spin S on the fixed axis. Since both the directions are equivalent, two solutions with the nonzero S have equal energies. It can be verified that at any g < 0, these solutions provide the minimum of the ground state en ergy, while the solution Φ0 with S = 0, the maximum. Indeed, calculations yield the gain in energy of the so
0,3
0,0
0,3
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.15 f=0.1
(a)

+
ξ +,
(p) p/pF
(b)
+n+,(p)
(c)
p/pF
s(p)
FIG. 9: The same as in Fig. 3 for f = 0.1, g = −0.15.
0,3
0,0
0,3
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.45 f=0.3
(a) 
+
ξ +,
(p) p/pF
(b)
+n+,(p)
(c)
p/pF
s(p)
FIG. 10: The same as in Fig. 3 for f = 0.3, g = −0.45.
lutions Φ1 and Φ ′ 1, as compared with the energy of the
Landau state E1/EL−1 ≃ −0.027, while for the solution Φ0, one obtains E0/EL − 1 ≃ −0.019. Thus, the reason for these alterations of the FC is the
occurrence of the spontaneous spin
S = ∑
p
s(p) = ∑
p
( n−(p)− n+(p)
) . (7)
arising at any g < 0. For illustration, the spin density s(p) is drawn in panel (c). We infer that the spin density s(p) differs from 0 only in the region [pi, pf ]. Therefore
5 0,3
0,0
0,3
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.45 f=0.1
(a) 
+
ξ +,
(p) p/pF
(b)
+ n
+, (p)
(c)
p/pF
s(p)
FIG. 11: The same as in Fig. 3 for f = 0.1, g = −0.45.
the total spin S depends primarily on the value of pa rameter f . For the parameters f and g, given in Fig. 3, one obtains S/ρ ≃ 0.299.
For comparison, results for the other two sets of the parameters: f = 0.3, g = −0.01 and f = 0.1, g = −0.033 are drawn in Figs. 4 and 5, respectively. The striking feature of these results is that for the given value f = 0.3 the magnitude of the spontaneous spin S for g = −0.01, S/ρ ≃ 0.232, is of the same order as in the case of g = −0.1. The energy gain of the solution with the spontaneous spin, corresponding to f and g for Fig. 4, E1/EL − 1 ≃ −0.019. Fig. 5 shows that multiplication of both the parameters f and g by the factor 1/3 results in the scaling of the spectra, occupation numbers, and spontaneous spin by the same factor.
The situation drastically changes if the magnitude of the spinspin force becomes comparable to that of the scalar one. This is seen in Figs. 6 and 7 show ing results obtained for the sets f = 0.3, g = −0.3 and f = 0.1, g = −0.1. We see that the FC prac tically disappears, while the Landau state reappears. Now the occupation numbers n+(p) and n−(p) have the Migdal jump from 0 to 1 at the points pi and pf , re spectively, and the spontaneous spin attains the maxi mum value, which is equal to the phase volume of the spherical layer between pi and pf (S/ρ ≃ 0.448 and S/ρ ≃ 0.149, for the parameters of Figs. 6 and 7, re spectively). Fig. 8 illustrates the evolution of solutions for the set (6) versus ξ0
p in the case when the abso
lute value of the spinspin constant g is larger than the scalar one f . In this case, the set (6) has five solutions: Φ0, corresponding to the chain A−B0−C0−D0−E; Φ1 and Φ′1, corresponding to the chains A−B1−C0−D1−E and A−B′1−C0−D
′ 1−E, respectively; and a new couple
A
0 ξ 
(a)
0
ξ +
(b)
B
0
0
ξ 
ξ +
(c)
C
0
0
ξ 
ξ +
'D1 D1
D0
(d)
0
0 ξ 
ξ +
'E1
E1
(e)
0
0
ξ 
ξ +
(f) 'F1
F0 F1
0
0
ξ 
ξ +
(g)
G
0
0
ξ 
ξ +
(h)
H
0
0
ξ 
ξ +
(i)
I
0
0
ξ 
ξ +
FIG. 12: Graphical illustration of the dynamics of solutions of the set (9) with increasing ξ0p in the case of g < f , βH < g/2.
of solutions, Φ2 and Φ ′ 2, corresponding to the chains
A−B2−C2−D2−E and A−B ′ 2−C
′ 2−D
′ 2−E. The calcu
lation of the energies yields the solutions Φ2 and Φ ′ 2,
equal in energy, beating Φ0, Φ1 and Φ ′ 1. The solution
Φ2 is drawn in Figs. 9, 10, 11 for three sets of parame ters: f = 0.1, g = −0.15; f = 0.3, g = −0.45; f = 0.1, g = −0.45. These figures show that in this case the spon taneous spin is defined only by the parameter g, and the parameter f defines the magnitudes of the jumps in the spectra ξ+(p) and ξ−(p). It is worth noting that, while the flat portions of the spectra, corresponding to FC, disappear and the Landau state is recovered, the phan tom of the FC manifests itself in the emergence of the spontaneous spin.
The above analysis can be generalized to the case of an external magnetic field. In this case, the equations for the sp spectra read
ξ+(p) = ξ 0 p +
1
2 f(n+(p) + n−(p))
+ 1
2 g(n+(p)− n−(p)) + βH ,
ξ−(p) = ξ 0 p +
1
2 f(n+(p) + n−(p))
− 1
2 g(n+(p)− n−(p))− βH , (8)
where H is the effective magnetic field acting on the fermion spin and β is the magnetic moment of the fermion. Upon rewriting Eqs. (8) to the form convenient
6 ξ 
ξ 
ξ +
ξ +
ξ +
ξ 
(a)
A
ξ 
0
0
ξ +
(b)
0
B
ξ 
0
ξ +
0
C
(c)
0
0
D
(d) 0
E 0
(e) 0
FIG. 13: The same as in Fig. 12 for the case of βH > g/2.
����
���
���
��� ��� ���
�
�
��� ��� ���
�
�
β+ ���� J ���� I ���
�D� �
�ξ ���
�S�
S�S )
�
�
�E�
��
��Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 14: The same as in Fig. 3 for f = 0.3, g = −0.1, βH = 0.04.
for the graphical analysis, we obtain
ξ+ = ( 1−
a
b
) ξ0 p + a
b ξ− +
( b−
a2
b
) n− +
( 1+
a
b
) βH,
ξ− = ( 1−
a
b
) ξ0 p + a
b ξ+ +
( b−
a2
b
) n+ −
( 1+
a
b
) βH.(9)
Solutions of Eqs. (9) are represented by the broken lines with the same slopes and jumps as in the case of
����
���
���
��� ��� ��� ��� ���
�
�
��� ��� ��� ��� ���
�
�
β+ ��� J ���� I ���
�D� �
�
ξ ���
�S� S�S )
�
�
�E�
��
��Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 15: The same as in Fig. 3 for f = 0.3, g = −0.1, βH = 0.1.
zero magnetic field but shifted by the value (1 + a/b)βH along the ξ+ and ξ− axes. For g < f , two cases can be distinguished. In the case of βH < g/2 (βH is in units of ε0F ), drawn in Fig. 12, two segments, vertical and horizontal, lying on the axes, intersect each other and the situation with the three intersection points, dis cussed above, holds. In the case βH > g/2, the single intersection point of the two lines remains at any ξ0 (see Fig. 13). The sp spectra and occupation numbers are shown in Fig. 14 for f = 0.3, g = −0.1, βH = 0.04 and in Fig. 15 for the same values of f and g but at βH = 0.1. In both the cases, the magnetic field promotes splitting the sp spectra ξ+(p) and ξ−(p) but does not influence their flat parts. In conclusion, in the Noziereslike model of a Fermi sys
tem with scalar and spindependent longrange forces, it is shown that the fermion condensation, occurring in the vicinity of a phase transition, results in the emergence of the weak magnetization which precedes the ferromag netic transition.
Acknowledgments
We thank A. Lichtenstein, G. Kotliar and V. Yakovenko for many valuable discussions. This research was supported in part by the National
Science Foundation under Grants No. PHY9907949 and No. PHY0140316, by the McDonnell Center for the Space Sciences at Washington University, and by Grant No. NS1885.2003.2 from the Russian Ministry of Indus try and Science (VAK and MVZ). One of the authors (VAK) thanks the University of California (Santa Bar
7 bara) for the kind hospitality.
[1] S. Doniach and S. Engelsberg, Phys. Rev. Lett. 17, 750, 1966.
[2] V. A. Khodel and V. R. Shaginyan, JETP Lett. 51, 553, 1990.
[3] L. D. Landau, JETP 30, 1058, 1956; 35, 97, 1958. [4] P. Nozieres, J. Phys. France I 2, 443, 1992. [5] G. E. Volovik, JETP Lett. 53, 222, 1991.
[6] V. A. Khodel, J. W. Clark and V. R. Shaginyan, Solid St. Comm. 96, 353, 1995.
[7] J. Dukelsky, V. A. Khodel, P. Schuck and V. R. Shaginyan, Z. Phys. B 102, 245, 1997.
[8] M. V. Zverev and M. Baldo, JETP 87, 1129, 1998; J. Phys. Cond. Mat. 11, 2059, 1999.
iv :c
on d
m at
/0 30
83 09
v1 [
co nd
m at.
str el
] 15
A ug
20 03
Spin degrees of freedom and flattening of the spectra of singleparticle
excitations in strongly correlated Fermi systems
V. A. Khodel1, P. Schuck2, and M. V. Zverev1 1 Kurchatov Institute, Russian Research Center,
Moscow, 123182 Russia 2 Groupe de Physique The´orique,
Institut de Physique Nucle´aire,
F01406 Orsay Cedex, France
(Dated: August 13, 2013)
The impact of longrange spinspin correlations on the structure of a flat portion in singleparticle spectra ξ(p), which emerges beyond the point, where the Landau state loses its stability, is studied. We supplement the wellknown Nozieres model of a Fermi system with limited scalar longrange forces by a similar longrange spindependent term and calculate the spectra versus its strength g. It is found that Nozieres results hold as long as g > 0. However, with g changing its sign, the spontaneous magnetization is shown to arise at any nonzero g. The increase of the strength g is demonstrated to result in shrinkage of the domain in momentum space, occupied by the flat portion of ξ(p), and, eventually, in its vanishing.
The investigation of flattening of singleparticle (sp) spectra of Fermi liquids is dated back to Ref. [1], where longrange correlations, enhanced in the vicinity of an impending ferromagnetic phase transition, were shown to result in the divergence of the effective mass M∗ at the transition point. Later in [2] an idea of the socalled fermion condensation, i.e. a rearrangement of the Lan dau state, occurring beyond a critical point in strongly correlated Fermi systems with longrange effective inter actions, was suggested. A striking feature of this rear rangement is ”swelling” of the Fermi surface (FS), i.e. the occurrence of a completely flat portion ξ(p) = 0, called the fermion condensate (FC), in the spectrum ξ(p), mea sured from the FS. To gain insight into the problem of fermion condensa
tion, let us turn to the Dyson equation, rewriting it in the form
ξ(p) = ξ0 p +Σ(p, ξ(p)) , (1)
appropriate for finding the FC solutions ξ(p) = 0. The sp mass operator Σ is usually determined by the formula
Σ(p, ε) =
∫ Wp, ε,p1, ε1)G(p1, ε1)
d4p1 (2pi)4i
, (2)
whereW is an effective interaction between particles, and G is the sp Green function. It is worth noting that the imaginary part of Σ(p, ε) vanishes at ε = 0. Therefore, in searching for the FC solutions, only the real part of Σ is relevant. For a long time, it was reckoned that in homogeneous
Fermi systems, there exists a onetoone correspondence between the momentum p and the sp energy ξ, at least, close to the FS, i.e. the derivative (dξ/dp)F is always positive, a postulate, virtually being a cornerstone of the Landau theory of Fermi liquid [3]. However, in sys tems with longrange forces it might be incorrect. Com pelling evidence for that is provided by a phenomenologi cal model, suggested by Nozieres [4], whereW is taken as
0(ξξp)/f
0
1
n(ξ)
ξ
FIG. 1: Graphical illustration of the solution of Eq. (3).
a constant in the coordinate space. Thus, in this model Eq. (1) takes the form
ξ − ξ0 p = fn(ξ) , (3)
where f is an effective coupling constant, the value of which is taken to be positive, while n(ξ) is the Landau quasiparticle momentum distribution, being 0 at positive ξ, and 1, at negative ones. This equation is easily solved, but we concentrate on its graphical solution (see Fig. 1). Let us draw both the r.h.s. and the l.h.s. of Eq. (3) as functions of ξ, taking p as an input parameter. The l.h.s. of (3), depicted at different p, provides a set of par allel straight lines, while the r.h.s. forms a kink, the ver tical segment of which is located at ξ = 0, no matter what is the input. Crossing points yield the sp spectrum ξ(p). If they lie on any of the two horizontal pieces of the kink, ξ(p) does coincide with the sp spectrum of ideal Fermi gas. An unconventional situation occurs when straight lines cross the steep section of the kink. In this case, the intersection point remains 0 in a finite momentum inter val [pi, pf ]. This plateau ξ(p) = 0, lying exactly at the FS, does form the fermion condensate. Since ξ(p) = 0 in the finite volume of momentum space, the density of states ρ(ε) acquires a singular term ρ(ε) ∼ δ(ε). Another salient feature of the fermion condensation phenomenon
2 is that in the FC domain, the quasiparticle momentum distribution n0(p) differs from the Landau one. Indeed, upon setting ξ = 0 in Eq.(3) one finds
n0(p) = − ξ0 p
f , pi < p < pf , (4)
in contrast to the conventional step function nF (p) = θ(pF − p). Thus, the FS does swell, whereas the basic assertion of Landau theory fails. In dealing with the FC problem, attention is usually
paid to the spinindependent part of the effective inter action W (see e.g. [2, 4, 5, 6, 7, 8]). Here we investigate effects, associated with longrange spinspin components of W , which involve the spinup quasiparticle distribu tion n+(p) and spindown one n−(p). It is instructive to start such an analysis with a generalization of the Nozieres model [4], supplementing it by longrange spin spin terms. As a result, one obtains two equations
ξ+(p) = ξ 0 p +
1
2 f (n+(p) + n−(p))
+ 1
2 g (n+(p)− n−(p)) ,
ξ−(p) = ξ 0 p +
1
2 f (n+(p) + n−(p))
− 1
2 g (n+(p)− n−(p)) , (5)
with a new constant g, specifying the longrange spin spin component of the model effective interaction W . It is worth noting that in the case of spin fluctuations with the nonzero critical momentum qc ≪ pF , presumably relevant to twodimensional liquid He3, the constants in Eq. (5) are related with each other: g = −f/3. In conventional Fermi liquids, the spontaneous spin
S arises only if strength of the spinspin interaction obeys the Pomeranchuk condition Gρ(0) > 1, where ρ(0) = pFM
∗/pi2 is the density of states at the FS in the Landau theory. In the case at issue, the density of states is infinite, so that one can expect the emergence of the spontaneous magnetization at any g < 0. To facilitate the solution of the problem, we recast the
system (5) to the form
ξ+ = ( 1−
a
b
) ξ0 p + a
b ξ− +
( b−
a2
b
) n− ,
ξ− = ( 1−
a
b
) ξ0 p + a
b ξ+ +
( b−
a2
b
) n+ , (6)
where a = (f + g)/2 and b = (f − g)/2. Let us now draw the plot ξ+(ξ−), proceeding from the first of Eqs.(6) and treating p as an input parameter. At ξ− > 0 the function ξ+(ξ−) is given by the straight line ξ+ = (1 − a/b)ξ
0 p +
aξ−/b, while at ξ− < 0 this straight line is slightly shifted: ξ+ = (1−a/b)ξ
0 p +(b2−a2)/b+aξ−/b. The drawing of its
counterpart ξ− from the second of Eqs.(6) yields another couple of straight lines: ξ− = (1 − a/b)ξ
0 p + aξ+/b at
ξ+ > 0, and ξ− = (1 − a/b)ξ 0 p − (b2 − a2)/b + aξ−/b at
ξ 
ξ 
ξ +
ξ +
ξ +
ξ 
(a)
A
ξ 
0
0
ξ +
(b)
0
B0
B1
'B1
ξ 
0
ξ +
0
C1
'C1 (c)
C0
0
0
D1
(d)
'D1
D0
0
E 0
(e) 0
FIG. 2: Graphical illustration of the dynamics of solutions of the set (6) with increasing ξ0p in the case of g < f .
ξ+ < 0. As we shall see, both these curves cross each other either at one or at three points, providing several solutions for the spectra ξ±(p).
Some elucidation of the procedure is ensured by Fig. 2, where the evolution of solutions of the set (6) versus the input parameter ξ0
p is shown. The solid line corresponds
to the dependence ξ+(ξ−), given by the first equation of the set (6), while the dashed one shows the same depen dence resulting from the second equation of the set. The intersection points are indicated by letters. Five panels of Fig. 2 show five different cases referring to different ξ0 p . In the first case (panel (a)), two zigzag lines have
one intersection point A, corresponding to the single so lution where both ξ+(p) and ξ−(p) are negative. With increasing ξ0
p , the intersection point moves to the origin.
At a certain value of this variable, depending on the pa rameter f only, bifurcation emerges, and there appear three intersection points (B0, B1 and B
′ 1 on the panel
(b)). One of them, B0, lies at the origin, while the other two, B1 and B
′ 1, on the ξ+ and ξ− axes, respectively. A
similar situation occurs in the panels (c) and (d). The intersection at the origin persists, while the other two crossing points lie inside the quadrants on the panel (c) (C1 and C
′ 1) and on the ξ− and ξ+ axes on the panel (d)
(D1 and D ′ 1). The last (e) panel shows the fifth case with
one intersection point at positive values of ξ+ and ξ−.
The chain A−B0−C0−D0−E is associated with the solution (denoted as Φ0), for which the sp spectra coin cide: ξ+(p) = ξ−(p) = 0 within the interval [pi, pf ], the length of which is determined by the parameter f only, as if there was no spinspin interaction at all. The analysis
3 0,1
0,0
0,1
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.1 f=0.3
(a) 
+
ξ +,
(p) p/pF
(b)
+n+,(p)
(c)
p/pF
s(p)
FIG. 3: Singleparticle spectra ξ±(p) in units of ε0F = p 2 F/2M
(panel (a)), occupation numbers n±(p) (panel (b)) and spon taneous spin s(p) = n−(p)− n+(p) (panel (c)), calculated for f = 0.3, g = −0.1 in units of ε0F .
0,1
0,0
0,1
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.01 f=0.3
(a)

+
ξ +,
(p) p/pF
(b)
 +
n +,
(p)
(c)
p/pF
s(p)
FIG. 4: The same as in Fig. 3 for f = 0.3, g = −0.01.
shows that, as long as the constant g remains positive, it is the solution Φ0 that has the lowest energy compared with the others. As for solutions, related to the chains
A−B1−C1−D1−E and A−B ′ 1−C
′ 1−D
′ 1−E (Φ1 and
Φ′1, respectively), the situation is different. The solution Φ1, corresponding to the first chain, is shown in Fig. 3. We see that the FC region, the boundaries of which were insensitive to g at any g > 0, is destroyed: the spectra
0,1
0,0
0,1
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.033 f=0.1
(a)

+
ξ +,
(p) p/pF
(b)
+ n
+, (p)
(c)
p/pF
s(p)
FIG. 5: The same as in Fig. 3 for f = 0.1, g = −0.033.
����
���
���
��� ��� ���
�
�
��� ��� ���
�
�
J ���� I ���
�D� �
�
ξ ���
�S�
S�S )
�
�
�E�
��
�� Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 6: The same as in Fig. 3 for f = 0.3, g = −0.3.
ξ+(p) and ξ−(p) repell each other. As the momentum p moves from the lower point, where these spectra simultaneously attain the FS, only the FC plateau at the spectrum ξ+(p) survives, with the energy splitting between the spectra ξ+(p) and ξ−(p) growing linearly with the p increase. Attaining the maximum at the point, where the plateau at ξ+(p) vanishes, this splitting ceases to increase and remains constant until the point, where a new FC plateau emerges at the spectrum ξ−(p), and, finally, both the spectra once again merge at the point, where the FC disappears forever.
4 ����
���
���
��� ��� ���
�
�
��� ��� ���
�
�
J ���� I ���
�D� �
�
ξ ���
�S�
S�S )
�
�
�E�
��
�� Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 7: The same as in Fig. 3 for f = 0.1, g = −0.1.
ξ 
ξ 
ξ +ξ+
ξ +
ξ 
(a)
A
ξ 
0
0
ξ +
(b)
0
B0 B2
'B2
B1
'B1
ξ 
0
ξ +
0
C2
'C2 (c)
C0
0
0
D2 D1
(d)
'D2 'D1
D0
0
E
0
(e) 0
FIG. 8: The same as in Fig. 2 for the case of g > f .
The second chain refers to the solution Φ′1, symmetric to the previous one, with changing± by ∓ and, thus, cor responds to the opposite sign of the projection of sponta neous spin S on the fixed axis. Since both the directions are equivalent, two solutions with the nonzero S have equal energies. It can be verified that at any g < 0, these solutions provide the minimum of the ground state en ergy, while the solution Φ0 with S = 0, the maximum. Indeed, calculations yield the gain in energy of the so
0,3
0,0
0,3
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.15 f=0.1
(a)

+
ξ +,
(p) p/pF
(b)
+n+,(p)
(c)
p/pF
s(p)
FIG. 9: The same as in Fig. 3 for f = 0.1, g = −0.15.
0,3
0,0
0,3
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.45 f=0.3
(a) 
+
ξ +,
(p) p/pF
(b)
+n+,(p)
(c)
p/pF
s(p)
FIG. 10: The same as in Fig. 3 for f = 0.3, g = −0.45.
lutions Φ1 and Φ ′ 1, as compared with the energy of the
Landau state E1/EL−1 ≃ −0.027, while for the solution Φ0, one obtains E0/EL − 1 ≃ −0.019. Thus, the reason for these alterations of the FC is the
occurrence of the spontaneous spin
S = ∑
p
s(p) = ∑
p
( n−(p)− n+(p)
) . (7)
arising at any g < 0. For illustration, the spin density s(p) is drawn in panel (c). We infer that the spin density s(p) differs from 0 only in the region [pi, pf ]. Therefore
5 0,3
0,0
0,3
0,9 1,0 1,1
0
1
0,9 1,0 1,1 0
1
g=0.45 f=0.1
(a) 
+
ξ +,
(p) p/pF
(b)
+ n
+, (p)
(c)
p/pF
s(p)
FIG. 11: The same as in Fig. 3 for f = 0.1, g = −0.45.
the total spin S depends primarily on the value of pa rameter f . For the parameters f and g, given in Fig. 3, one obtains S/ρ ≃ 0.299.
For comparison, results for the other two sets of the parameters: f = 0.3, g = −0.01 and f = 0.1, g = −0.033 are drawn in Figs. 4 and 5, respectively. The striking feature of these results is that for the given value f = 0.3 the magnitude of the spontaneous spin S for g = −0.01, S/ρ ≃ 0.232, is of the same order as in the case of g = −0.1. The energy gain of the solution with the spontaneous spin, corresponding to f and g for Fig. 4, E1/EL − 1 ≃ −0.019. Fig. 5 shows that multiplication of both the parameters f and g by the factor 1/3 results in the scaling of the spectra, occupation numbers, and spontaneous spin by the same factor.
The situation drastically changes if the magnitude of the spinspin force becomes comparable to that of the scalar one. This is seen in Figs. 6 and 7 show ing results obtained for the sets f = 0.3, g = −0.3 and f = 0.1, g = −0.1. We see that the FC prac tically disappears, while the Landau state reappears. Now the occupation numbers n+(p) and n−(p) have the Migdal jump from 0 to 1 at the points pi and pf , re spectively, and the spontaneous spin attains the maxi mum value, which is equal to the phase volume of the spherical layer between pi and pf (S/ρ ≃ 0.448 and S/ρ ≃ 0.149, for the parameters of Figs. 6 and 7, re spectively). Fig. 8 illustrates the evolution of solutions for the set (6) versus ξ0
p in the case when the abso
lute value of the spinspin constant g is larger than the scalar one f . In this case, the set (6) has five solutions: Φ0, corresponding to the chain A−B0−C0−D0−E; Φ1 and Φ′1, corresponding to the chains A−B1−C0−D1−E and A−B′1−C0−D
′ 1−E, respectively; and a new couple
A
0 ξ 
(a)
0
ξ +
(b)
B
0
0
ξ 
ξ +
(c)
C
0
0
ξ 
ξ +
'D1 D1
D0
(d)
0
0 ξ 
ξ +
'E1
E1
(e)
0
0
ξ 
ξ +
(f) 'F1
F0 F1
0
0
ξ 
ξ +
(g)
G
0
0
ξ 
ξ +
(h)
H
0
0
ξ 
ξ +
(i)
I
0
0
ξ 
ξ +
FIG. 12: Graphical illustration of the dynamics of solutions of the set (9) with increasing ξ0p in the case of g < f , βH < g/2.
of solutions, Φ2 and Φ ′ 2, corresponding to the chains
A−B2−C2−D2−E and A−B ′ 2−C
′ 2−D
′ 2−E. The calcu
lation of the energies yields the solutions Φ2 and Φ ′ 2,
equal in energy, beating Φ0, Φ1 and Φ ′ 1. The solution
Φ2 is drawn in Figs. 9, 10, 11 for three sets of parame ters: f = 0.1, g = −0.15; f = 0.3, g = −0.45; f = 0.1, g = −0.45. These figures show that in this case the spon taneous spin is defined only by the parameter g, and the parameter f defines the magnitudes of the jumps in the spectra ξ+(p) and ξ−(p). It is worth noting that, while the flat portions of the spectra, corresponding to FC, disappear and the Landau state is recovered, the phan tom of the FC manifests itself in the emergence of the spontaneous spin.
The above analysis can be generalized to the case of an external magnetic field. In this case, the equations for the sp spectra read
ξ+(p) = ξ 0 p +
1
2 f(n+(p) + n−(p))
+ 1
2 g(n+(p)− n−(p)) + βH ,
ξ−(p) = ξ 0 p +
1
2 f(n+(p) + n−(p))
− 1
2 g(n+(p)− n−(p))− βH , (8)
where H is the effective magnetic field acting on the fermion spin and β is the magnetic moment of the fermion. Upon rewriting Eqs. (8) to the form convenient
6 ξ 
ξ 
ξ +
ξ +
ξ +
ξ 
(a)
A
ξ 
0
0
ξ +
(b)
0
B
ξ 
0
ξ +
0
C
(c)
0
0
D
(d) 0
E 0
(e) 0
FIG. 13: The same as in Fig. 12 for the case of βH > g/2.
����
���
���
��� ��� ���
�
�
��� ��� ���
�
�
β+ ���� J ���� I ���
�D� �
�ξ ���
�S�
S�S )
�
�
�E�
��
��Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 14: The same as in Fig. 3 for f = 0.3, g = −0.1, βH = 0.04.
for the graphical analysis, we obtain
ξ+ = ( 1−
a
b
) ξ0 p + a
b ξ− +
( b−
a2
b
) n− +
( 1+
a
b
) βH,
ξ− = ( 1−
a
b
) ξ0 p + a
b ξ+ +
( b−
a2
b
) n+ −
( 1+
a
b
) βH.(9)
Solutions of Eqs. (9) are represented by the broken lines with the same slopes and jumps as in the case of
����
���
���
��� ��� ��� ��� ���
�
�
��� ��� ��� ��� ���
�
�
β+ ��� J ���� I ���
�D� �
�
ξ ���
�S� S�S )
�
�
�E�
��
��Q ���
�S�
�
�F�
S�S )
V�S�
��
�
FIG. 15: The same as in Fig. 3 for f = 0.3, g = −0.1, βH = 0.1.
zero magnetic field but shifted by the value (1 + a/b)βH along the ξ+ and ξ− axes. For g < f , two cases can be distinguished. In the case of βH < g/2 (βH is in units of ε0F ), drawn in Fig. 12, two segments, vertical and horizontal, lying on the axes, intersect each other and the situation with the three intersection points, dis cussed above, holds. In the case βH > g/2, the single intersection point of the two lines remains at any ξ0 (see Fig. 13). The sp spectra and occupation numbers are shown in Fig. 14 for f = 0.3, g = −0.1, βH = 0.04 and in Fig. 15 for the same values of f and g but at βH = 0.1. In both the cases, the magnetic field promotes splitting the sp spectra ξ+(p) and ξ−(p) but does not influence their flat parts. In conclusion, in the Noziereslike model of a Fermi sys
tem with scalar and spindependent longrange forces, it is shown that the fermion condensation, occurring in the vicinity of a phase transition, results in the emergence of the weak magnetization which precedes the ferromag netic transition.
Acknowledgments
We thank A. Lichtenstein, G. Kotliar and V. Yakovenko for many valuable discussions. This research was supported in part by the National
Science Foundation under Grants No. PHY9907949 and No. PHY0140316, by the McDonnell Center for the Space Sciences at Washington University, and by Grant No. NS1885.2003.2 from the Russian Ministry of Indus try and Science (VAK and MVZ). One of the authors (VAK) thanks the University of California (Santa Bar
7 bara) for the kind hospitality.
[1] S. Doniach and S. Engelsberg, Phys. Rev. Lett. 17, 750, 1966.
[2] V. A. Khodel and V. R. Shaginyan, JETP Lett. 51, 553, 1990.
[3] L. D. Landau, JETP 30, 1058, 1956; 35, 97, 1958. [4] P. Nozieres, J. Phys. France I 2, 443, 1992. [5] G. E. Volovik, JETP Lett. 53, 222, 1991.
[6] V. A. Khodel, J. W. Clark and V. R. Shaginyan, Solid St. Comm. 96, 353, 1995.
[7] J. Dukelsky, V. A. Khodel, P. Schuck and V. R. Shaginyan, Z. Phys. B 102, 245, 1997.
[8] M. V. Zverev and M. Baldo, JETP 87, 1129, 1998; J. Phys. Cond. Mat. 11, 2059, 1999.
Comments