A Detailed Study of the Role of Fermi Energy in Determining Properties of Superconducting NbN

Affiliation(s)

^{1}
School of Environmental Sciences, Jawaharlal Nehru University, New Delhi, India.

^{2}
Present Address: B 208 Sushant Lok 1, Gurgaon, Haryana, India.

ABSTRACT

The recent concern with the role of Fermi energy (*E*_{F}) as a determinant of the properties of a superconductor (SC) led us to present new *E*_{F}-dependent equations for the effective mass (m*) of superconducting electrons, their critical velocity, number density, and critical current density, and also the results of the calculations of these parameters for six SCs the *T*_{c}*s* of which vary between 3.72 and 110 K. While this work was based on, besides an idea due to Pines, equations for Tc and the gap at *T* = 0 that are explicitly *E*_{F}-dependent, it employed an equation for the dimensionless construct that depends on *E*_{F} only implicitly; *k* in this equation is the Boltzmann constant, θ is the Debye temperature, and P0 is the critical momentum of Cooper pairs. To meet the demand of consistency, we give here derivation of an equation for y that is also explicitly* E*_{F}-dependent. The resulting framework is employed to (a) review the previous results for the six SCs noted above and (b) carry out a study of NbN which is the simplest composite SC that can shed further light on our approach. The study of NbN is woven around the primary data of Semenov et al. For the additional required inputs, we appeal to the empirical data of Roedhammer et al. and of Antonova et al.

The recent concern with the role of Fermi energy (

1. Introduction

Some of the recent studies [1] - [7] concerned with high-T_{c} superconductors (SCs) have been motivated by the belief that Fermi energy (E_{F}) plays an important role in determining their T_{c}s and gap-structures. These studies make it natural to ask: why not incorporate E_{F} (equivalently, chemical potential μ) into the equations for the T_{c} and the gap of an SC, and then treat it as an independent variable? This is a departure from the usual practice because these parameters are conventionally calculated via equations sans E_{F} because of the assumption

(1)

where k is the Boltzmann constant and is the Debye temperature.

The proposed approach requires, besides the values of T_{c} and, another property of the SC in order to determine E_{F}. Upon choosing critical current density of the SC, new equations for both elemental and composite SCs valid at T = 0 were recently presented in [8] for j_{0} and the following of their properties: m*, , and n_{s}, which denote, respectively, the effective mass of superconducting electrons, their critical velocity at which vanishes, and the density of superconducting electrons. While the results of such a study for Sn, Pb, MgB_{2}, YBCO, Bi-2212, and Tl-2212 were also reported in [8] , it was based on, unlike the equations for and T_{c}, an equation for the dimensionless construct, defined below, that is dependent on E_{F} only implicitly.

(2)

where m*, P_{0}, and E_{F} are in units of electron volts.

To meet the demand of consistency, we present here the derivation of a new equation for that also contains E_{F} explicitly―to put it on par with the equations for T_{c} and. While this leads us to review our earlier results, we also undertake here a detailed study of the superconducting properties of NbN because:

(i) It is the simplest composite SC different samples of which (a) have been fabricated by the same method of preparation, (b) are geometrically similar, but (c) differ in size (e.g., film thickness), and for which (d) data in the form are available, where n_{e} is the density of conduction electrons. This is unlike the composite SCs dealt with earlier, which were not necessarily fabricated by the same method of preparation and for which the values of and n_{e} were not available. We were then led to estimate the values of for these SCs from the data at T = 4.2 K. Given the values of T_{c} and n_{e} for NbN, we can now also shed light on the ratio as a function of T_{c}.

(ii) Since the value of the highest T_{c} reported for it in [9] is 15.25 K, it is the simplest composite SC for which we believe one-phonon exchange mechanism (OPEM) to be operative. This is unlike, e.g., MgB_{2} for which, given its T_{c}, we need to invoke the two- phonon exchange mechanism (TPEM).

(iii) The above features make NbN the simplest testing ground for some key steps of our approach, such as the procedure followed for resolving θ_{NbN} into θ_{Nb} and θ_{N}.

The paper is organized as follows. In Section 2 are reproduced from [8] those equations that constitute our framework in the OPEM scenario, which may be defined as one in which the T_{c} of an SC can be accounted for by a value of the interaction parameter that satisfies the Bogoliubov constraint, i.e., λ < 0.5. Section 3 is devoted to derivation of the new equation for. The study of NbN is taken up in Section 4. A review of our earlier results is taken up in Section 5. The final two sections are devoted to a discussion and conclusions, respectively.

2. E_{F}-Incorporated Equations for Various Properties of an SC

Recalled below from [8] are some of the equations that we need for NbN. In these equations is to be identified with. Further, the equations have been written by assuming that, E_{F}, and have the same values at T = 0 and T = T_{c}, which is in accord with a tenet of the BCS theory. In the following we use and E_{F} interchangeably because they will be seen to differ negligibly. The modified equation for will be derived in the next section.

Equation for:

(3)

where

(4)

and

(5)

Equation for T_{c}:

(6)

where

(7)

and

(8)

In the above equations

(9)

After has been determined via (7) with the input of θ, T_{c}, and any assumed value of, the corresponding value of E_{F} can be determined by the following equation

(10)

Equation for y:

(11)

This equation has been obtained by assuming that

(12)

where

(13)

Equation for j_{0}(E_{F}):

(14)

where

(15)

3. The Modified Equation for y in the OPEM Scenario

Equation (11) has been derived in [10] (pp. 115-120) by assuming Inequality (1). In order to do away with this inequality, we begin here with the following equation for moving CPs because the present derivation differs from the earlier one only beyond it.

(16)

In this equation

(17)

(18)

(19)

(20)

(21)

Equation (16) was obtained via a Bethe-Salpeter equation. It seems interesting to point out that when, it reduces to the well known criterion of superconductivity derived by Thouless via the t-matrix approach, as can be seen from [11] and, in greater detail, in [12] .

The equation for the critical momentum at any temperature follows from (16) by putting W = 0. In terms of we then have

(22)

where

(23)

(24)

(25)

and we have used (9), (13) and (19). Besides, justification to follow, we have dropped E_{3} everywhere except in the denominator of (25) in order to avoid the singularity at. Compared with the earlier equation for, the new feature of (22) is that it has the additional factor of in each of its constituents.

In order to obtain the version of (22), we split both and into two parts: into and for which the limits of integration are and respectively, and into and, where the former is integrated from and the latter from It is then seen that, when T = 0, for and and (+1) for the remaining parts.

Because the constituents of both and differ from one another only in the matter of limits and an overall sign, we now consider the following indefinite integral:

(26)

where we have used (25), put and whence

(27)

Therefore, for (as will be seen to be so), we obtain

(28)

Taking into account the overall sign of, (28) yields

(29)

where Since we replace in the above equation by in order to make contact with (11)., and can be similarly calculated. For the sake of compactness, we define

(30)

(31)

Then substituting and into (22), we obtain

(32)

where

and has been defined for later convenience. Obtained by retaining the factor in and, (32) for y is the equation we had set out to

obtain. It generalizes (11) which was obtained without this factor. While we could earlier solve (11) in the OPEM scenario with the input of alone, solution of (32) requires the additional input of θ and E_{F}. In order to carry out a quick consistency check of (32), we recall that upon solving (6) for Sn (θ = 195 K, T_{c} = 3.72 K,), we had earlier obtained λ = 0.2466. The solution of (11) then led to. This is precisely the value we now obtain by solving (32) with the same inputs for, , and

4. Study of NbN Based on E_{F}-Incorporated Equations

4.1. Outline of Procedure

Working in the OPEM scenario, we

(A) Solve. (6) with the input of θ and T_{c} to determine for different assumed values of.

(B) Solve (32) to obtain the values of corresponding to each pair of values obtained above.

(C) Calculate via (14) for each triplet of values till it is found to agree with its experimental value.

As predictions, this process also yields the values of m*, n_{s}, and via equations derived in [8] and noted in Table 3. As a further check, we calculate via (3) by employing the values of and that led in (C) to the experimental value of.

Before we can proceed as above, we need to fix the Debye temperature of the ions that cause pairing in NbN, i.e., θ_{Nb}.

4.2. Debye Temperature of Nb Ions in NbN

θ_{NbN} is not quoted in [9] . The reported values for it vary in the range 250 - 335 K [13] [14] [15] [16] . We begin by adopting [13]

(33)

We now need to resolve θ_{NbN} into θ_{Nb} and θ_{N}, which must be different because masses of Nb and N ions are different. As in [8] , we do so via the following equations

(34)

(35)

where () is the atomic mass of N (Nb). While the first of the above equations has been routinely used for binaries, the second equation has been derived [10] by assuming that the constituents of the binary simulate weakly coupled oscillations of a double pendulum. The equations above have been written by assuming that Nb is the upper bob of the double pendulum. With, , and θ_{NbN} as in (33), the solutions of these equations yield

(36)

(37)

the corresponding values for θ_{N} being 272.2 and 564.3 K (which we do not need). In the following we shall perform all calculations with both the above values of θ_{Nb}.

4.3. Choosing the Values of T_{c} for Which the Data in [9] Are Addressed

In [9] , while values of T_{c} varying between 9.87 and 15.25 K have been reported for 13 samples of NbN for which the values of lie in range 2.92 - 13.30 MA∙cm^{−2}, the values of have been reported at only three values of, which are 10.72, 14.02, and 15.17 K. Hence we limit the scope of this paper to these values of only.

4.4. A Consistency Check of (6)

If we solve the usual BCS equation for T_{c} (i.e., the equation sans E_{F}) with θ = 105.7 (397.8 K) and T_{c} = 10.72 K, we obtain λ = 0.4142 (0.2682). These are precisely the values we obtain via (6) for the same values of T_{c} and θ and the additional input of μ (or E_{F}) = 100 kθ for each value of θ being considered. Note that manifestly satisfies constraint (1). It is hence seen that (6) incorporating is a valid generalization of the usual equation sans, and may therefore be used for arbitrary values of.

4.5. Fixing Additional Required Inputs

Having fixed the values of θ_{Nb} and T_{c}, we can carry out steps (A) and (B) spelled out in Section 4.1; to carry out step (C) we additionally need the values of γ and the cell parameters of different samples of NbN, which are not given in [9] . We fix these by appealing to the data in [13] . A summary of all the inputs required for this study is given in Table 1. Based on the data in [17] , this table includes the estimated values of at each of the T_{c}s under consideration

4.6. Results

For each of the three values of T_{c} and both the values of θ_{Nb} noted above, we carried out steps (A)-(C) noted in Section (4.1) for. For the sake of brevity, presented in Table 2 are the results corresponding to θ_{Nb} = 105.7 K for only those values of for which the calculated values of are in close agreement with their experimental values noted in Table 1. In obtaining these results we have assumed that θ_{NbN} and hence θ_{Nb} does not change significantly with T_{c}―as is seen from the data in [13] . Thus, up to this stage, having fixed the value of θ_{Nb} as 105.7 K, we have shown that each subset of the experimental values can be accounted for by a corresponding set of values. Since it is pertinent to ask if we could have achieved similar agreement by adopting a different value of θ_{Nb}, we observe that (i) (3) and (6) can be employed only for values of―otherwise we run into complex values because of the factor; (ii) for μ as any multiple of kθ_{Nb}, the value of calculated via either of these equations must be less than 0.5 in order to satisfy the Bogoliubov constraint, and (iii) for any value of increases as μ is increased.

Table 1. Experimental values of [9] , [13] , and [17] employed for the study of NbN in this paper.

Table 2. Results of calculations for θ_{Nb} = 105.7 K. The value of against each T_{c} is the one that led―via the values of E_{F}, , , and (the gram-atomic volume of NbN)―to a value of in close agreement with its experimental value noted in Table 1. was calculated with the input of a_{0} from Table 1 and the atomic masses of the N_{b} and N, as in [8] .

We now take up the results following from θ_{Nb} = 397.8 K. The least permissible value of μ corresponding to it, i.e., led to MA/cm^{2} for T_{c} = 10.72 K and MA/cm^{2} for T_{c} = 14.02 K. Since both these values are greater than their experimental counterparts, in the light of observation (iii) above, one might attempt to employ lower values of―which is ruled out because of (i). In fact the value 105.7 K seems like the upper limit for θ_{Nb} because we had to employ the least value of corresponding to it in order to achieve agreement between the calculated and the experimental values of at T_{c} = 10.72 K. As a concrete example in support of this statement, we note that θ_{Nb} = 125 K, led via the least permissible value of corresponding to it to the following results:

Since this value of exceeds the experimental value, we need to employ a lower value of―which is impermissible because we have already employed for it the lowest allowed value.

Our considerations so far have been based on the derived values of θ_{Nb} from θ_{NbN} = 335 K. In order to find if there is a lower limit on the value of θ_{Nb}, we now report our findings based on the values of θ_{Nb} derived from the lowest value of θ_{NbN} that was noted above, i.e., 250 K. This value leads via (34) and (35) to θ_{Nb} = 296.8 (Nb as the upper bob) and θ_{Nb} = 78.9 K (Nb as the lower bob). Since the former of these values exceeds the upper limit noted above, we did not pursue it any further. For the latter value, we obtained for any assumed value of

Because both these values of are in conflict with the Bogoliubov constraint, we conclude that θ_{Nb} cannot be as low as 78.9 K. The value closest to it that yields values of satisfying the Bogoliubov constraint at both the T_{c}s is θ_{Nb} = 100 K, for which, e.g.,

Above considerations raise the question: Could for NbN? If so, it would put NbN in the category of heavy-fermion SCs [18] . Since there is no compelling reason to believe that this may be so, we did not pursue this idea.

Given in Table 3 are the predicted values of various parameters concomitant with

Table 3. With θ_{Nb} = 105.7 K, predicted values of various parameters of NbN that are concomitant with the calculated values of given against each T_{c} in Table 2.

Notes: (i) The equations employed for the calculation of the above parameters have been derived in [8] and are as follows: (ii) The product [n_{s}(E_{F}) e v_{0}] at each T_{c} yields the same value for j_{0} as was calculated via (14) and given in Table 2. (ii) The values of are, and for T_{c} = 10.72, 14.02 and 15.17 K, respectively, which justify the approximation made in obtaining (32).

the experimental values of T_{c} and of NbN at the three T_{c}s. Among these, the values of are in reasonably good agreement with their experimental counterparts (see Table 1), considering that the latter are estimated values based on the data of [17] .

5. A Review of the Results Obtained in [8] in View of the Modified Equation for y

For Sn and Pb, all our earlier results remain unchanged because solution of (32) for these elements yields the same values for that were obtained via (11). Since the values of that were needed for these elements are rather large, 55 for Sn and 108 for Pb, this result was to be expected; it also establishes that (32) is a valid generalization of (11). To bring out the extent to which the solutions of the two equations differ for low values of, we note that if we erroneously employ (11) for Sn for, θ = 195 K and λ = 0.2516 (these values are consistent with of the SC), then we obtain y = 20.083 [8] ; employment of (11) in this case is erroneous because the equation was obtained by assuming that. On the other hand, solution of (32) for this case leads to.

For each of the high-T_{c} SCs dealt with in [8] , there are two―say, and― and two in the problem. The E_{F}-dependent equation for that we now need to employ is

(38)

where was defined in (32), θ being the Debye temperature of the SC, and E_{2} was defined in (13). It is hence seen that―as it ought to be. Without the multipliers and, would denote in the first term on the RHS of (38) and in second term, whereas with the multipliers has the same definition (i.e.,) for both the terms.

Equation (38) generalizes (32) to the TPEM scenario; because it explicitly contains E_{F} as a variable, it is also a generalized version of equation (30) in [8] . Upon solving (38) with the input of, and E_{F} as given in [8] for any of the high-T_{c} SCs, we obtain the same value for that we had obtained earlier. Notwithstanding the fact that all our results reported in [8] remain unchanged is fortuitous―for lower values of than were required in [8] , we ought to employ the more accurate (38) rather than equation (30) in [8] .

6. Discussion

In connection with fixing θ_{Nb}, we recall that Debye temperature is just another way to specify Debye frequency; it is not to be confused with thermodynamic temperature. We now note that, based on neutron powder diffraction experiments, different values of Debye temperature for the constituents of anisotropic LCO have been reported [19] . This lends support to the idea that the Debye temperature of a composite SC needs to be “resolved.” The results reported here depend only on the value of θ_{Nb}, for the identification of which we have simply employed (34) and (35) as a vehicle.

Among the five variables that determine―see Equation (14) ―seems to stand alone. We draw attention to a discussion of this variable in [8] .

7. Conclusions

The main results of this paper are: (i) a new E_{F}-dependent equation for the dimensionless construct defined in (2) has been derived, (ii) it has been shown that the experimental values of T_{c}, , and of NbN are explicable in the OPEM scenario by a value of θ_{Nb} in the range 100 - 106 K, (iii) predictions have been made about the values of m*, n_{s}, and that are concomitant with the T_{c} and values of NbN, (iv) the greater the value of the ratio, the greater is the value of T_{c}, and (v) it has been pointed out that we need to employ the new equations for presented here when.

The work reported here is continuation of an attempt to find via theory tangible clues about raising the T_{c}s of composite SCs. The role of experiment in this quest can hardly be over-emphasized. While huge amounts of such data about hundreds of SCs are now available, we have not come across a single composite SC for which all the relevant parameters identified here, i.e., θ, T_{c}, , , m*, , n_{e}, n_{s}, , and, have been reported.

We conclude by noting that the derivations of most of the equations employed in this paper and the concepts on which they are based, e.g., multiple Debye temperatures, superpropagator, and the Bogoliubov constraint, can be found at one place in [10] .

Acknowledgements

The author thanks Dr. A. Semenov for kindly responding to his queries concerned with the experimental data reported in [9] , and Professor D.C. Mattis for encouragement.

Cite this paper

Malik, G. (2017) A Detailed Study of the Role of Fermi Energy in Determining Properties of Superconducting NbN.*Journal of Modern Physics*, **8**, 99-109. doi: 10.4236/jmp.2017.81009.

Malik, G. (2017) A Detailed Study of the Role of Fermi Energy in Determining Properties of Superconducting NbN.

References

[1] Lee, D. (2012) Nature Physics, 8, 364-365.

[2] Zhang, Y., et al. (2012) Nature Physics, 8, 371-375.

[3] Allan, M.P., et al. (2012) Science, 336, 563-567.

https://doi.org/10.1126/science.1218726

[4] Lin, X, Zhu, Z., Fauqué, B. and Behnia, K. (2013) Physical Review X, 3, 021002.

https://doi.org/10.1103/PhysRevX.3.021002

[5] Alexandrov, A.S. (2001) Physica C: Superconductivity, 363, 231-236.

https://doi.org/10.1016/S0921-4534(01)01095-4

[6] Jarlborg, T. and Bianconi, A. (2013) Physical Review B, 87, 054514.

https://doi.org/10.1103/PhysRevB.87.054514

[7] van der Marel, D., van Mechelen, J.L.M. and Mazin, I.I. (2011) Physical Review B, 84, 205111.

https://doi.org/10.1103/PhysRevB.84.205111

[8] Malik, G.P. (2016) Journal of Superconductivity and Novel Magnetism, 29, 2755-2764.

[9] Semenov, A., et al. (2009) Physical Review B, 80, 054510.

https://doi.org/10.1103/PhysRevB.80.054510

[10] Malik, G.P. (2016) Superconductivity: A New Approach Based on the Bethe-Salpeter Equation in the Mean Field Approach. World Scientific, Singapore.

https://doi.org/10.1142/9868

[11] Thouless, D.J. (1960) Perturbation Theory in Statistical Mechanics and the Theory of Superconductivity. Annals of Physics, 10, 553-588.

https://doi.org/10.1016/0003-4916(60)90122-6

[12] Nozières, P. and Schmitt-Rink, S. (1985) Bose Condensation in an Attractive Fermion Gas: From Weak to Strong Coupling Superconductivity. Journal of Low Temperature Physics, 59, 195-211.

https://doi.org/10.1007/BF00683774

[13] Roedhammer, P., Gmelin, E., Weber, W. and Remeika, J.P. (1977) Observation of Phonon Anomalies in NbCxN1-x Alloys. Physical Review B, 15, 711-717.

https://doi.org/10.1103/PhysRevB.15.711

[14] Geballe, T.H., et al. (1966) High-Temperature SP-Band Superconductors. Physics, 2, 293-310.

[15] Pellan, Y., Dousselin G., Pinel, J. and Sohn, Y.U. (1990) Temperature and Magnetic Field Dependence of NbN Film Resistivity: 3D Weak Localization Effects. Journal of Low Temperature Physics, 78, 63-77.

https://doi.org/10.1007/BF00682110

[16] Geibel, C., Rietschel, H., Junod, A., Pelizzone M. and Muller, J. (1985) Electronic Properties, Phonon Densities of States and Superconductivity in Nb1-xVxN. Journal of Physics F: Metal Physics, 15, 405-416.

https://doi.org/10.1088/0305-4608/15/2/014

[17] Antonova, E.A., Dzhuraev, D.R., Motulevich, G.P. and Sukhov, V.A. (1981) Superconducting Energy Gap of Niobium Nitride. Journal of Experimental and Theoretical Physics, 53, 1270-1271.

[18] Malik, G.P. (2015) A Study of Heavy-Fermion Superconductors via BCS Equations Incorporating Chemical Potential. Journal of Modern Physics, 6, 1233-1241.

https://doi.org/10.4236/jmp.2015.69128

[19] Kwei, G.H., Lawson, A.C. and Mostoller, M. (1991) Vibrational Properties and Atomic Debye Temperature for La2CuO4 from Neutron Powder Diffraction. Physica C: Superconductivity, 175, 135-142.

https://doi.org/10.1016/0921-4534(91)90244-S

[1] Lee, D. (2012) Nature Physics, 8, 364-365.

[2] Zhang, Y., et al. (2012) Nature Physics, 8, 371-375.

[3] Allan, M.P., et al. (2012) Science, 336, 563-567.

https://doi.org/10.1126/science.1218726

[4] Lin, X, Zhu, Z., Fauqué, B. and Behnia, K. (2013) Physical Review X, 3, 021002.

https://doi.org/10.1103/PhysRevX.3.021002

[5] Alexandrov, A.S. (2001) Physica C: Superconductivity, 363, 231-236.

https://doi.org/10.1016/S0921-4534(01)01095-4

[6] Jarlborg, T. and Bianconi, A. (2013) Physical Review B, 87, 054514.

https://doi.org/10.1103/PhysRevB.87.054514

[7] van der Marel, D., van Mechelen, J.L.M. and Mazin, I.I. (2011) Physical Review B, 84, 205111.

https://doi.org/10.1103/PhysRevB.84.205111

[8] Malik, G.P. (2016) Journal of Superconductivity and Novel Magnetism, 29, 2755-2764.

[9] Semenov, A., et al. (2009) Physical Review B, 80, 054510.

https://doi.org/10.1103/PhysRevB.80.054510

[10] Malik, G.P. (2016) Superconductivity: A New Approach Based on the Bethe-Salpeter Equation in the Mean Field Approach. World Scientific, Singapore.

https://doi.org/10.1142/9868

[11] Thouless, D.J. (1960) Perturbation Theory in Statistical Mechanics and the Theory of Superconductivity. Annals of Physics, 10, 553-588.

https://doi.org/10.1016/0003-4916(60)90122-6

[12] Nozières, P. and Schmitt-Rink, S. (1985) Bose Condensation in an Attractive Fermion Gas: From Weak to Strong Coupling Superconductivity. Journal of Low Temperature Physics, 59, 195-211.

https://doi.org/10.1007/BF00683774

[13] Roedhammer, P., Gmelin, E., Weber, W. and Remeika, J.P. (1977) Observation of Phonon Anomalies in NbCxN1-x Alloys. Physical Review B, 15, 711-717.

https://doi.org/10.1103/PhysRevB.15.711

[14] Geballe, T.H., et al. (1966) High-Temperature SP-Band Superconductors. Physics, 2, 293-310.

[15] Pellan, Y., Dousselin G., Pinel, J. and Sohn, Y.U. (1990) Temperature and Magnetic Field Dependence of NbN Film Resistivity: 3D Weak Localization Effects. Journal of Low Temperature Physics, 78, 63-77.

https://doi.org/10.1007/BF00682110

[16] Geibel, C., Rietschel, H., Junod, A., Pelizzone M. and Muller, J. (1985) Electronic Properties, Phonon Densities of States and Superconductivity in Nb1-xVxN. Journal of Physics F: Metal Physics, 15, 405-416.

https://doi.org/10.1088/0305-4608/15/2/014

[17] Antonova, E.A., Dzhuraev, D.R., Motulevich, G.P. and Sukhov, V.A. (1981) Superconducting Energy Gap of Niobium Nitride. Journal of Experimental and Theoretical Physics, 53, 1270-1271.

[18] Malik, G.P. (2015) A Study of Heavy-Fermion Superconductors via BCS Equations Incorporating Chemical Potential. Journal of Modern Physics, 6, 1233-1241.

https://doi.org/10.4236/jmp.2015.69128

[19] Kwei, G.H., Lawson, A.C. and Mostoller, M. (1991) Vibrational Properties and Atomic Debye Temperature for La2CuO4 from Neutron Powder Diffraction. Physica C: Superconductivity, 175, 135-142.

https://doi.org/10.1016/0921-4534(91)90244-S