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

Molecular tail chemistry controls thermal transport in fullerene films

Herein we report on the thermal conductivities of alkyl- and indene-group functionalized fullerene derivative thin films as measured via time domain thermoreflectance and steady state thermoreflectance. The thermal conductivities vary from 0.064 ± 0.007 W m – 1 K – 1 for [6,6]-phenyl C 61 -butyric acid methyl ester (PCBM) to 0.15 ± 0.017 W m – 1 K – 1 for indene- C 60 bisadduct at room temperature and do not exhibit significant temperature dependence from 300 to 375 K. In comparison to the thermal conductivity of PCBM, increasing the length of the alkyl chain, as in the case of [6,6]-phenyl C 61 -butyric acid butyl ester, and [6,6]-phenyl C 61 -butyric acid octyl ester leads to higher thermal conductivities. Likewise, increasing the number of alkyl chains attached to the fullerenes as in the case of bisadduct PCBM leads to a higher thermal conductivity compared to that of PCBM. We present atomistic insights into the role of chemical functionalization on the overall heat transfer in these fullerene derivatives by conducting molecular dynamics simulations and lattice dynamics calculations. The thermal conductivities predicted via our atomistic simulations qualitatively agree with the experimental trends for our fullerene derivatives. Lattice dynamics calculations reveal that one of the main factors dictating the ultralow thermal conductivities in fullerene derivatives is the large reduction in modal diffusivities in the molecular crystals as calculated from the Allen-Feldman model, thus providing an explanation for their largely reduced thermal conductivities as compared to that of C 60 crystals. The low diffusivities result from high degrees of localization of Einstein-like vibrations in fullerene derivatives due to the molecular side chains, providing the ability to dial-in the properties of these low thermal conductivity solids via molecular engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Impacts of vacancy-induced polarization and distortion on diffusion in Solid Electrolyte Li 3 OCl

We report that Lithium-rich oxychloride antiperovskites are promising solid electrolytes for enabling next-generation batteries. Here, we report a comprehensive study varying Li + concentrations in Li 3 OCl using ab initio molecular dynamics simulations. The simulations accurately capture the complex interactions between Li + vacancies (V$^'_{Li}$), the dominant mobile species in Li 3 OCl. The V$^'_{Li}$ polarize and distort the host lattice, inducing additional non-vacancy-mediated diffusion mechanisms and correlated diffusion events that reduce the activation energy barrier at concentrations as low as 1.5% V$^'_{Li}$. Our analyses of discretized diffusion events in both space and time illustrate the critical interplay between correlated dynamics, polarization and local distortion in promoting ionic conductivity in Li 3 OCl.

25 ENERGY STORAGE↗
DOE OSTI2023

Fermionic quantum processing with programmable neutral atom arrays

Simulating the properties of many-body fermionic systems is an outstanding computational challenge relevant to material science, quantum chemistry, and particle physics.-5.4pc]Please note that the spelling of the following author names in the manuscript differs from the spelling provided in the article metadata: D. González-Cuadra, D. Bluvstein, M. Kalinowski, R. Kaubruegger, N. Maskara, P. Naldesi, T. V. Zache, A. M. Kaufman, M. D. Lukin, H. Pichler, B. Vermersch, Jun Ye, and P. Zoller. The spelling provided in the manuscript has been retained; please confirm. Although qubit-based quantum computers can potentially tackle this problem more efficiently than classical devices, encoding nonlocal fermionic statistics introduces an overhead in the required resources, limiting their applicability on near-term architectures. In this work, we present a fermionic quantum processor, where fermionic models are locally encoded in a fermionic register and simulated in a hardware-efficient manner using fermionic gates. We consider in particular fermionic atoms in programmable tweezer arrays and develop different protocols to implement nonlocal gates, guaranteeing Fermi statistics at the hardware level. We use this gate set, together with Rydberg-mediated interaction gates, to find efficient circuit decompositions for digital and variational quantum simulation algorithms, illustrated here for molecular energy estimation. Finally, we consider a combined fermion-qubit architecture, where both the motional and internal degrees of freedom of the atoms are harnessed to efficiently implement quantum phase estimation as well as to simulate lattice gauge theory dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2026

Interplay between Ultrafast Electronic and Librational Dynamics in Liquid Nitrobenzene Probed with Two-Color Four-Wave Mixing

Here, we present an experimental and theoretical study of the interplay between ultrafast electron dynamics and librational dynamics in liquid nitrobenzene. A femtosecond ultraviolet pulse and two femtosecond near-infrared pulses interact with nitrobenzene molecules, generating a four-wave mixing nonlinear signal measured in the Optical Kerr Effect geometry. The signal is measured to be nonzero only at negative time delays, corresponding to the near-infrared pulses arriving before the ultraviolet pulse. We perform time-dependent Quantum Master Equation calculations with classical libration to simulate the experiment. The simulations support the conclusion that the near-infrared pulses launch librational motion while creating electronic coherences resulting in a libration-modulated electronic nonlinear response. The analysis of the phase-matched four-wave mixing signals suggests a nonparametric process leaving the molecules in an excited electronic state, providing new insight into ultrafast nonlinear optical interactions in liquids and advancing toward probing ultrafast electronic coherences in complex molecular liquids.

Shivaram, Niranjan [Lawrence Berkeley National Lab↗
DOE OSTI2022

Chemistry-mediated Ostwald ripening in carbon-rich C/O systems at extreme conditions

Abstract There is significant interest in establishing a capability for tailored synthesis of next-generation carbon-based nanomaterials due to their broad range of applications and high degree of tunability. High pressure (e.g., shockwave-driven) synthesis holds promise as an effective discovery method, but experimental challenges preclude elucidating the processes governing nanocarbon production from carbon-rich precursors that could otherwise guide efforts through the prohibitively expansive design space. Here we report findings from large scale atomistically-resolved simulations of carbon condensation from C/O mixtures subjected to extreme pressures and temperatures, made possible by machine-learned reactive interatomic potentials. We find that liquid nanocarbon formation follows classical growth kinetics driven by Ostwald ripening (i.e., growth of large clusters at the expense of shrinking small ones) and obeys dynamical scaling in a process mediated by carbon chemistry in the surrounding reactive fluid. The results provide direct insight into carbon condensation in a representative system and pave the way for its exploration in higher complexity organic materials. They also suggest that simulations using machine-learned interatomic potentials could eventually be employed as in-silico design tools for new nanomaterials.

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

Coherent band-edge oscillations and dynamic longitudinal-optical phonon mode splitting as evidence for polarons in perovskites

The coherence of collective modes, such as phonons and polarons, and their modulation of electronic states is long sought in complex systems, which is a crosscutting issue in photovoltaics and quantum electronics. In photovoltaic cells and lasers based on metal halide perovskites, the presence of polarons, i.e., photocarriers dressed by the macroscopic motion of charged lattice, assisted by terahertz (THz) longitudinal-optical (LO) phonons, has been intensely studied yet is still debated. This may be key for explaining the remarkable properties of the perovskite materials, e.g., defect tolerance, long charge lifetimes, and diffusion lengths. In this work we use the intense single-cycle THz pulse with peak electric field up to E THz = 1000 kV/cm to drive coherent polaronic band-edge oscillations at room temperature in CH 3 NH 3 PbI 3 (MAPbI 3 ). We reveal the oscillatory behavior is dominated by a specific quantized lattice vibration mode at $ω$ LO ~ 4 THz, which is both dipole and momentum forbidden. THz-driven coherent polaron dynamics exhibits distinguishing features: room temperature coherent oscillations at ω LO longer than 1 ps in both single crystals and thin films, mode-selective modulation of different band-edge states assisted by electron-phonon interaction, and dynamic mode splitting at low temperature due to entropy and anharmonicity of organic cations. Our results demonstrate intense THz-driven coherent band-edge modulation is a powerful probe of electron-lattice coupling phenomena and polaronic quantum control in perovskites.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum Squeezing and Sensing with Pseudo-Anti-Parity-Time Symmetry

The emergence of parity-time (PT) symmetry has greatly enriched our study of symmetry-enabled non-Hermitian physics, but the realization of quantum PT symmetry faces an intrinsic issue of unavoidable symmetry-breaking Langevin noises. Here we construct a quantum pseudo-anti-PT (pseudo-APT) symmetry in a two-mode bosonic system without involving Langevin noises. We show that the spontaneous pseudo-APT symmetry breaking leads to an exceptional point, across which there is a transition between different types of quantum squeezing dynamics; i.e., the squeezing factor increases exponentially (oscillates periodically) with time in the pseudo-APT-symmetric (broken) region. Such dramatic changes of squeezing factors and quantum dynamics near the exceptional point are utilized for ultraprecision quantum sensing. These exotic quantum phenomena and sensing applications can be experimentally observed in two physical systems: spontaneous wave mixing nonlinear optics and atomic Bose-Einstein condensates. Furthermore, our Letter offers a physical platform for investigating exciting APT symmetry physics in the quantum realm, paving the way for exploring fundamental quantum non-Hermitian effects and their quantum technological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Time-resolved site-selective imaging of predissociation and charge transfer dynamics: the CH 3 I B-band

The predissociation dynamics of the 6s (B 2 E) Rydberg state of gas-phase CH 3 I were investigated by time-resolved Coulomb-explosion imaging using extreme ultraviolet (XUV) free-electron laser pulses. Inner-shell ionization at the iodine 4d edge was utilized to provide a site-specific probe of the ensuing dynamics. The combination of a velocity-map imaging (VMI) spectrometer coupled with the pixel imaging mass spectrometry (PImMS) camera permitted three-dimensional ionic fragment momenta to be recorded simultaneously for a wide range of iodine charge states. In accord with previous studies, initial excitation at 201.2 nm results in internal conversion and subsequent dissociation on the lower-lying A-state surface on a picosecond time scale. Examination of the time-dependent yield of low kinetic energy iodine fragments yields mechanistic insights into the predissociation and subsequent charge transfer following multiple ionization of the iodine products. The effect of charge transfer was observed through differing delay-dependencies of the various iodine charge states, from which critical internuclear distances for charge transfer could be inferred and compared to a classical over-the-barrier model. Time-dependent photofragment angular anisotropy parameters were extracted from the central slice of the Newton sphere, without Abel inversion, and highlight the effect of rotation of the parent molecule before dissociation, as observed in previous works. Our results demonstrate the ability to perform three-dimensional ion imaging at high event rates and showcase the potential benefits of this approach, particularly in relation to further time-resolved studies at free-electron laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

High-fidelity, hyper-accurate, and evolved mutants rewire atomic-level communication in CRISPR-Cas9

The high-fidelity (HF1), hyper-accurate (Hypa), and evolved (Evo) variants of the CRISPR-associated protein 9 (Cas9) endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in recognition subdomain 3 (Rec3) mediate specificity in these variants are poorly understood. Here, solution nuclear magnetic resonance and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic His-Asn-His (HNH) domain. The variants remodel the communication sourcing from the Rec3 α helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Molecular-level insights into structure and dynamics in ionic liquids and polymer gel electrolytes

We report designing new electrolytes requires a better understanding of the correlation between their transport properties and their molecular structure. In this work, we present a detailed study of ionic liquids and polymer gel electrolytes probing their structure and dynamics by nuclear magnetic resonance (NMR) spectroscopy. In particular, ammonium- and phosphonium-based ionic liquids (ILs) are combined with different LiTFSI concentrations, and then with different ratios of poly(methylmethacrylate) polymer. The temperature dependence of self-diffusion coefficients of mobile species, D Li+ , D TFSI- and D P4441+ or D N4441+ , measured by pulsed field gradient (PFG) NMR spectroscopy obeys the Arrhenius equation. Diffusivity of [P4441][TFSI] ILs is found to be greater than those of the [N4441][TFSI] ILs in both liquid and gel electrolytes. Solid state NMR experiments including 13 C and 19 F MAS, 13 C{ 19 F}, 13 C{ 1 H} CPMAS probe the local structure and molecular-level interactions between ions and polymer in gel electrolytes, particularly for samples with higher PMMA content (≥25 wt%). Finally, the fast field cycling relaxometry has been used to unveil the rotational and translational dynamics of P4441 + or N4441 + and TFSI - by measuring 1 H and 19 F R 1 relaxation rate profiles at different temperatures. A comprehensive NMR analysis including relaxation studies at low magnetic field provides decisive new insights regarding the formation of ionic clusters and the interaction of ions with the polymer chain in the case of the gel electrolytes.

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

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

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

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

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

Helfrich-Hurault elastic instabilities driven by geometrical frustration

Here, the Helfrich-Hurault (HH) elastic instability is a well-known mechanism behind patterns that form as a result of strain upon liquid crystal systems with periodic ground states. In the HH model, layered structures undulate and buckle in response to local, geometric incompatibilities in order to maintain the preferred layer spacing. Classic HH systems include cholesteric liquid crystals under electromagnetic field distortions and smectic liquid crystals under mechanical strains, where both materials are confined between rigid substrates. However, richer phenomena are observed when undulation instabilities occur in the presence of deformable interfaces and variable boundary conditions. Understanding how the HH instability is affected by deformable surfaces is imperative for applying the instability to a broader range of materials. In this review, the HH mechanism is reexamined and special focus is given to how the boundary conditions influence the response of lamellar systems to geometrical frustration. Lamellar liquid crystals confined within a spherical shell geometry are used as the model system. Made possible by the relatively recent advances in microfluidics within the past 15 years, liquid crystal shells are composed entirely of fluid interfaces and have boundary conditions that can be dynamically controlled at will. Past and recent work that exemplifies how topological constraints, molecular anchoring conditions, and boundary curvature can trigger the HH mechanism in liquid crystals with periodic ground states is examined. The review ends by identifying similar phenomena across a wide variety of materials, both biological and synthetic. The fact that the HH mechanism is a generic and often overlooked response of periodic materials to geometrical frustration is highlighted.

42 ENGINEERING↗
DOE OSTI2024

Dynamic population balance in molecular-level simulations of hypersonic flows

This report summarizes the work towards developing stochastic weighted particle methods (SWPM) for future application in hypersonic flows. Extensive changes to Sandia’s direct simulation Monte Carlo (DSMC) solver, SPARTA (Stochastic Particle Real Time Analyzer), were made to enable the necessary particle splitting and reduction capabilities for SWPM. The results from one-dimensional Couette and Fourier flows suggest that SWPM can reproduce the correct transport for a large range of Knudsen numbers with adequate accuracy. The associated velocity and temperature profiles are in good agreement with DSMC. An issue with particle placement during particle number reduction, is identified, to which, a simple but effective solution based on minimizing the center of mass error is proposed. High Mach wheel flows are simulated using the SWPM and DSMC methods. SWPM is capable of providing nearly an order of magnitude increase in efficiency over DSMC while retaining high accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Differentiating Three-Dimensional Molecular Structures Using Laser-Induced Coulomb Explosion Imaging

Coulomb explosion imaging (CEI) with x-ray free electron lasers has recently been shown to be a powerful method for obtaining detailed structural information of gas-phase planar ring molecules [R. Boll et al., X-ray multiphoton-induced Coulomb explosion images complex single molecules, Nat. Phys. 18, 423 (2022).]. In this Letter, we investigate the potential of CEI driven by a tabletop laser and extend this approach to differentiating three-dimensional structures. We study the static CEI patterns of planar and nonplanar organic molecules that resemble the structures of typical products formed in ring-opening reactions. Here, our results reveal that each molecule exhibits a well-localized and distinctive pattern in three-dimensional fragment-ion momentum space. We find that these patterns yield direct information about the molecular structures and can be qualitatively reproduced using a classical Coulomb explosion simulation. Our findings suggest that laser-induced CEI can serve as a robust method for differentiating molecular structures of organic ring and chain molecules. As such, it holds great promise as a method for following ultrafast structural changes, e.g., during ring-opening reactions, by tracking the motion of individual atoms in pump-probe experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron–nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional “virtual” nuclear wave packets whose interference with the “real” (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

74 ATOMIC AND MOLECULAR PHYSICS↗