Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics

We report nonadiabatic dynamics, which goes beyond the Born–Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.

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

Nuclear quantum effects on the quasiparticle properties of the chloride anion aqueous solution within the GW approximation

Photoelectron spectroscopy experiments in ionic solutions reveal important electronic structure information, in which the interaction between hydrated ions and water solvent can be inferred. Based on many-body perturbation theory with GW approximation, we theoretically compute the quasiparticle electronic structure of chloride anion solution, which is modeled by path-integral abinitio molecular dynamics simulation by taking account the nuclear quantum effects (NQEs). The electronic levels of hydrated anion as well as water are determined and compared with the recent experimental photoelectron spectra. It is found that NQEs improve the agreement between theoretical prediction and experiment because NQEs effectively weaken the hybridization of the between the Cl – anion and water. Furthermore, our study indicates that NQEs plays a small but non-negligible role in predicting the electronic structure of the aqueous solvation of ions of the Hofmeister series.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Biochemical Process Modeling and Simulation (BPMS)

The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design. We work closely with experimental projects to identify problems and iterate with experiments to find and refine solutions. By working with experimentalists, we decide on problems that can be solved with simulation that could otherwise not be solved or would take too long with experiment alone to reach BETO's targets. Over the years, we have produced solutions that have resulted in determining the most likely fatty-acid derivative for passive transport out of bacteria that upgrade biomass, and we have also designed enzyme mutations for enhanced lignin upgrading. Metabolic models have been developed to tune the activity of 2,3 butanediol production for the 2030 target. A computational method to deliver understanding of how complex omics data can be interpreted in the metabolic pathways of organisms used in the Agile Biofoundry. We have found methods to overcome specific barriers and continue to develop those methods. Our reactor studies have guided the design of both the microbes and reactors for aerobic and micro-aerobic production at all scales and have been instrumental in improving the accuracy of techno-economic analysis models. This project is essential in the process of selecting the final processes for 2030 SAF production targets. More specifically, recently, we have: 1) Predicted the strength of the basic structural interactions in commodity plastics to provide guidance for plastics upcycling strategies. 2) Developed computational tool to improve the characterization of lignin-derived compounds 3) Developed new methodologies to enable Machine Learning-based Directed Evolution for protein engineering. 4) Developed Machine Learning methods to predict protein promiscuity and mutations to further improve microbial and enzymatic driven processes and demonstrated the utility of ML approaches to engineering proteins from sparse experimental datasets. 5) Developed new methods to enable high-fidelity simulation of aerobic fermentation at industrial scale and resolving mismatch of time scales through subcycling/operator splitting 7) Identified the difficulty in preventing local high-oxygen conditions in industrial bubble columns, which leads to less-desirable acetoin production, suggesting future research directions in alternative reactor configurations (e.g loop reactors, shallow-channel reactors).

BIOMASS FUELS↗
DOE OSTI2022

Electron and Electron-Nuclear Dynamics of Open Quantum Subsystems

Leveraging an open-subsystem formulation of Density Functional Theory (DFT) we aim at describing periodic and molecular systems alike, including their electronic and nuclear dynamics. Subsystem DFT enables first principles simulations to approach realistic time- and length-scales, and most importantly sheds light on the dynamical behavior of complex systems. Taking subsystem DFT to the time domain allows us to inspect the electron dynamics of condensed-phase systems in real time. In liquids and interfaces, we observe all the relevant regimes proper of non-Markovian open quantum system dynamics, such as electronic energy transfer, and screening. In addition, the ab-initio modeling of system-bath interactions brought us to observe and justify the holographic time-dependent electron density theorem. Contrary to interactions between molecular (finite) systems, when molecules interact with metal or semiconductor surfaces the electron dynamics is strongly non-Markovian with dramatic repercussions to the molecule’s response to external perturbations. Metals and semiconductors typically have large polarizabilities, and even in a regime of low coupling their effect on impinging molecular species is significant–line broadening, peak shift, and intensity borrowing are observed, characterized, and explained in terms of inter-subsystem dynamical interactions and a many-body decomposition of the system’s density-density response function in a way that transcends the canons of Fermi Golden Rule.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Substitution of distal and active site residues reduces product inhibition of E1 from Acidothermus Cellulolyticus

Abstract Cellulases are largely afflicted by inhibition from their reaction products, especially at high-substrate loading, which represents a major challenge for biomass processing. This challenge was overcome for endoglucanase 1 (E1) from Acidothermus cellulolyticus by identifying a large conformational change involving distal residues upon binding cellobiose. Having introduced alanine substitutions at each of these residues, we identified several mutations that reduced cellobiose inhibition of E1, including W212A, W213A, Q247A, W249A and F250A. One of the mutations (W212A) resulted in a 47-fold decrease in binding affinity of cellobiose as well as a 5-fold increase in the kcat. The mutation further increased E1 activity on Avicel and dilute-acid treated corn stover and enhanced its productivity at high-substrate loadings. These findings were corroborated by funnel metadynamics, which showed that the W212A substitution led to reduced affinity for cellobiose in the +1 and +2 binding sites due to rearrangement of key cellobiose-binding residues.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Chemical diversity in angiosperms − monoterpene synthases control complex reactions that provide the precursors for ecologically and commercially important monoterpenoids

SUMMARY Monoterpene synthases (MTSs) catalyze the first committed step in the biosynthesis of monoterpenoids, a class of specialized metabolites with particularly high chemical diversity in angiosperms. In addition to accomplishing a rate enhancement, these enzymes manage the formation and turnover of highly reactive carbocation intermediates formed from a prenyl diphosphate substrate. At each step along the reaction path, a cationic intermediate can be subject to cyclization, migration of a proton, hydride, or alkyl group, or quenching to terminate the sequence. However, enzymatic control of ligand folding, stabilization of specific intermediates, and defined quenching chemistry can maintain the specificity for forming a signature product. This review article will discuss our current understanding of how angiosperm MTSs control the reaction environment. Such knowledge allows inferences about the origin and regulation of chemical diversity, which is pertinent for appreciating the role of monoterpenoids in plant ecology but also for aiding commercial efforts that harness the accumulation of these specialized metabolites for the food, cosmetic, and pharmaceutical industries.

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

Manipulating and Probing Ultrafast Atomic and Molecular Dynamics (Final Report)

This program has focused on the exploration and control of dynamics in atoms, small molecules, clusters, and micro- or nano-structures following ultrafast excitation and/or driven by strong oscillating fields. In our principal experiments, gas phase atoms and molecules were exposed to short pulses of electromagnetic radiation with frequencies ranging from THz into the XUV spectral regions, and with durations from picoseconds to attoseconds. The pulses coherently induce and/or probe dynamics in the electronic, vibrational, and rotational degrees of freedom of the targets. The overarching goal of the various projects has been to exploit strong-field processes and/or implement novel ultrafast techniques for manipulating and probing coherent electronic and nuclear motion. Typically, the relevant system evolution occurs over attosecond to picosecond timescales. In doing so, we have gained a more complete picture of correlated multi-particle dynamics in atoms, molecules, and other complex systems and expanded the toolbox of useful techniques for manipulating and probing those dynamics. Grant funds were primarily used to: (i) pay salaries and fringe benefits for the PI and postdocs and stipends, tuition, and health care costs for graduate students; (ii) cover travel expenses for the PI, postdocs, and students to attend and participate in professional meetings, conferences, and workshops; and (iii) purchase research materials, supplies, and equipment required for experiments. The research resulted in the publication of 24 peer reviewed articles (including 1 in Nature Physics, 2 in Nature Communications, 7 in Physical Review Letters, and 12 in Physical Review A), and 4 PhD dissertations. The specific results from different funding periods are summarized separately in the main report.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Two-dimensional ring polymer molecular dynamics determination of the MnO + + H 2 /D 2 reaction rates on a Δ-machine learned potential energy surface

In this work, we investigate the impact of nuclear quantum effects in the kinetics of the MnO + + H 2 reaction, a prototypical system for gas-phase H 2 activation by transition metal oxide ions. The DFT based potential energy surfaces (PESs) for the lowest-lying quintet and septet spin states reported in our previous work are improved by 2953 newly calculated CCSD(T)/AVDZ points using a delta-machine learning (Δ-ML) method. To examine nuclear quantum effects, the rate coefficients are computed using ring-polymer molecular dynamics (RPMD). Due to the mechanistic complexity of the reaction, two reaction coordinates are necessary to map out the free-energy surface and an extended RPMD rate theory is developed. The calculated RPMD rate coefficients on the new PES are in better agreement with experimental data. The calculated kinetic isotope effects (KIEs) range from 1.6 to 1.8, also in good agreement with the experimental values. Additionally, the converged RPMD rate coefficients are 2.3 to 3.1 times higher than their classical counterparts across the studied temperature range, suggesting the presence of moderate nuclear quantum effects in this reaction.

Liu, Yang [University of New Mexico, Albuquerque, ↗
DOE OSTI2019

Phonon-Mediated Ultrafast Hole Transfer from Photoexcited CdSe Quantum Dots to Black Dye

We have applied non-adiabatic molecular dynamics (NAMD) and Fewest Switching surface hopping (FSSH) method to simulate the energy dissipation and charge transfer processes in the CdSe quantum dot (QD) functionalized by the Black dye. Our calculations reveal ultrafast hole transfer from the photoexcited QD to the dye, despite a significant energy splitting between QD’s and dye’s occupied states. The mechanism of such ultrafast QD-to-dye hole transfer is rationalized by strong non-adiabatic couplings between the delocalized surface orbitals of the QD and the high-frequency vibrations of isothiocyanate ligands of the Black dye. Here, this mechanism is also confirmed by a very rapid lose of coherency of the excitation due to strong electronic couplings with vibrational modes of the dye.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

On the Dynamics of the Carbon–Bromine Bond Dissociation in the 1-Bromo-2-Methylnaphthalene Radical Anion

This paper studies the mechanism of electrochemically induced carbon–bromine dissociation in 1-Br-2-methylnaphalene in the reduction regime. In particular, the bond dissociation of the relevant radical anion is disassembled at a molecular level, exploiting quantum mechanical calculations including steady-state, equilibrium and dissociation dynamics via dynamic reaction coordinate (DRC) calculations. DRC is a molecular-dynamic-based calculation relying on an ab initio potential surface. This is to achieve a detailed picture of the dissociation process in an elementary molecular detail. From a thermodynamic point of view, all the reaction paths examined are energetically feasible. The obtained results suggest that the carbon halogen bond dissociates following the first electron uptake follow a stepwise mechanism. Indeed, the formation of the bromide anion and an organic radical occurs. The latter reacts to form a binaphthalene intrinsically chiral dimer. This paper is respectfully dedicated to Professors Anny Jutand and Christian Amatore for their outstanding contribution in the field of electrochemical catalysis and electrosynthesis.

1-bromo-2-methylnaphthalene↗
DOE OSTI2024

Exciting DeePMD: Learning excited-state energies, forces, and non-adiabatic couplings

We extend the DeePMD neural network architecture to predict electronic structure properties necessary to perform non-adiabatic dynamics simulations. While learning the excited state energies and forces follows a straightforward extension of the DeePMD approach for ground-state energies and forces, how to learn the map between the non-adiabatic coupling vectors (NACV) and the local chemical environment descriptors of DeePMD is less trivial. Most implementations of machine-learning-based non-adiabatic dynamics inherently approximate the NACVs, with an underlying assumption that the energy-difference-scaled NACVs are conservative fields. We overcome this approximation, implementing the method recently introduced by Richardson [J. Chem. Phys. 158, 011102 (2023)], which learns the symmetric dyad of the energy-difference-scaled NACV. Furthermore, the efficiency and accuracy of our neural network architecture are demonstrated through the example of the methaniminium cation CH 2 NH 2 + .

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

Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics

Practical methods for simulating photo-induced dynamics of molecules often rely on a trajectory-description of the nuclear dynamics, incorporating approximate electron-nuclear feedback. Through approximating electron–nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. Here, a set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.

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

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

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

Dissecting subcycle interference in photoelectron holography

Multipath holographic interference in strong-field quantum tunnel ionization is key to revealing subangstrom attosecond dynamics for molecular movies. This critical subcycle motion is often obscured by longer timescale effects such as ring-shaped patterns that appear in above-threshold ionization (ATI). In the present work, we overcome this problem by combining two techniques in theory and experimental analysis: unit-cell averaging and time-filtering data and simulations. Together these suppress ATI rings and enable an unprecedented highly detailed quantitative match between strong-field ionization experiments in argon and the Coulomb-quantum orbit strong-field approximation (CQSFA) theory. Velocity map images reveal fine modulations on the holographic spiderlike interference fringes that form near the polarization axis. CQSFA theory traces this to the interference of three types of electron pathways. The level of agreement between experiment and theory allows sensitive determination of quantum phase differences and symmetries, providing an important tool for quantitative dynamical imaging in quantum systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Importance of nuclear quantum effects on the hydration of chloride ion

Path-integral ab initio molecular dynamics (PI-AIMD) calculations have been employed to probe the nature of chloride ion solvation in aqueous solution. Nuclear quantum effects (NQEs) are shown to weaken hydrogen bonding between the chloride anion and the solvation shell of water molecules. As a consequence, the disruptive effect of the anion on the solvent water structure is significantly reduced compared to what is found in the absence of NQEs. Additionally, the chloride hydration structure obtained from PI-AIMD agrees well with information extracted from neutron scattering data. In particular, the observed satellite peak in the hydrogen-chloride-hydrogen triple angular distribution serves as a clear signature of NQEs. The present results suggest that NQEs are likely to play a crucial role in determining the structure of saline solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantifying the Long‐Range Coupling of Electronic Properties in Proteins with ab initio Molecular Dynamics**

Abstract The delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.

Yang, Zhongyue↗