Superconductivity in high-entropy alloy system containing Th | Scientific Reports

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  13, Article number: 16317 (2023 ) Cite this article Tablets

Superconductivity in high-entropy alloy system containing Th | Scientific Reports

Th-containing superconducting high entropy system with the nominal composition (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) was synthesized. Its structural and physical properties were investigated by X-ray diffraction, scanning electron microscopy, energy dispersive X-ray spectroscopy, specific heat, resistivity and magnetic measurements. Two main phases of alloy were observed: major bcc structure and minor fcc. The experimental results were supported by numerical simulation by the DFT Korringa-Kohn-Rostoker method with the coherent potential approximation (KKR-CPA).

High entropy alloys (HEA) can be defined as solid solutions of five or more mixed elements in non-negligible amounts (higher than 5 at.% each), which are characterized by well-defined, ordered crystal structure with high chemical disorder, i.e. random distribution of elements at equivalent crystallographic positions. Instead of forming binary, ternary or other intermetallic compounds, they retain simple, closed packed structures similar to those of simple metals. The name ’high entropy’ comes from the large change in configurational entropy during the synthesis of a multi-element alloy, as described by Yeh et al.1.

Among other things, the HEAs have become well known as alloys with unusually high mechanical strength and corrosion resistance2,3,4. However, superconducting alloys with high entropy have attracted much attention in recent years. Most of them consists of tantalum-niobium matrices doped with other transition metals, such as Ti, Zr and Hf5,6,7. It has been shown, that chemical composition and molar ratio of the constituent elements strongly affect the properties of the superconducting state (SC) in terms of critical temperature \(T_\text{c}\) , critical field \(H_\text{c}\) , and SC transition width8,9,10,11,12.

So far, very little is known about actinides-containing high entropy alloys. Currently, HEA with uranium are mainly being studied in the development of advanced high strength materials13,14. Recently, however, a superconducting state has also been found in one of the alloys, namely [TaNb]\(_{0.31}\) (TiUHf)\(_{0.69}\) 8, but to the best of our knowledge there are no known superconducting HEAs containing thorium. The main goal of our research is to investigate new high entropy systems with higher critical parameters. Thorium metal would be an interesting addition to Nb-rich HEA, as it has been reported that eutectic Nb-Th system can show much enhanced superconductivity compared to isolated metals15. It would also be interesting, from a fundamental research point of view, to study the magnetic behavior of the known HEA superconducting systems when Th is introduced. Thorium possesses 5f-electron states, which can cause resulting alloy to exhibit complex magnetic behavior16,17,18. Motivated by this state of affairs, we have undertook a search for the first Th-based superconducting HEA. We chose the recently described superconducting alloy (NbTa)\(_{0.67}\) (MoHfW)\(_{0.33}\) 19 as our starting point. In addition to the synthesis of a new alloy, followed by full structural and physical characterization, we also undertook numerical studies of the electron structure of the obtained alloy. The DFT calculation were also used to explain the formation of more than one phase in studied sample.

X-ray powder diffraction pattern of (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) along with results of Rietveld refinement of its crystal structure. Red circles and black line represent the experimental data and fit curve, respectively, and green line is the difference between the two. Symbols \(*\) , \(\blacklozenge\) , \(+\) and \(\swarrow\) mark Bragg reflections coming from phases found in the sample and a pie chart shows atomic content of those phases.

SEM micrographs with EDXS elemental mapping of sample surface of (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) . Individual panels are labeled with symbols the content of which was analyzed; panel without the indicated element name shows the sample surface.

Fig. 1 presents powder X-ray diffraction (XRD) patterns obtained for the synthesized alloy (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) . Analysis of the experimental data by the Rietveld method showed that the sample consist of four phases: (1) the expected alloy crystallizing in a cubic, body-centered (bcc) structure (87.2 at.%, assuming the atomic distribution from EDS measurements - discussed later in this paper), (2) the expected alloy, but with a face-centred (fcc) structure (11.6 at.%), (3) thorium oxide (0.4 at.%), and (4) tungsten (0.8 at.%). The thorium oxide (ThO\(_2\) or thoria) is a refractory oxide with 3573 K melting point, which is a common impurity in thorium-containing compounds and alloys. Its presence in the studied sample was therefore difficult to avoid and, so to speak, expected. The case is similar for tungsten, which has the highest melting points among metals. Therefore, the presence of a small amount of unreacted tungsten is not surprising even with repeated melting of the sample. Fortunately, neither thorium oxide nor free tungsten superconduct and order magnetically (at least in the temperature range studied)20, so such small amounts of these impurities had no significant effect on determining the physical properties of the majority phase. However, the influence of the presence of the fcc phase will be discussed in the description of individual properties.

The lattice parameter estimated for the bcc phase was found to be a = 3.2725(3) Å. It differs from the lattice parameters of any of the pure metals used in the synthesis, but is close to those of Nb (a = 3.3033 Å  according to database Powder Diffraction Data - PDF - ref. no. 34-370) and Ta (a = 3.3058 Å  according to PDF 4-788), the two majority elements of the alloy studied. In the case of the fcc phase, a is 5.0955(6) Å  which is close to the lattice parameter of pure thorium (a = 5.0722 Å  according to PDF 1-920), which also crystallizes in the fcc structure.

Analysis of the chemical composition of the surface of the synthesized sample (Fig. 2) clearly showed presence of two phases with slightly different composition: large grains of the primary (major) phase separated by thin layers of the secondary phase. A closer look at the distribution of individual elements in the sample reveals that the primary phase contains less Th and more Ta and W than the secondary one, while the content of Nb and Mo in both phases is much more similar.

The EDXS analysis of selected spots is shown in Table 1.

Although EDSX is a semi-quantitative method that provides only approximate data on chemical composition of particular phases (especially when they are not well separated), the data we obtained show that the secondary phase is indeed an alloy, and not just inclusions of pure thorium. In other words, our sample probably contains two HEA alloys, which is consistent with the results of crystal structure refinement mentioned above. Numerical simulation of several structures with various atomic compositions using Density Functional Theory (DFT) methods sheds more light on the possible stoichiometry of the detected alloys (it is discussed later in this paper).

(a) Temperature dependence of mass magnetization \(\sigma\) of (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) in an external magnetic field \(\mu _0 H\) , measured at low temperature in both zero-field cooling (ZFC) and field-cooling (FC) regimes; solid curves serve as guides for the eye, and \(T_\text{c}\) and \(T_\text{c, onset}\) mark superconducting transition in the alloy (see the text for details). (b) Field variation of \(\sigma\) measured at various temperatures in the superconducting state.

(a) Temperature variation of specific heat \(C_P\) of (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) ; insert shows \(C_P/\) T vs. \(T^2\) measured in various external magnetic fields \(\mu _0 H\) . Black solid curves are fits of Eqs. (1) (main panel) and (2) (insert) to the experimental data. (b) Electrical resistivity \(\rho\) of the alloy measured as a function of temperature in zero magnetic field; insert shows \(\rho (T)\) in various fields. Solid curves serve as guides for the eye.

Measurements of magnetic properties of (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) revealed that it has a weak, featureless and nearly temperature-independent magnetic susceptibility at least down to about 7 K. This means that both the bcc and fcc phase are Pauli-like itinerant paramagnets.

Magnetic measurements in zero field cooling (ZFC) regime (Fig. 3a) shown that sample undergoes transition from paramagnetic state to diamagnetic one at the low temperatures. This transition begins with the onset at about 7 K which starts to develop into sharper decline of mass magnetization \(\sigma\) at 5.64(2) K (\(T_c\) ). The sharp decline ends at 3.2 K but diamagnetic signal is progressively getting stronger with cooling. The measurement conducted in field cooling (FC) regime showed no diamagnetic transition, probably due to strong vortex pinning in type-II HEA superconductor. Figure 3b is presenting field dependent magnetisation curves measured at different temperatures. One can see that when temperature rises and disrupt fully developed Meissner state (showed in Fig. 3a as being below 3.2 K), value of \(\sigma\) drastically decreases but diamagnetic signal disappears only after heating above 7 K (T\(_{c, onset}\) ).

The \(C_p(T)\) dependence of (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) is presented in Fig. 4a. It bears no significant feature above 7 K and it follows Dulong- Petit law at RT. Experimental points of the curve measured in nominal field \(\mu _0 H\) =0 K were easily described by the conventional formula:

in which the first term describes the electron contribution to the specific heat according to the Sommerfeld model, and the second term is the phonon contribution to \(C_P(T)\) in terms of the Debye model, with the Sommerfeld coefficient \(\gamma\) and Debye temperature \(\Theta _\text{D}\) , as fitting parameters; R denotes the universal gas constant and r is the number of atoms in the formula unit, which in this case is equal to 1. Results show simple metallic behaviour of this alloy at the normal state. Insert in Fig. 4a presents \(C_p/\) T measured at several selected fields at low temperatures in function of \(T^2\) , reviling two inclines corresponding to T\(_{c, onset}\) and \(T_{c}\) previously observed in Fig. 3a. The level of incline is relatively small in comparison to \(\lambda\) -shaped superconductive transition below \(T_{c}\) .

Below 7 K function \(C_P/T\) of (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) in the normal state can be described by the \(T^3\) -Debye law (bold black line in Fig. 4a).:

with \(\gamma\)  = 3.31(4) mJ K\(^{-2}\)  mol\(^{-1}\) and \(\beta\)  = 0.1739(5) mJ K\(^{-4}\)  mol\(^{-1}\) . Due to the presence of a non-negligible amount of fcc phase both of those parameters should be treated as an average of two phases strongly dominated by the bcc phase. Using lattice specific heat coefficient \(\beta\) as fitting parameter to equation:

specific Debye temperature \({\Theta _\text{D}^\text{LT}}\) =223.59(2) K was found. The \(T_c\) and \({\Theta _\text{D}^\text{LT}}\) were used to calculate electron-phonon coupling from McMillan’s equation21:

where \(\mu ^{*}\) is the Coulomb repulsion constant. Taking \(\mu ^{*}\)  = 0.125 (a value commonly used for systems containing mostly d-electron elements), we obtained \(\lambda _\mathrm{el-ph}\)  = 0.72(4). That classifies this alloy in the range of intermediate electron-phonon coupled superconductors. To estimate conductive electron density of states at the Fermi level \(N(E_\text{F})\) = 1.4, Sommerfeld coefficient \(\gamma\) was used in following equation:

Then, density of non-interacting electrons \(N(E_\text{F})^*\) = 0.8(1)  states eV\(^{-1}\)  f.u.\(^{-1}\) was calculated based on following equation22:

Fig. 4b depicts electrical resistivity  \(\rho\) measured in nominal zero applied field. As shown, the electrical resistivity decreases with decreasing temperature in a manner characteristic of metals. At low temperatures the resistivity drops to zero, indicating superconducting properties of the compound under investigation. Here we define \(T_c\) at 10% of residual resistivity, which in nominal field is 7 K. The RRR ratio is equal to 2.54, and its low value is a result of the polycrystalline nature and high degree of structural disorder in the sample. Insert in Fig. 4b shows temperature dependent curves measured in various fields. Transition to zero resistance is visible even under applied nominal field higher than used in other methods.

The experimental curves of \(\mu _0 H_\text{c2}(T)\) obtained for (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) can be described by the Ginzburg-Landau (G-L) equation:

yielding \(\mu _0 H_\text{c2}^\text{mag}(0) = 0.73(4)\)  T (for data from Fig. 3b) and \(\mu _0 H_\text{c2}^\text{res}(0) = 0.98(2)\)  T (for data from Fig. 4b) as the least-squares fitting parameters (see the dashed and dotted lines in Fig. 5). The \(\mu _0 H_\text{c2}^\text{mag}(0)\) was also estimated using the full Werthamer-Helfand-Hohenberg (WHH) formalism for isotropic-gap BCS superconductors in the dirty limit, incorporating the spin-paramagnetic effect via the Maki parameter \(\alpha _\text{M}\) and the spin-orbit scattering constant \(\lambda _\text{SO}\) 23,24,25:

where \(\gamma \equiv \sqrt{(\alpha _\text{M} \overline{h})^2 - (\frac{1}{2} \lambda _{\text{SO}})^2}\) , \(\overline{h} = \frac{4}{\pi ^2} \frac{H_{\text{c2}}}{-dH_{\text{c2}}/dT}\) and \(t=\frac{T}{T_C}\) . Best fit to the experimental data was found for \(\lambda _\text{SO}\) equal to zero, what suggest that spin-orbit scattering has very weak effect on the value of the upper critical field. The value of Maki parameter in aforementioned fit was 0.07(1), and the \(\mu _0 H_\text{c2}^\text{mag}(0)\) obtained from WHH model was equal to 0.67(1) T In order to determine the pair-braking mechanism, the orbital upper critical field \(\mu _0 H_\text{c2}^\text{orb}\) =0.67(1) T (virtually equal to \(\mu _0 H_\text{c2}^\text{mag}(0)\) from WHH model) was estimated using \(\mu _0 H_\text{c2}^\text{mag}(0)\) data and the formula developed for type-II superconductors in the dirty limit scenario24,25:

where the value of \(\left[ \text{d}(\mu _0 H_\text{c2}(T))/\text{d}T \right] _{T = T_\text{c}}\) is given by the formula:

The Pauli limiting field \(\mu _0 H_\text{P}\) = 13.4(2) T was estimated from the equation26:

According to Eq. (2) electronic contribution \(\textit{C}_{el}\) to the specific heat can be estimated by simply subtracting \(\beta\) \(\textit{T}^3\) - the lattice contribution. The resulting \(\textit{C}_{el}\) (T) dependency can be described using the formula from the Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity:

where \(\Delta _{0}\) is a superconducting energy gap, and A is a constant. The least-squares fit of this equation to the experimental data yielded \(\Delta _{0}/k_\text{B}\)  = 7.4(3) K. Using the value of superconductive gap and Sommerfeld coefficient per volume (derived from the known mass of the sample and the estimated density of the alloy equal to 12903 kgm\(^{-3}\) ) \(\gamma _{v}\) =304.6(3.7) J m\(^{-3}\) K\(^{-2}\) , it is possible27 to calculate the thermodynamic critical field \(\mu _{0}H_\text{c}(0)\) :

as equal to 28.4(1) mT.

Using the values derived from Eqs. (9) and (13) it is possible to estimate the Ginzburg-Landau coherence length \(\xi _\text{GL}\) 28:

and the Ginzburg-Landau penetration depth \(\lambda _\text{GL}\) :

For the studied alloy, these values are \(\xi _\text{GL}\)  = 25.13 nm and \(\lambda _\text{GL}\)  = 326.52 nm, respectively. From them the Ginzburg-Landau parameter \(\kappa _\text{GL}\) , defined as:

The lower critical field \(\mu _0 H_\text{c1}(0)\) was estimated to be 3.4 mT and 4.0 mT respectively, using the equations29:

where \(\alpha\) is calculated from29.

Phase diagram of superconducting (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) obtained from magnetization and resistivity measurements. Dotted and dashed lines are obtained from Eq. (7) while solid line is an experimental fit to the WHH model.

From the collected data regarding superconducting state it is evident that this alloy posses two superconductive phases - one with \(T_c\) around 7 K and other with \(T_c\) at 5.6 K. The latter is most likely corresponding to the bulk, Th-depleted, bcc phase as its transition temperature is close to the other HEA of similar composition. Also depth of Meissner state as well as phase transition peak in \(C_p/\) T curve indicates that this is the signal of dominant phase. Overlapping of the phases signals makes it impossible to determine the exact \(T_c\) of bcc phase from \(C_p\) /T curve. On the other hand the diamagnetic signal of the phase with higher \(T_c\) depicted in Fig. 3 has only minimal incline so \(T_c\) of the main phase could be described as above. In the \(C_p\) /T high \(T_c\) phase transition is extremely broad. It is hard to attribute this to Th-rich phase crystallizing in fcc structure as one would expect sharper transition due to sheer amount of this phase. Instead, based on SEM-EDS observations, these signals could originate from minuscule phases with disturbed stochiometry that could be formed on the interphase between well crystallized, superconductive bcc phase and nonsuperconductive fcc phase as both of them are not well separated, but merging one into another. That would explain weak but broad transitions seen in Figs. 3 and 4a but also robust electrical resistance in Fig. 4b. Unusual width of transition in \(C_p\) /T can sometimes be related to a specific chemical composition9,19 but the inhomogeneity of the phases (seen in SEM-EDS micrographs) is most likely the case in this system. A link between mesoscopic inhomogeneity and such broad transitions have recently been reported in Ti-Hf-Nb-Ta-Re30 and Sc-Hf-Nb-Ta-Ti-Zr31 systems. However, the formation of nanostrucutres does not seem to affect the superconductive properties of HEA32. To further analyze the two phases in regards of their stability and superconductive properties, the DFT numerical simulations were performed.

Free energy changes calculated for both bcc and fcc phases in as a function of Th concentration; insert shows the equilibrium lattice parameter determined from the relaxation of the unit cell.

The experimental data evidently proved that our (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) sample contains two HEA. The major Th-depleted phase crystallizes in bcc structure while the secondary Th-rich phase crystallizes in fcc. In order to study the relative stability of bcc vs. fcc phases, the total energy computations were done using KKR-CPA method for these two phases varying the concentrations of the constituent elements of the alloy around their experimentally observed stechiometry. For each structure and composition, the equilibrium lattice parameter were derived from the relaxation of the unit cell. The free energy F at temperature T for bcc and fcc phases was computed using equation33:

where the ground state formation energy (\(E_\text{form}\) ) is the difference between the total energy of a given alloy \(E_\text{alloy}\) and sum of total energies of elemental bulks \(E_i\) , weighted over concentrations \(c_i\) , while S is the total entropy. In the case of HEA, S is dominated by the configurational entropy \(S_\text{conf}\) which in random solid system containing n elements is given by combinatorial formula1:

The calculations of F at \(T=0,300, 1000\) and 2000 K were performed for over 1000 bcc and 1000 fcc HEA structures in which the atomic concentration of each element in the NbTaMoWTh alloy varies in the range of 10% - 46% for Ta, 10% - 49% for Nb, 0.2% - 22% for Mo, 0% - 17% for W and 0% - 35% for Th. Compositions of considered phases were constructed randomly in a process similar to the method of acceptance-rejection sampling, where the atomic concentrations are drawn randomly for each element in the ranges mentioned above. If the concentrations summed up to 100% (with 1% margin) CPA-KKR calculations were performed. Insert in the Fig. 6 shows the equilibrium lattice parameters and F values computed for all taken into account HEA structures as function of Th concentration (\(c_\text{Th}\) ). The results presented here are for \(T=1000\) K, for other temperatures see the online supplementary information Figs. S1–S3. Since the crystal stability of the system is related to the minimum of the free energy, the obtained results show the preference of bcc structures for \(c_\text{Th}<5\%\) , while fcc one for \(c_\text{Th}\) between 5% and 32%. These findings are in agreement with XRD and SEM results which reveal the coexistence of bcc HEA phase with Th content close to 4 at.% and fcc HEA phase with about 25 at.% of Th. For higher Th concentrations (\(32\%<c_\text{Th} \le 35\%\) ), the difference of F derived for both considered structures becomes relatively small indicating that the bcc and fcc phases may coexist together (heterogeneous equilibrium). Alternative approach to the computation of the formation energy is described in the supplementary information. Results are presented on Supplementary Figs. S4–S7.

Electron density of states of of (a) bcc Nb\(_{0.38}\) Ta\(_{0.37}\) Mo\(_{0.11}\) W\(_{0.10}\) Th\(_{0.04}\) and (b) fcc Nb\(_{0.31}\) Ta\(_{0.25}\) Mo\(_{0.12}\) W\(_{0.07}\) Th\(_{0.25}\) alloys calculated by the KKR-CPA method. Total and partial atomic densities, color coded and weighted by their atomic concentrations.

The electronic band structures for all considered bcc and fcc NbTaMoWTh alloys were calculated using the KKR-CPA method. In the first step, we choose two particular HEA compositions (bcc Nb\(_{0.38}\) Ta\(_{0.37}\) Mo\(_{0.11}\) W\(_{0.10}\) Th\(_{0.04}\) and fcc Nb\(_{0.31}\) Ta\(_{0.25}\) Mo\(_{0.12}\) W\(_{0.07}\) Th\(_{0.25}\) ) which are close to those estimated using EDXS analysis. Total and partial atomic densities of states (DOS) of these two alloys are shown in Fig. 7. In the case of bcc alloy the shape of the total density of states (TDOS) is similar to this reported for the superconducting bcc HEA (NbTa)\(_{0.67}\) (MoHfW)\(_{0.33}\) 19. The Fermi level lies 0.54 eV above the nearest TDOS maximum while the largest contribution to TDOS in the Fermi-level region (\(E - E_\text{F}\) from -4 to 2 eV) comes from Nb and Ta atoms, due to their highest atomic concentrations in the alloy. For fcc Nb\(_{0.31}\) Ta\(_{0.25}\) Mo\(_{0.12}\) W\(_{0.07}\) Th\(_{0.25}\) , one can note larger TDOS at the Fermi energy \(N^\text{KKR}(E_\text{F})\) = 1.57  states eV\(^{-1}\)  f.u.\(^{-1}\) than 1.09  states eV\(^{-1}\)  f.u.\(^{-1}\) for the bcc Nb\(_{0.38}\) Ta\(_{0.37}\) Mo\(_{0.11}\) W\(_{0.10}\) Th\(_{0.04}\) . Since the DFT calculations which are performed in this work take into account only the non-interacting electrons, the obtained \(N^\text{KKR}(E_\text{F})\) values can be compared with experimentally determined \(N(E_\text{F})^*\) = 0.8(1)  states eV\(^{-1}\)  f.u.\(^{-1}\) . As one can notice, the theoretical value of \(N^\text{KKR}(E_\text{F})\) for bcc HEA is slightly higher than experimental one. In the case of fcc HEA, \(N^\text{KKR}(E_\text{F})\) is much higher than \(N(E_\text{F})^*\) and is similar to \(N(E_\text{F})\) . Therefore, taking into account Eq. (6) this result suggests that for fcc phase, \(\lambda _\mathrm{el-ph}\) could be close to 0.

In the next step, taking the calculated \(N^\text{KKR}(E_\text{F})\) values for all considered bcc and fcc NbTaMoWTh alloys, it is possible to predict the theoretical value of superconducting critical temperature (\(T_\text{c}^\text{KKR})\) . For each alloy \(T_\text{c}^\text{KKR}\) was estimated using McMillan’s formula (Eq. (4)). The Coulomb repulsion constant \(\mu ^{*}\) was obtained from the empirical relation34:

\({\Theta _\text{D}^\text{LT}}\) was calculated from the relation given by Moruzzi et al.35:

where B is the bulk modulus evaluated at the equilibrium Wigner-Seitz sphere radius \(S_0\) and M is the atomic mass of given HEA. Finally, \(\lambda _\mathrm{el-ph}\) parameter was estimated as follows19,22,36:

where the theoretical specific heat coefficient \(\gamma _\text{th}\) was obtained by inserting \(N^\text{KKR}(E_\text{F})\) value into relation Eq. (5) and the experimental specific heat coefficient \(\gamma _\text{exp}\) was approximated by a sum of specific heat coefficients of elemental bulks \(\gamma _i\) , weighted over concentrations \(c_i\) . The computed values of \(T_\text{c}^\text{KKR}\) for over 1000 bcc and 1000 fcc NbTaMoWTh structures containing up to 35 at.% of Th are presented in Fig. 8. It should be noted here that in the case of superconducting HEA, the critical temperatures calculated using KKR-CPA method are generally two time higher than the experimental ones. In particular, for Ta\(_{0.34}\) Nb\(_{0.33}\) Hf\(_{0.08}\) Zr\(_{0.14}\) Ti\(_{0.11}\) HEA, the experimental \(T_\text{c}\) = 7.3 K while theoretically predicted \(T_\text{c}^\text{KKR}\) = 15 K36. For (NbTa)\(_{0.67}\) (MoHfW)\(_{0.33}\) , \(T_\text{c}\) = 4.3 K and \(T_\text{c}^\text{KKR}\) = 8.8 K19. Since that, the results presented in Fig. 8 should be treated as a crude estimate of real critical temperatures. However, some general trends can be discussed. Firstly, up to \(c_\text{Th}=30\%\) , \(T_\text{c}^\text{KKR}\) calculated for bcc HEA are much higher than for fcc alloys with the same atomic composition. Secondly, above \(c_\text{Th}=10\%\) , \(T_\text{c}^\text{KKR}\) values rapidly decrease with Th concentration. This tendency is much more pronounced in the case of bcc HEA. Connecting these results with experimental data, one can state that the broad superconducting transition observed in magnetic and specific heat measurements is caused be inhomogeneous distribution of Th atoms in the bcc system which leads to a certain distribution of the \(T_\text{c}\) values. Finally, dividing \(T_\text{c}^\text{KKR}\) values by factor of 2, one can obtain the mean critical temperature for bcc HEA with \(c_\text{Th}<10\%\) close to 7-8 K which is comparable with \(T_\text{c, onset}\) observed experimentally, while for fcc HEA with \(c_\text{Th} \approx 25\%\) the mean critical temperature is close to 0 K and again this result is in agreement with experimental data.

Theoretical predictions of superconducting critical temperature (\(T_\text{c}^\text{KKR}\) ) for various NbTaMoWTh alloys containing up to 35 at.% of Th.

The synthesis and physical properties of the high-entropy system (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) were described. Refinement of the crystal structure as well as SEM imaging and mapping shows segregation on two major phases, one rich in Th (crystallizing in fcc structure) and other Th-depleted (crystallizing in bcc structure) as well as W and ThO\(_2\) impurities. Specific heat, magnetic and resistivity measurements let ascribe superconductive state with \(T_c\) = 5.64(2) K to bcc structure, while \(T_{c,onset}\) around 7 K as originating from minuscule, distorted phases forming on an interphase of fcc and bcc structures. This interpretation was supported by DFT numerical simulations which reveal that SC state in fcc phase would be very unlikely.

A polycristalline (NbTa)\(_{0.67}\) (MoWTh)\(_{0.33}\) sample was synthesised by conventional arc melting technique where the stochiometric amounts of pure elements and Ti-gettered Ar atmosphere were used. Process was repeated eight times. To ensure minimal amount of ThO\(_2\) in thorium metal, the raw material was polished with metal file and remelted few times before the synthesis of the alloy. The total weight loss after the synthesis was less than 0.3\(\%\) .

Crystal structure of the product was studied by powder X-ray diffraction (XRD) using a PANalytical X’pert Pro diffractometer with CuK\(\alpha\) radiation. The experimental XRD pattern was analyzed by the Rietveld method using the HighScore Plus software.

Chemical composition and phase composition of the sample were verified by energy dispersive X-ray spectroscopy (EDXS) using a FESEM FEI Nova NanoSEM 230 scanning electron microscope equipped with an EDAX Genesis XM4 spectrometer on polished surface of the cut specimen.

Magnetic properties of the alloy were studied in temperature range 1.72–300 K and in magnetic fields up to 20 kOe using a commercial Quantum Design MPMS-XL magnetometer. Heat capacity and electrical resistivity were measured from room temperature down to 1.8 K using a Quantum Design PPMS platform.

Free energy and electronic structure calculations were performed using the Korringa-Kohn-Rostoker (KKR) method, which in the case of disordered systems was implemented together with the coherent potential approximation (CPA), where a random arrangements of all elements is replaced by the ordered lattice representing an average over all possible configurations of the disordered lattice within the simple unit cell (bcc or fcc)37,38,39,40,41. In this study, the KKR-CPA method implemented in the AkaiKKR (MACHIKANEYAMA) package was used42,43,44. The Perdew-Burke-Ernzerhof exchange-correlation functional (PBE) was applied to construct the muffin-tin crystal potential in the semirelativistic approach45,46. The cutoff for the angular momentum was set to \(l_\text{max}\) = 3 and 256 k points were used to sample the irreducible part of the Brillouin zone during the self-consistent cycle and density of states calculations. Additionally, values up to 5216 k points were tested yielding virtually the same results. Atomic sphere approximation (ASA) was utilized in all calculations. For each considered HEA structure and composition, the lattice parameter a, bulk modulus B and its derivative \(B'\) of the crystal were derived using the Murnaghan equation of state47.

The data presented in this study are openly available in OSF repository at DOI 10.17605/OSF.IO/G4N6B (subsection: Superconductivity in high-entropy alloy system containing Th).

Yeh, J.-W. et al. Nanostructured high-entropy alloys with multiple principal elements: Novel alloy design concepts and outcomes. Adv. Eng. Mater. 6(5), 299–303. (2004).

Youssef, K. M., Zaddach, A. J., Niu, C., Irving, D. L. & Koch, C. C. A novel low-density, high-hardness, high-entropy alloy with close-packed single-phase nanocrystalline structures. Mater. Res. Lett. 3(2), 95–99. (2015).

Gludovatz, B. et al. A fracture-resistant high-entropy alloy for cryogenic applications. Science 345(6201), 1153–1158. (2014).

Article  ADS  CAS  PubMed  Google Scholar 

Chen, Y. Y., Hong, U. T., Yeh, J. W. & Shih, H. C. Selected corrosion behaviors of a Cu0.5NiAlCoCrFeSi bulk glassy alloy in 288\(^\circ\) C high-purity water. Scripta Materialia 54(12), 1997–2001. (2006).

Ishizu, N. & Kitagawa, J. New high-entropy alloy superconductor Hf21Nb25Ti15V15Zr24. Results Phys. 13, 102275. (2019).

Kim, G. et al. Strongly correlated and strongly coupled s-wave superconductivity of the high entropy alloy Ta1/6Nb2/6Hf1/6Zr1/6Ti1/6 compound. Acta Mater. 186, 250–256. (2020).

Article  ADS  CAS  Google Scholar 

Koželj, P. et al. Discovery of a superconducting high-entropy alloy. Phys. Rev. Lett. 113, 107001. (2014).

Article  ADS  CAS  PubMed  Google Scholar 

Nelson, W. et al. Superconductivity in a uranium containing high entropy alloy. Sci. Rep. 10, 4717. (2020).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

von Rohr , F. , Winiarski , MJ , Tao , J. , Klimczuk , T. & Cava , RJ Effect of electron count and chemical complexity in the ta-nb-hf-zr-ti high-entropy alloy superconductor .Proc.Natl.Acad.Sci.Rev. 113(46), 7144–7150. (2016).

Sun, L. & Cava, R. J. High-entropy alloy superconductors: Status, opportunities, and challenges. Phys. Rev. Mater. 3, 090301. (2019).

Gabovich, A. M., Li, M. S., Pekała, M., Szymczak, H. & Voitenko, A. I. Heat capacity of mesoscopically inhomogeneous superconductors: theory and applications to MgB2. Physica C 405(3), 187–211. (2004).

Article  ADS  CAS  Google Scholar 

Liao, T. F., Sung, H. H., Syu, K. J. & Lee, W. H. Alloying effects of Y on Tc in superconducting (La1-xY x)NiC2. Solid State Commun. 149(11), 448–452. (2009).

Article  ADS  CAS  Google Scholar 

Aizenshtein, M. et al. High entropy uranium-based alloys: Thermodynamics, characterization and mechanical properties. J. Nucl. Mater. 558, 153378. (2022).

Zhang, W., Li, Y., Liaw, P. K. & Zhang, Y. A strategic design route to find a depleted uranium high-entropy alloy with great strength. Metals (2022).

Cline, H. E., Rose, R. M. & Wulff, J. Niobium-thorium eutectic alloy as a high-field, high-current superconductor. J. Appl. Phys. 34(6), 1771–1774. (2004).

Moore, K. T. & van der Laan, G. Nature of the \(5f\) states in actinide metals. Rev. Mod. Phys. 81, 235–298. (2009).

Article  ADS  CAS  Google Scholar 

Maple, M., Baumbach, R., Butch, N., Hamlin, J. & Janoschek, M. Non-fermi liquid regimes and superconductivity in the low temperature phase diagrams of strongly correlated d- and f-electron materials. J. Low Temp. Phys. 161, 4–54. (2010).

Article  ADS  CAS  Google Scholar 

Pfleiderer, C. Superconducting phases of \(f\) -electron compounds. Rev. Mod. Phys. 81, 1551–1624. (2009).

Article  ADS  CAS  Google Scholar 

Sobota, P. et al. Superconductivity in the high-entropy alloy \({(\rm NbTa )}_{0.67}{(\rm MoHfW )}_{0.33}\) . Phys. Rev. B 106, 184512. (2022).

Article  ADS  CAS  Google Scholar 

Johnson, R. T., Vilches, O. E., Wheatley, J. C. & Gygax, S. Superconductivity of tungsten. Phys. Rev. Lett. 16, 101–104. (1966).

Article  ADS  CAS  Google Scholar 

McMillan, W. L. Transition temperature of strong-coupled superconductors. Phys. Rev. 167, 331–344. (1968).

Article  ADS  CAS  Google Scholar 

Tari, A. The Specific Heat of Matter at Low Temperatures (Imperial College Press, 2003).

Maki, K. Effect of pauli paramagnetism on magnetic properties of high-field superconductors. Phys. Rev. 148, 362–369. (1966).

Article  ADS  CAS  Google Scholar 

Werthamer, N. R., Helfand, E. & Hohenberg, P. C. Temperature and purity dependence of the superconducting critical field, \({H}_{c2}\) . iii. Electron spin and spin-orbit effects. Phys. Rev. 147, 295–302. (1966).

Article  ADS  CAS  Google Scholar 

Helfand, E. & Werthamer, N. R. Temperature and purity dependence of the superconducting critical field, \({H}_{c2}\) . ii. Phys. Rev. 147, 288–294. (1966).

Article  ADS  CAS  MATH  Google Scholar 

Clogston, A. M. Upper limit for the critical field in hard superconductors. Phys. Rev. Lett. 9, 266–267. (1962).

Parks, R. D. Superconductivity (Taylor and Francis, 1969).

Tinkham, M. Introduction to Superconductivity. Dover Books on Physics Series (Dover Publications, 2004).

Brandt, E. H. Properties of the ideal ginzburg-landau vortex lattice. Phys. Rev. B 68, 054506. (2003).

Article  ADS  CAS  Google Scholar 

Hattori, T. et al. Metallurgy, superconductivity, and hardness of a new high-entropy alloy superconductor Ti-Hf-Nb-Ta-Re. J. Alloys Metall. Syst. 3, 100020. (2023).

Krnel, M. et al. The effect of scandium on the structure, microstructure and superconductivity of equimolar sc-hf-nb-ta-ti-zr refractory high-entropy alloys. Materials (2022).

Article  PubMed  PubMed Central  Google Scholar 

Vrtnik, S. et al. Superconductivity in thermally annealed ta-nb-hf-zr-ti high-entropy alloys. J. Alloy. Compd. 695, 3530–3540. (2017).

Ikeda, Y., Grabowski, B. & Körmann, F. Ab initio phase stabilities and mechanical properties of multicomponent alloys: A comprehensive review for high entropy alloys and compositionally complex alloys. Mater. Charact. 147, 464–511. (2019).

Douglass, D. H. Superconductivity in D- and F- Band Metals Vol. Metals (Springer, 1976).

Moruzzi, V. L., Janak, J. F. & Schwarz, K. Calculated thermal properties of metals. Phys. Rev. B 37, 790–799. (1988).

Article  ADS  CAS  Google Scholar 

Jasiewicz, K., Wiendlocha, B., Korbeń, P., Kaprzyk, S. & Tobola, J. Superconductivity of Ta34Nb33Hf8Zr14Ti11 high entropy alloy from first principles calculations. Physica Status Solidi (RRL) Rapid Res. Lett. 10(5), 415–419. (2016).

Article  ADS  CAS  Google Scholar 

Soven, P. Coherent-potential model of substitutional disordered alloys. Phys. Rev. 156, 809–813. (1967).

Article  ADS  CAS  Google Scholar 

Gaspari, G. D. & Gyorffy, B. L. Electron-phonon interactions, \(d\) resonances, and superconductivity in transition metals. Phys. Rev. Lett. 28, 801–805. (1972).

Article  ADS  CAS  Google Scholar 

Butler, W. H. Theory of electronic transport in random alloys: Korringa-Kohn-Rostoker coherent-potential approximation. Phys. Rev. B 31, 3260–3277. (1985).

Article  ADS  CAS  Google Scholar 

Kaprzyk, S. & Bansil, A. Green’s function and a generalized Lloyd formula for the density of states in disordered muffin-tin alloys. Phys. Rev. B 42, 7358–7362. (1990).

Article  ADS  CAS  Google Scholar 

Bansil, A. & Kaprzyk, S. First-principles treatment of disorder effects in complex alloys: A study of \({\text{ Ba }}_{\mathit{x}}\) \({\text{ K }}_{1{-}\mathit{x}}\) \({\text{ BiO }}_{3}\) and \({\text{ BaPb }}_{1{-}\mathit{x}}\) \({\text{ Bi }}_{\mathit{x}}\) \({\text{ O }}_{3}\) . Phys. Rev. B. 43, 10335–10339. (1991).

Article  ADS  CAS  Google Scholar 

Akai, H. Electronic structure ni-pd alloys calculated by the self-consistent kkr-cpa method. J. Phys. Soc. Jpn. 51(2), 468–474. (1982).

Article  ADS  CAS  Google Scholar 

Akai, H. Fast Korringa-Kohn-Rostoker coherent potential approximation and its application to FCC Ni-Fe systems. J. Phys. Condens. Matter 1(43), 8045–8064. (1989).

Kotani, T. & Akai, H. KKR-ASA method in exact exchange-potential band-structure calculations. Phys. Rev. B 54, 16502–16514. (1996).

Article  ADS  CAS  Google Scholar 

Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868. (1996).

Article  ADS  CAS  PubMed  Google Scholar 

Perdew, J. P., Burke, K. & Wang, Y. Generalized gradient approximation for the exchange-correlation hole of a many-electron system. Phys. Rev. B 54, 16533–16539. (1996).

Article  ADS  CAS  Google Scholar 

Birch, F. Finite elastic strain of cubic crystals. Phys. Rev. 71, 809–824. (1947).

Article  ADS  CAS  MATH  Google Scholar 

R.T. acknowledges the support of Dioscuri program initiated by the Max Planck Society, jointly managed with the National Science Centre (Poland), and mutually funded by the Polish Ministry of Science and Higher Education and the German Federal Ministry of Education and Research. This research was funded in whole by National Centre of Science Grant number 2020/39/B/ST5/01782. For the purpose of Open Access, the author has applied a CC-BY public copyright licence to any Author Accepted Manuscript (AAM) version arising from this submission.

Institute of Experimental Physics, University of Wrocław, pl. M. Borna 9, 50-204, Wrocław, Poland

Piotr Sobota, Rafał Topolnicki, Tomasz Ossowski & Rafał Idczak

Institute of Low Temperature and Structure Research, Polish Academy of Sciences, ul. Okólna 2, 50-422, Wrocław, Poland

Piotr Sobota, Daniel Gnida & Adam Pikul

Dioscuri Center in Topological Data Analysis, Institute of Mathematics, Polish Academy of Sciences, ul. Śniadeckich 8, 00-656, Warsaw, Poland

Institute of Electronics and Information Technology, Lublin University of Technology, ul. Nadbystrzycka 38A, 20-618, Lublin, Poland

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

P.S.: Investigation, Formal analysis, Writing - Original Draft, Review and Editing; R.T.: Investigation, Formal analysis; T.O.: Investigation, Formal analysis; T.P.: Formal analysis; D.G.: Investigation, Formal analysis; R.I.: Formal analysis, Writing - Original Draft, Review and Editing, Funding acquisition; A.P.: Writing - Review and Editing.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Sobota, P., Topolnicki, R., Ossowski, T. et al. Superconductivity in high-entropy alloy system containing Th. Sci Rep 13, 16317 (2023).


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Superconductivity in high-entropy alloy system containing Th | Scientific Reports

Functional Materials Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.