Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR DYNAMICS METHOD”

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

Two-dimensional phase cartography for high-harmonic spectroscopy

Extremely nonlinear spectroscopy based on high-order-harmonic generation has become a powerful investigation method for attosecond dynamics in gas and solid targets. In particular, the phase of harmonic emission was shown to carry profound insight into atomic and molecular structure and dynamics. However, current techniques offer phase measurements only along specific directions, thus providing partial characterization. Here we report on a new approach combining optical and quantum interferometers measuring along two dimensions the intensity and phase of harmonic emission from aligned molecules in the exact same experimental conditions. This two-dimensional cartography technique measures the phase with no arbitrary offset and no uncertainty on its sign. Measurements along different dimensions can be combined in two ways: either a single mapping or a redundant mapping allowing high-precision phase recovery using a Shack–Hartmann-like algorithm. We demonstrate both methods in a nitrogen test case, which allows disentangling structural and dynamical effects. Two-dimensional phase cartography paves the way to high-resolution high-harmonic spectroscopy for applications such as quantum orbital tomography and attosecond charge migration in molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Simplified Ab Initio Molecular Dynamics-Based Raman Spectral Simulations

We describe a simplified approach to simulating Raman spectra from ab initio molecular dynamics (AIMD) calculations. The protocol relies on on-the-fly calculations of approximate molecular polarizabilities using the well-known sum over orbitals (as opposed to states) method. This approach bypasses the more accurate but computationally expensive approach to calculating molecular polarizabilities along AIMD trajectories, i.e., solving the coupled perturbed Hartree–Fock/Kohn–Sham equations. We demonstrate the advantages and limitations of our method through a few case studies targeting molecular systems of interest to surface- and/or tip-enhanced Raman spectroscopy practitioners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light–matter hybrid states, so-called polaritons, are formed due to quantum light–matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light–matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light–matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule–cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

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

Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful?

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born–Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. In this study, we address this issue by modeling the electronic energy transfer in a donor–acceptor–donor (D–A–D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are “sufficiently good” to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

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

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
NASA NTRS1989

Pressure broadening and line coupling in bending bands of CO2

The pressure broadening and line coupling cross sections in the Fano-Ben Reuven theory of line shapes are calculated for bending bands of CO2 in a bath of He atoms. Molecular collision dynamics are simplified by invoking the infinite order sudden (IOS) approximation for molecular rotational and vibrational angular momentum in a manner similar to but not identical with the method developed by Clary and shown to be accurate for CO2-He. Numerical values are obtained using a pairwise additive interaction potential developed by Clary. Predictions are in good accord with data for various infrared bands and pure rotational Raman spectra. It is found that all pressure broadening and state-to-state cross sections depend on only a few dynamical factors (generalized IOS cross sections) and are therefore closely interrelated. Results are used to assess models developed previously to analyze line shapes in this and similar systems.

Green, Sheldon↗
DOE OSTI2022

Review of the Effects of Polymer Binder Properties on Microstructure and Irreversible Volume Growth of Plastic Bonded Explosives Formulations

The rational design of effective polymeric binders for the formulation of plastic-bonded explosives (PBX) is challenging due to their inherent compositional complexity. The composites comprise irregularly shaped energetic material (EM) powders coated with low weight fractions of polymer via non-equilibrium processes such as slurry coating. Defects can deleteriously affect PBX stability and performance: nano- to micrometer-scale voids can act as loci for hot spots, lowering deflagration and detonation temperatures in unpredictable ways. Furthermore, some nominally desirable polymer properties are at odds with each other: e. g. good flow characteristics are desirable for coating and adhesion, but mechanical stiffness is needed to prevent deformation and cracking of PBX under mechanical stress. Good binder adhesion is critical, but the best means to predict and measure adhesion in PBX is not obvious. Experimental methods of determining binder adhesion on model surfaces may not capture polymer structural configurations relevant to deposition during coating. Molecular dynamics-based computational models have predicted key observables in PBX formulation, suggesting that they may be powerful tools for binder selection. In this review, primarily recent (~2006 and later) literature on polymeric binders for insensitive HE (IHE) is surveyed. We focus on how binder properties influence observable PBX properties as resistance to irreversible volume growth and void formation in PBX formulations mainly (but not exclusively) featuring 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the EM (e. g. PBX 9502, LX-17 and others). Conclusions from these studies yield useful guidelines for choosing HE binder candidates, as well as for the general design of highly-filled polymer composite materials. Finally, studies describing challenges in PBX formulation with the newer and more energetically dense high explosive, LLM-105, will be discussed.

adhesion↗
DOE OSTI2023

Probing pH-Dependent Dehydration Dynamics of Mg and Ca Cations in Aqueous Solutions with Multi-Level Quantum Mechanics/Molecular Dynamics Simulations

The dehydration of aqueous calcium and magnesium cations is the most fundamental process controlling their reactivity in chemical and biological phenomena, such as the formation of ionic solids or passing through ion channels. It holds particular relevance in light of recent advancements in the development of carbon capture techniques that rely on mineralization for long-term carbon storage. Specifically, dehydration of Ca 2+ and Mg 2+ is a key step in proposed carbon capture processes aiming to exploit the relatively high concentration of dissolved carbon dioxide in seawater via the formation of carbonate minerals from solvated Ca 2+ and Mg 2+ cations for sequestration and storage. Nevertheless, atomic-scale understanding of the dehydration of aqueous Ca 2+ and Mg 2+ cations remains limited. Here, we utilize rare event sampling via density functional theory molecular dynamics and embedded wavefunction theory calculations to elucidate the dehydration dynamics of aqueous Ca 2+ and Mg 2+ . Emphasis is placed on the investigation of the effect pH has on the stability of the different coordination environments. Here, our results reveal significant differences in the dehydration dynamics of the two cations and provide insight into how they may be modulated by pH changes.

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

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores. This method makes it easier to interpret the movement of the protein complex from heterogeneity analysis and to compare the structural dynamics observed from CryoEM data with results from other experimental and simulation techniques.

59 BASIC BIOLOGICAL SCIENCES↗
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 (↗
NASA NTRS1986

Electron-Diffraction Analysis of Growth of GaAs

Report describes experiments that used reflection high-energy-electron diffraction (RHEED) to investigate behavior of GaAs surfaces during and after growth by molecular-beam epitaxy (MBE). Experimental results show dynamic RHEED measurements useful both as probes of surface and growth kinetics and as methods for determining and reproducing surface and growth conditions.

Lewis, B. F.↗
DOE OSTI2024

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Ultrafast molecular frame electronic coherences from lab frame scattering anisotropies

Electronic coherences in molecules are ultrafast charge oscillations on the molecular frame (MF) and their direct observation and separation from electronic population dynamics is challenging. Here we present a valence shell lab frame (LF) scattering method suited to probing electronic coherences in isolated systems. MF electronic coherences lead to LF electronic anisotropies observable by ultrafast angle-resolved scattering. Moment analysis of the LF anisotropy completely separates electronic coherences from population dynamics, demonstrated in excited state NH 3 using ultrafast time-energy-angle-resolved photoelectron spectroscopy. Overall, this general approach applies equally to attosecond/femtosecond electronic coherences in isolated systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

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

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Stochastic Constrained Extended System Dynamics for Solving Charge Equilibration Models

In this work, we present a new stochastic extended Lagrangian molecular dynamics solution to charge equilibration that eliminates self-consistent field (SCF) calculations, thus eliminating the computational bottleneck in solving the charge distribution with standard SCF solvers. By formulating both charges and chemical potential as latent variables and introducing a holonomic constraint that satisfies charge conservation, the SC-XLMD method accurately reproduces thermodynamic, dynamic, and structural properties within the framework of ReaxFF for a bulk water system and highly reactive RDX molecules simulated at high temperature. The SC-XLMD method shows excellent computational performance and is available in the publicly available LAMMPS package.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗