^{1}

^{*}

^{2}

^{*}

MgB_{2} with T_{c} ≈ 40 K, is a record-breaking compound among the s-p metals and alloys. It appears that this material is a rare example of the two band electronic structures, which are weakly connected with each other. Experimental results clearly reveal that boron sub-lattice conduction band is mainly responsible for superconductivity in this simple compound. Experiments such as tunneling spectroscopy, specific heat measurements, and high resolution spectroscopy show that there are two superconducting gaps. Considering a canonical two band BCS Hamiltonian containing a Fermi Surface of π- and σ-bands and following Green’s function technique and equation of motion method, we have shown that MgB_{2} possess two superconducting gaps. It is also pointed out that the system admits a precursor phase of Cooper pair droplets that undergoes a phase locking transition at a critical temperature below the mean field solution. Study of specific heat and density of states is also presented. The agreement between theory and experimental results for specific heat is quite convincing. The paper is organized in five sections: Introduction, Model Hamiltonian, Physical properties, Numerical calculations, Discussion and conclusions.

The surprising discovery of superconductivity in the novel system MgB_{2} with T_{c} = 39 K by Nagamatsu et al. [_{c} for similar compounds.

The crystal structure of MgB_{2} is very simple. It is composed of layers of boron and magnesium, alternating along the c-axis. Each boron layer has a hexagonal lattice similar to that of graphite. The magnesium atoms are arranged between the boron layers in the centers of the hexagons. This has allowed to perform consistent calculations of its electronic structure. Band structure calculations of MgB_{2} show that there are at least two types of nearly separated bands with two superconducting gaps in the excitation spectrum at the Fermi surface. The first one is a heavy hole band, built up of boron σ orbitals. The second one is the broader band with a smaller effective mass, built up mainly of π boron orbitals [3-7].

It is now well established that MgB_{2} is an anisotropic two-gap superconductor [

The Fermi surface consists of four sheets: two three dimensional sheets form the π bonding and antibonding bands, and two nearly cylindrical sheets form the two-dimensional σ-band [4,8]. There is a large difference in the electron-phonon coupling on different Fermi surface sheets and this leads to multiband description of superconductivity. The average electronphonon coupling strength is found to have small values [9-11]. Ummarino et al. [_{2} is a weak coupling two band phononic system where the Coulomb pseudopotential and the interchannel paring mechanism are key terms to interpret the superconducting state. Garland [_{2} the Coulomb effect cannot be considered to explain the reduction of isotope exponent.

It is quite natural to describe a two-gap superconductor by means of a two-band model with interband coupling [14,15]. For MgB_{2}, an approach of such kind is also directly proposed by the nature of the electron spectrum mentioned. There is a number of two band type approaches for superconductivity in MgB_{2} [_{c} cuprate superconductivity [16,17].

Liu et al. [_{2}. In the present study, we use σ-π interband coupling with a strong σ-interband contribution of electron-phonon and Coulobmic nature. Following Liu et al. [

Using two band models, we study the basic MgB_{2} superconductivity characteristics, specific heat and density of states and compare the theoretical results qualitatively with the available experimental data.

The model Hamiltonian has the form [

where

and

Here p and d are momentum labels in the π- and σ- bands respectively with energies and, μ is the common chemical potential. Each band has its proper pairing interaction and, while the pair interchange between the two bands is assured by term.

We have assumed, and we define the following quantities

Further we define

Now in Equation (1) read as

Final Hamiltonian can be written as

We study the Hamiltonian (6) with the Green’s function technique and equation of motion method.

In order to study the physical properties, we define the following normal and anomalous Green’s functions [18- 28]:

Following equation of motion method, we obtain Green’s functions as follows. In obtaining Green’s functions, we have assumed

and

and

Then

1) Green’s functions for π-band

2) Green’s functions for σ-band

Using the following relation [23-27],

and employing the following identity,

we obtain the correlation functions for the Green’s functions given by Equations (9) and (10) as:

where

and are Fermi functions.

Similarly correlation functions for Green’s functions (11) and (12) for σ holes are obtained.

One can define the two superconducting order parameters related to the correlation functions corresponding to Green’s functions and for π- and σ-bands respectively. In a similar manner electronic specific heat can also be defined related to both π- and σ-bands.

Gap parameter is the superconducting order parameter, which can be determined self consistently from the gap equations

In a matrix form, the order parameter for the superconducting state is given by [

where is the pairing interaction constant and function G’s are defined as

Here and are density of states for π- and σ-bands respectively at the Fermi level.

There are two superconducting gaps corresponding to π- and σ-bands in this interband model. One can write the equations for superconducting gaps corresponding to π- and σ-bands as follows

where and is pairing interaction for π- and σ- bands respectively, while the pair interchange between the two bands is assured by the term. The quantity has been supposed to be operative and constant in the energy interval for higher band and lower band, keeping in mind the integration ranges, the gap parameter satisfy the system if the interband interactions are missing, i.e., the transition is solely induced by the interband interaction [

Using Equations (24), we can write the simultaneous equation as

The electronic specific heat per atom of a superconductor is determined from the following relation [3,23-28]1) For π-band

where is the energy of π-band and μ is the common chemical potential.

Substituting from (10) and changing the summation over p into an integration by using the relation, we obtain

where and are given by Equation (16).

2) For σ-band Similarly one can write the expression for electronic specific heat for σ-band, as

where and are

Electronics specific heat for π-band and σ-band are given by Equations (28) and (29) respectively.

The density of states is an important function. This helps in the interpretation of several experimental data e.g. many processes that could occur in crystal but are forbidden because they do not conserve energy. Some of them nevertheless take place, if it is possible to correct the energy imbalance by phonon-assisted processes, which will be proportional to [

where, is the density of state function for π- band. For σ-band we have

where is one particle Green function for π- and σ-bands, defined by Equations (9) and (11) respectively. We have the Green’s function Equation (9),

where. Now solving Equation (32) and using partial fraction method, we obtain

Now substituting the Green function from Equation (33) in Equation (30) and using the delta function property,

we obtain

Changing the summation into integration and after simplification, one obtains

Similarly, for σ-band

Values of various parameters appearing in equations obtained in the previous section are given in _{2}.

For the study of superconducting order parameter for MgB_{2} system within two band models, one finds the following situations1) The SC order parameter for π- and σ-bands. Using Equation (25), one can write

and

Changing the variables as, , and taking, we obtain Solving Equation (37) numerically, one can study the variation of superconductivity order parameters and with temperature corresponding π- and σ-bands. The behavior of superconducting order parameters corresponding to π- and σ-bands with temperature is shown in

2) SC order parameter in the presence of both π- and σ-bands.

The superconducting order parameter for combined π- and σ-bands can be studied by taking a simple sum of both the parameters. Taking the sum of order parameters, one can obtain the values by solving numerically. A comparison of with BCS type curve is shown in

1) For π-band Using Equation (28) and putting a_{1}, a_{2} and, after simplification, we obtain

Changing the variables as, , and using parameters from

2) For σ-band Similarly, we can write expression for specific heat for σ-band using Equation (29)

_{2} system.

The variation electronic specific heat (C_{es}) with temperature (T) for π- and σ-band is shown in

Density of states function for the π-band is given by Equation (35). Now using the following values of y, x_{1} for and x_{2} for, and taking, one obtains,

Similarly, using Equation (36), density of states for σ- band is obtained as

The above two expressions of density of states function for π- and σ-bands are similar, hence we have evaluated the values with different values of x_{1} and x_{2} for π- and σ-bands. The behavior of density of states function for both π- and σ-bands is shown in