Critical role of slags in pitting corrosion of additively manufactured stainless steel in simulated seawater.

Pitting corrosion in seawater is one of the most difficult forms of corrosion to identify and control. A workhorse material for marine applications, 316L stainless steel (316L SS) is known to balance resistance to pitting with good mechanical properties. The advent of additive manufacturing (AM), particularly laser powder bed fusion (LPBF), has prompted numerous microstructural and mechanical investigations of LPBF 316L SS; however, the origins of pitting corrosion on as-built surfaces is unknown, despite their utmost importance for certification of LPBF 316L SS prior to fielding. Here, we show that Mn-rich silicate slags are responsible for pitting of the as-built LPBF material in sodium chloride due to their introduction of deleterious defects such as cracks or surface oxide heterogeneities. In addition, we explain how slags are formed in the liquid metal and deposited at the as-built surfaces using high-fidelity melt pool simulations. Our work uncovers how LPBF changes surface oxides due to rapid solidification and high-temperature oxidation, leading to fundamentally different pitting corrosion mechanisms.

Quantifying disorder one atom at a time using an interpretable graph neural network paradigm.

Quantifying the level of atomic disorder within materials is critical to understanding how evolving local structural environments dictate performance and durability. Here, we leverage graph neural networks to define a physically interpretable metric for local disorder, called SODAS. This metric encodes the diversity of the local atomic configurations as a continuous spectrum between the solid and liquid phases, quantified against a distribution of thermal perturbations. We apply this methodology to four prototypical examples with varying levels of disorder: (1) grain boundaries, (2) solid-liquid interfaces, (3) polycrystalline microstructures, and (4) tensile failure/fracture. We also compare SODAS to several commonly used methods. Using elemental aluminum as a case study, we show how our paradigm can track the spatio-temporal evolution of interfaces, incorporating a mathematically defined description of the spatial boundary between order and disorder. We further show how to extract physics-preserved gradients from our continuous disorder fields, which may be used to understand and predict materials performance and failure. Overall, our framework provides a simple and generalizable pathway to quantify the relationship between complex local atomic structure and coarse-grained materials phenomena.

Hydrogen Storage in Partially Exfoliated Magnesium Diboride Multilayers.

Metal boride nanostructures have shown significant promise for hydrogen storage applications. However, the synthesis of nanoscale metal boride particles is challenging because of their high surface energy, strong inter- and intraplanar bonding, and difficult-to-control surface termination. Here, it is demonstrated that mechanochemical exfoliation of magnesium diboride in zirconia produces 3-4 nm ultrathin MgB2 nanosheets (multilayers) in high yield. High-pressure hydrogenation of these multilayers at 70 MPa and 330 °C followed by dehydrogenation at 390 °C reveals a hydrogen capacity of 5.1 wt%, which is ≈50 times larger than the capacity of bulk MgB2 under the same conditions. This enhancement is attributed to the creation of defective sites by ball-milling and incomplete Mg surface coverage in MgB2 multilayers, which disrupts the stable boron-boron ring structure. The density functional theory calculations indicate that the balance of Mg on the MgB2 nanosheet surface changes as the material hydrogenates, as it is energetically favorable to trade a small number of Mg vacancies in Mg(BH4 )2 for greater Mg coverage on the MgB2 surface. The exfoliation and creation of ultrathin layers is a promising new direction for 2D metal boride/borohydride research with the potential to achieve high-capacity reversible hydrogen storage at more moderate pressures and temperatures.

Flexible machine-learning interatomic potential for simulating structural disordering behavior of Li7La3Zr2O12 solid electrolytes.

Batteries based on solid-state electrolytes, including Li7La3Zr2O12 (LLZO), promise improved safety and increased energy density; however, atomic disorder at grain boundaries and phase boundaries can severely deteriorate their performance. Machine-learning (ML) interatomic potentials offer a uniquely compelling solution for simulating chemical processes, rare events, and phase transitions associated with these complex interfaces by mixing high scalability with quantum-level accuracy, provided that they can be trained to properly address atomic disorder. To this end, we report the construction and validation of an ML potential that is specifically designed to simulate crystalline, disordered, and amorphous LLZO systems across a wide range of conditions. The ML model is based on a neural network algorithm and is trained using ab initio data. Performance tests prove that the developed ML potential can predict accurate structural and vibrational characteristics, elastic properties, and Li diffusivity of LLZO comparable to ab initio simulations. As a demonstration of its applicability to larger systems, we show that the potential can correctly capture grain boundary effects on diffusivity, as well as the thermal transition behavior of LLZO. These examples show that the ML potential enables simulations of transitions between well-defined and disordered structures with quantum-level accuracy at speeds thousands of times faster than ab initio methods.

Heteroatom-Doped Graphenes as Actively Interacting 2D Encapsulation Media for Mg-Based Hydrogen Storage.

Nanoencapsulation using graphene derivatives enables the facile fabrication of two-dimensional (2D) nanocomposites with unique microstructures and has been generally applied to many fields of energy materials. Particularly, metal hydrides such as MgH2 encapsulated by graphene derivatives have emerged as a promising hybrid material for overcoming the disadvantageous properties of Mg-based hydrogen storage. Although the behavior of the graphene-Mg nanoencapsulation interface has been studied for many composite materials, the direct modification of graphene with nonmetal foreign elements for changing the interfacial behavior has been limitedly reported. In this regard, using B-doped graphene and N-doped graphene as nanoencapsulation media for tuning the interfacial behavior of graphene derivative-Mg nanoparticles, we present altered hydrogen storage kinetics of heteroatom-doped (B and N) graphene-Mg composites. The effect of heteroatom doping is studied in terms of bonding configurations and heteroatom doping concentrations. The enhancement in hydrogen uptake was observed for all of the heteroatom-doped graphene-Mg nanocomposites. On the other hand, a few samples exhibit significantly low activation energy at the early stage of desorption, which can be related to the facilitated nucleus formation. Density functional theory calculation indicates that B-doping and N-doping accelerate hydrogen absorption kinetics in different ways, aiding charge transfer and inducing surface deformation of Mg nanoparticles, respectively. Their effects can be augmented in the presence of structural defects on graphene, such as vacancies, pores, or graphene edges. These results demonstrate that hydrogen storage kinetics of Mg-based systems can be altered by utilizing heteroatom-doped graphene oxide derivatives as 2D nanoencapsulation media, suggesting that the addition of a nonmetal doping element can also be applied to Mg-based hydrogen storage by modifying the nanoencapsulation interface without forming Mg alloy phases.

Understanding Hydrogenation Chemistry at MgB2 Reactive Edges from Ab Initio Molecular Dynamics.

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB2 without a priori assumption of reaction pathways. Focusing on highly reactive (101̅0) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH4- units and key chemical intermediates, involving H2 dissociation, generation of functionalities and molecular complexes containing BH2 and BH3 motifs, and B-B bond breaking. The genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B-H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB2-based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

First-Principles Elucidation of Initial Dehydrogenation Pathways in Mg(BH4)2.

Complex borohydrides such as Mg(BH4)2 offer one of highest capacities to chemically store hydrogen for onboard applications; however, it suffers greatly from kinetic constraints that prevent realization of full capacity and reversibility. Understanding these kinetic limitations solely from experiments is extremely challenging due to the unusual complexity of various competing elemental reaction steps involved during the de/rehydrogenation reaction. This work aims to map out the energetics associated with initial dehydrogenation of Mg(BH4)2 from first-principles simulations and to identify the preferred reaction pathways. Our calculations suggest the rate-limiting step during BH4--B3H8- conversion is the formation of the B2H7- intermediate. We further emphasize and clarify that the B3H8- and H- intermediates, formed during initial Mg(BH4)2 decomposition, appear as molecular species that are embedded in the Mg-BH4-Mg matrix as evidenced in the nuclear magnetic resonance measurements and not as bulk MgH2 and Mg(B3H8)2 as previously assumed in theoretical predictions of the thermodynamics.

Beyond Thermodynamics: Assessing the Dynamical Softness of Hydrated Ions from First Principles.

Ion (de)hydration is a key rate-determining step in interfacial processes from corrosion to electrochemical energy storage. However, predicting the kinetics of ion (de)hydration remains challenging, prompting the use of static proxies such as hydration energy and valence. While useful for assessing thermodynamic preferences, such descriptors cannot fully capture the dynamical softness of the hydration shell that dictates kinetics. Accordingly, we use first-principles molecular dynamics to analyze hydration shell softness for a diverse set of metal cations. Three dynamic metrics are introduced to intuitively describe the bond rigidity, shape deformability, and exchange fluidity of the solvation shell. Together, these metrics capture the relevant physics in the static descriptors, while offering a far more complete and efficient representation for the overall propensity for (de)hydration. Application to the hydrated ion set demonstrates a weak connection between dynamical softness and hydration energy, confirming that dynamical descriptors of hydration are key for correctly describing ion transfer processes.

Spontaneous dynamical disordering of borophenes in MgB2 and related metal borides.

Layered boron compounds have attracted significant interest in applications from energy storage to electronic materials to device applications, owing in part to a diversity of surface properties tied to specific arrangements of boron atoms. Here we report the energy landscape for surface atomic configurations of MgB2 by combining first-principles calculations, global optimization, material synthesis and characterization. We demonstrate that contrary to previous assumptions, multiple disordered reconstructions are thermodynamically preferred and kinetically accessible within exposed B surfaces in MgB2 and other layered metal diborides at low boron chemical potentials. Such a dynamic environment and intrinsic disordering of the B surface atoms present new opportunities to realize a diverse set of 2D boron structures. We validated the predicted surface disorder by characterizing exfoliated boron-terminated MgB2 nanosheets. We further discuss application-relevant implications, with a particular view towards understanding the impact of boron surface heterogeneity on hydrogen storage performance.

Paradigms of frustration in superionic solid electrolytes.

Superionic solid electrolytes have widespread use in energy devices, but the fundamental motivations for fast ion conduction are often elusive. In this Perspective, we draw upon atomistic simulations of a wide range of superionic conductors to illustrate some ways frustration can lower diffusion cation barriers in solids. Based on our studies of halides, oxides, sulfides and hydroborates and a survey of published reports, we classify three types of frustration that create competition between different local atomic preferences, thereby flattening the diffusive energy landscape. These include chemical frustration, which derives from competing factors in the anion-cation interaction; structural frustration, which arises from lattice arrangements that induce site distortion or prevent cation ordering; and dynamical frustration, which is associated with temporary fluctuations in the energy landscape due to anion reorientation or cation reconfiguration. For each class of frustration, we provide detailed simulation analyses of various materials to show how ion mobility is facilitated, resulting in stabilizing factors that are both entropic and enthalpic in origin. We propose the use of these categories as a general construct for classifying frustration in superionic conductors and discuss implications for future development of suitable descriptors and improvement strategies. This article is part of the Theo Murphy meeting issue 'Understanding fast-ion conduction in solid electrolytes'.

Impacts of vacancy-induced polarization and distortion on diffusion in solid electrolyte Li3OCl.

Lithium-rich oxychloride antiperovskites are promising solid electrolytes for enabling next-generation batteries. Here, we report a comprehensive study varying Li+ concentrations in [Formula: see text] using ab initio molecular dynamics simulations. The simulations accurately capture the complex interactions between Li+ vacancies ([Formula: see text]), the dominant mobile species in [Formula: see text]. The [Formula: see text] polarize and distort the host lattice, inducing additional non-vacancy-mediated diffusion mechanisms and correlated diffusion events that reduce the activation energy barrier at concentrations as low as 1.5% [Formula: see text]. Our analyses of discretized diffusion events in both space and time illustrate the critical interplay between correlated dynamics, polarization and local distortion in promoting ionic conductivity in [Formula: see text]. This article is part of the Theo Murphy meeting issue 'Understanding fast-ion conduction in solid electrolytes'.

Evaluating the stability and activity of dilute Cu-based alloys for electrochemical CO2 reduction.

Cu-based catalysts currently offer the most promising route to actively and selectively produce value-added chemicals via electrochemical reduction of CO2 (eCO2R); yet further improvements are required for their wide-scale deployment in carbon mitigation efforts. Here, we systematically investigate a family of dilute Cu-based alloys to explore their viability as active and selective catalysts for eCO2R through a combined theoretical-experimental approach. Using a quantum-classical modeling approach that accounts for dynamic solvation effects, we assess the stability and activity of model single-atom catalysts under eCO2R conditions. Our calculations identify that the presence of eCO2R intermediates, such as CO*, H*, and OH*, may dynamically influence the local catalyst surface composition. Additionally, we identify through binding energy descriptors of the CO*, CHO*, and OCCO* dimer intermediates that certain elements, such as group 13 elements (B, Al, and Ga), enhance the selectivity of C2+ species relative to pure Cu by facilitating CO dimerization. The theoretical work is corroborated by preliminary testing of eCO2R activity and selectivity of candidate dilute Cu-based alloy catalyst films prepared by electron beam evaporation in a zero-gap gas diffusion electrode-based reactor. Of all studied alloys, dilute CuAl was found to be the most active and selective toward C2+ products like ethylene, consistent with the theoretical predictions. We attribute the improved performance of dilute CuAl alloys to more favorable dimerization reaction energetics of bound CO species relative to that on pure Cu. In a broader context, the results presented here demonstrate the power of our simulation framework in terms of rational catalyst design.

Defying Thermodynamics: Stabilization of Alane Within Covalent Triazine Frameworks for Reversible Hydrogen Storage.

The highly unfavorable thermodynamics of direct aluminum hydrogenation can be overcome by stabilizing alane within a nanoporous bipyridine-functionalized covalent triazine framework (AlH3 @CTF-bipyridine). This material and the counterpart AlH3 @CTF-biphenyl rapidly desorb H2 between 95 and 154 °C, with desorption complete at 250 °C. Sieverts measurements, 27 Al MAS NMR and 27 Al{1 H} REDOR experiments, and computational spectroscopy reveal that AlH3 @CTF-bipyridine dehydrogenation is reversible at 60 °C under 700 bar hydrogen, >10 times lower pressure than that required to hydrogenate bulk aluminum. DFT calculations and EPR measurements support an unconventional mechanism whereby strong AlH3 binding to bipyridine results in single-electron transfer to form AlH2 (AlH3 )n clusters. The resulting size-dependent charge redistribution alters the dehydrogenation/rehydrogenation thermochemistry, suggesting a novel strategy to enable reversibility in high-capacity metal hydrides.

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields.

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. We have focused our workflow on TiH2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

Reversing the Irreversible: Thermodynamic Stabilization of LiAlH4 Nanoconfined Within a Nitrogen-Doped Carbon Host.

A general problem when designing functional nanomaterials for energy storage is the lack of control over the stability and reactivity of metastable phases. Using the high-capacity hydrogen storage candidate LiAlH4 as an exemplar, we demonstrate an alternative approach to the thermodynamic stabilization of metastable metal hydrides by coordination to nitrogen binding sites within the nanopores of N-doped CMK-3 carbon (NCMK-3). The resulting LiAlH4@NCMK-3 material releases H2 at temperatures as low as 126 °C with full decomposition below 240 °C, bypassing the usual Li3AlH6 intermediate observed in bulk. Moreover, >80% of LiAlH4 can be regenerated under 100 MPa H2, a feat previously thought to be impossible. Nitrogen sites are critical to these improvements, as no reversibility is observed with undoped CMK-3. Density functional theory predicts a drastically reduced Al-H bond dissociation energy and supports the observed change in the reaction pathway. The calculations also provide a rationale for the solid-state reversibility, which derives from the combined effects of nanoconfinement, Li adatom formation, and charge redistribution between the metal hydride and the host.

Neural network potential from bispectrum components: A case study on crystalline silicon.

In this article, we present a systematic study on developing machine learning force fields (MLFFs) for crystalline silicon. While the main-stream approach of fitting a MLFF is to use a small and localized training set from molecular dynamics simulations, it is unlikely to cover the global features of the potential energy surface. To remedy this issue, we used randomly generated symmetrical crystal structures to train a more general Si-MLFF. Furthermore, we performed substantial benchmarks among different choices of material descriptors and regression techniques on two different sets of silicon data. Our results show that neural network potential fitting with bispectrum coefficients as descriptors is a feasible method for obtaining accurate and transferable MLFFs.

Structural Anomalies and Electronic Properties of an Ionic Liquid under Nanoscale Confinement.

Ionic liquids (ILs) promise far greater electrochemical performance compared to aqueous systems, yet key physicochemical properties governing their assembly at interfaces within commonly used graphitic nanopores remain poorly understood. In this work, we combine synchrotron X-ray scattering with first-principles molecular dynamics simulations to unravel key structural characteristics of 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([TFSI]-) ionic liquids confined in carbon slit pores. X-ray scattering reveals selective pore filling due to size exclusion, while filled pores exhibit disruption in the IL intermolecular structure, the extent of which increases for narrower slit pores. First-principles simulations corroborate this finding and quantitatively describe how perturbations in the local IL structure, particularly the hydrogen-bond network, depend strongly on the degree of confinement. Despite significant deviations in structure under confinement, electrochemical stability remains intact, which is important for energy storage based on nanoporous carbon electrodes (e.g., supercapacitors).

Nanoconfinement of Molecular Magnesium Borohydride Captured in a Bipyridine-Functionalized Metal-Organic Framework.

The lower limit of metal hydride nanoconfinement is demonstrated through the coordination of a molecular hydride species to binding sites inside the pores of a metal-organic framework (MOF). Magnesium borohydride, which has a high hydrogen capacity, is incorporated into the pores of UiO-67bpy (Zr6O4(OH)4(bpydc)6 with bpydc2- = 2,2'-bipyridine-5,5'-dicarboxylate) by solvent impregnation. The MOF retained its long-range order, and transmission electron microscopy and elemental mapping confirmed the retention of the crystal morphology and revealed a homogeneous distribution of the hydride within the MOF host. Notably, the B-, N-, and Mg-edge XAS data confirm the coordination of Mg(II) to the N atoms of the chelating bipyridine groups. In situ 11B MAS NMR studies helped elucidate the reaction mechanism and revealed that complete hydrogen release from Mg(BH4)2 occurs as low as 200 °C. Sieverts and thermogravimetric measurements indicate an increase in the rate of hydrogen release, with the onset of hydrogen desorption as low as 120 °C, which is approximately 150 °C lower than that of the bulk material. Furthermore, density functional theory calculations support the improved dehydrogenation properties and confirm the drastically lower activation energy for B-H bond dissociation.

Toward Engineering of Solution Microenvironments for the CO2 Reduction Reaction: Unraveling pH and Voltage Effects from a Combined Density-Functional-Continuum Theory.

Engineering the electrolyte microenvironment represents an attractive route to tuning the selectivity of electrocatalytic reactions beyond catalyst composition and morphology. However, harnessing the full potential of this approach requires understanding the interplay between voltage, electrolyte composition, and adsorbate binding within the electrical double layer, which is absent from the usual theoretical approaches. In this work, we apply a recently developed density functional theory (DFT)-continuum approach based on the effective screening medium method and reference interaction site model (ESM-RISM) to explore electrolyte effects with an enhanced description of the electrochemical interface. Applying this method to the binding of CO adsorbates in potassium-containing electrolytes on copper, a problem of direct relevance to CO2 electroreduction to value-added products, we show that the interdependence of voltage and pH leads to an unexpected change in adsorption site preference on Cu(001) terraces. Our findings highlight the often-overlooked importance of the electrical double-layer structure for predicting catalyst operation.

A Mechanistic Analysis of Phase Evolution and Hydrogen Storage Behavior in Nanocrystalline Mg(BH4)2 within Reduced Graphene Oxide.

Magnesium borohydride (Mg(BH4)2, abbreviated here MBH) has received tremendous attention as a promising onboard hydrogen storage medium due to its excellent gravimetric and volumetric hydrogen storage capacities. While the polymorphs of MBH-alpha (α), beta (ß), and gamma (γ)-have distinct properties, their synthetic homogeneity can be difficult to control, mainly due to their structural complexity and similar thermodynamic properties. Here, we describe an effective approach for obtaining pure polymorphic phases of MBH nanomaterials within a reduced graphene oxide support (abbreviated MBHg) under mild conditions (60-190 °C under mild vacuum, 2 Torr), starting from two distinct samples initially dried under Ar and vacuum. Specifically, we selectively synthesize the thermodynamically stable α phase and metastable ß phase from the γ-phase within the temperature range of 150-180 °C. The relevant underlying phase evolution mechanism is elucidated by theoretical thermodynamics and kinetic nucleation modeling. The resulting MBHg composites exhibit structural stability, resistance to oxidation, and partially reversible formation of diverse [BH4]- species during de- and rehydrogenation processes, rendering them intriguing candidates for further optimization toward hydrogen storage applications.