Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerated molecular 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A generative artificial intelligence framework based on a molecular diffusion model for the design of metal-organic frameworks for carbon capture

Metal-organic frameworks (MOFs) exhibit great promise for CO 2 capture. However, finding the best performing materials poses computational and experimental grand challenges in view of the vast chemical space of potential building blocks. Here, we introduce GHP-MOFassemble, a generative artificial intelligence (AI), high performance framework for the rational and accelerated design of MOFs with high CO 2 adsorption capacity and synthesizable linkers. GHP-MOFassemble generates novel linkers, assembled with one of three pre-selected metal nodes (Cu paddlewheel, Zn paddlewheel, Zn tetramer) into MOFs in a primitive cubic topology. GHP-MOFassemble screens and validates AI-generated MOFs for uniqueness, synthesizability, structural validity, uses molecular dynamics simulations to study their stability and chemical consistency, and crystal graph neural networks and Grand Canonical Monte Carlo simulations to quantify their CO 2 adsorption capacities. We present the top six AI-generated MOFs with CO 2 capacities greater than 2m mol g -1 , i.e., higher than 96.9% of structures in the hypothetical MOF dataset.

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

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

FOXO3 targets are reprogrammed as Huntington's disease neural cells and striatal neurons face senescence with p16 INK4a increase

Neurodegenerative diseases (ND) have been linked to the critical process in aging-cellular senescence. However, the temporal dynamics of cellular senescence in ND conditions is unresolved. Here, we show senescence features develop in human Huntington's disease (HD) neural stem cells (NSCs) and medium spiny neurons (MSNs), including the increase of p16 INK4a , a key inducer of cellular senescence. We found that HD NSCs reprogram the transcriptional targets of FOXO3, a major cell survival factor able to repress cell senescence, antagonizing p16 INK4a expression via the FOXO3 repression of the transcriptional modulator ETS2. Additionally, p16 INK4a promotes cellular senescence features in human HD NSCs and MSNs. These findings suggest that cellular senescence may develop during neuronal differentiation in HD and that the FOXO3-ETS2-p16 INK4a axis may be part of molecular responses aimed at mitigating this phenomenon. Our studies identify neuronal differentiation with accelerated aging of neural progenitors and neurons as an alteration that could be linked to NDs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Double Paddle‐Wheel Enhanced Sodium Ion Conduction in an Antiperovskite Solid Electrolyte

Abstract Antiperovskite structure compounds (X 3 AB, where X is an alkali cation and A and B are anions) have the potential for highly correlated motion between the cation and a cluster anion on the A or B site. This so‐called “paddle‐wheel” mechanism may be the basis for enhanced cation mobility in solid electrolytes. Through combined experiments and modeling, the first instance of a double paddle‐wheel mechanism, leading to fast sodium ion conduction in the antiperovskite Na 3− x O 1− x (NH 2 ) x (BH 4 ), is shown. As the concentration of amide (NH 2 − ) cluster anions is increased, large positive deviations in ionic conductivity above that predicted from a vacancy diffusion model are observed. Using electrochemical impedance spectroscopy, powder X‐ray diffraction, synchrotron X‐ray diffraction, neutron diffraction, ab initio molecular dynamics simulations, and NMR, the cluster anion rotational dynamics are characterized and it is found that cation mobility is influenced by the rotation of both NH 2 − and BH 4 − species, resulting in sodium ion conductivity a factor of 10 2 higher at x = 1 than expected for the vacancy mechanism alone. Generalization of this phenomenon to other compounds could accelerate fast ion conductor exploration and design.

25 ENERGY STORAGE↗
DOE OSTI2020

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

X-ray scattering signatures of early-time accelerations in iodine dissociation

Ultrafast time-resolved x-ray scattering (TRXS) from a photoexcited molecular ensemble measures a distribution $S(Q, τ)$ of the x-ray momentum transfer Q and pump-probe delay τ in which all modes of motion induced by the excitation overlap. Frequency-resolved x-ray scattering (FRXS) based on $\tilde {S}\left(Q,\omega \right)$ separates each oscillation and dissociation channel in the TRXS data, enabling measurements of vibrational frequencies and phases, and dissociation velocities and time shifts. In this work, we extend FRXS analysis to study early-time accelerations as well. We show how these appear as diffuse scattering patterns with characteristic phase evolution in FXRS and we develop a set of transformations that isolate individual channels to measure the early-time accelerated motion. This procedure is used to analyze diatomic iodine x-ray scattering data with multiple dissociations, and the ability of this technique to characterize early-time accelerations of one dissociation channel even in the presence of another dissociation is demonstrated.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Two-component dynamics in supercritical $\text {CO}_2$ from inelastic X-ray scattering

Supercritical fluids are characterized by unique thermodynamic properties. One of these properties is the existence of two-component dynamics that is associated with distinct low-frequency and high-frequency vibrational responses of the fluid. However, the origin of this behavior remains unknown. By combining inelastic X-ray scattering and molecular dynamics simulations, we show that this behavior can be connected to density heterogeneities arising from molecular clusters. Analyses of measurements and molecular trajectories suggest that the two-component dynamics emerges due to distinct momentum fluctuations of clustered and unbound molecules. This connection between clusters and two-component dynamics highlights the importance of molecular-structural heterogeneities in supercritical fluids, colloids, and condensed-matter systems.

biophysics↗
DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2025

Control of Permanent Porosity in Type 3 Porous Liquids via Solvent Clustering

Porous liquids (PLs) are an exciting new class of materials for carbon capture due to their high gas adsorption capacity and ease of industrial implementation. They are composed of sorbent particles suspended in a nonadsorbed solvent, forming a liquid with permanent porosity. While PLs have a vast number of potential compositions based on the number of solvents and sorbent materials available, most of the research has been focused on the selection of the sorbent rather than the solvent. Therefore, PL design criteria on the supramolecular structures of the solvent are explored to create a fundamental understanding of how the solvent enables PL formation for rapid discovery of new PL compositions. Atomistic molecular dynamics simulation of eight solvents with a range of molecular sizes, shapes, and intramolecular bonding was performed, identifying that the shape and size of molecular clusters formed in the solvent are the driving predictor of PL formation rather than the size of the individual solvent molecule. The results demonstrate a significant departure from common approaches to PL formation based on the steric exclusion of solvent molecules from the sorbent via the size of the pore aperture. A modeling and experimental validation study further supports these findings. In conclusion, through this computational material design study, a previously unexplored mechanism in PL formation, solvent–solvent clustering, is identified as a critical factor for the accelerated discovery of liquid phase carbon capture materials.

Carbon capture↗
DOE OSTI2026

Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR

Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic–inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. Here, to understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31–calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ∼2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31–calcite interactions. While this work provides the first experimental confirmation that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31–calcite interface has advanced our fundamental understanding of the atomic interactions at the organic–inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.

Close, Emily G. S. [Pacific Northwest National Lab↗
DOE OSTI2025

DataSet for Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR

Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic-inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. To understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31-calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ~2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31-calcite interactions. These studies provide the first experimental evidence to confirm that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template. However, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31-calcite interface has advanced our fundamental understanding of the atomic interactions at the organic-inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.

Saccuzzo Close, Emily Grace [Pacific Northwest Nat↗
DOE OSTI2026

Gauging the Impact of Cosmic-ray Feedback on the Stellar Initial Mass Function

Cosmic rays (CRs) drive ionization and influence gas dynamics in molecular clouds (MCs), potentially impacting the resulting star formation outcomes. Although previous simulations of individual star formation have included methods for CR transport (CRT), none have been large enough to resolve the stellar initial mass function (IMF). We conduct numerical simulations following the collapse of a 20,000 M ⊙ MC and the subsequent star formation including CRT, both with and without CRs accelerated by winds from the young massive stars, and compare against a non-CRT simulation. We show that after the first massive stars form, the cavity produced by feedback is more pronounced in the CRT simulations because the external CRs are able to propagate inward and compress the gas into higher-density structures. This increases the subsequent star formation in the cloud; by the end of the simulation, the star formation efficiency (SFE) in the CRT simulation including stellar wind CRs is 43% higher than the non-CRT simulation. The IMF is also top-heavy in comparison, with a slope above 1 M ⊙ that is shallower by ∼20%. These effects are also present in the simulation without wind-accelerated CRs, but they are not as pronounced; the SFE is only 16% higher than the non-CRT simulation, and the IMF high-mass slope is shallower by ∼10%. These results may explain some of the observed top-heavy IMFs, which typically occur in high-CR environments such as the Galactic center.

Fitz Axen, Margot [Oak Ridge National Laboratory (↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effect of differently oriented interlayer phases on the radiation damage of Inconel-$\mathrm{Ni}$ multimetallic layered composite

Multimetallic layered composites (MMLCs) have shown an excellent potential for application under extreme environments, e.g., accident-tolerant fuel cladding, because of their low oxidation tendency and high corrosion resistance. Interfacial phases or complexions in nanocrystalline materials accelerate the annihilation of defects and enhance the radiation resistance of materials, making MMLCs with engineered interlayer phases compelling to deploy in extreme conditions. However, implementation of MMLCs in full capacity remained a challenge due to a lack of fundamental understanding of the underlying mechanisms governing the characteristics of the interface between the metallic layers. The precise role of interlayer phases in MMLCs and their interaction with defects, specifically under extreme conditions, is still unexplored. Pursuing atomistic simulations for various Inconel-Ni MMLCs model materials, we revealed accelerated defect mobility in interlayers with larger crystalline misorientation and the inverse relationship between the interface sink strength to the misorientation angle. Furthermore, we found a linear relation between interlayer misorientation angle with the density of radiation-induced defects and radiation enhanced displacements. So, our results indicate that radiation-induced material degradation is accelerated by the higher defect formation tendency of MMLCs with a high-angle interlayer interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Segmental relaxation of sequence defined polymers

For this work, the dynamical behavior of sequence defined polymers, P(C n EG 4 ), was studied using dielectric spectroscopy showing one segmental relaxation in addition to a secondary relaxation. In case of segmental relaxation, the relaxation times strongly depend on the molecular weight at low temperatures, while at higher temperatures, unlike to linear homo polymers, this effect levels out. With increasing length of C-units, the segmental relaxation accelerates. This is also reflected in the glass transition temperature, extracted from dielectric spectroscopy. With increasing length of C-units the glass transition temperature decreases, compatible with the accelerated segmental relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The time-resolved atomic, molecular and optical science instrument at the Linac Coherent Light Source

The newly constructed time-resolved atomic, molecular and optical science instrument (TMO) is configured to take full advantage of both linear accelerators at SLAC National Accelerator Laboratory, the copper accelerator operating at a repetition rate of 120 Hz providing high per-pulse energy as well as the superconducting accelerator operating at a repetition rate of about 1 MHz providing high average intensity. Both accelerators power a soft X-ray free-electron laser with the new variable-gap undulator section. With this flexible light source, TMO supports many experimental techniques not previously available at LCLS and will have two X-ray beam focus spots in line. Thereby, TMO supports atomic, molecular and optical, strong-field and nonlinear science and will also host a designated new dynamic reaction microscope with a sub-micrometer X-ray focus spot. The flexible instrument design is optimized for studying ultrafast electronic and molecular phenomena and can take full advantage of the sub-femtosecond soft X-ray pulse generation program.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Bidentate Lewis bases are preferred for passivation of MAPbI 3 surfaces: A time-domain ab initio analysis

We report chemical passivation is a predominant approach to inhibit intrinsic defects responsible for electron–hole recombination in perovskite solar cells. Using time-domain density functional theory combined with nonadiabatic molecular dynamics, we demonstrate defect passivation by separation of electrons and holes in metal halide perovskites (MHPs), and show that bidentate ligands exhibit the best performance. Defects in traditional semiconductors create deep midgap states, and passivation eliminates these states. In contrast, common defects produce no deep midgap states in MHPs. Instead, defects localize electrons and holes around defect sites, enhancing electron–hole interaction. Defect passivation in MHPs acts to separate charges, decreasing electron–hole and charge–phonon interactions, and increasing charge lifetimes. Bidentate ligands work best, because they can passivate both unsaturated chemical bonds created due to vacancy defects. Bidentate ligands with spatially separated binding sites are preferred, since they provide better match to the sparse inorganic lattice of MHPs. Similar to the Lewis base ligands, water also acts as a ligand and extends charge lifetimes, with the oxygen atom donating its lone electron pair to the defect site. However, water accelerates chemical degradation of MHPs. The binding energy of most ligands is larger than that of water, and therefore, the ligands displace water, increase MHP stability and prolong carrier lifetimes. The established physical mechanism of defect passivation, and the specific principles guiding the choice of passivating molecules advance our understanding of the exceptional properties of MHPs and suggest routes for further improvement of MHP performance in solar energy and optoelectronic applications.

36 MATERIALS SCIENCE↗