Open AccessArticle Unveiling Superior Fracture Toughness in MnCoSb Half-Heusler Alloy: A First-Principles Guide for Designing Damage-Tolerant Functional Materials 1 College of Electronic and Information Engineering, Anshun Uniνersity, Anshun 561000, China 2 School of Physics and Electronic Science, Zunyi Normal University, Zunyi 563000, China 3 Institute of Advanced Optoelectronic Materials and Technology, College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, China * Author to whom correspondence should be addressed. Molecules 2026, 31(12), 1994; https://doi.org/10.3390/molecules31121994 (registering DOI) Submission received: 7 May 2026 / Revised: 27 May 2026 / Accepted: 4 June 2026 / Published: 7 June 2026 Abstract In this study, the stability, electronic, structural, and fracture toughness, and mechanical properties of the Half-Heusler(HH) alloys MnCoSb, MnCoAs, MnCoP, and MnNiSb were comprehensively investigated using first-principles calculations based on density functional theory (DFT). The calculated results reveal that all four alloys exhibit half-metallic characteristics, characterized by the presence of a substantial band gap in the spin-down channel. The phonon spectra and negative formation energies confirm that these alloys possess both dynamic and thermodynamic stability. The Born criteria further validate the structural stability in terms of mechanical properties. Three-dimensional representations of the Young’s modulus, bulk modulus, and shear modulus for the four alloys indicate that MnCoP exhibits the most pronounced anisotropy. The overall fracture toughness of the alloys ranges from 1.58 MPa·m 1/2 to 2.63 MPa·m 1/2, which falls within the typical range for half-metallic materials, albeit at the lower end, attributable to the relatively ductile nature of the four alloys. Although the two methods yield different absolute values, the explicit crack model (Method I) is considered more reliable for anisotropic systems because it directly simulates crack propagation and accounts for local relaxations, while the empirical formula (Method II) provides a useful reference for high-throughput screening. Among the alloys, MnCoSb demonstrates a superior mechanical performance, with KIC values of 2.63 MPa·m 1/2 and 1.58 MPa·m 1/2 and brittleness indices M of 8.97 and 14.94, indicating excellent damage tolerance compared to the other three alloys. In contrast, MnCoP exhibits higher brittleness and lower mechanical reliability, with KIC values of 2.00 MPa·m 1/2 and 1.63 MPa·m 1/2 and higher M values of 13.83 and 16.99. This study provides quantitative predictions of fracture toughness and establishes a relationship between microscopic and mechanical properties. These findings offer a theoretical foundation for the application of damage-tolerant HH alloys in fields such as spintronics and magnetism. 1. Introduction Materials are increasingly integral to daily life, with diverse applications in aerospace, high-end energy manufacturing equipment, shipbuilding, automotive engineering, and transportation infrastructure. The demand for advanced materials is gradually shifting toward multifunctional integration. Heusler alloys have attracted the attention of the scientific community for decades [ 1] due to their tailorable physicochemical properties, making them a subject of extensive research. Their unique characteristics render them highly suitable for advanced applications in spintronics, thermoelectrics, and optoelectronics [ 2, 3, 4, 5, 6]. Early research on Heusler alloys has been summarized by F. Dahmane et al. [ 7]. A foundational discovery by German mining engineer Friedrich Heusler showed that specific combinations of non-ferromagnetic elements, such as Mn and Cu with Sn, Al, or Sb, could exhibit ferromagnetic behavior [ 8]. This phenomenon has since spurred extensive research on Heusler alloys, leading to the discovery of a range of remarkable properties [ 9], including ferromagnetism, significant spin-polarization effects [ 10], and promising semiconducting thermoelectric characteristics [ 11, 12]. Structurally, HH compounds, which combine excellent functional properties (e.g., thermoelectric performance) with good high-temperature strength, are currently a research focus. Among the thousands of known HH compounds, this class of materials exhibits both superior thermoelectric conversion efficiency and outstanding mechanical reliability, with particular attention paid to their fracture resistance under coupled high-temperature and stress conditions. As a new generation of functional materials, HH compounds have gained growing attention due to their unique magnetic and electronic properties, as well as excellent performance in terms of structural stability, mechanical strength, and thermal transport. Furthermore, their half-metallic character endows them with great potential for cutting-edge fields like spintronics [ 13, 14], making them a key focus of both fundamental and applied research. The structural integrity and service life of HH alloys in practical applications are critically dependent on mechanical properties such as fracture toughness and hardness. A significant advantage of these materials is their capacity to tailor mechanical performance through rational design. Chemical substitution and structural tuning, two strategic modification approaches, constitute a fundamental method for the intentional regulation of material property profiles, enabling the targeted adjustment of material characteristics by altering their basic composition and core structure, thus facilitating their incorporation into durable devices requiring high mechanical stability [ 15]. Additionally, HH compounds have emerged as promising, cost-effective, and eco-friendly materials for thermoelectric (TE) applications. They are widely used in power generators for waste heat recovery systems and cooling devices utilizing the Peltier effect, which realizes energy conversion by exploiting temperature gradients [ 16, 17]. The inherent tunability and multifunctionality of HH compounds have become the main driving force for the ongoing exploration of novel HH systems, with the optimization of TE properties as a central research objective [ 17]. Beyond superior functional properties, excellent mechanical characteristics, especially fracture toughness, are essential to ensure the long service life of advanced materials. However, current research on Heusler alloys has mainly focused on their excellent thermoelectric and optical properties, with relatively limited investigations into their fracture toughness. This research gap hinders the accurate prediction of the damage tolerance of such materials in practical service. Along with the continuous development of novel alloy systems, considerable research efforts are also devoted to improving the mechanical properties of mature thermoelectric (TE) materials, as their practical application and commercial scalability require viable candidates to meet key prerequisites, including facile synthesis, high structural stability, and tunable electronic and phonon properties [ 17]. Although previous studies have extensively investigated the thermoelectric and optical properties of HH compounds [ 18, 19, 20], research on their fracture toughness, a critical parameter determining the reliability of structural applications, remains notably scarce. To the best of our knowledge, only a small number of studies have reported the evaluation of fracture toughness for HH alloys. 2. Methods and Calculation Details All first-principles calculations in this study were performed using the VASP code, with the exchange-correlation functional treated by the PBE generalized gradient approximation (GGA) [ 26]. The plane-wave cutoff energy was set to 500 eV, and a 12 × 12 × 12 k-point mesh was employed for Brillouin zone integration [ 27]. Convergence tests performed with respect to cutoff energy and k-point mesh (using MnCoSb) showed that increasing the cutoff to 600 eV or the k-mesh to 16 × 16 × 16 changes the total energy by less than 0.5 meV/atom, confirming that the chosen parameters are well converged (see Supplementary Tables S1 and S2). The convergence threshold for electronic self-consistency was set to 10 −6 eV, and the force tolerance for ionic relaxation was set to 0.0001 eV/Å. The projector augmented wave (PAW) method was adopted to describe the electron–ion interactions within the system [ 28]. Geometric structure optimization, including key structural parameters such as lattice constants and atomic positions, was carried out using the BFGS algorithm [ 29]. Phonon spectra were calculated using density functional perturbation theory (DFPT) as implemented in VASP. A 2 × 2 × 2 supercell containing 96 atoms was constructed, and atomic displacements of 0.015 Å were applied to compute the force constants. A 6 × 6 × 6 q-point mesh was used for integration over the Brillouin zone. The phonon dispersion curves were then interpolated along the high-symmetry path Γ–X–M–Γ–R–X using the Phonopy code. Elastic constants were determined using the stress–strain method. Six independent strain patterns were applied with a maximum amplitude of ±0.5% in steps of 0.1%. Cubic symmetry constraints were enforced, allowing only the three independent elastic constants C11, C12, and C44 to be extracted from the linear regime of the stress–strain curves. 3. Results and Discussions 3.1. Crystal Structure and Stability The ferromagnetic HH alloys MnNiSb and MnCoZ (Z = Sb, P, As) were selected as the subjects of this investigation. The intrinsic phases of MnNiSb, MnCoSb, MnCoP, and MnCoAs all crystallize in the cubic space group F-43m (#216), with their unit cell containing a total of 12 atoms. The alloy structure follows the XYZ-type configuration. Three distinct elements (X, Y, Z) form a cubic crystal system with a uniquely ordered arrangement ( Figure 1a–d). Specifically, X, Y, and Z atoms are located at the face-centered cubic (FCC) lattice sites, tetrahedral interstitial sites, and octahedral interstitial sites, respectively. This configuration ensures the structural integrity as illustrated. In the MnCoZ (Z = Sb, P, As) system, the Sb, P, or As atoms occupy the 4a Wyckoff position (0, 0, 0), Co occupies the 4b position (1/4, 1/4, 1/4), and Mn occupies the 4c position (1/2, 1/2, 1/2) [ 30]. The MnCoP and MnCoA compounds were modeled by substituting the Sb atom in MnCoSb, whereas in MnNiSb, Ni occupies the 4a site (0, 0, 0), Mn the 4b site (1/4, 1/4, 1/4), and Sb the 4c site (1/2, 1/2, 1/2). The HH compounds investigated in this study feature a characteristic atomic arrangement: each elemental component forms an independent FCC Bravais lattice, creating a structure analogous to NaCl. As detailed in Table 1, the lattice parameters derived from geometric optimization show close consistency with available reference data [ 16, 21, 22, 31, 32], thereby validating our computational approach. This distinctive crystal structure underpins the excellent physicochemical properties observed, indicating significant potential for applications in thermoelectric and magnetic fields. Here, Etot represents the total energy of the compound system, while Ei and ni denote the energy of an isolated atom of the i-th element and its stoichiometric content in the unit cell, respectively. Table 1 lists all the calculated enthalpies of formation, which are consistently negative, confirming the excellent thermodynamic stability of these materials. Further analysis of the enthalpy data indicates that the calculated formation energies ( Table 1) are negative for MnCoSb, MnCoAs, MnCoP, and MnNiSb, indicating that they are synthesizable and exhibit good environmental stability. Based on this characteristic, we conducted an in-depth investigation of the electronic band characteristics of the studied systems using first-principles calculations. As shown in Figure 3a–h, MnCoSb, MnCoAs, MnCoP, and MnNiSb exhibit half-metallic properties, which is consistent with previous reports on such materials [ 16, 22, 32]. The spin-up channel displays metallic characteristics, with its valence and conduction bands intersecting at the Fermi level, indicating the presence of a high concentration of mobile charge carriers. In contrast, the spin-down channel exhibits a distinct band gap at the Fermi level, suggesting semiconductor or insulator behavior for the spin-down states. This asymmetry between the spin-up and spin-down channels is a hallmark of half-metallicity and confirms the presence of ferromagnetism in all four alloys. As illustrated in Figure 3a–h, the HH alloys MnCoSb, MnCoAs, MnCoP, and MnNiSb are predicted to exhibit 100% spin polarization at the Fermi level. Consequently, these thermodynamically stable half-metallic alloys are potential candidates for spintronic devices such as spin valves, magnetic random-access memories, and spin injectors. Table 1. Key computational findings, including lattice parameters, formation energy, and assessments of excellent dynamic stability under various loading conditions. Table 1. Key computational findings, including lattice parameters, formation energy, and assessments of excellent dynamic stability under various loading conditions. Crystal Structure Types a = b = c (Å) EFormation (eV/atom) Dynamical Stability MnCoSb 5.812 −3.54 YES 5.825 [ 34] -- -- 5.807 [ 31] -- -- 5.872 [ 22] -- -- 5.835 [ 16] MnCoAs 5.532 −5.987 YES 5.55 [ 22] -- -- MnCoP 5.351 −25.556 YES 5.37 [ 22] -- -- MnNiSb 5.906 −4.324 YES 5.93 [ 32] -- -- 5.91 [ 21] -- -- 3.2. Mechanical Properties Having confirmed the structural stability of these alloys, we further evaluated their macroscopic mechanical properties, which are essential for their potential applications as structural materials. Elastic constants [ 32], defined as manifestations of structural integrity, govern a material’s resistance to deformation, dictate its stress response, and establish fundamental criteria for judging the mechanical stability of the lattice [ 41]. Among the core mechanical characteristics, Young’s modulus ( E), the shear modulus, and Poisson’s ratio ( v) are crucial for assessing mechanical performance. Young’s modulus ( E) quantitatively describes a material’s resistance to tensile and compressive deformation, with its magnitude directly reflecting the inherent stiffness of the material. The shear modulus is another critical elastic parameter; it specifically characterizes the material’s intrinsic resistance to shear deformation, thereby quantifying its dimensional stability under shear stress. Poisson’s ratio, a parameter central to elasticity theory, defines the extent of transverse deformation relative to axial deformation under tensile or compressive loading. Collectively, these elastic constants form a constitutive framework that comprehensively characterizes the fundamental mechanical behavior of a material. The application of the data presented in Table 2 to the stability criteria confirms compliance with all necessary conditions. The calculated elastic stiffness constants for the material were determined to be the MnNiSb and MnCoZ (Z = Sb, P, As) alloys, which fulfill the theoretical requirements of the Born rule, demonstrating the mechanical robustness of the intrinsic structures of all four materials. Within the framework of the Voigt–Reuss–Hill approximation scheme [ 44, 45, 46], the bulk modulus and shear modulus [ 46] are computed from the derived data using the following relationships: B V = B R = C 11 + 2 C 12 3 (3) G V = C 11 − C 12 + 3 C 44 5 (4) G R = 5 C 44 ( C 11 − C 12 ) 4 C 44 + 3 C 11 − C 12 (5) Among the physical quantities, such as BV, BR, GV, and GR, the subscripts “V” and “R” represent the Voigt and Reuss approximations, respectively, the former establishing the theoretical upper bound and the latter the lower bound for the effective elastic moduli of polycrystalline materials of the Voigt–Reuss–Hill approximation method. By taking the arithmetic average of the values derived from the Voigt and Reuss models, one obtains the Hill model [ 38]: B H = B V + B R 2 (6) G H = G V + G R 2 (7) The values of E and ν are determined through the rigorous application of foundational theoretical relations connecting them to the parameters B and G [ 47]: E = 9 G B G + 3 B (8) v = 3 B − 2 G 6 B + 2 G (9) Here, k = G / B . The Cauchy stress-defined pressure ( CP), a variable widely used in the intersection of materials science and solid mechanics, functions as a key parameter for predicting the ductile versus brittle behavior of a material. For intermetallic compounds and complex alloys, the Cauchy index serves as an indicator for the ductile versus brittle behavior and provides significant practical utility by offering a metric to predict the mechanical response, particularly the ductile versus brittle character: C P = C 12 − C 44 [ 50]. Substituting the C12 and C44 values of the alloy’s materials from Table 2 into this formula yields positive Cauchy pressure values. The observed enhancement in ductility implies that deformation is primarily accommodated by dislocation glide, suppressing brittle fracture. The calculated mechanical parameters, including the bulk modulus ( B), shear modulus ( G), Young’s modulus ( E), Poisson’s ratio ( ν), Cauchy pressure, B/ G ratio, Cauchy pressure ( Cp), and Vickers hardness ( Hv), are systematically summarized for the alloy materials in Table 3. Based on (DFT) calculations, the intrinsic elastic moduli and related physical properties of MnNiSb and MnCoZ (Z = Sb, P, As) Heusler alloys were obtained, with the calculated results summarized in Table 3. As two key macroscopic mechanical parameters, the bulk modulus ( B) and shear modulus ( G) are directly determined from the set of fundamental elastic constants. They serve as core parameters for evaluating a material’s resistance to volumetric and shear deformation, where the bulk modulus represents the resistance to volume change under elastic deformation, and the shear modulus reflects the resistance to shape change under plastic deformation. Among the alloys investigated, MnCoP exhibits the highest bulk modulus (175.87 GPa), indicating its excellent compressive resistance and high stiffness. In contrast, MnCoSb shows the lowest bulk modulus (126.91 GPa). Regarding the shear modulus ( G), which reflects the material’s resistance to plastic deformation, MnNiSb demonstrates a value (39.64 GPa) considerably higher than those of the other three compounds, suggesting that dislocation motion in MnNiSb may be more difficult, potentially affecting its plastic deformation mechanism. Additionally, Young’s modulus ( E), an important parameter for assessing tensile stiffness, varies significantly among these alloys, consistent with the overall trends observed for the bulk modulus ( B) and shear modulus ( G). As illustrated in Figure 4, three-dimensional surface representations of the shear ( Figure 4a), bulk ( Figure 4b), and Young’s ( Figure 4c) modulus are presented. MnCoSb exhibits a moderate bulk modulus (126.91 GPa) but high mechanical anisotropy, as reflected by its three-dimensional modulus surfaces, whereas MnCoAs and MnNiSb exhibit progressively less mechanical anisotropy, and MnCoP shows the most pronounced mechanical anisotropy, consistent with the data summarized in Table 3. The combination of a high bulk modulus and a high Young’s modulus implies that these materials can maintain structural integrity when employed as structural components or subjected to stress conditions. Furthermore, a comprehensive understanding of the anisotropy of Young’s modulus is crucial for the future design of microelectronic or functional devices based on these materials, as device performance may be influenced by grain orientation and the direction of applied stress fields. This favorable mechanical stability, together with the thermodynamic stability and half-metallic electronic characteristics discussed above, collectively establishes the potential of these HH alloys for applications in advanced spintronic devices. The ductile or brittle behavior of materials can be empirically assessed using Poisson’s ratio ( v) and the B/ G ratio. According to Pugh’s criterion, a B/ G ratio greater than 1.75 typically indicates ductile behavior, whereas a lower value suggests brittleness. In this study, the B/ G ratios of all four alloys significantly exceed 1.75, ranging from 1.93 to 2.50, and their Poisson’s ratios ( v) are all greater than 0.3. These results indicate that MnNiSb, MnCoSb, MnCoP, and MnCoAs are expected to exhibit good ductility. This finding is supported by the fact that their Cauchy pressures ( Cp) are all significantly positive, further confirming their ductile behavior and the dominance of metallic bonding in these alloys. However, it is worth noting that although MnNiSb has a B/ G ratio of 1.97 and a Poisson’s ratio of 0.38, its relatively high shear modulus suggests that it may exhibit a stronger tendency toward brittleness compared to the other three alloys, possibly implying mixed metallic–covalent bonding characteristics. It should be emphasized that all four alloys remain within the ductile regime according to the Pugh criterion; the differences among them represent a spectrum within ductility rather than a transition to brittle behavior. MnNiSb is simply less ductile (i.e., relatively more brittle) than the others, but it is still ductile overall. In terms of stiffness, the calculated Young’s modulus ( E) shows that MnNiSb possesses the highest value (109.48 GPa), indicating its strongest resistance to elastic deformation, whereas MnCoP exhibits the lowest stiffness (80.25 GPa). Regarding Vickers hardness ( Hv), the values for the four alloys fall within a narrow range (23.61–27.7 GPa). Overall, the theoretical calculations suggest that MnCoSb, MnNiSb, MnCoP, and MnCoAs are all predicted to exhibit good ductility. 3.3. Electronic Structure The differences in macroscopic mechanical properties essentially originate from the microscopic electronic structure and chemical bonding characteristics of the materials. To further understand the aforementioned mechanical behavior, particularly the differences in the ductility–brittleness tendency between MnNiSb and the other alloys, we calculated and analyzed the density of states (DOS) of all investigated systems. The variations in mechanical properties are primarily attributed to the influence of different elemental compositions on the electronic structure and interatomic bonding characteristics [ 51, 52, 53]. As confirmed by the electronic DOS results of the alloys shown in Figure 5, the fundamental reason for the differences in macroscopic mechanical properties lies in the unique electronic structures and chemical bonding characteristics arising from their distinct chemical compositions. The DOS analysis reveals that all systems exhibit residual DOS at the Fermi level, indicating their metallic nature. Specifically, for MnCoSb ( Figure 5a,b), the DOS distribution near the Fermi level is relatively steep, with the Mn- d and Sb- p orbitals ( Figure 6a) showing strong hybridization peaks and a high degree of energy overlap. This suggests the presence of strong directional covalent bonds between atoms, which is closely related to the alloy’s high elastic modulus and inherent brittleness. In contrast, MnNiSb, MnCoP, and MnCoAs exhibit a flatter DOS distribution near the Fermi level, indicating greater electron delocalization [ 54, 55] and thus more pronounced metallic bonding characteristics. A correlation analysis was subsequently conducted between the electronic configurations and macroscopic mechanical properties of the HH alloys (MnNiSb, MnCoAs, MnCoSb, MnCoP). The element-resolved DOS results shown in Figure 6a–d demonstrate that different elemental compositions have a decisive influence on the distribution of electronic states near the Fermi level. The Mn- d orbitals ( Figure 6a–d) exhibit significant spin polarization in the energy range from −2 to 1 eV, with the spin-up and spin-down d orbitals distinctly separated, reflecting the crystal field splitting effect, and thus the electronic origin of the magnetic anisotropy in these materials [ 56, 57]. Notably, the difference between the transition metals Co ( Figure 6a–c) and Ni ( Figure 6d) plays a significant role in modulating the electronic structure: the Co-3 d orbitals give rise to new electronic states near the Fermi level, shifting and broadening the DOS peaks. The d- d coupling likely enhances the metallic character of the material, which is consistent with theoretical trends in electrical conductivity [ 58]. Furthermore, the differences in spin-resolved DOS at the Fermi level are particularly pronounced for the P-based system ( Figure 6c), suggesting that certain compositions may exhibit half-metallicity, where electrons of one spin orientation display metallic behavior while those of the opposite spin exhibit insulating characteristics. This property holds great potential for spintronic applications and further confirms that MnCoP and MnCoAs, obtained by substituting Sb in MnCoSb, retain the half-metallic nature of the parent compound. The difference in the ductility–brittleness tendency of the MnNiSb alloy is directly reflected in its DOS plots ( Figure 5g,h): a prominent peak originating from the Ni-3d orbitals is observed at 1–2 eV above the Fermi level. This corresponds to the occupation of antibonding states arising from the additional d electrons of Ni compared to Co [ 59]. The occupation of these antibonding states directly weakens the net bond strength of the Mn-Ni and Ni-Sb bonds, leading to changes in elastic moduli, higher stiffness, and a tendency toward brittleness in terms of the ductility brittleness tendency [ 54, 55]. MnCoP exhibits the most distinctive mechanical behavior, with the lowest strength and most pronounced anisotropy, but potentially superior plastic deformability. Its DOS plots show a moderate total DOS at the Fermi level, and the d-band centers of Mn and Co are located farther from the Fermi level compared to the other systems. Meanwhile, the hybridization peaks between the P p-orbitals, Figure 6c, and the transition metal d-orbitals are broader and shallower. This electronic structure suggests a bonding nature with less covalency and a stronger tendency toward metallicity. The weaker bond directionality likely allows dislocations to nucleate under lower stresses, while the metallic character facilitates dislocation slip and multiplication. Consequently, more significant plastic flow during deformation can dissipate energy, potentially contributing to greater fracture strain [ 60, 61]. Magnetic properties. Since all four alloys exhibit half-metallic ferromagnetism, we further calculated their total magnetic moments and atomic-resolved contributions. The total magnetic moments per formula unit are MnCoSb = 4.00 μ B, in agreement with the literature [ 16, 17, 22]; MnCoAs = 4.00 μ B [ 22]; MnCoP = 3.95 μ B [ 22]; and MnNiSb = 4.00 μ B [ 32]. The magnetic moments are mainly contributed by the Mn atoms (≈3.2–3.5 μ B), while Co and Ni atoms provide small opposite contributions (≈−0.2 μ B); the p-block elements exhibit negligible magnetic moments. These values follow the Slater–Pauling rule ( Zt − 18 = 4 μ B) [ 32]. To assess the influence of magnetism on elastic constants, we performed a control calculation for MnCoSb by switching off spin polarization (non-magnetic state). The resulting elastic constants C11 = 170.32 GPa, C12 = 105.88 GPa, and C44 = 35.40 GPa change by less than 5%, and the B/ G ratio remains above 1.75. Hence, magnetism does not dominate the mechanical behavior; the ductile nature and fracture toughness reported herein are robust against magnetic ordering. A similar weak dependence of elastic constants on magnetism has been observed in other HH systems [ 32, 62]. In summary, the range of macroscopic mechanical properties (strength, plasticity) of these HH alloys is precisely governed by the DOS distribution near the Fermi surface and the intensity and characteristics of p- d orbital hybridization. Strong hybridization and a high DOS at the Fermi level promote high strength and high stiffness, potentially at the expense of ductility, whereas moderate hybridization combined with enhanced metallic character favors a better strength–ductility synergy. The electronic origin of these property differences ultimately determines the material’s resistance to crack propagation, i.e., its fracture toughness. 3.4. Fracture Toughness The fracture resistance of a material is crucial to its service safety and structural reliability. For load-bearing alloys, their ability to resist unstable crack propagation is commonly characterized by the plane-strain fracture toughness ( KIC), i.e., the critical stress intensity factor. The theoretical basis of this parameter stems from Griffith’s energy balance theory [ 63], which states that the significant discrepancy between the practical and theoretical fracture strengths of materials originates from stress concentration at the tips of pre-existing microcracks under external loading. Crack propagation becomes unstable when the elastic strain energy released by the system due to crack extension is sufficient to compensate for the energy required to create new fracture surfaces, a critical condition known as the Griffith criterion. From an engineering perspective, fracture toughness KIC is an important mechanical property that quantifies a material’s resistance to rapid fracture containing defects or cracks. It is defined as the critical stress intensity factor, KIC. A higher KIC value indicates a greater ability to suppress crack growth, thereby enhancing component safety [ 64, 65]. In this study, the specific crystallographic direction (001) was selected as a potential cleavage plane, and a crack model was constructed accordingly. Specifically, a 1 × 1 × 1 supercell was used, with a vacuum layer of 15 Å added perpendicular to the crack plane to avoid periodic interactions. All atoms within a 5 Å radius of the crack tip were fully relaxed until the residual forces were less than 0.01 eV/Å. The crack formation energy Ecrack was then calculated as the energy difference between the cracked and pristine supercells, divided by the crack area Acrack. G I C = E c r a c k − E i n i t i a l A c r a c k (12) where Ecrack represents the crack formation energy in electron volts (eV), defined as the energy difference between the cracked and initial systems; Einitial denotes the total energy of the pristine crystal structure in eV; and Acrack indicates the cross-sectional area of the crack in square angstroms (Å 2). The empirical formula (Method II) is as follows [ 68]: K IC = V 0 1 / 6 G ( B / G ) 1 / 2 (13) Here, V0 is the volume of the primitive cell ( V0 = a × b × c), and B and G are the bulk and shear moduli. To evaluate damage tolerance, the brittle index M is calculated as [ 39, 68] M = H V / K IC (14) The damage tolerance of a material can be assessed by the brittle index M, a parameter that includes the Vickers hardness, with a lower M value denoting superior tolerance. As summarized in Table 4, the key parameters required for calculating the mode I fracture toughness ( KIC) via two different methods are presented, along with the final results. The two methods yield different absolute values (e.g., a difference of about 40% for MnCoSb). Therefore, the explicit crack model (Method I) is preferred for quantitative assessment of fracture toughness in anisotropic HH alloys, whereas the empirical formula (Method II) should be used only as a rough comparative reference or for preliminary screening. In the study of HH alloys, accurate evaluation of fracture toughness ( KIC) and damage tolerance is crucial for assessing their potential as structural–functional integrated materials. This study presents two theoretical methods for estimating KIC based on intrinsic material parameters (atomic volume V0, bulk modulus B, and shear modulus G), establishing a bridge between electronic/atomic-scale calculations and macroscopic fracture properties. To further quantify the brittle ductile characteristics, the brittleness index M is introduced, where lower values indicate superior damage tolerance and enhanced resistance to brittle fracture while maintaining high hardness. As summarized in Table 4 for the MnNiSb and MnCoZ (Z = Sb, P, As) alloy systems, the calculated results indicate that among these four alloys, both models confirm that MnCoSb exhibits the most outstanding comprehensive performance. It possesses the lowest brittleness index ( M = 8.97 I, 14.94 II) and relatively high fracture toughness ( KIC = 2.63 I MPa·m 1/2, 1.58 II MPa·m 1/2), suggesting excellent damage tolerance. In contrast, MnCoP exhibits the lowest KIC values (2.00 I MPa·m 1/2,1.63 II MPa·m 1/2) and the highest M values among the four half-metals (13.85 I, 16.99 II), implying a greater tendency toward brittle fracture and potentially inferior mechanical reliability. The properties of MnNiSb ( KIC = 2.55 I MPa·m 1/2,1.90 II MPa·m 1/2; M = 10.54 I, 14.14 II) and MnCoAs ( KIC = 2.28 I MPa·m 1/2, 1.62 II MPa·m 1/2; M = 11.61 I, 16.32 II) lie between these two alloys. Overall, the fracture toughness values of the four half-metallic materials are consistent with the expected calculated results. The comparison between the two methods, despite differences in absolute values, confirms that the explicit crack model (Method I) is more reliable for anisotropic systems because it directly simulates crack propagation, accounts for local atomic relaxations, and does not rely on isotropic assumptions. In contrast, Method II serves as a rough comparative reference due to its simplicity. The ability to quantitatively predict key mechanical properties of XYZ-type intermetallic compounds through DFT calculations provides critical parameters and theoretical insights. This enables the rational design and strategic screening of novel alloys with excellent damage tolerance based on fundamental theoretical principles. 4. Conclusions A systematic evaluation of the HH alloys MnCoSb, MnCoAs, MnCoP, and MnNiSb using DFT confirms their excellent thermodynamic and mechanical stability. The negative formation energies and the absence of imaginary frequencies in the phonon spectra verify their synthesizability and dynamic stability under ambient conditions. Furthermore, their elastic constants satisfy the Born stability criteria for cubic crystals, thereby ensuring mechanical stability. The fracture toughness of all alloys lies in the range of 1.58 MPa·m 1/2–2.63 MPa·m 1/2, and their Poisson’s ratios and B/ G ratios exceed the critical value of 1.75. Notably, MnNiSb exhibits a relatively high tendency toward brittleness. Nevertheless, all alloys satisfy the B/G > 1.75 criterion, confirming that they are all ductile; the observed differences reflect a spectrum within the ductile regime. The higher tendency toward brittleness of MnNiSb is attributed to its relatively high shear modulus (39.64 GPa) and the significant covalent bonding character revealed by its electronic structure. Electronic structure characterization confirms the half-metallic nature of all compounds, with a theoretical spin polarization of 100%. For fracture toughness, the explicit crack model (Method I) is physically more realistic for anisotropic HH alloys, while the empirical formula (Method II) serves as a convenient reference. Although the two methods give different absolute values, the comparison confirms that Method I is preferred for quantitative assessment. The results indicate that MnCoSb exhibits the most outstanding overall performance, demonstrating the lowest brittleness index M (8.97, 14.94) and relatively high KIC values (2.63 MPa·m 1/2, 1.58 MPa·m 1/2), suggesting its excellent damage tolerance. Among the other three half-metals, MnCoAs shows KIC values of 2.28 MPa·m 1/2 and 1.62 MPa·m 1/2 and M values of 11.61 and 16.32; MnCoP exhibits KIC values of 2.00 MPa·m 1/2 and 1.63 MPa·m 1/2 and the highest M values of 13.85 and 16.99; and MnNiSb shows KIC values of 2.55 MPa·m 1/2 and 1.90 MPa·m 1/2 and M values of 10.54 and 14.14. The fracture toughness results for the four half-metallic materials indicate that MnCoP exhibits a greater tendency toward brittleness compared to the other three half-metals, while MnCoSb demonstrates the most favorable overall properties. This study achieves a quantitative prediction of the fracture toughness of HH alloys through first-principles calculations, revealing the excellent fracture toughness and damage-tolerant nature of the MnCoSb HH alloy, thereby establishing a theoretical foundation for the design of functional materials with excellent damage tolerance. The proposed methodology and model provide a paradigm for the rational design of high-performance materials integrating structural and functional properties. Supplementary Materials The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules31121994/s1, Table S1: Convergence of total energy (in eV) with respect to plane wave cutoff energy for the four HH alloys, with the k point mesh fixed at 12 × 12 × 12. The optimized lattice constants a = b = c (Å) are also listed; Table S2: Convergence of total energy (in eV) with respect to k point mesh for the four HH alloys, with the plane wave cutoff energy fixed at 500 eV. The optimized lattice constants a = b = c (Å) are also listed. Author Contributions A.Q.: writing—review and editing, supervision, investigation. S.-B.C.: writing—review and editing, supervision, investigation, conceptualization. L.-Z.T.: writing—review and editing, writing—original draft, investigation. J.-H.W.: writing—review and editing, supervision. W.-J.Y.: writing—review and editing, investigation, formal analysis. T.G.: writing—review and editing, investigation, software. K.-M.G.: writing—review and editing, investigation, formal analysis. J.Z.: writing—review and editing. All authors have read and agreed to the published version of the manuscript. Funding This research received support from the Training Program of Innovation and Entrepreneurship for College Students (Grant Nos. 202510667481 and S2024106670020); the Scientific Research Training Program for College Students (Grant No. asxysrt202327); the National Natural Science Foundation of China (Grant No. 12364017); the Guizhou Provincial Department of Education Open Recruitment and Leadership Scientific and Technological Attack Project (Guizhou Education Technology [2024] No. 003); Anshun City’s 2025 Basic Research Project of the Science and Technology Bureau (Grant No. [2025]01); and the Guizhou Provincial Science and Technology Plan Project (Subject) (No.: Qian Ke He Ji Chu MS (2026) 011). Institutional Review Board Statement Not applicable. Informed Consent Statement Not applicable. Data Availability Statement This manuscript has associated data in a data repository. The raw data supporting the conclusions of this article will be made available by the authors on request. Conflicts of Interest The authors declare no conflicts of interest. References Palmstrøm, C.J. Heusler compounds and spintronics. Prog. Cryst. Growth Charact. Mater. 2016, 62, 371–397. [ Google Scholar] [ CrossRef] Ouahrani, T.; Otero-de-la-Roza, A.; Khenata, R.; Luaña, V.; Amrani, B. Structural and thermodynamic properties of SbAsGa 2 and SbPGa 2 chalcopyrites. Comp. Mater. Sci. 2010, 47, 655–659. [ Google Scholar] [ CrossRef] El Krimi, Y.; Masrour, R.; Jabar, A. Structural, electronic and magnetic properties of full-Heusler alloy Co 2CrAl. Inorg. Chem. Commun. 2020, 121, 108207. [ Google Scholar] [ CrossRef] Idrissi, S.; Ziti, S.; Labrim, H.; Bahmad, L.; El Housni, I.; Khalladi, R.; Mtougui, S.; El Mekkaoui, N. Half-metallic behavior and magnetic proprieties of the quaternary Heusler alloys YFeCrZ (Z = Al, Sb and Sn). J. Alloys Compd. 2020, 820, 153373. [ Google Scholar] [ CrossRef] Alam, J.; Bran, C.; Chiriac, H.; Lupu, N.; Óvári, T.A.; Panina, L.V.; Rodionova, V.; Varga, R.; Vazquez, M.; Zhukov, A. Cylindrical micro and nanowires: Fabrication, properties and applications. J. Magn. Magn. Mater. 2020, 513, 167074. [ Google Scholar] [ CrossRef] Tian, X.; Xu, H.; Yang, J.; Zhao, L.; Wang, X.; Na, S.; Jin, X.; Li, J.; Zhao, W.; Tan, C. New all-d-metal Co-V-Ti Heusler alloys: Superior mechanical strength and excellent elastocaloric effect. Mater. Sci. Eng. A 2025, 926, 147965. [ Google Scholar] [ CrossRef] Dahmane, F.; Mogulkoc, Y.; Doumi, B.; Tadjer, A.; Khenata, R.; Bin Omran, S.; Rai, D.P.; Murtaza, G.; Varshney, D. Structural, electronic and magnetic properties of Fe 2-based full Heusler alloys: A first principle study. J. Magn. Magn. Mater. 2016, 407, 167–174. [ Google Scholar] [ CrossRef] Moussali, A.; Amina, M.B.; Fassi, B.; Ameri, I.; Ameri, M.; Al-Douri, Y. First-principles calculations to investigate structural and thermodynamic properties of Ni 2LaZ (Z = As, Sb and Bi) Heusler alloys. Indian J. Phys. 2020, 94, 1733–1747. [ Google Scholar] [ CrossRef] Idrissi, S.; Labrim, H.; Ziti, S.; Bahmad, L. Structural, electronic, magnetic properties and critical behavior of the equiatomic quaternary Heusler alloy CoFeTiSn. Phys. Lett. A 2020, 384, 126453. [ Google Scholar] [ CrossRef] Jiang, E.-H.; Zhu, X.-F.; Chen, L.-F. First-principles study of the electronic structure, magnetism, and spin-polarization in Heusler alloy Co 2MnAl(100) surface. Acta Phys. Sin. 2015, 64, 147301. [ Google Scholar] [ CrossRef] Xu, J.-L.; Jia, L.-Y.; Jin, X.-Q.; Hao, X.-N.; Ma, L.; Hou, D.-L. Structure and half-metallic ferromagnetism of quaternary Heusler compounds CoMnZnZ. Acta Phys. Sin. 2019, 68, 157501-1–157501-07. [ Google Scholar] [ CrossRef] Kulkova, S.E.; Eremeev, S.V.; Kulkov, S.S.; Skripnyak, V.A. Ab initio investigations of magnetic properties of thin film Heusler alloys. Mater. Sci. Eng. A 2008, 481–482, 209–213. [ Google Scholar] [ CrossRef] Zeier, W.G.; Schmitt, J.; Hautier, G.; Aydemir, U.; Gibbs, Z.M.; Felser, C.; Snyder, G.J. Engineering half-Heusler thermoelectric materials using Zintl chemistry. Nat. Rev. Mater. 2016, 1, 16032. [ Google Scholar] [ CrossRef] Kawasaki, J.; Chatterjee, S.; Canfield, P. Full and half-Heusler compounds. MRS Bull. 2022, 47, 555–558. [ Google Scholar] [ CrossRef] Zhang, Y.; Kumar, S.V.; Xiang, W.; Wu, Z.; Sun, Z. Martensitic transformation and mechanical properties of grain refined Ni–Co–Mn–Sn Heusler alloys via Cr doping. Mater. Sci. Eng. A 2021, 804, 140777. [ Google Scholar] [ CrossRef] Hori, A.; Minami, S.; Saito, M.; Ishii, F. First-principles calculation of lattice thermal conductivity and thermoelectric figure of merit in ferromagnetic half-Heusler alloy CoMnSb. Appl. Phys. Lett. 2020, 116, 242408. [ Google Scholar] [ CrossRef] Gupta, V.K.; Kumar, R.; Kumar, U. Unraveling the multifunctional capabilities of MnCoSb half-Heusler alloy: A computational exploration of structural, mechanical, optical, and thermoelectric properties. J. Phys. Chem. Solids 2026, 208, 113071. [ Google Scholar] [ CrossRef] Chen, S.-B.; Yan, W.-J.; Gao, T.-H. Electronic, elastic, optic, thermal dynamic, and thermoelectric properties of the half-Heusler compounds ZrCoBi and TiCoBi. Eur. Phys. J. Plus 2023, 138, 966. [ Google Scholar] [ CrossRef] Li, W.; Ghosh, S.; Liu, N.; Poudel, B. Half-Heusler thermoelectrics: Advances from materials fundamental to device engineering. Joule 2024, 8, 1274–1311. [ Google Scholar] [ CrossRef] Wang, J.-H.; Feng, C.-J.; Zhang, L.; Chen, Y.; Chen, S.-B. Effect of Co, Ni doping on the thermodynamic properties of half-Heusler alloy FeNbSb: A first-principles study. J. Mater. Sci. 2025, 60, 9327–9343. [ Google Scholar] [ CrossRef] Wang, Y.-M.; Pei, H.-X.; Ding, J.; Wen, L.-W. First-principles study of magnetism and electronic structureof Sb-containing half-Heusler alloys. Acta Phys. Sin. 2011, 60, 047110. [ Google Scholar] [ CrossRef] Xu, B.; Zhang, M. Effects of structural change on half-metallicity of CoMnZ (Z = P, As, Sb): A first-principle study. J. Magn. Magn. Mater. 2011, 323, 939–944. [ Google Scholar] [ CrossRef] Li, G.; An, Q.; Duan, B.; Borgsmiller, L.; Malki, M.A.; Agne, M.; Aydemir, U.; Zhai, P.; Zhang, Q.; Morozov, S.I.; et al. Fracture toughness of thermoelectric materials. Mater. Sci. Eng. 2021, 144, 100607. [ Google Scholar] [ CrossRef] Zou, M.; Bao, L.; Yang, A.; Duan, Y.; Peng, M.; Cao, Y.; Li, M. Structural, elastic, mechanical, electronic, damage tolerance, fracture toughness, and optical properties of Cr–Al–B MAB phases studied by first-principles calculations. J. Mater. Res. 2023, 38, 1396–1409. [ Google Scholar] [ CrossRef] Wang, B.L.; Li, J. Fracture behavior and thermal shock resistance analysis of thermoelectric material plates and shells under thermal and electric shocks. Eng. Fract. Mech. 2020, 225, 106130. [ Google Scholar] [ CrossRef] Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [ Google Scholar] [ CrossRef] [ PubMed] Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [ Google Scholar] [ CrossRef] Gajdoš, M.; Hummer, K.; Kresse, G.; Furthmüller, J.; Bechstedt, F. Linear optical properties in the projector-augmented wave methodology. Phys. Rev. B 2006, 73, 045112. [ Google Scholar] [ CrossRef] Fischer, T.H.; Almlof, J. General methods for geometry and wave function optimization. J. Phys. Chem. Solids 1992, 96, 9768–9774. [ Google Scholar] [ CrossRef] Jain, A.; Ong, S.P.; Hautier, G.; Chen, W.; Richards, W.D.; Dacek, S.; Cholia, S.; Gunter, D.; Skinner, D.; Ceder, G.; et al. Commentary: The Materials Project: A materials genome approach to accelerating materials innovation. Appl. Mater. 2013, 1, 011002. [ Google Scholar] [ CrossRef] V, A.S.; Sarwan, M.; Singh, S. Ground State Structural, Elastic, Electronic Properties and Pressure-Induced Structural Phase Transition of XCoSb (X = Sc, Ti, V, Cr and Mn). J. Supercond. Nov. Magn. 2020, 33, 1821–1829. [ Google Scholar] [ CrossRef] Toual, Y.; Mouchou, S.; Rani, U.; Azouaoui, A.; Hourmatallah, A.; Masrour, R.; Rezzouk, A.; Bouslykhane, K.; Benzakour, N. Probing electronic, magnetic and thermal properties of NiMnSb Half-Heusler alloy for spintronics as an environmentally friendly energy resource: A DFT + U and Monte Carlo study. Mater. Today Commun. 2024, 38, 108064. [ Google Scholar] [ CrossRef] Toual, Y.; Mouchou, S.; Azouaoui, A.; Fakrach, B.; Masrour, R.; Bouslykhane, K.; Rezzouk, A.; Benzakour, N.; Hormatallah, A. DFT analysis of defect-induced changes in half-metallicity of the perfect structure α-CoMnSb half-Heusler for smart spintronic device applications. Comput. Condens. Matter 2025, 44, e01069. [ Google Scholar] [ CrossRef] Baral, M.; Ganguli, T.; Chakrabarti, A. Investigation of structural, magnetic and electronic properties of CoMnSb superstructure: A DFT study. Comput. Mater. Sci. 2022, 210, 111441. [ Google Scholar] [ CrossRef] Uto, O.T.; Adebambo, P.O.; Akinlami, J.O.; Adebayo, G.A. Predicting the stable Type-I phase of XMnSb (X = Co, Fe, Os) compounds and its thermodynamic, electronic and magnetic properties from first-principles calculations. Solid State Sci. 2020, 105, 106208. [ Google Scholar] [ CrossRef] Wu, Y.; Wu, B.; Wei, Z.; Zhou, Z.; Zhao, C.; Xiong, Y.; Tou, S.; Yang, S.; Zhou, B.; Shao, Y. Structural, half-metallic and elastic properties of the half-Heusler compounds NiMnM (M = Sb, As and Si) and IrMnAs from first-principles calculations. Intermetallics 2014, 53, 26–33. [ Google Scholar] [ CrossRef] Wei, X.-P.; Zhang, X.; Shen, J.; Chang, W.-L.; Tao, X. Gilbert damping, electronic and magnetic properties for quaternary Heusler alloys CrYCoZ: First-principles and Monte Carlo studies. Comp. Mater. Sci. 2022, 210, 111453. [ Google Scholar] [ CrossRef] Chen, S.-B.; Guo, S.-D.; Lv, B.; Ang, Y.S. Electronic, optical, elastic, and thermal properties of the half-Heusler-derived Ti 2FeNiSb 2: A DFT study. Phys. B Condens. Matter. 2024, 689, 416201. [ Google Scholar] [ CrossRef] Li, R.; Duan, Y.; Ma, L.; Yang, A.; Zheng, S.; Su, L.; Su, L.; Li, M. Elastic properties, fracture toughness, electronic properties and thermal conductivity of In-doped η’-Cu 6Sn 5 in Sn-Cu solders using first-principles calculations. Mater. Today Commun. 2025, 43, 111634. [ Google Scholar] [ CrossRef] Chen, S.-B.; Wang, J.-H.; Liu, W.-J.; Yan, W.-J.; Gao, T.-H. First-principles investigation of mechanical properties and fracture toughness in doped FeNbSb half-Heusler alloys. iScience 2025, 28, 114189. [ Google Scholar] [ CrossRef] Hearmon, R.F.S. The Elastic Constants of Anisotropic Materials. Rev. Mod. Phys. 1946, 18, 409–440. [ Google Scholar] [ CrossRef] Rogge, S.M.J.; Waroquier, M.; Van Speybroeck, V. Reliably Modeling the Mechanical Stability of Rigid and Flexible Metal–Organic Frameworks. Acc. Chem. Res. 2018, 51, 138–148. [ Google Scholar] [ CrossRef] Mouhat, F.; Coudert, F.-X. Necessary and sufficient elastic stability conditions in various crystal systems. Phys. Rev. B 2014, 90, 224104. [ Google Scholar] [ CrossRef] Villalobos-Portillo, E.E.; Fuentes-Montero, L.; Montero-Cabrera, M.E.; Burciaga-Valencia, D.C.; Fuentes-Cobas, L.E. Polycrystal piezoelectricity: Revisiting the Voigt-Reuss-Hill approximation. Mater. Res. Express 2019, 6, 115705. [ Google Scholar] [ CrossRef] Man, C.-S.; Huang, M. A Simple Explicit Formula for the Voigt-Reuss-Hill Average of Elastic Polycrystals with Arbitrary Crystal and Texture Symmetries. J. Elast. 2011, 105, 29–48. [ Google Scholar] [ CrossRef] Hill, R. The Elastic Behaviour of a Crystalline Aggregate. Proc. Phys
