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 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↗
DOE OSTI2024

Origins of enhanced oxygen reduction activity of transition metal nitrides

Transition metal nitride (TMN-) based materials have recently emerged as promising non-precious-metal-containing electrocatalysts for the oxygen reduction reaction (ORR) in alkaline media. However, the lack of fundamental understanding of the oxide surface has limited insights into structure–(re)activity relationships and rational catalyst design. Here, in this work, we demonstrate how a well-defined TMN can dictate/control the as-formed oxide surface and the resulting ORR electrocatalytic activity. Structural characterization of MnN nanocuboids revealed that an electrocatalytically active Mn 3 O 4 shell grew epitaxially on the MnN core, with an expansive strain along the [010] direction to the surface Mn 3 O 4 . The strained Mn 3 O 4 shell on the MnN core exhibited an intrinsic activity that was over 300% higher than that of pure Mn 3 O 4 . A combined electrochemical and computational investigation indicated/suggested that the enhancement probably originates from a more hydroxylated oxide surface resulting from the expansive strain. This work establishes a clear and definitive atomistic picture of the nitride/oxide interface and provides a comprehensive mechanistic understanding of the structure–reactivity relationship in TMNs, critical for other catalytic interfaces for different electrochemical processes.

electrocatalysis↗
DOE OSTI2025

Origin of proton irradiation-induced deep acceptors in Al 0.70 Ga 0.30 N

Deep level defect introduction and carrier removal were characterized using steady-state photocapacitance (SSPC), deep level transient spectroscopy (DLTS) and lighted capacitance–voltage for proton irradiated n-type Al 0.70 Ga 0.30 N Schottky diodes grown by metal-organic vapor phase epitaxy on AlN-on-sapphire templates. SSPC observed deep levels in the as-grown diode with zero-phonon transition energies of 2.20, 2.65, 3.10, 3.40, and 4.65 eV relative to the conduction band minimum (E c ), and an additional deep level emerged at 1.20 eV with irradiation. Lighted capacitance–voltage measurements quantified the deep level concentration (N t ) of states detected by SSPC, and it was observed that N t increased with proton fluence only for the 1.2 and 4.65 eV levels. Carrier removal was much larger than the increase in N t of the 1.2 and 4.65 eV deep levels, suggesting that radiation-induced deep level compensators existed beyond what was detected with SSPC. DLTS detected additional, proton-induced deep acceptors at E c —0.55, 0.82, and 1.16 eV, the latter of which is likely the same 1.20 eV deep state observed by SSPC. The concentration of the E c —0.82 eV defect state was large enough to reconcile carrier removal with total deep level introduction. Comparing the E c —0.82 and 1.16 eV deep acceptor levels to previous experimental and theoretical reports suggests that their atomistic origins could be the nitrogen vacancy (V N ) and oxygen substituting on the nitrogen sub-lattice (O N ), respectively. This defect behavior contrasts starkly with GaN, where V N and O N are shallow donors, and demonstrates that the electronic properties of defects can evolve drastically within the AlGaN alloy system ranging from wide bandgap GaN to ultra-wide bandgap AlN.

Armstrong, Andrew Michael [Sandia National Laborat↗
DOE OSTI2021

Understanding copper diffusion in CuInSe 2 with first-principles based atomistic and continuum models

We investigate the diffusion of copper in CuInSe 2 using thermodynamic and kinetic models based on density functional theory calculations, attempting to reconcile large differences in reported experimental diffusivities. We find that observations of rapid chemical diffusion can be explained by large thermodynamic factors, which we calculate using a compositionally constrained model of intrinsic point defect formation. We further characterize how copper diffusion coefficients depend on material synthesis conditions and exhibit their variation across the CuInSe 2 secondary phase diagram. In doing so, we identify stable off-stoichiometries that are dominated by either vacancy- or interstitial-mediated diffusion mechanisms. These results are employed in the development of a continuum reaction–diffusion model, which we use to simulate experimental depth profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe–N–C Materials

The systematic improvement of Fe–N–C materials for fuel cell applications has proven challenging, due in part to an incomplete atomistic understanding of the oxygen reduction reaction (ORR) under electrochemical conditions. Herein, a multilevel computational approach, which combines ab initio molecular dynamics simulations and constant potential density functional theory calculations, is used to assess proton-coupled electron transfer (PCET) processes and adsorption thermodynamics of key ORR intermediates. These calculations indicate that the potential-limiting step for ORR on Fe–N–C materials is the formation of the Fe III –OOH intermediate. They also show that an active site model with a water molecule axially ligated to the iron center throughout the catalytic cycle produces results that are consistent with the experimental measurements. In particular, reliable prediction of the ORR onset potential and the Fe(III/II) redox potential associated with the conversion of Fe III –OH to Fe II and desorbed H 2 O requires an axial H 2 O co-adsorbed to the iron center. The observation of a five-coordinate rather than four-coordinate active site has significant implications for the thermodynamics and mechanism of ORR. Importantly, these findings highlight the importance of solvent–substrate interactions and surface charge effects for understanding the PCET reaction mechanisms and transition-metal redox couples under realistic electrochemical conditions.

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

Directing reaction pathways via in situ control of active site geometries in PdAu single-atom alloy catalysts

Abstract The atomic scale structure of the active sites in heterogeneous catalysts is central to their reactivity and selectivity. Therefore, understanding active site stability and evolution under different reaction conditions is key to the design of efficient and robust catalysts. Herein we describe theoretical calculations which predict that carbon monoxide can be used to stabilize different active site geometries in bimetallic alloys and then demonstrate experimentally that the same PdAu bimetallic catalyst can be transitioned between a single-atom alloy and a Pd cluster phase. Each state of the catalyst exhibits distinct selectivity for the dehydrogenation of ethanol reaction with the single-atom alloy phase exhibiting high selectivity to acetaldehyde and hydrogen versus a range of products from Pd clusters. First-principles based Monte Carlo calculations explain the origin of this active site ensemble size tuning effect, and this work serves as a demonstration of what should be a general phenomenon that enables in situ control over catalyst selectivity.

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

Facet-dependent structure and dissociation of water at pristine IrO 2 /water interfaces

Understanding the microscopic structure of water at metal oxide interfaces is crucial for advancing electrocatalysis. IrO 2 , specifically, has shown exceptional activity for electrochemical water oxidation, but we currently lack a fundamental understanding of how the surface structure of IrO 2 impacts water reactivity. In this work, we developed a machine learning potential trained to first-principles accuracy for modeling IrO 2 /water interfaces across different facets: (110), (100), (101), and (001). Using extensive machine learning molecular dynamics simulations, we investigated the spontaneous dissociation of water molecules at these interfaces. Our results reveal a distinct dissociation probability trend: (110) > (100) ≈ (101) > (001), which we attribute primarily to the reaction thermodynamics of surface water dissociation. A strong correlation is observed between the surface Ir–O bond distances and the dissociation probabilities, highlighting the role of surface geometry in modulating reactivity. As a consequence, the interfacial solvation structures and hydrogen bonding environments are dynamically tuned by the varying water dissociation capabilities across facets. This work elucidates how water dissociation energetics depend on surface orientation and interfacial structure, offering atomistic insights into manipulating reaction chemistry at electrocatalytic interfaces.

organic↗
DOE OSTI2024

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

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

Metastable piezoelectric group-IV monochalcogenide monolayers with a buckled honeycomb structure

Multiple two-dimensional materials are being naïvely termed stable on the grounds of displaying phonon dispersions with no negative frequencies and of not collapsing on molecular dynamics calculations at fixed volume. But, if these phases do not possess the smallest possible structural energy, how does one understand and establish their actual meta stability? To answer this question, twelve two-dimensional group-IV monochalcogenide monolayers (SiS, SiSe, SiTe, GeS, GeSe, GeTe, SnS, SnSe, SnTe, PbS, PbSe, and PbTe) with a buckled honeycomb atomistic structure—belonging to symmetry group P3m1—displaying an out-of-plane intrinsic electric polarization are shown to be metastable by three independent methods. First, we uncover a coordination-preserving structural transformation from the low-buckled honeycomb structure onto the lower-energy Pnm2 1 (or Pmmn for PbS, PbSe, and PbTe) phase to estimate energy barriers E B that must be overcome during such structural transformation. Using the curvature of the local minima and E B as inputs to Kramers escape formula, large escape times are found, implying the structural metastability of the buckled honeycomb phase (with the exception of PbS and PbSe, these phases display escape times ranging from 700 years to multiple times the age of the universe and can be considered “stable” for practical purposes in that relative sense). The second demonstration is provided by phonon dispersion relations that include the effect of long-range Coulomb forces and display no negative vibrational modes. The third and final demonstration of structural metastability is furnished by room-temperature ab initio molecular dynamics for selected compounds. Here, the magnitude of the electronic band gap evolves with chemical composition. Different from other binary two-dimensional compounds such as transition metal dichalcogenide monolayers and hexagonal boron nitride monolayers which only develop an in-plane piezoelectric response, the twelve group-IV monochalcogenide monolayers with a buckled honeycomb structure also display out-of-plane piezoelectric properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Revealing Local Structures through Machine-Learning-Fused Multimodal Spectroscopy

Atomistic structures of materials offer valuable insights into their functionality. Determining these structures remains a fundamental challenge in materials science, especially for systems with defects. While both experimental and computational methods exist, each has limitations in resolving nanoscale structures. Core-level spectroscopies, such as X-ray absorption (XAS) or electron energy-loss spectroscopies (EELS), have been used to determine the local bonding environment and structure of materials. Recently, machine learning (ML) methods have been applied to extract structural and bonding information from XAS/EELS data. However, frameworks relying solely on a single data stream, defined as characterization data derived from a single element using one technique, are often insufficient because multiple local environments can yield similar spectral features, making it challenging to differentiate between competing structural hypotheses. Here, in this work, we address this challenge by integrating multimodal ab initio simulations, experimental data acquisition, and ML techniques for structure characterization. Our goal is to determine local structures and properties using EELS and XAS data from multiple elements and edges. To showcase our approach, we use various lithium nickel manganese cobalt (NMC) oxide compounds which are used for lithium ion batteries, including those with oxygen vacancies and antisite defects, as the sample material system. We successfully inferred local element content, ranging from lithium to transition metals, with quantitative agreement with experimental data. Beyond local element inference, we find that ML model based on multimodal spectroscopic data is able to determine whether local defects such as oxygen vacancy and antisites are present, a task which is impossible for single mode spectra or other experimental techniques. Furthermore, our framework is able to provide physical interpretability, bridging spectroscopy with the local atomic and electronic structures.

battery↗
DOE OSTI2023

An open database of computed bulk ternary transition metal dichalcogenides

Abstract We present a dataset of structural relaxations of bulk ternary transition metal dichalcogenides (TMDs) computed via plane-wave density functional theory (DFT). We examined combinations of up to two chalcogenides with seven transition metals from groups 4–6 in octahedral (1T) or trigonal prismatic (2H) coordination. The full dataset consists of 672 unique stoichiometries, with a total of 50,337 individual configurations generated during structural relaxation. Our motivations for building this dataset are (1) to develop a training set for the generation of machine and deep learning models and (2) to obtain structural minima over a range of stoichiometries to support future electronic analyses. We provide the dataset as individual VASP xml files as well as all configurations encountered during relaxations collated into an ASE database with the corresponding total energy and atomic forces. In this report, we discuss the dataset in more detail and highlight interesting structural and electronic features of the relaxed structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗
DOE OSTI2024

Neural network interatomic potential-driven analysis of phase stability in Ti–V alloys at the atomistic scale

The evolution of the ω phase in titanium–vanadium (Ti–V) alloys is critical for their mechanical properties, particularly in aerospace and biomedical applications. Here, this study employs a Rapid Artificial Neural Network (RANN) potential to model the ω phase evolution at the atomistic level, demonstrating a high degree of consistency with experimental observations, unlike the Modified Embedded Atom Method (MEAM), which fails to capture this phase transformation accurately. RANN simulations replicate key phenomena such as the nucleation of α precipitates at ω/β interfaces and accurate lattice orientations, enhancing our understanding of phase stability and transformation kinetics. The findings affirm that RANN potentials can significantly improve the prediction accuracy of complex material behaviors, offering a powerful tool for designing advanced materials with tailored properties such as solute effect in various stacking fault energies. This approach not only bridges the gap between theoretical predictions and empirical data but also sets a new direction for future research in materials science, emphasizing the integration of machine learning techniques in the development and optimization of new alloys.

36 MATERIALS SCIENCE↗