A Fermi Energy-Incorporated Framework for Dealing with the Temperature- and Magnetic Field-Dependent Critical Current Densities of Superconductors and Its Application to Bi-2212

Show more

1. Introduction

The critical current density (j_{c}) of a superconductor (SC) is an important parameter because the greater its value, the greater is the practical use to which it can be put. Most of the plethora of formulae available in the literature for calculating this parameter may be mainly categorized as following from the framework of the Londons’ equations or the Ginsburg-Landau (GL) equations. The salient features of such approaches are that they are based on diverse criteria such as the type of SC being dealt with (type I or II) and its geometry [1] [2] [3] [4] [5]. Another limitation of both the Londons’ and the GL theories is that they work for the j_{c} of an SC only when its temperature T is close to its critical temperature T_{c}. In contrast with these phenomenological approaches, there is also available in the literature a smaller body of work dealing with the j_{c} of an SC on the basis of the microscopic theory of superconductivity. Notable among these is the often-used approach of Kupriyanov and Lukichev [6] based on a simplification of the Eilenberger theory, which in turn is derived from the original Gor’kov theory under the assumption that
${\rho}_{F}l\gg 1$, where
${\rho}_{F}$ is the electrical resistivity of electrons at the Fermi surface and l their mean free path.

An overview of the study reported herein is as follows. Guided by a substantial body of recent work, e.g., [6] - [12], which suggests that low values of the Fermi energy (E_{F}) play a pivotal role in determining properties of high-T_{c} SCs, we have been following a course where E_{F} is directly incorporated into the generalized BCS equations for the gaps (Δs) and the T_{c}s of both elemental and composite SCs [13] [14]. This is a departure from the usual practice since these parameters are conventionally calculated via equations independent of E_{F} because of the assumption that
${E}_{F}/k\theta \gg 1$, where k is the Boltzmann constant and θ the Debye temperature of the SC. As a supplement to this framework, we reported in [15] and [16] the results of an exercise that includes an E_{F}-incorporated equation for j_{c}, leading to a unified framework for dealing with the Δ_{s}, T_{c}s and j_{c}s of both elemental and composite SCs at T = 0 and H = 0, where H is the applied external magnetic field. Since the j_{c}-values of an SC are generally reported for non-zero values of both T and H, we present here a framework to deal with such a situation. In essence, the present work is concerned with a generalization of: 1) the work reported in [17] to include E_{F} in the pairing equations that already incorporate T and H and 2) the work reported in [18] to include H in the equations that already incorporate T, E_{F} and the momentum P of the pairs. We believe that presented herein is the first attempt that brings the T-and H-dependent j_{c} of an SC under the purview of an E_{F}-incorporated microscopic theory of superconductivity.

The paper is organized as follows. In Section 2, we show how, without making the usual approximation
${E}_{F}\gg k\theta $, the dynamical equation for a Cooper pair interacting via the model BCS interaction can be generalized to include T and H via the Matsubara prescription and the Landau quantization scheme, respectively. We thus obtain the pairing equation incorporating E_{F}, T and H corresponding to the 1-phonon exchange mechanism (1PEM) and P = 0, which is appropriate for dealing with the situation when j = 0. In order to deal with the situation when j ≠ 0, we need to do away with the P = 0 constraint of this section. This is done in Section 3, assuming that the dimensionless interaction parameter remains unchanged as we move out of the center-of-mass frame to the lab frame. Based on the equations thus obtained, in Section 4 are derived expressions for the density of the superconducting electrons n_{s}, their critical velocity v_{c}, j_{c} and s ≡ m^{*}/m_{e}, where m^{*} is the effective mass of an electron and m_{e} the free electron mass. Application of these equations to deal with the empirical data of Bi_{2}Sr_{2}CaCu_{2}O_{8} (Bi-2212) for non-zero values of T and H is taken up in Section 5.1 (for j = 0) and Section 5.2 (for j ≠ 0). In Section 5.3 we present the solutions of these equations in the scenario of 1PEM due to each of the ion species: Ca, Sr, and Bi. Sections 6 and 7 are devoted, respectively, to a Discussion of our findings and Conclusions. Finally, for the reader’s convenience, we have given in an Appendix the conversion factors needed to go over from the natural system of units to the more familiar units employed in the BCS theory.

2. Pairing Equation Incorporating Fermi Energy, Temperature and Applied Field When the Momentum of Cooper Pairs Is Zero

Incorporating T and H into the pairing equation is a two-step process [17]: 1) temperature-generalization of the T = 0 Bethe-Salpeter Equation (BSE) via the Matsubara recipe and 2) further generalization of the equation obtained in (1) to include H via the Landau quantization scheme.

For elemental SCs for which pairing arises from the one-phonon exchange mechanism (1PEM) due to a single species of ions, Step (1) leads to

$1=\frac{V}{{\left(2\pi \right)}^{3}}\frac{1}{2}{\displaystyle {\int}_{{E}_{F}-k\theta}^{{E}_{F}+k\theta}{\text{d}}^{3}p\frac{\mathrm{tanh}\left[\left(\beta /2\right)\left({p}^{2}/2m-{E}_{F}-W/2\right)\right]}{{p}^{2}/2m-{E}_{F}-W/2}},$ (1)

where (−V) ≠ 0 in a narrow region of $\pm k\theta $ around the Fermi surface is the model BCS interaction parameter, $\beta =1/kT$, k is the Boltzmann constant, θ the Debye temperature of the ions, and W is one-half the binding energy of a pair which is to be identified with Δ. The units employed are: eV, $\hslash =c=1$, generally. However, for the convenience of the reader, in all the final equations that are actually employed in our calculations, the factors of $\hslash $ and c have been made explicit.

If the field H is applied in the z-direction, then Step (2) consists of making the following substitutions in (1):

$\begin{array}{l}{\displaystyle \int \text{d}{p}_{x}\text{d}{p}_{y}}=2\pi eH{\displaystyle \underset{n}{\sum}\text{\hspace{0.05em}}},\text{}\frac{{p}_{x}^{2}}{2m}+\frac{{p}_{y}^{2}}{2m}=\left(n+1/2\right)\hslash {\Omega}_{1}\left(H\right)\\ {\Omega}_{1}\left(H\right)=\frac{eH}{{m}^{*}c}={\Omega}_{0}H/s,\\ {\Omega}_{0}=e/{m}_{e}c=1.7588\times {10}^{7}\text{\hspace{0.17em}}\text{rad}\cdot {\mathrm{sec}}^{-1}\cdot {\text{G}}^{-1},\text{}\left(s={m}^{*}/{m}_{e}\right)\end{array}$ (2)

where Ω_{0}(H) is the cyclotron frequency corresponding to the free electron mass and we assume that m^{*} = m_{e} when j_{c} = 0. The transverse components of momentum are thus quantized into Landau levels and we have (1) as

$1=\frac{eHV}{8{\pi}^{2}}{\displaystyle {\int}_{-L}^{L}\text{d}{p}_{z}}{\displaystyle \underset{n=0}{\overset{{n}_{m}}{\sum}}\frac{\mathrm{tanh}\left[\left(\beta /2\right)\left({p}_{z}^{2}/2{m}^{*}-{E}_{F}+\left(n+1/2\right)\hslash {\Omega}_{1}\left(H\right)-W/2\right)\right]}{{p}_{z}^{2}/2{m}^{*}-{E}_{F}+\left(n+1/2\right)\hslash {\Omega}_{1}\left(H\right)-W/2}}$ (3)

where L and n_{m} are usually
$\infty $. Since the energy of an electron in our problem is constrained to lie in a narrow shell around the Fermi surface, we fix L and n_{m} by appealing to the law of equipartition of energy and split the region in which V ≠ 0 as

$\begin{array}{l}-\frac{2k\theta}{3}\le \frac{{p}_{x}^{2}+{p}_{y}^{2}}{2{m}^{*}}=\left(n+1/2\right)\hslash {\Omega}_{1}\left(H\right)\le \frac{2k\theta}{3}\\ -\frac{k\theta}{3}\le \frac{{p}_{z}^{2}}{2{m}^{*}}-{E}_{F}\le \frac{k\theta}{3}.\end{array}$ (4)

In terms of the variable $\xi ={p}_{z}^{2}/2{m}^{*}-{E}_{F}$, $\text{d}{p}_{z}=\left(1/2\right)\sqrt{2{m}^{*}/\left({E}_{F}+\xi \right)}\text{d}\xi $ (3) is given by

$\begin{array}{c}1=\frac{eHV}{16{\pi}^{2}}\sqrt{\frac{2{m}^{*}}{{E}_{F}}}{\displaystyle {\int}_{-{L}_{1}\left(\theta \right)}^{{L}_{1}\left(\theta \right)}\frac{\text{d}\xi}{\sqrt{1+\xi /{E}_{F}}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times {\displaystyle \underset{n=0}{\overset{{n}_{m}\left(\theta ,H\right)}{\sum}}\frac{\mathrm{tanh}\left[\left(\beta /2\right)\left(\xi +\left(n+1/2\right)\hslash {\Omega}_{1}\left(H\right)-W/2\right)\right]}{\xi +\left(n+1/2\right)\hslash {\Omega}_{1}\left(H\right)-W/2}}\end{array}$ (5)

where

${L}_{1}\left(\theta \right)=k\theta /3$

and

${n}_{m}\left(\theta ,H\right)=\text{floor}\left\{\frac{2k\theta}{3\hslash {\Omega}_{1}\left(H\right)}-\frac{1}{2}\right\}.$ (6)

Putting W = 0 and rewriting (5) in terms of T and the dimensionless variable
$x=\xi /\hslash {\Omega}_{1}\left({H}_{c}\right)$, we obtain the equation for the critical field H_{c} at temperature T (equivalently, the equation for the critical temperature T_{c} at field H) as

$\begin{array}{c}1={\lambda}_{m}{\displaystyle {\int}_{-{L}_{2}\left(\theta ,{H}_{c}\right)}^{{L}_{2}\left(\theta ,{H}_{c}\right)}\frac{\text{d}x}{\sqrt{1+\hslash {\Omega}_{1}\left({H}_{c}\right)x/{E}_{F1}}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times {\displaystyle \underset{n=0}{\overset{{n}_{m2}\left(\theta ,{H}_{c}\right)}{\sum}}\frac{\mathrm{tanh}\left[\left(\hslash {\Omega}_{1}\left({H}_{c}\right)/2k{T}_{c}\right)\left(x+n+1/2\right)\right]}{x+n+1/2}}\end{array}$ (7)

where E_{F} has been relabeled as E_{F}_{1} in order to distinguish it from the E_{F} that occurs in the equations when P ≠ 0,

${L}_{2}\left(\theta ,{H}_{c}\right)=k\theta /\left[3\hslash {\Omega}_{1}\left({H}_{c}\right)\right],$

${n}_{m2}\left(\theta ,{H}_{c}\right)=\text{floor}\left[2{L}_{2}\left(\theta ,{H}_{c}\right)-1/2\right]$

and

${\lambda}_{m}\left({H}_{c},V\right)=\frac{e{H}_{c}V}{16{\pi}^{2}}\sqrt{\frac{2{m}^{*}}{{E}_{F1}}}.$ (8)

Incorporating T, H and E_{F}, (7) is the equation we had set out to obtain; it reduces to Equation (16) of [17] when
${E}_{F}\gg k\theta $. The integrand in (7) is manifestly dimensionless. In Appendix A, we provide the necessary conversion factors which show that λ_{m} too is dimensionless and enable one to employ in our framework the more familiar BCS units, i.e., eV-cm^{3} for V and Gauss for H.

3. Pairing Equation Incorporating Fermi Energy, Temperature and Applied Field When the Momentum of Cooper Pairs Is Non-Zero

It is a tenet of the BCS theory that the same interaction parameter λ occurs in both the equation for {T = T_{c}, Δ = 0} and the equation for {T = 0, Δ_{0} ≠ 0}, where Δ_{0} is the gap at T = 0. Similarly, we assume here that λ_{m} remains unchanged when we go over from the P = 0 to P ≠ 0 equations. We now draw attention to the fact that in [15] where we dealt with the Δ_{0}s and j_{c}s of various SCs, it was assumed that E_{F} too has the same value in both the equation for Δ_{0} (H = 0) and the equation for P_{c}/j_{c}--a plausible justification for which being that we were dealing with the j_{c}s too at T = 0 and H = 0. For the sake of generality, we now assume that the value of E_{F} when P ≠ 0 is different from its value when P = 0. Hence, E_{F} is labelled as E_{F}_{2} in the present section.

The T- and P-dependent equation for pairing in the 1PEM scenario is [18]:

$1=\frac{V}{16{\pi}^{2}}{\displaystyle \underset{\mathcal{l}}{\overset{u}{\int}}{\text{d}}^{3}p\frac{\mathrm{tanh}\left[\beta C\left(p\right)/2\right]+\mathrm{tanh}\left[\beta D\left(p\right)/2\right]}{C\left(p\right)+D\left(p\right)}},$ (9)

where

$\begin{array}{l}C\left(p\right)=\frac{{\left(P/2+p\right)}^{2}}{2{m}^{*}}-E/2=\frac{{P}^{2}}{8{m}^{*}}+\frac{\left|P\right|\left|p\right|\mathrm{cos}\left(P,p\right)}{2{m}^{*}}+\frac{{p}^{2}}{2{m}^{*}}-{E}_{F2}-W/2\\ D\left(p\right)=\frac{{\left(P/2-p\right)}^{2}}{2{m}^{*}}-E/2=\frac{{P}^{2}}{8{m}^{*}}-\frac{\left|P\right|\left|p\right|\mathrm{cos}\left(P,p\right)}{2{m}^{*}}+\frac{{p}^{2}}{2{m}^{*}}-{E}_{F2}-W/2\end{array}$

V ≠ 0 in the region

${E}_{F2}-k\theta \le \frac{{\left(P/2+p\right)}^{2}}{2{m}^{*}},\frac{{\left(P/2-p\right)}^{2}}{2{m}^{*}}\le {E}_{F2}+k\theta ,$ (10)

and the total energy of a pair is $E=2{E}_{F2}+W$. When P = 0, (10) fixes the limits $\mathcal{l}$ and u in (9) as $\mathcal{l}={E}_{F2}-k\theta $ and $u={E}_{F2}+k\theta $. When P ≠ 0, the lower limit in (9) follows from:

${E}_{F2}-k\theta \le \frac{{\left({P}_{c}/2-p\right)}^{2}}{2{m}^{*}}={p}^{2}/2{m}^{*}-\left|{P}_{c}\right|{p}_{z}/2{m}^{*},$ (11)

where P_{c} is the critical momentum corresponding to W = 0,
${p}_{z}=\left|p\right|\mathrm{cos}\left({P}_{c},p\right)$ and the
${P}_{c}^{2}$ -term in the expansion of
${\left({P}_{c}/2\pm p\right)}^{2}$ has been neglected (as justified in [18]). Note that (11) automatically ensures that
${E}_{F2}-k\theta <{\left({P}_{c}/2+p\right)}^{2}/2{m}^{*}$. Appealing to the law of equipartition of energy as earlier, we now assume that

${p}_{z}^{2}/2{m}^{*}\cong {E}_{F2}/3,\text{i}\text{.e}\text{.,}{p}_{z}=\sqrt{2{m}^{*}{E}_{F2}/3},$

whence we have (11) as

$\begin{array}{l}-{E}_{1}\left(\theta \right)+{E}_{2}\left({E}_{F2}\right)\le \left[\left({p}_{z}^{2}/2{m}^{*}-{E}_{F2}\right)+\left({p}_{x}^{2}+{p}_{y}^{2}\right)/2{m}^{*}\right]\\ \left({E}_{1}\left(\theta \right)=k\theta ,\text{}{E}_{2}\left({E}_{F2}\right)={P}_{c}\sqrt{{E}_{F2}/6{m}^{*}}\right).\end{array}$ (12)

which is now the lower limit in (9). Working out the upper limit similarly, we have the limits of (9) as:

$\begin{array}{c}-{E}_{1}\left(\theta \right)+{E}_{2}\left({E}_{F2}\right)\le \left[\left({p}_{z}^{2}/2{m}^{*}-{E}_{F2}\right)+\left({p}_{x}^{2}+{p}_{y}^{2}\right)/2{m}^{*}\right]\\ \le {E}_{1}\left(\theta \right)-{E}_{2}\left({E}_{F2}\right).\end{array}$ (13)

Employing substitutions (2) in order to introduce H into (9), we obtain

$\begin{array}{c}1=\frac{e{H}_{c}V}{8{\pi}^{2}}{\displaystyle \underset{-{\mathcal{l}}_{1}\left(\theta ,{E}_{F2}\right)}{\overset{{\mathcal{l}}_{1}\left(\theta ,{E}_{F2}\right)}{\int}}\text{d}{p}_{z}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times {\displaystyle \underset{n=0}{\overset{{n}_{u}{}_{1}\left(\theta ,{E}_{F2},{H}_{c}\right)}{\sum}}\frac{\mathrm{tanh}\left[{A}_{1}\left({T}_{c},{H}_{c},{E}_{F2},{p}_{z}\right)\right]+\mathrm{tanh}\left[{B}_{1}\left({T}_{c},{H}_{c},{E}_{F2},{p}_{z}\right)\right]}{2\left[\left({p}_{z}^{2}/2{m}^{*}-{E}_{F2}\right)+\left(n+1/2\right)\hslash {\Omega}_{2}\left({H}_{c}\right)\right]}}\end{array}$ (14)

where

${\Omega}_{2}\left({H}_{c}\right)={\Omega}_{0}{H}_{c}/s$ (15)

$\begin{array}{l}{A}_{1}\left({T}_{c},{H}_{c},{E}_{F},{p}_{z}\right)=\left(\frac{1}{2k{T}_{c}}\right)\left[\left(\frac{{p}_{z}^{2}}{2{m}^{*}}-{E}_{F2}\right)+\left(n+\frac{1}{2}\right)\hslash {\Omega}_{2}\left({H}_{c}\right)+{E}_{2}\left({E}_{F2}\right)\right]\\ {B}_{1}\left({T}_{c},{H}_{c},{E}_{F},{p}_{z}\right)=\left(\frac{1}{2k{T}_{c}}\right)\left[\left(\frac{{p}_{z}^{2}}{2{m}^{*}}-{E}_{F2}\right)+\left(n+\frac{1}{2}\right)\hslash {\Omega}_{2}\left({H}_{c}\right)-{E}_{2}\left({E}_{F2}\right)\right]\end{array}$ (16)

and

$\begin{array}{l}{\mathcal{l}}_{1}\left(\theta ,{E}_{F2}\right)=\left(1/3\right)\left[{E}_{1}\left(\theta \right)-{E}_{2}\left({E}_{F2}\right)\right]\\ {n}_{u}{}_{1}\left(\theta ,{E}_{F2},{H}_{c}\right)=\text{floor}\left[2{\mathcal{l}}_{1}\left(\theta ,{E}_{F2}\right)/\hslash {\Omega}_{2}\left({H}_{c}\right)-1/2\right].\end{array}$

We note that the factors of (1/3) in the limits of the integral and (2/3) in the upper limit of the sum occur because, appealing to the equipartition law, we have split the inequalities (13) as

$\begin{array}{l}\left(1/3\right)\left[-{E}_{1}\left(\theta \right)+{E}_{2}\left({E}_{F2}\right)\right]\le \left({p}_{z}^{2}/2{m}^{*}-{E}_{F2}\right)\le \left(1/3\right)\left[{E}_{1}\left(\theta \right)-{E}_{2}\left({E}_{F2}\right)\right]\\ \left(2/3\right)\left[-{E}_{1}\left(\theta \right)+{E}_{2}\left({E}_{F2}\right)\right]\le \left[\left({p}_{x}^{2}+{p}_{y}^{2}\right)/2{m}^{*}=\left(n+1/2\right)\hslash {\Omega}_{2}\left({H}_{c}\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\le \left(2/3\right)\left[{E}_{1}\left(\theta \right)-{E}_{2}\left({E}_{F2}\right)\right].\end{array}$

With

$\xi ={p}_{z}^{2}/2{m}^{*}-{E}_{F2},\text{d}{p}_{z}=\frac{1}{2}\sqrt{2{m}^{*}/\left({E}_{F2}+\xi \right)}\text{d}\xi $

(14) is transformed into

$\begin{array}{c}1=\frac{{\lambda}_{m}\left(V,{H}_{c}\right)}{2}{\displaystyle \underset{-{\mathcal{l}}_{1}\left(\theta ,{E}_{F2}\right)}{\overset{{\mathcal{l}}_{1}\left(\theta ,{E}_{F2}\right)}{\int}}\frac{\text{d}\xi}{\sqrt{1+\xi /{E}_{F2}}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times {\displaystyle \underset{n=0}{\overset{{n}_{u1}\left(\theta ,{E}_{F2},{H}_{c}\right)}{\sum}}\frac{\mathrm{tanh}\left[{A}_{2}\left(T,{H}_{c},{E}_{F2},\xi \right)\right]+\mathrm{tanh}\left[{B}_{2}\left(T,{H}_{c},{E}_{F2},\xi \right)\right]}{\xi +\left(n+1/2\right)\hslash {\Omega}_{2}\left({H}_{c}\right)}},\end{array}$ (17)

where λ_{m}(V, H_{c}) is given by (8), A_{2}(..) and B_{2}(..) are given by A_{1}(..) and B_{1}(..), respectively, with
$\left({p}_{z}^{2}/2{m}^{*}-{E}_{F2}\right)$ replaced by ξ. Employing the dimensionless variables
$x=\xi /\hslash {\Omega}_{2}\left({H}_{c}\right)$ and

$y={E}_{1}\left(\theta \right)/{E}_{2}\left({E}_{F2}\right)=\frac{k\theta}{{P}_{c}}\sqrt{6{m}^{*}/{E}_{F2}},$ (18)

we finally obtain the desired equation incorporating T, H and E_{F2} when P ≠ 0 in terms of y as

$\begin{array}{c}1=\frac{{\lambda}_{m}\left(V,{H}_{c}\right)}{2}{\displaystyle \underset{-{\mathcal{l}}_{2}\left(\theta ,{H}_{c},r,y\right)}{\overset{{\mathcal{l}}_{2}\left(\theta ,{H}_{c},r,y\right)}{\int}}\frac{\text{d}x}{\sqrt{1+\hslash {\Omega}_{2}\left({H}_{c}\right)x/{E}_{F2}}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times {\displaystyle \underset{n=0}{\overset{{n}_{u2}\left(\theta ,{H}_{c},r,y\right)}{\sum}}\frac{\mathrm{tanh}\left[{A}_{3}\left(\mathrm{...}\right)\right]+\mathrm{tanh}\left[{B}_{3}\left(\mathrm{...}\right)\right]}{x+n+1/2}},\end{array}$ (19)

where

${\mathcal{l}}_{2}\left(\theta ,{H}_{c},r,y\right)=\frac{{E}_{1}\left(\theta \right)\left(1-1/ry\right)}{3\hslash {\Omega}_{2}\left({H}_{c}\right)}$ (20)

${n}_{u2}\left(\theta ,{H}_{c},r,y\right)=\text{floor}\left[2{\mathcal{l}}_{2}\left(\theta ,{H}_{c},r,y\right)-1/2\right]$ (21)

${A}_{3}\left(\mathrm{...}\right)={A}_{3}\left(\theta ,{T}_{c},{H}_{c},r,y,x,n\right)=\left[\frac{\hslash {\Omega}_{2}\left({H}_{c}\right)}{2k{T}_{c}}\left\{x+n+1/2+\frac{{E}_{1}\left(\theta \right)}{ry\hslash {\Omega}_{2}\left({H}_{c}\right)}\right\}\right]$ (22)

${B}_{3}\left(\mathrm{...}\right)={B}_{3}\left(\theta ,{T}_{c},{H}_{c},r,y,x,n\right)=\left[\frac{\hslash {\Omega}_{2}\left({H}_{c}\right)}{2k{T}_{c}}\left\{x+n+1/2-\frac{{E}_{1}\left(\theta \right)}{ry\hslash {\Omega}_{2}\left({H}_{c}\right)}\right\}\right].$ (23)

Although the multiplier r of y in the above equations is unity for an elemental SC, it has been introduced for later convenience when we deal with a composite SC characterized by multiple θs. Taking the example of Bi-2212, we note that if the Debye temperature of the SC is θ_{0} and pairing takes place via the Ca or the Sr ions, then r_{1} = θ_{Ca}/θ_{0} = 1 (because θ_{Ca} = θ_{0}) for the former case and r_{2} = θ_{Sr}/θ_{0} ≠ 1 (because θ_{Sr} differs from θ_{0}) for the latter.

4. Equations for the Density of Superconducting Electrons, their Critical Velocity, Effective Mass, and Critical Current Density

If we assume that (7) corresponding to a set of {E_{F}_{1}, T_{c}_{1}, H_{c}_{1}, j_{c}_{1} = 0}-values and (19) corresponding to a set of {E_{F}_{2}, T_{c}_{2}, H_{c}_{2}, j_{c}_{2} ≠ 0}-values have provided a value of y, then we are enabled to obtain the following equations for the density of superconducting electrons n_{s}, P_{c}, the critical velocity v_{c} and j_{c}, where the requisite factors of
$\hslash $ and c have been inserted for the convenience of the reader.

$\begin{array}{l}{n}_{s}\left(s,{E}_{F2}\right)=\frac{1}{3{\pi}^{2}}{\left(\frac{2{m}^{*}{E}_{F2}}{{\left(\hslash c\right)}^{2}}\right)}^{3/2}={C}_{1}{\left(s{E}_{F2}\right)}^{3/2}\\ \text{}\left({C}_{1}=\frac{{\left(2{m}_{e}{c}^{2}\right)}^{3/2}}{3{\pi}^{2}{\left(\hslash c\right)}^{3}}\cong 4.541\times {10}^{21}{\left(\text{eV}\right)}^{-\text{3/2}}{\text{cm}}^{-\text{3}}\right)\end{array}$ (24)

${P}_{c}\left(\theta ,{E}_{F2},s,y\right)=k{\left(6s{m}_{e}\right)}^{1/2}\frac{\theta}{y}{E}_{F2}^{-1/2},$

whence

$\begin{array}{l}{v}_{c}\left({\theta}_{0},{E}_{F2},s,y\right)=\frac{c{P}_{c}}{2s{m}_{e}}={C}_{2}\frac{{\theta}_{0}}{y}{\left(s{E}_{F2}\right)}^{-1/2}\\ \text{}\left({C}_{2}=\frac{{6}^{1/2}ck{m}_{e}^{-1/2}}{2}\cong 4.426\times {10}^{3}{\text{eV}}^{\text{1/2}}{\text{K}}^{-}{}^{\text{1}}{\text{cms}}^{-\text{1}}\right)\end{array}$ (25)

$\begin{array}{l}{j}_{c}\left({\theta}_{0},{E}_{F2},s,y\right)=e{n}_{s}{v}_{c}={C}_{3}\frac{{\theta}_{0}}{y}s{E}_{F2}\\ \text{}\left({C}_{3}={\text{eC}}_{\text{1}}{\text{C}}_{\text{2}}=\frac{2\sqrt{3}eck{m}_{e}}{3{\pi}^{2}{\left(\hslash c\right)}^{3}}\cong 3.22\times {10}^{6}{\text{CK}}^{-\text{1}}{\text{(eV)}}^{-\text{1}}{\text{cm}}^{-\text{2}}{\text{s}}^{-\text{1}}\right)\end{array}$ (26)

and

$s\left(\theta ,{E}_{F2},{j}_{c},r,y\right)=\frac{{j}_{c}ry}{{C}_{3}\theta {E}_{F2}}.$ (27)

5. Bi_{2}Sr_{2}CaCu_{2}O_{8} (Bi-2212)

5.1. Pairing Equation in the 2PEM Scenario for Non-Zero Values of Temperature and Applied Field when the Critical Current Density Is Zero

We now undertake a study of Bi-2212 for which T_{c} = 86 K when H = 0. When H ≠ 0, it is characterized by [19]

${T}_{c1}=65\text{K},\text{}{H}_{c1}=36\text{Tesla},\text{}{j}_{c1}=0.$ (28)

Among the 14 different values of j_{c} that have been reported for the same SC at different values of T and H, we focus here on the following two which are particularly interesting because they correspond to the same values of T and H:

${T}_{c2}=4.2\text{K},\text{}{H}_{c2}=12\text{Tesla},\text{}{j}_{c2}=2.4\times {10}^{5}\text{A}/{\text{cm}}^{\text{2}}$ (29)

${T}_{c3}={T}_{c2},\text{}{H}_{c3}={H}_{c2},\text{}{j}_{c3}=1.0\times {10}^{6}\text{A}/{\text{cm}}^{\text{2}}.$ (30)

For the data in (28), we need to employ (7) which corresponds to j_{c} = 0, whereas for the data in (29) and (30) we need (19). In either case, first of all, we require the value of the Debye temperature θ_{0} of the SC, which is 237 K. Since we are now dealing with a composite SC [(7) was derived for an elemental SC], we have the following multiple Debye temperatures at play because of the anisotropy of the SC [20]:

$\begin{array}{l}\theta \left(\text{Bi-2212}\right)={\theta}_{0}=237\text{K}\\ {\theta}_{\text{Ca}}=237,\text{}{\theta}_{\text{Sr}}=286\text{K},\text{}{\theta}_{\text{Bi}}=269\text{K}.\end{array}$ (31)

Note that y in (18) is defined in terms of the Debye temperature of the SC (not in terms of the Debye temperature of any of the constituents of the SC if it consists of several elements). Sticking to the same definition, we have

$y\left(\text{Bi-2212}\right)={E}_{1}\left({\theta}_{0}\right)/{E}_{2}\left({E}_{F2}\right)=\frac{k{\theta}_{0}}{{P}_{c}}\sqrt{\frac{6s{m}_{e}}{{E}_{F2}}},$ (32)

and, in all the four equations following (19), we need to have [15] r = θ_{Ca}/θ_{0}, θ_{Sr}/θ_{0} and θ_{Bi}/θ_{0} for pairing via the Ca, Sr and Bi ions, respectively.

Following our study dealing with the T_{c}s, Δ_{0}s and j_{0}s of Bi-2212 reported in [20], it seems imperative that we generalize (7) and (19) to the case of pairing via the 2PEM scenario. To this end, it is convenient to define

$\begin{array}{l}{L}_{3}\left(\theta ,{H}_{c}\right)=\frac{k\theta}{3\hslash {\Omega}_{1}\left({H}_{c}\right)},\text{}{n}_{u3}\left(\theta ,{H}_{c}\right)=\text{floor}\left[2{L}_{3}\left(\theta ,{H}_{c}\right)-1/2\right]\\ {\varphi}_{1}\left({H}_{c},{E}_{F},x\right)=1/\sqrt{1+\hslash {\Omega}_{1}\left({H}_{c}\right)x/{E}_{F}}\end{array}$

${F}_{1}\left(\theta ,{T}_{c},{H}_{c},x\right)={\displaystyle \underset{n=0}{\overset{{n}_{u3}\left(\theta ,{H}_{c}\right)}{\sum}}\frac{\mathrm{tanh}\left[\frac{\hslash {\Omega}_{1}\left({H}_{c}\right)}{2k{T}_{c}}\left(x+n+1/2\right)\right]}{x+n+1/2}}$

${G}_{1}\left(\theta ,{T}_{c},{H}_{c},{E}_{F1}\right)={\displaystyle \underset{-{L}_{3}\left(\theta ,{H}_{c}\right)}{\overset{{L}_{3}\left(\theta ,{H}_{c}\right)}{\int}}dx\text{}{\varphi}_{1}\left({H}_{c},{E}_{F},x\right){F}_{1}\left(\theta ,{T}_{c},{H}_{c},x\right)},$

which lead to the following equation in the 2PEM scenario as a generalization of (7) for pairing via the Ca and Sr ions when j_{c} = 0:

$\begin{array}{l}Eq1\left({\theta}_{1},{\theta}_{2},{T}_{c1},{H}_{c1},{E}_{F1},{\lambda}_{m1},{\lambda}_{m2}\right)\\ \equiv \mathrm{Re}\left[1-{\lambda}_{m1}{G}_{1}\left({\theta}_{1},{T}_{c1},{H}_{c1},{E}_{F1}\right)-{\lambda}_{m2}{G}_{1}\left({\theta}_{2},{T}_{c1},{H}_{c1},{E}_{F1}\right)\right]=0.\end{array}$ (33)

For the θ-values in (31), we note that (i) solution of (33) with inputs from (28) leads to (a) the value of λ_{m}_{1} for pairing in the 1PEM scenario via the Ca ions by putting
${\theta}_{1}={\theta}_{\text{Ca}}$, λ_{m}_{2} = 0 and (b) the value of λ_{m}_{2} for pairing in the 1PEM scenario via the Sr ions by putting λ_{m}_{1} = 0,
${\theta}_{2}={\theta}_{\text{Sr}}$ and (ii) the need to invoke 2PEM will arise only if either of the so-obtained λ_{m}s exceeds the Bogoliubov limit of 0.5, because beyond this value the system becomes unstable.

5.2. Pairing Equation in the 2PEM Scenario for Non-Zero Values of Temperature, Applied Field and Critical Current Density

If we obtain the generalized version of (19) to deal with pairing in the 2PEM scenario following the same procedure as above, then with the inputs from (29) and (30) we have two equations besides (33). These three equations are not sufficient to deal with at least the six unknowns in our problem: λ_{m}_{1}, λ_{m}_{2}, E_{F}_{1}, E_{F}_{2}, E_{F}_{3} and y. Therefore, ab-initio calculation of the j_{c}s via this approach will involve making ad hoc assumptions about some of the unknowns. To avoid such a situation, we adopt the strategy of incorporating in the generalized version of (19) the empirical j_{c}-values noted in (29) or (30). Even though the number of equations in this approach continues to be three, we show below how it sheds light on several important features of the data in (29) and (30).

In order to incorporate j_{c} into (19), we employ (27) for s(θ, E_{F}, j_{c}, r, y), whence we have (15) as

${\Omega}_{2}\left({H}_{c}\right)\to {\Omega}_{2}\left(\theta ,{H}_{c},{E}_{F},{j}_{c},r,y\right)=\frac{{C}_{3}{\Omega}_{0}{H}_{c}\theta {E}_{F}}{{j}_{c}ry}$

and (20)-(23) as

${\mathcal{l}}_{4}\left({H}_{c},{E}_{F},{j}_{c},r,y\right)=\frac{k\text{}{j}_{c}\left(ry-1\right)}{3{C}_{3}\hslash \text{}{\Omega}_{0}{H}_{c}{E}_{F}}$ (34)

${n}_{u}{}_{4}\left({H}_{c},{E}_{F},{j}_{c},r,y\right)=\text{floor}\left[2{\mathcal{l}}_{4}\left({H}_{c},{E}_{F},{j}_{c},r,y\right)-1/2\right]$ (35)

$\begin{array}{l}{A}_{3}\left(\mathrm{...}\right)={A}_{3}\left(\theta ,{T}_{c},{H}_{c},{E}_{F},{j}_{c},r,y,x,n\right)\\ =\left[\frac{\hslash {\Omega}_{2}\left(\theta ,{H}_{c},{E}_{F},{j}_{c},r,y\right)}{2k{T}_{c}}\left\{x+n+\frac{1}{2}+\frac{k\theta}{\hslash {\Omega}_{2}\left(\theta ,{H}_{c},{E}_{F},{j}_{c},r,y\right)}\right\}\right]\end{array}$ (36)

$\begin{array}{l}{B}_{3}\left(\mathrm{...}\right)={B}_{3}\left(\theta ,{T}_{c},{H}_{c},{E}_{F},{j}_{c},r,y,x,n\right)\\ =\left[\frac{\hslash {\Omega}_{2}\left(\theta ,{H}_{c},{E}_{F},{j}_{c},r,y\right)}{2k{T}_{c}}\left\{x+n+1/2-\frac{k\theta}{\hslash {\Omega}_{2}\left(\theta ,{H}_{c},{E}_{F},{j}_{c},r,y\right)}\right\}\right].\end{array}$ (37)

With the definitions

${\varphi}_{2}\left(\theta ,{H}_{c},{E}_{F},{j}_{c},x,r,y\right)=\frac{1}{\sqrt{1+\hslash {\Omega}_{2}\left(\theta ,{H}_{c},{E}_{F},{j}_{c},r,y\right)x/{E}_{F}}}$

${F}_{2}\left(\theta ,{T}_{c},{H}_{c},{E}_{F},{j}_{c},x,r,y\right)={\displaystyle \underset{n=0}{\overset{{n}_{u4}\left({H}_{c},{E}_{F},{j}_{c},r,y\right)}{\sum}}\frac{{A}_{3}\left(\mathrm{...}\right)+{B}_{3}\left(\mathrm{...}\right)}{x+n+1/2}},$

where ${\mathcal{l}}_{4}\left(\mathrm{...}\right)$, ${n}_{u4}\left(\mathrm{...}\right)$, ${A}_{3}\left(\mathrm{...}\right)$ and ${B}_{3}\left(\mathrm{...}\right)$ are given by (34), (35), (36) and (37), respectively, and

${G}_{2}\left(\theta ,{T}_{c},{H}_{c},{E}_{F},{j}_{c},r,y\right)={\displaystyle \underset{-{\mathcal{l}}_{4}\left(\mathrm{...}\right)}{\overset{{\mathcal{l}}_{4}\left(\mathrm{...}\right)}{\int}}\text{d}x\text{}}{\varphi}_{2}\left(\mathrm{...}\right){F}_{2}\left(\mathrm{...}\right),$ (38)

we obtain, in the 2PEM scenario as a generalization of (19) for pairing via two species of ions and a given value of j_{c}, the following equation:

$\begin{array}{l}Eq2\left({\theta}_{1},{\theta}_{2},{T}_{c},{H}_{c},q,{E}_{F},{j}_{c},{r}_{1},{r}_{2},y,{\lambda}_{m1},{\lambda}_{m2}\right)\\ \equiv \mathrm{Re}\left[1-{\lambda}_{m1}{G}_{2}\left({\theta}_{1},{T}_{c},{H}_{c},{E}_{F},{j}_{c},{r}_{1},y\right)-{\lambda}_{m2}{G}_{2}\left({\theta}_{2},{T}_{c},{H}_{c},{E}_{F},{j}_{c},{r}_{2},y\right)\right]=0.\end{array}$ (39)

5.3. Solutions in the Scenario of 1PEM Due to Either of the Ca, Bi and the Sr Ions

In order to unravel the empirical features of Bi-2212 noted in (28), (29) and (30) in the above framework, we proceed as follows:

1) We first deal with the data in (28) via (33) where, in the latter equation, E_{F1} is an independent variable. If we solve this equation for different assumed values of E_{F}_{1}= ρkθ_{0} and λ_{m}_{2} = 0, we obtain the corresponding values of λ_{m}_{1} for pairing via the Ca ions. We thus find that for ρ = 10, λ_{m}_{1} = 0.31283 and for all values of ρ > 10, λ_{m}_{1} has the same value up to four significant digits. For the same value of ρ and λ_{m}_{1} = 0, the corresponding value of λ_{m}_{2}(Bi) and λ_{m}_{2}(Sr) via (33) are 0.22043 and 0.20778, respectively Similarly obtained values of λ_{m}_{1}(Ca), λ_{m}_{2}(Bi) and λ_{m}_{2}(Sr) for some select values of
$\rho \le 10$ are given in Table 1.

2) Since none of the values of λ_{m} in Table 1 exceeds the Bogoliubov limit of 0.5, we now employ (39) also in the 1PEM scenario in order to determine E_{F}_{2} corresponding to the data in (29). To this end for ρ = 10, for pairing via the Ca ions, we use the following values in (39):

Table 1. Obtained via (33), the interaction parameters λ_{m}_{1}(Ca), λ_{m}_{2}(Bi) and λ_{m}_{2} (Sr) for pairing via the Ca, Bi, and Sr ions corresponding to T_{c}_{1} = 65 K, H_{c}_{1} = 36 × 10^{4} G with the additional inputs {λ_{m}_{2} = 0, θ_{1} = 237}, {λ_{m}_{1} = 0, θ_{2} = 269 K} and {λ_{m}_{1} = 0, θ_{2} = 286 K}, respectively.

$\begin{array}{l}{\lambda}_{m2}=0,\text{}{\lambda}_{m1}=0.31283;\text{}{T}_{c}={T}_{c2},{H}_{c}={H}_{c2},\text{}{j}_{c}={j}_{c2};\\ r={r}_{1}=1,\text{}{E}_{F}=q{E}_{F1}\text{}\left({E}_{F1}=\rho k{\theta}_{0}\right),\end{array}$ (40)

which leaves out y which is yet to be specified. Since y too is an independent variable, we need to solve (39) for a range of values of ry > 1; that ry must be greater than unity follows from (34). We find that for each such value of y, the solution of the transcendental Equation (39) yields, in general, multiple roots each of which corresponds not only to j_{c}_{2} (which is an in-built feature of our formalism), but also, to the accuracy with which they are quoted in Table 2 and Table 3, the same values for n_{s} and v_{c} via (24) and (25), respectively. For this reason, corresponding to each y, given in Table 2 are only the greatest roots of (39), which are found by first plotting (39) against q to determine the range outside which the equation does not have any roots.

In the context of the critical current densities, such a plot given in Figure 1 seems to be unusual and will be discussed below. The largest value of q corresponding to any value of y at which the function being plotted crosses zero is then found more accurately using numerical root-finding methods and is given in Table 2 where similar results for ρ = 5 and 1 for both pairings via the Ca and the Sr ions are also given.

The parameters, the values of which are different for different roots, corresponding to the same value of y are n_{u}_{4}, s and E_{F}_{2} where n_{u}_{4} and s are determined via (35) and (27), respectively, and E_{F}_{2} = qE_{F}_{1}. Of these, the values of s dictate the range of y relevant for the data under consideration. To elaborate, for y ≤ 2.40, it is found that s ≥ 37. Since, contrary to the known features of the SC under consideration, such values of s would put it in the category of heavy-fermion SCs [13], they must be excluded and so we set the lower limit of y at a value which leads to s ≈ 10. The value of s decreases as y is progressively increased. This feature enables us to set the upper limit on the value of y as one for which s ≈ 1.

3) For the data in (30), we also give in Table 3 the results that follow when the operative 1PEM is due to the Sr ions.

Table 2. The largest values of q ≡ E_{F}_{2}/E_{F}_{1} corresponding to y_{1} and y_{2} obtained by solving (39) for the set S_{1} = {T_{c}, H_{c}, j_{c}, ion species} as noted and different pairs of S_{2} = {ρ, λ_{m}} values taken from Table 1. n_{u}_{4}, s, n_{s} and v_{c} corresponding to q_{1} and q_{2} are calculated via (35), (27), (24) and (25), respectively.

Table 3. The largest values of q ≡ E_{F}_{2}/E_{F}_{1} corresponding to y_{1} and y_{2} obtained by solving (39) for the set S_{1} = {T_{c}, H_{c}, j_{c}, ion species} as noted and different pairs of S_{2} = {ρ, λ_{m}} values taken from Table 1. n_{u}_{4}, s, n_{s} and v_{c} corresponding to q_{1} and q_{2} are calculated via (35), (27), (24) and (25), respectively.

4) Pairing in Bi-2212 can of course also take place via the Bi ions in the 1PEM scenario. From the results obtained via the Ca and the Sr ions and the value of θ_{Bi}, which lies between θ_{Ca} and θ_{Sr} [vide (31)], one would expect that λ_{m}(Bi) for any of the value of ρ in either of the two Tables should lie between the corresponding values of λ_{m}(Ca) and λ_{m}(Sr) and that, but for this change, the overall results should be substantially similar to those already obtained. This is indeed found to be so, as can be seen from the following results corresponding to T = 4.2 K, H_{c} = 12 × 10^{4} G and j_{c} = 2.4 × 10^{5} A/cm^{2} as an example:

Figure 1. Plot of Eq 2 (θ_{1}, θ_{2}, T_{c}_{2}, H_{c}_{2}, q, E_{F}_{1}, j_{c}_{2}, r_{1}, r_{2}, y, λ_{m}_{1}, 0)vide (39) as a function of q showing multiple roots of (39) in the indicated range corresponding to θ_{1} = 237, r_{1} = 1, ρ = 10, λ_{m}_{1} = 0.31283, y = 2.60 and {T_{c}_{2}, H_{c}_{2}, j_{c}_{2}} as in (29).

$\begin{array}{l}\rho =10,\text{}\lambda =0.22043;\text{}2.54\le y\le 3.30\\ 4.0226\times {10}^{-4}\le q\le 5.4521\times {10}^{-3},\text{}70\ge {n}_{u}\ge 7,\text{}9.72\ge s\ge 0.93\\ 1.03\times {10}^{17}\le {n}_{s}\le 1.52\times {10}^{17},\text{}1.46\times {10}^{7}\ge {v}_{c}\ge 9.87\times {10}^{6}\\ 0.82\times {10}^{-4}\le {E}_{F2}\le 11.1\times {10}^{-4},\end{array}$

where the units for n_{s}, v_{c }and E_{F} are cm^{-3}, cm/s and eV, respectively.

6. Discussion

To give a perspective of the approach followed in this paper, it is pertinent to point out that, conventionally, the j_{c} of an SC is determined via one or the other critical state models; for characteristic equations of nine critical state models, see ( [19], p61). It is postulated in such models that for low applied fields or currents the outer part of the sample is in the so-called critical state which is characterized by particular values of j_{c} and H and that the interior of the SC is shielded from these fields and currents. For Bi-2212, the most commonly employed model is Bean’s model where its j_{c} is determined via the geometry of the sample and the magnetization width ΔM of the M(H) hysteresis loop; see, e.g. [21], which gives values of j_{c}(H) of the melt quenched and the non-melted samples of the SC at 5 K. Following the conventional approach, the significantly different j_{c }(H)-values of the samples in [21] are attributed to material properties of the samples such as their cell parameters, alignment and inter-connectivity of the grains and the grain boundaries. It is hence seen that the approach followed in our paper differs radically from the conventional approach.

A remark about the operator Re in (33): The function G_{1}(θ_{1}, T_{c}_{1}, H_{c}_{1}, E_{F}_{1}) in this equation becomes pure imaginary for all values of
$x<-{E}_{F}/\hslash {\Omega}_{1}\left({H}_{c}\right)$ because of
${\varphi}_{1}\left({H}_{c1},{E}_{F1},x\right)$. This is a situation which also occurs in several other problems, e.g., while dealing with heavy-fermion SCs [13] and BCS-BEC crossover [22]. In order then to obtain real solutions, one alternative is to manually shift the lower limit of integration. This becomes cumbersome if one is simultaneously dealing with more than one such equation. The operator Re provides a much simpler, one-step, alternative, as was also noted in the context of Fe-based SCs [23]. This remark also applies to (39).

As was noted above, the solution of (39) leads to multiple roots for any value of y for which s falls in the range of our interest. These are shown in Figure 1 for a particular value of y. A notable feature of this figure is its saw-tooth appearance (a series of “Vs”), which is attributable to a combination of the fluctuations in Fermi energy and the floor function employed in our formalism. The unmarked vertical limbs of the Vs are discontinuities which occur when the summation index n changes discretely from one integral value to another due to the floor function. We note that no root is found even when such limbs cross zero, as also that the saw-tooth behavior is not seen when (33) is solved despite the fact that it too employs the floor function, which is so because it is solved for a fixed value of E_{F}.

Another feature of (39) is that as y is progressively increased, the number of its roots keeps decreasing, as an example of which we note that as one goes from y = 2.60 to 3.10 in the upper-half of Table 2, the number of roots changes from 7 to 2. This is reminiscent of Melde’s experiment with stationary waves on a taut string clamped at both ends which are therefore nodal points. In-between these, the string displays a varying number of nodes depending on the frequency of the tuning fork with which it is induced to vibrate. One is thus led to surmise that: 1) (39) embodies content that leads to behavior akin to a Melde’s string, 2) the variation of y in the present case is equivalent to changing the frequency of the tuning fork and 3) the values of q at which (39) vanishes (i.e., the roots of the equation) are equivalent to the nodes of the Melde’s string.

Regarding the employment of Ω_{0} (which corresponds to free electron mass) in obtaining solutions of (33), we note that there is no loss of generality in making this assumption in so far as the solutions of this equation are employed to shed light on why the same SC at the same values of T and H sustains different values of j_{c}—which is the chief objective of this study. This assertion follows from the fact that different pairs of {ρ, λ_{m}} in both Table 2 and Table 3 lead to essentially similar results for each of the following parameters:
${n}_{u4},s,\text{}{n}_{s},{v}_{c}$ and
${E}_{F}$.

In so far as the n_{u} values given in Table 2 and Table 3 are concerned, we recall from [17] that the radius of the largest Landau orbit is given by

${r}_{n}\left(H\right)={r}_{0}\sqrt{2{n}_{u}\left(H\right)+1}$

where
${r}_{0}=\sqrt{\hslash c/e{H}_{0}}$ and H_{0} is the critical applied field at T = 0 and that n_{u}(H_{0}) for the elemental SCs Hg, In, Tl and Sn, for example, has the values 1062, 1900, 2201 and 3172, respectively. Since r_{n}(H) and therefore n_{u} may also be regarded as a measure of the coherence length of the SC, the low values of the latter in the two Tables signify that Bi-2212 has a much smaller coherence length than the elemental SCs, which is in accord with the known empirical facts.

If, for both the values of j_{c}, i.e., j_{c}_{2} = 2.4 × 10^{5} (for an Ag-sheathed tape) and j_{c3} = 1.0 × 10^{6} A/cm^{2} (for a multilayer tape), the SC is assumed to have nearly the same value of s, which is unlikely, then the results in Table 2 and Table 3 show that 1)
${n}_{s}\left({j}_{c2}\right)<{n}_{s}\left({j}_{c3}\right)$, 2)
${v}_{c}\left({j}_{c2}\right)>{v}_{c}\left({j}_{c3}\right)$ and 3)
${E}_{F}\left({j}_{c2}\right)<{E}_{F}\left({j}_{c3}\right)$. All of these results remain valid even if s(j_{c2}) ≈ 10 and s(j_{c3}) ≈ 1. However, if s(j_{c2}) ≈ 1 and s(j_{c}_{3}) ≈ 10 then, while the first two results still hold, we have E_{F}(j_{c}_{2}) > E_{F}(j_{c3}). In so far as the absolute values of E_{F} in Table 2 and Table 3 are concerned, the lowest among them are of the order of 0.3 meV or less; these are interpretable as being near the nodes or the node lines on the Fermi surface.

7. Conclusions

1) Based on the microscopic BSE customized to deal with superconductivity, we have given here derivation of equations that incorporate E_{F} (equivalently, the chemical potential μ), T, H, and P.

2) Among the main results of this paper are (33) and (39). The former generalizes the equation given in [17] that incorporated T and H, but not E_{F}; the latter generalizes the equation given in [18] that incorporated T, E_{F} and P, but not H.

3) Another notable result of the approach followed here is that it sheds light on why the cuprate that we have dealt with has much smaller coherence length than elemental SCs.

4) A novel aspect of our work is that it incorporates j_{c}, which is a sample-specific property, into the dynamical equations that govern pairing in the SC.

5) As is well-known, the j_{c} of an SC depends on several factors such as its size, shape, and how it is doped and prepared. Based on the premise that the E_{F} of an SC subsumes all of these features, we have given here a framework for testing it, and applied it for a detailed study of Bi-2212.

6) The upshot of the present work is that for greater substantiation of the above premise, there is a need not only to monitor via experiment, insofar as it may be feasible, the following parameters for Bi-2212: ${n}_{u4},s,{n}_{s}$ and ${v}_{c}$, but also to carry out similar studies for other SCs. Unfortunately, none of the conventionally employed critical state models sheds light on any of these parameters.

We conclude by noting that a detailed exposition of most of concepts of the BSE-based approach employed in this paper can be found in [24].

Acknowledgements

GPM thanks Dr. V.P.S. Awana for a valuable discussion of the conventional critical state models and the approach followed in this paper for determining the j_{c} of an SC. It is gratifying to note that invoking the electronic structure of an SC in this context, as done here, seemed to be a refreshing new approach to him.

Appendix A

In order to see that λ_{m} defined in (8) in the natural system of units is dimensionless, we need the following conversion factors:

$\begin{array}{l}1\text{g}={a}_{1}eV{c}^{-2},\text{}1\text{cm}={a}_{2}e{V}^{-1}\left(\hslash c\right),\text{}e={b}_{1}{\left(\hslash c\right)}^{1/2},\\ 1\text{s}={a}_{3}e{V}^{-1}\text{}\hslash ,\text{}1\text{G}={\text{g}}^{1/2}\cdot {\text{cm}}^{-1/2}\cdot {\text{s}}^{-1}={b}_{2}e{V}^{2}{\left(\hslash c\right)}^{-3/2},\end{array}$

where

$\begin{array}{l}{a}_{1}=5.60958616\times {10}^{32},\text{}{a}_{2}=5.06772886\times {10}^{4},\text{}{a}_{3}=1.51926689\times {10}^{15}\\ {b}_{1}={\left(137.03599895\right)}^{-1/2},\text{}{b}_{2}={a}_{1}^{1/2}{a}_{2}^{-1/2}{a}_{3}^{-1}\text{}=\text{}6.92507774\times {10}^{-2}.\end{array}$

Hence, we have

$e={b}_{1}{\left(\hslash c\right)}^{1/2}={b}_{1}\text{}\left(\text{if}\text{\hspace{0.17em}}\hslash =c=1\right)$ (A1)

${H}_{c}\left(BCS,G\right)={b}_{2}\left|{H}_{c}\right|e{V}^{2}{\left(\hslash c\right)}^{-3/2}={b}_{2}\left|{H}_{c}\right|e{V}^{2}\text{}\left(\text{if}\text{\hspace{0.17em}}\hslash =c=1\right)$ (A2)

${V}_{BCS}=\left|V\right|eVc{m}^{3}={a}_{2}^{3}\left|V\right|e{V}^{-2}{\left(\hslash c\right)}^{3}={a}_{2}^{3}\left|V\right|e{V}^{-2}\text{}\left(\text{if}\text{\hspace{0.17em}}\hslash =c=1\right)$ (A3)

It follows from the above equations that when
$\hslash \ne 1$ and
$c\ne 1$, V_{BCS} has to be divided by
${\left(\hslash c\right)}^{3}$ in order for it to play the required role of a propagator in the BSE-based approach and that, after inserting a factor of
$\hslash c$ on the RHS of the first equation in the set of Equations (2), we obtain

${\lambda}_{m}=\frac{{a}_{2}^{3}{b}_{1}{b}_{2}\left|H\right|\left|V\right|}{16{\pi}^{2}}\sqrt{\frac{2{m}^{*}{c}^{2}}{{E}_{F}}},$ (A4)

which is dimensionless, since both |H| and |V| are pure numbers.

We note that (A4) is needed only when one wishes to decouple the constituent V of λ_{m} after the latter has been determined corresponding to a given value of H.

References

[1] Tinkham, M. (1975) Introduction to Superconductivity. McGraw Hill, New York.

[2] Ibach, H. and Lüth, H. (1996) Solid State Physics. Springer, Berlin.

https://doi.org/10.1007/978-3-642-88199-2

[3] Bardeen, J. (1962) Critical Fields and Currents in Superconductors. Reviews of Modern Physics, 34, 667-681.

https://doi.org/10.1103/RevModPhys.34.667

[4] Bean, C.P. (1964) Magnetization of High-Field Superconductors. Reviews of Modern Physics, 36, 31-39.

https://doi.org/10.1103/RevModPhys.36.31

[5] Kim, Y.B., Hempstead, C.F. and Strand, A.R. (1963) Magnetization and Critical Supercurrents. Physical Review, 129, 528-535.

https://doi.org/10.1103/PhysRev.129.528

[6] Kupriyanov, M.Y. and Lukichev, V.F. (1980) Temperature Dependence of Pair-Breaking Current in Superconductors. Soviet Journal of Low Temperature Physics, 6, 210.

[7] Lee, D. (2012) Iron-Based Superconductors: Nodal Rings. Nature Physics, 8, 364-365.

https://doi.org/10.1038/nphys2301

[8] Zhang, Y., et al. (2012) Nodal Superconducting Gap-Structure in Ferropnictide Superconductor Ba2Fe2(As0.7P0.3)2. Nature Physics, 8, 371-375.

https://doi.org/10.1038/nphys2248

[9] Allan, M.P., et al. (2012) Anisotropic Energy Gaps of Iron-Based Superconductivity from Intraband Quasiparticle Interference in LiFeAs. Science, 336, 563-567.

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

[10] Lin, X., et al. (2013) Fermi Surface of the Most Dilute Superconductor. Physical Review X, 3, Article ID: 021002.

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

[11] Alexandrov, A.S. (2001) Nonadiabatic Polaronic Superconductivity in MgB2 and Cuprates. Physica C: Superconductivity, 363, 231-236.

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

[12] Jarlborg, T. and Bianconi, A. (2013) Fermi Surface Reconstruction of Superoxygenated La2CuO4 Superconductors with Ordered Oxygen Interstitials. Physical Review B, 87, Article ID: 054514.

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

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

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

[14] Malik, G.P. and Varma, V.S. (2015) A Study of Superconducting La2CuO4 via Generalized BCS Equations Incorporating Chemical Potential. World Journal of Condensed Matter Physics, 5, 148-159.

https://doi.org/10.4236/wjcmp.2015.53017

[15] Malik, G.P. (2016) On the Role of Fermi Energy in Determining Properties of Superconductors: A Detailed Comparative Study of Two Elemental Superconductors (Sn and Pb), a Non-Cuprate (MgB2) and Three Cuprates (YBCO, Bi-2212 and Tl-2212). Journal of Superconductivity and Novel Magnetism, 29, 2755-2764.

https://doi.org/10.1007/s10948-016-3637-5

[16] Malik, G.P. (2018) Correction to: On the Role of Fermi Energy in Determining Properties of Superconductors: A Detailed Comparative Study of Two Elemental Superconductors (Sn and Pb), a Non-Cuprate (MgB2) and Three Cuprates (YBCO, Bi-2212 and Tl-2212). Journal of Superconductivity and Novel Magnetism, 31, 941-941.

https://doi.org/10.1007/s10948-017-4520-8

[17] Malik, G.P. (2010) On Landau Quantization of Cooper Pairs in a Heat Bath. Physica B, 405, 3475-3481.

https://doi.org/10.1016/j.physb.2010.05.026

[18] Malik, G.P. (2013) On a New Equation for Critical Current Density in Terms of the BCS Interaction Parameter, Debye Temperature and the Fermi Energy of the Superconductor. World Journal of Condensed Matter Physics, 3, 103-110.

https://doi.org/10.4236/wjcmp.2013.32017

[19] Poole, C.P. (2000) Handbook of Superconductivity. Academic Press, San Diego.

[20] Malik, G.P. and Varma, V.S. (2019) Generalized BCS Equations: A Review and a Detailed Study of the Superconducting Features of Ba2Sr2CaCu2O8. In: Superconductivity and Superfluidity, IntechOpen, London, 21 p.

https://doi.org/10.5772/intechopen.84340

[21] Kumar, J., Sharma, D., Ahluwalia, P.K. and Awana, V.P.S. (2013) Enhanced Superconducting Performance of Melt Quenched Bi2Sr2CaCu2O8 (Bi-2212) Superconductor. Materials Chemistry and Physics, 139, 681-688.

https://doi.org/10.1016/j.matchemphys.2013.02.016

[22] Malik, G.P. (2014) BCS-BEC Crossover without Appeal to Scattering Length Theory. International Journal of Modern Physics B, 8, Article ID: 1450054.

https://doi.org/10.1142/S0217979214500544

[23] Malik, G.P. (2017) On the S±-Wave Superconductivity in the Iron-Based Superconductors: A Perspective Based on a Detailed Study of Ba0.6K0.4Fe2As2 via the Generalized Bardeen-Cooper-Schrieffer Equations Incorporating Fermi Energy. Open Journal of Composite Materials, 7, 130-145.

https://doi.org/10.4236/ojcm.2017.73008

[24] Malik, G.P. (2016) Superconductivity: A New Approach Based on the Bethe-Salpeter Equation in the Mean-Field Approximation. Volume 21 in the Series on Directions in Condensed Matter Physics. World Scientific, Singapore.

https://doi.org/10.1142/9868