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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantum Simulation of Polarized Light-Induced Electron Transfer with a Trapped-Ion Qutrit System

Electron transfer within and between molecules is crucial in chemistry, biochemistry, and energy science. Here this study describes a quantum simulation method that explores the influence of light polarization on electron transfer between two molecules. By implementing precise and coherent control among the quantum states of trapped atomic ions, we can induce quantum dynamics that mimic the electron-transfer dynamics in molecules. We use three-level systems (qutrits), rather than traditional two-level systems (qubits), to enhance the simulation efficiency and realize high-fidelity simulations of electron-transfer dynamics. We treat the quantum interference between the electron coupling pathways from a donor with two degenerate excited states to an acceptor and analyze the transfer efficiency. We also examine the potential error sources that enter the quantum simulations. The trapped-ion systems have favorable scalings with system size compared to those of classical computers, promising access to richer electron-transfer simulations.

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

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗
DOE OSTI2023

Inherent-state melting and the onset of glassy dynamics in two-dimensional supercooled liquids

Below the onset temperature T o , the equilibrium relaxation time of most glass-forming liquids exhibits glassy dynamics characterized by a super-Arrhenius temperature dependence. In this supercooled regime, the relaxation dynamics also proceeds through localized elastic excitations corresponding to hopping events between inherent states, i.e., potential-energy-minimizing configurations of the liquid. Despite its importance in distinguishing the supercooled regime from the high-temperature regime, the microscopic origin of T o is not yet known. Here, we construct a theory for the onset temperature in two dimensions and find that an inherent-state melting transition, described by the binding–unbinding transition of dipolar elastic excitations, delineates the supercooled regime from the high-temperature regime. The corresponding melting transition temperature is in good agreement with the onset temperature found in various two-dimensional (2D) atomistic models of glass formers and an experimental binary colloidal system confined to a water–air interface. Additionally, we find the predictions for the renormalized elastic moduli to agree with the experimentally observed values for the latter 2D colloidal system. We further discuss the predictions of our theory on the displacement and density correlations at supercooled conditions, which are consistent with observations of the Mermin–Wagner fluctuations in experiments and molecular simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Two-Dimensional Projected-Momentum Covariance Mapping for Coulomb Explosion Imaging

We introduce projected-momentum covariance mapping, an extension of recoil-frame covariance mapping for 2D ion imaging studies. By considering the two-dimensional projection of the ion momenta as recorded by the detector, one opens the door to a complex suite of analysis tools adapted from three-dimensional momentum imaging studies. This includes the use of different frames of reference to unravel the dynamics of fragmentation and the application of fragment momentum constraints to isolate specific fragmentation channels. The technique is demonstrated on data from a two-dimensional ion imaging study of the Coulomb explosion of the cis and trans isomers of 1,2-dichloroethene, following strong-field ionization by an intense near-infrared femtosecond laser pulse. Classical simulations are used to guide the interpretation of projected-momentum covariance maps. The results offer a detailed insight into the distinct Coulomb explosion dynamics for this pair of isomers and lay the groundwork for future time-resolved studies of photoisomerization dynamics in this molecular system.

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

Observation of a highly conductive warm dense state of water with ultrafast pump–probe free-electron-laser measurements

The electrical conductivity of water under extreme temperatures and densities plays a central role in modeling planetary magnetic fields. Experimental data are vital to test theories of high-energy-density water and assess the possible development and presence of extraterrestrial life. These states are also important in biology and chemistry studies when specimens in water are confined and excited using ultrafast optical or free-electron lasers (FELs). Here we utilize femtosecond optical lasers to measure the transient reflection and transmission of ultrathin water sheet samples uniformly heated by a 13.6 nm FEL approaching a highly conducting state at electron temperatures exceeding 20 000 K. The experiment probes the trajectory of water through the high-energy-density phase space and provides insights into changes in the index of refraction, charge carrier densities, and AC electrical conductivity at optical frequencies. At excitation energy densities exceeding 10 MJ/kg, the index of refraction falls to n = 0.7, and the thermally excited free-carrier density reaches ne = 5 × 1027 m−3, which is over an order of magnitude higher than that of the electron carriers produced by direct photoionization. Significant specular reflection is observed owing to critical electron density shielding of electromagnetic waves. The measured optical conductivity reaches 2 × 104 S/m, a value that is one to two orders of magnitude lower than those of simple metals in a liquid state. At electron temperatures below 15 000 K, the experimental results agree well with the theoretical calculations using density-functional theory/molecular-dynamics simulations. With increasing temperature, the electron density increases and the system approaches a Fermi distribution. In this regime, the conductivities agree better with predictions from the Ziman theory of liquid metals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonradiative relaxation dynamics in the [Au 25- n Ag n (SH) 18 ] -1 ( n = 1, 12, 25) thiolate-protected nanoclusters

Evaluation of the electron-nuclear dynamics and relaxation mechanisms of gold and silver nanoclusters and their alloys is important for future photocatalytic, light harvesting, and photoluminescence applications of these systems. In this work, the effect of silver doping on the nonradiative excited state relaxation dynamics of the atomically precise thiolate-protected gold nanocluster [Au 25- n Ag n (SH) 18 ] –1 ( n = 1, 12, 25) is studied theoretically. Time-dependent density functional theory is used to study excited states lying in the energy range 0.0–2.5 eV. The fewest switches surface hopping method with decoherence correction was used to investigate the dynamics of these states. The HOMO–LUMO gap increases significantly upon doping of 12 silver atoms but decreases for the pure silver nanocluster. Doped clusters show a different response for ground state population increase lifetimes and excited state population decay times in comparison to the undoped system. The ground state recovery times of the S 1 –S 6 states in the first excited peak were found to be longer for [Au 13 Ag 12 (SH) 18 ] –1 than the corresponding recovery times of other studied nanoclusters, suggesting that this partially doped nanocluster is best for preserving electrons in an excited state. The decay time constants were in the range of 2.0–20 ps for the six lowest energy excited states. Among the higher excited states, S 7 has the slowest decay time constant although it occurs more quickly than S 1 decay. Overall, these clusters follow common decay time constant trends and relaxation mechanisms due to the similarities in their electronic structures.

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

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics↗
DOE OSTI2023

Scattering solution of an interacting Hamiltonian for the electronic control of molecular spin qubits

We theoretically study how a scattered electron can entangle molecular spin qubits (MSQs). This requires solving the inelastic transport of a single electron through a scattering region described by a tight-binding interacting Hamiltonian. We accomplish this using a Green's-function solution. We can model realistic physical implementations of MSQs by parametrizing the tight-binding Hamiltonian with first-principles descriptions of magnetic anisotropy and exchange interactions. As a result, we find that, for two-MSQ systems with inversion symmetry, projective measurement of the spin degree of freedom of the scattered electron offers probabilistic control of the degree of entanglement between the MSQs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

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

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2024

How Rigid Are Anthranilamide Molecular Electrets?

As important as molecular electrets are for electronic materials and devices, conformational fluctuations strongly impact their macrodipoles and intrinsic properties. Herein, we employ molecular dynamics (MD) simulations with the polarizable charge equilibrium (PQEq) method to investigate the persistence length (L P ) of molecular electrets composed of anthranilamide (Aa) residues. The PQEq-MD dissipates the accepted static notions about Aa macromolecules, and L P represents the shortest Aa rigid segments. The classical model with a single L P value does not describe these oligomers. Introducing multiple L P values for the same macromolecule follows the observed trends and discerns the enhanced rigidity in their middle sections from the reduced stiffness at their terminal regions. Furthermore, L P distinctly depends on solvent polarity. The Aa oligomers maintain extended conformations in nonpolar solvents with L P exceeding 4 nm, while in polar media, increased conformational fluctuations reduce L P to about 2 nm. These characteristics set key guidelines about the utility of Aa conjugates for charge-transfer systems within organic electronics and energy engineering.

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

Water transport in reverse osmosis membranes is governed by pore flow, not a solution-diffusion mechanism

We performed nonequilibrium molecular dynamics (NEMD) simulations and solvent permeation experiments to unravel the mechanism of water transport in reverse osmosis (RO) membranes. The NEMD simulations reveal that water transport is driven by a pressure gradient within the membranes, not by a water concentration gradient, in marked contrast to the classic solution-diffusion model. We further show that water molecules travel as clusters through a network of pores that are transiently connected. Permeation experiments with water and organic solvents using polyamide and cellulose triacetate RO membranes showed that solvent permeance depends on the membrane pore size, kinetic diameter of solvent molecules, and solvent viscosity. This observation is not consistent with the solution-diffusion model, where permeance depends on the solvent solubility. Motivated by these observations, we demonstrate that the solution-friction model, in which transport is driven by a pressure gradient, can describe water and solvent transport in RO membranes.

42 ENGINEERING↗
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↗