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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

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

Charge accumulation and solvation in $β$-NiOOH: Surface chemistry of an OER catalyst from ML-aided simulations

Electrochemical water splitting is a key technology for a sustainable energy transition, providing a route to store surplus electricity from renewable sources. A central bottleneck is the sluggish oxygen evolution reaction (OER), which drives the search for catalysts that are active, stable, and inexpensive enough for large-scale deployment. Within this context, pure and doped NiO x H y combine high activity with low cost, making them prime candidates for alkaline OER. Yet, despite extensive study, the atomistic structure of NiOOH under operando conditions and the associated reaction mechanisms remain debated. Here, we investigate the structural complexity of pure β-NiOOH, the scaffold for its doped derivatives. We systematically investigate the oxidation of the surface adsorbates via proton-coupled electron transfer steps across relevant facets and sites, identifying the most probable sequence of deprotonation events. Our results reveal asymmetric charge accumulation on Wulff-relevant surfaces and show how applied potential can promote morphological restructuring. Explicit solvation is included through machine-learning interatomic potential molecular dynamics of the NiOOH/water interface, which allows us to resolve the hydrophobic and hydrophilic character of different surfaces and the associated interfacial water structure. Together, these insights demonstrate how surface chemistry and solvation jointly govern the stability of NiOOH and the accumulation of surface charge, with possible implications for catalytic performance.

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

The role of oxygen in automotive grade lithium-ion battery cathodes: an atomistic survey of ageing

The rising demand for high-performance lithium-ion batteries, pivotal to electric transportation, hinges on key materials like the Ni-rich layered oxide LiNi x Co y Al z O 2 (NCA) used in cathodes. The present study investigates the redox mechanisms, with particular focus on the role of oxygen in commercial NCA electrodes, both fresh and aged under various conditions (aged cells have performed >900 cycles until a cathode capacity retention of ~80%). Our findings reveal that oxygen participates in charge compensation during NCA delithiation, both through changes in transition metal (TM)–O bond hybridization and formation of partially reversible O 2 , the latter occurs already below 3.8 V vs. Li/Li + . Aged NCA material undergoes more significant changes in TM–O bond hybridization when cycling above 50% SoC, while reversible O 2 formation is maintained. Nickel is found to be redox active throughout the entire delithiation and shows a more classical oxidation state change during cycling with smaller changes in the Ni–O hybridization. By contrast, Co redox activity relies on a stronger change in Co–O hybridization, with only smaller Co oxidation state changes. The Ni–O bond displays an almost twice as large change in its bond length on cycling as the Co–O bond. The Ni–O 6 octahedra are similar in size to the Co–O 6 octahedra in the delithiated state, but are larger in the lithiated state, a size difference that increases with battery ageing. These contrasting redox activities are reflected directly in structural changes. The NCA material exhibits the formation of nanopores upon ageing, and a possible connection to oxygen redox activity is discussed. The difference in interaction of Ni and Co with oxygen provides a key understanding of the mechanism and the electrochemical instability of Ni-rich layered transition metal oxide electrodes. Our research specifically highlights the significance of the role of oxygen in the electrochemical performance of electric-vehicle-grade NCA electrodes, offering important insights for the creation of next-generation long-lived lithium-ion batteries.

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

Nanoscale-correlated octahedral rotations in BaZrO 3

BaZrO 3 is reportedly one of few perovskites that retains an ideal cubic structure down to zero Kelvin even though its ground state as predicted by first-principles calculations has remained contentious. Here, we combine electron diffraction with total neutron scattering measured on ceramic samples at cryogenic temperatures to demonstrate that below 80 K, this compound undergoes a structural change associated with the onset of correlated out-of-phase rotations of [ZrO 6 ] octahedra. Here, the change is manifested in the appearance of weak but relatively sharp superlattice reflections at ½hkl (h, k, l = odd), which are readily detectable by electron diffraction. Atomistic structural refinements from the total neutron scattering data confirmed that these reflections are associated with octahedral tilting. The observed rotations are consistent with the a 0 a 0 c - type and feature a coherence length of about 3 nm. The average structure of BaZrO 3 , as seen by Bragg diffraction, remains cubic. According to the electron diffraction data, Nb-doped BaZrO 3 undergoes a similar structural transition but in this case, the tilting couldn't be recovered from the neutron scattering, suggesting that the distortions are smaller compared to those in the pure compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ethanol-Induced Condensation and Decondensation in DNA-Linked Nanoparticles: A Nucleosome-like Model for the Condensed State

Inspired by the conventional use of ethanol to induce DNA precipitation, ethanol condensation has been applied as a routine method to dynamically tune “bond” lengths (i.e., the surface-to-surface distances between adjacent nanoparticles that are linked by DNA) and thermal stabilities of colloidal crystals involving DNA-linked nanoparticles. However, the underlying mechanism of how the DNA bond that links gold nanoparticles changes in this class of colloidal crystals in response to ethanol remains unclear. Here, we conducted a series of all-atom molecular dynamic (MD) simulations to explore the free energy landscape for DNA condensation and decondensation. Our simulations confirm that DNA condensation is energetically much more favorable under 80% ethanol conditions than in pure water, as a result of ethanol’s role in enhancing electrostatic interactions between oppositely charged species. Moreover, the condensed DNA adopts B-form in pure water and A-form in 80% ethanol, which indicates that the higher-order transition does not affect DNA’s conformational preferences. We further propose a nucleosome-like supercoiled model for the DNA condensed state, and we show that the DNA end-to-end distance derived from this model matches the experimentally measured DNA bond length of about 3 nm in the fully condensed state for DNA where the measured length is 16 nm in water. Altogether, this study provides an atomistic understanding of the mechanism underlying ethanol-induced condensation and water-induced decondensation, while our proposed nucleosome-like model allows the design of new strategies for interpreting experimental studies of DNA condensation.

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

Scalable Simulation of Pressure Gradient-Driven Transport of Rarefied Gases in Complex Permeable Media Using Lattice Boltzmann Method

Accurate representations of slip and transitional flow regimes present a challenge in the simulation of rarefied gas flow in confined systems with complex geometries. In these regimes, continuum-based formulations may not capture the physics correctly. This work considers a regularized multi-relaxation time lattice Boltzmann (LB) method with mixed Maxwellian diffusive and halfway bounce-back wall boundary treatments to capture flow at high Kn. The simulation results are validated against atomistic simulation results from the literature. We examine the convergence behavior of LB for confined systems as a function of inlet and outlet treatments, complexity of the geometry, and magnitude of pressure gradient and show that convergence is sensitive to all three. The inlet and outlet boundary treatments considered in this work include periodic, pressure, and a generalized periodic boundary condition. Compared to periodic and pressure treatments, simulations of complex domains using a generalized boundary treatment conserve mass but require more iterations to converge. Convergence behavior in complex domains improves at higher magnitudes of pressure gradient across the computational domain, and lowering the porosity deteriorates the convergence behavior for complex domains.

42 ENGINEERING↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Role of nanoscale coherent precipitates on the thermo-mechanical response of martensitic materials (Final Report)

Coherent second phases can have a profound effect on the properties and performance of martensitic materials, and we lack a comprehensive, mechanistic understanding of the underlying processes. This limits our ability to rationally design second phases to tailor the response of martensitic materials. To address this situation, this project sought to develop a mechanistic, predictive understanding of the thermal and mechanical response of martensitic materials with nanoscale coherent heterogeneities. A synergistic combination of atomistic simulations and experiments was used to relate the local properties of each phase (in particular, their free energy landscape) and the material nanostructure (volume fraction, shape and size of second phase precipitates, and defects that may disrupt coherency) to the overall materials response. The focus was on second phases expected to induce phenomena or properties not otherwise achievable: second order martensitic transformation in SMAs, ultra-low stiffness, increased control of transition temperatures, and fatigue resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Influence of tungsten doping on nonradiative electron–hole recombination in monolayer MoSe 2 with Se vacancies

Two-dimensional transition metal dichalcogenides (TMDs) are receiving significant attention due to their excellent electronic and optoelectronic properties. The material quality is greatly affected by defects that are inevitably generated during material synthesis. Focusing on chalcogenide vacancies, which constitute the most common defect, we use the state-of-the-art simulation methodology developed in our group to demonstrate that W doping of MoSe 2 with Se vacancies reduces charge carrier losses by two mechanisms. First, W doping makes the formation of double Se vacancies unfavorable, while it is favorable in undoped MoSe 2 . Second, if a Se vacancy is present, the charge carrier lifetimes are extended in the W-doped MoSe 2 . Combining ab initio real-time time-dependent density functional theory with nonadiabatic molecular dynamics, the simulations show that the nonradiative carrier losses in the presence of Se vacancies proceed by sub-10 ps electron trapping and relaxation down the manifold of trap states, followed by a 100 ps recombination of trapped electrons with free holes. The electron–vibrational energy exchange is driven by both in-plane and out-of-plane vibrational motions of the MoSe 2 layer. Additionally, the atomistic studies advance our understanding of the influence of defects on charge carrier properties in TMDs and guide improvements of material quality and development of TMD applications.

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

Atomic-Scale Dynamic Mechanisms of Embedded MoS 2 Wires

Nanowires composed of a 1:1 stoichiometry of transition metals and chalcogen ions can be fabricated from two-dimensional transition metal dichalcogenides (TMDs) by using electron beam irradiation. Wires fabricated through in situ experiments can be geometrically connected to TMD sheets in various ways, and their physical properties can vary accordingly. Understanding the structural transformation caused by electron beams is critical for designing wire-sheet structures for nanoelectronics. In this study, we report the behavior of nanowires formed inside a monolayer MoS 2 sheet by combining phase-contrast images and large-scale atomistic modeling. Here we investigate the effect of vacancies on the dynamic evolution of wires, such as rotations with different edge structures and breaking, by considering the interactions between MoS wires and MoS 2 nanosheets. The obtained insights can be applied to other monolayer TMDs to guide the behavior of TMD wires and fabricate favorable geometries for various applications.

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

Role of solvation site segmental dynamics on ion transport in ethylene-oxide based side-chain polymer electrolytes

Ion conducting capability is often imparted to polymeric materials through short polyether side-chains, and yet the impact of this graft polymer architecture on ion solvation and conduction has not been fully explored. Here, we use a combination of impedance spectroscopy, vibrational spectroscopy, and atomistic molecular dynamics (MD) to compare the conductivity, ionic interactions, and polymer dynamics in a series of graft polyether electrolytes. We find that in poly[(oligo ethylene oxide)methyl ether methacrylate] (POEM), a widely used graft polymer electrolyte, the ionic conductivity drops more than an order of magnitude as the side-chain length is decreased from nine ethylene oxide (EO) units to three. This difference in conductivity is unexplained by differences in the calorimetric glass transition temperature (T g ), which varies only slightly with side-chain length. Furthermore, through vibrational spectroscopy and MD simulations we demonstrate that both linear and graft polyethers solvate Li⁺ ions effectively and dissociate them from large counterions, irrespective of side-chain length. Li⁺ ions do, however, show preferential solvation by EO units far from the methacrylate backbone. Similarly, EO units far from the backbone show enhanced segmental dynamics, while those near the immobile methacrylate group move substantially more slowly, as quantified by bond vector autocorrelation relaxation times. This heterogeneity in both ion solvation and local segmental relaxation explains variation in ion conductivity where material-averaged properties such as T g and number of free ions fail to do so. Importantly, the ionic conductivity is dictated primarily by the segmental mobility of the EO units which form effective solvation sites, rather than system-wide dynamics.

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

Hidden phonon highways promote photoinduced interlayer energy transfer in twisted transition metal dichalcogenide heterostructures

Vertically stacked van der Waals (vdW) heterostructures exhibit unique electronic, optical, and thermal properties that can be manipulated by twist-angle engineering. However, the weak phononic coupling at a bilayer interface imposes a fundamental thermal bottleneck for future two-dimensional devices. Using ultrafast electron diffraction, we directly investigated photoinduced nonequilibrium phonon dynamics in MoS 2 /WS 2 at 4° twist angle and WSe 2 /MoSe 2 heterobilayers with twist angles of 7°, 16°, and 25°. We identified an interlayer heat transfer channel with a characteristic timescale of ~20 picoseconds, about one order of magnitude faster than molecular dynamics simulations assuming initial intralayer thermalization. Atomistic calculations involving phonon-phonon scattering suggest that this process originates from the nonthermal phonon population following the initial interlayer charge transfer and scattering. Our findings present an avenue for thermal management in vdW heterostructures by tailoring nonequilibrium phonon populations.

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

Skyrmion formation in Ni-based Janus dihalide monolayers: Interplay between magnetic frustration and Dzyaloshinskii-Moriya interaction

Here we present a comprehensive theory of the magnetic phases in monolayer nickel-based Janus dihalides (NiIBr, NiICl, NiBrCl) through a combination of first-principles calculations and atomistic simulations. Phonon band structure calculations and finite-temperature molecular dynamics simulations show that Ni-dihalide Janus monolayers are stable. The parameters of the interacting spin Hamiltonian extracted from ab initio calculations show that nickel-dihalide Janus monolayers exhibit varying degrees of magnetic frustration together with a large Dzyaloshinskii-Moriya interaction due to their inherent inversion symmetry breaking. The atomistic simulations reveal a delicate interplay between these two competing magnetic interactions in giving rise to skyrmion formation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Allosterically coupled conformational dynamics in solution prepare the sterol transfer protein StarD4 to release its cargo upon interaction with target membranes

Complex mechanisms regulate the cellular distribution of cholesterol, a critical component of eukaryote membranes involved in regulation of membrane protein functions directly and through the physiochemical properties of membranes. StarD4, a member of the steroidogenic acute regulator-related lipid-transfer (StART) domain (StARD)-containing protein family, is a highly efficient sterol-specific transfer protein involved in cholesterol homeostasis. Its mechanism of cargo loading and release remains unknown despite recent insights into the key role of phosphatidylinositol phosphates in modulating its interactions with target membranes. We have used large-scale atomistic Molecular dynamics (MD) simulations to study how the dynamics of cholesterol bound to the StarD4 protein can affect interaction with target membranes, and cargo delivery. We identify the two major cholesterol (CHL) binding modes in the hydrophobic pocket of StarD4, one near S136&S147 (the Ser-mode), and another closer to the putative release gate located near W171, R92&Y117 (the Trp-mode). We show that conformational changes of StarD4 associated directly with the transition between these binding modes facilitate the opening of the gate. To understand the dynamics of this connection we apply a machine-learning algorithm for the detection of rare events in MD trajectories (RED), which reveals the structural motifs involved in the opening of a front gate and a back corridor in the StarD4 structure occurring together with the spontaneous transition of CHL from the Ser-mode of binding to the Trp-mode. Further analysis of MD trajectory data with the information-theory based NbIT method reveals the allosteric network connecting the CHL binding site to the functionally important structural components of the gate and corridor. Mutations of residues in the allosteric network are shown to affect the performance of the allosteric connection. These findings outline an allosteric mechanism which prepares the CHL-bound StarD4 to release and deliver the cargo when it is bound to the target membrane.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

Tensor Network Path Integral Study of Dynamics in B850 LH2 Ring with Atomistically Derived Vibrations

The recently introduced multisite tensor network path integral (MS-TNPI) allows simulation of extended quantum systems coupled to dissipative media. We use MS-TNPI to simulate the exciton transport and the absorption spectrum of a B850 bacteriochlorophyll (BChl) ring. The MS-TNPI network is extended to account for the ring topology of the B850 system. Accurate molecular-dynamics-based description of the molecular vibrations and the protein scaffold is incorporated through the framework of Feynman–Vernon influence functional. To relate the present work with the excitonic picture, an exploration of the absorption spectrum is done by simulating it using approximate and topologically consistent transition dipole moment vectors. Comparison of these numerically exact MS-TNPI absorption spectra are shown with second-order cumulant approximations. Finally, the effect of temperature on both the exact and the approximate spectra is also explored.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effects of Al:Si and (Al+Na):Si Ratios on the Properties of the International Simple Glass, Part II: Structure

High-alumina containing high-level waste (HLW) will be vitrified at the Waste Treatment Plant at the Hanford Site. The resulting glasses, high in alumina, will have distinct composition-structure-property (C-S-P) relationships compared to previously studied HLW glasses. These C-S-P relationships determine the processability and product durability of glasses and therefore must be understood. The main purpose of this study is to understand the detailed structural changes caused by Al:Si and (Al+Na):Si substitutions in a simplified nuclear waste model glass (ISG, International Simple Glass) by combining experimental structural characterizations and molecular dynamics (MD) simulations. The structures of these two series of glasses were characterized by neutron total scattering and 27Al, 23Na, 29Si, and 11B solid-state nuclear magnetic resonance (NMR) spectroscopy. Additionally, MD simulations were used to generate atomistic structural models of the borosilicate glasses and simulation results were validated by the experimental structural data. Short-range (e.g., bond distance, coordination number, etc.) and medium-range (e.g., oxygen speciation, network connectivity, polyhedral linkages) structural features of the borosilicate glasses were systematically investigated as a function of the degree of substitution. The results show that bond distance and coordination number of the cation-oxygen pairs are relatively insensitive to Al:Si and (Al+Na):Si substitutions with the exception of the B-O pair. Additionally, the Al:Si substitution results in an increase of tri-bridging oxygen species, while (Al+Na):Si substitution creates non-bridging oxygen species. Charge compensator preferences were found for Si-[NBO] (Na+), [3]B-[NBO] (Na+), [4]B (mostly Ca2+), [4]Al (nearly equally split Na+ and Ca2+), and [6]Zr (mostly Ca2+). The network former-BO-network former linkages preferences were also tabulated; Si-O-Al and Al-O-Al were preferred at the expense of lower Si-O-[3]B and [3]B-O-[3]B linkages. These results provide insights on the structural origins of property changes such as glass transition temperature caused by the substitutions, providing a basis for future improvements of theoretical and computer simulation models.

Lu, Xiaonan↗
DOE OSTI2021

Surface localized magnetism in transition metal doped alumina

Alumina is a structural ceramic that finds many uses in a broad range of applications. It is widely employed in the aerospace and biomedical sectors due to its stability at high temperatures and in harsh chemical environments. Here, we show that magnetism can be induced at alumina surfaces by doping with 3d transition metals. We analyze the electronic structure, spin magnetic moments, and spin density of α-Al 2 O 3 as a function of both dopant species (Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu) and depth using first principles calculations. Our results show that all dopants, with the exception of Sc, produce magnetic moments that are concentrated to the surface of alumina with varying degrees of delocalization. It is seen that all of the dopants are at least meta-stable on the surface and must overcome an energy barrier of 0.19–1.14 eV in order to diffuse from the surface into the bulk. As a result of judiciously doping with select 3$\textit{d}$ transition metals the surface of alumina can be made magnetic. This could lead to novel applications in data storage, catalysis, and biomedical engineering through an added surface functionality.

36 MATERIALS SCIENCE↗