# Quasi-Particles, Thermodynamic Consistency, and the Gap Equation.

1. IntroductionThe Bogoliubov transformation is one of the main discoveries in mathematical physics. In this article it is shown how to compute the thermodynamic potentials--functions of the (T, V, N) variables--in superconductivity by means of the Bogoliubov transformation. It is demonstrated that the Bogoliubov transformation avoids the problem of the thermodynamic consistency of the quasi-particle approach. The quasi-particles vacuum energy and the dependence of the chemical potential on the electron density are derived.

In the low temperature limit, the thermodynamic behaviour of any macroscopic physical system is determined by the structure of the energy levels of the system which belong to a neighbourhood of the ground state energy. The experimental data show that, usually, these energy levels are organized in such a way that one can give quite a good description of their thermodynamic properties by means of a gas of quasi-particles [1,2]. The energy spectrum of a single quasi-particle is specified by its dispersion relation, which is determined by the nature of the interactions among the "atoms" or the "elementary constituents" of the system.

Let us consider the framework in which the gas of quasi-particles describes the thermal fluctuations of the system above its ground state. In this formalism, the number of quasi-particles is not conserved--in fact, it vanishes at zero temperature--and therefore the quasi-particles chemical potential is zero. Even if the system is composed of strongly interacting atoms--like in a quantum liquid--the quasi-particles are weakly interacting, and then one can determine, with a good approximation, the thermodynamic potentials of the system. Indeed, in the low temperature limit the density of quasi-particles becomes very low, so that their mutual interactions produce minor effects. Thus a crucial part of the information on the thermal properties of the whole system is mainly contained in the form of the dispersion relation of a single quasi-particle. Let e(p) be the energy of a single quasi-particle in the state labelled by the value p of the momentum. The function [epsilon](p) can be a rather complicated function of the momentum, which can be determined by the experimental data, or it can be deduced--in some approximated form-- from the Hamiltonian of the system. In superconductors, for instance, the dispersion relation for a single quasi-particle [3, 4] assumes the form

[mathematical expression not reproducible] (1)

where [mu] denotes the chemical potential for the electrons. The gap function [[DELTA].sub.p] can be different from zero in a neighbourhood of the surface of the Fermi sphere and satisfies

[mathematical expression not reproducible] (2)

in which [w.sub.D] indicates the Debye frequency of the ions lattice. The gap [DELTA](T) is not vanishing when the temperature T takes values below the critical value [T.sub.c], and [DELTA](T) = 0 for T [greater than or equal to] [T.sub.c]. As a result, the energy [epsilon] of a single quasi-particle nontrivially depends on the temperature, [epsilon] = [epsilon] (p, T). But when the dispersion relation [epsilon] (p, T) nontrivially depends on the temperature, the standard expressions of the thermodynamic potentials--derived in statistical mechanics--for a gas of weakly interacting constituents do not satisfy the thermodynamic relations (or Maxwell relations).

More precisely, if A represents a standard thermodynamic potential which is computed in statistical mechanics when the atomic energy [[epsilon].sub.0] (p) does not depend on the temperature, let us denote by [A.sub.q] the corresponding potential for the gas of quasi-particles which is obtained from A by means of the substitution [[epsilon].sub.0] (p) [right arrow] [epsilon](p, T). For instance, the free energy [F.sub.q] and the internal energy [U.sub.q] of a gas of noninteracting quasi particles with energy [epsilon](p, T), satisfying the Fermi-Dirac statistics, are given by

[mathematical expression not reproducible]. (3)

where g denotes the spin degeneration factor. When the thermodynamic relations are satisfied, the entropy S is given by S = -[partial derivative]F(T,V)/[partial derivative]T; and since F = U - TS, one has U = F - T([partial derivative]F/[partial derivative]T). However, if [[partial derivative][epsilon](p, T)/[partial derivative]T] [not equal to] 0, one finds

[U.sub.q] [not equal to] [F.sub.q] -T [([partial derivative][F.sub.q]/[partial derivative]T).sub.v] (4)

and therefore expressions (3) do not satisfy the thermodynamic relations.

In addition to their agreement with experiments, the thermodynamic relations codify to some extent the laws of thermodynamics and are necessary for the sake of logic. The requirement of validity of the thermodynamic relations will be called the thermodynamic consistency. The thermodynamic consistency implies that, in general, the thermodynamic potentials of the macroscopic systems cannot be exactly equal to the potentials [A.sub.q] of the quasi-particles. The presence of adjusting terms [delta][A.sub.q] is needed in order to rectify the [A.sub.q] expressions and make the final combinations [A.sub.q] + [delta][A.sub.q] correct. As far as the computation of SA is concerned, the knowledge of the quasi-particles distribution is of no help, because [delta][A.sub.q] is essentially determined by the quasi-particles vacuum energy.

One of the purposes of the present article is to produce the correction terms to the quasi-particles potentials in the case of superconductivity by means of the Bogoliubov-Valatin formalism. It is shown that the accurate determination of the quasi-particles vacuum energy completes the construction of the thermodynamic potentials satisfying the Maxwell relations. It is verified that, in agreement with the Landau principle, the value of the gap corresponds to a stationary point of the grand potential [OMEGA](T, V, [mu]). Then the expressions of the thermodynamically consistent potentials as functions of the (T, V' N) variables are obtained, and how the chemical potential is related to the electron density is determined. As a final check, a derivation of the superconducting condensation energy from [OMEGA](T, V, [mu]) is presented.

The thermodynamic consistency in the presence of nontrivial medium-dependent dispersion relation has been considered, for instance, by Shanenko et al. [5, 6] in the case of clustering matter and by Gorenstein and Yang [7] in the case of a gluon plasma. This issue has been elaborated also in [8-12]. The Landau principle, which can be interpreted as an equilibrium condition for ordinary thermal states, implies the thermodynamic consistency. It has been argued [12] that, conversely, the validity of the thermodynamic consistency does not necessarily imply the Landau principle. The results for the superconducting electrons are in agreement with the analysis of [7].

In superconductivity, the dependence of the thermodynamic potentials on the statistical mean value N of the number of electrons is a rather nontrivial issue because when T < [T.sub.c], the U(1) symmetry which is related to the number of electrons in superconductors is broken [3, 4]. The solution to this problem which is presented in the following chapters and the outcomes of the computed thermodynamic potential [OMEGA](T, V [mu]) are in agreement with the known results [3, 4, 13, 14] on superconductivity. In the present work it is shown that the identification of the variables which are associated with the quasi-particles leads to a precise determination of the quasi-particles vacuum energy. The method, which is presented here to solve this question, can find applications also in the study of topological superconductivity and in the new developments on the topological states of matter [15-24].

The basic concepts which are connected with the use of the quasi-particles distribution are briefly recalled in Section 2. Starting from the BCS Hamiltonian, in Section 3 the Bogoliubov-Valatin formalism [25, 26] is used to determine the grand potential Q, according to the procedure envisaged, for instance, by Rickayzen [27]. By means of a Bogoliubov transformation acting on the creation and annihilation operators of the electron fluctuations around the Fermi sphere, the full expression of [OMEGA](T, V, [mu]) satisfying the Maxwell relations is obtained. The thermodynamic consistency of the result is discussed in Section 4, where a new derivation of the gap equation is described. In Section 5 the relation connecting the chemical potential with the electron density is derived, and the consistent expressions of various thermodynamic potentials--functions of the (T, V, N) variables--of the superconducting electrons are produced in the low temperature limit.

2. Quasi-Particles Distribution

Let us briefly recall the basic notions which are related to the use of the quasi-particles formalism. The description of the low temperature physics by means of quasi-particles represents a phenomenological approach in which the dispersion relation of a single quasi-particle nontrivially depends on the macroscopic variables of the material [epsilon](p) = [epsilon](p , T, [mu], ...). Thus, in addition to the temperature, the failure of the thermodynamic consistency for the quasi-particle potentials [A.sub.q] really concerns several state variables. In certain circumstances the dependence of [epsilon] on these variables is weak enough so that, in a limited range of variability, one can assume [28] that [epsilon] only depends on the momentum. But in general--and in particular when a phase transition occurs--one cannot neglect the dependence of [epsilon] on the macroscopic variables; this dependence will be denoted by [epsilon] = [epsilon](p , T).

As it has been suggested by Landau [2], a possible way to connect the thermodynamic potentials of the system with the quasi-particles functions makes use of the quasi-particles distribution. The only thermodynamic potential whose absolute value is fixed--and cannot be modified by any additive constant--is the entropy S:

S = k

* ln {number of microstates with fixed macro conditions}. (5)

Since the quantum states of the system are described (in the low temperature limit) precisely by the gas of quasi-particles, the entropy of the system is equal to the entropy of the quasi-particles gas. For noninteracting quasi-particles, the number of microstates in a given macrostate can be determined by means of the distribution [n.sub.p] of the quasi-particles. Indeed, in the quantum case of Fermi statistics, the analogue of the Boltzmann H-functional is given [2] by

[mathematical expression not reproducible]. (6)

Let us recall that, in the case of particles with vanishing chemical potential, the entropy represents the thermodynamic potential for the variables (U, V). This means that, with fixed (U,V), the stable thermal state of the system corresponds to a maximum of S. Therefore the distribution [n.sub.p] can be determined by the requirement that, for fixed (U, V), the variation of S with respect to a generic fluctuation [delta][n.sub.p] must vanish. The associated (conditioned) variational principle takes the form

[partial derivative]S/[partial derivative][n.sub.p] - [lambda] [partial derivative]U/[partial derivative][n.sub.p] = 0, (7)

where the value of the Lagrange multiplier [lambda] is given by [lambda] = 1/T as a consequence of the Maxwell relation TdS - dU PdV = 0. Because of the explicit dependence (6) of S on V, there is no need to introduce a Lagrange multiplier for the volume. Since a modification of [n.sub.p] causes the following change in the energy

[delta]U = Vq/[h.sup.3] [integral] [d.sup.3] p[epsilon] (p' T) [delta][n.sub.p]' (8)

the solution of (7) is given precisely by the Fermi-Dirac distribution

[mathematical expression not reproducible] (9)

By using (9) several variables can be computed. Unfortunately, the knowledge of the distribution (9)--or of S(T, V)--alone is not enough to determine all the potentials; for instance, the complete expressions of the free energy F and of the internal energy U cannot be obtained from (9).

3. Computation of the Grand Potential

The low energy BCS Hamiltonian [3] for conducting electrons in superconductor metals can be written as

[mathematical expression not reproducible] (10)

where [b.sup.[??].sub.ps] and [b.sub.ps] denote the creation and annihilation operators for one electron in the state |p , s), where s = [+ or -] refers to the value of one component of the spin:

[mathematical expression not reproducible] (11)

The interaction kernel [U.sub.p,q] is related to the amplitude of the electron-electron scattering. [U.sub.p,q] does not depend on the values of the electron thermodynamic variables, as the temperature T, the volume V, and the chemical potential u. The main effects of the interactions between electrons in the superconducting materials--which are relevant for the superconducting phase transition--are found when [mathematical expression not reproducible], and in this region [U.sub.p,q] can be approximated by -[U.sub.0]/V, where [U.sub.0] is a positive constant. The grand partition function [??], which is defined by

[mathematical expression not reproducible], (12)

can be interpreted as the partition function of a system with total Hamiltonian

[mathematical expression not reproducible] (13)

in which the "effective" kinetic energy [eta](p) of one electron (where p = [absolute value of p]) is given by

[eta](p) = [p.sup.2]/2m - [mu] (14)

Since trace (12) cannot be evaluated exactly, one can introduce a self-consistent approximation (some kind of a mean field approximation) which permits proceeding with the computation. In the grand canonical ensemble, let us introduce the statistical mean values

[mathematical expression not reproducible] (15)

In order to simplify the exposition, let us assume that [X.sup.*.sub.p] = [X.sub.p] (in fact, it turns out [3, 4,13] that one can always choose the phases of the creation and annihilation operators in such a way that the mean values (15) are real). The couple of operators [mathematical expression not reproducible] can be written as the sum of their means values plus a fluctuation term:

[mathematical expression not reproducible]. (16)

Now one can assume [13, 27] that the effects of the fluctuations are sufficiently small so that, inserting identities (16) into expression (13), the resulting term which is quadratic in the fluctuations can be neglected. One then finds

[mathematical expression not reproducible] (17)

where

[[DELTA].sub.p] = [[summation].sub.q] [U.sub.p,q] [X.sub.q]. (18)

In order to recover the quasi-particles description of the low temperature behaviour of the system, it is convenient to introduce new creation and annihilation operators. The corresponding procedure is composed of two steps:

(1) Introduction of the operators [d.sup.[??].sub.ps] and [d.sub.ps] which create and annihilate quasi-particles in the case of free electrons.

(2) Introduction of the operators [a.sup.[??].sub.ps] and [a.sub.ps]--which are obtained from [d.sup.[??].sub.ps] and [d.sub.ps] by means of a Bogoliubov transformation--that diagonalize Hamiltonian (17).

Step 1. The operators [d.sup.[??].sub.ps] and [d.sub.ps] should describe the one-particle fluctuations with respect to the configuration of the Fermi sphere. Then [d.sup.[??].sub.ps] and [d.sub.ps] are chosen in such a way that their vacuum state coincides with the free electrons ground state corresponding to the Fermi sphere. In the T[right arrow]^0 limit, the electron states with p = [absolute value of p] < [p.sub.0] are occupied whereas the states with p = [absolute value of p] [p.sub.0] are empty, where

[p.sub.0] = [square root of 2mu] (19)

Therefore [d.sup.[??].sub.ps] and [d.sub.ps] are defined [25] by

[mathematical expression not reproducible] (20)

The creation operator [d.sup.[??].sub.ps] of a quasi-particle corresponds to the creation of a one-electron hole in the Fermi sphere if [absolute value of p] < [p.sub.0] and the creation of a one-electron occupied state if [absolute value of p] > [p.sub.0]. As a result of the creation of a hole in the Fermi sphere, which corresponds to an electron state with momentum -p and spin component -s, the total momentum of the system increases by p and the total spin of the system increases by s. For this reason, when p < [p.sub.0] the quantum numbers of the operators [d.sup.[??].sub.ps] and [d.sub.ps] are opposite with respect to the quantum numbers of [b.sup.[??].sub.ps] and [b.sub.ps]. The operators [d.sup.[??].sub.ps] and [d.sub.ps] satisfy the canonical anticommutation relations

[mathematical expression not reproducible] (21)

The operator number of quasi-particles [[??].sub.qp] is given by [mathematical expression not reproducible] counts the number of holes inside the Fermi sphere plus the number of occupied states outside the Fermi sphere. The number of quasi-particles is not conserved, and thus quasi-particles have vanishing chemical potential. In the T [right arrow] 0 limit, the number of quasi-particles vanishes. By using the relation [mathematical expression not reproducible] and definition (20), one finds

[mathematical expression not reproducible] (22)

in which

[mathematical expression not reproducible] (23)

The gas of quasi-particles, which are associated with operators (20), can also be used to compute the thermodynamic variables for a degenerate Fermi gas (as an alternative to the Sommerfeld expansion). Note that, differently from the effective kinetic energy [eta](p) of one electron, the kinetic energy w(p) of one quasi-particle satisfies w(p) [greater than or equal to] 0.

Step 2. By using the operators [d.sup.[??].sub.ps] and [d.sub.ps], Hamiltonian (17) reads

[mathematical expression not reproducible] (24)

Since [mathematical expression not reproducible] one can write

[mathematical expression not reproducible] (25)

Thus operator (24) can be written in the form

[mathematical expression not reproducible] (26)

where

[mathematical expression not reproducible] (27)

By means of the diagonalizing matrix [M.sub.p] [epsilon] SU(2) that satisfies the relation

[mathematical expression not reproducible] (28)

in which

[mathematical expression not reproducible] (29)

one can introduce new creation and annihilation operators [a.sup.[??].sub.ps] and [a.sub.ps] according to

[mathematical expression not reproducible] (30)

One then obtains

[mathematical expression not reproducible] (31)

By using the operators [a.sup.[??].sub.ps] and [a.sub.ps], the effective Hamiltonian (26) assumes the diagonal form

[mathematical expression not reproducible] (32)

where

[mathematical expression not reproducible] (33)

The operators [a.sup.[??].sub.ps] and [a.sub.ps] verify

[mathematical expression not reproducible] (34)

and they represent the creation and annihilation operators for quasi-particles (the so-called bogolons) with dispersion relation [epsilon](p). The grand partition function (12) is given by

[mathematical expression not reproducible] (35)

and then the grand potential [OMEGA] for the system of conducting electrons turns out to be

[mathematical expression not reproducible] (36)

The introduction of Step 1 clarifies the nature of the bogolons and the resulting formalism is consistent with the vanishing of the bogolon chemical potential. Note that if the Bogoliubov transformation is directly applied to the creation and annihilation operators [b.sup.[??].sub.ps] and [b.sub.ps] of the electrons, one finds sign ambiguities [27] in the energy spectrum. The method that has been presented here has the advantage of avoiding these ambiguities. By the way, the correct solution of these ambiguities is rather nontrivial, because if the electron effective energy [eta](p) takes the place of w(p), the thermodynamic consistency requires that the right sign of the dispersion relation should be +[epsilon](p) when [absolute value of p] > [p.sub.0] and -[epsilon](p) when [absolute value of p] < [p.sub.0]. The correct value (33) of the energy [E.sub.0] of the vacuum of the quasi-particles is a fundamental ingredient in the following discussions in Sections 4 and 5. Group properties of the quasi-particles vacuum and features of the solution of the gap equation have been discussed, for instance, in [29-32].

In order to complete the derivation of [OMEGA], one has to impose the statistical self-consistency of approximation (17), in which the mean values (15) have been used. Let us briefly recall the method which is usually presented in literature for this purpose. The mean values (15) can be related to the mean values of the corresponding operators defined in terms of [a.sup.[??].sub.ps] and [a.sub.ps]. Since the quasi-particles created and annihilated by [a.sup.[??].sub.ps] and [a.sub.ps] behave as free particles, the only nontrivial mean value is given by ([a.sup.[??].sub.ps] and [a.sub.ps] = np, where the quasi- particle distribution [n.sub.p] is shown in(9). In this way, by using definition (18), one can obtain the gap equation [3,4, 27, 28, 33]

1 = [U.sub.0]/V [summation over P)] 1/2[epsilon](p) Th ([epsilon](p)/2kT], (37)

where [[summation].sub.p] denotes the integral in the neighbourhood of the surface of the Fermi sphere which is specified by the condition [mathematical expression not reproducible]. In the next section, an alternative derivation of the gap equation will be produced.

4. Thermodynamic Consistency and the Gap Equation

From definition (18) and the experimental fact that the superconducting electrons belong to a neighbourhood of the Fermi surface, where [U.sub.p,q] [equivalent] [U.sub.0]/V, it follows that one can put

[mathematical expression not reproducible] (38)

and then

- [[summation].sub.P] [[DELTA].sub.p] [X.sub.p] = V/[U.sub.0] [[DELTA].sup.2] (39)

Let us now consider expression (36) and let us assume, for the moment, that the value [DELTA] of the gap is a free undeterminate which is not specified by the gap equation; one has

[mathematical expression not reproducible], (40)

where

[mathematical expression not reproducible]. (41)

Let us denote by B the neighbourhood of the Fermi surface which is defined by the condition [mathematical expression not reproducible]. Then

[mathematical expression not reproducible] (42)

Since [OMEGA] is the thermodynamic potential associated with the variables (T, V, u), for fixed (T,V,p) the system tends to minimize Q. So, in the equilibrium state, the value A of the gap can be obtained by imposing the variational condition

[mathematical expression not reproducible] (43)

Since

[mathematical expression not reproducible]; (44)

one finds

[mathematical expression not reproducible] (45)

Moreover

[mathematical expression not reproducible] (46)

Equations (45)-(46) imply

[mathematical expression not reproducible] (47)

By using the identity

[mathematical expression not reproducible] (48)

the minimum condition (43) turns out to be equivalent to (apart from the trivial solution)

[mathematical expression not reproducible] (49)

which coincides precisely with the gap equation (37). In what follows, the value of the gap satisfying (49) will be simply denoted by [DELTA] = [DELTA](T, u).

To sum up, the grand potential [OMEGA] for the electrons system takes the form

[OMEGA] (T, V, u) = [F.sub.q] (T, V; A (T, u)) + [[epsion].sub.0] (V, u [DELTA] (T, u)), (50)

where [F.sub.q] denotes the free energyfor the gas of quasi-particles, with dispersion relation e(p) shown in (42) and [E.sub.0] of (41) represents the vacuum energy of the quasi-particles. The value [DELTA](T,[mu]) of the gap corresponds to a minimum of [OMEGA] for fixed (T, V, [mu]).

Let us now consider the issue of the thermodynamic consistency of result (50). The starting expression (12) is thermodynamically consistent, because the dependence of the grand partition function on the macroscopic variables (T, V, u) is the standard dependence of statistical physics. However, when the dependence of the gap on the variables T and u has been introduced, [DELTA] = [DELTA](T, u), the standard rules of statistical physics have been modified. Now, the dependence of [OMEGA] on the variables (T, V, [mu]) consists of two parts:

(i) The (thermodynamically consistent) explicit dependence on (T, V, [mu]) displayed in expressions (40)-(42)

(ii) The indirect dependence on (T, [mu]) through [DELTA] = [DELTA](T, [mu]).

Let us distinguish these two possibilities by means of the notation

[OMEGA] = [OMEGA]T, V, [mu]; [DELTA] (T, [mu])). (51)

When computing the derivatives of [OMEGA] one finds, for instance,

[mathematical expression not reproducible](52)

But since the value of the gap minimizes [OMEGA], the term [([partial derivative][OMEGA]/[partial derivative][DELTA]A).sub.T] vanishes, and then the derivative dA/dT does not contribute to [partial derivative][OMEGA]/[partial derivative]T. Therefore, only the thermodynamically consistent dependence of [OMEGA] on (T, V, u) really contributes to the derivates. This means that the Maxwell relations associated with [OMEGA] are satisfied.

Let us now consider a generic thermodynamic potential which is derived from [OMEGA] by means of appropriate Legendre transformations. As far as the variable [DELTA] is concerned, the so-called theorem of small increments [34] takes the form

[mathematical expression not reproducible] (53)

and it states that [DELTA] really corresponds to a stationary point of each potential (provided the corresponding set of thermodynamic variables is kept fixed). Consequently, when constructing the thermodynamic relations by means of the first derivatives of the potentials, [DELTA] effectively behaves like a constant term and does not alter the Maxwell relations. Thus, because of the validity of the gap equation, the quasi-particles description of the systems, which is given--in the low temperature limit--by the grand potential (50), is perfectly consistent with the validity of the standard thermodynamic relations.

This example of superconducting electrons is in agreement with the Landau principle, which states that the value of any order parameter entering a given thermodynamic potential actually corresponds to a stationary point of the same potential, with the appropriate thermodynamic variables held fixed.

5. Thermodynamic Potentials

Quite often, the complete expression (50) of the grand potential and, in particular, the complete expression (41) of the energy [E.sub.0] are not explicitly displayed in literature. In certain cases, the combination [[summation].sub.p][w(p) - [epsilon](p) - [[DELTA].sub.p] [X.sub.p]] is reported; however the term [E.sub.F](V, u) is missing. Consequently, in these cases the partial derivatives of [OMEGA] with respect to V and u are invalidated and thus the deduction of the thermodynamic potentials as functions of the variables (T, V, N) becomes problematic. Usually, the explicit derivation of the general dependence of the superconducting thermodynamic potentials on (T, V, N) represents a rather laborious task. So, let us now derive from (50) a few thermodynamic potentials as functions of (T, V, N) when T < [T.sub.c] in the low temperature limit. Let us compute, in particular, how the chemical potential gets modified by the superconducting phase transition and how [mu] depends on the electron density. The dependence of [mu], on (T, V, N) is indeed the fundamental ingredient for the determination of the consistent thermodynamic potentials.

In the T [right arrow] 0 limit, the quasi-particles distribution is mainly concentrated around the minimum of the energy spectrum. The energy (29) of a single quasi-particle can be written as [epsilon] = [[epsilon].sub.w] = [square root of [w.sup.2] + [[DELTA].sup.2] and, in a neighbourhood of the minimum w = 0, it can be approximated as

[[epsilon].sub.w] = [square root of [w.sup.2] + [DELTA].sup.2] [equivalent] [DELTA] + [w.sup.2]/2[DELTA]. (54)

Accordingly, in the computations of [F.sub.q] and of [E.sub.0] the integral in momentum space is effectively dominated by the integration in a neighbourhood of the surface of the Fermi sphere and it can be written as

[mathematical expression not reproducible](55)

The grand potential (50) reads

[mathematical expression not reproducible]. (56)

In the thermodynamic limit, the statistical mean value N of the number of electrons is specified by the relation N = -[partial derivative][OMEGA]/[partial derivative]u, which takes the form

[mathematical expression not reproducible] (57)

In the low temperature limit one has [DELTA]/[[epsilon].sub.F] [much less than] < 1 and (kT/[[epsilon].sub.F]) [much less than] 1, where [[epsilon].sub.F] denotes the Fermi energy

[mathematical expression not reproducible](58)

So, the chemical potential can be determined perturbatively in powers of (kT/[[epsilon].sub.F]) [much less than] 1 and [DELTA]/ [[epsilon].sub.F] [much less than] 1 from (57); to lowest nontrivial orders one finds

[mathematical expression not reproducible]. (59)

This relation gives the desired expression for the chemical potential. It should be noted that, in addition to the quasi-particles contributions to p, the last term on the right-hand-side of (59) describes precisely the effects of the quasi-particles vacuum rearrangement--which takes place at the critical point--on the chemical potential. The free energy F is given by F = [OMEGA] + N[mu]:

[mathematical expression not reproducible] (60)

The entropy S = -[partial derivative]F/[partial derivative]T turns out to be

[mathematical expression not reproducible] (61)

and it coincides with the entropy of the gas of quasi-particles, as it must be. Finally, the internal energy U = F + TS is given by

[mathematical expression not reproducible] (62)

The consequences of expressions (56)-(62) are in agreement with the known results [3,4,13,27,35,36] on superconductivity. As a check, let us compute the condensation energy. The energy of the quasi-particles gas can be obtained by means of the Gaussian approximation

[mathematical expression not reproducible] (63)

this energy contribution vanishes in the T [right arrow] 0 limit. Since [DELTA]/[??][w.sub.D] [much less than] 1, one finds

[mathematical expression not reproducible] (64)

In the T [right arrow] 0 limit, the value of [DELTA] is specified by the gap equation (49); in particular, when (3N[U.sub.0]/4V[[epsilon].sub.F]) [much less than] 1, one obtains [3, 4,13]

[mathematical expression not reproducible]. (65)

By means of relation (65), one gets

[mathematical expression not reproducible] (66)

Within the considered approximations, the value [E.sub.C] of the condensation energy is obtained by subtracting the energy of a free fermions gas at T = 0 from the system energy U(T = 0, V, N). Since the energy of a free electrons gas in the Fermi sphere is equal to (3/5)N[[epsilon].sub.F], the condensation energy [E.sub.C] turns out to be

[mathematical expression not reproducible]. (67)

which is in agreement with the results of the alternative derivations presented in [3, 4, 13].

6. Summary and Conclusions

In the Bogoliubov-Valatin approach, it has been shown that the relevant variables--which must be considered for the computation of the thermodynamic potentials--are associated with the energy fluctuations of the system around the Fermi sphere configuration. By means of a Bogoliubov transformation in these variables one gets the expression of the quasi-particles vacuum energy, which is the fundamental ingredient in the construction of the grand potential [OMEGA](T, V, u) for the superconducting electrons. The value of the gap corresponds to a minimum of [OMEGA] and the gap equation precisely represents the stationarity condition. When the gap equation is satisfied, [OMEGA] and the thermodynamic potentials that are derived from [OMEGA] satisfy the Maxwell relations. The expressions of the thermodynamic potentials as functions of the (T, V, N) variables have been derived. The dependence of the chemical potential on the electron density and the effects of the superconducting phase transition on the chemical potential have been computed in the low energy limit. As a final check, a derivation of the condensation energy from [OMEGA](T, V, [mu]) has been presented.

https://doi.org/10.1155/2017/6265427

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author wishes to thank E. Vicari for discussions.

References

[1] L. D. Landau, "On the theory of superfluidity of helium II," J. Phys. U.S.S.R, vol. 11, p. 91, 1947.

[2] L. D. Landau, "The theory of a Fermi liquid," Soviet Physics. JETP, vol. 3, pp. 920-925, 1957.

[3] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, "Theory of superconductivity," Physical Review Letters, vol. 108, pp. 1175-1204, 1957.

[4] L. D. Landau and E. M. Lifshitz, Course of Theoretical Physics, vol. 9 of Statistical Physics, Part 2, Theory of the Condensed State, Pergamon Press, 1981.

[5] A. A. Shanenko, E. P. Yukalova, and V. I. Yukalov, "Statistical model of quark-hadron matter," Il Nuovo Cimento A Series 11, vol. 106, no. 9, pp. 1269-1281, 1993.

[6] A. A. Shanenko, E. P. Yukalova, and V. I. Yukalov, "Statistical models of clustering matter," Physica A: Statistical Mechanics and its Applications, vol. 197, no. 4, pp. 629-666, 1993.

[7] M. I. Gorenstein and S. N. Yang, "Gluon plasma with a medium-dependent dispersion relation," Physical Review D, vol. 52, no. 9, pp. 5206-5212, 1995.

[8] T. S. Biro, A. A. Shanenko, and V. D. Toneev, "Toward thermodynamic consistency of quasiparticle picture," Physics of Atomic Nuclei, vol. 66, no. 5, pp. 982-996, 2003.

[9] S. Y. Yin and R. K. Su, "The intrinsic degree of freedom for quasiparticle in thermodynamics with medium effects," https:// arxiv.org/abs/0709.0179.

[10] V. M. Bannur, "Comments on quasiparticle models of quarkgluon plasma," Physics Letters B, vol. 647, p. 271, 2007.

[11] F. G. Gardim and F. M. Steffens, "Thermodynamics of quasiparticles," Nuclear Physics A, vol. 797, no. 1-2, pp. 50-66, 2007.

[12] F. G. Gardim and F. M. Steffens, "Thermodynamics of quasiparticles at finite chemical potential," Nuclear Physics A, vol. 825, no. 3-4, pp. 222-244, 2009.

[13] M. Tinkham, Introduction to Superconductivity, Dover Publications Inc., New York, NY, USA, 2nd edition, 2004.

[14] H. Rogalla and P. H. Kes, Eds., 100 Years of Superconductivity, CRC Press, Taylor & Francis Group, New York, NY, USA, 2012.

[15] A. Y. Kitaev, "Unpaired Majorana fermions in quantum wires," Usp. Fiz. Nauk (Suppl.), vol. 171, p. 131, 2001.

[16] A. P. Mackenzie and Y. Maeno, "The superconductivity of Sr2RuO4 and the physics of spin-triplet," Review of Modern Physics, vol. 75, p. 657, 2003.

[17] L. Fu and C. L. Kane, "Superconducting proximity effect and Majorana fermions at the surface of a topological insulator," Physical Review Letters, vol. 100, Article ID 96407, 2008.

[18] X.-L. Qi and S.-C. Zhang, "Topological insulators and superconductors," Review of Modern Physics, vol. 83, p. 1057, 2011.

[19] A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, and H. Shtrikman, "Zero-bias peaks and splitting in an Al-InAs nanowire topological superconductor as a signature of Majorana fermions," Nature Physics, vol. 8, no. 12, pp. 887-895, 2012.

[20] S. Nadjj-Perge, I. K. Drozdov, J. Li et al., "Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor," Science, vol. 346, no. 6209, pp. 602-607, 2014.

[21] Y. Qiu, K. N. Sanders, J. Dai et al., "Time reversal symmetry breaking superconductivity in topological materials," https:// arxiv.org/abs/1512.03519.

[22] C. Kallin and J. Berlinsky, "Chiral superconductors," Reports on Progress in Physics, vol. 79, Article ID 054502, 2016.

[23] C. Beenakker and L. Kouwenhoven, "A road to reality with topological superconductors," Nature Physics, vol. 12, no. 7, pp. 618-621, 2016.

[24] M. Asorey, "Space, matter and topology," Nature Physics, vol. 12, no. 7, pp. 616-618, 2016.

[25] N. N. Bogoljubov, "On a new method in the theory of superconductivity," Il Nuovo Cimento Series 10, vol. 7, no. 6, pp. 794-805, 1958.

[26] J. G. Valatin, "Comments on the theory of superconductivity," Il Nuovo Cimento Series 10, vol. 7, no. 6, pp. 843-857, 1958.

[27] G. Rickayzen, "The theory of Bardeen, Cooper and Schrieffer," in Superconductivity, R. D. Parkes, Ed., vol. 1, Marcel Dekker, Inc., New York, NY, USA, 1969.

[28] D. Pines and P. Noziers, The Theory of Quantum Liquids, vol. 1, 2, W A. Benjamin, Inc., New York, NY, USA, 1966.

[29] K. Takayanagi, "Utilizing group property of Bogoliubov transformation," Nuclear Physics A, vol. 808, no. 1-4, pp. 17-26, 2008.

[30] P. Billard and G. Fano, "An existence proof for the gap equation in the superconductivity theory," Communications in Mathematical Physics, vol. 10, p. 274, 1968.

[31] C. Hainzl and R. Seiringer, "The Bardeen-Cooper-Schrieffer functional of superconductivity and its mathematical properties," Journal of Mathematical Physics, vol. 57, no. 2, Article ID 021101, 2016.

[32] A. Deuchert, A. Geisinger, C. Hainzl, and M. Loss, "Persistence of translational symmetry in the BCS model with radial pair interaction," https://arxiv.org/abs/1612.03303.

[33] L. P. Gor'kov, "On the energy spectrum of superconductors," Soviet Phys. JETP, vol. 34, p. 505, 1958.

[34] L. D. Landau and E. M. Lifshitz, Statistical Physics, Part 1, Pergamon Press Ltd., 1980.

[35] A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, Methods ofquantum field theory in statistical physics, Dover Publications, Inc., New York, NY, USA, 1975.

[36] A. Zagoskin, Quantum theory of many-body systems, Graduate Texts in Physics, Springer, 2nd edition, 2014.

Enore Guadagnini

Dipartimento di Fisica "E. Fermi", Universita di Pisa and INFN Sezione di Pisa, Largo B. Pontecorvo 2, 56127Pisa, Italy

Correspondence should be addressed to Enore Guadagnini; enore.guadagnini@unipi.it

Received 18 April 2017; Revised 29 May 2017; Accepted 6 June 2017; Published 10 July 2017

Academic Editor: Xavier Leoncini

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Guadagnini, Enore |

Publication: | Advances in Mathematical Physics |

Date: | Jan 1, 2017 |

Words: | 6030 |

Previous Article: | Stability Monitoring of Batch Processes with Iterative Learning Control. |

Next Article: | Intrinsic Optimal Control for Mechanical Systems on Lie Group. |

Topics: |