Cooper Pairs Distribution function for bcc Niobium under pressure from first-principles

In this paper, we report Cooper Pairs Distribution function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{cp}(\omega ,{T}_{c})$$\end{document}Dcp(ω,Tc) for bcc Niobium under pressure. This function reveals information about the superconductor state through the determination of the spectral regions for Cooper-pairs formation. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{cp}(\omega ,{T}_{c})$$\end{document}Dcp(ω,Tc) is built from the well-established Eliashberg spectral function and phonon density of states, calculated by first-principles. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{cp}(\omega ,{T}_{c})$$\end{document}Dcp(ω,Tc) for Nb suggests that the low-frequency vibration region \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left(\omega <6 \,{\text{meV}}\right)$$\end{document}ω<6meV is where Cooper-pairs are possible. From \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{cp}(\omega ,{T}_{c})$$\end{document}Dcp(ω,Tc), it is possible to obtain the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${N}_{cp}$$\end{document}Ncp parameter, which is proportional to the total number of Cooper-Pairs formed at a temperature \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{c}$$\end{document}Tc. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${N}_{cp}$$\end{document}Ncp parameter allows an approach to the understanding of the Nb \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{c}$$\end{document}Tc anomalies, measured around 5 and 50 GPa.

www.nature.com/scientificreports/ The population of vibrational states (bosons) at temperature T(g n p ), with energy between ω and ω + dω , is determined by the density of vibrational states, N ph (ω) , times Bose-Einstein distribution at temperature T, where β = k B T and k B is the Boltzman constant. The distribution for an additional state with energy ω at temperature T is given by 11 , For electrons ( g n e ), the occupied states are described by electronic state density N(ε) , times Fermi factor at temperature T, A corresponding distribution to vacant electronic states is given by, An electron and a phonon in a crystal may or may not interact between them. The electron-phonon coupling probability is associated with the Eliashberg spectral function α 2 F(ω) . A measure of this probability is defined by the quotient α 2 (ω) = α 2 F(ω)/N ph (ω) . This quantity has been calculated for Nb at zero pressure by several authors 12,13 . Now, the probability that a pair of electrons is coupled by a phonon with energy between ω and ω + dω is obtained by simultaneous likelihood, summed over all electronic states, that one electron with energy ε interacts with the lattice and transfers to it an energy ω , taking the ε − ω state, and a second electron with energy ε ′ , by interaction with the lattice, absorbs the phonon of energy ω and goes to ε ′ + ω state, at the temperature T c , ergo: where D cp (ω, T c ) , namely Cooper Pairs Distribution function 9,10 , tells us the range of vibrational spectrum where the Cooper pairs are formed. Since electrons that form the Cooper pairs are near to the Fermi level E F , interacting through the lattice phonons, it implies that these are into E F − ω c < ε < E F + ω c energy interval. This cutoff phonon energy, ω c , should be the higher in the whole crystal-vibrational-spectrum (ω max ) . However, for T > 0 there are electrons occupying states with energy beyond Fermi level, then we could have ε > E F + ω max or ε < E F − ω max . Accordingly, we choose ω c such that to ω > ω c , electronic energies do not have an appreciable contribution to D cp (ω, T c ).
Further, we can get an estimate of the total number of Cooper pairs formed at the temperature T c through a quantity proportional to it, the N cp parameter, Next, we show the results of D cp (ω, T c ) and N cp evaluated at both pressure intervals that involve the anomalies at corresponding experimental T c values 3 .

Results and discussion
The Cooper Pairs Distribution functions D cp (ω, T c ) obtained for bcc Niobium under pressure (see Fig. 1) situates to Cooper-pairs formation in the 0-6 meV interval. It is observed that these D cp (ω, T c ) have a shape that mimics a Gaussian distribution centered around 1.5 meV. In general, the spectral Eliashberg function α 2 F(ω) moves towards high energies as pressure increases (see Fig. 2). However, in a narrow interval (0 to 4 meV) α 2 F(ω) spectra go back to low energy, as can be seen in Figs. 1 and 2. This behavior is in accordance with the work of Wierzbowska et al. 5,6 . By and large, our α 2 F(ω) spectra are in good agreement with previous theoretical reports [12][13][14] .
At the range pressures of 0−10 GPa , the superconducting critical temperature Niobium undergoes a decrease with an ulterior raising 3 , namely Nb anomaly. The tendency change happens at 4.0 GPa , approximately 3,15 . Nevertheless, electronic band structure and Fermi surface do not suggest any evident relationship with the Nb anomaly 4,16-24 . When we look at Eliashberg functions ( Fig. 2), we do not find a visible variation that we can associate with the Nb anomaly, except for the narrow interval from 0 to 4 meV. On a first view, this interval has a subtle or no contribution to Nb superconductivity. But this interval claims its importance in superconducting properties when is observed that Cooper pair formation (around 1.5 meV) coincides with this one (see Fig. 1).
Theoretical and experimental studies about Nb Fermi surface reported by several authors 4,16-23 do not detect any appreciable change at low pressure. However, Wierzbowska et al. 5,6 , by nesting factor, found slight features of the Nb Fermi surface. This could be associate with the Nb anomaly. On the other hand, Struzhkin et al. 3 and one previous study 9 showed that the electron-phonon parameter , has a good correlation with the critical www.nature.com/scientificreports/ temperature behavior. However, ours D cp (ω, T c ) specifically show that a small region of the whole vibrational spectrum contributes to the Cooper pair formation (Fig. 1). The D cp (ω, T c ) at pressures 0 ≤ P ≤ 10 GPa , are located around 1.5 meV, precisely in the spectral range where the corresponding α 2 F(ω) undergoes softening. Then, we could associate this demeanor with the Nb anomaly. Further, D cp (ω, T c ) marks the change of the tendency of T c (P) around 4.0 GPa . The D cp (ω, T c ) at 0 < P ≤ 4.5 GPa are always under D cp (ω, T c ) calculated at zero pressure (see the arrow down in Fig. 1a), while at 4.5 < P ≤ 10 GPa D cp (ω, T c ) are over it (see arrow up in Fig. 1b). This behavior allows us to infer that there is a change in physical conditions at ~ 4.0 GPa , which could be inducing the Nb anomaly.
On the other hand, the second Nb anomaly occurs around 50 GPa , where T c begins to decrease significantly 3 . This behavior has been explained by the electronic band structure and Fermi level variations and the trend of weakening electron-phonon interaction under pressure 3,4,7 . Our results show that Eliashberg functions harden at the pressure range of 45−75 GPa (Fig. 3a). A detailed view of α 2 F(ω) at low energies (Fig. 3b) displays a clear diminution of the electron-phonon interaction. In this case, the pressure effects are more evident than at low pressures.
All D cp (ω, T c ) calculated at the pressure range of 45−75 GPa reveal a Gaussian shape and are in the 0-5 meV interval (Fig. 3b). A significant decrease in D cp (ω, T c ) with increasing pressure is observed, mainly between 45 and 55 GPa where the Nb anomaly has been measured. Due to α 2 F(ω) decreasing as pressure increases (at low  www.nature.com/scientificreports/ energies), the Cooper pairs formation falls, and in that way, T c does too. So, as at low pressure (0−10 GPa) , the low energy phonons take remarkable importance in the behavior of T c . From D cp (ω, T c ) , it was possible to obtain the N cp parameter, which is proportional to the total number of Cooper-Pairs formed at a temperature T c . We found that the N cp parameters calculated (for both pressure intervals) have a significant correlation with experimental data for T c . It is observed that N cp as a function of pressure reproduce well the trend of T c (P) measured (see Fig. 4).
According to our results (Fig. 4), we can infer that the pressure induces modifications on the physical conditions (at low-energies) that lead to generation (or reduction) of the number of Cooper pairs, which cause the increase (or decrease) of the T c .
Finally, the analysis from D cp (ω, T c ) bring us to suggest an experiment where low-energy phonon could be stimulated in superconducting samples to research T c improvement. These results validate the use of D cp (ω, T c ) as a theoretical tool for the study of conventional superconductors.

Conclusions
Here, we presented the Cooper Pair Distribution functions D cp (ω, T c ) for Nb (bcc), calculated by first-principles. These results enabled us to broaden the understanding of the anomalous T c (P) behavior measured in the Nb, in the pressure ranges of 0−10 GPa and 45−75 GPa . D cp (ω, T c ) showed that the Cooper-pairs formation energy intervals are located at low-energies (ω < 6 meV) . At low-pressures (P < 10 GPa) , the low-energy region of  www.nature.com/scientificreports/ α 2 F(ω) moves towards lower frequencies in the same region of the Cooper pairs formation energy interval, which leads to the recognition of the importance of low-energy phonons in the superconducting behavior of Nb under pressure, despite their slight spectral weight. From D cp (ω, T c ) it was possible to obtain the N cp parameter, which is proportional to the total number of Cooper pairs formed at different T c and pressures. N cp as a function of pressure achieved an adequate reproduction of the trend of T c (P) measured, that is to say, the Nb anomalies around 5 and 50 GPa. From our D cp (ω, T c ) findings, it is expected that electron-phonon interaction specifically at low-energies (ω < 6 meV) contributes to superconducting properties, despite the α 2 F(ω) having a stronger spectral value at high vibrational energy. All of our results validate the use of D cp (ω, T c ) as a theoretical tool for the study of conventional superconductors.

Methods
In order to determine the Cooper Pair Distribution function D cp (ω, T c ) , we require electronic density states, vibrational density states, and Eliashberg function. To do these ab initio calculations, we first relax the internal degrees of freedom and the lattice vectors of the Nb structure using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton algorithm [25][26][27][28] at each pressure to get the corresponding lattice constants. From these relaxed structure configurations, we calculated the electronic and phonon band structures, electron (DOS) and phonon (PHDOS) densities of states, and the Eliashberg function α 2 F(ω) . We used a kinetic energy cut-off of 80Ry for the expansion of the wave function into plane waves and 320Ry for the density. To integrate over the Brillouin zone (BZ), we used for the electronic integration a k-grid of 32 × 32 × 32 and for the phononic integration a q-grid of 8 × 8 × 8 according to the Monkhorst-Pack scheme 29 . The calculations were done with the pseudopotential plane-wave (PW) method of Perdew et al. 30 , using the generalized gradient approximation (GGA) and the Troullier and Martins norm-conserving pseudopotential 31 . The cut-off and grids were chosen big enough as to obtain a good precision in α 2 F(ω) calculated within the density-functional perturbation theory (DFPT) frame 32,33 . We used the Quantum Espresso code 34 for all these calculations.