Search NASA⌕ Search

SEARCH · Search NASA

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Highly tail-asymmetric lipids interdigitate and cause bidirectional ordering

Phospholipids form structurally and compositionally diverse membranes. A less studied type of compositional diversity involves phospholipid tail variety. Some phospholipids contain two acyl tails which differ in length. These tail-asymmetric lipids are shown to contribute to temperature sensitivity, oxygen adaptability, and membrane fluidity. Membranes of a highly virulent intracellular bacterium, Francisella tularensis, contain highly tail-asymmetric 1-lignoceroyl-2-decanoyl-sn-glycero-3-phosphatidylethanolamine (XJPE) lipids which were previously shown to inhibit inflammatory responses in host cells. XJPE tails have unusually high asymmetry, and how they contribute to membrane properties on a molecular level is unknown. Here, we use small angle X-ray scattering and molecular dynamics simulations to investigate how varying XJPE ratios alters properties of simple membranes. Our results demonstrate that at high concentration they promote liquid-to-gel transition in otherwise liquid membranes, while at low concentration they are tolerated well, minimally altering membrane properties. In liquid membranes, XJPE lipids dynamically adopt two main conformations; with the long tail extended into the opposing leaflet or bent-back residing in its own leaflet. When added to both leaflets XJPE primarily adopts an extended confirmation, while asymmetric addition results in more bent-back orientations. The former increases tail ordering and the latter decreases it. XJPE tails adopt different conformations that induce composition- and leaflet-dependent bidirectional effect on membrane fluidity and this suggests that Francisella tularensis could use tail asymmetry to facilitate vesicle fusion and destabilize host cells. The effect of tail-asymmetric lipids on complex membranes should be further investigated to reveal the regulatory roles of high tail asymmetry.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Ab Initio Bulk Free Energy Surface of Proper Ferroelectrics

We report a systematic and accurate approach for deriving the bulk free energy surface (FES), a function of temperature, polarization, and strain, from the first-principles density functional theory (DFT) of proper ferroelectrics. The core of our approach is the metadynamics algorithm that extracts the polarization dependence of the FES from all-atom molecular dynamics simulations without an a priori ansatz. The rest of the FES is derived from the metadynamics trajectories that span the relevant phase space. We demonstrate our approach in the case of lead titanate. The errors across the phase transition, due to DFT numerics, all-atom molecular dynamics, and free energy evaluation by enhanced sampling, can be systematically controlled and are of the order of 1 meV/atom. The accuracy of the resulting ab initio FES is only limited by the adopted functional approximation of DFT.

Xie, Pinchen [Lawrence Berkeley National Laborator↗
DOE OSTI2025

Spectra of magnetic fields from electroweak symmetry breaking

We characterize magnetic fields produced during electroweak symmetry breaking by nondynamical numerical simulations based on the Kibble mechanism. The generated magnetic fields were thought to have an energy spectrum ∝ k3 for small wave numbers k, but here we show that it is actually a spectrum ∝ k4 along with characteristic fluctuations in the magnetic helicity. Using scaling results from magnetohydro- dynamics simulations for the evolution and assuming that the initial magnetic field is coherent on the electroweak Hubble scale, we estimate the magnetic field strength to be ∼10−13 G on kpc scales at the present epoch for nonhelical fields. For maximally helical fields we obtain ∼10−10 G on Mpc scales. We also give scalings of these estimates for partially helical fields.

Vachaspati, Tanmay [Arizona State University]↗
DOE OSTI2025

Spin Dynamics of Triple-𝐐 Magnetic Orderings in a Triangular Lattice: Implications for Multi-𝐐 Orderings in General Two-Dimensional Lattices

Multi-𝐐 magnetic structures on two-dimensional (2D) lattices provide a key route to realizing topological physics in 2D magnetism. A major experimental challenge is to unambiguously confirm their formation by excluding the possibility of topologically trivial multidomain single- or double-𝐐 magnetic orders, which cannot be distinguished using conventional diffraction techniques. Here, we propose that long-wavelength spin dynamics offers a universal diagnostic for triangular lattices: Triple-𝐐 orders that preserve rotational symmetry and single- or double-𝐐 orders that break it exhibit qualitatively distinct anisotropies in their Goldstone-mode velocities, stemming from fundamental differences in their underlying spin configurations. We validate this concept using the metallic triangular-lattice antiferromagnet Co 0.325 ⁢TaS 2 , which hosts both a stripe-type single-𝐐 state and a triple-𝐐 tetrahedral ordering at different temperatures. Using inelastic neutron-scattering and spin dynamics simulations, we first refine the spin Hamiltonian by fitting the paramagnetic excitation spectra, allowing us to develop an unbiased model independent of magnetic ordering. We then show that the observed velocity profiles of the Goldstone modes agree with the high-temperature model’s predictions: markedly anisotropic for the single-𝐐 phase and near isotropic for the triple-𝐐 phase. Importantly, this contrast persists across various exchange parameters, highlighting its model-independent nature and suggesting potential applicability to other 2D lattice systems. Beyond the long-wavelength regime, we present a substantial discrepancy between the measured and simulated magnon spectra exclusively in the triple-𝐐 phase. We attribute this discrepancy to magnon energy renormalization arising from order-of-magnitude-enhanced magnon-magnon interactions in the triple-𝐐 phase, due to its noncollinear configuration. This work provides universal insight into the dynamical properties of topological multi-𝐐 magnetic orderings in 2D lattice structures, offering a broadly applicable diagnostic to distinguishing them from topologically trivial single- or double-𝐐 counterparts. The unequivocal confirmation of the triple-𝐐 structure in Co 0.325 ⁢TaS 2 further establishes it as a prominent material platform for exploring topological spin textures in the genuine 2D limit.

Park, Pyeongjae [Oak Ridge National Laboratory (OR↗
DOE OSTI2025

Relationship between microporous structure and light gas transport through glassy polymeric membranes revealed by molecular simulations

Microporous glassy polymers are attractive materials for gas separation membranes, due to their high permeability and tailorable selectivity, provided physical aging can be delayed. The archetypal microporous glassy polymer, PTMSP, can be blended with a hyper-crosslinked isatin–triptycene porous polymer network (PPN) to delay physical aging. However, while PPN is effective at reducing physical aging, it also affects the permeability of light gases through PTMSP. Molecular dynamics simulations were used here to shed fundamental light on the mechanisms responsible for these effects. Atomistic models are developed that satisfactorily reproduce experimental observations such as matrix density and cavity size distributions for neat PTMSP as well as for PTMSP–PPN blends. Analysis of the simulation results suggests that physical aging is delayed because the PPN inclusions slow down PTMSP relaxation while reducing the connectivity between free volume pockets. To understand how PPN inclusions affect light gas permeability, the atomistic models developed are used to probe CO 2 and CH 4 diffusion and sorption. These simulations are conducted for PTMSP matrices exhibiting varying density, towards reproducing experimental permeability data. Interrogating the simulation trajectories, it is discovered that while CH 4 travels preferentially through the free volume cavities, CO 2 preferentially interacts with the available surfaces, especially in the presence of PPN. These differences help interpret experimental observations. The transport of C 2 H 6 and H 2 S through PTMSP matrices was also investigated. These gases, and in particular H 2 S, were found to absorb within PTMSP, yielding very low diffusion coefficients. The ability to predict differences in diffusion pathways for various natural gas components unveils the possibility of engineering membranes to control permeability and selectivity towards large-scale applications.

Bao Le, Tran Thi [Univ. of Oklahoma, Norman, OK (U↗
DOE OSTI2025

Biphasic solvents for post-combustion CO 2 capture from natural gas flue Gas

Fossil fuel fired power plants are generally expected to remain one of the most significant global sources of electricity for decades to come. Consequently, carbon management technologies are needed to reduce or eliminate ongoing emissions from these sources. Amongst the many techniques for carbon capture, aqueous amine-based absorbents (monoethanolamine, MEA, in particular) are, presently, considered the leading technology for post-combustion CO 2 point-source capture. These technologies are nevertheless limited by their high capital and regeneration energy costs. Biphasic solvents have been identified as an attractive alternative to traditional MEA based absorbents due to their potential energy savings. Thus far, however, the research on biphasic solvents has largely focused on their performance in coal flue gas while more dilute natural gas flue gas applications have received relatively little attention. Here, this work examines the performances of two novel biphasic solvent blends, diethylenetriamine (DETA) and triethylenetetramine (TETA), in CO 2 capture from a natural gas flue gas simulant. Across several regeneration tests, both solvents achieved considerable energy savings over the benchmark MEA solution. Specifically, the energy consumption per mol CO 2 recovered for the DETA-based and TETA-based solvents was 46 % and 35 % less than that of the benchmark MEA solution, respectively. Molecular dynamics simulations were also performed to gain a deeper understanding of the phase separation phenomena that occur as a consequence of CO 2 absorption. These simulations indicated that phase change was driven by the strong interaction between the absorption products and water, while the degree of separation depended on the CO 2 loading.

Biphasic solvents↗
DOE OSTI2025

Molecular-resolution imaging of ice crystallized from liquid water by cryogenic liquid-cell TEM

Despite the ubiquity of ice, a molecular-resolution image of nanoscopic defects or microstructures in ice crystallized from liquid water has never been obtained. This is mainly due to the difficulties in preparing and preserving crystalline ice samples that can survive under high-resolution imaging conditions. Here, we report the stabilization and Å-resolution electron imaging of ice Ih crystallized from liquid water by developing cryogenic liquid-cell transmission electron microscopy (CRYOLIC-TEM). We combine lattice mapping with molecular dynamics simulations to reveal that ice formation is highly tolerant to nanoscale defects such as misoriented subdomains and trapped gas bubbles, which are stabilized by molecular-scale structural motifs. Importantly, bubble surfaces adopt low-energy nanofacets and create negligible strain fields in the surrounding crystal. These bubbles can dynamically nucleate, grow, migrate, dissolve, and coalesce under electron irradiation and be monitored in situ near a steady state. This work improves our understanding of water crystallization behaviors at a molecular spatial resolution.

Materials science↗
DOE OSTI2025

Molecular dynamics study of grain boundaries as defect sinks under irradiation in LiAlO 2 and LiAl 5 O 8

Lithium aluminate ceramics, LiAlO 2 and LiAl 5 O 8 , show promise in nuclear environments due to their excellent radiation tolerance. Molecular dynamics simulations investigate grain boundaries (GB) and their role in defect evolution. Results reveal that GBs act as efficient defect sinks, with Li and Al atoms exhibiting distinct behaviors during displacement cascades. Tritium migration in LiAlO 2 is also studied, showing rapid diffusion and stable configurations with oxygen, corroborated by ab initio simulations from the literature. The calculated tritium diffusion coefficient of 1.33 × 10 - ¹⁴ m²/s aligns with the literature, validating the model. LiAl 5 O 8 demonstrates superior defect healing compared to LiAlO 2 , attributed to enhanced atomic transfer between grains and GBs. These findings reveal key insights into defect dynamics, providing essential insights for their application in tritium-producing burnable absorber rods (TPBARs).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Leveraging prior mean models for faster Bayesian optimization of particle accelerators

Tuning particle accelerators is a challenging and time-consuming task that can be automated and carried out efficiently using suitable optimization algorithms, such as model-based Bayesian optimization techniques. One of the major advantages of Bayesian algorithms is the ability to incorporate prior information about beam physics and historical behavior into the model used to make control decisions. In this work, we examine incorporating prior accelerator physics information into Bayesian optimization algorithms by utilizing fast executing, neural network models trained on simulated or historical datasets as prior mean functions in Gaussian process models. We show that in ideal cases, this technique substantially increases convergence speed to optimal solutions in high-dimensional tuning parameter spaces. Additionally, we demonstrate that even in non-ideal cases, where prior models of beam dynamics do not exactly match experimental conditions, the use of this technique can still enhance convergence speed. Finally, we demonstrate how these methods can be used to improve optimization in practical applications, such as transferring information gained from beam dynamics simulations to online control of the LCLS injector, and transferring knowledge gained from experimental measurements across different operating modes, such as accelerating different ion species at the ATLAS heavy ion accelerator.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Investigation of Low-Energy Lattice Dynamics and Their Role in Superionic Na Diffusion and Ultralow Thermal Conductivity of Na 3 PSe 4 as a Solid-State Electrolyte

The atomic dynamics of Na 3 PSe 4 were investigated using a combination of neutron scattering experiments and ab initio and machine-learned molecular dynamics simulations to probe the interplay of fast ionic diffusion with atomic vibrations (phonons) of the host lattice. Our results reveal the existence of low-energy vibrational modes, simultaneously involving motions of Na + ions and framework polyanion subunits, and show that these modes become strongly overdamped in the superionic regime as they couple with the Na + hopping process. In particular, the Na+ migration energy landscape is strongly impacted by low-energy phonons derived from a soft acoustic branch of the host lattice, which modulates the diameter of the Na + diffusion channel at the bottleneck. We find that an additional factor for the enhanced Na + conductivity in Na 3 PSe 4 is the presence of Na-vacancies, which also affect the low-frequency dynamics and thermal vibration amplitudes, pointing to an interplay between Na + vacancies and host dynamics, jointly enhancing ionic diffusivity. Finally, we investigate the origin of ultralow thermal conductivities in Na 3 PSe 4 and Na 3 PS 4 using Green–Kubo simulations and find that low-energy acoustic phonon modes of the overall crystal framework provide a dominant contribution to the thermal conductivity.

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

Transport coefficient approach for characterizing nonequilibrium dynamics in soft matter

Nonequilibrium states in soft condensed matter require a systematic approach to characterize and model materials, enhancing predictability and applications. Among the tools, X-ray photon correlation spectroscopy (XPCS) provides exceptional temporal and spatial resolution to extract dynamic insight into the properties of the material. However, existing models might overlook intricate details. We introduce an approach for extracting the transport coefficient, denoted as $J(t)$, from the XPCS studies. This coefficient is a fundamental parameter in nonequilibrium statistical mechanics and is crucial for characterizing transport processes within a system. Our method unifies the Green–Kubo formulas associated with various transport coefficients, including gradient flows, particle–particle interactions, friction matrices, and continuous noise. We achieve this by integrating the collective influence of random and systematic forces acting on the particles within the framework of a Markov chain. We initially validated this method using molecular dynamics simulations of a system subjected to changes in temperatures over time. Subsequently, we conducted further verification using experimental systems reported in the literature and known for their complex nonequilibrium characteristics. The results, including the derived $J(t)$ and other relevant physical parameters, align with the previous observations and reveal detailed dynamical information in nonequilibrium states. This approach represents an advancement in XPCS analysis, addressing the growing demand to extract intricate nonequilibrium dynamics. Further, the methods presented are agnostic to the nature of the material system and can be potentially expanded to hard condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2025

Ion density waves driving the formation of filamentary dust structures

The PlasmaKristall-4 (PK-4) experiment on the International Space Station allows for the study of the 3-dimensional interaction between plasma and dust particles. Previous simulations of the PK-4 environment have discovered fast moving ionization waves in the dc discharge. These ionization waves vary the plasma parameters by up to an order of magnitude, which may affect the mechanisms responsible for the self-organization of chains seen in the PK-4 experiment. Here, we adapt a molecular dynamics simulation to employ temporally varying plasma conditions in order to investigate the effect on the dust charging and electrostatic potential. In order to describe the differences between the average of the plasma conditions and the time-varying plasma condition, we present a model to reproduce the interaction that takes into account the negative potential from the dust grain and the positive potential from the ion wake.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Elucidating the Solvation Structures and Dynamics in Iron-Based TFSI – Aqueous Systems

This study investigates the solvation structures and dynamics of bis(trifluoromethanesulfonyl)imide (TFSI – )-based aqueous electrolytes, focusing on Fe(TFSI) 2 and Fe(TFSI) 3 . Using advanced characterization techniques, including small-angle X-ray scattering, molecular dynamics simulations, Raman spectroscopy, Fourier-transform infrared spectroscopy, and nuclear magnetic resonance, we elucidate how the electrolyte concentration influences ion association, solvation structures, and transport properties. Furthermore, our findings show that higher electrolyte concentrations promote the formation of contact ion pairs and anion networks, leading to reduced ion mobility and altered hydrogen-bonding dynamics. These insights provide a deeper understanding of solvation phenomena in TFSI-based electrolytes and contribute to the development of efficient and environmentally friendly iron electrodeposition processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Shapiro steps observed in a two-dimensional Yukawa solid modulated by a one-dimensional vibrational periodic substrate

Depinning dynamics of a two-dimensional (2D) solid dusty plasma modulated by a one-dimensional (1D) vibrational periodic substrate are investigated using Langevin dynamical simulations. As the uniform driving force increases gradually, from the overall drift velocity varying with the driving force, four significant Shapiro steps are discovered. The data analysis results indicate that, when the ratio of the frequency from the drift motion over potential wells to the external frequency from the modulation substrate is close to integers, dynamic mode locking occurs, corresponding to the discovered Shapiro steps. In conclusion, around both termini of the first and fourth Shapiro steps, the transitions are found to be always continuous, however, the transition between the second and third Shapiro steps is discontinuous, probably due to the different arrangements of particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Large scale polymer toughening of two-dimensional materials revealed by in situ TEM fracture tests and multiscale simulations

Two-dimensional (2D) materials offer significant potential for applications in energy-harvesting devices, batteries, sensors, and transistors. However, their intrinsic brittleness makes them prone to mechanical failure, limiting their practical use. In this work, we perform in situ transmission electron microscopy (TEM) fracture tests on monolayer MoSe2 and uncover an extrinsic toughening effect induced by an ultrathin adsorbed polystyrene adlayer. This adlayer substantially enhances the fracture resistance of the 2D flakes. Through a combination of molecular dynamics simulations and finite element analysis, we elucidate the molecular mechanism behind this toughening effect. Further, it arises from the active crack-bridging behavior of entangled polymer chains and the formation of a fracture process zone that stabilizes crack propagation and increases the energy required for crack extension. The proposed toughening mechanism offers a pathway to improving the mechanical reliability of 2D material-based devices by mitigating the risk of sudden failure.

2D materials↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗