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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High Explosive Ignition through Chemically Activated Nanoscale Shear Bands

Shock initiation and detonation of high explosives is considered to be controlled through hot spots, which are local regions of elevated temperature that accelerate chemical reactions. Using classical molecular dynamics, we predict the formation of nanoscale shear bands through plastic failure in shocked TATB high explosive crystal. By scale-bridging with quantum-based molecular dynamics, we show that shear bands exhibit lower reaction barriers. While shear bands quickly cool,they remain chemically activated and support increased reaction rates without the local heating typically evoked by the hot spot paradigm. We describe this phenomenon as chemical activation through shear banding

Chemistry - Inorganic↗
DOE OSTI2024

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

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

Accelerated Sequence Design of Star Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning

Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatibilization, chemical transformations, and separations. s-BCPs have chain architectures where three or more linear diblock copolymer arms comprised of two chemically distinct linear polymers, e.g., solvophobic and solvophilic chains, are covalently joined at one point. The chemical composition of each of the subunit polymer chains comprising the arms, their molecular weights, and the number of arms can be varied to tailor the surface and interfacial activity of these architecturally unique molecules. Further, this makes identification of the optimal s-BCP design nontrivial as the total number of plausible s-BCP architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with a reinforcement learning-based Monte Carlo tree search (MCTS) to identify s-BCP designs that minimize the interfacial tension between polar and nonpolar solvents. We first validate the MCTS approach for the design of small- and medium-sized s-BCPs and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified by using the configurations obtained from MD simulations. Chemical insights into the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provides important groundwork for future experimental investigation of s-BCPs for various applications.

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

Synthetic Control of Exciton Dynamics in Bioinspired Cofacial Porphyrin Dimers

Understanding how the complex interplay among excitonic interactions, vibronic couplings, and reorganization energy determines coherence-enabled transport mechanisms is a grand challenge with both foundational implications and potential payoffs for energy science. We use a combined experimental and theoretical approach to show how a modest change in structure may be used to modify the exciton delocalization, tune electronic and vibrational coherences, and alter the mechanism of exciton transfer in covalently linked cofacial Zn-porphyrin dimers (meso-beta linked AB m-β and meso-meso linked AA m-m ). While both AB m-β and AA m-m feature zinc porphyrins linked by a 1,2-phenylene bridge, differences in the interporphyrin connectivity set the lateral shift between macrocycles, reducing electronic coupling in AB m-β and resulting in a localized exciton. Pump-probe experiments show that the exciton dynamics is faster by almost an order of magnitude in the strongly coupled AA m-m dimer, and two-dimensional electronic spectroscopy (2DES) identifies a vibronic coherence that is absent in AB m-β . Theoretical studies indicate how the interchromophore interactions in these structures, and their system-bath couplings, influence excitonic delocalization and vibronic coherence-enabled rapid exciton transport dynamics. Real-time path integral calculations reproduce the exciton transfer kinetics observed experimentally and find that the linking-modulated exciton delocalization strongly enhances the contribution of vibronic coherences to the exciton transfer mechanism, and that this coherence accelerates the exciton transfer dynamics. These benchmark molecular design, 2DES, and theoretical studies provide a foundation for directed explorations of nonclassical effects on exciton dynamics in multiporphyrin assemblies.

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

Exciton-vibrational dynamics induces efficient self-trapping in a substituted nanoring

Cycloparaphenylenes, being the smallest segments of carbon nanotubes, have emerged as prototypes of the simplest carbon nanohoops. Their unique structure–dynamics–optical properties relationships have motivated a wide variety of synthesis of new related nanohoop species. Studies of how chemical changes, introduced in these new materials, lead to systems with new structural, dynamics and optical properties, expand their functionalities for optoelectronics applications. Herein, we study the effect that conjugation extension of a cycloparaphenylene through the introduction of a satellite tetraphenyl substitution has on its structural and dynamical properties. Our non-adiabatic excited state molecular dynamics simulations suggest that this substitution accelerates the electronic relaxation from the high-energy band to the lowest excited state. This is partially due to efficient conjugation achieved between specific phenyl units as introduced by the tetraphenyl substitution. We observe a particular exciton redistribution during relaxation, in which the tetraphenyl substitution plays a significant role. As a result, an efficient inter-band energy transfer takes place. Besides, the observed phonon-exciton interplay induces a significant exciton self-trapping. Finally, our results encourage and guide the future studies of new phenyl substitutions in carbon nanorings with desired optoelectronic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Physics through the 1990s: Atomic, molecular and optical physics

The volume presents a program of research initiatives in atomic, molecular, and optical physics. The current state of atomic, molecular, and optical physics in the US is examined with respect to demographics, education patterns, applications, and the US economy. Recommendations are made for each field, with discussions of their histories and the relevance of the research to government agencies. The section on atomic physics includes atomic theory, structure, and dynamics; accelerator-based atomic physics; and large facilities. The section on molecular physics includes spectroscopy, scattering theory and experiment, and the dynamics of chemical reactions. The section on optical physics discusses lasers, laser spectroscopy, and quantum optics and coherence. A section elucidates interfaces between the three fields and astrophysics, condensed matter physics, surface science, plasma physics, atmospheric physics, and nuclear physics. Another section shows applications of the three fields in ultra-precise measurements, fusion, national security, materials, medicine, and other topics.

Source record↗
DOE OSTI2026

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗
DOE OSTI2025

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗
DOE OSTI2022

Studying Performance Portability of LAMMPS across Diverse GPU-based Platforms

The molecular dynamics simulation software, LAMMPS, utilizes the Kokkos acceleration library to port computation to a diverse set of architectures including those based on GPU accelerators. In addition to Kokkos, LAMMPS contains a vast code base that leverages the CUDA application programming interface using library functions such as cuFFT, CUDA’s fast-fourier transform (FFT) library, and, more recently, also support for AMD’s Heterogeneous Interface for Portability (HIP) that is rapidly growing. While preparing LAMMPS tests for the AMD GPU-based test system precursors to Frontier, we investigated several strategies for accelerating LAMMPS on AMD GPUs, using the AMD Instinct MI100 and MI250X. In this work, we integrated the HIP FFT library, hipFFT, into the particle-particle particle-mesh (PPPM) long-range solver, which allowed the porting of PPPM calculations to the GPUs. Kokkos behavior on the MI100 and MI250X was also investigated through the package kokkos command of LAMMPS, targeting com- munication, memory usage, and particle grid decomposition. The Tersoff, Reax, Lennard-Jones (LJ), EAM, Granular, and PPPM potentials were investigated in this effort, and results from these experiments are provided. The selected potentials were run on Spock (AMD Instinct MI100), Crusher (AMD Instinct MI250X), AFW HPC11 (NVIDIA A100) and Summit (NVIDIA V100), for comparison. Operational roofline models were constructed and analyzed for the Tersoff, Reax, and Lennard-Jones potentials on Crusher and Summit.

Hagerty, Nick↗
DOE OSTI2021

Spectral Denoising for Accelerated Analysis of Correlated Ionic Transport

Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

The heterogeneous nature of mechanically accelerated grain growth

Abstract While grain growth is traditionally viewed as a purely thermally driven process, nanocrystalline metals can undergo grain growth under mechanical loads, even at room temperature. We performed a detailed atomistic study of the heterogeneous nature of mechanically accelerated grain growth in a polycrystalline Pt nanowire. Using molecular dynamics simulations, we compared the grain-growth behavior of individual grains during tensile and shear cyclic loading, for three different equivalent strain levels, and at two temperatures. Pure thermal grain growth with no mechanical loading provided a baseline reference case. On average, grains that were already susceptible to thermal grain growth were stimulated to grow faster with mechanical loading, as expected. However, when analyzed on a grain-by-grain basis, the results were far more complex: grains that grew fastest under one stimuli were less accelerated under other stimuli. Even when the magnitude of loading changed, the relative growth of individual grains was distorted. We interpret this complexity from the perspective of superimposed growth mechanisms.

Chen, Elton Y.↗
DOE OSTI2022

HPC Molecular Simulation Tries Out a New GPU: Experiences on Early AMD Test Systems for the Frontier Supercomputer

Molecular simulation is an important tool for nu- merous efforts in physics, chemistry, and the biological sciences. Simulating molecular dynamics requires extremely rapid cal- culations to enable sufficient sampling of simulated temporal molecular processes. The Hewlett Packard Enterprise (HPE) Cray EX Frontier supercomputer installed at the Oak Ridge Leadership Computing Facility (OLCF) will provide an exascale resource for open science, and will feature graphics processing units (GPUs) from Advanced Micro Devices (AMD). The future LUMI supercomputer in Finland will be based on an HPE Cray EX platform as well. Here we test the ports of several widely used molecular dynamics packages that have each made substantial use of acceleration with NVIDIA GPUs, on Spock, the early Cray pre-Frontier testbed system at the OLCF which employs AMD GPUs. These programs are used extensively in industry for pharmaceutical and materials research, as well as academia, and are also frequently deployed on high-performance computing (HPC) systems, including national leadership HPC resources. We find that in general, performance is competitive and installation is straightforward, even at these early stages in a new GPU ecosystem. Our experiences point to an expanding arena for GPU vendors in HPC for molecular simulation.

Sedova, Ada↗
DOE OSTI2023

The impact of stacking and phonon environment on energy transfer in organic chromophores: computational insights

Energy transfer in organic materials is extensively studied due to many applications in optoelectronics. The electronic and vibrational relaxations within molecular assemblies can be influenced by stacking arrangements or additions of a backbone that unites them. Here, we present the computational study of the photoexcitation dynamics of a perylene diimide monomer, and face-to-face stacked dimer and trimer. By using non-Energy transfer in organic materials is extensively studied due to many applications in optoelectronics. The electronic and vibrational relaxations within molecular assemblies can be influenced by stacking arrangements or additions of a backbone that unites them. In this work, we present the computational study of the photoexcitation dynamics of a perylene diimide monomer, and face-to-face stacked dimer and trimer. By using non-adiabatic excited-state molecular dynamics simulations, we show that the non-radiative relaxation is accelerated with the number of stacked molecules. This effect is explained by differences in the energy splitting between states that impacts their corresponding nonadiabatic couplings. Additionally, our analysis of the vibronic dynamics reveals that the passage through the different conical intersections that participate in the relaxation of the stacked systems, activate a positive feedback mechanism. This effect involves a narrow set of vibrational normal modes that accelerate the process by increasing the efficiency of its vibronic dynamics. In contrast, an addition of a biologically inspired backbone slows down the relaxation rate due to its participation in the vibronic dynamics of the molecular stacking arrangements. Our results suggest the stacking arrangements and common backbones as strategies to modulate the efficiency of electronic and vibrational relaxation of diimide-based systems and other molecular aggregates.

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

Quantum-Based Molecular Dynamics Simulations Using Tensor Cores

Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency but with the caveat that they were designed for tensor contractions (matrix–matrix multiplications) using only low-precision floating-point operations. Despite this perceived limitation, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born–Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born–Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low-precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.

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

AI-powered exploration of molecular vibrations, phonons, and spectroscopy

The vibrational dynamics of molecules and solids play a critical role in defining material properties, particularly their thermal behaviors. However, theoretical calculations of these dynamics are often computationally intensive, while experimental approaches can be technically complex and resource-demanding. Recent advancements in data-driven artificial intelligence (AI) methodologies have substantially enhanced the efficiency of these studies. This review explores the latest progress in AI-driven methods for investigating atomic vibrations, emphasizing their role in accelerating computations and enabling rapid predictions of lattice dynamics, phonon behaviors, molecular dynamics, and vibrational spectra. Key developments are discussed, including advancements in databases, structural representations, machine-learning interatomic potentials, graph neural networks, and other emerging approaches. Compared to traditional techniques, AI methods exhibit transformative potential, dramatically improving the efficiency and scope of research in materials science. The review concludes by highlighting the promising future of AI-driven innovations in the study of atomic vibrations.

Han, Bowen [Oak Ridge National Laboratory (ORNL), ↗
DOE OSTI2021

Ligand modulation of the conformational dynamics of the A 2A adenosine receptor revealed by single-molecule fluorescence

G protein-coupled receptors (GPCRs) are the largest class of transmembrane proteins, making them an important target for therapeutics. Activation of these receptors is modulated by orthosteric ligands, which stabilize one or several states within a complex conformational ensemble. The intra- and inter-state dynamics, however, is not well documented. Here, we used single-molecule fluorescence to measure ligand-modulated conformational dynamics of the adenosine A 2A receptor (A 2A R) on nanosecond to millisecond timescales. Experiments were performed on detergent-purified A 2R in either the ligand-free (apo) state, or when bound to an inverse, partial or full agonist ligand. Single-molecule Förster resonance energy transfer (smFRET) was performed on detergent-solubilized A 2A R to resolve active and inactive states via the separation between transmembrane (TM) helices 4 and 6. The ligand-dependent changes of the smFRET distributions are consistent with conformational selection and with inter-state exchange lifetimes ≥ 3 ms. Local conformational dynamics around residue 229 6.31 on TM6 was measured using fluorescence correlation spectroscopy (FCS), which captures dynamic quenching due to photoinduced electron transfer (PET) between a covalently-attached dye and proximal aromatic residues. Global analysis of PET-FCS data revealed fast (150–350 ns), intermediate (50–60 μs) and slow (200–300 μs) conformational dynamics in A 2A R, with lifetimes and amplitudes modulated by ligands and a G-protein mimetic (mini-G s ). Most notably, the agonist binding and the coupling to mini-G s accelerates and increases the relative contribution of the sub-microsecond phase. Molecular dynamics simulations identified three tyrosine residues (Y112, Y288 7.53 , and Y290 7.55 ) as being responsible for the dynamic quenching observed by PET-FCS and revealed associated helical motions around residue 229 6.31 on TM6. As a result, this study provides a quantitative description of conformational dynamics in A 2A R and supports the idea that ligands bias not only GPCR conformations but also the dynamics within and between distinct conformational states of the receptor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Light by design: emerging frontiers in ultrafast photon sciences and light–matter interactions

Photon sciences and technologies establish the building blocks for myriad scientific and engineering frontiers in life and energy sciences. Because of their overarching functionality, the developmental roadmap and opportunities underpinned by photonics are often semiotically mediated by the delineation of subject areas of application. In this perspective article, we map current and emerging linkages between three intersecting areas of research stewarded by advanced photonics technologies, namely light by design, outlined as (a) quantum and structured photonics, (b) light–matter interactions in accelerators and accelerator-based light sources, and (c) ultrafast sciences and quantum molecular dynamics. In each section, we will concentrate on state-of-the-art achievements and present prospective applications in life sciences, biochemistry, quantum optics and information sciences, and environmental and chemical engineering, all founded on a broad range of photon sources and methodologies. We hope that this interconnected mapping of challenges and opportunities seeds new concepts, theory, and experiments in the advancement of ultrafast photon sciences and light–matter interactions. Additionally, through this mapping, we hope to inspire a critically interdisciplinary approach to the science and applications of light by design.

42 ENGINEERING↗