Hydrogen generation and related energy applications heavily rely on the hydrogen evolution reaction (HER), which faces challenges of slow kinetics and high overpotential. Efficient electrocatalysts, particularly single-atom catalysts (SACs) on two-dimensional (2D) materials, are essential. This study presents a few-shot machine learning (ML) assisted high-throughput screening of 2D septuple-atomic-layer Ga2CoS4−x supported SACs to predict HER catalytic activity. Initially, density functional theory (DFT) calculations showed that 2D Ga2CoS4 is inactive for HER. However, defective Ga2CoS4−x (x = 0–0.25) monolayers exhibit excellent HER activity due to surface sulfur vacancies (SVs), with predicted overpotentials (0–60 mV) comparable to or lower than commercial Pt/C, which typically exhibits an overpotential of around 50 mV in the acidic electrolyte, when the concentration of surface SV is lower than 8.3%. SVs generate spin-polarized states near the Fermi level, making them effective HER sites. We demonstrate ML-accelerated HER overpotential predictions for all transition metal SACs on 2D Ga2CoS4−x. Using DFT data from 18 SACs, an ML model with high prediction accuracy and reduced computation time was developed. An intrinsic descriptor linking SAC atomic properties to HER overpotential was identified. This study thus provides a framework for screening SACs on 2D materials, enhancing catalyst design.

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 development of advanced catalysts with innovative nanoarchitectures is critical for addressing energy and environmental challenges such as the electrochemical CO₂ reduction reaction (CO₂ RR). Herein, the synthesis of an innovative copper–sulfur planar structure, Cu–S–BDC, within a metal–organic framework (MOF) catalyst is presented, which demonstrates 100% selectivity toward formate as the sole carbon product. Structural analysis and surface characterizations reveal that Cu–S–BDC exhibits quasi-2D inorganic building units, with Cu bonded to two S-CH (Formula presented.) groups and one BDC linker, while carboxylate groups adopt a bridging coordination mode. This unique arrangement not only imparts remarkable structural stability but also enhances the electronic properties of the MOF, as evidenced by a narrow bandgap of 1.203 eV that facilitates efficient charge transfer and increased electrochemical current density in CO (Formula presented.) RR. Notably, it offers a Faradaic efficiency of 92% for formate at an overpotential as low as −0.4 V versus the reversible hydrogen electrode (RHE) in an aqueous electrolyte of 1 m KOH, as well as a current density of −25.8 mA cm² at −0.9 V versus RHE, averaged over 24 h of electrolysis. This study highlights a fresh perspective in the field of MOF electrocatalysts by demonstrating that engineering the metal coordination environment can significantly enhance the electronic properties and consequently improve the electrocatalytic performance of these materials.

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 conductivity predictions of KOH(aq) are crucial for electrolysis applications. OH– is transferred in water by the Grotthuss transfer mechanism, thereby increasing its mobility compared to that of other ions. Classical and ab initio molecular dynamics struggle to capture this enhanced mobility due to limitations in computational costs or in capturing chemical reactions. Most studies to date have provided only qualitative descriptions of the structure during Grotthuss transfer, without quantitative results for the transfer rate and the resulting transport properties. Here, machine learning molecular dynamics is used to investigate 50,000 transfer events. Analysis confirmed earlier works that Grotthuss transfer requires a reduction in accepted and a slight increase in donated hydrogen bonds to the hydroxide, indicating that hydrogen-bond rearrangements are rate-limiting. The computed self-diffusion coefficients and electrical conductivities are consistent with experiments for a wide temperature range, outperforming classical interatomic force fields and earlier AIMD simulations.

Recent advances in machine learning, combined with the generation of extensive density functional theory (DFT) datasets, have enabled the development of universal machine learning interatomic potentials (uMLIPs). These models offer broad applicability across the periodic table, achieving first-principles accuracy at a fraction of the computational cost of traditional DFT calculations. In this study, we demonstrate that state-of-the-art pretrained uMLIPs can effectively replace DFT for accurately modeling complex defects in a wide range of metals and alloys. Our investigation spans diverse scenarios, including grain boundaries and general defects in pure metals, defects in high-entropy alloys, hydrogen-alloy interactions, and solute-defect interactions. Remarkably, the latest EquiformerV2 models achieve DFT-level accuracy on comprehensive defect datasets, with root mean square errors below 5 meV atom−1 for energies and 100 meV Å−1 for forces, outperforming specialized machine learning potentials such as moment tensor potential and atomic cluster expansion. We also present a systematic analysis of accuracy versus computational cost and explore uncertainty quantification for uMLIPs. A detailed case study of tungsten (W) demonstrates that data on pure W alone is insufficient for modeling complex defects in uMLIPs, underscoring the critical importance of advanced machine learning architectures and diverse datasets, which include over 100 million structures spanning all elements. These findings establish uMLIPs as a robust alternative to DFT and a transformative tool for accelerating the discovery and design of high-performance materials.

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.

Rechargeable lithium–sulfur batteries (LiSBs) assembled with earth-abundant and safe Li anodes are less prone to form dendrites on the surface, and sulfur-containing cathodes offer considerable potential for achieving high energy densities. Nevertheless, suitable sulfur host materials and their interaction with electrolytes are at present key factors that retard the commercial introduction of these batteries. Here we propose a two-dimensional metallic carbon phosphorus framework, namely, 2D CP3, as a promising sulfur host material for inhibiting the shuttle effect and improving electronic conductivity in high-performance Li–S batteries. The good electrical conductivity of CP3 eliminates the insulating nature of most sulfur-based electrodes. The dissolution of lithium polysulfides (LiPSs) into the electrolyte is largely prevented by the strong interaction between CP3 and LiPSs. In addition, the deposition of Li2S on CP3 facilitates the kinetics of the LiPS redox reaction. Therefore, the use of CP3 for Li–S battery cathodes is expected to suppress the LiPS shuttle effect and to improve the overall performance, which is ideal for the practical application of Li–S batteries.

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.

Photocatalytic water splitting represents a promising approach for sustainable hydrogen production, with two-dimensional Janus materials offering unique advantages through intrinsic electric fields that enhance charge separation. We present a comprehensive first-principles investigation of Janus AlXY₂ (X = Ga, In; Y = S, Se, Te) monolayers using density functional theory and ab initio molecular dynamics simulations. All six systems exhibit excellent structural, thermal, and mechanical stability with HSE06 bandgaps of 2.029–2.969 eV suitable for UV-light absorption. The asymmetric structure generates strong intrinsic electric fields of 5.391–6.437 V perpendicular to the monolayer plane, significantly enhancing photogenerated charge carrier separation. While pristine monolayers show poor hydrogen evolution reaction (HER) activity with Gibbs free energies of 1.937–2.371 eV, strategic introduction of metal vacancies dramatically improves performance, reducing ΔG_H values to −0.371 to ＋0.607 eV and approaching optimal catalytic conditions. These findings demonstrate the potential of defect-engineered 2D Janus AlXY₂ materials for efficient photocatalytic hydrogen production.

In this study, we explore the mechanisms underlying the exceptional intrinsic strength of face-centered cubic (FCC) Multi-Principal Element Alloys (MPEAs) using a multifaceted approach. Our methods integrate atomistic simulations, informed by both embedded-atom model and neural network potentials, with first-principles calculations, stochastic Peierls-Nabarro (PN) modeling, and symbolic machine learning. We identify a consistent, robust linear correlation between the strength of MPEAs and the standard deviation of the maximum stacking-fault restoring force (τ_max,sd) across various potentials. This finding is substantiated by comparing the experimental strengths of Cantor alloys’ subsystems and Ni_62.5V_37.5 against τ_max,sd values from high-throughput first-principle calculations. Our theoretical insights are derived from integrating the stochastic Peierls-Nabarro model with a shearable precipitation hardening framework, demonstrating that lattice distortion alone does not directly enhance intrinsic strength. Instead, τ_max,sd emerges as a critical determinant, capable of boosting the strength of MPEAs by up to tenfold. Our analysis reveals the critical role of the exponential form of the PN model in achieving substantial strength improvement by transforming the Gaussian-like distribution of τ_max into an exponential-like distribution of local Peierls stress. Additionally, using an advanced symbolic machine learning technique, the sure independence screening and sparsifying operator (SISSO) method, we derive interpretable relationships between MPEA strength, elastic properties, and τ_max statistics, offering new insights into the design and optimization of advanced MPEAs. These findings highlight that the nonlinear physics and atomic fluctuations characterizing MPEAs not only underpin their unconventional intrinsic strength but also contribute to other complex properties such as sluggish diffusion and cocktail effect.

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.

Purpose: Prednisolone, a synthetic corticosteroid drug, is extensively utilized to treat inflammatory diseases and regulates metabolism and the immune system in cancer treatment. However, these drugs are toxic and cause severe side effects if administrated for long durations and in large doses. This work intends to study the atomistic interactions of popular polymeric carrier like PLGA with the drug and thereby provide insights into achieving better loading and a sustained release. Methods: Molecular dynamics (MD) simulations of prednisolone (drug) encapsulated in Poly Lactic-co-Glycolic acid (PLGA) are performed in this study. Grand Canonical Monte Carlo (GCMC) simulations with MD simulations are conducted to determine the water penetration in PLGA polymer and polymer stability in water. The investigations from this study encompasses structural and dynamical parameters, including the end-to-end distance, radius of gyration of polymer chains, interaction energy, and diffusion coefficient of the drug. Results: The polymer-drug interactions are studied and identified from the simulation data of PLGA(75:25) and PLGA(50:50) polymers with prednisolone in an aqueous medium for optimal drug carrying capacity and effective drug release. Also, the polymeric systems of PLGA(75:25) and PLGA(50:50) are analyzed with the water penetrant loading using the Grand Canonical Monte Carlo (GCMC) and MD simulations. Water loading analysis revealed that PLGA(75:25) has the highest swelling compared to PLGA(50:50). Conclusion: This study highlights the characteristics and critical parameters for developing an optimal drug delivery system by investigating polymer-drug interactions, drug encapsulation, and water uptake in polymers using MD and GCMC simulations.

Efficiently designing lightweight alloys with combined high corrosion resistance and mechanical properties remains an enduring topic in materials engineering. Due to the inadequate accuracy of conventional stress-strain machine learning (ML) models caused by corrosion factors, a novel reinforcement self-learning ML algorithm combined with calculated features (accuracy R² >0.92) is developed. Based on the ML models, calculated work functions and mechanical moduli, a Computation Designed Corrosion-Resistant Al alloy is fabricated and verified. The performance (elongation reaches ∼30 %) is attributed to the H trapping Al-Sc-Cu phases (-1.44 eV H⁻¹) and Cu-modified η/η' precipitates inside the grain boundaries (GBs).

Vapor-Liquid Equilibria (VLE) of hydrogen (H₂) and aqueous electrolyte (KOH and NaCl) solutions are central to numerous industrial applications such as alkaline electrolysis and underground hydrogen storage. Continuous fractional component Monte Carlo simulations are performed to compute the VLE of H₂ and aqueous electrolyte solutions at 298-423 K, 10-400 bar, 0-8 mol KOH/kg water, and 0-6 mol NaCl/kg water. The densities and activities of water in aqueous KOH and NaCl solutions are accurately modeled (within 2% deviation from experiments) using the non-polarizable Madrid-2019 Na⁺/Cl⁻ ion force fields for NaCl and the Madrid-Transport K⁺ and Delft Force Field of OH⁻ for KOH, combined with the TIP4P/2005 water force field. A free energy correction (independent of pressure, salt type, and salt molality) is applied to the computed infinite dilution excess chemical potentials of H₂ and water, resulting in accurate predictions (within 5% of experiments) for the solubilities of H₂ in water and the saturated vapor pressures of water for a temperature range of 298-363 K. The compositions of water and H₂ are computed using an iterative scheme from the liquid phase excess chemical potentials and densities, in which the gas phase fugacities are computed using the GERG-2008 equation of state. For the first time, the VLE of H₂ and aqueous KOH/NaCl systems are accurately captured with respect to experiments (i.e., for both the liquid and gas phase compositions) without compromising the liquid phase properties or performing any refitting of force fields.

The technological landscape for industrial processes handling asphaltene is evolving at a rapid pace due to the increase in the extraction of heavy crude oil. The main underlying challenges in this regard are the flow assurance, the recovery of the spent solvent, and the sophisticated extractor setup required to develop the process to an industrial scale. The number of studies focused on the handling of the asphaltene at the atomic and molecular scales is growing enormously in order to identify new sustainable solvents for the effective extraction of asphaltene from heavy crude oil or oil-bearing sands. This Perspective focuses on the importance of density functional theory and molecular dynamics simulations to explore the broader range of asphaltene inhibitors, e.g., nanoparticles, ionic liquids, and deep eutectic solvents, to prevent asphaltene precipitation. We provide a concise overview of the major accomplishments, analyze the aspects that require attention, and highlight the path-breaking studies having a significant impact on the process of chemical enhanced oil recovery from heavy crude oil reservoirs primarily based on atomistic and molecular simulations.

This study aims to identify photo-/electrocatalysts that can enhance the oxygen evolution reaction (OER), hydrogen evolution reaction (HER), and oxygen reduction reaction (ORR), which are of utmost importance in electro-/photochemical energy systems, such as solar energy, fuel cells, water electrolyzers, or metal-air batteries. Our study focused on investigating the 2D Ge₂Se₂P₄ monolayer and found that it exhibits a bifunctional photocatalyst with a very high solar-to-hydrogen efficiency. The two-dimensional (2D) Ge₂Se₂P₄ monolayer has superior HER activity compared to that of most 2D materials, and it also outperforms the reference catalysts IrO₂(110) and Pt(111) in terms of low overpotential values for ORR and OER mechanisms. Such superior catalytic performance in the 2D Ge₂Se₂P₄ monolayer can be attributed to its electron states, charge transfer process, and suitable band alignments referring to normal hydrogen electrodes. Overall, the study suggests that the Ge₂Se₂P₄ monolayer could be an excellent bifunctional catalyst for advancing photo-/electrochemical energy systems.

Non-polarizable force fields fail to accurately predict free energies of aqueous electrolytes without compromising the predictive ability for densities and transport properties. A new approach is presented in which (1) TIP4P/2005 water and scaled charge force fields are used to describe the interactions in the liquid phase and (2) an additional Effective Charge Surface (ECS) is used to compute free energies at zero additional computational expense. The ECS is obtained using a single temperature-independent charge scaling parameter per species. Thereby, the chemical potential of water and the free energies of hydration of various aqueous salts (e.g., NaCl and LiCl) are accurately described (deviations less than 5% from experiments), in sharp contrast to calculations where the ECS is omitted (deviations larger than 20%). This approach enables accurate predictions of free energies of aqueous electrolyte solutions using non-polarizable force fields, without compromising liquid-phase properties.

The introduction of pharmaceuticals into aquatic ecosystems can lead to the generation of antibiotic-resistant bacteria. This paper employed molecular dynamics simulations to examine the interactions between cationic/anionic surfactants and two antibiotics or drugs, namely, ciprofloxacin and azithromycin. The analysis focused on many factors to elucidate the mechanism by which the surfactant bilayer molecular structure affects the selected antibiotics. These factors include the tilt angle, rotational angle of the surfactants, electrostatic potential, and charge density along the bilayers. Our molecular-level investigation of the adsorption mechanisms of hydrophobic (azithromycin) and hydrophilic (ciprofloxacin) drugs on the cationic/anionic surfactant bilayer offers a crucial understanding for comprehending the optimal selection of surfactants for effectively separating antibiotics.