1. Introduction
The rapid progress of industrialization worldwide has led to the extensive utilization of hard materials in various fields, including weaponry, aviation, machinery, construction, geological exploration, and mining [
1,
2]. Traditionally, hard materials are composed of compounds featuring shorter covalent bonds centered around lighter elements with smaller atomic radii, such as boron, carbon, nitrogen, and oxygen. However, due to the limitations of traditional hard materials in terms of mechanical and electrical properties, there is a growing interest in exploring novel hard materials comprising compounds formed from light elements and transition metals [
3,
4]. Due to the high density and large bulk modulus of transition metals, the doping of light elements can be carried out with the
d-orbital of transition metals to form quasi-covalent bonds [
5], significantly enhancing the hardness. Moreover, the distribution of the
d-orbitals of transition metals in the outer space of the nucleus is more complex compared to the traditional superhard materials. Therefore, these materials frequently form diverse structures, exhibiting exceptional properties and potential benefits for substituting conventional hard materials in domains such as electronics, information technology, and aerospace [
6].
Significant advancements have been achieved in recent years regarding the investigation of transition metal borides. In addition to high hardness, a high melting point, and strong corrosion resistance similar to traditional superhard materials, certain compounds also exhibit properties such as superconductivity [
7], magnetism, electrical conductivity [
8], and catalysis [
9]. Transition metal borides have emerged as extensively investigated among the transition metals, owing to their finely dispersed powder form and their proclivity to readily react with metals under relatively moderate conditions. This reactivity leads to the formation of a diverse range of stable compounds. ZrB
12, as a noteworthy representative of this category, has garnered significant attention due to its remarkable physicochemical properties. The experiments and first-principles calculations on ZrB
12 show that the symmetrical B-B covalent network can form delocalized
π-bonds, and the interaction of the
4d-orbitals of Zr atoms with these
π-bonds creates electronic channels, thereby maintaining the material’s ultrahardness while also providing good electrical conductivity [
10,
11,
12]. Notably, heavily doped ZrB
12 has a calculated hardness of up to 40 GPa [
13,
14], while its superconducting transition temperature approaches 6 K [
15]. Singularly, ZrB
12 exhibits a range of behaviors that deviate significantly from the typical characteristics observed in conventional superconductors. Several theoretical models have been proposed to elucidate its superconducting properties, each offering a distinct perspective on the nature of this intriguing material. One such model characterizes ZrB
12 as a traditional superconductor with amplified surface effects that could be responsible for its unique behavior [
16,
17]. Another model suggests that ZrB
12 possesses a strongly coupled
s-wave gap structure [
18], while another intriguing proposal posits that it may exhibit
d-wave superconductivity. The temperature-dependent behavior of the magnetic penetration depth in ZrB
12 reveals a pronounced two-band structure feature [
19], which implies that the interaction between these bands is relatively weak [
20]. Recent band structure calculations have determined that the Fermi surface of ZrB
12 consists of both an open sheet and a closed sheet [
21], which is a characteristic that could underpin the observed multiband nature of its superconductivity. Additionally, specific heat and vortices measurements consistent with both type 1 and type 2 superconducting behaviors [
22,
23]. Zhang et al. provided further insight by confirming, through the use of X-ray diffraction techniques, that ZrB
12 belongs to the category of type 1.5 superconductivity [
24,
25]. These results suggest that ZrB
12 may represent a distinct class of superconductors, inherently capable of simultaneously exhibiting both type 1 and type 2 characteristics simultaneously. Therefore, the contribution of electron–electron interactions in ZrB
12 is expected to be substantial. However, recent research on ZrB
12 has predominantly focused on the superconducting type, with limited attention given to investigating the electron–electron interactions in the superconducting character. This intriguing possibility has piqued our interest and motivated us to further explore the superconducting properties of ZrB
12, particularly through the application of high-pressure techniques.
Applying pressure can often effectively change the electronic structure of compounds, and even induce phase transitions, without introducing additional impurities [
26,
27]. In this study, we conducted resistance characterization under pressure to investigate electron–electron interactions in the superconductor ZrB
12. Encouragingly, ZrB
12 exhibits exceptional stability with minimal structural impact even at pressures as high as 44 GPa, thereby preserving its phase integrity [
28,
29]. Our research reveals that the superconductivity of ZrB
12 exhibits anisotropic characteristics, which persist even under high-pressure conditions. Furthermore, under a pressure of 1.5 GPa, ZrB
12 still exhibits superconductivity at approximately 5.7 K. With increasing pressure, the superconducting transition temperature decreases. ZrB
12 exhibits metallic behavior under pressures up to 31.4 GPa. The diminishing superconductivity observed in ZrB
12 can be attributed to the intensified effect of Zr’s movement on phonon dispersion, ultimately causing a reduction in carrier concentration.
3. Results
ZrB
12 crystallizes in a cubic crystal system characterized by the space group
Fm-3m, which is assigned the number 225 in the International Tables for Crystallography. Within this symmetrical framework, the Zr atoms are located at specific Wyckoff positions, which are designated sites within the unit cell that can be occupied by atoms or ions. The Zr atoms are surrounded by clusters of boron atoms, forming B
12 icosahedra. These B
12 clusters form a rigid, three-dimensional network that contributes to the exceptional properties of ZrB
12, including a high melting point and excellent thermal conductivity. The cubic arrangement with B
12 clusters renders ZrB
12 a material of significant interest for various high-temperature and high-strength applications [
30]. The Zr atoms in ZrB
12 have an electron configuration of
4s24p64d35s1 and the B atoms have
2s22p1, which serve as valence electrons. These contribute to covalent bonding, thereby fulfilling the characteristics of ultra-hard materials, including a high average valence electron count per atom, a maximal number of covalent bonds per unit volume, and a high bond energy. In order to investigate the electrical transport properties under external pressure on ZrB
12, we measured temperature dependence resistance from 2 to 300 K for both [110] and [100] orientations of a ZrB
12 single crystal, as shown in
Figure 1a. Here, the
Tconset is determined as the interception between two straight lines below and above the superconducting transition, while
Tczero is defined as the zero-resistivity temperature. As the temperature decreases, the
Tconset for the [100] and [110] is 5.8 K and 5.7 K, whereas the
Tczero is 4.6 K and 3.4 K, respectively [
31]. Moreover, the resistivity has a certain difference, and the resistivity in the [100] orientation is significantly greater than that of the [110] orientation. These findings demonstrate the anisotropic electrical behavior of ZrB
12, which is consistent with prior research reports [
13]. The specific heat measurements of the [100] orientation below 10 K reveal a peak heat capacity within the weak magnetic field near 0 T, as shown in
Figure 1b. An exothermic process around 5.4 K was confirmed, providing further evidence for a superconducting transition around this temperature [
15,
32]. In ZrB
12, which is classified as a two-gap superconductor exhibiting dynamic charge stripes, the interplay between these fluctuating charges and magnetic vortices can occur even under relatively low magnetic fields. This interaction induces a pronounced anisotropy in the superconducting properties of the material. Specifically, the behavior of the superconducting characteristics, such as the critical temperature and the critical current density, becomes highly dependent on the orientation of the applied magnetic field. This anisotropy is a direct consequence of the dynamic charge stripe order within the superconductor, which modulates the electronic structure and thereby influences the pinning and motion of vortices within the material. Bolotina et al. show a stepped singularity of the specific heat capacity of the [100] and [110] crystal orientations at around 5.8 K, and another hump of the [110] orientation appeared at around 5.7 K, demonstrating the band-gap anisotropy of ZrB
12 single crystals [
33]. Generally, the resistivity of ZrB
12 will immediately approach 0 upon reaching the first superconducting temperature. Therefore, in our resistivity experiment, we can only observe the tendency of the resistivity drop to zero and cannot observe the second peak of the [110] orientation. The presence of such a pronounced anisotropy in the superconducting behavior of ZrB
12 highlights the intricate relationship between the charge stripe fluctuations and the vortex dynamics, and it underscores the complexity of the superconducting state in materials with this unique electronic structure [
33]. Consequently, the change in the heat capacity under the magnetic field, the magnitude of the resistivity, and the differences in the superconducting temperature collectively confirm the anisotropy of ZrB
12.
Furthermore, we characterized the temperature dependence of resistance for a ZrB
12 single crystal with the [100] orientation at 1.5 GPa under an external magnetic field ranging from 0 to 1 T (as shown in
Figure 2). As shown in
Figure 2a, the
Tconset at 0.05, 0.1, 0.15, 0.2, 0.25, and 0.3 T is 5.4, 5.0, 4.7, 4.3, 3.9, and 3.5 K, while the corresponding
Tczero values are 3.2, 2.8, 2.1 K, respectively.
TC decreases with an increasing applied magnetic field H, and this variation is illustrated in
Figure 2b. Using the Werthamer–Helfand–Hohenberg (WHH) formula,
where
Tc0 is the superconducting transition temperature under zero magnetic field. At 1.5 GPa, the upper critical field
Hc2(0) was estimated to be 0.64 T. From
Figure 1b and
Figure 2a, it is apparent that the intensification of the magnetic field to 1 Tesla completely suppresses superconductivity, indicating that the [100] orientation of the ZrB
12 single crystal belongs to the classification of type 2 superconductivity. Based on the Ginzburg–Landau theory, Abrikosov categorized superconductors into two distinct types: type 1 and type 2, according to the magnitude of the Ginzburg–Landau parameter. ZrB
12 has been proposed to exhibit a two band-effect and related type 1.5 superconductivity (containing type 1 and type 2) by using X-ray diffraction and
μSR studies [
26,
34]. The observation of type 2 superconductivity presented in
Figure 2b does not contradict the conclusion of type 1.5 superconductivity, as the type 1.5 superconductivity of ZrB
12 can only be observed below
Tc, whereas our discussion primarily focuses on the relationship between
Tc and the magnetic field H. We take into account the anisotropy along the [100] direction, which restricts our observations solely to type 2 superconductivity. Furthermore, as illustrated in
Figure 2a, at relatively low magnetic fields and below
Tc, the superconducting transition exhibits a remarkably sharp profile. As the magnetic field intensifies, the transition gradually broadens, which may be attributed to various vortex phase transitions [
26].
Figure 2b shows the characteristic of the upper critical field
Hc2(0) in type 2 superconductors, which provides insight into the microscopic mechanism of ZrB
12. Generally, two independent mechanisms govern the suppression of superconductivity by magnetic fields. One involves the orbital pair breaking of Cooper pairs in the superconducting state, which is associated with screening currents that exclude the external field which is commonly referred to as the orbital limit. The other mechanism is a spin effect resulting from Zeeman splitting, which is applicable only to singlet pairings and is widely known as the Pauli paramagnetic limit. Evidently, ZrB
12 aligns with the former theoretical explanation.
In order to comprehensively investigate the electrical transport properties within a high-pressure environment, it is imperative to ensure and validate the structural stability of the ZrB
12 compound under compressive forces, without experiencing any detrimental alterations or deformations. This validation is crucial because the stability of ZrB
12’s crystal lattice directly impacts its ability to effectively conduct electricity under such extreme conditions, thereby influencing the reliability and accuracy of measurements pertaining to electrical transport properties. A high-pressure synchrotron X-ray diffraction (λ = 0.6199 Å) study conducted in a DAC was used to investigate the phase stability of ZrB
12.
Figure 3a presents the XRD patterns obtained during ambient-temperature compression, revealing that the ZrB
12 single crystal sustains its structural stability up to 43 GPa, without any observable phase transition. The lattice parameters of ZrB
12 in
Figure 3b exhibit a nearly linear decrease with increasing pressure. In
Figure S3, the pressure–volume data for ZrB
12 demonstrate a gradual reduction in both the lattice parameters and volume as pressure increases. This can be attributed to the compression of atomic distances under elevated pressures, resulting in enhanced repulsive interactions among atoms and consequently higher crystal incompressibility. The Zr atomic radius is significantly larger than that of B, the interplanar spacing of ZrB
12 gradually decreases under pressure, and the influence on the orbital of Zr is much greater than that of B. Therefore, on the premise of ensuring the structural stability of ZrB
12 under pressure, our research on the high-pressure behavior of ZrB
12 primarily focuses on investigating alterations in the Zr-
d atomic orbital. The third-order Birch–Murnaghan equation of state (EOS) was employed to fit the ZrB
12 unit cell, resulting in a bulk modulus B = 221.1 GPa, as illustrated in
Figure S3, which demonstrates its compressibility. Therefore, it can be concluded that ZrB
12 maintains its structural stability under high pressure. Subsequently, we investigated the electrical transport characteristics of ZrB
12 single crystals under pressures ranging from 4.4 to 44 GPa, as depicted in
Figure 4. Specifically,
Figure 4a exhibits the [100] orientation resistance–temperature curves from 2 to 300 K across 4.4 to 20 GPa, with an inset figure presenting these curves normalized between 2 and 10 K. The obtained data clearly demonstrate that as pressure increases, ZrB
12 consistently maintains its superconducting state, while exhibiting a decrease in
Tc. Furthermore, under pressures ranging from 31.4 to 44 GPa, ZrB
12 exhibits metallic behavior, as illustrated in
Figure 4b.
Figure 5 presents the phase diagram of the superconducting transition temperature (SC) dependence pressure of ZrB
12, illustrating a region of gradual decrease in
Tconset and
Tczero from 0 to 44 GPa. The relationship between the superconducting transition temperature and pressure exhibits a gradual decline with a slope of approximately 0.236 K/GPa. From the phase diagram, we can easily visualize the evolution and intimated correlations between the metal and SC as a function of pressure. As pressure increases gradually,
Tconset is suppressed, accompanied by an initial enhancement of
Tc with a broad superconducting transition width, as displayed in
Figure 4a. Interestingly, our observations under applied pressure indicate that despite the absence of structural phase transition in ZrB
12 under elevated pressures, the superconducting transition temperature disappears at 20 GPa, coinciding with an increase in pressure. This phenomenon exhibits a metallic behavior, which is strongly associated with the electron–electron interaction within the single crystal.
4. Discussion
The exceptional hardness exhibited by ZrB
12 can be largely attributed to the presence of strong covalent bonds between the boron atoms within its crystalline structure. These B-B covalent bonds are a direct consequence of the unique integration and interaction of boron elements within the framework of transition metals, specifically Zr. The robustness of these bonds significantly contributes to the overall hardness of the ZrB
12 crystal, rendering it a material of interest in various high-performance applications where resistance to wear and deformation is crucial. The existence and positioning of Zr atoms play a crucial and indispensable role in facilitating and enhancing the electrical conductivity and superconducting capabilities of this specific material. Furthermore, under conditions of zero pressure, it is observed that there is a complete absence of energetic interaction or overlap between the optical phonon modes and the acoustic phonon modes [
35]. This intriguing phenomenon can be attributed to the distinct and varying atomic masses of Zr atoms and B atoms as they are arranged within the crystal lattice unit cell. The disparity in their masses ensures that the vibrational modes associated with these atoms do not intersect energetically, thereby maintaining a clear separation between the two types of phonon modes. This separation is essential for the material to exhibit its unique electrical properties, particularly its ability to conduct electricity and transition into a superconducting state under specific conditions [
35]. As pressure rises, ZrB
12 exhibits remarkable resistance to compression, which is attributed to the robust B-B covalent bond. The band structure of ZrB
12 remains largely unaltered, with the primary variations originating from the Zr atomic orbital. Consequently, our focus lies in exploring the disappearance of superconductivity resulting from the electronic–electron interactions in ZrB
12 under pressure, along with the corresponding changes in the Zr atomic orbital. The electronic density of states (DOS) and the partial density of states (PDOS) of ZrB
12 at 0 GP demonstrated that the DOS at the Fermi level is predominantly dominated by Zr-
d states. This indicates that the electron carriers in ZrB
12 originate primarily from the Zr-
d states. In addition, the changes in cell parameters and volume with pressure derived from our experimental values (as shown in
Figure 3 and
Figure S3) are basically consistent with the theoretical calculated values of Li, X. et al. [
35]. According to the calculation results of ZrB
12 by Li, X. et al., the DOS of ZrB
12 at 0, 20, 50 GPa was 1.4679 eV, 1.4001 eV, and 1.3571 eV, respectively. When the pressure is below 50 GPa, the carrier concentration gradually decreases under the action of pressure and reaches the minimum value at 50 GPa, indicating that ZrB
12 exhibits limited compressibility due to its high bulk modulus. In general, in the case of internal stability of ZrB
12 the change in carrier concentration and resistance are correlated. As can be seen from
Figure 4a, the resistance of our sample gradually decreases under the applied pressure. However, from
Figure 4b, when the pressure is further increased to 44 GPa, ZrB
12 exhibits significant resistance to compression, with a resistance level similar to that observed at 31.4 GPa, indicating that the change in carrier concentration during this stage is relatively small and gradually reaches a minimum value. These results are consistent with the calculated value. As the pressure increases, the contribution of Zr movement to phonon dispersion becomes more significant [
35]. With the growing influence of Zr’s movement on phonon dispersion, it may result in a reduction in carrier concentration. Therefore, within the pressure range of 0 to 44 GPa, the carrier concentration of ZrB
12 steadily decreases as pressure rises. This reduction in carrier concentration leads to a decrease in its superconducting transition temperature. Lower carrier concentrations require even lower temperatures for achieving superconducting transition. Notably, at 20 GPa, the superconducting transition was undetectable, suggesting orbital pair breaking of Cooper pairs in the superconducting state of ZrB
12, thereby indicating a transition to a metallic state.