Previous Download PDF Next
Download PDF
Research Article Open Access 8 Sep 2026

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Views:327 Downloads:19 Cited: 0
J. Mater. Inf. 2026, 6, 44. 10.20517/jmi.2026.16
Article Notes

Graphical Abstract

Abstract

High-entropy alloys (HEAs) exhibit exceptional stability in extreme environments, yet their expansive design space presents a “curse of dimensionality” for traditional discovery methods. While machine learning (ML) offers a data-driven paradigm for material screening, the scarcity of experimental data often results in overfitting and limited physical interpretability. To address these challenges, this study proposes a hybrid physics-informed machine learning (Hybrid PIML) framework for accelerated hardness prediction. By integrating classical solid solution strengthening theory with a residual learning artificial neural network (ANN), the model explicitly embeds the physical coupling of shear modulus and lattice distortion (G·δr2/3) as prior knowledge. This approach ensures predictions adhere to metallurgical principles while significantly outperforming benchmark algorithms, achieving a coefficient of determination (R2) of 0.976 and reducing the root mean square error (RMSE) by approximately 43%. SHapley Additive exPlanations (SHAP) analysis confirms that physics-enhanced features dominate the decision-making process, validating the model’s internalization of strengthening mechanisms. Furthermore, the research elucidates a phase-dependent non-linear correlation between hardness and yield strength, correcting the failure of the classical Tabor formula in work-hardening face-centered cubic (FCC) alloys. Finally, a high-throughput virtual screening funnel based on this framework successfully identified optimized non-equiatomic candidates within the refractory Co-Cr-Ti-Mo-W system. This work establishes a precise, physically consistent pathway for inverse material design under data-constrained conditions.

Keywords

High-entropy alloysphysics-informed machine learninghardness predictionsolid solution strengtheninghigh-throughput virtual screeningresidual learning
Reprints
Download PDF

INTRODUCTION

Since the Bronze Age, human alloy design strategies have long followed the “single principal element” paradigm, basing materials on one metal matrix and adding minor alloying elements to fine-tune properties. However, in the early 21st century, the concept of high-entropy alloys (HEAs) proposed by Yeh et al. completely overturned this tradition[1]. HEAs are typically defined as alloy systems composed of five or more principal elements in equiatomic or near-equiatomic ratios[2-4]. This novel design philosophy has pushed the frontier of materials research and development from the corners of phase diagrams to their vast central regions[5,6].

The scientific value of HEAs stems primarily from their unique “four core effects”, which endow them with mechanical properties distinct from traditional alloys: High Entropy Effect: The high configurational entropy generated by multi-component mixing significantly lowers the system’s Gibbs free energy at high temperatures, suppressing the formation of brittle intermetallic compounds (IMCs) and promoting the stability of simple solid solution phases [e.g., face-centered cubic (FCC), body-centered cubic (BCC)][2,7]. Severe Lattice Distortion: Elements with different atomic sizes randomly occupy lattice sites, resulting in a severe, non-uniform elastic stress field within the lattice. This distortion significantly increases resistance to dislocation motion and is the primary physical origin of the high hardness and yield strength observed in HEAs[8-10]. Sluggish Diffusion: In complex multi-component environments, cooperative atomic diffusion becomes difficult, leading to elevated activation energy for lattice diffusion. This characteristic imparts HEAs with excellent high-temperature creep resistance and microstructural thermal stability[11-13]. Cocktail Effect: This is a composite synergistic effect, implying that the overall performance of the alloy is not merely a simple summation of its constituents’ properties but includes additional gains arising from element-element interactions[12,14,15].

Despite the promising prospects of HEAs, their design faces a formidable challenge: the “curse of dimensionality”. Assuming a selection of 30 common metal elements from the periodic table, the number of combinations for equiatomic alloys composed of just 5 elements exceeds 142,000. If non-equiatomic variations (e.g., the continuous variation of AlxCoCrFeNi) and complex processing heat treatments are considered, the potential design space explodes exponentially, potentially reaching millions or even billions of combinations[16-18].

Traditional “trial-and-error” experimental workflows - encompassing melting, casting, cutting, polishing, and characterization - are lengthy and costly, making them wholly inadequate for such massive screening tasks[19-21]. Even high-throughput computational methods based on CALPHAD are often limited by a lack of accurate multi-component thermodynamic databases, leading to reduced prediction accuracy in unknown compositional regions[22-24]. Therefore, there is an urgent need for an efficient screening tool to navigate this vast uncharted territory[21,25].

In recent years, machine learning (ML), as a data-driven “fourth paradigm”, has been widely applied in materials science. Algorithms such as random forest (RF), support vector machines (SVM), and gradient boosting decision trees [Light Gradient Boosting Machine (LightGBM)] have demonstrated powerful capabilities in establishing “composition-property” mappings[7,26].

However, in the field of HEAs, the application of ML faces a specific “small data” dilemma. Unlike image recognition fields with millions of labeled data points, high-quality data for HEA hardness or yield strength is often limited to a few hundred or a thousand entries[2,23]. This data scarcity introduces core problems:

Overfitting Risk: Models are prone to memorizing noise in the training set rather than learning the intrinsic laws of the material, leading to inflated performance on validation sets[23,27]. Poor Extrapolation: Pure data-driven models are essentially interpolators. When predicting compositions outside the coverage of training data, model predictions often fail completely, sometimes yielding results that violate physical common sense[28,29]. Lack of Interpretability: “Black box” models cannot explain why a specific composition possesses high hardness, making it difficult to guide new material design at the scientific level[30,31].

To overcome the aforementioned bottlenecks, physics-informed machine learning (PIML) has emerged as a transformative approach. In recent years, interpretable ML has proven invaluable for uncovering complex, non-linear structure-property relationships that defy simple empirical assumptions, such as the counterintuitive role of defect morphology in additive manufacturing[32]. Concurrently, PIML frameworks have emerged to solve the “small data” dilemma. For instance, physics-constrained neural networks (PCNNs) have been developed to integrate empirical physics equations into dual-loss functions, successfully predicting microstructure-dependent corrosion behaviors through inverse parameter identification[33]. However, while solid solution strengthening (SSS) in HEAs provides a classical linear baseline, the extreme lattice distortions and complex interatomic interactions introduce significant non-linear deviations. Therefore, a novel architecture is required to bridge macroscopic physical laws with microscopic non-linear complexities. The core idea of PIML is to embed domain knowledge into ML models, using known physical laws to constrain or guide the learning process, thereby reducing dependence on data volume and improving generalization ability[34,35].

In HEA hardness prediction, while SSS theory provides a solid physical foundation, its classical application poses a significant challenge. In traditional dilute alloys, SSS is often treated as a relatively simple linear or analytical problem based on the “average effective medium” assumption (e.g., the classical Labusch model)[9]. However, the extreme local lattice distortions, multi-principal element synergistic interactions, short-range order (SRO), and phase structural mutations inherently present in HEAs introduce massive non-linear deviations. Consequently, simple linear SSS models completely fail to capture these complex behaviors in expansive HEA compositional spaces. This fundamental limitation provides the critical motivation for introducing ML. A novel architecture is required to bridge the macroscopic linear physical laws with microscopic non-linear complexities. In our Hybrid PIML framework, classical analytical models (combining atomic size mismatch and modulus mismatch) are not discarded but utilized as an accurate “first-order linear approximation” or “physical baseline”, leaving the neural network to capture the highly complex non-linear residuals.

Beyond constructing a high-precision hardness prediction model, this work addresses the disconnect in current research by bridging the gap between numerical predictions and the understanding of the non-linear yield strength-hardness relationship. Additionally, it demonstrates the utility of this model for accelerating material discovery in complex systems like Refractory HEAs. A Hybrid PIML framework is presented herein, which integrates SSS theory by using a physics-based SSS model as a baseline and training the ML model to learn the ‘residual’ between experimental and theoretical values. This methodology seamlessly merges the robustness of physical principles with the flexibility of data-driven models, ensuring high-precision, physically consistent, and data-efficient predictions of HEA mechanical properties.

MATERIALS AND METHODS

Constructing a robust ML model requires, first and foremost, high-quality data infrastructure. Since raw compositional ratios cannot directly reflect internal physical interactions, Feature Engineering is essential to transform compositions into physical descriptors capable of describing phase stability and strengthening mechanisms.

Data acquisition and curation

Data quality and representativeness are the bedrock of high-precision models, especially for PIML frameworks relying on physical laws. The data must include not only chemical compositions but also key thermodynamic and physical parameters.

Dataset compilation

A dataset comprising 205 experimentally verified HEA samples was compiled for this study. The data were primarily curated from peer-reviewed literature and public materials databases[2,36-38]. This dataset encompasses diverse alloy systems, predominantly including 3d Transition Metal Systems (e.g., Al-Co-Cr-Fe-Ni and its derivatives, spanning FCC/BCC/Multi-phase structures) and Refractory HEA Systems (e.g., W-Mo-Ta-Nb-V-Ti systems, typically displaying a BCC structure)[39,40]. Each data entry records the complete chemical composition (in atomic percentage) alongside the corresponding experimental Vickers hardness (HV). To facilitate the training of the physics-informed model, the feature space was not limited to raw elemental ratios. Instead, 12 key physical and thermodynamic descriptors were calculated for each alloy utilizing CALPHAD methods and empirical formulas[8,23]. Core features include: Atomic Size Difference (δr), which quantifies the degree of lattice distortion and serves as a critical parameter for SSS; Mixing Enthalpy (ΔHmix) and Mixing Entropy (ΔSmix), which characterize thermodynamic phase stability; valence electron concentration (VEC), a key electronic parameter for predicting FCC/BCC phase stability; and Shear Modulus (G), reflecting the material’s resistance to elastic deformation.

Data cleaning and preprocessing

To eliminate noise and ensure training stability, a strict cleaning process was executed: Missing and Outlier Treatment: Samples with missing physical parameters or hardness records were removed. Outliers significantly deviating from physical common sense were excluded based on the 3σ principle[12,41]. Physical Consistency Check: All physical descriptors were verified to ensure inputs complied with the basic assumptions of the SSS theory.

Feature standardization

Given the immense dimensional differences among physical descriptors (e.g., ΔH is roughly 101 kJ/mol, while δr is 10-2), direct input into a neural network would cause convergence difficulties or weight bias. Therefore, Z-Score Standardization was applied[3,42]:

$$ x'=\frac{x-\mu }{\sigma } $$

where μ is the mean and σ is the standard deviation. This step is crucial for the artificial neural network (ANN) model and distance-based SHapley Additive exPlanations (SHAP) analysis, ensuring that the weight contribution of each physical feature in the loss function is of the same order of magnitude.

Dataset splitting

To objectively evaluate generalization, the cleaned dataset was randomly split into a Training Set (80%) and a Test Set (20%)[5,42]. To ensure reproducibility, the random seed was fixed (Random_State = 42). The training set was used for fitting the physical baseline and updating residual network parameters, while the test set was reserved strictly for final evaluation and was never involved in gradient updates.

Definition of model inputs and output

To ensure a rigorous mapping relationship within the framework, the inputs and output are explicitly defined as follows: Model Inputs (Xaug): A total of 17 physics-augmented descriptors are utilized. These include 14 fundamental chemical and thermodynamic features (e.g., VEC, ΔHmix, ΔSmix, δr, and shear modulus G) calculated via CALPHAD and empirical formulas, alongside 3 critical physics-based terms, with the LSterm (G·δr2/3) being the primary prior. Model Output (Y): The primary and sole direct output of the Hybrid PIML model is the HV. Derived Property: The yield strength (σy) is not a direct output of the ML model but is subsequently derived from the predicted HV based on the phase-dependent physical correlation analysis established in this study.

Descriptor construction

To enable the ML model to truly “comprehend” the physical mechanisms underlying the materials, rather than merely memorizing elemental symbols, this study constructed a multi-dimensional descriptor system encompassing Geometric, Thermodynamic, Electronic Structure, and Physics-Based Strengthening factors. These descriptors not only determine the phase formation rules of HEAs but also serve as core input variables for the SSS model.

Geometric and thermodynamic descriptors

As SSS originates from the lattice elastic stress field induced by solute atoms, the atomic size mismatch (δr) is identified as the primary geometric parameter. In parallel, Mixing Enthalpy (ΔHmix), Mixing Entropy (ΔSmix), and the Ω parameter were selected to assess the thermodynamic stability of solid solution phases against IMCs. The corresponding calculation formulas are detailed below[23,43]:

Atomic size mismatch (δr):

$$ \delta _r=\sqrt{\sum_{i=1}^nc_i(1-\frac{r_i}{\bar{r}})^2} $$

where ci denotes the atomic fraction of element i, ri is the atomic radius, and r = Σciri represents the average atomic radius. The parameter δr directly quantifies the severity of lattice distortion and serves as a core input for the subsequent physics-based strengthening model.

Mixing Enthalpy (ΔHmix) and Mixing Entropy (ΔSmix):

$$ \mathrm{\Delta H_{mix}=\sum_{i < j}4\Delta H_{ij}^{mix}c_{i}c_{j}} $$

$$ \mathrm{\Delta S_{mix}=-R\sum_{i=1}^{n}c_{i}lnc_{i}} $$

where ΔHijmix represents the binary mixing enthalpy between elements i and j, and R is the gas constant.

Thermodynamic Parameter (Ω):

$$ \Omega =\frac{T_m\Delta S_{mix}}{|\Delta H_{mix}|} $$

This parameter integrates the competition between entropic and enthalpic effects and is typically used to determine the stability of solid solutions (where Ω > 1.1 favors solid solution formation)[44,45]. To visualize the coverage and phase stability trends of the dataset, the distribution of the 205 compiled HEA samples in the δr - ΔHmix parameter space is illustrated in Figure 1.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 1. Distribution of 205 HEAs in the δr - ΔHmix parameter space. HEAs: High-entropy alloys; FCC: face-centered cubic; BCC: body-centered cubic.

Electronic descriptors

Experimental studies indicate that the crystal structure of HEAs exhibits a strong dependence on the VEC. Since different lattice structures possess varying numbers of slip systems, their hardness baselines differ vastly (BCC is generally harder than FCC)[3]. Therefore, introducing VEC assists the model in implicitly learning the “Phase Structure-Hardness” correlation[46]:

$$ \mathrm{VEC=\sum_{i=1}^{n}c_i(VEC)_i} $$

Where (VEC)i is the number of valence electrons for the i-th element.

Physics-based strengthening descriptors

This section highlights the key distinction between the proposed Hybrid PIML framework and traditional “composition-property” models. To allow the neural network to explicitly learn SSS laws, this study calculated mechanical modulus parameters directly related to dislocation motion resistance based on the SSS theory:

Average Shear Modulus (G):

$$ \mathrm{G} =\sum_{i=1}^{n}c_iG_i $$

The shear modulus determines the magnitude of dislocation line tension and serves as the “substrate” for measuring material stiffness[31].

Physics-Based Strengthening Factor (LSterm): To provide direct physical guidance to the Hybrid PIML model, a theoretical strengthening term was pre-constructed. Based on the classical Labusch model and its subsequent extensions for multi-principal element alloys (e.g., by Senkov and Toda-Caraballo), the SSS increment is highly dependent on the coupled effect of lattice distortion and modulus mismatch. Accordingly, we explicitly embedded this physical prior by formulating the composite feature:

$$ \mathrm{LS_{term}} =\propto \mathrm{G\cdot \delta _r^{2/3}} $$

This composite feature physically couples lattice distortion (δr) with matrix stiffness (G)[11,38], capturing the theoretical baseline of resistance to dislocation motion. In the subsequent model training, this feature is input into the neural network as the most critical prior knowledge, enabling the model to perform residual correction directly on the basis of physical laws, thereby substantially enhancing the physical consistency of predictions[9].

Feature selection strategy

Although a comprehensive descriptor system was established, an excessive number of input features - particularly those exhibiting high dimensionality and strong correlation - can lead to the ‘curse of dimensionality’ and increase the risk of overfitting. To balance the retention of physical information with the reduction of model complexity, this study implemented a Physics-Guided Hybrid Feature Selection strategy.

Multicollinearity analysis

First, a Pearson correlation coefficient matrix (r) was employed to evaluate the pairwise correlations among all candidate descriptors, as illustrated in Figure 2. In HEAs, certain physical parameters (e.g., atomic radius difference δr and lattice mismatch parameters) may exhibit high mathematical collinearity (|r| > 0.9)[9,22]. The screening principle adopted was as follows: for highly correlated feature pairs, priority was given to retaining descriptors with clearer physical significance and tighter linkage to the SSS theory[23]. For instance, when the simple atomic size difference (δr) was found to be highly correlated with the theoretical strengthening factor (LSterm), the latter was retained because LSterm incorporates shear modulus information, thereby better characterizing the comprehensive resistance to dislocation slip.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 2. Pearson correlation coefficient heatmap of physical descriptors. VEC: Valence electron concentration.

Tree-based recursive feature elimination

Following the removal of redundant features, a recursive feature elimination (RFE) algorithm based on RF was applied[19,43]. Utilizing the “Gini Impurity” metric inherent to tree models, the remaining features were preliminarily ranked by importance. Noise features making negligible contributions to hardness prediction (Importance < 0.01) - such as the atomic number features of certain non-critical elements - were eliminated.

Physics prior retention

Representing the crux of the feature selection strategy, this step introduces domain knowledge intervention, distinguishing the approach from purely data-driven methods[36]. Irrespective of RFE rankings, the retention of three core physical descriptors was enforced for the final model input: LSterm(G·δr2/3) as the theoretical baseline for SSS[38]; VEC as the indicator of phase structural stability; and ΔHmix as the thermodynamic measure of interatomic binding forces. This three-step process reduced the high-dimensional feature space to 17 key descriptors. Consequently, the computational load for ANN training was minimized while ensuring the input information retained clear metallurgical meaning, a prerequisite for the subsequent SHAP mechanistic analysis.

PIML methodology

The PIML framework proposed in this study aims to address the generalization challenges faced by pure data-driven models on “small-sample” HEA data, facilitating the high-throughput development of novel materials[47]. Unlike traditional “black-box” learning, this framework adopts a “theory-guided + data-corrected” hybrid residual learning architecture. This architecture explicitly injects prior knowledge (Inductive Bias) from physical metallurgy into the machine learning process, ensuring that the model’s predictions not only align with physical common sense but also achieve extremely high data-fitting accuracy.

Theoretical baseline: SSS model

As the physical cornerstone of the PIML framework, an analytical baseline model based on classical SSS theory was first constructed. The theoretical foundation originates from the Labusch model[38], which describes the strengthening effect arising from the elastic interaction between dislocations and solute atoms in binary alloys. For multi-principal element systems such as HEAs, recent studies have extended this framework by incorporating the coupled effects of lattice distortion and modulus mismatch[9,36]. Following the approach established by Wen et al.[36] and Huang et al.[9], the theoretical hardness increment (Hphy) can be approximated by a simplified linear expression incorporating the atomic size mismatch (δr) and shear modulus (G):

$$ \mathrm{H_{phy}} =\mathrm{A\cdot G\cdot \delta _r^{2/3}+B} $$

where A and B are fitting constants obtained by linear regression against the experimental hardness data in the training set, representing an average effective-medium approximation of the strengthening behavior. This baseline model thus provides a deterministic first-order prediction. However, the parameters A and B are empirical in nature and inherently embed the simplifying assumptions of a homogeneous, single-phase matrix; they do not capture sample-to-sample deviations arising from phase transitions, short-range ordering, or the differential work-hardening behavior of FCC vs. BCC matrices. The uncertainty introduced by these simplifications is the primary source of prediction error for the baseline model [as seen in its high root mean square error (RMSE) of 212.61 HV in Table 1]. Consequently, the residual term FML in the Hybrid PIML framework is specifically tasked with learning and correcting these complex, physically un-modeled deviations.

Table 1

Performance comparison of different hardness prediction models

Model R2 score RMSE (HV) MAE (HV)
Physical baseline model (Labusch-type) 0.2134 212.61 166.61
RF 0.9257 65.36 38.86
SVR 0.9274 64.60 41.37
XGBoost/LightGBM 0.9452 56.13 37.01
ANN (standard MLP) 0.9702 39.83 26.90
Hybrid PIML (physics-augmented) 0.9762 37.01 25.74

Hybrid architecture: residual learning strategy

To correct the deviations of the physical model, this study introduced a Residual Learning mechanism. In this residual architecture, the primary objective (Output) is the HV. The ANN is specifically tasked with learning the non-linear residual correction term FML - representing the deviation between Hexp and the theoretical baseline Hphy - thereby enabling precise total hardness prediction[34]. The prediction logic of the entire framework is as follows:

$$ \mathrm{H_{final}=H_{phy}+F_{ML}(X_{aug})} $$

Where:
Hphy: The physical baseline value obtained by linear regression fitting of the physics-based baseline formula[9,38]. FML: The non-linear residual correction term constructed by an ANN. Xaug: The physics-augmented input feature vector.

The core advantage of this architecture lies in “divide and conquer” [Figure 3]: allowing the physical formula to handle macroscopic linear laws while enabling the neural network to focus on capturing microscopic, non-linear complex interactions (such as phase structural mutations caused by VEC). This significantly reduces the learning difficulty for the neural network, allowing it to converge quickly and avoid overfitting even with small-sample data[2,27]. Unlike conventional PCNNs that embed empirical physical equations directly into a dual-loss function to regularize network weights[33], our framework adopts a structural “theory-guided + data-corrected” Residual Learning architecture. By explicitly decoupling the macroscopic physical baseline (the Labusch-type model) from the machine learning residual term (FML), the neural network is unburdened from learning the fundamental linear laws from scratch. Instead, it focuses exclusively on capturing the complex non-linear phase-dependent deviations. This structural decoupling not only enhances training efficiency and prevents optimization conflicts common in complex multi-loss formulations but also provides a more transparent mechanism for tracking physical consistency.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 3. Architecture of the proposed Hybrid PIML framework. Hybrid PIML: Hybrid physics-informed machine learning; HEA: high-entropy alloy; VEC: valence electron concentration; ANN: artificial neural network.

Physics-augmented input space

To further enhance the model’s sensitivity to physical laws, the training of the residual network (FML) moved beyond simple raw elemental ratios to construct a physics-augmented input space. Supplementing the fundamental descriptors previously established, this study incorporated the physics-based strengthening term (LSterm = G·δr2/3) as an independent, explicit feature within the neural network. This integration implies that the network not only perceives ‘which elements are present’ but also directly ‘observes’ the theoretically predicted strengthening potential. Such Feature-level Fusion enables the model to keenly identify outliers that violate physical laws and dynamically adjust theoretical values based on other thermodynamic parameters (e.g., automatically lowering hardness predictions to account for work-hardening effects when VEC indicates an FCC phase)[3,39].

Algorithm selection

To construct the optimal PIML framework, this study first conducted extensive benchmarking on various mainstream regression algorithms. The goal was to screen for the algorithm best suited to capture the non-linear mapping relationship between HEA composition and properties as the “Residual Learner” of the framework[7,41]. Five representative classes of models were selected for comparative study: Physical Baseline Model (Labusch-type): Based solely on the classical SSS formula (in linear regression form), serving as a benchmark for evaluating pure theoretical prediction capability. Ensemble Learning Models: Including RF and gradient boosting decision trees [eXtreme Gradient Boosting (XGBoost)/LightGBM], which excel in handling tabular data and feature interactions. Kernel Methods: support vector regression (SVR), utilizing the radial basis function (RBF) kernel to handle non-linear relationships. Deep Learning Models: multi-layer perceptron (ANN/MLP), possessing powerful universal function approximation capabilities. Hybrid PIML framework: The architecture proposed in this study, which explicitly augments the ANN input with the physics-based strengthening term. The specific performance of each model is presented in Table 1. To ensure absolute transparency and mathematical rigor for the results presented in Table 1, a strict machine learning evaluation protocol was utilized. The dataset was randomly partitioned into an 80% training set and a 20% unseen test set using a fixed random seed (Random_State = 42) to guarantee reproducibility. For all ML algorithms (including RF, SVR, XGBoost, and the standard ANN), hyperparameter tuning was conducted utilizing 5-fold cross-validation strictly within the training set to prevent data leakage. The performance metrics reported in Table 1 - namely R2, RMSE, and mean absolute error (MAE) - were calculated exclusively based on the predictions made on the unseen 20% test set. This strict protocol guarantees that the reported values reflect the true generalization and extrapolation capabilities of the models on out-of-sample HEA compositions, rather than mere training data memorization. First, the fundamental limitations of treating HEA SSS as a simple linear problem were unequivocally exposed by the pure physical model (Labusch Model). Its R2 was an abysmal 0.213, with an RMSE as high as 212.61 HV. This catastrophic empirical failure directly refutes the assumption that HEA strengthening can be resolved by simple linear combinations of lattice distortion and shear modulus. It demonstrates that the hardness mechanism of HEAs is governed by profound non-linearities, such as structural phase transitions induced by VEC and complex strain-hardening behaviors. This directly motivates the absolute necessity of introducing deep data-driven methods[9,38]. By utilizing the ANN as a residual learner to complement the linear physical baseline[27,48], the Hybrid PIML model successfully captured these deep non-linear laws, jumping to an R2 of 0.976 and dropping the RMSE to 37.01 HV.

RESULTS AND DISCUSSION

Model performance benchmarking

This paper subjected all models to rigorous testing using the same training and test sets (80:20 split, Random_State = 42). The quantitative performance metrics for each model are summarized in Table 1, while the corresponding parity plots comparing predicted versus experimental hardness are illustrated in Figure 4. The Hybrid PIML model achieved the best performance among all models by introducing physics-enhanced features.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 4. Comparison of experimental vs. predicted HV across different models. (A) Physical baseline model, (B) RF, (C) SVR, (D) XGBoost, (E) ANN and (F) the proposed Hybrid PIML model. The dashed line indicates perfect prediction (y = x). The Hybrid PIML model demonstrates significantly tighter scattering around the diagonal, validating its superior accuracy. HV: Vickers hardness; RF: random forest; SVR: support vector regression; XGBoost: eXtreme Gradient Boosting; ANN: artificial neural network; Hybrid PIML: hybrid physics-informed machine learning; R2: coefficient of determination; RMSE: root mean square error; MLP: multi-layer perceptron.

The Hybrid PIML model reduced the RMSE to 37.01 HV, significantly outperforming both the pure theoretical model and pure ML models. This confirms the effectiveness of the residual learning strategy: the physics-based baseline model handles the primary SSS contribution, greatly reducing the learning difficulty for the ML model, allowing it to focus solely on learning the complex non-linear effects not covered by the theory (such as phase mutation). An R2 approaching 0.98 indicates that the framework explains almost all sources of variance in the dataset.

Mechanistic interpretability: SHAP analysis

To open the “black box” of the machine learning model and verify whether it has learned hardness mechanisms consistent with physical metallurgy principles, this study employed the SHAP method to conduct an in-depth interpretability analysis of the XGBoost, Standard ANN, and Hybrid PIML models. Figures 5-7 display the SHAP Summary Plots, SHAP Feature Importance Bar Plots, and Univariate Dependence Plots for the three models on the test set.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 5. Mechanistic interpretability analysis of the baseline XGBoost model. (A) SHAP summary plot and (B) feature importance ranking derived from the standard XGBoost algorithm; (C) SHAP dependence plot for the top-ranked feature. XGBoost: eXtreme Gradient Boosting; SHAP: SHapley Additive exPlanations; VEC: valence electron concentration.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 6. Baseline interpretability analysis of the standard ANN model. (A) SHAP summary plot illustrating the feature impact distribution for the pure data-driven ANN; (B) Global feature importance ranking. While the model identifies key descriptors, the ranking relies solely on statistical correlations within the training data; (C) SHAP dependence plot for the top-ranked feature. ANN: Artificial neural network; SHAP: SHapley Additive exPlanations; VEC: valence electron concentration.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 7. Mechanistic interpretability analysis of the Hybrid PIML framework. (A) SHAP summary beeswarm plot. The physics-derived strengthening parameter (LSterm) exhibits a distinct positive correlation with hardness, validating the model’s adherence to SSS principles; (B) Global feature importance ranking. The LSterm is identified as the dominant predictor among all features, confirming that the PIML architecture successfully prioritizes physical mechanisms over raw compositional data; (C) SHAP dependence plot for LSterm. The curve reveals a clear monotonic increasing relationship between the theoretical strengthening factor and predicted hardness, demonstrating the model’s capability to capture the underlying non-linear physical laws. Hybrid PIML: Hybrid physics-informed machine learning; SHAP: SHapley Additive exPlanations; SSS: solid solution strengthening; VEC: valence electron concentration.

Limitations of traditional ML models (XGBoost and standard ANN)

As seen in Figures 5A and 6A, while these two pure data-driven models can identify some key features, their decision logic relies primarily on statistical correlation rather than physical causality. Spurious Statistical Correlation: In the XGBoost model, atomic number or the content of specific elements (e.g., Co or Ni) are often assigned excessively high weights. Although these elements are associated with phase structure to some extent, they are not direct determinants of hardness in SSS theory. Ambiguity of Physical Mechanisms: Although the standard ANN model captured the importance of VEC, correctly reflecting VEC’s control over BCC/FCC phase stability, it lacked sensitivity to the interaction between lattice distortion (δr) and shear modulus (G). Consequently, when predicting high-distortion alloys (such as refractory HEAs), the model often infers results by “memorizing” composition data rather than “understanding” the lattice stress field, limiting its extrapolation capability in unknown compositional spaces.

Physical consistency of the Hybrid PIML model

Conversely, the Hybrid PIML model, augmented with the physics-based feature (LSterm = G·δr2/3), demonstrates superior physical consistency. SHAP analysis highlights the following mechanisms: (1) Dominance of Physical Descriptors: The engineered LS_term dominates the feature importance ranking [Figure 7B], surpassing individual elemental features. This confirms a shift in the model’s attention from statistical composition to physical fundamentals, identifying the coupling of shear modulus (G) and atomic size mismatch (δr) as the core determinant of HEA hardness; (2) Theoretically Aligned Strengthening Direction: Figure 7A shows that high LS_term values correlate positively with SHAP values, indicating that increased lattice distortion and modulus stiffness lead to higher predicted hardness. This behavior is consistent with classical SSS theory (H∝G·δr2/3), validating that the model has internalized physical causality rather than relying solely on data fitting; (3) Non-linear Threshold Effect: Figure 7C illustrates a non-linear contribution of the physical term to hardness. While the relationship is linear in low-distortion regions, a slope change is observed in high-distortion regions (e.g., refractory HEAs). This deviation implies complex coupling between SSS and phenomena such as SRO or nano-cluster strengthening under high lattice stress - non-linear details effectively captured by the PIML framework.

Summary

In summary, the comparative analysis shows that while pure data-driven models (XGBoost/ANN) perform acceptably in numerical precision, they suffer from the defect of “knowing the what but not the why” regarding physical interpretability. The Hybrid PIML framework, however, not only achieved the highest prediction accuracy (R2 = 0.976) but, more critically, established a decision-making process consistent with physical metallurgy principles. This “physics-aware” capability makes the PIML model more reliable in new material discovery, effectively avoiding erroneous predictions that violate thermodynamic common sense.

Correlation analysis: yield strength vs. hardness

In classical physical metallurgy, the relationship between yield strength (σy) and HV is a long-established, fundamental concept. Derived from the Tabor relation, the plastic constraint developed beneath the indenter yields a well-known approximate empirical scaling of σy (MPa) ≈ C × HV, where C is a universally recognized conversion factor typically ranging from 2.5 to 3.5. We do not claim this classical scaling as a new finding; rather, we utilize this foundational metallurgical literature as a physical baseline.

The novelty of our correlation analysis lies in revealing that this widely accepted classical approximation critically fails in complex HEA systems. Statistical analysis of the 205 HEA samples in this study explicitly demonstrates that a universal empirical constant C cannot be applied across the broad HEA compositional space. Instead, the hardness-strength scaling exhibits a profound Phase-Structure Dependency that requires machine learning to accurately capture and correct.

Figure 8 illustrates the correlation between HV and yield strength for HEAs with distinct phase structures. By segregating the data into BCC, FCC, and multi-phase categories, significant mechanistic divergences are observed: BCC Alloys (Quasi-linear Behavior): As depicted by the red data points, BCC HEAs (e.g., W-Mo-Ta-Nb systems) closely follow the traditional approximation line of σy ≈ 3 × HV (Pearson r ≈ 0.88). This alignment indicates that BCC alloys, characterized by high Peierls stresses and limited work-hardening capacity, behave as near-ideal elastic-plastic materials where the yield strength is the dominant factor controlling indentation resistance. FCC Alloys (Work-Hardening Deviation): Conversely, FCC alloys (blue data points) exhibit a phenomenon of “pseudo-high hardness”, falling significantly below the theoretical reference line of σy = 3HV (i.e., for a given hardness, the actual yield strength is lower than predicted, often σy ≈ 1.5 × HV). This deviation is attributed to the exceptional work-hardening capability of FCC HEAs. The Vickers indentation process induces substantial plastic strain (~8%), causing the material to work-harden significantly (σflow >> σy). Consequently, the measured hardness reflects the flow stress at high strain rather than the initial yield strength, leading to an overestimation of σy if the standard coefficient of 3 is used. Adaptive Correction by PIML: The PIML model’s ability to capture this non-linear “pseudo-high hardness” in FCC alloys highlights the critical necessity of ML beyond simple linear solid-solution models. Just as interpretable ML has revealed counterintuitive, non-linear relationships between pore irregularity and yield strength - challenging the assumption that all defects are universally detrimental - our framework successfully identifies fundamental deviations from the standard macroscopic scaling law (σy ≈ 3 × HV). Furthermore, regarding the integration of microstructural features, while recent advancements have successfully incorporated explicit microstructural parameters (e.g., grain size, boundary misorientation) into PCNNs for localized corrosion prediction via inverse parameter identification[33], acquiring such exhaustive experimental microstructural data across vast, uncharted HEA compositional spaces remains experimentally prohibitive for high-throughput screening. Our framework elegantly circumvents this experimental bottleneck. By utilizing physically informed thermodynamic and electronic descriptors (e.g., VEC, δr) as robust, computable proxies for phase stability and microstructural tendencies, the model enables rapid and reliable virtual screening across tens of thousands of candidates under data-constrained conditions. The observed discrepancy underscores that a single conversion coefficient is invalid for all HEAs. The proposed Hybrid PIML framework addresses this by incorporating phase-stability descriptors as inputs. This allows the model to implicitly classify the alloy system and dynamically adjust the hardness-strength mapping - applying a lower coefficient for work-hardenable FCC alloys and a standard coefficient for BCC alloys - thereby achieving precise predictions across the entire compositional space.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 8. Phase-dependent mapping relationship between hardness and yield strength. FCC: Face-centered cubic; BCC: body-centered cubic; HV: Vickers hardness.

High-throughput virtual screening strategy based on PIML

The ultimate objective of constructing the Hybrid PIML framework in this study is to address the ‘curse of dimensionality’ in HEA design. Leveraging the model’s superior performance in both accuracy and physical interpretability, a high-throughput virtual screening (HTVS) strategy was proposed to guide the accelerated discovery of potential superhard HEAs.

Workflow design

In the virtual screening funnel, the ‘optimization performance’ is defined by two hierarchical criteria: Maximization of HV (Target Output): Identifying compositions with predicted HV > 900 HV; Phase Stability Constraint: Filtering for BCC-dominant solid solution structures while maintaining a balance between high hardness and the avoidance of brittle intermetallic phases. Traditional experimental trial-and-error methods are often cost-prohibitive and limited to the central regions of phase diagrams. Leveraging the Hybrid PIML model, a Screening Funnel was established to navigate the expanded design space, as illustrated in Figure 9. First, the target system (e.g., the refractory W-Mo-Ta-Nb-V-Ti system) is defined with a step size of 5 at.%, generating approximately 53,000 potential non-equiatomic combinations. Next, physics-based physical descriptors (e.g., δr, ΔHmix, VEC) are calculated for each virtual composition. The trained Hybrid PIML model is then deployed to predict the hardness of all candidates, ranking them by predicted value. Finally, phase formation laws implicitly learned by the model (as shown in the SHAP analysis) are applied to filter out compositions with extremely low predicted hardness - implying the formation of soft FCC phases - or abnormal physical parameters that suggest a risk of amorphous formation.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 9. HTVS workflow based on the Hybrid PIML framework. HTVS: High-throughput virtual screening; Hybrid PIML: hybrid physics-informed machine learning; VEC: valence electron concentration; FCC: face-centered cubic.

Case study and validation: compositional optimization of refractory HEAs

To validate the screening strategy, the Co-Cr-Ti-Mo-W system was subjected to retrospective analysis. While the traditional equiatomic design (Co20Cr20Ti20Mo20W20) exhibits limited performance, the Hybrid PIML model successfully identified optimized non-equiatomic regions [Figure 10]. The analysis suggests that reducing Ti/W content and increasing the Co/Cr matrix fraction optimizes the trade-off between SSS and brittle phase suppression, specifically within the range of (Co + Cr + Mo) ≈ 0.9 and (W + Ti) ≈ 0.1. These predictions are strongly supported by recent experimental literature[11,41,49]. Evidence shows that alloys in this optimized interval surpass the equiatomic benchmark, maintaining ultra-high hardness (> 900 HV) while mitigating brittle intermetallic phases (e.g., σ or Laves phase) precipitation induced by severe lattice distortion[50]. This finding accords with Wen et al. on the synergistic enhancement of strength and ductility via non-equiatomic design[51]. Ultimately, this case study demonstrates that the Hybrid PIML model offers both robust interpolation accuracy and the extrapolation capability required to navigate complex compositional spaces and locate local optima.

Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

Figure 10. Compositional optimization heatmap of the Co-Cr-Mo-W-Ti system predicted by the Hybrid PIML model. Hybrid PIML: Hybrid physics-informed machine learning; HV: Vickers hardness.

To conclusively address the critical necessity for experimental validation regarding phase constitutions (e.g., confirming the solid-solution-only condition) and macroscopic strengthening mechanisms, we cross-verified our PIML-optimized compositions against extensive experimental efforts documented in recently published, independent literature.

To address the critical necessity for experimental validation regarding both macroscopic hardness and microscopic phase constitutions, we cross-verified our PIML-optimized compositions against recent, independent high-throughput experimental efforts on the exact CoxCryTizMouWv system[52]. While direct in-house synthesis and microstructural characterizations [scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS)] are beyond the scope of this computational methodology study, the rigorous experimental data provided by Liu et al. offers robust out-of-sample validation[52]. Crucially, their experimental microstructural analysis [X-ray diffraction (XRD)/SEM] of equimolar CoCrTiMoW revealed the precipitation of Co-Cr-Ti IMCs alongside the BCC solid-solution dendrites due to excessive Ti content. This experimental finding perfectly justifies the optimization logic internalized by our Hybrid PIML framework. As revealed by our SHAP analysis, the model heavily relies on the solid-solution strengthening term (LSterm). To maximize this solid-solution contribution and suppress brittle second-phase precipitations, our model correctly predicted that the atomic ratios of severe lattice distorters must be significantly constrained [i.e., (W + Ti) ≈ 0.1 with an elevated Co/Cr matrix]. Furthermore, Liu et al. experimentally confirmed that non-equiatomic alloys designed following similar compositional constraints successfully maintain structural integrity and exhibit ultra-high HV exceeding 900 HV (up to 971 HV)[52]. This exceptional agreement between independent experimental observations and our mathematical outputs not only validates the predicted hardness values but also explicitly proves that the Hybrid PIML framework successfully captured the true physical metallurgy principles governing solid-solution optimization in complex HEA design spaces. Future work will entail targeted in-house synthesis to further map the phase boundaries of these optimized candidates.

CONCLUSIONS

To address the “curse of dimensionality” and data scarcity challenges in the design of HEAs, this study developed and validated a Hybrid PIML framework integrating classical SSS theory. By synergizing a physics-based baseline with a residual learning neural network, the framework achieves high-precision hardness prediction while successfully internalizing fundamental metallurgical principles such as dislocation motion resistance. Furthermore, the model accurately captures the phase-dependent non-linear mapping between yield strength and hardness, demonstrating superior physical consistency and generalization across diverse crystallographic structures. The high-throughput screening strategy and design methodology established herein provide an efficient and universal computational platform for accelerating the discovery of superhard, high-performance structural materials within expansive compositional spaces.

DECLARATIONS

Authors’ contributions

Investigation, data curation and analysis, conceptualization, figure preparation, writing - original draft: Gao, A.

Data processing, figure preparation: Yan, Y.

Supervision, validation, funding acquisition, methodology, writing - review and editing: Liu, X.; Chong, X.

Discussion; project administration, writing - review and editing: Tao, Q.; Yu, J.; Xi, S.; Li, G.

Availability of data and materials

Raw data that support the findings are available from the corresponding author upon reasonable request.

AI and AI-assisted tools statement

Not applicable.

Financial support and sponsorship

The authors would like to thank the financial support from the Major R&D Project of Yunnan Province (Nos. 202402AB080007).

Conflicts of interest

Liu, X. is the Executive Editor-in-Chief, and Chong, X. is a Editorial Board Member of the journal Journal of Materials Informatics, but were not involved in any steps of editorial processing, notably including reviewer selection, manuscript handling, and decision making, while the other authors have declared that they have no conflicts of interest.

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Copyright

© The Author(s) 2026.

REFERENCES

1. Yeh, J.; Chen, S.; Lin, S.; et al. Nanostructured high‐entropy alloys with multiple principal elements: novel alloy design concepts and outcomes. Adv. Eng. Mater. 2004, 6, 299-303.

2. Liu, C.; Meng, M.; Luo, X. Composition design and property prediction for AlCoCrCuFeNi high-entropy alloy based on machine learning. Metals 2025, 15, 733.

3. Chen, D.; Guo, J.; Guo, L.; et al. A phase prediction strategy of high entropy alloys manufactured by selective laser melting based on feature selection and machine learning. J. Mater. Res. Technol. 2025, 37, 79-88.

4. Kandavalli, M.; Agarwal, A.; Poonia, A.; Kishor, M.; Ayyagari, K. P. R. Design of high bulk moduli high entropy alloys using machine learning. Sci. Rep. 2023, 13, 20504.

5. Arun, S.; Radhika, N. Phase formation in vacuum arc melted high entropy alloys: insights from machine learning and experimental validation. IEEE. Access. 2025, 13, 146772-85.

6. Zhao, S.; Jiang, B.; Song, K.; et al. Machine learning assisted design of high-entropy alloys with ultra-high microhardness and unexpected low density. Mater. Design. 2024, 238, 112634.

7. Paturi, U. M. R.; Ishtiaq, M.; Lakshmi Narayana, P.; Maurya, A. K.; Choi, S.; Reddy, N. G. S. Evaluating machine learning models for predicting hardness of AlCoCrCuFeNi high-entropy alloys. Crystals 2025, 15, 404.

8. Ding, S.; Zhang, Y.; Lei, S.; et al. Data-driven machine learning with lattice distortion and thermodynamic parameters guided strength optimization of refractory high-entropy alloys. J. Mater. Res. Technol. 2025, 38, 133-42.

9. Huang, X.; Jin, C.; Zhang, C.; Zhang, H.; Fu, H. Machine learning assisted modelling and design of solid solution hardened high entropy alloys. Mater. Design. 2021, 211, 110177.

10. Kim, G.; Diao, H.; Lee, C.; et al. First-principles and machine learning predictions of elasticity in severely lattice-distorted high-entropy alloys with experimental validation. Acta. Mater. 2019, 181, 124-38.

11. Shen, L.; Li, Y.; Zhang, W.; et al. Machine learning-assisted design of strong and ductile BCC high-entropy alloys. Mater. Res. Lett. 2025, 13, 1260-8.

12. Li, J.; Xiong, W.; Zhang, T.; et al. Machine learning and explainable AI-guided design and optimization of high-entropy alloys as binder phases for WC-based cemented carbides. Comput. Mater. Contin. 2025, 84, 2189-216.

13. Ghouchan Nezhad Noor Nia, R.; Jalali, M.; Mail, M.; Ivanisenko, Y.; Kübel, C. Machine learning approach to community detection in a high-entropy alloy interaction network. ACS. Omega. 2022, 7, 12978-92.

14. Wu, Y.; Yan, G.; Yu, W.; Shen, S. Investigating nanostructure-property relationship of WTaVCr high-entropy alloy via machine learning optimized reactive potential. J. Mater. Res. Technol. 2024, 32, 2624-37.

15. Mazitov, A.; Springer, M. A.; Lopanitsyna, N.; Fraux, G.; De, S.; Ceriotti, M. Surface segregation in high-entropy alloys from alchemical machine learning. J. Phys. Mater. 2024, 7, 025007.

16. Sato, K.; Hayashi, G.; Ogushi, K.; et al. Computational materials design of high-entropy alloys based on full potential Korringa-Kohn-Rostoker coherent potential approximation and machine learning techniques. Mater. Trans. 2023, 64, 2174-8.

17. Giles, S. A.; Sengupta, D.; Broderick, S. R.; Rajan, K. Machine-learning-based intelligent framework for discovering refractory high-entropy alloys with improved high-temperature yield strength. npj. Comput. Mater. 2022, 8, 926.

18. Zhang, L.; Qian, K.; Schuller, B. W.; Shibuta, Y. Prediction on mechanical properties of non-equiatomic high-entropy alloy by atomistic simulation and machine learning. Metals 2021, 11, 922.

19. Zhong, K.; Wu, S.; Liu, S.; Qi, J.; Cheng, L.; Ai, D. Machine learning-driven optimization of hardness and toughness in high-entropy alloy coatings based on composition and descriptor. J. Appl. Phys. 2025, 137, 195305.

20. Elgack, O.; Almomani, B.; Syarif, J.; Elazab, M.; Irshaid, M.; Al-Shabi, M. Molecular dynamics simulation and machine learning-based analysis for predicting tensile properties of high-entropy FeNiCrCoCu alloys. J. Mater. Res. Technol. 2023, 25, 5575-85.

21. Gupta, K. K.; Barman, S.; Dey, S.; Mukhopadhyay, T. Explainable machine learning assisted molecular-level insights for enhanced specific stiffness exploiting the large compositional space of AlCoCrFeNi high entropy alloys. Mach. Learn. Sci. Technol. 2024, 5, 025082.

22. He, S.; Wang, Y.; Zhang, Z.; et al. Interpretable machine learning workflow for evaluation of the transformation temperatures of TiZrHfNiCoCu high entropy shape memory alloys. Mater. Design. 2023, 225, 111513.

23. Berry, J.; Snell, R.; Anderson, M.; et al. Design and selection of high entropy alloys for hardmetal matrix applications using a coupled machine learning and calculation of phase diagrams methodology. Adv. Eng. Mater. 2024, 26, 2302064.

24. Klimenko, D.; Stepanov, N.; Li, J.; Fang, Q.; Zherebtsov, S. Machine learning-based strength prediction for refractory high-entropy alloys of the Al-Cr-Nb-Ti-V-Zr system. Materials 2021, 14, 7213.

25. Wu, M.; Guan, B.; Wang, J.; et al. Accelerated design of eutectic high-entropy alloys using high-throughput phase diagram calculations and machine learning. Mater. Design. 2025, 254, 114125.

26. Abdul Salam, M. Y.; Ogunmuyiwa, E. N.; Manisa, V. K.; Yahya, A.; Badruddin, I. A. Enhancing phase characterization of AlCuCrFeNi high entropy alloys using hybrid machine learning models: a comprehensive XRD analysis. J. Mater. Res. Technol. 2025, 36, 592-605.

27. Akhlaghi, Z. A.; Shahmir, H.; Sharafat, A. R. Machine learning for obtaining values of thermomechanical processing parameters of high-entropy alloys with desirable strength and ductility. Intermetallics 2025, 185, 108892.

28. Liu, S.; Bocklund, B.; Diffenderfer, J.; et al. A comparative study of predicting high entropy alloy phase fractions with traditional machine learning and deep neural networks. npj. Comput. Mater. 2024, 10, 1335.

29. Thampiriyanon, J.; Khumkoa, S. Machine learning–based prediction of complex combination phases in high-entropy alloys. Metals 2025, 15, 227.

30. He, J.; Li, Z.; Zhao, P.; et al. Machine learning-assisted design of high-entropy alloys with superior mechanical properties. J. Mater. Res. Technol. 2024, 33, 260-86.

31. Zheng, S.; Yang, L.; Fang, L.; et al. Exploration of ductility for refractory high entropy alloys via interpretive machine learning. J. Mater. Res. Technol. 2025, 37, 1243-56.

32. Wang, J.; Lee, S.; Kim, Y. W.; Kim, K. T.; Park, J. M.; Raabe, D. Data-driven analysis relates mechanical properties to pore morphology in laser powder bed fusion. Acta. Mater. 2026, 304, 121751.

33. Rahman, M. W.; Zhao, B.; Pan, S.; Qu, Y. Microstructure-informed machine learning for understanding corrosion resistance in structural alloys through fusion with experimental studies. Comput. Mater. Sci. 2025, 248, 113624.

34. Liu, S.; Lee, K.; Balachandran, P. V. Integrating machine learning with mechanistic models for predicting the yield strength of high entropy alloys. J. Appl. Phys. 2022, 132, 105105.

35. Thiercelin, L.; Peltier, L.; Meraghni, F. Physics-informed machine learning prediction of the martensitic transformation temperature for the design of “NiTi-like” high entropy shape memory alloys. Comput. Mater. Sci. 2024, 231, 112578.

36. Wen, C.; Zhang, Y.; Wang, C.; et al. Machine learning assisted design of high entropy alloys with desired property. Acta. Mater. 2019, 170, 109-17.

37. Li, S.; Li, S.; Liu, D.; Zou, R.; Yang, Z. Hardness prediction of high entropy alloys with machine learning and material descriptors selection by improved genetic algorithm. Comput. Mater. Sci. 2022, 205, 111185.

38. Ren, W.; Zhang, Y.; Wang, W.; Ding, S.; Li, N. Prediction and design of high hardness high entropy alloy through machine learning. Mater. Design. 2023, 235, 112454.

39. He, J.; Li, Z.; Lin, J.; et al. Machine learning-assisted design of refractory high-entropy alloys with targeted yield strength and fracture strain. Mater. Design. 2024, 246, 113326.

40. Bhandari, U.; Rafi, M. R.; Zhang, C.; Yang, S. Yield strength prediction of high-entropy alloys using machine learning. Mater. Today. Commun. 2021, 26, 101871.

41. Pan, Y.; Xu, B.; Wang, J.; et al. A machine learning approach to predict coefficient of friction of laser cladding high-entropy alloy coatings. J. Mater. Res. Technol. 2025, 36, 8240-53.

42. Sivaraman, S.; Radhika, N.; Abubaker Khan, M. Machine learning-driven prediction of wear rate and phase formation in high entropy alloy coatings for enhanced durability and performance. IEEE. Access. 2025, 13, 33956-75.

43. Otieno, R.; Odhong, E. V.; Ondieki, C. Prediction of phases and mechanical properties of magnesium-based high-entropy alloys using machine learning. J. King. Saud. Univ. Sci. 2024, 36, 103456.

44. Hou, S.; Li, Y.; Bai, M.; et al. Phase prediction of high-entropy alloys by integrating criterion and machine learning recommendation method. Materials 2022, 15, 3321.

45. Halpren, E.; Yao, X.; Chen, Z. W.; Singh, C. V. Machine learning assisted design of BCC high entropy alloys for room temperature hydrogen storage. Acta. Mater. 2024, 270, 119841.

46. Singh, S.; Katiyar, N. K.; Goel, S.; Joshi, S. N. Phase prediction and experimental realisation of a new high entropy alloy using machine learning. Sci. Rep. 2023, 13, 4811.

47. Mooraj, S.; Chen, W. A review on high-throughput development of high-entropy alloys by combinatorial methods. J. Mater. Inf. 2023, 3, 4.

48. Veeresham, M.; Sake, N.; Lee, U.; Park, N. Unraveling phase prediction in high entropy alloys: a synergy of machine learning, deep learning, and ThermoCalc, validation by experimental analysis. J. Mater. Res. Technol. 2024, 29, 1744-55.

49. Wang, S.; Li, D.; Xiong, J. Prediction of elastic properties of face-centered cubic high-entropy alloys by machine learning. Trans. Nonferrous. Met. Soc. China. 2023, 33, 518-30.

50. Salam, M. Y. A.; Ogunmuyiwa, E. N.; Manisa, V. K.; Yahya, A.; Khader, A.; Badruddin, I. A. Investigation of heat treatment effects on microstructure and mechanical properties of AlCuCrFeNiTi high-entropy alloys employing hybrid model and feature-adaptive machine learning. Case. Stud. Therm. Eng. 2025, 76, 107325.

51. Wen, C.; Zhang, Y.; Wang, C.; et al. Machine-learning-assisted compositional design of refractory high-entropy alloys with optimal strength and ductility. Engineering 2025, 46, 214-23.

52. Liu, Y.; Wang, J.; Xiao, B.; Shu, J. Accelerated development of hard high-entropy alloys with data-driven high-throughput experiments. J. Mater. Inf. 2022, 2, 3.

Cite This Article

Research Article
Open Access
Physics-informed machine learning framework integrating solid solution strengthening theory for accelerated hardness prediction in high-entropy alloys

How to Cite

Download Citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click on download.

Export Citation File

Type of Import

Tips on Downloading Citation

This feature enables you to download the bibliographic information (also called citation data, header data, or metadata) for the articles on our site.

Citation Manager File Format

Use the radio buttons to choose how to format the bibliographic data you're harvesting. Several citation manager formats are available, including EndNote and BibTex.

Type of Import

If you have citation management software installed on your computer your Web browser should be able to import metadata directly into your reference database.

Direct Import: When the Direct Import option is selected (the default state), a dialogue box will give you the option to Save or Open the downloaded citation data. Choosing Open will either launch your citation manager or give you a choice of applications with which to use the metadata. The Save option saves the file locally for later use.

Indirect Import: When the Indirect Import option is selected, the metadata is displayed and may be copied and pasted as needed.

Data & Comments

Data

Views
327
Downloads
19
Citations
0
Comments
0
0

Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at [email protected].

Journal of Materials Informatics
ISSN 2770-372X (Online)
Follow Us

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/