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

Light-induced electron spin qubit coherences in the purple bacteria reaction center protein

Photosynthetic reaction center proteins (RCs) provide ideal model systems for studying quantum entanglement between multiple spins, a quantum mechanical phenomenon wherein the properties of the entangled particles become inherently correlated. Following light-generated sequential electron transfer, RCs generate spin-correlated radical pairs (SCRPs), also referred to as entangled spin qubit (radical) pairs (SQPs). Understanding and controlling coherence mechanisms in SCRP/SQPs is important for realizing practical uses of electron spin qubits in quantum sensing applications. The bacterial RC (bRC) provides an experimental system for exploring quantum effects in the SCRP P 865 + Q A − , where P 865 , a special pair of bacteriochlorophylls, is the primary donor, and Q A is the primary quinone acceptor. In this study, we focus on understanding how local molecular environments and isotopic substitution, particularly deuteration, influence spin coherence times (T M ). Using high-frequency electron paramagnetic resonance (EPR) spectroscopy, we observed that the local environment surrounding P 865 and Q A plays a significant role in determining T M . Our findings show that while deuteration led to a modest increase in T M , particularly at low temperatures, but the effect was substantially smaller than predicted by classical nuclear spin diffusion alone. This result is in contrast to our previous study of the photosystem I (PSI) RC, where no increase in T M was observed upon deuteration. Theoretical modeling identified several methyl groups at key distances from the spin centers of both bRC and PSI, and methyl group tunneling at low temperatures has been previously suggested as a mechanism for enhanced spin decoherence. Additionally, our study revealed a strong dependence of spin coherence on the orientation of the external magnetic field, highlighting the influence of the protein microenvironment on spin dynamics. In conclusion, these results offer new insights for optimizing coherence times in quantum system design for quantum information science and sensing applications.

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

Breaking the curse of dimensionality: Solving configurational integrals for crystalline solids by tensor networks

Accurately evaluating configurational integrals for dense solids remains a central and difficult challenge in the statistical mechanics of condensed systems. Here, we present a tensor network approach that reformulates the high-dimensional configurational integral for identical-particle crystals into a sequence of computationally efficient summations. We represent the integrand as a high-dimensional tensor and apply tensor-train (TT) decomposition together with a custom TT-cross interpolation. This approach circumvents the need to explicitly construct the full tensor. We introduce tailored rank-1 and rank-2 schemes optimized for sharply peaked Boltzmann probability densities, typical for identical-particle crystals. When applied to the calculation of internal energy and pressure-temperature curves for crystalline Cu and Ar at high (GPa) pressures, as well as the alpha-to-beta phase transition diagram of Sn, our method accurately reproduces molecular dynamics simulation results using tight-binding, machine learning, hierarchical interacting particle–neural network, and modified embedded atom method potentials,all within seconds of computation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Sensing ultrashort electronic coherent beating at conical intersections by single-electron pulses

In a theoretical study, we present an ultrafast technique for probing time-dependent molecular charge densities. An ultrafast optical pump first brings the molecule into an electronic nonstationary state. This is followed by coherent inelastic scattering of a broadband single-electron probe pulse with a variable delay T, which is detected spectrally. The technique is applied to reveal phase-sensitive background-free coherent electron beating in the conical intersection passage in uracil and reveals the otherwise elusive coherent beating of strongly coupled electrons and nuclei.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Adaptive pruning-based optimization of parameterized quantum circuits

Abstract Variational hybrid quantum–classical algorithms are powerful tools to maximize the use of noisy intermediate-scale quantum devices. While past studies have developed powerful and expressive ansatze, their near-term applications have been limited by the difficulty of optimizing in the vast parameter space. In this work, we propose a heuristic optimization strategy for such ansatze used in variational quantum algorithms, which we call ‘parameter-efficient circuit training (PECT)’. Instead of optimizing all of the ansatz parameters at once, PECT launches a sequence of variational algorithms, in which each iteration of the algorithm activates and optimizes a subset of the total parameter set. To update the parameter subset between iterations, we adapt the Dynamic Sparse Reparameterization scheme which was originally proposed for training deep convolutional neural networks. We demonstrate PECT for the Variational Quantum Eigensolver, in which we benchmark unitary coupled-cluster ansatze including UCCSD and k -UpCCGSD, as well as the Low-Depth Circuit Ansatz (LDCA), to estimate ground state energies of molecular systems. We additionally use a layerwise variant of PECT to optimize a hardware-efficient circuit for the Sycamore processor to estimate the ground state energy densities of the one-dimensional Fermi-Hubbard model. From our numerical data, we find that PECT can enable optimizations of certain ansatze that were previously difficult to converge and more generally can improve the performance of variational algorithms by reducing the optimization runtime and/or the depth of circuits that encode the solution candidate(s).

Physics↗
DOE OSTI2025

Neuromorphic heat transport effects in a molecular junction

Understanding energy transport at the nanoscale is an open and fundamental challenge in the molecular sciences with direct implications for the design of new electronics, computing devices, and materials. While nanoscale energy transport under steady-state conditions has been studied extensively, there is much less known about energy transport under time-dependent driving forces, particularly in the far-from-equilibrium regime. In this work, we use nonequilibrium molecular dynamics simulations and stochastic thermodynamics to investigate energy transport in a well-studied nanoscale system—a molecular junction—subjected to a time-periodic temperature gradient. The primary observation is that molecular junctions can exhibit heat transport hysteresis, a phenomenon in which the heat flux through a system depends not only on the instantaneous value of a time-dependent temperature bias but also on the temporal history of that bias. The presented findings illustrate that molecular junctions can exhibit the specific memory effect—heat transport hysteresis—that is essential for the design of thermal neuromorphic computers. This work elucidates a potential pathway toward the realization of such devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ultrafast infrared nano-imaging of far-from-equilibrium carrier and vibrational dynamics

Ultrafast infrared nano-imaging has demonstrated access to ultrafast carrier dynamics on the nanoscale in semiconductor, correlated-electron, or polaritonic materials. However, mostly limited to short-lived transient states, the contrast obtained has remained insufficient to probe important long-lived excitations, which arise from many-body interactions induced by strong perturbation among carriers, lattice phonons, or molecular vibrations. Here, we demonstrate ultrafast infrared nano-imaging based on excitation modulation and sideband detection to characterize electron and vibration dynamics with nano- to micro-second lifetimes. As an exemplary application to quantum materials, in phase-resolved ultrafast nano-imaging of the photoinduced insulator-to-metal transition in vanadium dioxide, a distinct transient nano-domain behavior is quantified. In another application to lead halide perovskites, transient vibrational nano-FTIR spatially resolves the excited-state polaron-cation coupling underlying the photovoltaic response. These examples show how heterodyne pump-probe nano-spectroscopy with low-repetition excitation extends ultrafast infrared nano-imaging to probe elementary processes in quantum and molecular materials in space and time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scaling of silicon nanoparticle growth in low temperature flowing plasmas

Low temperature plasmas are an emerging method to synthesize high quality nanoparticles (NPs). An established and successful technique to produce NPs is using a capacitively coupled plasma (CCP) in cylindrical geometry. Although a robust synthesis technique, optimizing or specifying NP properties using CCPs, is challenging. In this paper, results from a computational investigation for the growth of silicon NPs in flowing inductively coupled plasmas (ICPs) using Ar/SiH 4 gas mixtures of up to a few Torr are discussed. ICPs produce more locally constrained and quiescent plasma potentials. These positive plasma potentials produce an electrostatic trap for negatively charged NPs, which can significantly extend the residence time of NPs in the plasma, which in turn provides a controllable period for particle growth. The computational platforms used in this study consist of a two-dimensional plasma hydrodynamics model, a three-dimensional nanoparticle growth and trajectory tracking model, and a molecular dynamics simulation for deriving reactive sticking coefficients of silane radicals on Si NPs. Trends for the nanoparticle growth as a function of SiH 4 inlet fraction, gas residence time, energy deposition per particle, pressure, and reactor diameter are discussed. The general path for particle synthesis is the trapping of small NPs in the positive electrostatic potential, followed by entrainment in the gas flow upon reaching a critical particle size. Optimizing or controlling NP synthesis then depends on the spatial distribution of plasma potential, the density of growth species, and the relative time that particles spend in the electrostatic trap and flowing through higher densities of growth species upon leaving the trap.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Multiscale Microtubule Dynamics in Active Nematics

In microtubule-based active nematics, motor-driven extensile motion of microtubule bundles powers chaotic large-scale dynamics. We quantify the interfilament sliding motion both in isolated bundles and in a dense active nematic. The extension speed of an isolated microtubule pair is comparable to the molecular motor stepping speed. In contrast, the net extension in dense 2D active nematics is significantly slower; the interfilament sliding speeds are widely distributed about the average and the filaments exhibit both contractile and extensile relative motion. These measurements highlight the challenge of connecting the extension rate of isolated bundles to the multimotor and multifilament interactions present in a dense 2D active nematic. Furthermore, they also provide quantitative data that is essential for building multiscale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fragmentation Dynamics of Fluorene Explored Using Ultrafast XUV-Vis Pump-Probe Spectroscopy

We report on the use of extreme ultraviolet (XUV, 30.3 nm) radiation from the Free-electron LASer in Hamburg (FLASH) and visible (Vis, 405 nm) photons from an optical laser to investigate the relaxation and fragmentation dynamics of fluorene ions. The ultrashort laser pulses allow to resolve the molecular processes occurring on the femtosecond timescales. Fluorene is a prototypical small polycyclic aromatic hydrocarbon (PAH). Through their infrared emission signature, PAHs have been shown to be ubiquitous in the universe, and they are assumed to play an important role in the chemistry of the interstellar medium. Our experiments track the ionization and dissociative ionization products of fluorene through time-of-flight mass spectrometry and velocity-map imaging. Multiple processes involved in the formation of each of the fragment ions are disentangled through analysis of the ion images. The relaxation lifetimes of the excited fluorene monocation and dication obtained through the fragment formation channels are reported to be in the range of a few tens of femtoseconds to a few picoseconds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Sequence-selective dynamic covalent assembly of information-bearing oligomers

Relatively robust dynamic covalent interactions have been employed extensively to mediate molecular self-assembly reactions; however, these assembly processes often do not converge to a thermodynamic equilibrium, instead yielding mixtures of kinetically-trapped species. Here, we report a dynamic covalent self-assembly process that mitigates kinetic trapping such that multiple unique oligomers bearing covalently coreactive pendant groups are able to undergo simultaneous, sequence-selective hybridization with their complementary strands to afford biomimetic, in-registry molecular ladders with covalent rungs. Analogous to the thermal cycling commonly employed for nucleic acid melting and annealing, this is achieved by raising and lowering the concentration of a multi-role reagent to effect quantitative dissociation and subsequently catalyze covalent bond rearrangement, affording selective assembly of the oligomeric sequences. The hybridization specificity afforded by this process further enabled information encoded in oligomers to be retrieved through selective hybridization with complementary, mass-labeled sequences.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Forces that control self-organization of chemically-propelled Janus tori

Control of the individual and collective behavior of self-propelled synthetic micro-objects has immediate application for nanotechnology, robotics, and precision medicine. Despite significant progress in the synthesis and characterization of self-propelled Janus (two-faced) particles, predictive understanding of their behavior remains challenging, especially if the particles have anisotropic forms. Here, by using molecular simulation, we describe the interactions of chemically-propelled microtori near a wall. The results show that a torus hovers at a certain distance from the wall due to a combination of gravity and hydrodynamic flows generated by the chemical activity. Moreover, electrostatic dipolar interactions between the torus and the wall result in a spontaneous tilt and horizontal translation, in a qualitative agreement with experiment. Simulations of the dynamics of two and four tori near a wall provide evidence for the formation of stable self-propelled bound states. Our results illustrate that self-organization at the microscale occurs due to a combination of multiple factors, including hydrodynamic, chemical, electrostatic and steric interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗
DOE OSTI2025

A spline-based method to obtain spatially dependent viscosity in confined flows

Coupling chemical physics to continuum theories is a critical step to understanding multi-scale phenomena. This paper will connect non-equilibrium molecular dynamics simulations to a continuum-based Navier-Stokes equation that has relaxed the assumption of spatial uniformity in viscosity. Using a form for viscosity based on spline interpolation, viscosity as a function of position is obtained from the least squares fit of the velocity profile measured from molecular simulations of flow in a nanochannel. Viscosity can vary widely, particularly near the channel boundaries, indicating that uniform viscosity is no longer appropriate. Variations of the viscosity near the channel surfaces imply that considering solution and surface chemistry could be necessary to rigorously understand molecular-scale flows in nanochannels.

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

Viscosity of the magnetized strongly coupled one-component plasma

The viscosity tensor of the magnetized one-component plasma, consisting of five independent shear viscosity coefficients, a bulk viscosity coefficient, and a cross coefficient, is computed using equilibrium molecular dynamics simulations and the Green-Kubo relations. A broad range of Coulomb coupling and magnetization strength conditions are studied. Magnetization is found to strongly influence the shear viscosity coefficients when the gyrofrequency exceeds the Coulomb collision frequency. Three regimes are identified as the Coulomb coupling strength and magnetization strength are varied. Here, the Green-Kubo relations are used to separate kinetic and potential energy contributions to each viscosity coefficient, showing how each contribution depends upon the magnetization strength. The shear viscosity coefficient associated with the component of the pressure tensor parallel to the magnetic field, and the two coefficients associated with the component perpendicular to the magnetic field, are all found to merge to a common value at strong Coulomb coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Influence of solution phase environmental heterogeneity and fluctuations on vibronic spectra: Perylene diimide molecular chromophore complexes in solution

Ensembles of ab initio parameterized Frenkel-exciton model Hamiltonians for different perylene diimide dimer systems are used, together with various dissipative quantum dynamics approaches, to study the influence of the solvation environment and fluctuations in chromophore relative orientation and packing on the vibronic spectra of two different dimer systems: a π-stacked dimer in aqueous solution in which the relative chromophore geometry is strongly confined by a phosphate bridge and a side-by-side dimer in dichloromethane involving a more flexible alkyne bridge that allows quasi-free rotation of the chromophores relative to one another. These entirely first-principles calculations are found to accurately reproduce the main features of the experimental absorption spectra, providing a detailed mechanistic understanding of how the structural fluctuations and environmental interactions influence the vibronic dynamics and spectroscopy of solutions of these multi-chromophore complexes.

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

Geometric Phase in Anisotropic Kepler Problem: Perspective for Realization in Rydberg Atoms

We predict a gyroscopic effect that can be demonstrated with Rydberg atoms following the dynamics of a Kepler Hamiltonian with an additional uniaxial anisotropy induced by optical ponderomotive force. This effect is analogous to the rotation of the Foucault pendulum in response to the Earth’s rotation. Furthermore, we argue that in Rydberg states with a large principal quantum number a similar geometric angle can be generated by mechanical rotations of an atomic-optical setup on timescales between 1 μ⁢s and 1 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗