Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of pore connectivity and tortuosity on the dynamics of fluids confined in sub-nanometer pores

Dynamical behavior of fluids under nano-pore confinement is studied extensively as it has important implications for several industrial as well as geological processes. Pore network in many porous materials exhibits a varied degree of inter connections. The extent of this pore connectivity may affect the structural and dynamical behavior of the confined fluid. However, studies of fluid confinement addressing these effects systematically are lacking. In this paper we report molecular dynamics simulation studies addressing the effects of pore connectivity on the dynamics of two representative fluids – CO 2 and ethane in silicalite by systematically varying the degree of pore connectivity through selectively blocking some pore space with immobile methane molecules. By selectively turning off the pore spaces in the shape of straight, or tortuous zigzag channels, we also probe the effects of pore tortuosity. In general, pore connectivity is found to facilitate both the translational as well as rotational dynamics of both fluids, while the intermolecular modes of vibration in both fluids remain largely unaffected. The effects of providing connections between a set of straight or zigzag channel-like pores are however more nuanced. Pore tortuosity facilitates the rotational motion, but suppresses the translational motion of CO 2 , while its effects on the rotational and translational motion of ethane are less pronounced. The intermolecular vibrational modes of both fluids shift to higher energies with an increase in the number of tortuous pores. The results reported here provide a detailed molecular level understanding of the effects of pore connectivity on the dynamics of fluids and thus have implications for applications like fluid separation.

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

Phase-matching analysis in high-order harmonic generation with nonzero orbital angular momentum Laguerre-Gaussian beams

Through high-order harmonic generation driven by intense ultrashort vortex infrared or midinfrared lasers, a nonzero orbital angular momentum can be imprinted onto extreme ultraviolet (XUV) or soft-x-ray (SXR) light pulses. Here we simulate the generation of vortex XUV harmonics in the gas medium as well as their propagation in vacuum till reaching the far field. In this work, we find that the intensity and phase of generated high harmonics are very sensitive to the position of gas jet with respect to the laser focus. The topological charge of the $\textit{q}$th harmonic is found to be $\textit{q}$ times that of the driving Laguerre-Gaussian beam. Each harmonic in the far field appears as a single ring in the transverse plane with an invariant diameter which is scalable with the fundamental topological charge only when the gas jet is placed after the laser focus. The underlying phase-matching mechanism is analyzed by examining the spatial map of the coherence length and by calculating the evolution of harmonic emission in the medium. We anticipate this work to stimulate interest in generating intense vortex XUV or SXR attosecond pulses for probing dynamics of molecules where special molecular features are difficult to be detected with linear or circular XUV or SXR pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

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

Accelerating Ab Initio Simulation via Nested Monte Carlo and Machine Learned Reference Potentials

As a corollary of the rapid advances in computing, ab initio simulation is playing an increasingly important role in modeling materials at the atomic scale. Two strategies are possible, ab initio Monte Carlo (AIMC) and molecular dynamics (AIMD) simulation. The former benefits from exact sampling from the correct thermodynamic distribution, while the latter is typically more efficient with its collective all-atom coordinate updates. In this study, using a relatively simple test model comprised of helium and argon, we show that AIMC can be brought up to, and even above, the performance levels of AIMD via a hybrid nested sampling/machine learning (ML) strategy. Here, ML provides an accurate classical reference potential (up to three-body explicit interactions) that can pilot long collective Monte Carlo moves that are accepted or rejected in toto à la nested Monte Carlo (NMC); this is in contrast to the single move nature of a naive implementation. Our proposed method only requires a small up front expense from evaluating the ab initio energies and forces of (100) random configurations for training. Importantly, our method does not totally rely on the trained, assuredly imperfect, interaction. We show that high performance and exact sampling at the desired level of theory can be realized even when the trained interaction has appreciable differences from the ab initio potential. Remarkably, at the highest levels of performance realized via our approach, a pair of statistically uncorrelated atomic configurations can be generated with (1) ab initio calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Accurate Modeling of Bromide and Iodide Hydration with Data-Driven Many-Body Potentials

Ion–water interactions play a central role in determining the properties of aqueous systems in a wide range of environments. However, a quantitative understanding of how the hydration properties of ions evolve from small aqueous clusters to bulk solutions and interfaces remains elusive. Here, we introduce the second generation of data-driven many-body energy (MBnrg) potential energy functions (PEFs) representing bromide–water and iodide–water interactions. The MB-nrg PEFs use permutationally invariant polynomials to reproduce two-body and three-body energies calculated at the coupled cluster level of theory, and implicitly represent all higher-body energies using classical many-body polarization. Further, a systematic analysis of the hydration structure of small Br - (H 2 O) n and I - (H 2 O) n clusters demonstrates that the MBnrg PEFs predict interaction energies in quantitative agreement with “gold standard” coupled cluster reference values. Importantly, when used in molecular dynamics simulations carried out in the isothermal-isobaric ensemble for single bromide and iodide ions in liquid water, the MB-nrg PEFs predict extended X-ray absorption fine structure (EXAFS) spectra that accurately reproduce the experimental spectra, which thus allows for characterizing the hydration structure of the two ions with high level of confidence.

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

Enhancing Graphene Plasmonic Device Performance via its Dielectric Environment

Graphene plasmons provide a compelling avenue toward chip-scale dynamic tuning of infrared light. Dynamic tunability emerges through controlled alterations in the optical properties of the system defining graphene’s plasmonic dispersion. Typically, electrostatic induced alterations of the carrier concentration in graphene working in conjunction with mobility have been considered the primary factors dictating plasmonic tunability. We find here that the surrounding dielectric environment also plays a primary role, dictating not just the energy of the graphene plasmon but so too the magnitude of its tuning and spectral width. To arrive at this conclusion, poles in the imaginary component of the reflection coefficient are used to efficiently survey the effect of the surrounding dielectric on the tuning of the graphene plasmon. By investigating several common polar materials, optical phonons (i.e., the Reststrahlen band) of the dielectric substrate are shown to appreciably affect not only the plasmon’s spectral location but its tunability, and its resonance shape as well. In particular, tunability is maximized when the resonances are spectrally distant from the Reststrahlen band, whereas sharp resonances (i.e., high-Q) are achievable at the band’s edge. Overall, these observations both underscore the necessity of viewing the dielectric environment in aggregate when considering the plasmonic response derived from two-dimensional materials and provide heuristics to design dynamically tunable graphene-based infrared devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Active liquid crystals powered by force-sensing DNA-motor clusters

Significance Single-molecule techniques have elucidated how isolated molecular motors generate piconewton forces with unprecedented detail. However, in diverse biological and synthetic settings, force-generating proteins collectively power nonequilibrium dynamics, including continuous large-scale rearrangements and persistent fluid flows. Characterizing motor-generated forces in these dense and dynamical environments remains a challenge. We assembled a reversible DNA-based force-sensing probe that, by an optical readout, reveals the molecular arrangements and the force loads experienced by kinesin motors. These probes provide insight into motor-generated forces that collectively power the unique dynamics of microtubule-based active nematics, a noteworthy example of an internally driven active matter system. DNA-based force probes can be extended to study forces and stresses in various synthetic systems as well as diverse cellular environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2024

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. Here, we use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

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

Second-harmonic generation in zinc blende crystals under combined action of femtosecond optical and strong terahertz fields

The influence of a short intense (with an electric field strength up to 250 kV cm{sup −1}) terahertz (THz) pulse on the generation of second harmonic (SH) of Ti : sapphire laser radiation in crystals of zinc blende type (InAs and GaAs), characterised by nonzero bulk quadratic susceptibility, is investigated. It is experimentally shown for InAs(100) that, in the case of s-polarised first and second harmonics, an application of s-polarised THz field changes significantly the SH signal. The THz field-induced azimuthal dependence of the SH signal is in good agreement with the results of theoretical calculation within a phenomenological approach. The dependence of the SH signal on the delay time between the optical and THz pulses is investigated. This dependence for the GaAs crystal repeats the envelope of the THz pulse intensity, whereas in the case of InAs crystal there is a significant discrepancy, caused by the nonlinear dynamics of strong THz field in InAs. (nonlinear optical phenomena)

36 MATERIALS SCIENCE↗
DOE OSTI2021

Drive dependence of the Hall angle for a sliding Wigner crystal in a magnetic field

We numerically examine the depinning and sliding dynamics of a Wigner crystal in the presence of quenched disorder and a magnetic field. In the disorder-free limit, the Wigner crystal Hall angle is independent of crystal velocity, but when disorder is present, we find that the Hall angle starts near zero at the depinning threshold and increases linearly with increasing drive before reaching a saturation close to the disorder free value at the highest drives. The drive dependence is the result of a side jump effect produced when the charges move over pinning sites. The magnitude of the side jump is reduced at the higher velocities. Here, the drive dependent Hall angle is robust for a wide range of disorder parameters and should be a generic feature of classical charges driven in the presence of quenched disorder and a magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Recombination of Autodissociated Water Ions in a Nanoscale Pure Water Droplet

The recombination of water ions has diverse scientific and practical implications, ranging from acid-base chemistry and biological systems to planetary environments and applications in fuel cell and carbon conversion technologies. While spatial confinement affects the physicochemical properties of water dynamics, its impact on the recombination process has rarely been studied. In this work, we investigate the dynamics of water, the water ion distribution, and the ion recombination process in water droplets as a function of droplet size through molecular dynamics simulations and adaptive quantum mechanical/molecular mechanical calculations. We compare the dynamics of recombination in water droplet sizes ranging from 100 to 18 000 waters, both in their interiors and on their surfaces. We found that the self-diffusion of water dramatically decreases in droplets with a diameter below 2.2 nm. Using a classical RexPoN force-field, we found that the ions in 1000 H2O's spend almost 50% of the time on the surface and 0.5 nm beneath it with a slight preference for OH- ion to reside longer on the surface. We estimate that, on average, recombination in these drops occurs at 400 ps in 1000 H2O's and 1 ns in 3000 H2O's. We also found that recombination is not limited by the local structure of the surface or the size of the droplet but can be influenced by the geometry of the water wire connecting the ions as they approach each other, which can often prevent recombination. Our results provide insights to the reaction microenvironments presented by nanoscopic water droplets.

Kwon, Soonho↗
DOE OSTI2020

Understanding the Interplay between the Nonvalence and Valence States of the Uracil Anion upon Monohydration

In this work we present an ab initio investigation into the effect of monohydration on the interaction of uracil with low energy electrons. Electron attachment and photodetachment experimental studies have previously shown dramatic changes in uracil upon solvation with even a single water molecule, due to an inversion of the character of the ground state of the anion. Here we explore the interplay between the nonvalence and valence states of the uracil anion, as a function of geometry and site of solvation. Our model provides unambiguous interpretation of previous photoelectron studies, reproducing the binding energies and photoelectron images for bare uracil and a single isomer of the U•(H 2 O) 1 cluster. Furthermore, the results of this study provide insight into how electrons may attach to hydrated nucleobases. These results lay the foundations for further investigations into the effect of microhydration on the electronic structure and electron capture dynamics of nucleobases.

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

Dissociation slowdown by collective optical response under strong coupling conditions

We consider an ensemble of diatomic molecules resonantly coupled to an optical cavity under strong coupling conditions at normal incidence. Photodissociation dynamics is examined via direct numerical integration of the coupled Maxwell–Schrödinger equations with molecular rovibrational degrees of freedom explicitly taken into account. It is shown that the dissociation is significantly affected (slowed down) when the system is driven at its polaritonic frequencies. The observed effect is demonstrated to be of transient nature and has no classical analog. An intuitive explanation of the dissociation slowdown at polaritonic frequencies is proposed.

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

The Birth of Minerals: From Single Step to Multiple Step Mechanisms

Mineral formation from ions in aqueous solutions begins with complex initial stages, where amorphous and liquid-like precursors play pivotal roles before crystalline growth occurs. Both classical and non-classical nucleation and growth theories, introduced in previous chapters, offer explanations, each with their own strengths and limitations, for the complex intermediate phases observed in experimental research. Analytical techniques play a critical role in detecting and characterizing precursor phases, offering valuable insights into nucleation and growth mechanisms across various temporal and spatial scales. In conclusion, molecular dynamics and modelling provide in-depth perspectives on these phases, allowing for a closer examination of their nucleation and growth mechanisms at the molecular level, and revealing the intricate processes that govern their behaviour.

amorphous↗