Nitrogen alloying improves the mechanical performance of martensitic stainless steel, while tempering is required to mitigate brittleness and enhance processability. However, tempering-induced microstructural changes markedly influence the semiconducting properties of the passive film, thereby affecting corrosion resistance. This work examines the passive film of high-nitrogen martensitic stainless steel (HNMSS) tempered at different temperatures. Tempering at 200 °C yields the highest pitting potential (461.28 mV), attributed to the transformation of CrN to Cr₂O₃, which enhances film protection. At higher temperatures, abundant Cr-rich nitride/precipitate (M₂N) depletes N and Cr in the matrix. The increase of precipitate/matrix interfaces and defect density in the passivation film impairs pitting corrosion resistance. First-principles calculations and quasi-in-situ scanning Kelvin probe force microscopy reveal that M₂N precipitates exhibit a higher Volta potential and the highest work function than those of the matrix, acting as cathodes to accelerate localized matrix dissolution, which reduces nitrogen incorporation into the passive film. These findings clarify the relationship of tempering, microstructure, and corrosion in the HNMSS.

Universal machine-learning interatomic potentials (uMLIPs) are emerging as foundation models for atomistic simulation, offering near-ab initio accuracy at far lower cost. Their safe, broad deployment is limited by the absence of reliable, general uncertainty estimates. We present a unified, scalable uncertainty metric, U, built from a heterogeneous ensemble that reuses existing pretrained MLIPs. Across diverse chemistries and structures, U strongly tracks true prediction errors and robustly ranks configuration-level risk. Using U, we perform uncertainty-aware distillation to train system-specific potentials with far fewer labels: for tungsten, we match full density-functional-theory (DFT) training using 4% of the DFT data; for MoNbTaW, a dataset distilled by U supports high-accuracy potential training. By filtering numerical label noise, the distilled models can in some cases exceed the accuracy of the MLIPs trained on DFT data. This framework provides a practical reliability monitor and guides data selection and fine-tuning, enabling cost-efficient, accurate, and safer deployment of foundation models.

Carbide particles act as effective hydrogen traps that can mitigate hydrogen embrittlement in steels, with their trapping efficacy governed by the binding energy of hydrogen at atomic-scale sites within the carbide and at the carbide/matrix interface. Within the carbide, carbon vacancies are established trap sites, while at the carbide/matrix interface, multiple distinct trap sites may form depending upon the interface character. State-of-the-art experimental methods cannot attribute binding energies to individual sites unambiguously, particularly at the carbide/matrix interface, where distinct trap sites may exhibit similar binding energies. Density functional theory (DFT) calculations have partially addressed this limitation, but their computational cost prohibits accurate treatment of semicoherent and incoherent interfaces, for which hydrogen binding energies remain poorly characterised. In this work, we develop a DFT-accurate machine-learning Fe–TiC–H interatomic potential that enables simulations of hydrogen energetics and kinetics at various TiC/Fe interfaces in the Baker–Nutting orientation. Predictions of structural, dynamic and energetic properties of both bulk phases (i.e., bcc Fe and rocksalt TiC) and the coherent interface with and without hydrogen are in good agreement with DFT. We apply the potential to examine the atomic structure and hydrogen trapping at semicoherent and incoherent interfaces. We find that the binding energy of hydrogen at misfit dislocations is similar to that at coherent interfaces, with a difference of no more than 0.1 eV. However, the misfit strain field reduces carbon vacancy formation energies at the interface, where hydrogen can be trapped with similar binding energies as that in a bulk C vacancy (approximately –1 eV). Lower vacancy formation energies are also found at the incoherent interface, where hydrogen binding energies span a broader range, depending sensitively on the local atomic environment. Overall, the results enable interpretation of experimental observations of hydrogen segregation around semicoherent and incoherent TiC interfaces.

Enhancing the hydrogen embrittlement (HE) resistance of alloys caters to the urgent needs of engineering safety and long-distance hydrogen transportation. Highly dense precipitates in the alloys act as H traps, however, some of them cannot strongly trap H thus failing to prevent its accumulation at the critical regions. Experimentally, it is challenging to expeditiously identify and generate phases causing strengthening and acting as strong H traps. Here, we demonstrate a computation-based design strategy to generate precipitates strongly trapping H. Based on the quantum machine learning Al-Sc-Cu potential, the optimal processing parameters of strong H trapping phases are determined, even though they are metastable in nature. Elemental mapping in electron microscope and atom probe tomography confirms the presence of Cu in Al₃Sc and its strong interaction with H. Hence, we envisage the proposed strategy will accelerate the design of HE-resistant microstructures of various technologically relevant materials via identification of desirable phases.

Iron oxidation is a complex process involving critical atomistic events, such as atomic adsorption, diffusion, and surface reconstruction, understanding of which is significant for both surface science and coating technology. Atomistic simulation serves as an useful tool to investigate the processes, where description of interatomic interactions is required. However, selecting appropriate force field or interatomic potential is not only difficult, but also essential for getting accurate result. In this work, we present a detailed benchmark of reactive force fields (ReaxFFs) and universal machine learning interatomic potentials (uMLIPs) against density functional theory (DFT) calculations of oxygen adsorption on various α-iron surfaces, which is the first yet crucial step towards oxidation. The comparisons show the coverage-dependent performance and improvable accuracy of both ReaxFFs and uMLIPs at reproducing DFT results, with ReaxFFs outperforming uMLIPs. Subsequently, iron oxidation is simulated using ReaxFF and uMLIP. The results reveal the strong capability of ReaxFF and poor stability of uMLIP for describing reactive process, i.e., the formation of iron oxide. This may be attributed to the suitable functional form of ReaxFF for the description of bond changes. The insights presented here not only provide an example of benchmarking force field or interatomic potential for system of interest, but also highlight the applicability of ReaxFF and scopes of improvement of uMLIP.

Machine-learning interatomic potentials (MLIPs) with broad chemical flexibility are essential for atomistic simulations of compositionally complex alloys, but their deployment in large-scale molecular dynamics requires a balance among accuracy, efficiency, stability, transferability, and uncertainty quantification. Here, we benchmark two chemically scalable MLIP frameworks, neuroevolution potential (NEP) and graph atomic cluster expansion (GRACE), for 16 elemental metals and their multicomponent alloys. GRACE-FS shows higher training efficiency and generally better average accuracy, chemical transferability, and finite-temperature robustness, whereas UNEP-v1 provides substantially higher inference speed and remains competitive in selected stress and large-error metrics. We further show that chemical transferability is closely linked to high-temperature MD stability in highly multicomponent environments and that ensemble-based uncertainty provides a more reliable error indicator than D-optimality for the heterogeneous systems considered here. Finally, 3×106-atom shock simulations demonstrate that UNEP-v1, combined with ensemble uncertainty, enables uncertainty-aware simulations under extreme dynamic conditions, yielding robust global spall-strength predictions while revealing model sensitivity in local damage pathways. These results provide practical guidelines for selecting and deploying MLIPs in large-scale simulations of multicomponent alloys.

Multi-principal element alloys (MPEAs) are renowned for their enhanced mechanical strength relative to their constituent metals, as evidenced by various experimental techniques such as tension/compression tests and instrumental indentation. Nevertheless, atomistic simulations sometimes produce conflicting results, casting doubt on the consistently superior mechanical properties of MPEAs. In this study, machine-learning interatomic potentials (MLIPs) with first-principles accuracy were developed for body-centered cubic refractory MoNbTaW MPEAs, enabling systematic atomistic simulations under various deformation scenarios. The new MLIPs are supported by a comprehensive dataset encompassing extensive defects, and the established embedded-atom model (EAM) potential was benchmarked against both this dataset and the new MLIP. Simulations covering diverse compositions confirm that both MLIPs and EAM accurately capture the critical strengthening mechanisms in MoNbTaW MPEAs. It is revealed that MPEAs generally exhibit superior mechanical strength compared to their constituent metals in macroscale specimens, primarily due to solid solution strengthening during dislocation motion. However, at the nanoscale—where plasticity is predominantly governed by dislocation nucleation and grain boundary deformation—the constituent metals may outperform MPEAs. A critical length scale is identified above which MPEAs demonstrate enhanced mechanical strength relative to their constituent elements; below this scale, the advantage diminishes, underscoring a significant size-dependent effect that is crucial for optimizing MPEA applications, particularly at the nanoscale.

Understanding atomic hydrogen (H) diffusion in multi-principal element alloys (MPEAs) is crucial for enhancing hydrogen transport and storage technologies. However, the vast compositional space and complex chemical environments of MPEAs pose significant challenges. We develop highly accurate machine learning force field and neural network-driven kinetic Monte Carlo simulations to investigate H diffusion in body-centered cubic (BCC) MoNbTaW MPEAs. H diffusion exhibits super-Arrhenius behavior in MPEAs, dominated by the low percentile of the H solution energy spectrum. Robust analytical models are derived via machine learning symbolic regression to predict H diffusivity across general BCC MPEAs. Additionally, it is revealed that chemical short-range order (SRO) generally does not impact H diffusion in MoNbTaW MPEAs, except it enhances diffusion when H-favoring elements are present in low concentrations. These insights not only deepen our understanding of H diffusion dynamics in MPEAs but also guide the strategic development of advanced MPEAs for hydrogen-related applications by manipulating element type, composition, and SRO.

Titanium dioxide (TiO₂) has been widely used as a photocatalyst in CO₂ reduction reaction (CO₂RR) due to its low cost, high stability, and strong absorption in the close-to-visible ultra-violet (UV) range. However, TiO₂ films suffer from poor selectivity in CO₂ reduction due to their unfavorable electronic properties. In this work, we address this challenge by fabricating ultra-thin (14 nm) defective TiO₂ films (TiO₂-DTF) to enhance the selectivity of CO₂RR towards formate. TiO₂ sol was prepared using a facile and reproducible sol-gel method and directly deposited onto the surface of the electrode, forming a uniform, ultra-thin TiO₂ layers with a high number of defects. The activity of the TiO₂-DTF catalyst was studied in both photochemical and photoelectrochemical CO₂RR, indicating that the applied potential increases both the yield and selectivity of CO₂RR to formate. The TiO₂-DTF photocathode exhibited remarkable formate production during CO₂ reduction, achieving exceptional Faradaic efficiencies of up to 45 %. To elucidate the mechanism of photoelectrochemical CO₂RR on TiO₂-DTF, an in-situ attenuated total reflection Fourier-transform infrared spectroscopy (in-situ ATR-FTIR) was used and experimental results were supported by density functional theory (DFT) calculations. This study demonstrates that ultra-thin highly defective TiO₂ film, prepared using the cost-effective and environmentally friendly sol-gel method, can be used as photoelectrocatalyst for CO₂ reduction.

Accurate prediction of thermodynamic properties of hydrocarbons is essential for chemical process modelling. Conventional group contribution methods often are used to predict these properties. However, these methods often require extensive parameter sets to handle structural complexities. A refined group contribution method for predicting thermodynamic properties of hydrocarbon isomers with reduced complexity and improved accuracy is presented and discussed. By combining the structural framework of Constantinou and Gani (CG94) with a sensitivity-based selection of second-order groups, a reduced yet highly effective set of twelve second-order groups is identified. This reduced set retains the predictive power comparable to more complex models while significantly reducing the number of parameters. Linear regression is applied to model enthalpies and Gibbs free energies of formation for a wide temperature range. To test broader applicability, the model is further extended to properties that require nonlinear regression, including critical temperatures, critical pressures, acentric factors, and liquid densities. For all cases, the proposed model achieves high predictive accuracy, demonstrating its robustness and generalizability. This methodology balances interpretability, efficiency, and performance, making it suitable for both research and industrial thermodynamic modelling.

Extended defects such as dislocation networks and general grain boundaries are ubiquitous in metals, and accurate modeling these extensive defects is crucial to elucidate their deformation mechanisms. However, existing machine learning interatomic potentials (MLIPs) often fall short in adequately describing these defects, as their large characteristic scales exceed the computational limits of first-principles calculations. To address this challenge, we present a computational framework combining a defect genome constructed via empirical interatomic potential-guided sampling, with an automated reconstruction technique that enables accurate first-principles modeling of general defects by converting atomic clusters into periodic configurations. The effectiveness of this approach was validated through simulations of nanoindentation, tensile deformation, and fracture in BCC tungsten. This framework enhances the modeling accuracy of extended defects in crystalline materials and provides a robust foundation for advancing MLIP development by leveraging defect genomes strategically.

The high thermal stability of a thermoelectric material, which maintains a stable conversion efficiency under prolonged heat exposure, is essential for sustainable thermoelectric applications. Despite the well-known relationship between thermal degradation and microstructural evolution, their underlying interplay remains unclear, with contradictory findings reported in the literature owing to the complex dependence of microstructural changes on the material composition. Herein, the effect of Sb doping on the thermal stability of NbCoSn half-Heusler compounds is investigated in detail by comprehensively analyzing their microstructural evolution. The results reveal that introducing 3.3 at.% Sb into NbCoSn markedly enhances the thermal stability, by preserving the lattice thermal conductivity after heat exposure. Advanced techniques, including atom probe tomography, scanning transmission electron microscopy, and neutron diffraction, show that this improvement is driven by the evolution of Sb-induced complementary point defects. Although heat exposure significantly reduces lattice disorder in intrinsic NbCoSn, NbCoSn_0.9Sb_0.1 retains its lattice disorder by forming alternative point defects, thereby maintaining its lattice thermal conductivity. This detailed experimental work, corroborated by ab initio calculations, highlights the pivotal role of the point defect dynamics in achieving robust thermoelectric performances in half-Heusler compounds for high-temperature applications.

Group contribution methods (GCMs) provide a practical and computationally efficient approach for predicting thermodynamic properties of hydrocarbons, especially when experimental data are scarce. This review evaluates the evolution of GCMs from classical first-order schemes (e.g. Lydersen method, Joback method) to more advanced second-order frameworks (e.g. CG94, Sharma method), hybrid extensions, and emerging machine learning integrations. While first-order models are simple and widely used, these models struggle with branched and long-chain molecules. Second-order approaches significantly improve structural sensitivity and predictive accuracy, achieving deviations below 2–3% for critical properties and within 1 kcal/mol for formation enthalpies of branched alkanes. Nevertheless, challenges remain in extrapolating to highly complex molecules, underrepresented functional groups, and extreme conditions. Promising directions include reinforcement of second-order GCMs with molecular theory, systematic expansion of experimental and quantum-based datasets, and hybrid GCM–machine learning models that retain interpretability while improving generalisability. We recommend prioritising models that balance accuracy, robustness, simplicity, and transferability to accelerate sustainable process and product designs, particularly in applications such as fuel upgrading including hydroisomerisation, separation processes, and green chemical development.

One of the most promising energy carriers for transport applications are hydrogen-based energy carriers. NaBH₄ is a hydrogen energy carrier and produces hydrogen bubbles when it is dissolved in water. The formation of hydrogen bubbles hinders experimental measurements of the thermodynamic and transport properties of aqueous NaBH₄ solutions at elevated temperatures. Accurate knowledge of these properties is essential for the NaBH₄ hydrolysis reactor modeling and design. Molecular dynamics (MD) simulations provide the option to study the thermodynamic and transport properties of NaBH₄ aqueous solutions without hindering hydrogen bubble formation. In this work, a new force field is developed for BH₄^-, namely, the Delft force field of BH₄^- (DFF/BH₄^-), which, combined with additional force fields, can accurately describe experimental densities and viscosities of 0 to 5 m (mol salt/kg water) NaBH₄, 0 to 3 m NaB(OH)₄, and 1 m NaOH aqueous solutions at 295 K within 1.8% and 10.8% maximum deviation, respectively. Empirical fitting correlations are created for densities, viscosities, and self-diffusivities obtained from the MD simulations of 0 to 5 m NaBH₄, 0 to 5 m NaB(OH)₄, and 0 to 1 m NaOH aqueous solutions at 295-363 K for NaBH₄ hydrolysis reactor modeling and design purposes.