Abstract

The thermodynamic potentials of superconducting electrons are derived by means of the Bogoliubov-Valatin formalism. The thermodynamic potentials can be obtained by computing the free energy of a gas of quasi-particles, whose energy spectrum is conditional on the gap function. However, the nontrivial dependence of the gap on the temperature jeopardises the validity of the standard thermodynamic relations. In this article, it is shown how the thermodynamic consistency (i.e., the validity of the Maxwell relations) is recovered, and the correction terms to the quasi-particles potentials are computed. It is shown that the Bogoliubov-Valatin transformation avoids the problem of the thermodynamic consistency of the quasi-particle approach; in fact, the correct identification of the variables, which are associated with the quasi-particles, leads to a precise calculation of the quasi-particles vacuum energy and of the dependence of the chemical potential on the electron density. The stationarity condition for the grand potential coincides with the gap equation, which guarantees the thermodynamic consistency. The expressions of various thermodynamic potentials, as functions of the variables, are produced in the low temperature limit; as a final check, a rederivation of the condensation energy is presented.

1. Introduction

The 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 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 be the energy of a single quasi-particle in the state labelled by the value of the momentum. The function 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 where denotes the chemical potential for the electrons. The gap function can be different from zero in a neighbourhood of the surface of the Fermi sphere and satisfies in which indicates the Debye frequency of the ions lattice. The gap is not vanishing when the temperature takes values below the critical value , and for . As a result, the energy of a single quasi-particle nontrivially depends on the temperature, . But when the dispersion relation 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 represents a standard thermodynamic potential which is computed in statistical mechanics when the atomic energy does not depend on the temperature, let us denote by the corresponding potential for the gas of quasi-particles which is obtained from by means of the substitution . For instance, the free energy and the internal energy of a gas of noninteracting quasi-particles with energy , satisfying the Fermi-Dirac statistics, are given bywhere denotes the spin degeneration factor. When the thermodynamic relations are satisfied, the entropy is given by ; and since , one has . However, if , one finds 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 of the quasi-particles. The presence of adjusting terms is needed in order to rectify the expressions and make the final combinations correct. As far as the computation of is concerned, the knowledge of the quasi-particles distribution is of no help, because 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 . Then the expressions of the thermodynamically consistent potentials as functions of the 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 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 [812]. 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 of the number of electrons is a rather nontrivial issue because when , the 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 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 [1524].

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 , 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 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 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 . Thus, in addition to the temperature, the failure of the thermodynamic consistency for the quasi-particle potentials really concerns several state variables. In certain circumstances the dependence of on these variables is weak enough so that, in a limited range of variability, one can assume [28] that only depends on the momentum. But in general—and in particular when a phase transition occurs—one cannot neglect the dependence of on the macroscopic variables; this dependence will be denoted by .

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 :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 of the quasi-particles. Indeed, in the quantum case of Fermi statistics, the analogue of the Boltzmann -functional is given [2] by Let us recall that, in the case of particles with vanishing chemical potential, the entropy represents the thermodynamic potential for the variables . This means that, with fixed , the stable thermal state of the system corresponds to a maximum of . Therefore the distribution can be determined by the requirement that, for fixed , the variation of with respect to a generic fluctuation must vanish. The associated (conditioned) variational principle takes the form where the value of the Lagrange multiplier is given by as a consequence of the Maxwell relation . Because of the explicit dependence (6) of on , there is no need to introduce a Lagrange multiplier for the volume. Since a modification of causes the following change in the energy the solution of (7) is given precisely by the Fermi-Dirac distribution By using (9) several variables can be computed. Unfortunately, the knowledge of the distribution (9)—or of —alone is not enough to determine all the potentials; for instance, the complete expressions of the free energy and of the internal energy 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 where and denote the creation and annihilation operators for one electron in the state , where refers to the value of one component of the spin: The interaction kernel is related to the amplitude of the electron-electron scattering. does not depend on the values of the electron thermodynamic variables, as the temperature , the volume , and the chemical potential . The main effects of the interactions between electrons in the superconducting materials—which are relevant for the superconducting phase transition—are found when and , and in this region can be approximated by , where is a positive constant. The grand partition function , which is defined bycan be interpreted as the partition function of a system with total Hamiltonian in which the “effective” kinetic energy of one electron (where ) is given by 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 In order to simplify the exposition, let us assume that (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 and can be written as the sum of their means values plus a fluctuation term: 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 where 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 and which create and annihilate quasi-particles in the case of free electrons.(2)Introduction of the operators and —which are obtained from and by means of a Bogoliubov transformation—that diagonalize Hamiltonian (17).

Step 1. The operators and should describe the one-particle fluctuations with respect to the configuration of the Fermi sphere. Then and are chosen in such a way that their vacuum state coincides with the free electrons ground state corresponding to the Fermi sphere. In the limit, the electron states with are occupied whereas the states with are empty, whereTherefore and are defined [25] by The creation operator of a quasi-particle corresponds to the creation of a one-electron hole in the Fermi sphere if and the creation of a one-electron occupied state if . As a result of the creation of a hole in the Fermi sphere, which corresponds to an electron state with momentum and spin component , the total momentum of the system increases by and the total spin of the system increases by . For this reason, when the quantum numbers of the operators and are opposite with respect to the quantum numbers of and . The operators and satisfy the canonical anticommutation relations The operator number of quasi-particles is given by ; 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 limit, the number of quasi-particles vanishes. By using the relation and definition (20), one finds in whichThe 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 of one electron, the kinetic energy of one quasi-particle satisfies .

Step 2. By using the operators and , Hamiltonian (17) reads Since , one can writeThus operator (24) can be written in the form whereBy means of the diagonalizing matrix that satisfies the relation in which one can introduce new creation and annihilation operators and according toOne then obtainsBy using the operators and , the effective Hamiltonian (26) assumes the diagonal form whereThe operators and verify and they represent the creation and annihilation operators for quasi-particles (the so-called bogolons) with dispersion relation . The grand partition function (12) is given by and then the grand potential for the system of conducting electrons turns out to beThe 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 and 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 takes the place of , the thermodynamic consistency requires that the right sign of the dispersion relation should be when and when . The correct value (33) of the energy 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 [2932].
In order to complete the derivation of , 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 and . Since the quasi-particles created and annihilated by and behave as free particles, the only nontrivial mean value is given by , where the quasi-particle distribution is shown in (9). In this way, by using definition (18), one can obtain the gap equation [3, 4, 27, 28, 33] where denotes the integral in the neighbourhood of the surface of the Fermi sphere which is specified by the condition . 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 , it follows that one can put and then Let us now consider expression (36) and let us assume, for the moment, that the value of the gap is a free undeterminate which is not specified by the gap equation; one has where Let us denote by the neighbourhood of the Fermi surface which is defined by the condition . ThenSince is the thermodynamic potential associated with the variables , for fixed the system tends to minimize . So, in the equilibrium state, the value of the gap can be obtained by imposing the variational condition Sinceone finds Moreover Equations (45)-(46) imply By using the identity the minimum condition (43) turns out to be equivalent to (apart from the trivial solution) which coincides precisely with the gap equation (37). In what follows, the value of the gap satisfying (49) will be simply denoted by .

To sum up, the grand potential for the electrons system takes the form where denotes the free energy for the gas of quasi-particles, with dispersion relation shown in (42) and of (41) represents the vacuum energy of the quasi-particles. The value of the gap corresponds to a minimum of for fixed .

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 is the standard dependence of statistical physics. However, when the dependence of the gap on the variables and has been introduced, , the standard rules of statistical physics have been modified. Now, the dependence of on the variables consists of two parts:(i)The (thermodynamically consistent) explicit dependence on displayed in expressions (40)–(42)(ii)The indirect dependence on through .

Let us distinguish these two possibilities by means of the notation When computing the derivatives of one finds, for instance, But since the value of the gap minimizes , the term vanishes, and then the derivative does not contribute to . Therefore, only the thermodynamically consistent dependence of on really contributes to the derivates. This means that the Maxwell relations associated with are satisfied.

Let us now consider a generic thermodynamic potential which is derived from by means of appropriate Legendre transformations. As far as the variable is concerned, the so-called theorem of small increments [34] takes the form and it states that 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, 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 are not explicitly displayed in literature. In certain cases, the combination is reported; however the term is missing. Consequently, in these cases the partial derivatives of with respect to and are invalidated and thus the deduction of the thermodynamic potentials as functions of the variables becomes problematic. Usually, the explicit derivation of the general dependence of the superconducting thermodynamic potentials on represents a rather laborious task. So, let us now derive from (50) a few thermodynamic potentials as functions of when in the low temperature limit. Let us compute, in particular, how the chemical potential gets modified by the superconducting phase transition and how depends on the electron density. The dependence of on is indeed the fundamental ingredient for the determination of the consistent thermodynamic potentials.

In the 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 and, in a neighbourhood of the minimum , it can be approximated as Accordingly, in the computations of and of 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 asThe grand potential (50) reads In the thermodynamic limit, the statistical mean value of the number of electrons is specified by the relation , which takes the form In the low temperature limit one has and , where denotes the Fermi energy So, the chemical potential can be determined perturbatively in powers of and from (57); to lowest nontrivial orders one finds This relation gives the desired expression for the chemical potential. It should be noted that, in addition to the quasi-particles contributions to , 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 is given by : The entropy turns out to be and it coincides with the entropy of the gas of quasi-particles, as it must be. Finally, the internal energy is given by 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 this energy contribution vanishes in the limit. Since , one findsIn the limit, the value of is specified by the gap equation (49); in particular, when , one obtains [3, 4, 13] By means of relation (65), one getsWithin the considered approximations, the value of the condensation energy is obtained by subtracting the energy of a free fermions gas at from the system energy . Since the energy of a free electrons gas in the Fermi sphere is equal to , the condensation energy turns out to be 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 for the superconducting electrons. The value of the gap corresponds to a minimum of and the gap equation precisely represents the stationarity condition. When the gap equation is satisfied, and the thermodynamic potentials that are derived from satisfy the Maxwell relations. The expressions of the thermodynamic potentials as functions of the 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 has been presented.

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.