Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic dynamics”

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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Effects of breaking gravity waves on the chemical composition of the mesosphere and lower thermosphere

Eddy diffusion and momentum forcing in the mesosphere and lower thermosphere arise principally from the effects of breaking gravity and waves. The propagation and dissipation of gravity waves depend strongly on season and latitude because of its relationship to the background zonal wind. These processes have significant effects on the transport of chemical species in that region, so that observations of variations in chemical constituents provide useful tracers for dynamical theory. A parameterization of gravity wave propagation and dissipation was incorporated in the coupled dynamical-chemical model to study these effects. Several easily observable airglow features of the mesosphere-thermosphere region are shown to provide sensitive indications of gravity wave influence. In particular, theoretical predictions of large seasonal and latitudinal variations in the atomic oxygen green line, the Meinal bands of excited OH, and the O3 densities inferred from O2 (delta supra 1 sub g) emission at altitudes from 75 to 90 km are shown to be in remarkably good agreement with observations. Implications of these sensitive chemical tracers for the dynamics of the mesosphere-lower thermosphere region and its relationship to gravity waves are discussed.

Solomon, S.↗
DOE OSTI2026

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗
DOE OSTI2020

Broadening the Gas Separation Utility of Monolayer Nanoporous Graphene Membranes by an Ionic Liquid Gating

Ultrathin two-dimensional (2D) monolayer atomic crystal materials offer great potential for extending the field of novel separation technology due to their infinitesimal thickness and mechanical strength. One difficult and ongoing challenge is to perforate the 2D monolayer material with subnanometer pores with atomic precision for sieving similarly sized molecules. Here, we demonstrate the exceptional separation performance of ionic liquid (IL)/graphene hybrid membranes for challenging separation of CO 2 and N 2 . Notably, the ultrathin ILs afford dynamic tuning of the size and chemical affinity of nanopores while preserving the high permeance of the monolayer nanoporous graphene membranes. The hybrid membrane yields a high CO 2 permeance of 4000 GPU and an outstanding CO 2 /N 2 selectivity up to 32. Finally, this rational hybrid design provides a universal direction for broadening gas separation capability of atomically thin nanoporous membranes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Tailoring the Local Environment of Platinum in Single‐Atom Pt 1 /CeO 2 Catalysts for Robust Low‐Temperature CO Oxidation

Abstract A single‐atom Pt 1 /CeO 2 catalyst formed by atom trapping (AT, 800 °C in air) shows excellent thermal stability but is inactive for CO oxidation at low temperatures owing to over‐stabilization of Pt 2+ in a highly symmetric square‐planar Pt 1 O 4 coordination environment. Reductive activation to form Pt nanoparticles (NPs) results in enhanced activity; however, the NPs are easily oxidized, leading to drastic activity loss. Herein we show that tailoring the local environment of isolated Pt 2+ by thermal‐shock (TS) synthesis leads to a highly active and thermally stable Pt 1 /CeO 2 catalyst. Ultrafast shockwaves (>1200 °C) in an inert atmosphere induced surface reconstruction of CeO 2 to generate Pt single atoms in an asymmetric Pt 1 O 4 configuration. Owing to this unique coordination, Pt 1 δ+ in a partially reduced state dynamically evolves during CO oxidation, resulting in exceptional low‐temperature performance. CO oxidation reactivity on the Pt 1 /CeO 2 _TS catalyst was retained under oxidizing conditions.

Jiang, Dong↗
DOE OSTI2022

Accurate Fe–He machine learning potential for studying He effects in BCC-Fe

Nanostructured Ferritic Alloys (NFAs) containing oxide nanoparticles are prospective advanced structural materials in future fusion reactors. Lacking fusion neutron sources with an adequate flux, modeling is critical to understand radiation damage in the NFAs, including helium effects. Machine learning interatomic potentials (MLPs) have been demonstrated to be capable of describing atomic interactions in chemically complex systems with comparable accuracy to density functional theory (DFT) but with much less computational costs. Hence, MLPs are chosen to study He bubble accumulation in chemically complex NFAs. As a preliminary step, we develop a Fe–He potential, based on ~10,000 atomic configurations. The developed MLP accurately predicts the bulk properties of BCC-Fe compared to DFT, including the lattice constant (2.832 Å), elastic constants (c 11 = 271 GPa, c 12 = 141, c 44 = 93), and phonon frequencies (maximum error < 3.5%). The MLP predicts He binding energies in He n V bubbles and He n clusters with a mean absolute error (MAE) of only ~70 meV, which is ~3 to ~5 times smaller than the MAE of binding energies estimated using empirical potentials. The MLP also accurately predicts the migration barrier of interstitial He (64 meV). Furthermore, the MLP correctly predicts the maximum number of He atoms (critical size) at which a He n V bubble (He n cluster) can kick-out a self-interstitial atom (SIA), forming a bigger He n V 2 bubble (He n V bubble). Subsequently, we use the MLP to perform molecular dynamics simulations to investigate the effect of temperature on the critical size of He n V bubbles and He n clusters. We find that the critical size decreases with increasing temperature — indicating that smaller bubbles and clusters can kick-out an SIA at higher temperatures, thereby escalating the He bubble growth. We report the current Fe–He MLP can be further developed to include all chemical interactions in NFAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Localized Phonon Densities of States at Grain Boundaries in Silicon

Since it is now possible to record vibrational spectra at nanometer scales in the electron microscope, it is of interest to explore whether extended defects in crystals such as dislocations or grain boundaries will result in measurable changes of the phonon densities of states (dos) that are reflected in the spectra. Phonon densities of states were calculated for a set of high angle grain boundaries in silicon. The boundaries are modeled by supercells with up to 160 atoms, and the vibrational densities of states were calculated by taking the Fourier transform of the velocity–velocity autocorrelation function from molecular dynamics simulations with larger supercells doubled in all three directions. In selected cases, the results were checked on the original supercells by comparison with the densities of states obtained by diagonalizing the dynamical matrix calculated using density functional theory. Near the core of the grain boundary, the height of the optic phonon peak in the dos at 60 meV was suppressed relative to features due to acoustic phonons that are largely unchanged relative to their bulk values. This can be attributed to the variation in the strength of bonds in grain boundary core regions where there is a range of bond lengths.

Materials Science↗
DOE OSTI2023

Incommensurate Spiral Spin Order in CaMn 2 Bi 2 Observed via High-Pressure Neutron Diffraction

Here, high-pressure neutron diffraction is employed to investigate the magnetic behavior of CaMn 2 Bi 2 in extreme conditions. In contrast to antiferromagnetic ordering on Mn atoms reported at ambient pressure, our results reveal that at high pressure, incommensurate spiral spin order emerges due to the interplay between magnetism on the Mn atoms and strong spin–orbit coupling on the Bi atoms: sinusoidal spin order is observed at pressures as high as 7.4 GPa. First-principles calculations with a noncollinear spin orientation demonstrate band crossing behavior near the Fermi level as a result of strong hybridization between the d orbitals of Mn and the p orbitals of Bi atoms. Competing antiferromagnetic order is observed at different temperatures in the partially frustrated lattice. Theoretical models have been developed to investigate spin dynamics. This research provides a unique toolbox for conducting experimental and theoretical magnetic and spin dynamics studies of magnetic quantum materials via high-pressure neutron diffraction.

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

Multipole Expansion of Atomic Electron Density Fluctuation Interactions in the Density-Functional Tight-Binding Method

The accuracy of the density-functional tight-binding (DFTB) method in describing noncovalent interactions is limited due to its reliance on monopole-based spherical charge densities. In this study, we present a multipole-extended second-order DFTB (mDFTB2) method that takes into account atomic dipole and quadrupole interactions. Furthermore, we combine the multipole expansion with the monopole-based third-order contribution, resulting in the mDFTB3 method. To assess the accuracy of mDFTB2 and mDFTB3, we evaluate their performance in describing noncovalent interactions, proton transfer barriers, and dipole moments. Here, our benchmark results show promising improvements even when using the existing electronic parameters optimized for the original DFTB3 model. Both mDFTB2 and mDFTB3 outperform their monopole-based counterparts, DFTB2 and DFTB3, in terms of accuracy. While mDFTB2 and mDFTB3 perform comparably for neutral and positively charged systems, mDFTB3 exhibits superior performance over mDFTB2 when dealing with negatively charged systems and proton transfers. Overall, the incorporation of the multipole expansion significantly enhances the accuracy of the DFTB method in describing noncovalent interactions and proton transfers.

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

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗
DOE OSTI2021

Effects of interlayer confinement and hydration on capacitive charge storage in birnessite

Nanostructured birnessite exhibits high specific capacitance and nearly ideal capacitive behaviour in aqueous electrolytes, rendering it an important electrode material for low-cost, high-power energy storage devices. The mechanism of electrochemical capacitance in birnessite has been described as both Faradaic (involving redox) and non-Faradaic (involving only electrostatic interactions). To clarify the capacitive mechanism, we characterized birnessite’s response to applied potential using ex situ X-ray diffraction, electrochemical quartz crystal microbalance, in situ Raman spectroscopy and operando atomic force microscope dilatometry to provide a holistic understanding of its structural, gravimetric and mechanical responses. These observations are supported by atomic-scale simulations using density functional theory for the cation-intercalated structure of birnessite, ReaxFF reactive force field-based molecular dynamics and ReaxFF-based grand canonical Monte Carlo simulations on the dynamics at the birnessite–water–electrolyte interface. We show that capacitive charge storage in birnessite is governed by interlayer cation intercalation. We conclude that the intercalation appears capacitive due to the presence of nanoconfined interlayer structural water, which mediates the interaction between the intercalated cation and the birnessite host and leads to minimal structural changes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermodynamic and kinetic properties of layered-CaCo2O4 for the Ca-ion batteries: a systematic first-principles study

One of the more promising directions in multivalent energy storage is systems based on Ca ion intercalation due to the potential for high voltage and capacity. A major challenge for enabling such a battery is to find cathode materials capable of fast ionic diffusion and reversible insertion of Ca ions. Here, on the basis of first-principles calculations, we have demonstrated that layered CaCo 2 O 4 exhibits favorable thermodynamic and kinetic properties that should enable topotactic Ca ion intercalation reactions. The P3-type layered Ca x Co 2 O 4 (0 < x < 1) with either of space groups of P 1 or P 2 1 / m are stable at multiple Ca concentrations and show a smooth voltage plateau higher than 3 V up to X = 0.5. The energy barriers of the single Ca ion migration are as low as 0.36 eV and 0.27 eV at the dilute and high vacancy concentration limits, respectively. Therefore, although varying the vacancy environments of the diffusing atom influences the migration barriers, they do not exceed 0.6 eV. Stochastic analysis of Ca hopping events performed by ab initio molecular dynamics (AIMD) simulation has shown that the migration barriers are lower than 0.32 eV. Therefore, the Ca diffusivity at room temperature extrapolated from the AIMD results is comparable to Li diffusivity (>10 -10 cm 2 s -1 ) in conventional Li cathode materials, suggesting the feasibility of layered Ca x Co 2 O 4 as multivalent cathode materials. Finally, the structural factors that enable fast diffusion are discussed.

25 ENERGY STORAGE↗
DOE OSTI2024

Unconventional gas-phase synthesis of biphenyl and its atropisomeric methyl-substituted derivatives

The biphenyl molecule (C 12 H 10 ) acts as a fundamental molecular backbone in the stereoselective synthesis of organic materials due to its inherent twist angle causing atropisomerism in substituted derivatives and in molecular mass growth processes in circumstellar environments and combustion systems. Here, we reveal an unconventional low-temperature phenylethynyl addition–cyclization–aromatization mechanism for the gas-phase preparation of biphenyl (C 12 H 10 ) along with ortho-, meta-, and para-substituted methylbiphenyl (C 13 H 12 ) derivatives through crossed molecular beams and computational studies providing compelling evidence on their formation via bimolecular gas-phase reactions of phenylethynyl radicals (C 6 H 5 CC, X 2 A 1 ) with 1,3-butadiene-d 6 (C 4 D 6 ), isoprene (CH 2 C(CH 3 )CHCH 2 ), and 1,3-pentadiene (CH 2 CHCHCHCH 3 ). The dynamics involve de-facto barrierless phenylethynyl radical additions via submerged barriers followed by facile cyclization and hydrogen shift prior to hydrogen atom emission and aromatization to racemic mixtures (ortho, meta) of biphenyls in overall exoergic reactions. Furthermore, these findings not only challenge our current perception of biphenyls as high temperature markers in combustion systems and astrophysical environments, but also identify biphenyls as fundamental building blocks of complex polycyclic aromatic hydrocarbons (PAHs) such as coronene (C 24 H 12 ) eventually leading to carbonaceous nanoparticles (soot, grains) in combustion systems and in deep space thus affording critical insight into the low-temperature hydrocarbon chemistry in our universe.

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

Element-specific ultrafast lattice dynamics in FePt nanoparticles

Light–matter interaction at the nanoscale in magnetic alloys and heterostructures is a topic of intense research in view of potential applications in high-density magnetic recording. While the element-specific dynamics of electron spins is directly accessible to resonant x-ray pulses with femtosecond time structure, the possible element-specific atomic motion remains largely unexplored. We use ultrafast electron diffraction (UED) to probe the temporal evolution of lattice Bragg peaks of FePt nanoparticles embedded in a carbon matrix following excitation by an optical femtosecond laser pulse. The diffraction interference between Fe and Pt sublattices enables us to demonstrate that the Fe mean square vibration amplitudes are significantly larger that those of Pt as expected from their different atomic mass. Both are found to increase as energy is transferred from the laser-excited electrons to the lattice. Contrary to this intuitive behavior, we observe a laser-induced lattice expansion that is larger for Pt than for Fe atoms during the first picosecond after laser excitation. This effect points to the strain-wave driven lattice expansion with the longitudinal acoustic Pt motion dominating that of Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Optical properties of Platonic clusters of plasmonic nanoparticles

In the framework of the dipole approximation, we have developed a model of optical properties of a meta-atom consisting of spherical nanoparticles located at the vertices of Platonic solids. Based on the model, we have found and analysed the dynamics of changes in the optical spectra with a change in the length of the edge of a polyhedron. We have observed strong hybridisation and splitting of initially degenerate modes of individual nanoparticles. The obtained results can be used as the basis for the development of an optical nanosensor, which can determine the change in the chemical and biological composition of the environment. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Driving and characterizing nucleation of urea and glycine polymorphs in water

Crystal nucleation is relevant across the domains of fundamental and applied sciences. However, in many cases, its mechanism remains unclear due to a lack of temporal or spatial resolution. To gain insights into the molecular details of nucleation, some form of molecular dynamics simulations is typically performed; these simulations, in turn, are limited by their ability to run long enough to sample the nucleation event thoroughly. To overcome the timescale limits in typical molecular dynamics simulations in a manner free of prior human bias, here, we employ the machine learning-augmented molecular dynamics framework “reweighted autoencoded variational Bayes for enhanced sampling (RAVE).” We study two molecular systems—urea and glycine—in explicit all-atom water, due to their enrichment in polymorphic structures and common utility in commercial applications. From our simulations, we observe multiple back-and-forth nucleation events of different polymorphs from homogeneous solution; from these trajectories, we calculate the relative ranking of finite-sized polymorph crystals embedded in solution, in terms of the free-energy difference between the finite-sized crystal polymorph and the original solution state. We further observe that the obtained reaction coordinates and transitions are highly nonclassical.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗