Interplay of Chemical Disorder, Vortex Pinning, and Electron Correlations in the Mechanically Robust Medium-Entropy Alloy Superconductor $\text{Nb}_{0.45}\text{Ta}_{0.25}\text{Ti}_{0.15}\text{Hf}_{0.15}$

Introduction

Superconducting materials capable of sustaining high critical current densities while maintaining excellent structural and mechanical integrity are essential for next-generation technologies, including fusion reactors, high-field magnets, particle accelerators, quantum technologies, and advanced power transmission systems [1, 2, 3]. Although conventional superconductors have enabled many technological advances, achieving an optimal balance between superconducting performance, mechanical robustness, and environmental stability remains a persistent challenge. In this context, the emergence of complex entropy-stabilized alloys (high- and medium-entropy alloys) has transformed materials design by demonstrating that chemical complexity can yield combinations of properties that are difficult to achieve in conventional alloys [4, 5, 6]. Unlike conventional alloys that rely on one or two principal elements, high entropy alloys (HEAs) and medium entropy alloys (MEAs) contain four or more elements mixed in equiatomic ratios [5, 7, 8]. Thermodynamically, the stability of HEAs and MEAs anticipated by the high entropy of mixing according to the Gibbs free energy equation:

where T is the temperature, $\Delta H_{\mathrm{mix}}$ is enthalpy of mixing and $\Delta S_{\mathrm{mix}}$ is configurational entropy [9]. The configurational entropy increases logarithmically with an increase in constituent elements as $\Delta S_{\mathrm{mix}}^{\mathrm{conf}} = -R\sum_{i=1}^{N} c_i \ln c_i$, where $R$ is the gas constant, $N$ the number of components, and $c_i$ the atomic fraction of the $i$-th element. At elevated temperatures, the resulting $T\Delta S_{\mathrm{mix}}$ term becomes large enough to overrule the mixing enthalpy, thereby stabilizing simple solid-solution phases over intermetallic or complex phases according to Hume-Rothery criteria. Conventionally, alloys with $0.69R \le\Delta S_{\mathrm{mix}} < 1.61R$ are known as MEAs, whereas those with $\Delta S_{\mathrm{mix}} \ge1.61R$ are categorized as HEAs [10]. Most of the reported HEAs and MEAs crystallize in BCC, CsCl, HCP, and $\alpha$-Mn-type phases [11, 12, 13].

A defining characteristic of entropy-stabilized alloys is their exceptional mechanical performance, including high yield strength, high fracture toughness, superior microhardness, and remarkable fatigue resistance, which significantly outperform those of conventional alloys [14].

Moreover, their high damage tolerance and structural stability under intense particle bombardment make them attractive candidate for application in the severe radiation and harsh environments associated with nuclear reactors [15].

These exceptional structural features, combined with the intrinsic chemical complexity of entropy-stabilized alloys, have sparked growing interest in exploring superconducting properties, motivated by the need for mechanically robust superconducting materials for operating under extreme conditions. The first superconductivity was discoverd in $\mathrm{Ta}_{34}\mathrm{Nb}_{33}\mathrm{Hf}_{8}\mathrm{Zr}_{14}\mathrm{Ti}_{11}$, which is a type-II superconductor with $T_C = 7.3\,\mathrm{K}$ and upper critical field $= 8.2\,\mathrm{T}$ [16]. This pioneering work established that superconductivity can persist despite extreme chemical disorder observed in these systems. Since then, several refractory HEAs and MEAs have been reported, focusing on exploring new materials with improved performance and the potential for unconventional superconductivity [16, 17, 18, 19]. For instance, an equimolar ScTiVNbCu HEA, which exhibits a high $J_c$ of $0.18 \times10^{6}\,\mathrm{A}/\mathrm{cm}^{2}$ and a high upper critical field of $\mu_{0}H_{c2}$ of $12.9\,\mathrm{T}$ [18]. In eutectic NbScTiZr $J_c$ also exceeds $10^{6}\,\mathrm{A}/\mathrm{cm}^{2}$ while in the BCC $\mathrm{Nb}_{2/6}\mathrm{Ta}_{1/6}\mathrm{Ti}_{1/6}\mathrm{Zr}_{1/6}\mathrm{Hf}_{1/6}$ $J_c$ is of the order of $10^{5}\,\mathrm{A}/\mathrm{cm}^{2}$ at $2\,\mathrm{K}$ [20, 21, 22]. Thin films of NbTiTaHfZr exhibit a high $J_c$ exceeding $10^{6}\,\mathrm{A}/\mathrm{cm}^{2}$ at $4.2\,\mathrm{K}$ [23]. Equiatomic $\mathrm{Nb}_{0.25}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.25}\mathrm{Zr}_{0.25}$ shows multiband feature with high $J_c$ value of $\approx0.13 \times10^{6}\,\mathrm{A}/\mathrm{cm}^{2}$ [19]. Also, Ti-rich HEAs have been reported to exhibit a high $\mu_{0}H_{c2}(0) = 13.5\,\mathrm{T}$ [24].

Recently, exceptional fracture toughness even at cryogenic temperature has been observed in MEA $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$ [25, 75].This material highlights fracture toughness exceeding $253\,\mathrm{MPa}\cdot\mathrm{m}^{1/2}$ over a wide temperature range of $77$–$1473\,\mathrm{K}$, along with room-temperature yield strength of $596\,\mathrm{MPa}$. The remarkable deformation behavior in $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$ is attributed to the kink band formation driven by the coordinated slip of $\langle111\rangle$ edge dislocations across intersecting $\{110\}$ and $\{112\}$ slip systems, while edge dislocation pile up on $\{112\}$ planes, initiates the dislocation motion that causes kink band formation. The repeated nucleation and propagation of these kink bands accommodate heterogeneous deformation, redistribute local stresses away from crack tips, and mitigate strain localization, thereby establishing an intrinsic micromechanical hardening mechanism and enabling exceptional ductility and fracture resistance in this refractory BCC MEA [25].

Motivated by the need for mechanically robust superconducting materials, we investigate the superconducting properties and local distortion effects of mechanically robust $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$. In entropy-stabilized alloys, the observed severe lattice distortion, spatial fluctuations in the local atomic environment, and associated defects not only enhance mechanical robustness but can also serve as effective vortex-pinning centers, thereby suppressing vortex motion and improving the critical current density. In this study, we synthesized the MEA $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$ and conducted a systematic experimental and theoretical investigation to understand the phase stability, normal- and superconducting-state properties, and lattice-distortion effects in this material.

Methods

Experimental Methods

Sample Synthesis

An entropy-stabilized polycrystalline alloy of composition $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$ was synthesized using the arc melting technique. The constituent elements were weighed in stoichiometric ratios and then melted in a high-purity argon atmosphere in the presence of a Ti getter. The alloy was melted to form an ingot, then flipped and remelted several times to ensure chemical homogeneity of the sample. The sample formed was hard and exhibits negligible weight loss after melting .

Structural and Microstructural Analysis

The structural and microstructural analysis of the alloy was done using a Rigaku X-Ray Diffractometer (Rigaku SmartLab SE) with Cu-K$\alpha1$ radiation ($\lambda= 1.5406\,\mathrm{\AA}$) and a ZEISS GEMINI field-emission scanning electron microscope (FESEM), respectively. Elemental distribution was determined using energy-dispersive X-ray spectroscopy (EDS). The Le Bail refinement of the resulting X-ray diffraction (XRD) patterns was performed using the FULLPROF Suite.

Photoemission Spectroscopy

Photoemission spectroscopic measurements (PES), including X-ray Photo-electron spectroscopy (XPS), were performed on Thermoscientific Nexsa G2 series using a spot size of $400\,\mu\mathrm{m}$, whereas ultraviolet Photo-electron spectroscopy measurement (UPS) was done at AIPES beamline of the Indus-1 synchrotron radiation facility at Raja Ramanna Center for Advanced Technology (RRCAT), Indore, using a photon energy of $40\,\mathrm{eV}$ at room temperature. Prior to photoemission measurements, the alloy surface was etched by using $\mathrm{Ar}^{+}$ ions at 3000 eV for 200 s to remove the oxide layer and surface contaminations. XPS provides information on core-level binding energies and oxidation states, while UPS provides valence-band features near the Fermi level.

Image not extracted

Magnetic and Physical Properties

The magnetic properties were analyzed using a Superconducting Quantum Interference Device (SQUID) measurement system (MPMS-XL7). The electrical transport and specific heat measurements were carried out using the Physical Property Measurement System (PPMS) by Quantum Design. The conventional four-probe method is used for electrical transport measurements.

Computational Methods

To provide a microscopic interpretation of the experimentally observed crystal structure, metallic conductivity, and stability of the MEA Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$, density functional theory (DFT)-based first-principles calculations were carried out [26, 27]. Since the alloy exhibits intrinsic substitutional disorder arising from the random occupation of lattice sites by multiple constituent elements, an explicit treatment of chemical disorder is essential for obtaining a realistic description of its electronic structure. Hence, for a realistic description of this random occupation of the body-centered cubic (bcc) lattice while retaining periodic boundary conditions, we adopted the Special Quasirandom Structure (SQS) [28] approach. This approach is known to reproduce the most important correlation functions of an ideal random alloy and is therefore considered a well-established framework for investigating the electronic structure and lattice dynamics of chemically disordered high- and medium-entropy alloys. The SQS was generated using the mcsqs module of the Alloy Theoretic Automated Toolkit (ATAT), which employs a Monte Carlo simulated annealing algorithm to optimize the atomic configuration [29, 30, 31, 32]. During the optimization, the atomic arrangement was chosen to reproduce as closely as possible the pair and higher-order correlation functions of an ideal random solid solution. A 20-atom supercell was found to provide an optimal balance between computational efficiency and an accurate representation of the experimentally observed substitutional disorder. The optimized SQS model was subsequently employed for all electronic-structure and lattice-dynamical calculations. First-principles calculations were carried out based on the framework of density functional theory (DFT) [26, 27] using the projector augmented-wave (PAW) method [33, 34], as implemented in the Vienna $ab$-initio Simulation Package (VASP) [35, 36, 37]. The exchange-correlation effects were treated within the generalized gradient approximation (GGA) [38] using the Perdew-Burke-Ernzerhof (PBE) parametrization [39]. A plane-wave kinetic-energy cutoff of 550 eV was employed throughout the calculations to ensure well-converged total energies and electronic properties. Brillouin-zone integrations were carried out using an $8 \times7 \times7$ $\Gamma$-centered Monkhorst-Pack $k$-point mesh [40]. To obtain the equilibrium structure, both the lattice parameters and internal atomic coordinates were fully relaxed until the Hellmann-Feynman forces on every atom were smaller than $10^{-3}$ eV/Å, while the electronic self-consistency was achieved with a total-energy convergence criterion of $10^{-6}$ eV. The optimized structure served as the starting point for all subsequent calculations. The calculated equilibrium lattice parameter was directly compared with the X-ray diffraction measurements, thereby providing an independent validation of the structural model adopted for the chemically disordered alloy. To assess the dynamical stability of the considered structure, phonon calculations were carried out within the framework of density functional perturbation theory (DFPT) [41, 42, 43, 44] as implemented in VASP. The interatomic force constants and dynamical matrix were evaluated. The calculated phonon spectrum was further used to verify the absence of unstable vibrational modes, suggesting a dynamically stable phase of this structure. The calculations aim to establish a microscopic understanding of the experimentally observed structural and electronic properties by validating the stability of the chemically disordered bcc phase of $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$ and elucidating the nature of the low-energy electronic states near the Fermi level.

Results and Discussion

Structure and Microstruture

Figure 1(a) represents the X-ray diffraction pattern of $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$ collected at room temperature. The X-ray diffraction pattern of the sample confirms formation of a single phase alloy, with peak broadening due to the presence of disorder. The Le Bail refinement of recorded XRD data performed using Fullproof Suite reveals a single phase body-centered cubic (BCC) crystal structure corresponding to $\mathrm{Im}\bar{3}\mathrm{m}$ space group. The lattice constant is estimated as $a=b=c=3.3257(2)$ Å. Theoretically, the lattice parameter can be calculated by using Vegard’s law of mixtures, which is given as $a_{\mathrm{mix}}=\sum_i x_i a_i$, where $x_i$ and $a_i$ are the atomic fraction and lattice constant of the constituent elements present in the sample [45]. The evaluated theoretical value of the lattice parameter is 3.34 Å, close to the experimentally obtained value. Figure 1 (c) shows the FESEM image as well as the EDS mapping of the sample. The EDS mapping is performed at multiple positions across the sample to ensure chemical homogenity and the obtained average elemental composition is $\mathrm{Nb}_{0.46}\mathrm{Ta}_{0.27}\mathrm{Ti}_{0.12}\mathrm{Hf}_{0.16}$. The EDS data confirms the uniform distribution of all elements at the microscale in the sample. To assess the thermodynamic stability of the alloy, various phase formation parameters are calculated and summarized in Table I.

Parameters	Criteria (single phase)	$\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$
$\Delta S_{\mathrm{mix}}$ (J/K)	—	1.24R
$\Delta H_{\mathrm{mix}}$ (kJ/mol)	$\approx0$	3.29
$\Omega$	$1.1<\Omega<229.8$	78.03
$\Delta\chi$	—	0.10
$\delta(\%)$	$<6.6$	4.31
VEC	$<6.87$ (BCC)	4.73

Photoelectron Spectroscopy

To get a deeper insight into the electronic states and chemical environment of constituent elements in $\mathrm{Nb}_{0.45}\mathrm{Ta}_{0.25}\mathrm{Ti}_{0.15}\mathrm{Hf}_{0.15}$, PES measurements were performed, using Al K$\alpha$ (1486.6 eV) for XPS and a photon energy of 40 eV at room temperature for UPS. The core level spectral regions corresponding to Nb 3$d$, Ti 2$p$, Hf 4$f$, and Ta 4$f$ are shown in Figure 2 (a-d). The core levels of Nb, Ta, and Hf exhibit spin-orbit split doublets showing asymmetric Doniach-Sunjic line shape, which shows the metallic nature of the alloy. This asymmetric Doniach-Sunjic line shape arises due to final state interaction of the core hole and the conduction electrons, causing a tail toward higher binding energies [46, 47]. However, Ti shows only a single peak at 2$p_{3/2}$, but the Doniach-Sunjic type asymmetric line shape confirms its metallic nature. The peak at 2$p_{1/2}$ is suppressed due to surface contamination or oxide layer formation. The valence band UPS spectrum of the alloy, using a photon energy of 40 eV at room temperature, is shown in Figure 2 (e). The spectral weight in the vicinity of the Fermi level ($E_F$) was determined from the normalized intensity within the binding-energy range of 0–2 eV, providing a quantitative measure of the electronic states near $E_F$ [48]. The spectrum exhibits finite intensity at $E_F$, with the corresponding spectral weight estimated to be approximately 0.71. To illustrate the overall valence-band electronic structure, the spectrum over the wider binding-energy range of 0–10 eV is also presented. The observed non-zero spectral weight at $E_F$ further confirms the metallic nature of the alloy.

Normal and Superconducting-State Properties

The dc electrical resistivity ($\rho(T)$) as a function of temperature is shown in Figure 3 (a), is measured over 2–300 K under zero applied magnetic field to study the superconducting as well as normal state properties of the alloy. The inset shows the transition to the superconducting state as resistivity begins to drop at $T_c^{\mathrm{onset}}=7.81\ \mathrm{K}$.

Image not extracted

The residual resistivity ratio ($RRR$), given by $\rho_{300}/\rho_{10}$, is estimated to be 1.29. This low value of $RRR$ indicates strong disorder-dominated scattering due to high atomic-scale disorder in Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$, comparable to other reported HEAs and MEAs [18, 49]. The characteristic quadratic temperature dependence of electrical resistivity is predicted by Fermi-liquid theory. The zero-field low temperature resistivity is well fitted using the Fermi liquid theory, which is given as:

where $\rho_{0}$ is the residual resistivity resulting from temperature-independent scattering while $A$ is the electron scattering coefficient extracted from the low-temperature resistivity region and $n = 2$ for the Fermi-liquid behavior. The obtained values of $\rho_{0}$ and $A$ are 69.016 $\mu\Omega$ cm and 0.00081279 $\mu\Omega$ cm K$^{-2}$ respectively. The normal state resistivity of the sample was also fitted using the Bloch-Gruneisen (BG) model [50, 51]. The BG

formula is given as:

where C is the material-dependent constant, n identifies the primary dominant scattering mechanism in the material, and $\theta_{D}$ is the Debye temperature. The BG model best fitted using $n=3$, showing that electron-phonon scattering dominates the normal state electrical transport on this sample, and the obtained $\theta_{D}$ is 260 K. A similar trend is observed in various HEAs and MEAs [52, 53]. Figure 4 (a) shows the temperature dependent dc electrical resistivity over magnetic field ranging from 0–6 T. The measurements were carried over the temperature range of 2 to 12 K. The upper critical field ($\mu_{0}H_{c2}(0)$) is determined from temperature-dependent electrical resistivity. Figure 4 (c) shows temperature dependence of $\mu_0H_{c2}(0)$ fitted using GL equation which is given as:

Image not extracted

where $\mu_0H_{c2}(T)$ is temperature-dependent upper critical field. The estimated upper critical field $\mu_0H_{c2}(0)$ is 7.14 T by extrapolating the GL fit. Furthermore, the Ginzburg-Landau coherence length ($\xi_{GL}(0)$) is calculated using the following relation [54, 55]:

where $\phi_0$ is the magnetic flux quantum number ($h/2e$) and equals $2.07 \times10^{-15}\ \mathrm{T\,m^2}$. By using $\mu_0H_{c2}(0)$, the value of $\xi_{GL}(0)$ comes to be 6.79 nm. The superconducting state in type II superconductors are completely destroyed above $\mu_0H_{c2}(0)$. The suppression of superconductivity in type-II materials is governed by two different mechanisms: the orbital pair-breaking effect and the spin-paramagnetic pair-breaking effect. The orbital pair breaking mechanism arises when the kinectic energy of the Cooper pairs exceed the condensation energy while spin paramagnetic mechanism is due to the Zeeman splitting of the Copper pairs in the applied magnetic field.

The orbital limiting field $\mu_0H_{c2}(0)$ for BCS superconductors can be approximated by the WHH model [12, 54, 55] and given as:

In this equation, $\alpha$ is a dirty limit constant and equals 0.693. Here, we calculated the orbital pair breaking field by using the slope $-\frac{dH_{c2}}{dT}=-1.11$ at $T=T_C$ and obtained as $H_{c2}^{\mathrm{orbital}}(0)=5.92\ \mathrm{T}$. The Pauli paramagnetic field can be calculated using $\mu_0H_{c2}^{P}(0)=1.85T_C$ and equals to 14.25 T by taking $T_C=7.70\ \mathrm{K}$ (considering $T_C^{\mathrm{mid}}$). The relative strength of orbital and spin-paramagnetic effects are quantified using the Maki parameter $\alpha_M$ and given as: [56]:

The obtained $\alpha_M$ is 0.59, indicating superconductivity is limited by orbital limiting field. To examine the contributions of orbital and pauli-paramagnetic effect, the experimental $\mu_0H_{c2}(T)$ is modelled using the WHH model

[19, 54] and is given by using following equation:

where $t=T/T_{\mathrm{C}}$ is the reduced temperature and $\psi(x)$ is the digamma function.

where $\alpha_{\mathrm{M}}$ is the Maki parameter and $\lambda_{\mathrm{so}}$ is the spin-orbit interaction parameter. The WHH fit yields $\mu_0H_{\mathrm{c2}}(0)=5.61$ T without spin–orbit effect ($\lambda_{\mathrm{so}}=0$ and $\mu_0H_{\mathrm{c2}}(0)=6.43$ T by considering spin–orbit effect ($\lambda_{\mathrm{so}}=2.4$). The upper critical field derived from WHH considering the spin–orbit effect ($\lambda_{\mathrm{so}}=2.4$) value is comparable to the value obtained using eq. (5). The magnetic susceptibility measurements were performed in field-cooled warming (FCW) and zero-field-cooled warming (ZFCW) modes over the temperature range of 1.8–12 K, as shown in Figure 3 (b). The superconducting onset temperature $T_{\mathrm{c}}^{\mathrm{onset}}$ obtained from the ZFCW measurement is found to be 7.51 K. The separation of ZFCW and FCW red magnetization curves indicates irreversibility due to strong vortex pinning. To determine the lower critical field ($\mu_0H_{\mathrm{c1}}(0)$), isothermal magnetization measurements were performed over the temperature range 2 K–5.5 K. The corresponding magnetization curves are presented in the inset of Figure 4(a). The deviation from Meissner state is taken as $\mu_0H_{\mathrm{c1}}(T)$ at different temperatures. The estimated $\mu_0H_{\mathrm{c1}}(0)$ is 0.023 T, using the Ginzburg–Landau approximation as:

Image not extracted

Also, the Ginzburg–Landau penetration depth ($\lambda_{\mathrm{GL}}(0)$) is evaluated using relationship between lower critical field $\mu_0H_{\mathrm{c1}}(0)$ and coherence length $\xi_{\mathrm{GL}}(0)$ and given as follows [55, 57]:

The calculated value of $\lambda_{\mathrm{GL}}(0)$ from above eq. 10 is 151.5 nm, using $\mu_0H_{\mathrm{c1}}(0)=0.023$ T and $\xi_{\mathrm{GL}}(0)=6.79$ nm. In order to classify this material as a type I or type II superconductor, we calculated the Ginzburg–Landau parameter ($\kappa_{\mathrm{GL}}$), which is defined as $\kappa_{\mathrm{GL}}=\frac{\lambda_{\mathrm{GL}}(0)}{\xi_{\mathrm{GL}}(0)}$. The estimated value of $\kappa_{\mathrm{GL}}$ is 22.31, which is more than $\frac{1}{\sqrt{2}}$, confirming type II superconductivity. Also, the temperature dependence of the lower critical field was analyzed using the single-root Ginzburg–Landau expression, which is given as [58]:

where, $\lambda_{\mathrm{GL}}(0)$ is the Ginzburg–Landau penetration depth at zero temperature, $\xi_{\mathrm{GL}}(0)$ is the Ginzburg–Landau coherence length, and $T_{\mathrm{C}}$ is the superconducting transition temperature. The experimental $\mu_0H_{\mathrm{c1}}(T)$ data were fitted using eq. (11), with $\lambda_{\mathrm{GL}}(0)$ as free parameter, while $\xi_{\mathrm{GL}}(0)$ and $T_{\mathrm{C}}$ were fixed to the values obtained from the upper critical field analysis and resistivity measurements, respectively, as shown in Figure 4 (d). The best fit as shown in Figure 4 (d) yields $\lambda_{\mathrm{GL}}(0)=158.49$ nm. This fitted value is in close agreement with the value obtained using eq. 10.

Furthermore, heat capacity measurements are performed to elucidate the material’s superconducting nature. The heat capacity data show a jump from the superconducting to the normal state at $T_{\mathrm{c}}=7.6$ K, consistent with magnetization and resistivity data. Figure 5 (a) shows the specific heat data plotted as $C/T$ versus $T^2$. The normal state data is fitted using the following expression:

Image not extracted

where $\gamma_{\mathrm{n}}$ is the Sommerfeld coefficient related to the electronic contribution (normal state) to specific heat, and $\beta_1$ and $\beta_2$ corresponds to the phononic and anharmonic contribution to specific heat. The $\gamma_{\mathrm{n}}$, $\beta_1$ and $\beta_2$ are 4.22 $mJmol^{-1}K^{-2}$, 0.1257 $mJmol^{-1}K^{-4}$ and 0.00362 $\mu Jmol^{-1}K^{-6}$ respectively. The obtained value of $\beta_1$ is used to calculate Debye temperature $\Theta_{\mathrm{D}}$ according to the equation:

where R = 8.314 $Jmol^{-1}K^{-1}$ is the universal gas constant and N is the total number of atoms in formula unit [59]. The Debye temperature is 248.96 K, comparable to other HEAs and MEAs [12]. The electronic density of state ($D(E_{\mathrm{F}})$) at the Fermi level is calculated by using Sommerfeld coefficient ($\gamma_{\mathrm{n}}$) as follows:

where $k_{\mathrm{B}}$ is the Boltzmann’s constant ($1.38\times10^{-23}\,JK^{-1}$) and $\gamma_{\mathrm{n}}$ is 4.22 $mJmol^{-1}K^{-2}$[18]. The computed density of state at fermi level using eq. 14 is 1.79 states/eV f.u.. To analyze the nature of the superconducting state and pairing mechanism, electronic specific heat capacity ($C_{\mathrm{el}}$) is calculated by subtracting phononic contribution from measured total specific heat capacity and given as $C_{\mathrm{el}}(T)=$ $C_{\mathrm{tot}}(T)-\beta_1T^3-\beta_2T^5$. The electronic contribution ($C_{\mathrm{el}}$) is well fitted with isotropic s-wave BCS model as shown in Figure 5(b) and given as:

where f represents Fermi-Dirac distribution function and is given in terms of quasiparticle occupation index as $f(\xi)=\left[\exp\left(\frac{E(\xi)}{k_BT}\right)+1\right]^{-1}$, quasiparticle energy $E(\xi)=\sqrt{\xi^2+\Delta(t)^2}$, $y=\xi/\Delta(0)$, $t=T/T_C$ and $\Delta(t)=\Delta(0)\tanh\left[1.82\left[1.018\left(t^{-1}-1\right)\right]^{0.51}\right]$ is the defined as superconducting gap [49]. Using eq. (15), we have estimated the superconducting gap ratio $\frac{\Delta(0)}{k_BT_C}=2.01$. The obtained value exceeds the weak coupling BCS limit of 1.76, indicating the strong coupling effects in the present alloy. The obtained value is larger than the weak-coupling BCS value of 1.76, indicating the presence of strong-coupling effects. To further examine the thermodynamic consistency of the superconducting state, the strong-coupling formalism is employed [60]. The logarithmic average phonon frequency, $\omega_{\ln}$, is estimated using the experimentally determined normalized specific-heat jump, $\Delta C_{\mathrm{el}}/\gamma_nT_C=1.972$, through the relation

which yields $\omega_{\ln}=115$ K using $T_C=7.70$ K. Using $\omega_{\ln}=115$ K and $T_C=7.70$ K, $\frac{\Delta(0)}{k_BT_C}$ is further evaluated as:

resulting in $\frac{\Delta(0)}{k_BT_C}=1.96$. Further, the coupling strength of elctron-phonon interaction ($\lambda_{\mathrm{el-ph}}$) was estimated using the Allen–Dynes equation [60, 61] and given as:

where $T_C=7.70$ K, $\omega_{\ln}$ is 115 K, and $\mu^*$ is the Coulomb pseudopotential and taken as 0.13 [61]. The calculated value of the electron–phonon coupling constant $\lambda_{\mathrm{el-ph}}$, using eq. 18 is 1.07, indicating that the present system lies in the strong-coupling regime.

Critical Current Density and Pinning Mechanism

Figure 6 (a) shows the magnetic hysteresis loop (M-H) of Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$ at 2K. The magnetic hysteresis curve exhibits type II superconductivity. The $J_c$ is a key parameter for potential applications of superconductors and is calculated from the magnetization curve, using Bean’s model [62] and given as:

where $\Delta M$ is the vertical magnetization width and $a$ and $b$ are the sample dimensions such that $a<b$. The sample dimensions are 1 mm and 2 mm, respectively, in our case. The variation of $J_c$ with field is shown in Figure 6 (b). The obtained value of $J_c$ is $0.48 \times10^{6} A/cm^{2}$. This value of $J_c$ exceeds the benchmark of $10^{5} A/cm^{2}$, required for practical applications. Remarkably, the $J_c$ value surpasses the value obtained for other previously reported as-cast HEAs and MEAs. For instance, Nb$_{0.4}$Hf$_{0.2}$Zr$_{0.2}$Ti$_{0.2}$ shows $J_c$ of $1 \times10^{4} A/cm^{2}$ [63]. In 3d-rich Sc$_{0.2}$Ti$_{0.2}$V$_{0.2}$Nb$_{0.2}$Cu$_{0.2}$ shows $J_c = 0.18 \times10^{6} A/cm^{2}$ at 2 K [18]. In equiatomic Nb$_{0.25}$Ta$_{0.25}$Ti$_{0.25}$Zr$_{0.25}$ critical current density is $\approx0.13 \times10^{6} A/cm^{2}$ at temperatures ranging from 2 K to 4 K [19]. The eutectic NbScTiZr exhibits $J_c$ exceeding $10^{6} A/cm^{2}$ while, the BCC Nb$_{2/6}$Ta$_{1/6}$Ti$_{1/6}$Zr$_{1/6}$Hf$_{1/6}$ and (NbTa$_{0.7}$)(TiZrHf$_{0.5}$) alloys exhibit $J_c$ of the order of $10^{5} A/cm^{2}$ at 2 K [20]. Remarkably, the present material highlights the highest critical current density value observed in any reported pure-phase as-prepared HEAs and MEAs.

Image not extracted

In order to elucidate the flux pinning mechanism, the normalized pinning force $f_p = F_p/F_{p\max}$ is plotted as a function of the reduced magnetic field $h = H/H_{\mathrm{irr}}(T)$ and shown in the Figure 6 (c). Using $J_c$ value, the pinning force $F_p$ is given as $F_p \propto\mu_0 HJ_c$. The scaling behavior of $f_p$ is described by using the Dew-Hughes model as $f_p \propto h^p(1-h)^q$ [64], where p and q are fitting parameters to provide information about the flux pinning mechanism. The experimental data is fitted with the characteristic exponents $p = 0.5$, $q = 2$ for surface pinning and $p = 1$, $q = 2$ for point pinning. The fitting results indicate that the flux pinning behavior in the present alloy is mainly governed by the point pinning contribution. The severe lattice distortion discussed in the latter part of the manuscript could enhance the pinning effect at grain boundaries and therefore cause high $J_c$ in this sample.

Strongly Correlated Behaviour

The Kadowaki–Woods ratio (KWR) is calculated to investigate the correlation between the electron–electron scattering coefficient (A) and quasiparticle mass renormalization factor in terms of $\gamma_n$, thereby probing electronic correlations. Figure 7(a) represents the KWR plot of $A$ vs. $\gamma_n^2$ for various materials, where $A$ is the coefficient of $T^2$ obtained from the low-temperature resistivity ($\rho= \rho_0 + AT^2$) and $\gamma_n$ is obtained from the electronic specific heat. The KWR is defined as $A/\gamma_n^2$. For transition metals the characteristics value of KWR is $a_{\mathrm{TM}} = 0.4 \,\mu\Omega$ cm mol$^2$ K$^2$ J$^{-2}$, whereas in heavy fermion system it is significantly enhanced and reaching $a_{\mathrm{HF}} = 10 \,\mu\Omega$ cm mol$^2$ K$^2$ J$^{-2}$ [65, 66]. For Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$, the KWR is $45.48\mu\Omega$ cm mol$^2$ K$^2$ J$^{-2}$, placing the sample near heavy-fermion line [21, 57]. This large ratio of $A/\gamma_n^2$ indicates pronounced electron-electron correlations and classifies Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$ as a strongly correlated material. The large KWR value is likely to arise from enhanced disorder, and electron-phonon interactions [57]. Strongly correlated behavior is uncommon in alloy superconductors, and most of the time, they highlight anomalous superconducting behavior [69]. Therefore, in the current MEA, the emergence of strongly correlated behavior opens a promising opportunity to explore the interplay among superconductivity, electronic correlations, and extreme disorder [21, 57].

Image not extracted

Figure 7 (b) shows the variation of transition temperature $T_C$ with valence electron Count (VEC) for various crystalline metal alloys, amorphous alloys, and BCC HEAs. According to the Matthias rule, BCC HEAs exhibit the highest $T_C$ within the VEC range of 4.67–4.70. The present sample is consistent with the VEC range reported for other HEAs [67, 68].

Lattice Distortion

HEAs and MEAs exhibit severe lattice distortion due to the presence of multicomponent elements having different sizes [72, 73]. In order to quantify the local lattice distortion in BCC Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$, interatomic distances are calculated using 60 atoms SQS as shown in Figure 8 (c). The broadness of the interatomic distances plot indicates that atoms deviate significantly from their mean positions, causing pronounced lattice distortion. This peak broadness is attributed to size mismatch of constituent elements as $r_{\mathrm{Ti}} \leq r_{\mathrm{Nb}} \leq r_{\mathrm{Ta}} \leq r_{\mathrm{Hf}}$. The average interatomic distance for the first five NN pairs is 2.89, 3.32, 4.72, 5.50, and 5.73 Å, respectively. However, the peak position is still centered around their ideal interatomic distances, predicting crystalline behaviour is preserved over long-range. Figure 8 (a) shows a 60-atom supercell of the BCC structure using SQS. The atomic size mismatch is commonly used to quantify the effects of local lattice distortion (LLD) in HEAs and MEAs, using the equation:

Image not extracted

where $\bar{r}$ being the average atomic radius of the constituent elements and given as $\sum_{i=1}^{n} c_i r_i$, where $c_i$ and $r_i$ are atomic fraction and atomic radius of individual elements, respectively. This approach yields somewhat ambiguous results, as atomic radii depend on factors such as the nature of bonding, crystal structure, and external conditions [70]. Hence, for accurate measurement of local lattice distortion with atomic radii measured under different conditions, we measured the mean atomic displacement as follows:

where $(x_i,y_i,z_i)$ and $(x_i',y_i',z_i')$ are the reduced coordinates of the $i$ th atom before and after the structural relaxation. The obtained value of mean atomic displacement ($\Delta d$) is 0.086 Å, which indicates the displacement of atoms from their ideal positions. This result predicts that $\delta$, as a radius-dependent parameter, can overestimate the lattice distortion, while $\Delta d$ gives a more precise intrinsic structural response to chemical disorder. Figure 8 (b) shows a linear dependence of the $\delta$ and $\Delta d$. Typically, a larger atomic size mismatch ($\delta$) is associated with a greater local lattice distortion ($\Delta d$). Nevertheless, in some 3d HEAs and MEAs, $\delta$ tends to slightly overestimate the local lattice distortion (LLD) when it is calculated using the atomic radii $r_i$ of the individual elements [71].

To get deeper insight of microscopic origin of LLD, a pair-wise analysis of fist nearest neighbour (1NN) shell was carried out and presented in Figure 8 (d). The BCC alloy Nb$_{0.45}$Ti$_{0.15}$Ta$_{0.25}$Hf$_{0.15}$ composed of both homoatomic as well as heteroatomic pairs of atoms. The homoatomic pairs include Nb-Nb, Ta-Ta, Ti-Ti, Hf-Hf while heteroatomic pairs composed of Nb-Ti, Nb-Ta, Nb-Hf, Ti-Ta, Ti-Hf and Ta-Hf. The computed results yield comparable bond lengths for MEAs and HEAs. The mean and median corresponding to each pair are shown in Figure 8 (d). The pair-wise LLD can be evaluated using the formula:

where $d_i$ corresponds to the interatomic distance of the $i$ th atom in relaxed SQS supercell and $d_r$ is the reference distance in unrelaxed structure. Using eq. 22, the maximum distortion is shown by the Hf-Hf pair (3.45%) and the Ti-Ti pair (0.85%), respectively, reflecting the larger atomic radius and correspondingly stronger local lattice strain of Hf compared with Ti in the BCC environment. The higher value for Hf-Hf arises from the pronounced size-mismatch, driven expansion when two Hf atoms occupy nearest-neighbor sites, while the smaller Ti–Ti distortion indicates that Ti induces only a moderate local strain due to its relatively smaller size mismatch with the BCC lattice. A similar trend was observed in $(TaNb)_{0.67}(HfZrTi)_{0.33}$ [74]. Among heteroatomic pairs, the largest distortion is exhibited by Hf-Ti and Hf-Ta, which are respectively 3.13% and 2.40%. The atomic pair is compressed when the median bond length is shorter than corresponding reference bond length of the pair e.g., Nb-Nb pair. The heteroatomic pairs such as Nb-Ti and Nb-Ta shows less deviations from their neighbouring atoms, indicating a balanced atomic environment. The heteroatomic pairs, such as Nb with Ti and Ta, appear to be balanced with their neighboring atoms. The elongation or compression of bond lengths is governed by the local environment of a atomic pair, where certain neighboring atoms induce compression while others cause elongation of the bonds. Also, changes in bond length are related to chemical interactions among multiprincipal elements in MEAs. Further, the shortest and longest first nearest neighbor bond lengths for Nb-Nb and Ta-Ta pairs in SQS supercell is shown in Figure 8 (e-f).

First-Principles Studies of Structural Stability and Electronic Properties

To gain microscopic insight into the experimentally observed structural stability and metallic normal-state behavior of Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$, we performed first-principles calculations within the generalized gradient approximation (GGA). Since the alloy possesses intrinsic substitutional disorder arising from the random occupation of lattice sites by multiple constituent elements, the chemical disorder was explicitly modeled using the special quasirandom structure (SQS) approach. The calculated lattice dynamical and electronic properties are summarized in Figure 9. The calculated phonon dispersion Figure 9(a) exhibits no imaginary phonon frequencies throughout the Brillouin zone, confirming the dynamical stability of the chemically disordered body-centered cubic phase. The optimized equilibrium lattice parameter is also in excellent agreement with the experimentally determined value, demonstrating that the SQS model faithfully captures the disordered crystal structure realized experimentally. The phonon spectrum extends up to approximately 8 THz, with acoustic modes below ∼2 THz followed by relatively flat optical branches at higher frequencies. Such broad phonon features originate from the

Image not extracted

Properties	unit	Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$
$T_{\mathrm{C}}^{\mathrm{onset}}$ ($\rho$ - $T$)	K	7.8 (1)
$T_{\mathrm{C}}^{\mathrm{mid}}$ ($\rho$ - $T$)	K	7.7 (1)
$\mu_0H_{\mathrm{c1}}(0)$	T	0.0232 (1)
$\mu_0H_{\mathrm{c2}}(0)$	T	7.14 (9)
$\mu_0H_{\mathrm{c2}}(0)$ (WHH)	T	7.14 (9)
$\mu_0H_{\mathrm{c2}}^{\mathrm{P}}(0)$	T	14.33
$\xi_{\mathrm{GL}}(0)$	Å	67.9 (4)
$\lambda_{\mathrm{GL}}(0)$	Å	1515 (3)
$\kappa_{\mathrm{GL}}(0)$		22.31
$\gamma$	mJ mol$^{-1}$ K$^{-2}$	4.22(6)
$\beta_1$	mJ mol$^{-1}$ K$^{-4}$	0.1257(6)
$\beta_2$	mJ mol$^{-1}$ K$^{-4}$	0.0036(62)
$J_{\mathrm{C}}$	A/cm$^2$	0.48 × 10$^6$
$\theta_{\mathrm{D}}$	K	248.96 (4)
$KWR$	μΩ cm mol$^2$ K$^2$ J$^{-2}$	45.48(12)
$D(E_{\mathrm{F}})$	states/eV f.u	1.79 (3)
$\omega_{\mathrm{ln}}$	K	115 (6)
$\lambda_{\mathrm{el-ph}}$		1.07 (3)
$\Delta_0/k_{\mathrm{B}}T_{\mathrm{C}}$		2.01
$\Delta C_{\mathrm{el}}/\gamma_{\mathrm{n}}T_{\mathrm{C}}$		1.98

coexistence of constituent atoms with different masses and local bonding environments, reflecting the intrinsic chemical disorder characteristic of entropy-stabilized alloys. These results establish that the chemically disordered BCC phase is both structurally and dynamically stable, consistent with the experimentally observed single-phase crystal structure.

ordered BCC phase is both structurally and dynamically stable, consistent with the experimentally observed single-phase crystal structure.

The corresponding electronic structure is shown in Figure 9(b,Figure 9c). The total and atom-resolved density of states reveal a finite electronic density at the Fermi level, confirming the metallic ground state of Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$. This is fully consistent with the experimentally observed metallic resistivity above the superconducting transition. The electronic states near $E_{\mathrm{F}}$ are dominated by transition-metal $d$ orbitals, with Ti and Ta contributing the largest spectral weight, while Nb and Hf provide broader contributions over a wider energy range through strong orbital hybridization. The absence of a pseudogap or electronic instability near $E_{\mathrm{F}}$ indicates that the pronounced chemical disorder does not disrupt the itinerant electronic states required for superconductivity. The calculated band structure further elucidates the low-energy electronic states responsible for the metallic behavior. Multiple dispersive bands cross the Fermi level along the high-symmetry directions, giving rise to several electron and hole pockets and indicating an intrinsically multiband electronic structure. Rather than arising from an individual atomic species, these low-energy bands originate from the hybridized $d$-orbitals of Nb, Ta, Ti, and Hf, reflecting the cooperative electronic nature of the chemically disordered alloy. Such multiband electronic states naturally account for the finite density of states at the Fermi level and provide the microscopic electronic framework for the experimentally observed metallic transport and superconductivity.

Thus, our first-principles calculations establish a coherent microscopic picture of Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$. The absence of unstable phonon modes confirms the stability of the chemically disordered BCC phase, while the multiband $d$-electron metallic state explains the robust normal-state conductivity that precedes superconductivity. Together with the experimentally observed strong electron-phonon coupling and superconducting transition, these calculations demonstrate that superconductivity in Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$ emerges from a dynamically stable, chemically disordered metallic phase whose low-energy electronic structure is governed by strongly hybridized transition-metal $d$ states.

Conclusion

In summary, we have carried out a comprehensive experimental and first-principles investigation of the structural, superconducting, electronic, and local lattice-distortion properties of the mechanically robust medium-entropy alloy Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$. Structural characterization confirms the formation of a homogeneous single-phase body-centered cubic structure with a nearly stoichiometric elemental distribution. Low-temperature magnetization, transport, and specific-heat measurements establish bulk type-II superconductivity with a transition temperature of $T_{\mathrm{c}}^{\mathrm{onset}} \approx7.8$ K, an upper critical field of $\mu_0 H_{\mathrm{c2}}(0) = 7.14 \pm0.08$ T, and a strongly coupled isotropic $s$-wave superconducting state. Remarkably, Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$ exhibits a critical current density of $0.48 \times10^6$ A/cm$^2$, the highest reported among as-cast medium- and high-entropy superconducting alloys, highlighting its exceptional capability for sustaining large supercurrents. Photoemission spectroscopy confirms the metallic electronic character with finite spectral weight at the Fermi level, while GGA-based first-principles calculations establish the dynamical stability of the chemically disordered BCC phase and reveal a multiband metallic electronic structure dominated by hybridized Ti and Ta $d$ states near the Fermi level. The SQS calculations further quantify substantial disorder-induced local lattice distortion ($\Delta d = 0.086$ Å), providing a microscopic picture of the intrinsic structural disorder that is likely responsible for the enhanced vortex-pinning landscape and the observed high critical current density.

Beyond its outstanding superconducting performance, the unusually large Kadowaki-Woods ratio places Nb$_{0.45}$Ta$_{0.25}$Ti$_{0.15}$Hf$_{0.15}$ close to the heavy-fermion regime, suggesting that significant electronic correlations coexist with strong electron-phonon coupling despite the presence of extreme chemical disorder. The concurrent presence of high critical current density, strong vortex pinning, enhanced electronic correlations, and exceptional mechanical robustness positions this material as a highly promising candidate for next-generation superconductors for operation under extreme conditions.

Acknowledgments

S. M. acknowledges financial support from the ANRF, Government of India, for the Advanced Research Grant (ANRF/ARG/2025/004870/PS). S. M. also acknowledges the financial support from the SPARC, Ministry of Education, Govt. of India (Project no. P4139). S.K.P. acknowledges support from ANRF (previously SERB), Government of India, for the core research grant (CRG/2023/003063)