Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic simulation”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗
DOE OSTI2026

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Toward the Prediction and Control of Glass Transition Temperature for Donor–Acceptor Polymers

Semiconducting donor–acceptor (D–A) polymers have attracted considerable attention toward the application of organic electronic and optoelectronic devices. However, a rational design rule for making semiconducting polymers with desired thermal and mechanical properties is currently lacking, which greatly limits the development of new polymers for advanced applications. Here, polydiketopyrrolopyrrole (PDPP)-based D–A polymers with varied alkyl side-chain lengths and backbone moieties are systematically designed, followed by investigating their thermal and thin film mechanical responses. The experimental results show a reduction in both elastic modulus and glass transition temperature (Tg) with increasing side-chain length, which is further verified through coarse-grained molecular dynamics simulations. Informed from experimental results, a mass-per-flexible bond model is developed to capture such observation through a linear correlation between Tg and polymer chain flexibility. Using this model, a wide range of backbone Tg over 80 °C and elastic modulus over 400 MPa can be predicted for PDPP-based polymers. This study highlights the important role of side-chain structure in influencing the thermomechanical performance of conjugated polymers, and provides an effective strategy to design and predict Tg and elastic modulus of future new D–A polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interactions between Per- and Polyfluoroalkyl Substances (PFAS) at the Water–Air Interface

Per- and polyfluoroalkyl substances (PFAS)--so-called “forever chemicals”--contaminate the drinking water of about 100 million people in the U.S. alone and are inefficiently removed by standard treatment techniques. A key property of these compounds that underlies their fate and transport and the efficacy of several promising remediation approaches is that they accumulate at the water−air interface. This phenomenon remains incompletely understood, particularly under conditions relevant to natural and treatment systems where water−air interfaces often carry significant loads of other organic contaminants or natural organic matter. To understand the impact of organic loading on PFAS adsorption, we carried out molecular dynamics simulations of PFAS at varying interfacial densities. We find that adsorbed PFAS form strong mutual interactions (attraction between perfluoroalkyl chains and electrostatic interactions among charged head groups) that give rise to ordered interfacial coatings. Furthermore, these interactions often involve near-cancellation of hydrophobic attraction and Coulomb repulsion. Our findings explain an apparent paradox whereby PFAS adsorption isotherms often suggest minimal mutual interactions while simultaneously displaying a high sensitivity to the composition and density of interfacial coatings. Consideration of the compounds present with PFAS at the interface has the potential to allow for more accurate predictions of fate and transport and the design of more efficient remediation approaches.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Molecular Origins of Temperature-Dependent X-ray Absorption in YAG:Dy

The absorption and emission of X-rays in dysprosium-doped yttrium aluminum garnet (YAG:Dy) has produced unexpected thermographic behavior, which is investigated using a combination of finite temperature ab initio molecular dynamic simulations, structural characterization, and electronic structure calculations of X-ray characteristics. Calculated average peak X-ray absorption spectra (XAS) from simulations between 300 and 600 K result in peak intensity loss due to thermalization effects, matching experimentally measured behavior of YAG:Dy. Here, investigation of atomic snapshots indicates structural factors that correlated with the X-ray behavior, with the first Y–O coordination sphere identified as the primary structural feature unique to high XAS intensity as calculated by radial and pair distribution functions.

Ab initio molecular dynamics↗
DOE OSTI2024

A strain density function to analyze particle size effects during high velocity impacts of yttria

Here, the impact behavior of a single yttria (Y 2 O 3 ) nanoparticle onto a Y 2 O 3 substrate is studied as a function of particle velocity (300–1200 m/s) and diameter (12–50 nm) using molecular dynamics simulations. To analyze the results, a strain density function is developed to provide both quantitative and qualitative information about the deformation mechanisms that contribute to the final state of the system. This function provides both clear evidence of shear localization and insight into the conditions for which localization occurs during the impact of Y 2 O 3 particles as they approach experimentally relevant sizes. Further analysis shows that localization and fragmentation only occur for Y 2 O 3 impacts with diameters of at least 50 nm, while particles with diameters of 25 nm and smaller deform primarily by amorphization and viscous flow. Implications for the deposition of films of Y 2 O 3 and other rare-earth oxides via the micro cold spray process are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Processes at lithium-hydride/deuteride surfaces upon low energy impact of H/D

Sputtering, reflection, and retention processes at amorphous and crystalline lithium hydride surfaces due to impact of low energy (1–100 eV) hydrogen and deuterium atoms over the range of 0 o −85 o angle of incidence at 300 K surface temperature were investigated by atomistic computational methods. Classical molecular dynamics simulations were performed with improved reactive bond-order force field (ReaxFF) potentials that include long-range polarization effects. In addition to probabilities of surface processes, the energy and angular spectra of ejected particles were obtained. Comparison of these results with those previously obtained on pristine lithium surfaces indicates the importance of saturation of the Li surface and near-surface region with hydrogen. We show that such saturation, which is typical in both laboratory and fusion device experiments with lithium coating of the plasma-facing surfaces, significantly changes the surface processes with hydrogen irradiation in the understudied low-energy region of impact energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Revealing the structure and electronic characteristics of Te-rich threshold switching materials for high-density integration

Ovonic threshold switching (OTS) selectors play an important role in the integration of advanced three-dimensional memory. Selectors based on tellurium (Te)-containing materials exhibit significant promise due to their low threshold voltages and superior consistency. Here we have theoretically studied the structure and electronic properties of a typical OTS material, amorphous GeTe 6 , to explore the switching mechanisms using ab initio molecular dynamics simulations. The results indicate that Ge atoms tend to bond with Te atoms, forming stable chemical bonds. The Te-centered clusters are predominantly in the form of distorted octahedrons, while the Ge-centered clusters are in the form of both octahedrons and tetrahedrons. Notably, the proportion of tetrahedrons within the 4-coordinated Ge-centered clusters reaches an impressive 66.9%. These tetrahedrons are randomly dispersed throughout the simulated cell, leading to a stable amorphous configuration. The mid-gap state observed in the mobility bandgap originates from the atomic chain composed of both over-coordinated Ge and Te atoms. It is the inherent stability of the chemical environment within amorphous GeTe 6 that enables it to maintain its amorphous phase under a repeated threshold voltage, a characteristic that distinguishes it from non-OTS materials, such as amorphous tellurium. Furthermore, our findings provide an in-depth understanding of the structure and electronic characteristics of amorphous GeTe 6 , which can promote the design and application of the Te-rich threshold switching materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Intermolecular correlations of liquid and glassy CS 2 studied by synchrotron radiation x-ray diffraction

How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .

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

Visualizing Anisotropic Oxygen Diffusion in Ceria under Activated Conditions

Oxygen reactivity plays a key role in the performance of ceria-based catalysts. Aberration-corrected transmission electron microscopy and molecular dynamics simulations were used to study the oxygen atom diffusion in ceria under activated conditions. Reactive oxygen atom and its real-time diffusion were visualized. Additionally, the interplay between cerium and oxygen atoms originating from a Coulomb interaction was revealed by the out-of-plane buckling of cerium atoms associated with oxygen transport. Anisotropic oxygen atom diffusion that depends on crystal orientations was discovered, demonstrating a preferential [001] crystallographic diffusion pathway. These findings reveal prospects for applications of anisotropic orientation-relevant fluorite-structured oxides.

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

Bulk viscosity of the rigid rotor one-component plasma

Bulk viscosity of a plasma consisting of strongly coupled diatomic ions is computed using molecular dynamics simulations. The simulations are based on the rigid rotor one-component plasma, which is introduced as a model system that adds two degrees of molecular rotation to the traditional one-component plasma. It is characterized by two parameters: the Coulomb coupling parameter, Γ, and the bond length parameter, Ω. Results show that the long-range nature of the Coulomb potential can lead to long rotational relaxation times, which in turn yield large values for bulk viscosity. Here, the bulk-to-shear viscosity ratio is found to span from small to large values depending on the values of Γ and Ω. Although bulk viscosity is often neglected in plasma modeling, these results motivate that it can be large in molecular plasmas with rotational degrees of freedom.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the role of methyl groups in the molecular architectures of mesophase pitches

The role of methyl groups on the liquid–crystal structure of mesophase pitches was investigated by combining experimental characterizations and atomic-scale computational modeling, using three pitches synthesized from different precursors. Of the three pitches, C-9 alkyl benzene and naphthalene-based pitches have 13 and 7 methyl groups per average polyaromatic hydrocarbon, respectively. By contrast, mesophase produced from a coal-tar pitch has about one methyl group. The coal tar–based mesophase pitch is hydrogen deficient or more aromatic compared with C-9 alkyl benzene- and naphthalene-based pitches. Additionally, X-ray diffraction data showed that average coherent domain sizes of C-9 alkyl benzene (3.7 nm) and naphthalene-based (3.6 nm) pitches with more methyl groups are larger than that of coal tar–based mesophase (2.4 nm). Based on the identified features, the influence of the methyl group on the layering structures was investigated via molecular dynamics simulations. The results revealed that methyl groups are critical in mesophase layering in C-9 alkyl benzene- and naphthalene-based pitches, by reducing CH-π interaction. However, similar alignment could be achieved without the same degree of methyl substitutions for the coal tar-based pitch because of stronger π-π interaction than the other precursors. The insights from this study contribute to our understanding of the formation of conventional mesophase pitch and have implications for the processing of coal-derived materials. In conclusion, this knowledge is vital to produce valuable products like carbon fiber and graphite from pitches.

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

Exploring the relation between transonic dislocation glide and stacking fault width in FCC metals

Theory predicts limiting gliding velocities that dislocations cannot overcome. Computational and recent experiments have shown that these limiting velocities are soft barriers and dislocations can reach transonic speeds in high rate plastic deformation scenarios. In this paper we systematically examine the mobility of edge and screw dislocations in several face centered cubic (FCC) metals (Al, Au, Pt, and Ni) in the extreme large-applied-stress regime using molecular dynamics simulations. Our results show that edge dislocations are more likely to move at transonic velocities due to their high mobility and lower limiting velocity than screw dislocations. Importantly, among the considered FCC metals, the dislocation core structure determines the dislocation’s ability to reach transonic velocities. This is likely due to the variation in stacking fault width due to relativistic effects near the limiting velocities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect-Free Growth of Decagonal Quasicrystals around Obstacles

Quasicrystals are aperiodic solids with exotic physical properties. The mechanisms driving their growth are far from understood, particularly in the presence of rigid obstacles. Here, using in situ synchrotron x-ray tomography and molecular dynamics simulations, we investigate the interaction between decagonal quasicrystals and obstacles (shrinkage pores) during solidification of an Al 79 ⁢Co 6 ⁢Ni 15 alloy. The results reveal defect-free quasicrystal growth around the pores, independent of pore size and geometry. Our findings point to a universal feature of aperiodic solids: their ability to accommodate structural disruptions due to additional (phasonic) degrees of freedom. This Letter opens new possibilities for synthesizing large-scale, single quasicrystals from a liquid metal precursor.

alloys↗
DOE OSTI2023

Polarity-dependence of the nonlinear dielectric response in interfacial water

Molecular dynamics simulations are used to study the nonlinear dielectric responses of a confined aqueous film in a planar nanopore under perpendicular electric fields at varied voltages between confining graphene sheets. Dielectric saturation reminiscent of the bulk phase behavior is prevalent at very strong fields, whereas we observe a nonmonotonic permittivity dependence on the electric field at intermediate strengths where field-alignment and spontaneous polarization of interfacial water are of comparable magnitude. The coupling between the two effects results in distinct dielectric responses at opposite confinement walls. The normal component of both the differential dielectric constant and dielectric difference constant tensors averaged over the region closer to the wall under an incoming electric field (field pointing from the liquid to the solid phase) initially increases with the strength of the imposed field. The differential permittivity peaks at a field strength previously shown to offset the surface-induced orientation bias of hydration molecules at this wall. Further strengthening of the field results in a conventional saturation behavior. At the opposite wall (subject to outgoing field) and in the central region of the water slab, the nonlinear dielectric response resembles bulklike saturation. The conditions at the permittivity extremum coincide with the window of accelerated reorientation rates of interfacial water molecules under an incoming field we uncovered in earlier molecular dynamics analyses.

Chemistry↗
DOE OSTI2020

Mode-specific pressure effects on the relaxation of an excited nitromethane molecule in an argon bath

The vibrational and rotational mode-specific relaxations of CH 3 NO 2 with 50 kcal/mol of initial internal energy in an argon bath is computed at 300 K at pressures of 10-400 atm. This work uses archived information from our previously published [J. Chem. Phys. 142, 014303 (2015)] molecular dynamics simulations and employs our previous published [J. Chem. Phys. 151, 034303 (2019)] method for projecting time-dependent Cartesian velocities onto normal mode eigenvectors. The computed relaxations cover three types of energies: vibrational, rotational, and Coriolis. In general, rotational and Coriolis relaxations in all modes are initially fast followed by an orders of magnitude slower relaxation. For all modes, that slower relaxation rate is approximately comparable to the vibrational relaxation rate. For all three types of energies, there are small-scale mode-to-mode variations. Of particular prominence is the exceptionally fast relaxation shared in common by the external rotation about the C-N axis, the internal hindered rotation of the CH 3 group relative to the NO 2 group, and the symmetric stretch of the CH 3 group.

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

Elucidating the water–anatase TiO 2 (101) interface structure using infrared signatures and molecular dynamics

The structure and dynamics of water on solid surfaces critically affect the chemistry of materials in ambient and aqueous environments. Here, we investigate the hydrogen bonding network of water adsorbed on the majority (101) surface of anatase TiO 2 , a widely used photocatalyst, using polarization- and azimuth-resolved infrared spectroscopy combined with neural network potential molecular dynamics simulations. Our results show that one monolayer of water saturates the undercoordinated titanium (Ti 5c ) sites, forming one-dimensional chains of molecule hydrogen bonded to surface undercoordinated bridging oxygen (O 2c ) atoms. As the coverage increases, water adsorption on O 2c sites leads to significant restructuring of the water monolayer and the formation of a two-dimensional hydrogen bond network characterized by tightly bound pairs of water molecules on adjacent Ti 5c and O 2c sites. This structural motif likely persists at ambient conditions, influencing the reactions occurring there. In conclusion, the results reported here provide critical details of the structure of the water–anatase (101) interface that were previously hypothesized but unconfirmed experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗