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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Molecular dynamics simulations of energy accommodation between gases and polymers for ultra-low thermal conductivity insulation

Determining the energy accommodation between gases and solids is essential to developing porous thermal insulation materials with ultra-low effective thermal conductivity that reduce energy use, greenhouse gas emissions, and fossil fuel consumption. The energy accommodation coefficients of most gases, however, have been rarely studied, especially with respect to solids that have relatively high thermal resistivity, e.g., polymers. Here, by using all-atom nonequilibrium molecular dynamics simulations, we reveal the accommodation coefficients of He, Ar, N 2 , and O 2 with polymers, mainly polystyrene. We find that their values are around 0.51, 0.72, 0.79, and 0.90, respectively, suggesting a critical reexamination of the commonly used theoretical maximum value of 1. We have also conducted experiments and validated the value for air, which is about 0.81. Such a change in accommodation coefficients can lead to a reduction of about 70%, 50%, 35%, and 20% in the thermal conductivity of He, Ar, N 2 , and O 2 gases in nano pores (below 100 nm) or at low pressures (below 1 millibar). With these new accommodation coefficients, we find that in a 10 nm pore with ambient pressure at 300 K, the gas thermal conductivity of He, Ar, N 2 , and O 2 in porous polystyrene can be as low as 9.7 × 10 –4 , 3.4 × 10 –4 , 7.3 × 10 –4 , and 8.5 × 10 –4 W·m –1 ·K –1 , respectively, which are two to three orders of magnitude lower than their bulk values, promising higher thermal resistivity of insulation materials. This work reveals the fundamental energy exchange between gases and polymers, providing important guidance for designing high-performance thermal insulation materials for various applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Subdiffusive High-Pressure Hydrogen Gas Dynamics in Elastomers

We report elastomeric rubber materials serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. High-pressure exposure and subsequent rapid decompression often lead to cavitation and stress-induced damage of the elastomer due to localization of the hydrogen gas. Here, we use all-atom classical molecular dynamics simulations to assess the impact of compositional variations on gas diffusion within the commonly used elastomer ethylene–propylene–diene monomer (EPDM). With the aim to build a predictive understanding of precursors to cavitation and to motivate material formulations that are less sensitive to hydrogen-induced failure, we perform systematic simulations of gas dynamics in EPDM as a function of temperature, gas concentration, and cross-link density. Our simulations reveal anomalous, subdiffusive hydrogen motion at pressure and intermediate times. We identify two groups of gas with different mobilities: one group exhibiting high mobility and one group exhibiting low mobility due to their motion being impeded by the polymer. With decreasing temperatures, the low-mobility group shows increased gas localization, the necessary precursor for cavitation damage in these materials. At lower temperatures, increasing cross-link density led to greater hydrogen gas mobility and a lower fraction of caged hydrogen, indicating that increasing cross-link density may reduce precursors to cavitation. Finally, we use a two-state kinetic model to determine the energetics associated with transitions between these two mobility states.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Geometric Phase in Anisotropic Kepler Problem: Perspective for Realization in Rydberg Atoms

We predict a gyroscopic effect that can be demonstrated with Rydberg atoms following the dynamics of a Kepler Hamiltonian with an additional uniaxial anisotropy induced by optical ponderomotive force. This effect is analogous to the rotation of the Foucault pendulum in response to the Earth’s rotation. Furthermore, we argue that in Rydberg states with a large principal quantum number a similar geometric angle can be generated by mechanical rotations of an atomic-optical setup on timescales between 1 μ⁢s and 1 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Computational Modeling of Battery Materials

This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.

25 ENERGY STORAGE↗
DOE OSTI2020

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R 0 ) and fluorophore orientation factor (κ 2 ) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Machine learning bridges local static structure with multiple properties in metallic glasses

A long-standing challenge in the metallic glass (MG) community has been how to quantitatively gauge the influence of the intricate local packing environment on the response (such as the propensity for atomic rearrangement) of the atomic configuration to external stimuli. Here we establish this structure–property relation by representing the complex amorphous structure using a single, flexibility-orientated structural quantity. This structural flexibility (SF) couples to a bona fide structural representation, the pair distribution function (PDF) of individual atoms, through a weighting function that reflects what matters in the static atomic configuration to dynamic responses. Machine learning is used, employing microscopic flexibility volume as the supervisory signal, to establish via direct regression an optimized weighting vector, which is proven robust for all quenching rates, deformation conditions, and different compositions in a given (e.g., Cu x Zr 100-x ) alloy system. Additionally, the SF is evaluated solely from the particle positions (PDF), for any structure variation, from the atomic scale up to sample average. Strong correlations are demonstrated between SF and a broad range of properties, including vibrational, diffusional, as well as elastic and plastic relaxation responses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multichannel photodissociation dynamics in CS 2 studied by ultrafast electron diffraction

In this study, the structural dynamics of photoexcited gas-phase carbon disulfide (CS 2 ) molecules are investigated using ultrafast electron diffraction. The dynamics were triggered by excitation of the optically bright 1 B 2 ( 1 Σ u + ) state by an ultraviolet femtosecond laser pulse centred at 200 nm. In accordance with previous studies, rapid vibrational motion facilitates a combination of internal conversion and intersystem crossing to lower-lying electronic states. Photodissociation via these electronic manifolds results in the production of CS fragments in the electronic ground state and dissociated singlet and triplet sulphur atoms. The structural dynamics are extracted from the experiment using a trajectory-fitting filtering approach, revealing the main characteristics of the singlet and triplet dissociation pathways. Finally, the effect of the time-resolution on the experimental signal is considered and an outlook to future experiments provided.

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

Emergent interaction-induced topology in Bose-Hubbard ladders

We investigate the quantum many-body dynamics of bosonic atoms hopping in a two-leg ladder with strong on-site contact interactions. We observe that when the atoms are prepared in a staggered pattern with pairs of atoms on every other rung, singlon defects, i.e., rungs with only one atom, can localize due to an emergent topological model, even though the underlying model in the absence of interactions admits only topologically trivial states. This emergent topological localization results from the formation of a zero-energy edge mode in an effective lattice formed by two adjacent chains with alternating strong and weak hoping links (Su-Schrieffer-Heeger chains) and opposite staggering which interface at the defect position. Our findings open the opportunity to dynamically generate nontrivial topological behaviors without the need for complex Hamiltonian engineering. Published by the American Physical Society 2025

Wellnitz, David (ORCID:0000000349788938)↗
DOE OSTI2023

Bond clusters control rupture force limit in shear loaded histidine-Ni 2+ metal-coordinated proteins

Dynamic noncovalent interactions are pivotal to the structure and function of biological proteins and have been used in bioinspired materials for similar roles. Metal-coordination bonds, in particular, are especially tunable and enable control over static and dynamic properties when incorporated into synthetic materials. Despite growing efforts to engineer metal-coordination bonds to produce strong, tough, and self-healing materials, the systematic characterization of the exact contribution of these bonds towards mechanical strength and the effect of geometric arrangements is missing, limiting the full design potential of these bonds. In this work, we engineer the cooperative rupture of metal-coordination bonds to increase the rupture strength of metal-coordinated peptide dimers. Utilizing all-atom steered molecular dynamics simulations on idealized bidentate histidine-Ni 2+ coordinated peptides, we show that histidine-Ni 2+ bonds can rupture cooperatively in groups of two to three bonds. We find that there is a strength limit, where adding additional coordination bonds does not contribute to the additional increase in the protein rupture strength, likely due to the highly heterogeneous rupture behavior exhibited by the coordination bonds. Further, we show that this coordination bond limit is also found natural metal-coordinated biological proteins. Using these insights, we quantitatively suggest how other proteins can be rationally designed with dynamic noncovalent interactions to exhibit cooperative bond breaking behavior. Altogether, this work provides a quantitative analysis of the cooperativity and intrinsic strength limit for metal-coordination bonds with the aim of advancing clear guiding molecular principles for the mechanical design of metal-coordinated materials.

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

Activation of atomic transport via vibrational coupling-induced force fluctuations

Vibrational coupling, although well-explored in many fields, has seldom been investigated in the context of mass transport. In this Letter, we examine the impact of vibrational coupling on atomic transport using simple molecular dynamics simulations. Our study shows that the atomic transport can be activated when the natural frequency of the atomic slit is close to the natural frequency of the atom being transported. We uncover the presence of fluctuating forces induced by vibrational coupling with higher amplitudes observed when the coupling is strong. We show that the transport activation mechanism is due to the high force fluctuations that arise during strong vibrational coupling allowing the atom to temporarily surpass the transport barrier of the slit. Our findings will serve as a foundation for the continued examination of vibrational coupling in the realm of mass transport.

Physics↗
DOE OSTI2021

Two-dimensional overdamped fluctuations of the soft perovskite lattice in CsPbBr 3

Lead halide perovskites exhibit structural instabilities and large atomic fluctuations thought to impact their optical and thermal properties, yet detailed structural and temporal correlations of their atomic motions remain poorly understood. Here, these correlations are resolved in CsPbBr 3 crystals using momentum-resolved neutron and X-ray scattering measurements as a function of temperature, complemented with first-principles simulations. We uncover a striking network of diffuse scattering rods, arising from the liquid-like damping of low-energy Br-dominated phonons, reproduced in our simulations of the anharmonic phonon self-energy. Here, these overdamped modes cover a continuum of wave vectors along the edges of the cubic Brillouin zone, corresponding to two-dimensional sheets of correlated rotations in real space, and could represent precursors to proposed two-dimensional polarons. Further, these motions directly impact the electronic gap edge states, linking soft anharmonic lattice dynamics and optoelectronic properties. These results provide insights into the highly unusual atomic dynamics of halide perovskites, relevant to further optimization of their optical and thermal properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling the Atomic Mechanism of the Crystalline Phase‐Dependent Structural Features and Special Spectral Design of α‐, β‐, and Ɛ‐Ga₂O₃

Atomic‐scale phase transformations profoundly influence the functional properties of Ga₂O₃ polymorphs. By combining irradiation experiments with microstructure characterization and theoretical approaches, phase‐specific energy‐dissipation pathways in α‐, β‐, and ε‐Ga₂O₃ are uncovered and strategies for targeted property design are outlined. Competing antiphase boundaries (APBs) and twin domain boundaries (TDBs) promote irreversible α→ε interconversion through domain fragmentation. In β‐Ga₂O₃, defect‐induced stress gradients drive two distinct local transformations: surface Ga‐aggregated β→δ that stabilizes transient states, and latent‐track‐confined β→κ phase transition with recoverable distortions via cation reordering. Under electronic excitation, β‐Ga₂O₃ forms nanohillocks via robust GaO₆ octahedra (high density/strong Ga─O bonds), while α/ε‐Ga₂O₃ generates nanopores from tetrahedral Ga looseness (low bonding energy), highlighting phase‐dependent surface dynamics shaped by atomic packing and bonding anisotropy. Defect‐regulated recombination suppresses visible photoluminescence in α/β‐Ga₂O₃, whereas in ε‐Ga₂O₃ bandgap narrowing of ΔE: 0.30 eV is observed, enhancing emission. Linking phase‐dependent defect‐carrier interactions and metastable‐phase engineering in Ga₂O₃ enables property optimization for power‐electronics and optoelectronics devices.

electronic state configuration↗
DOE OSTI2021

Deep potential generation scheme and simulation protocol for the Li 10 GeP 2 S 12 -type superionic conductors

We report solid-state electrolyte materials with superior lithium ionic conductivities are vital to the next-generation Li-ion batteries. Molecular dynamics could provide atomic scale information to understand the diffusion process of Li-ion in these superionic conductor materials. Here, we implement the deep potential generator to set up an efficient protocol to automatically generate interatomic potentials for Li 10 GeP 2 S 12 -type solid-state electrolyte materials (Li 10 GeP 2 S 12 , Li 10 SiP 2 S 12 , and Li 10 SnP 2 S 12 ). The reliability and accuracy of the fast interatomic potentials are validated. With the potentials, we extend the simulation of the diffusion process to a wide temperature range (300 K–1000 K) and systems with large size (~1000 atoms). Important technical aspects such as the statistical error and size effect are carefully investigated, and benchmark tests including the effect of density functional, thermal expansion, and configurational disorder are performed. The computed data that consider these factors agree well with the experimental results, and we find that the three structures show different behaviors with respect to configurational disorder. Our work paves the way for further research on computation screening of solid-state electrolyte materials.

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

Quantum spin liquids bootstrapped from Ising criticality in Rydberg arrays

Arrays of Rydberg atoms constitute a highly tunable, strongly interacting venue for the pursuit of exotic states of matter. We develop a strategy for accessing a family of fractionalized phases known as quantum spin liquids in two-dimensional Rydberg arrays. We specifically use effective field theory methods to study arrays assembled from Rydberg chains tuned to an Ising phase transition that famously hosts emergent fermions propagating within each chain. This highly entangled starting point allows us to naturally access spin liquids familiar from Kitaev's honeycomb model—albeit from an entirely different framework. In particular, we argue that finite-range repulsive Rydberg interactions, which frustrate nearby symmetry-breaking orders, can enable coherent propagation of emergent fermions between the chains in which they were born. Delocalization of emergent fermions across the full two-dimensional Rydberg array yields a gapless Z 2 spin liquid with a single massless Dirac cone. Here, the Rydberg occupation numbers exhibit universal power-law correlations that provide a straightforward experimental diagnostic of this phase. We further show that explicitly breaking symmetries perturbs the gapless spin liquid into gapped, topologically ordered descendants: Breaking lattice symmetries generates toric-code topological order, whereas introducing Floquet-mediated chirality generates non-Abelian Ising topological order. In the toric-code phase, we analytically construct microscopic incarnations of non-Abelian defects, which can be created and transported by dynamically controlling the atom positions in the array. Our work suggests that appropriately tuned Rydberg arrays provide a cold-atoms counterpart of solid-state “Kitaev materials” and, more generally, it spotlights a different angle for pursuing experimental platforms for Abelian and non-Abelian fractionalization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Relevance of structural defects to the mechanism of mechanical deformation in metallic glasses

Abstract It is known that deformation in disordered materials such as metallic glasses and supercooled liquids occurs via the cooperative rearrangement of atoms or constituent particles at dynamical heterogeneities, commonly regarded as point-like defects. We show via molecular-dynamics simulations that there is no apparent relationship between atomic rearrangements and the local atomic environment as measured by the atomic-level stresses, kinetic and potential energies, and the per-atom Voronoi volume. In addition, there is only a weak correlation between atomic rearrangements and the largest and smallest eigenvalues of the dynamical matrix. Our results confirm the transient nature of dynamical heterogeneities and suggest that the notion of defects may be less relevant than that of a propensity for rearrangement.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Direct Visualization of Defect‐Controlled Diffusion in van der Waals Gaps

Abstract Diffusion processes govern fundamental phenomena such as phase transformations, doping, and intercalation in van der Waals (vdW) bonded materials. Here, the diffusion dynamics of W atoms by visualizing the motion of individual atoms at three different vdW interfaces: hexagonal boron nitride (BN)/vacuum, BN/BN, and BN/WSe 2 , by recording scanning transmission electron microscopy movies is quantified. Supported by density functional theory (DFT) calculations, it is inferred that in all cases diffusion is governed by intermittent trapping at electron beam‐generated defect sites. This leads to diffusion properties that depend strongly on the number of defects. These results suggest that diffusion and intercalation processes in vdW materials are highly tunable and sensitive to crystal quality. The demonstration of imaging, with high spatial and temporal resolution, of layers and individual atoms inside vdW heterostructures offers possibilities for direct visualization of diffusion and atomic interactions, as well as for experiments exploring atomic structures, their in situ modification, and electrical property measurements of active devices combined with atomic resolution imaging.

Chemistry↗
DOE OSTI2026

The Focus of an Ultraintense Laser Pulse

Back in 1885, von Lommel studied the diffraction of a circular opening. His study was one step further than the Airy diffraction model (describing patterns in the focal plane) because von Lommel’s model was able to describe what the light distribution before and after the focus is. Today’s extreme lasers have a flat-top profile whose diffraction pattern in the vicinity of the focus is very much reminiscent of von Lommel’s description. In the present Chapter, we will compute, based on a Huygens-Fresnel integral method, the laser field distribution around the focus. This modelling is particularly important because the intensity distribution before and after the focal plane (where the Airy disk exists) is very relevant to many contemporary experiments, and is key to the ionization of atoms and the dynamics of liberated electrons. For clarity, we are going to analyze three cases corresponding to three f-numbers: f/10, f/1.8, and f/1.1. Our computed results reveal the present-day significance of von Lommel’s pioneering results. Here, they help us understand the light distribution in the entire focal region of an ultraintense laser pulse including the full vectorial fields, not just the focal plane.

Roso, Luis [Universidad de Salamanca (Spain)]↗
DOE OSTI2025

Empirical evidence that glucan-interacting amino acid side chains within the transmembrane channel collectively facilitate cellulose synthase function

The fundamental mechanism of cellulose synthesis is widely conserved across Kingdoms and depends on cellulose synthases, which are processive, dual-function, family 2 glycosyltransferases (GT-2). These enzymes polymerize glucose on the cytoplasmic side of the plasma membrane and export the glucan chain to the cell surface through an integral transmembrane (TM) channel. Structural studies of active plant cellulose synthases (CESAs) have revealed interactions between the nascent glucan chain and the side chains of polar, charged, and aromatic amino acid residues that line the TM channel. However, the functional consequences of modifying these side chains have not been tested in vivo in CESAs or other processive GT-2s. To test this, we used an established in vivo assay based on genetic complementation of CESA5 in the moss, Physcomitrium patens. For accurate prediction of glucan-interacting amino acid residues, we generated a complete homotrimeric molecular model of PpCESA5 using a combination of homology and de novo modeling. All-atom molecular dynamics-based analyses of contact metrics and interaction energy identified 23 amino acid residues with high propensity to interact with the nascent glucan chain within the TM channel or on the apoplastic surface of PpCESA5. Mutating any one of 18 of these amino acid residues to alanine, thereby removing their side chains, abolished or impaired CESA function, with the strongest effects observed upon the loss of charged amino acid side chains. This provides direct evidence to support the hypothesis that multiple amino acid residues collectively maintain a smooth energy landscape within the TM channel to facilitate glucan translocation.

59 BASIC BIOLOGICAL SCIENCES↗