Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic transitions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Molecular Models of Atomically Dispersed Uranium at MoS 2 Surfaces Reveal Cooperative Mechanism of Water Reduction

Single atoms of uranium supported on molybdenum sulfide surfaces (U@MoS 2 ) have been recently demonstrated to facilitate the hydrogen evolution reaction (HER) through electrocatalysis. Theoretical calculations have predicted uranium hydroxide moieties bound to edge-sulfur atoms of MoS 2 as a proposed transition state involved in the HER process. However, the isolation of relevant intermediates involved in this process remains a challenge, rendering mechanistic hypotheses unverified. The present work describes the isolation and characterization of a uranium-hydroxide intermediate on molybdenum sulfide surfaces using [(Cp* 3 Mo 3 S 4 )UCp*], a molecular model of a reduced uranium center supported at MoS 2 . Mechanistic investigations highlight the metalloligand cooperativity with uranium involved in the water activation pathway. The corresponding uranium-oxo analogue, [(Cp* 3 Mo 3 S 4 )Cp*U(═O)], was also accessed from the hydroxide cluster via hydrogen atom transfer and from [(Cp* 3 Mo 3 S 4 )UCp*] through an alternative direct oxygen atom transfer. These results provide an atomistic perspective on the reactivity of low-valent uranium at molybdenum sulfide surfaces toward water, modeling key intermediates associated with the HER of U@MoS 2 catalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Mixed excitonic nature in water-oxidized BiVO 4 surfaces with defects

BiVO 4 is a promising photocatalyst for efficient water oxidation, with surface reactivity determined by the structure of active catalytic sites. Surface oxidation in the presence of oxygen vacancies induces electron localization, suggesting an atomistic route to improve the charge transfer efficiency within the catalytic cycle. In this paper, we study the effect of oxygen vacancies on the electronic and optical properties at BiVO 4 surfaces upon water oxidation. We use density functional theory and many-body perturbation theory to explore the change in the electronic and quasiparticle energy levels and to evaluate the electron-hole coupling as a function of the underlying structure. Further we show that while the presence of defects alters the atomic structure and largely modifies the wave-function nature, leading to defect-localized states at the quasiparticle gap region, the optical excitations remain largely unchanged due to the substantial hybridization of defect and nondefect electron-hole transitions. Our findings suggest that defect-induced surface oxidation supports improved electron transport, both through bound and tunable electronic states and via a mixed nature of the optical transitions, expected to reduce electron-hole defect trapping.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Subnanosecond phase transition dynamics in laser-shocked iron

Iron is one of the most studied chemical elements due to its sociotechnological and planetary importance; hence, understanding its structural transition dynamics is of vital interest. By combining a short pulse optical laser and an ultrashort free electron laser pulse, we have observed the subnanosecond structural dynamics of iron from high-quality x-ray diffraction data measured at 50-ps intervals up to 2500 ps. We unequivocally identify a three-wave structure during the initial compression and a two-wave structure during the decaying shock, involving all of the known structural types of iron (α-, γ-, and ε-phase). In the final stage, negative lattice pressures are generated by the propagation of rarefaction waves, leading to the formation of expanded phases and the recovery of γ-phase. Our observations demonstrate the unique capability of measuring the atomistic evolution during the entire lattice compression and release processes at unprecedented time and strain rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Segregation competition and complexion coexistence within a polycrystalline grain boundary network

Interfacial segregation can stabilize grain structures and even lead to grain boundary complexion transitions. However, understanding of the complexity of such phenomena in polycrystalline materials is limited, as most studies focus on bicrystal geometries. In this work, we investigate interfacial segregation and subsequent complexion transitions in polycrystalline Cu-Zr alloys using hybrid Monte Carlo/molecular dynamics simulations. No significant change in the grain size or structure is observed upon Zr dopant addition to a pure Cu polycrystal at moderate temperature, where grain boundary segregation is the dominant behavior. Segregation within the boundary network is inhomogeneous, with some boundaries having local concentrations that are an order of magnitude larger than the global value and others having almost no segregation, and changes to physical parameters such as boundary free volume and energy are found to correlate with dopant concentration. Further, another alloy sample is investigated at a higher temperature to probe the occurrence of widespread transitions in interfacial structure, where a significant fraction of the originally ordered boundaries transition to amorphous complexions, demonstrating the coexistence of multiple complexion types, each with their own distribution of boundary chemical composition. Overall, this work highlights that interfacial segregation and complexion structure can be diverse in a polycrystalline network. Here, the findings shown here complement existing computational and experimental studies of individual interfaces and help pave the way for unraveling the complexity of interfacial structure in realistic microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In-situ HRTEM observations of intermediate phase transformation in lattice reorientation in HCP rhenium

In-situ high-resolution transmission electron microscopy (HRTEM) is performed to investigate the deformation behavior of hexagonal close-packed rhenium (Re) which is compressed along the $\langle1\overline{1}00\rangle$ direction. Atomistic simulations are also conducted to better understand the deformation mechanisms. Two types of lattice reorientation are observed during compression. The first type involves the reorientation of one lattice by ~90° around $\langle1\overline{1}00\rangle$, which is accomplished by the formation of an intermediate face-center-cubic (FCC) phase at the interface. This transformation sequence can be described as {$1\overline{1}00$} matrix → {$111$} FCC → ($0001$) twin . In the second type, a new grain is formed but does not satisfy any known twin relationship with the matrix, and an intermediate FCC phase is also formed. The transformation sequence can be described as {$1\overline{1}00$} matrix → {$111$} FCC → ($0001$) grain . Mechanisms responsible for the observed lattice reorientation and sequential phase transitions are analyzed by conducting lattice correspondence analyses on the simulation results. Strain accommodation is also analyzed to explain the mechanisms for lattice reorientation and the intermediate phase transformations. In conclusion, the results provide new insight into the deformation behavior of HCP metals

36 MATERIALS SCIENCE↗
DOE OSTI2025

Efficient machine learning interatomic potentials robust for liquid and multiple solid polymorphs of NaF and KF

Achieving atomic-level understanding of crystallization of molten salts is of importance to a wide range of technological applications. Recent work [Fan et al., Proc. Natl. Acad. Sci. USA 122, e2425702122 (2025)] revealed that crystal nucleation in molten LiF salt is a multistage process according to the molecular-dynamics (MD) simulations based on an atomic cluster expansion (ACE) machine-learning interatomic potential (MLIP). In order to understand the influence of increasing cation size on nucleation pathways and nucleation rates of molten fluoride salts, here we develop two new ACE MLIPs for NaF and KF. The two ACE MLIPs feature DFT-SCAN-level accuracy for liquid and multiple solid polymorphs over a wide temperature (0–2000 K) and pressure (0–100 GPa) range, and also reproduce well a number of experimental data for solid and liquid equilibrium properties. The efficiency of the two ACE MLIPs enable million-atom-scale or microsecond-scale MD simulations. The two general-purpose ACE MLIPs are expected to be useful for atomistic simulations for different purposes, in addition to studying crystallization of molten NaF and KF salts.

Crystal melting↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗
DOE OSTI2025

Dynamic Interfacial Architectures: Cruciferin‐Stabilized Oil/Water Interfaces for Sustainable Emulsions

Stabilizing oil-water interfaces in emulsions by plant-based proteins provides sustainable and tunable ways for designing emulsions with specific properties, for food, healthcare, and pharmaceuticals. Cruciferin, a protein from rapeseed, has great potential as green emulsifier, but details about its structure and mobility at oil-water interfaces are largely unknown. Here, these properties are studied with small angle neutron and x-ray scattering, and neutron spin echo spectroscopy, analyzed by atomistic modelling of scattering curves and coarse-grained modelling, to gain insight into interface coverage, and molecular conformation and mobility at the interface. Cruciferin assumes trimeric conformations at the interface, as in solution, but with its protrusions from the central core of the subunits (“arms”) more compressed. Interfacial mobility is only marginally lower than in solution, indicating the arms still transiently extend and preserve a network, for the first time revealing the mechanism how cruciferin forms highly elastic 2d gel-like oil-water interfaces, as observed in macroscopic rheology. The high interfacial mobility may help in self-repairing non-stabilized interfacial fractions, reducing coalescence. These findings provide a deeper molecular level understanding of proteins at oil-water interfaces, which can stimulate development of new plant-based emulsion products, and contribute to the global protein transition.

cruciferin↗
DOE OSTI2025

Machine Learning Thermodynamics And Kinetics of Defects For Accelerated Materials Discovery

Atomistic defects play a pivotal role in functional and structural materials’ performance across a myriad of technology applications. Quantitative prediction of the thermodynamics and kinetics of defect formation and migration, respectively, typically requires accurate but expensive first-principles approaches, such as density functional theory (DFT). Their computational expense limits the throughput needed to perform high-throughput materials discovery/screening exercises or to perform materials modeling tasks relying on extensive sampling techniques. Therefore, in this Sandia National Laboratories Laboratory Directed Research and Development (LDRD) project (Project #229366), we developed a variety of machine learning techniques, trained on density functional theory calculations, to accelerate the discovery and modeling of materials in which vacancy and interstitial defects primarily dictate material performance. These include applications such as metal oxides for water-splitting or mixed ionic-electronic conduction, metal hydrides for hydrogen storage, and transition metal dichalcogenides for electronics, and the approaches developed herein can further be applied to many other domains that similarly depend on materials’ thermodynamic and kinetic defect properties for their desired functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2024

Tuning Two-Dimensional Phthalocyanine Dual Site Metal–Organic Framework Catalysts for the Oxygen Reduction Reaction

Metal-organic frameworks (MOFs) offer an interesting opportunity for catalysis, particularly for metal-nitrogen-carbon (M-N-C) motifs by providing an organized porous structural pattern and well-defined active sites for the oxygen reduction reaction (ORR), a key need for hydrogen fuel cells and related sustainable energy technologies. Here, in this work, we leverage electrochemical testing with computational models to study the electronic and structural properties in these systems and their relationship to ORR activity and stability based on dual transitional metal centers. These consists of two M1 metals with amine nodes coordinated to a single M2 metal with a phthalocyanine linker, where M1/M2 = Co, Ni, or Cu. Co-based metal centers, in particular Ni-Co, demonstrate the highest overall activity of all nine tested MOFs. Computationally, we identify the dominance of Co-sites, relative higher importance of the M2 site, and the role of layer M1 interactions on the ORR activity. Selectivity measurements indicate that M1 sites of MOFs, particularly Co, exhibits lowest (< 4%), and Ni demonstrates highest (>46%) two-electron selectivity, in good agreement with computational studies. Direct in-situ stability characterization, measuring dissolved metal ions, and calculations, using an alkaline stability metric, confirm that Co is the most stable metal in the MOF, while Cu exhibits notable instability at the M1. Overall, this study reveals how atomistic coupling of electronic and structural properties affects the ORR performance of dual site MOF catalysts and opens new avenues for tunable design and future development of these systems for practical electrochemical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A Universal Model of Cation Effects in Electrocatalysis

Electrolyte cations are conventionally viewed as inert spectators in electrocatalysis. However, a wealth of observations show that catalytic rates are often highly sensitive to cation identity. Despite their prevalence, these cation effects have resisted a unified mechanistic explanation, with different physical phenomena implicated across reaction chemistries, catalyst compositions, and choice of solvent. In this perspective, we describe a general framework for understanding cation effects in electrocatalysis based on electrostatics. We argue that cations influence reaction rates by modifying the strength of the electric field present at the catalyst surface, which alters the energetics of adsorbed intermediates and transition states according to their dipole moments and polarizabilities. The magnitude of this field depends on how cations arrange at the electrode surface, controlled by their size, shape, solvation, and packing efficiency. Cations that can arrange more densely result in a steeper potential drop at the electrode surface and consequently a stronger electric field. Our model further identifies two criteria for observing cation effects: (1) the operating potential must be negative of the electrode’s potential of zero total charge, ensuring that cations accumulate at the interface, and (2) the energetics of the kinetically relevant elementary step must be field sensitive. This framework reconciles previously inconsistent trends, including why cation effects appear only for some catalysts, why reaction selectivity is sensitive to cation identity, and why activity can increase with cation size on certain metals but decrease on others. Supported by kinetic measurements, spectroscopy, and atomistic simulations, the model provides both conceptual value for building intuition about catalysis at charged interfaces and predictive value for anticipating trends for new reactions, catalysts, and electrolytes. We conclude by highlighting the importance of electric fields across electrochemical, thermochemical, and biological catalysis and propose that considering the electrostatic environment around active sites offers new opportunities for improving activity and selectivity.

cation effects↗
DOE OSTI2025

Sliding dynamics of skyrmion molecular crystals

Using both atomistic and particle-based simulations, we investigate the current-driven dynamics of skyrmions on two-dimensional periodic substrates when there are multiple skyrmions per substrate minimum. At zero drive, the system forms pinned skyrmion molecular crystal states consisting of dimers, trimers, or dimer-trimer mixtures that have both positional and orientational order. On a square substrate lattice, the motion above depinning occurs via a running soliton that travels completely transverse to the applied current. This motion is generated by a torque from the Magnus force, which rotates the n-mer states perpendicular to the applied current. At higher drives, the flow becomes disordered while the Hall angle diminishes and gradually approaches the intrinsic value. In some cases, we also find directional locking where the Hall angle becomes locked to certain symmetry directions of the substrate over a range of currents. The transitions into and out of directionally locked states are accompanied by negative differential mobility in which the net velocity decreases as the drive increases. In conclusion, on a triangular substrate, we find no transverse mobility effects, but still observe directionally locked motion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shapiro steps and stability of skyrmions interacting with alternating anisotropy under the influence of ac and dc drives

Here we use atomistic simulations to examine the sliding dynamics of a skyrmion in a two-dimensional system containing a periodic one-dimensional stripe pattern of variations between low and high values of the perpendicular magnetic anisotropy. The skyrmion changes in size as it crosses the interface between two anisotropy regions. On applying combined dc and ac driving in either parallel or perpendicular directions, we observe a wide variety of Shapiro steps, Shapiro spikes, and phase-locking phenomena. The phase-locked orbits have two-dimensional dynamics due to the gyrotropic or Magnus dynamics of the skyrmions and are distinct from the phase-locked orbits found for strictly overdamped systems. Along a given Shapiro step when the ac drive is perpendicular to the dc drive, the velocity parallel to the ac drive is locked while the velocity in the perpendicular direction increases with increasing drive to form Shapiro spikes. At the transition between adjacent Shapiro steps, the parallel velocity jumps up to the next step value, and the perpendicular velocity drops. The skyrmion Hall angle shows a series of spikes as a function of increasing dc drive, where the jumps correspond to the transition between different phase-locked steps. At high drives, the Shapiro steps and Shapiro spikes are lost. When both the ac and dc drives are parallel to the stripe periodicity direction, Shapiro steps appear, while if the dc drive is parallel to the stripe periodicity direction and the ac drive is perpendicular to the stripe periodicity, then there are only two locked phases, and the skyrmion motion consists of a combination of sliding along the interfaces between the two anisotropy values and jumping across the interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Effects of N -glycans on the structure of human IgA2

The transition of IgA antibodies into clinical development is crucial because they have the potential to create a new class of therapeutics with superior pathogen neutralization, cancer cell killing, and immunomodulation capacity compared to IgG. However, the biological role of IgA glycans in these processes needs to be better understood. This study provides a detailed biochemical, biophysical, and structural characterization of recombinant monomeric human IgA2, which varies in the amount/locations of attached glycans. Monomeric IgA2 antibodies were produced by removing the N-linked glycans in the CH1 and CH2 domains. The impact of glycans on oligomer formation, thermal stability, and receptor binding was evaluated. In addition, we performed a structural analysis of recombinant IgA2 in solution using Small Angle X-Ray Scattering (SAXS) to examine the effect of glycans on protein structure and flexibility. Our results indicate that the absence of glycans in the Fc tail region leads to higher-order aggregates. SAXS, combined with atomistic modeling, showed that the lack of glycans in the CH2 domain results in increased flexibility between the Fab and Fc domains and a different distribution of open and closed conformations in solution. When binding with the Fcα-receptor, the dissociation constant remains unaltered in the absence of glycans in the CH1 or CH2 domain, compared to the fully glycosylated protein. These results provide insights into N-glycans’ function on IgA2, which could have important implications for developing more effective IgA-based therapeutics in the future.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. Here we isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗