# Superconducting Properties of 3D Low-Density Translation-Invariant Bipolaron Gas.

1. IntroductionThe theory of superconductivity is one of the finest and oldest subject matters of condensed matter physics which involves both macroscopic and microscopic theories as well as derivation of macroscopic equations of the theory from microscopic description [1]. In this sense the theory was thought to be basically completed and its further development was to have been concerned with further details and consideration of various special cases.

The situation changed when the high-temperature superconductivity (HTSC) was discovered [2]. Surprisingly, it was found that, in oxide ceramics, the correlation length is some orders of magnitude less than that in traditional metal superconductors while the ratio of the energy gap to the temperature of superconducting transition is much greater [3]. The current status of research can be found in books and reviews [4-18].

Today the main problem in this field is to develop a microscopic theory capable of explaining experimental facts which cannot be accounted for by the standard BCS theory. One might expect that development of such a theory would not affect the macroscopic theory based on phenomenological approach.

With all the variety of modern versions of HTSC microscopic descriptions: phonon, plasmon, spin, exciton, and other mechanisms, the central point of constructing the microscopic theory is the effect of electron pairing (Cooper effect). In what follows such a bosonization of electrons provides the basis for the description of their superconducting condensate.

The phenomenon of pairing, in a broad sense, is considered as arising of bielectron states, while, in a narrow sense, if the description is based on phonon mechanism, it is treated as formation of bipolaron states [19]. For a long time this view was hindered by a large correlation length or the size of Cooper pairs in BCS theory. For the same reason, over a long period, the superconductivity was not viewed as a boson condensate (see footnote at p. 1177 in [20]). A significant reason of this lack of understanding was a standard idea that bipolarons are very compact particles.

The most dramatic illustration is the use of the small-radius bipolaron (SRB) theory to describe HTSC [10, 21, 22]. It implies that a stable bound bipolaron state is formed at one node of the lattice and subsequently such small-radius bipolarons are considered as a gas of charged bosons (as a variant individual SRP are formed and then are considered within BCS of creation of the bosonic states). Despite the elegance of such a picture, its actual realization for HTSC comes up against inextricable difficulties caused by impossibility to meet antagonistic requirements. On the one hand, the constant of electron-phonon interaction (EPI) should be large for bipolaron states of small radius to form. On the other hand, it should be small for the bipolaron mass (on which the superconducting temperature depends [23-28]) to be small too. Obviously, the HTSC theory based on SRP concept which uses any other (nonphonon) interaction mechanism mentioned above will run into the same problems.

Alternatively, in describing HTSC one can believe that the role of a fundamental charged boson particle can be played by large-radius bipolarons (LRB) [30-34]. Historically just this assumption was made by Ogg [30] and Schafroth [35] long before the development of the SRP theory. When viewing Cooper pairs as a peculiar kind of large-radius bipolaron states, one might expect that the LRP theory should be used to solve the HTSC problem.

As pointed out above, the main obstacle to consistent use of the LRP theory for explaining high-temperature superconductivity was an idea that electron pairs are localized in a small region, the constant of electron-phonon coupling should be large, and, as a consequence, the effective mass of electron pairs should be large.

In the light of the latest advances in the theory of LRP and LRB, namely, in view of development of an all-new concept of delocalized polaron and bipolaron states, translation-invariant polarons (TI-polarons) and bipolarons (TI-bipolarons) [36-42], it seems appropriate to consider their role in the HTSC theory in a new angle.

We recall the main results of the theory of TI-polarons and bipolarons obtained in [36-42]. Notice that consideration of just electron-phonon interaction is not essential for the theory and can be generalized to any type of interaction.

In what follows we will deal only with the main points of the theory important for the HTSC theory. The main result of papers [36-42] is construction of delocalized polaron and bipolaron states in the limit of strong electron-phonon interaction. The theory of TI-bipolarons is based on the theory of TI-polarons [36, 37] and retains the validity of basic statements proved for TI-polarons. The chief of them is the theorem of analytic properties of the ground state of a TI-polaron (accordingly TI-bipolaron) depending on the constant of electron-phonon interaction a. The main implication of this statement is the absence of a critical value of the EPI constant [[alpha].sub.c], below which the bipolaron state becomes impossible since it decays into independent polaron states. In other words, if there exists a value of ac, at which the TI-state becomes energetically disadvantageous with respect to its decay into individual polarons, then nothing occurs at this point but for [alpha] < [[alpha].sub.c] and the state becomes metastable. Hence, over the whole range of [alpha] variation we can consider TI-polarons as charged bosons capable of forming a superconducting condensate.

Another important property of TI-bipolarons is the possibility of changing the correlation length over the whole range of [0, [infinity]] depending on the Hamiltonian parameters [39]. Hence, it can be both much larger (as is the case in metals) and much less than the characteristic size between the electrons in an electron gas (as happens with ceramics).

A detailed description of the theory of TI-polarons and bipolarons and description of their various properties is given in review [42].

An outstandingly important property of TI-polarons and bipolarons is the availability of an energy gap between their ground and excited states (Section 3).

The above-indicated characteristics can be used to develop a microscopic HTSC theory on the basis of TI-bipolarons.

The paper is arranged as follows. In Section 2 we take Pekar-Froehlich Hamiltonian for a bipolaron as an initial Hamiltonian. The results of three canonical transformations, such as Heisenberg transformation, Lee-Low-Pines transformation, and that of Bogolyubov-Tyablikov are briefly outlined. Equations determining the TI-bipolaron spectrum are derived.

In Section 3 we analyze solutions of the equations for the TI-bipolaron spectrum. It is shown that the spectrum has a gap separating the ground state of a TI-bipolaron from its excited states which form a quasicontinuous spectrum. The concept of an ideal gas of TI-bipolarons is substantiated.

With the use of the spectrum obtained, in Section 4, we consider thermodynamic characteristics of an ideal gas of TI-bipolarons. For various values of the parameters, namely, phononfrequencies, we calculate the values of critical temperatures of Bose condensation, latent heat of transition into the condensed state, heat capacity, and heat capacity jumps at the point of transition.

In Section 5 we discuss the nature of current states in Bose-condensate of TI-bipolarons. It is shown that the transition from a currentless state to a current one is sharp.

In Section 6 the results obtained are compared with the experiment.

In Section 7 we consider the problems of expanding the theory which would enable one to make a more detailed comparison with experimental data on HTSC materials.

In Section 8 we sum up the results obtained.

2. Pekar-Froehlich Hamiltonian: Canonical Transformations

Following [38-42], in describing bipolarons, we will proceed from Pekar-Froehlich Hamiltonian:

[mathematical expression not reproducible], (1)

where [mathematical expression not reproducible] are coordinates of the first and second electrons, respectively; [a.sup.+.sub.k], [a.sub.k] are operators of the birth and annihilation of the field quanta with energy [??][[omega].sup.0.sub.k] = [??][[omega].sub.0]; [m.sup.*] is the electron effective mass; the quantity U describes Coulomb repulsion between the electrons; [V.sub.k] is the function of the wave vector k:

[mathematical expression not reproducible], (2)

where e is the electron charge; [[epsilon].sub.[infinity]] and [[epsilon].sub.0] are high-frequency and static dielectric permittivities; [alpha] is the constant of electron-phonon interaction; V is the systems volume.

In the system of the center of mass Hamiltonian (1) takes the form:

[mathematical expression not reproducible]. (3)

In what follows in this section we will believe [??] = 1, [[omega].sup.0.sub.k] = 1, [M.sub.e] = 1 (accordingly [[mu].sub.e] = 1/4).

The coordinates of the center of mass [??] can be excluded from Hamiltonian (3) using Heisenberg's canonical transformation [43]:

[mathematical expression not reproducible]. (4)

Accordingly, the transformed Hamiltonian will be written as

[mathematical expression not reproducible]. (5)

From (5) it follows that the exact solution of the bipolaron function is determined by the wave function [psi](r), which contains only relative coordinates r and, therefore, is translation-invariant.

Averaging of [??] over [psi](r) yields the Hamiltonian [bar.H]:

[mathematical expression not reproducible]. (6)

Equation (6) suggests that the bipolaron Hamiltonian differs from the polaron one in that in the latter the quantity [V.sub.k] is replaced by [[bar.V].sub.k] and the constants [[bar.V].sub.T], [[bar.V].sub.U] are added.

With the use of Lee-Low-Pines canonical transformation [44]:

[mathematical expression not reproducible], (7)

where [f.sub.k] are variational parameters having the sense of the distance by which the field oscillators are displaced from their equilibrium positions:

[mathematical expression not reproducible], (8)

for Hamiltonian [??]:

[??] = [S.sup.-1.sub.2] [bar.H][S.sub.2], (9)

we get

[mathematical expression not reproducible] (10)

[mathematical expression not reproducible]. (11)

Hamiltonian [H.sub.1] contains linear, threefold, and fourfold terms in the birth and annihilation operators. Its explicit form is given in [36-38].

Then, as is shown in [36, 37], the use of Bogolyubov-Tyablikov canonical transformation [45] for passing on from operators [a.sup.+.sub.k], [a.sub.k] to new operators [[alpha].sup.+.sub.k], [[alpha].sub.k]:

[mathematical expression not reproducible] (12)

(in which [H.sub.0] is a diagonal operator), makes mathematical expectation of [H.sub.1] equal to zero.

In the new operators [[alpha].sup.+.sub.k], [[alpha].sub.k] Hamiltonian (11) takes on the form [??]:

[mathematical expression not reproducible] (13)

[mathematical expression not reproducible], (14)

where [DELTA][E.sub.r] is the so-called recoil energy. The general expression for [DELTA][E.sub.r] = [DELTA][E.sub.r]{[f.sub.k]} was obtained in [37]. Actually, calculation of the ground state energy [E.sub.bp] was performed in [41] by minimization of (14) in [f.sub.k] and in [psi].

Notice that in the theory of a polaron with broken symmetry a diagonalized electron-phonon Hamiltonian has the form of (13)-(14) [46]. This Hamiltonian can be treated as a Hamiltonian of a polaron and a system of its associated renormalized real phonons or as a Hamiltonian whose quasiparticle excitations spectrum is determined by (13)-(14) [47]. In the latter case excited states of a polaron are Fermi quasiparticles.

In the case of a bipolaron the situation is qualitatively different since a bipolaron is a boson quasiparticle whose spectrum is determined by (13)-(14). Obviously, a gas of such quasiparticles can experience Bose-Einstein condensation (BEC). Treatment of (13)-(14) as a bipolaron and its associated renormalized phonons does not prevent their BEC since maintenance of the number of particles required in this case takes place automatically due to commutation of the total number of renormalized phonons with Hamiltonian (13)-(14).

Renormalized frequencies [v.sub.k] involved in (13)-(14), according to [36, 37], are determined by the equation for s:

[mathematical expression not reproducible], (15)

solutions of which yield the spectrum of s = {[v.sup.2.sub.k]} values.

3. Energy Spectrum of a TI-Bipolaron

Hamiltonian (13)-(14) is conveniently presented in the form:

[mathematical expression not reproducible], (16)

[mathematical expression not reproducible], (17)

where in the case of a three-dimensional ionic crystal is

[mathematical expression not reproducible], (18)

where [mathematical expression not reproducible] is the number of atoms along the ith crystallographic axis.

Let us prove the validity of the expression for the spectrum (16), (17). Since operators [[alpha].sup.+.sub.n], [[alpha].sub.n] obey Bose commutation relations:

[mathematical expression not reproducible], (19)

they can be considered to be operators of birth and annihilation of TI-bipolarons. The energy spectrum of TI-bipolarons, according to (15), is determined by the equation

F(s) = 1, (20)

where

[mathematical expression not reproducible]. (21)

It is convenient to solve (20) graphically (Figure 1).

Figure 1 suggests that the frequencies [mathematical expression not reproducible] (index i is omitted) lie between the frequencies [mathematical expression not reproducible] and [mathematical expression not reproducible]. Hence, the spectrum [mathematical expression not reproducible] as well as the spectrum [mathematical expression not reproducible] are quasicontinuous: [mathematical expression not reproducible], which just proves the validity of (16), (17).

It follows that the spectrum of a TI-bipolaron has a gap between the ground state [E.sub.bp] and the quasicontinuous spectrum, equal to [[omega].sub.0].

Below we will consider the case of low concentration of TI-bipolarons in a crystal. Then they can adequately be considered as an ideal Bose gas, whose properties are determined by Hamiltonian (16).

4. Statistical Thermodynamics of Low-Density TI-Bipolaron Gas

Let us consider an ideal Bose gas of TI-bipolarons which represents a system of N particles occurring in some volume V. Let us write [N.sub.0] for the number of particles in the lower one-particle state and N' for the number of particles in higher states. Then

[mathematical expression not reproducible], (22)

[mathematical expression not reproducible]. (23)

In expression N' (23), we will perform integration over quasicontinuous spectrum (instead of summation) (16), (17) and assume [mu] = [E.sub.bp]. As a result, from (22), (23) we get an equation for determining the critical temperature of Bose condensation [T.sub.c]:

[mathematical expression not reproducible], (24)

[mathematical expression not reproducible], (25)

where n = N/V. Figure 2 shows a graphical solution of (24) for the values of parameters [M.sub.e] = 2[m.sup.*] = 2[m.sub.0], where [m.sub.0] is the mass of a free electron in vacuum, [[omega].sup.*] = 5 meV ([approximately equal to]58 K), n = [10.sup.21] [cm.sup.-3], and the values [mathematical expression not reproducible].

It is seen from Figure 2 that the critical temperature grows with increasing phonon frequency [[omega].sub.0]. The relations of critical temperatures [T.sub.ci]/[[omega].sub.0i] corresponding to the chosen parameter values are given in Table 1. Table 1 suggests that the critical temperature of a TI-bipolaron gas is always higher than that of ideal Bose gas (IBG). It is also evident from Figure 2 that an increase in the concentration of TI-bipolarons n will lead to an increase in the critical temperature, while a gain in the electron mass [m.sup.*] to its decrease. For [??] = 0 the results go over into the limit of IBG. In particular, (24), for [??] = 0, yields the expression for the critical temperature of IBG:

[mathematical expression not reproducible]. (26)

It should be stressed, however, that (26) involves [M.sub.e] = 2[m.sup.*], rather than the bipolaron mass. This resolves the problem of the low temperature of condensation which arises both in the SRP theory and in the LRP theory in which expression (26) involves the bipolaron mass [31-34]. Another important result is that the critical temperature [T.sub.c] for the parameter values considerably exceeds the gap energy [[omega].sub.0].

From (22), (23), it follows that

[mathematical expression not reproducible] (27)

Figure 3 shows temperature dependencies of the number of supracondensate particles N' and the number of particles [N.sub.0] occurring in the condensate for the above-indicated parameter values [[??].sub.i].

Figure 3 suggests that, as could be expected, the number of particles in the condensate grows as the gap [[omega].sub.i] increases.

The energy of a TI-bipolaron gas E is determined by the expression

[mathematical expression not reproducible]. (28)

With the use of (16), (17), and (28) the specific energy (i.e., the energy per one TI-bipolaron) [mathematical expression not reproducible] will be

[mathematical expression not reproducible], (29)

[mathematical expression not reproducible], (30)

where [mu] is determined from the equation

[mathematical expression not reproducible], (31)

[mathematical expression not reproducible]. (32)

Relation of [??] with the chemical potential of the system [mu] is given by the formula [??] = ([mu] - [E.sub.bp])/[[omega].sup.*]. From (29)- (31) expressions for the free energy [DELTA]F = -(2/3) [DELTA]E, [DELTA]F = F - [E.sub.bp]N, [DELTA]E = E - [E.sub.bp]N and entropy S = -[partial derivative]F/[partial derivative]T also follow.

Figure 4 illustrates temperature dependencies [DELTA]E = [??] - [[??].sub.bp] for the above-indicated parameter values [[omega].sub.i]. Break points on the curves [DELTA][E.sub.i]([??]) correspond to the values of critical temperatures [mathematical expression not reproducible].

The dependencies obtained enable us to find the heat capacity of a TI-bipolaron gas: [mathematical expression not reproducible]. With the use of (29) for [??] [less than or equal to] [[??].sub.c], we express [C.sub.v]([??]) as

[mathematical expression not reproducible], (33)

[mathematical expression not reproducible]. (34)

Expression (33) yields a well-known exponential dependence of the heat capacity at low temperatures [C.sub.v] ~ exp(-[[omega].sub.0]/T), caused by the availability of the energy gap [[omega].sub.0].

Figure 5 shows temperature dependencies of the heat capacity [C.sub.v]([??]) for the above-indicated parameter values [[??].sub.i]. Table 1 lists the values of the heat capacity jumps:

[mathematical expression not reproducible] (35)

at the transition points for the parameter values [[??].sub.i].

The dependencies obtained will enable one to find the latent heat of transition q = TS, where S is the entropy of supracondensate particles. At the point of transition this value is q = 2[T.sub.c][C.sub.v]([T.sub.c] - 0)/3, where [C.sub.v](T) is determined by formula (33). For the above-indicated parameter values [[omega].sub.i], it is given in Table 1.

5. Current States of a TI-Bipolaron Gas

In the foregoing we have considered equilibrium properties of a TI-bipolaron gas. The formation of Bose-condensate per se does not mean that it has superconducting properties. To demonstrate such a possibility let us consider the total momentum of a TI-bipolaron:

[mathematical expression not reproducible], (36)

where [[??].sub.1] and [[??].sub.2] are the momenta of the first and second electron, respectively. It is easy to check that [??] commutes with Hamiltonian (1) and therefore is a constant value, that is, c-number.

For this reason, to consider nonequilibrium properties and, particularly current states, we can use generalization of Heisenberg transformation (4) such that

[mathematical expression not reproducible]. (37)

The general expression for the functional of the total energy of a TI-bipolaron for [??] [not equal to] 0 is given in [36]. TI-bipolarons occurring in condensed state have a common wave function for the_whole condensate and do not thermalize in the state when [??] [not equal to] 0. In the supracondensate part, phonons whose wave vectors contribute to [??] (36) are thermalized and their total momentum will be equal to zero.

Hence, all the changes arising in considering the case of [??] [not equal to] 0 concern only the expression for the ground state energy which becomes dependent on [??], while the spectrum of excited states (16), (17) remains unchanged. It follows that, in passing over the critical point, the currentless state suddenly becomes current which is in agreement with the experiment.

6. Comparison with the Experiment

Figure 4 shows typical dependencies of E([??]). They suggest that at the point of transition the energy is a continuous function of [??]. This means that the transition per se occurs without energy expenditure being a phase transition of the second kind in complete agreement with the experiment. At the same time transition of Bose particles from a condensate state to a supracondensate one occurs with consumption of energy which is determined by the value q (Section 4, Table 1), determining the latent heat of transition of a Bose gas which makes it a phase transition of the 1st kind.

Byway of example let us consider HTSC Y[Ba.sub.2][Cu.sub.3][O.sub.7] with the temperature of transition 90 / 93K, volume of the unit cell 0, 1734 x [10.sup.-21] [cm.sup.-3], and concentration of holes n [approximately equal to] [10.sup.21] [cm.sup.-3]. According to estimates [48], Fermi energy is equal to [[member of].sub.F] = 0,37 eV. Concentration of TI-bipolarons in Y[Ba.sub.2][Cu.sub.3][O.sub.7] is found from (24):

[mathematical expression not reproducible] (38)

with [[??].sub.c] = 1, 6. Table 1 lists the values of [n.sub.bp,i] for the values of parameters [[??].sub.i] given in Section 4. It follows from Table 1 that [n.sub.bp,i] [much less than] n. Hence, only a small part of charge carriers is in a bipolaron state which justifies the approximation of a low-density TI-bipolaron gas used by us. This fact correlates with recent experiments [49] where it was shown that only a small part of electrons are in SC state. The energy levels of such TI-bipolarons lie near Fermi surface and are described by the wave function:

[mathematical expression not reproducible], (39)

which leads to replacement of [bar.T], involved in (6), by

[bar.T] = -2<[phi][absolute value of [[DELTA].sub.r]] [phi]> + 2[k.sup.2.sub.F], (40)

that is, reckoning of the energy from Fermi level (the last term in the right-hand side of (40) in dimensional units is equal to 2[[member of].sub.F], where [[member of].sub.F] = [h.sup.2][k.sup.2.sub.F]/2[m.sup.*]).

According to our approach, superconductivity arises when coupled states are formed. The condition for the formation of such states has the form

[E.sub.bp] < 0, (41)

where [E.sub.p] is the energy of a TI-polaron [41]. Condition (41) determines the value of a pseudogap:

[[DELTA].sub.1] = [absolute value of [E.sub.bp]]. (42)

For [[DELTA].sub.1] [much greater than] [[omega].sub.0] the value of a pseudogap can greatly exceed both [T.sub.c] and the energy of the gap (i.e., [[omega].sub.0]). The expression for the spectrums [E.sub.bp] and [E.sub.p] (16)-(17) suggests that the angular dependence of superconducting gap is completely determined by the symmetry of the isoenergetic surface of the [[omega].sub.0] ([??]). Earlier this conclusion was made by Bennet [50], who proved that the main source of anisotropy of superconducting properties is the angular dependence of the phonon spectrum, though some contribution is also made by the anisotropy of Fermi surface.

It follows from what has been said that formation of a pseudogap is a phase transition preceding the phase transition to the superconducting state. Recent experiments [51] also testify in favor of this statement.

In paper [39] correlation length for TI-bipolarons was calculated. According to [39], in HTSC materials its value can vary from several angstroms to several tens of angstroms, which is also in agreement with the experiment.

Of special interest is to determine the characteristic energy of phonons responsible for the formation of TI-bipolarons and superconducting properties of oxide ceramics. To do this let us compare the calculated values of the heat capacity jumps with experimental data.

As is known, in BCS theory, a jump in the heat capacity is equal to

[[C.sub.s] - [C.sub.n]/[C.sub.n]]|[T.sub.c] = 1, 43, (43)

where [C.sub.s] is the heat capacity in the superconducting phase and [C.sub.n] is the heat capacity in the normal phase and is independent of the parameters of the model Hamiltonian. As it follows from numerical calculations shown in Figure 5 and in Table 1, as distinct from the BCS theory, the value of the jump depends on the phonon frequency. Hence, the approach presented predicts the existence of an isotopy effect for the heat capacity jump.

As it is seen from Figure 6, the heat capacity jump calculated theoretically (Section 4) coincides with the experimental value in Y[Ba.sub.2][Cu.sub.3][O.sub.7] [29], for [??] = 1, 5, that is, for [omega] = 7, 5 meV. This corresponds to the concentration of TI-bipolarons equal to [n.sub.bp] = 2, 6 x [10.sup.18] [cm.sup.-3]. Hence, in contrast to the widespread notion that in oxide ceramics superconductivity is determined by high-energyphonons (with energy 70 / 80 meV [52]), actually, the superconductivity in HTSC materials should be determined by soft phonon modes.

Notice that in calculations of the temperature of transition it was believed that the effective mass [M.sub.e] in (24) is independent of the direction of the wave vector; that is, isotropic case was dealt with.

In the anisotropic case, choosing principal axes of vector k, as coordinate axes, we will get the quantity [([M.sub.ex][M.sub.ey][M.sub.ez]).sup.1/3] instead of the effective mass [M.sub.e]. In complex HTSC materials the values of effective masses lying in the plane of layers [M.sub.ex], [M.sub.ey] are close in value. Assuming in this case [M.sub.e] = [M.sub.ex] = [M.sub.ey] = [M.sub.[parallel]], [M.sub.ez] = [M.sub.[parallel]], we will get instead of [C.sub.bp], determined by (24), the value [[??].sub.bp] = [C.sub.bp]/[gamma]; [[gamma].sup.2] = [M.sub.[perpendicular to]]/[M.sub.[parallel]] is the parameter of anisotropy. Hence anisotropy of effective masses gives for the concentration [n.sub.bp] the value [[??].sub.bp] = [gamma][n.sub.bp]. Therefore taking account of anisotropy can, in an order of magnitude, enhance the estimate of the concentration of TI-bipolarons. If for Y[Ba.sub.2][Cu.sub.3][O.sub.7] we take the estimate [[gamma].sub.2] = 30 [52], then for the concentration of TI- bipolarons we will get [[??].sub.bp] = 1, 4 x [10.sup.19] [cm.sup.-3], which holds valid the general conclusion: in the case under consideration only a small number of charge carriers are in TI-bipolaron state. The situation can change if the anisotropy parameter is very large. Thus, for example, in layered HTSC Bi-Sr-Ca-Cu-O the anisotropy parameter is [gamma] > 100; accordingly, the concentration of TI-bipolarons in these compounds can have the same order of magnitude as the total concentration of charge carriers.

Another important conclusion emerging from taking account of the anisotropy of effective masses is that the temperature of the transition [T.sub.c] depends not on [n.sub.bp] and [M.sub.[parallel]] individually, but on their relation which straightforwardly follows from (24).

7. Essential Generalizations of the Theory

In the foregoing we considered the case of an ideal TI-bipolaron gas. At small concentration of TI-bipolarons their Coulomb interaction will be greatly screened which justifies the use of the model of an ideal gas.

If the concentration of TI-bipolarons is large (e.g., n = [10.sup.21] [cm.sup.-3], as was believed in Section 4), then such Bose gas can no longer be considered to be ideal. Taking account of Coulomb interaction between bipolarons becomes necessary. For this purpose we can use Bogolyubov theory [53] for weakly imperfect Bose gas, which implies that the spectrum of elementary excitations with the momentum k will be determined by the expression:

[mathematical expression not reproducible], (44)

where u(k) = [square root of [n.sub.B]V(k)/[m.sub.B]], V(k) is the Fourier component of the potential of pairwise interaction between charged bosons: V(k) = 4[pi][e.sup.2.sub.B]/[[member of].sub.0][k.sup.2], [e.sub.B] and [m.sub.B] are the charge and mass of a boson, and [n.sub.B] is the concentration of bosons.

Hence it follows that

[mathematical expression not reproducible]. (45)

[[omega].sub.p] is the plasma frequency. According to (45), the spectrum of excitations of quasiparticles of an imperfect gas is characterized by a finite energy gap which in the long-wavelength limit is equal to [[omega].sub.p]. The same result can be arrived at if the Coulomb interaction between the electrons is considered as a result of electron-plasmon interaction. According to [54], Hamiltonian of electron-plasmon interaction coincides in structure with Froehlich Hamiltonian which involves plasmon frequency instead of phonon frequency. This straightforwardly leads to the energy gap equal to the plasmon frequency [[omega].sub.p], if in (45) we put [e.sub.B] = 2e, [m.sub.B] = 2[m.sup.*], [n.sub.B] = n/2, where n is the electron concentration. Hence, for [[omega].sub.p] < [[omega].sub.0], the energy gap will be determined by [[omega].sub.p], while for [[omega].sub.p] > [[omega].sub.0] it will still be determined by [[omega].sub.0].

Actually, in real HTSC, there are not only phonon and plasmon branches, but also some other elementary excitations which can take part in electron pairing. An example is spin fluctuations. In generalization of the theory to the case of interaction with various brunches of excitations which will contribute to the ground state energy, the value of the gap and the dispersion law of a TI-bipolaron are a topical problem. Obviously, taking account of this interaction will lead to an increase in the coupling energy of both TI-bipolarons and TI-polarons. Therefore, a priori, without any particular calculations one cannot say anything of how the condition of stability of TI-bipolaron states determined by (41) will change.

Another important problem is generalization of the theory to the case of intermediate coupling of electron-phonon interaction. Formally, the expression for the functional of the ground state energy of a TI-bipolaron (14) is valid for any value of the electron-phonon coupling constant. For this reason such a calculation will not change the spectrum of a TI-bipolaron; however it will alter the criteria of fulfillment of the conditions of a TI-bipolaron stability (41).

8. Conclusive Remarks

In this paper we have presented conclusions emerging from consistent translation-invariant consideration of EPI. It implies that pairing of electrons, for any coupling constant, leads to a concept of TI-polarons and TI-bipolarons. Being bosons, TI-bipolarons can experience Bose condensation leading to superconductivity. Let us list the main results following from this approach. First and foremost the theory resolves the problem of the great value of the bipolaron effective mass (Section 4). As a consequence, formal limitations on the value of the critical temperature of the transition are eliminated too. The theory quantitatively explains such thermodynamic properties of HTSC-conductors as availability (Section 4) and value (Section 6) of the jump in the heat capacity lacking in the theory of Bose condensation of an ideal gas. The theory also gives an insight into the occurrence of a great ratio between the width of the pseudogap and [T.sub.c] (Section 6). It accounts for the small value of the correlation length [39] and explains the availability of a gap (Section 3) and a pseudogap (Section 6) in HTSC materials. The angular dependence of the gap and pseudogap gets a natural explanation (Section 6).

Accordingly, isotopic effect automatically follows from expression (24), where the phonon frequency [[omega].sub.0] acts as a gap. The conclusion of the dependence of the temperature of the transition [T.sub.c] on the relation [n.sub.hp]/[M.sub.[parallel]] (Section 6) correlates with Uemura law universal for all HTSC materials which implies that the temperature of the transition is related to the concentration of charge carriers divided by its effective mass [55, 56]. The TI-bipolaron theory of superconductivity developed in this paper gives an answer to the question of paper [49], namely, where most of electrons disappear in the superconductors analyzed. It lies in the fact that only a small part of electrons are paired. The results obtained suggest that in order to raise the critical temperature [T.sub.c] one should increase the concentration of bipolarons.

Application of the theory to 1D and 2D systems leads to qualitatively new results since the occurrence of a gap in the TI-bipolaron spectrum automatically removes divergences at small momenta, inherent in the theory of ideal Bose gas. These problems were considered by the author in [57]. The case of 1D-polaron was discussed in [58].

https://doi.org/10.1155/2018/1380986

Conflicts of Interest

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

Acknowledgments

The work was done with the support from the Russian Foundation for Basic Research (RFBR Project no. 16-0700305).

References

[1] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2: Theory of the Condensed State, Butterworth-Heinemann, 1st edition, 1980.

[2] J. G. Bednorz and K. A. Muller, "Possible high [T.sub.c] superconductivity in the Ba-La-Cu-O system," Zeitschrift fur Physik B, vol. 64, pp. 189-193, 1986.

[3] V. L. Ginzburg, "Superconductivity: the day before yesterday--yesterday --today--tomorrow," Physics-Uspekhi, vol. 43, p. 573, 2000.

[4] S. L. Kakani and S. Kakani, Superconductivity, Anshan, China, 2009.

[5] T. Tohyama, "Recent progress in physics of high-temperature superconductors," Japanese Journal of Applied Physics, vol. 51, Article ID 010004, 2012.

[6] S. Kruchinin, H. Nagao, and S. Aono, Modern Aspects of superconductivity. Theory of superconductivity, World Scientific Publishing Co. Pte. Ltd., River Edge, NJ, USA, 2011.

[7] K. P. Sinha and S. L. Kakani, "Fermion local charged boson model and cuprate superconductors," Proceedings--National Academy of Sciences, India. Section A, Physical Sciences, vol. 72, p. 153, 2002.

[8] K. H. Benneman and J. B. Ketterson, Superconductivity: Conventional and Unconventional Superconductors 1-2, Springer, New York, NY, UK, 2008.

[9] J. R. Schrieffer, Theory of Superconductivity, Westview Press, Oxford, UK, 1999.

[10] A. S. Alexandrov, Theory of Superconductivity from weak to strong oupling, IOP, publishing, Bristol, UK, 2003.

[11] N. M. Plakida, "High temperature cuprate superconductors: experiment," in Theory and Applications, vol. 166 of Springer Series in Solid-State Sciences, pp. 479-494, Springer, Berlin, Germany, 2010.

[12] I. Askerzade, "Physical Properties of Unconventional Superconductors," in Unconventional Superconductors, vol. 153 of Springer Series in Materials Science, pp. 1-26, Springer Berlin Heidelberg, Berlin, Germany, 2012.

[13] O. Gunnarsson and O. Rosch, "Interplay between electron-phonon and Coulomb interactions in cuprates," Journal of Physics, vol. 20, Article ID 043201, 2008.

[14] T. Moriya and K. Ueda, "Spin fluctuations and high temperature superconductivity," Advances in Physics, vol. 49, no. 5, pp. 556-606, 2000.

[15] D. Manske, Theory of Unconventional Superconductors, Springer, Heidelberg, Germany, 2004.

[16] R. Szczesniak and A. P. Durajski, "Anisotropy of the gap parameter in the hole-doped cuprates," Superconductor Science and Technology, vol. 27, no. 12, Article ID 125004, 2014.

[17] A. P. Durajski, "Anisotropic evolution of energy gap in Bi2212 superconductor," Frontiers of Physics, vol. 11, no. 5, Article ID 117408, 2016.

[18] R. Szczesniak, A. P. Durajski, and A. M. Duda, "Pseudogap in Eliashberg approach based on electron-phonon and electron-electron-phonon interaction," Annalen der Physik, vol. 529, Article ID 1600254, 2017.

[19] V. D. Lakhno, "Phonon interaction of electrons in the translation-invariant strong-coupling theory," Modern Physics Letters B, vol. 30, Article ID 1650031, 2016.

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

[21] A. S. Alexandrov and N. Mott, Polarons Bipolarons, World Sci, Pub. CO Inc, Singapore, Singapore, 1996.

[22] A. S. Alexandrov and A. B. Krebs, "Reviews of topical problems," Soviet Physics--Uspekhi, vol. 35, no. 5, pp. 345-383, 1992.

[23] P. W. Anderson, "Condensed matter: the continuous revolution," Physics World, vol. 8, no. 12, p. 37, 1995.

[24] N. F. Mott, "High-temperature superconductivity debate heats up," Physics World, vol. 9, no. 1, p. 16, 1996.

[25] P. W. Anderson, Physics World, vol. 9, no. 1, p. 16, 1996.

[26] E. V. de Mello and J. Ranninger, "Dynamical properties of small polarons," Physical Review B: Condensed Matter and Materials Physics, vol. 55, no. 22, pp. 14872-14885, 1997.

[27] Y. A. Firsov, V. V. Kabanov, E. K. Kudinov, and A. S. Alexandrov, "Comment on "Dynamical properties of small polarons"," Physical Review B: Condensed Matter and Materials Physics, vol. 59, p. 12132, 1999.

[28] E. V. de Mello and J. Ranninger, "Reply to "Comment on 'Dynamical properties of small polarons' "," Physical Review B: Condensed Matter and Materials Physics, vol. 59, no. 18, pp. 12135-12136, 1999.

[29] N. Overend, M. A. Howson, and I. D. Lawrie, "3D X-Y scaling of the specific heat of Y[Ba.sub.2][Cu.sub.3][O.sub.7-[delta]] single crystals," Physical Review Letters, vol. 72, p. 3238, 1994.

[30] R. A. Ogg, "Superconductivity in solid metal-ammonia solutions," Physical Review, vol. 70, p. 93, 1946.

[31] V. L. Vinetskii and E. A. Pashitskii, "Superfluidity of charged Bose-gas and bipolaron mechanism of superconductivity," Ukrainian Journal of Physics, vol. 20, p. 338, 1975.

[32] E. A. Pashitskii and V. L. Vinetskii, "Plasmon and bipolaron mechanisms of high-temperature superconductivity," JETP Letters, vol. 46, p. S104, 1987.

[33] D. Emin, "Formation, motion, and high-temperature superconductivity of large bipolarons," Physical Review Letters, vol. 62, no. 13, pp. 1544-1547, 1989.

[34] V. L. Vinetskii, N. I. Kashirina, and E. A. Pashitskii, "Bipolaron states in ion crystals and the problem of high temperature superconductivity," Ukrainian Journal of Physics, vol. 37, no. 76, 1992.

[35] M. R. Schafroth, "Superconductivity of a charged boson gas," Physical Review, vol. 96, no. 4, p. 1149, 1954.

[36] A. V. Tulub, "Recoil accounting in nonrelativistic quantum field theory," Vestnik Leningrad University, vol. 15, no. 22, p. 104, 1960.

[37] A. V. Tulub, "Slow electrons in polar crystals," Soviet Physics-JETP, vol. 14, no. 1301, 1962.

[38] V. D. Lakhno, "Energy and critical ionic-bond parameter of a 3D large-radius bipolaron," Journal of Experimental and Theoretical Physics, vol. 110, no. 5, pp. 811-815, 2010.

[39] V. D. Lakhno, "Translation-invariant bipolarons and the problem of high-temperature superconductivity," Solid State Communications, vol. 152, p. 621, 2012.

[40] N. I. Kashirina, V. D. Lakhno, and A. V. Tulub, "The virial theorem and the ground state problem in polaron theory," Journal of Experimental and Theoretical Physics, vol. 114, no. 5, pp. 867-869, 2012.

[41] V. D. Lakhno, "Translation invariant theory of polaron (bipolaron) and the problem of quantizing near the classical solution," Journal of Experimental and Theoretical Physics, vol. 116, no. 6, pp. 892-896, 2013.

[42] V. D. Lakhno, "Pekar's ansatz and the strong coupling problem in polaron theory," Physics-Uspekhi, vol. 58, p. 295, 2015.

[43] W. Heisenberg, "Die Selbstenergie des Elektrons," Zeitschriftfur Physik, vol. 65, no. 1-2, pp. 4-13, 1930.

[44] T. D. Lee, F. E. Low, and D. Pines, "The motion of slow electrons in a polar crystal," Physical Review, vol. 90, no. 2, pp. 297-302, 1953.

[45] S. V. Tyablikov, Methods of Quantum Theory of Magnetizm, Moscow, Russia, 1975.

[46] S. J. Miyake, in Polarons and Applications, V. D. Lakhno, Ed., Wiley, Chichester, 1994.

[47] I. B. Levinson and E. I. Rashba, "Threshold phenomena and bound states in the polaron problem," Soviet Physics--Uspekhi, vol. 16, no. 6, pp. 892-912, 1974.

[48] L. P. Gor'kov and N. B. Kopnin, "High-Tcsuperconductors from the experimental point of view," Soviet Physics--Uspekhi, vol. 31, no. 9, pp. 850-860, 1988.

[49] I. Bozovic, X. He, J. Wu, and A. T. Bollinger, "Dependence of the critical temperature in overdoped copper oxides on superfluid density," Nature, vol. 536, no. 7616, pp. 309-311, 2016.

[50] A. J. Bennett, "Theory of the anisotropic energy gap in superconducting lead," Physical Review, vol. 140, p. A1902, 1965.

[51] R.-H. He, M. Hashimoto, H. Karapetyan et al., "From a single-band metal to a high-temperature superconductor via two thermal phase transitions," Science, vol. 331, p. 1579, 2011.

[52] A. Marouchkine, Room-Temperature Superconductivity, vol. 7, Cambridge Int. Sci. Publ., Cambridge, UK, 2004.

[53] N. N. Bogolyubov, "On the theory of superfluidity," Journal of Physics USSR, vol. 11, p. 77, 1947.

[54] V. D. Lakhno, In Polarons and Applications, V D. Lakhno, Ed., Wiley, Chichester, 1994.

[55] Y. J. Uemura, G. M. Luke, B. J. Sternlieb et al., "Universal correlations between [T.sub.c] and [n.sub.s]/m (carrier density over effective mass) in High-[T.sub.c] cuprate superconductors," Physical Review Letters, vol. 62, p. 2317, 1989.

[56] Y. J. Uemura, L. P. Le, G. M. Luke et al., "Basic similarities among cuprate, bismuthate, organic, Chevrel-phase, and heavy-fermion superconductors shown by penetration-depth measurements," Physical Review Letters, vol. 66, p. 2665, 1991.

[57] V. D. Lakhno, "A translation invariant bipolaron in the Holstein model and superconductivity," SpringerPlus, vol. 5, p. 1277, 2016.

[58] V. D. Lakhno, "Large-radius Holstein polaron and the problem of spontaneous symmetry breaking," Progress of Theoretical and Experimental Physics, vol. 2014, no. 7, Article ID 073I01, 2014.

V. D. Lakhno (ID)

Keldysh Institute of Applied Mathematics, RAS, Miusskaya Sq. 4, Moscow 125047, Russia

Correspondence should be addressed to V. D. Lakhno; lak@impb.ru

Received 12 December 2017; Accepted 28 January 2018; Published 15 April 2018

Academic Editor: Jorg Fink

Caption: Figure 1: Graphical solution of (20).

Caption: Figure 2: Solutions of (24) with [C.sub.bp] = 331, 3 and [[??].sub.i] = {0, 2; 1; 2; 10; 15; 20}, which correspond to [mathematical expression not reproducible].

Caption: Figure 3: Temperature dependencies of the relative number of supracondensate particles N'/N and the particles occurring in the condensate [N.sub.0]/N = 1 - N'/N for the parameter values [[??].sub.i], given in Figure 2.

Caption: Figure 4: Temperature dependencies [mathematical expression not reproducible] for the parameter values [[??].sub.i] presented in Figures 2 and 3.

Caption: Figure 5: Temperature dependencies of the heat capacity for various values of the parameters [mathematical expression not reproducible].

Caption: Figure 6: Comparison of the theoretical (solid line) and experimental [29] (broken line) dependencies in the region of the heat capacity jump.

Table 1: Calculated characteristics of Bose gas of TI-bipolarons with concentration n = [10.sup.21] [cm.sup.-3]. i 0 1 [[??].sub.i] 0 0,2 [T.sub.ci]/[[omega].sub.oi] 136,6 [q.sub.i]/[T.sub.ci] 1,3 1,44 -[DELTA]([partial derivative] [C.sub.v,i]/[partial derivative][??]) 0,11 0,12 [C.sub.v,i] ([T.sub.c] - 0) 1,9 2,16 ([C.sub.s] - [C.sub.n])/ [C.sub.n] 0 0,16 [n.sub.bpi] x [cm.sup.3] 16 x [10.sup.19] 9,4 x [10.sup.18] i 2 3 [[??].sub.i] 1 2 [T.sub.ci]/[[omega].sub.oi] 30 16 [q.sub.i]/[T.sub.ci] 1,64 1,8 -[DELTA]([partial derivative] [C.sub.v,i]/[partial derivative][??]) 0,12 0,13 [C.sub.v,i] ([T.sub.c] - 0) 2,46 2,7 ([C.sub.s] - [C.sub.n])/ [C.sub.n] 0,36 0,52 [n.sub.bpi] x [cm.sup.3] 4,2 x [10.sup.18] 2,0 x [10.sup.18] i 4 5 [[??].sub.i] 10 15 [T.sub.ci]/[[omega].sub.oi] 4,2 3 [q.sub.i]/[T.sub.ci] 2,5 2,8 -[DELTA]([partial derivative] [C.sub.v,i]/[partial derivative][??]) 0,14 0,15 [C.sub.v,i] ([T.sub.c] - 0) 3,74 4,2 ([C.sub.s] - [C.sub.n])/ [C.sub.n] 1,23 1,53 [n.sub.bpi] x [cm.sup.3] 1,2 x [10.sup.17] 5,2 x [10.sup.14] i 6 [[??].sub.i] 20 [T.sub.ci]/[[omega].sub.oi] 2,5 [q.sub.i]/[T.sub.ci] 3 -[DELTA]([partial derivative] [C.sub.v,i]/[partial derivative][??]) 0,15 [C.sub.v,i] ([T.sub.c] - 0) 1,6 ([C.sub.s] - [C.sub.n])/ [C.sub.n] 1,8 [n.sub.bpi] x [cm.sup.3] 2,3 x [10.sup.13] [[??].sub.i] = [[omega].sub.i]/[[omega].sup.*], [[omega].sup.*] = 5 meV, [[omega].sub.i] is the energy of an optical phonon; [mathematical expression not reproducible] is the critical temperature of the transition, [q.sub.i] is the latent heat of the transition from condensate to supracondensate state; [mathematical expression not reproducible] is the jump in heat capacity during SC transition, [??] = T/[[omega].sup.*]; [C.sub.v,i] ([T.sub.c] - 0) is the heat capacity in the SC phase at the critical point; [C.sub.s] = [C.sub.v]([T.sub.c] - 0), [C.sub.n] = [C.sub.v]([T.sub.c] + 0). The calculations are carried out for the concentration of TI-bipolarons n = [10.sup.21] [cm.sup.-3] and the effective mass of a band electron [m.sup.*] = [m.sup.0]. The table also lists the values of concentrations of TI-bipolarons [n.sub.bpi] for HTSC [YBa.sub.2][Cu.sub.3][O.sub.7], based on the experimental value of the transition temperature [T.sub.c] = 93K (Section 6).

Printer friendly Cite/link Email Feedback | |

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

Author: | Lakhno, V.D. |

Publication: | Advances in Condensed Matter Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 7441 |

Previous Article: | Experimental Demonstration of Special-Shaped 32-Quadrature Amplitude Modulation Constellations for Visible Light Communications. |

Next Article: | Stabilization of Silver Nanoparticles with a Dithiocarbamate Ligand and Formation of Nanocomposites by Combination with Polythiophene Derivative... |

Topics: |