SEARCH · Engineering Papers
Results for “quantum molecular dynamics”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations
We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.
Charge carrier nonadiabatic dynamics in non-metal doped graphitic carbon nitride
Graphitic carbon nitride (GCN) has attracted significant attention due to its excellent performance in photocatalytic applications. Non-metal doping of GCN has been widely used to improve the efficiency of the material as a photocatalyst. Using a combination of time-domain density functional theory with nonadiabatic molecular dynamics, we study the charge carrier dynamics in oxygen and boron doped GCN systems. The reported simulations provide a detailed time-domain mechanistic description of the charge separation and recombination processes that are of fundamental importance while evaluating the photovoltaic and photocatalytic performance of the material. The appearance of smaller energy gaps due to the presence of dopant states improves the visible light absorption range of the doped systems. At the same time, the nonradiative lifetimes are shortened in the doped systems as compared to the pristine GCN. In the case of boron doped at a carbon (B–C–GCN), the charge recombination time is very long as compared to the other two doped systems owing to the smaller electron–phonon coupling strength between the valence band maximum and the trap state. The results suggest B–C–GCN as the most suitable candidate among three doped systems studied in this work for applications in photocatalysis. Finally, this work sheds light into the influence of dopants on quantum dynamics processes that govern GCN performance and, thus, guides toward building high-performance devices in photocatalysis.
Amplifying Magnetic Field Effects on Upconversion Emission via Molecular Qubit-Driven Triplet–Triplet Annihilation
Triplet-triplet annihilation (TTA) enables photon upconversion by combining two lower-energy triplet excitons to produce a higher-energy singlet exciton. This mechanism enhances light-harvesting efficiency for solar energy conversion and enables the use of lower-energy photons in bioimaging and photoredox catalysis applications. The magnetic modulation of such high-energy excitons presents an exciting opportunity to develop molecular quantum information technologies. While the spin dynamics of triplet exciton pairs are sensitive to external magnetic fields, the magnetic field effects (MFEs) associated with these pairs are generally limited by spin statistics to at most 10% at low fields (<1 Tesla), making them challenging to apply in technological advancements. In contrast, MFEs on spin-correlated radical pairs (SCRPs) can be significantly greater, surpassing those on triplet pairs. By using SCRPs-based molecular qubits as triplet sensitizers in the sensitized TTA scheme, we can magnetically modulate TTA and consequently, the delayed fluorescence of annihilators. In our current system, we have achieved more than 70% magnetic modulation of delayed fluorescence, effectively harnessing and even amplifying magnetic modulation within SCRPs to influence high-energy excitons. In conclusion, this work opens new opportunities for advancing spin-controlled chemical reactions and molecular quantum information technologies.
Quantum noise and the threshold of hearing
It is argued that the sensitivity of the ear reaches a limit imposed by the uncertainty principle. This is possible only if the receptor cell holds the detector elements in a special nonequilibrium state which has the same noise characteristics as a ground (T = 0 K) state. To accomplish this 'active cooling' the molecular dynamics of the system must maintain quantum mechanical coherence over the time scale of the measurement.
Phase transitions and dimensional cross-over in layered confined solids
The nature of solid phases and cross-over of order–disorder phase transitions from two-dimensional (2D) layers to three-dimensional (3D) bulk in confined atomic systems remain largely unexplained. To this end, we consider noble gases and aluminum confined between graphene sheets at different pressures and temperatures. Using crystal structure search methods and molecular dynamics based on machine-learned potentials with quantum-mechanical accuracy, we identify structures of multilayer confined solids that deviate from simple close packing. Upon heating, we find that confined 2D monolayers melt according to the two-step continuous Kosterlitz–Thouless–Halperin–Nelson–Young theory. However, multilayer solids transition continuously into an intermediate layered-hexatic phase before melting discontinuously into an isotropic liquid. This intermediate phase persists at least up to 12 layers studied here. This change can be qualitatively understood based on the cross-over from 2D topological defects toward 3D ones during melting as the number of layers increases.
Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states
Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. Furthermore, this work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.
Step-by-step state-selective tracking of fragmentation dynamics of water dications by momentum imaging
Abstract The double photoionization of a molecule by one photon ejects two electrons and typically creates an unstable dication. Observing the subsequent fragmentation products in coincidence can reveal a surprisingly detailed picture of the dynamics. Determining the time evolution and quantum mechanical states involved leads to deeper understanding of molecular dynamics. Here in a combined experimental and theoretical study, we unambiguously separate the sequential breakup via D + + OD + intermediates, from other processes leading to the same D + + D + + O final products of double ionization of water by a single photon. Moreover, we experimentally identify, separate, and follow step by step, two pathways involving the b 1 Σ + and a 1 Δ electronic states of the intermediate OD + ion. Our classical trajectory calculations on the relevant potential energy surfaces reproduce well the measured data and, combined with the experiment, enable the determination of the internal energy and angular momentum distribution of the OD + intermediate.
Unveiling the spatial distribution of molecular coherences at conical intersections by covariance X-ray diffraction signals
Significance Monitoring the real-space motion of the electrons directly involved in ultrafast molecular dynamics is a challenging goal, hindered by the strong contributions of the remaining inactive electrons. Molecular quantum coherences can provide a direct window into elementary electronic dynamics, but they could not be separately imaged so far. Here, we demonstrate a time- and frequency-resolved X-ray diffraction technique, which can employ existing stochastic X-ray pulses to separately image the coherence contributions. The technique can selectively follow the relevant electrons involved in the ultrafast dynamics and identify the confined molecular regions covered by their motion. It enables readily feasible measurements and can be straightforwardly extended to full three-dimensional imaging via X-ray and electron diffraction.
Environmental Contributions to Proton Sharing in Protein Low-Barrier Hydrogen Bonds
Hydrogen bonds (H-bonds) are central to biomolecular structure and dynamics. Although H-bonds are typically characterized by well-defined proton positions, proton delocalization has been proposed to play a role in facilitating enzyme catalysis and allostery in some systems. Experimentally locating protons is difficult, hampering the study of proton mobility in H-bonds. We used neutron crystallography, atomic resolution X-ray bond length analysis, and large quantum mechanics/molecular mechanics-Born–Oppenheimer molecular dynamics (QM/MM-BOMD) simulations to comprehensively characterize the shared proton/deuteron in a Glu–Asp low-barrier hydrogen bond (LBHB) in the bacterial protein YajL that is a conventional H-bond in the homologous disease-associated human protein DJ-1. X-ray bond length analysis of protiated and perdeuterated DJ-1 and YajL shows no significant effect of deuteron substitution on these carboxylic acid-carboxylate H-bonds but does reveal an effect at the active site glutamic acid near a cysteine thiolate. Residues in an H-bonded network that might favor LBHB formation in YajL were interrogated by the mutation of homologous residues in DJ-1. A distal DJ-1 substitution increases proton delocalization in the Glu–Asp H-bond, demonstrating that mutations within extended H-bond networks can modulate proton transfer barriers in carboxylic acid-carboxylate H-bonds. In addition, proton mobility in the H-bond is correlated with dimer-spanning motions in the QM/MM-BOMD simulations of YajL and DJ-1. Our results show that proton delocalization can be tuned using combined bioinformatic, structural, and computational information, opening the possibility of using engineered proton delocalization as a probe of H-bonding environments and as a tool to test hypotheses about LBHB function.
Atomic scale etching of diamond: insights from molecular dynamics simulations
Diamond is a promising material for multiple applications in quantum information processing and sensing as well as applications in microelectronics. However, diamond devices can be limited by surface defects that compromise charge stability and spin coherence, among others. Improved strategies in plasma etching of diamond could play an important role in minimizing or eliminating these defects. In this work, we explore plasma-assisted atomic scale etching of diamond using argon ions (Ar + ), hydrogen ions (H + ) and hydrogen atoms (H). We employ classical molecular dynamics (MD) simulations and test several interatomic potentials based on the Reactive Empirical Bond Order (REBO) form with comparisons to a variety of published experimental results. We performed MD simulations of low-energy hydrogen ($\leqslant$50 eV) and argon ( $\leqslant$200 eV) ion bombardment of diamond surfaces. Ar + bombardment can be used to locally smooth initially rough diamond surfaces via the formation of an amorphous C layer, the thickness of which increases with argon ion energy. Subsequent exposure with hydrogen ions (or fast neutrals) will selectively etch this amorphous C layer, leaving the underlying diamond layer mostly intact if the H energy is maintained below about 10 eV. The simulations suggest that combining Ar + smoothing with selective, near threshold energy H removal of amorphous C can be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices.
Control of Dipolar Dynamics by Geometrical Programming
We propose and theoretically analyze methods for quantum many-body control through geometric reshaping of molecular tweezer arrays. Dynamic rearrangement during entanglement is readily available due to the extended coherence times of molecular rotational qubits. We show how motional dephasing can be suppressed and enhanced spin squeezing can be achieved in an actively rearranged short-range XY model. We also analyze in detail a specific static geometry that significantly suppresses decoherence. These general methods as applied to programmable quantum systems offer robust control modalities that are well suited to molecules.
Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition
Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.
Quantum dynamics simulations of the 2D spectroscopy for exciton polaritons
We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton–polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light–matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton–phonon coupling among polariton states due to the strong light–matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity.
DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations
DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. Here, the flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.
T 1 Anisotropy Elucidates Spin Relaxation Mechanisms in an S = 1 Cr(IV) Optically Addressable Molecular Qubit
Paramagnetic molecules offer unique advantages for quantum information science owing to their spatial compactness, synthetic tunability, room-temperature quantum coherence, and potential for optical state initialization and readout. However, current optically addressable molecular qubits are hampered by rapid spin–lattice relaxation ( T 1 ) even at sub-liquid nitrogen temperatures. Here, we use temperature- and orientation-dependent pulsed electron paramagnetic resonance (EPR) to elucidate the negative sign of the ground state zero-field splitting (ZFS) and assign T 1 anisotropy to specific types of motion in an optically addressable S = 1 Cr(o-tolyl) 4 molecular qubit. The anisotropy displays a distinct sin 2 (2θ) functional form that is not observed in S = 1/2 Cu(acac) 2 or other Cu(II)/V(IV) microwave addressable molecular qubits. The Cr(o-tolyl) 4 T 1 anisotropy is ascribed to couplings between electron spins and rotational motion in low-energy acoustic or pseudoacoustic phonons. Our findings suggest that rotational degrees of freedom should be suppressed to maximize the coherence temperature of optically addressable qubits.
Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics
In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024
Hydrogen isotope population near dislocations in zirconium from molecular dynamics
Not Available