Engineering Papers⌕ Search

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.

At least 361 records · Page 20

Path sampling of recurrent neural networks by incorporating known physics

Recurrent neural networks have seen widespread use in modeling dynamical systems in varied domains such as weather prediction, text prediction and several others. Often one wishes to supplement the experimentally observed dynamics with prior knowledge or intuition about the system. While the recurrent nature of these networks allows them to model arbitrarily long memories in the time series used in training, it makes it harder to impose prior knowledge or intuition through generic constraints. In this work, we present a path sampling approach based on principle of Maximum Caliber that allows us to include generic thermodynamic or kinetic constraints into recurrent neural networks. We show the method here for a widely used type of recurrent neural network known as long short-term memory network in the context of supplementing time series collected from different application domains. These include classical Molecular Dynamics of a protein and Monte Carlo simulations of an open quantum system continuously losing photons to the environment and displaying Rabi oscillations. Our method can be easily generalized to other generative artificial intelligence models and to generic time series in different areas of physical and social sciences, where one wishes to supplement limited data with intuition or theory based corrections.

59 BASIC BIOLOGICAL SCIENCES↗

Fock-Space Schrieffer–Wolff Transformation: Classically-Assisted Rank-Reduced Quantum Phase Estimation Algorithm

We present an extension of many-body downfolding methods to reduce the resources required in the quantum phase estimation (QPE) algorithm. In this paper, we focus on the Schrieffer–Wolff (SW) transformation of the electronic Hamiltonians for molecular systems that provides significant simplifications of quantum circuits for simulations of quantum dynamics. We demonstrate that by employing Fock-space variants of the SW transformation (or rank-reducing similarity transformations (RRST)) one can significantly increase the locality of the qubit-mapped similarity-transformed Hamiltonians. The practical utilization of the SW-RRST formalism is associated with a series of approximations discussed in the manuscript. In particular, amplitudes that define RRST can be evaluated using conventional computers and then encoded on quantum computers. The SW-RRST QPE quantum algorithms can also be viewed as an extension of the standard state-specific coupled-cluster downfolding methods to provide a robust alternative to the traditional QPE algorithms to identify the ground and excited states for systems with various numbers of electrons using the same Fock-space representations of the downfolded Hamiltonian. The RRST formalism serves as a design principle for developing new classes of approximate schemes that reduce the complexity of quantum circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular dynamics discovery of an extraordinary ionic migration mechanism in dislocation-containing TlBr crystals

TlBr can surpass CZT as the leading semiconductor for γ- and X-radiation detection. Unfortunately, the optimum properties of TlBr quickly decay when an operating electrical field is applied. Quantum mechanical studies indicated that if this property degradation comes from the conventional mechanism of ionic migration of vacancies, then an unrealistically high vacancy concentration is required to account for the rapid aging of TlBr seen in experiments. Here, we have applied large scale molecular dynamics simulations to study the effects of dislocations on ionic migration of TlBr crystals under electrical fields. We found that electrical fields can drive the motion of edge dislocations in both slip and climb directions. These combined motions eject enormous vacancies in the dislocation trail. Both dislocation motion and a high vacancy concentration can account for the rapid aging of the TlBr detectors. These findings suggest that strengthening methods to pin dislocations should be explored to increase the lifetimes of TlBr crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Quantum Simulations of Vibrational Strong Coupling via Path Integrals

A quantum simulation of vibrational strong coupling (VSC) in the collective regime via thermostated ring-polymer molecular dynamics (TRPMD) is reported. For a collection of liquid-phase water molecules resonantly coupled to a single lossless cavity mode, the simulation shows that as compared with a fully classical calculation, the inclusion of nuclear and photonic quantum effects does not lead to a change in the Rabi splitting but does broaden polaritonic line widths roughly by a factor of 2. Moreover, under thermal equilibrium, both quantum and classical simulations predict that the static dielectric constant of liquid water is largely unchanged inside vs outside the cavity. Furthermore, this result disagrees with a recent experiment demonstrating that the static dielectric constant of liquid water can be resonantly enhanced under VSC, suggesting either limitations of our approach or perhaps other experimental factors that have not yet been explored.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photonic Interrogation and Control of Nano Processes

My research activities for the summer of 2003 consisted of two projects: One project was concerned with determining a method for predicting the static and dynamic assembly properties of nano-structures using laser tweezers. The other project was to investigate the generation of Laguerre-Gaussian modes using a spatial light modulator incorporated into an optical tweezers system. Concerning the first project, I initially pursued the approach suggested by my NASA colleague Dr. Art Decker. This approach involved mimicking the model of the structure of atomic nucleus for the assembly of 1 to 100 atoms using allowed quadruple transitions induced by orbital angular momentums of a Laguerre- Gaussian (Doughnut) laser mode. After realizing the inaptness of the nuclear model with the nanostructure model as far as the binding forces and transitions were concerned, I focused on using quantum dot modei. This model was not attuned also for the host lattice influences the electronic structure of the quantum dot. Thus one other option that I decided to pursue was the approach of molecular quantum mechanics. In this approach the nanostructure is treated as a large (10-100 nm) molecule constructed from single element or multi-elements. Subsequent to consultation with Dr. Fred Morales, a chemical engineer at NASA GRC, and Dr. David Ball, a computational chemist at Cleveland State University, I acquired a molecular-quantum computation software, Hyperchem 7.0. This software allows simulation of different molecular structures as far as their static and dynamic behaviors are concerned. The time that I spent on this project was about eight weeks. Once this suitable approach was identified, I realized the need to collaborate with a computational quantum chemist to pursue searching for stable nanostructures in the range of 10-100 nm that we can be assembled using laser tweezers. The second project was about generating laser tweezers that possess orbital angular momentum. As shown, we were able to generate laser tweezers modes of different orbital angular momentum using a spatial light modulator incorporated into a laser tweezers system. The motivation for investigating these types of modes stems from being able to spin particles at high speeds and also to orient two particles in separate traps and then join them together. Also, there has been recent intense interest on fundamental physics research on orbital angular momentum of light. The fact that circularly polarized light may have associated with it angular momentum that relates to the spin of individual photons (spin 0 for the plane polarized light, spin +1 for the right-circularly polarized light and spin -1 for the left-circularly polarized light) was first demonstrated by Beth in 1936. Orbital angular momentum is, however, distinct from spin in that the spin angular momentum of light is intrinsically linked to the behavior of the electric field in the light whereas orbital angular momentum is a consequence of inclined wavefronts. In 1992 L. Allen, et al showed that the Laguerre-Gaussian (LG) modes could possess well-defined orbital angular momentum that can exceed 1 planck's constant, i.e. l plancks constant per photon, where l is the azimuthal index of the mode.

Jassemnejad, Baha↗

Elucidating Quintet-State Dynamics in Singlet Fission Oligomers and Polymers with Tetracene Pendants

To unlock the potential of molecular engineering for practical quantum sensing and computing, it is essential to create and control pure magnetic states in molecular systems. Singlet fission (SF) in organic materials offers a promising approach by generating pairs of triplet excited states from photoexcited singlets. Here, in this work, we investigate SF in a polymer with strategically positioned tetracene pendant groups along a polynorbornene backbone and its oligomeric counterparts, facilitating intrapolymer through-space coupling. Using continuous-wave and pulsed time-resolved electron paramagnetic resonance (EPR) spectroscopy, we elucidate the spin dynamics and identify key intermediates, including the quintet state, that emerge during SF. Our findings reveal that exciton translational motion along the pendant groups enhances the dissociation of triplet pairs, with oligomer length playing a critical role in modulating spin state interconversion and exciton transport. Our results provide key insights into the SF mechanism in polymeric materials and highlight the role of oligomer length in modulating spin state interconversion and exciton transport. This work advances our understanding of SF in polymers, paving the way for their application in quantum information science and energy conversion technologies.

Bindra, Jasleen K. [Argonne National Laboratory (A↗

Cavity molecular dynamics simulations of vibrational polariton-enhanced molecular nonlinear absorption

We report recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry– Pérot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity-induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics simulations to numerically investigate the linear and the nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light–matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption that is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes—whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states—under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Phase Fullerene Nanotechnology: 1996

NASA has used exotic materials for spacecraft and experimental aircraft to good effect for many decades. In spite of many advances, transportation to space still costs about $10,000 per pound. Drexler has proposed a hypothetical nanotechnology based on diamond and investigated the properties of such molecular systems. These studies and others suggest enormous potential for aerospace systems. Unfortunately, methods to realize diamonoid nanotechnology are at best highly speculative. Recent computational efforts at NASA Ames Research Center and computation and experiment elsewhere suggest that a nanotechnology of machine phase functionalized fullerenes may be synthetically relatively accessible and of great aerospace interest. Machine phase materials are (hypothetical) materials consisting entirely or in large part of microscopic machines. In a sense, most living matter fits this definition. To begin investigation of fullerene nanotechnology, we used molecular dynamics to study the properties of carbon nanotube based gears and gear/shaft configurations. Experiments on C60 and quantum calculations suggest that benzyne may react with carbon nanotubes to form gear teeth. Han has computationally demonstrated that molecular gears fashioned from (14,0) single-walled carbon nanotubes and benzyne teeth should operate well at 50-100 gigahertz. Results suggest that rotation can be converted to rotating or linear motion, and linear motion may be converted into rotation. Preliminary results suggest that these mechanical systems can be cooled by a helium atmosphere. Furthermore, Deepak has successfully simulated using helical electric fields generated by a laser to power fullerene gears once a positive and negative charge have been added to form a dipole. Even with mechanical motion, cooling, and power; creating a viable nanotechnology requires support structures, computer control, a system architecture, a variety of components, and some approach to manufacture. Additional information is contained within the original extended abstract.

Globus, Al↗

Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature

It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. Here, we find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD)

Here, we present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-MD classical and nuclear quantum dynamical simulations in which the energies and forces of the full system are approximated by the 2-, 3-, and 4-body terms of the MBE expansion and compare the average potential and the vibrational density of states (VDOS) with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the 4-body term converges with near identical behavior to the one from the full simulation. The 3-body makes a substantial contribution (~20%) to the energy, whereas the 4- body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (~2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm -1 ) at the 4-body level, and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a 4-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N 4 ,N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Importance of Nuclear Quantum Effects on Aqueous Electrolyte Transport under Confinement in Ti 3 C 2 MXenes

Protons display a high chemical activity and strongly affect the charge storage capability in confined interlayer spaces of two-dimensional (2D) materials. As such, an accurate representation of proton dynamics under confinement is important for understanding and predicting charge storage dynamics in these materials. While often ignored in atomistic-scale simulations, nuclear quantum effects (NQEs), e.g., tunneling, can be significant under confinement even at room temperature. Here, using the thermostatted ring polymer molecular dynamics implementation of path integral molecular dynamics (PIMD) in conjunction with the ReaxFF force field, density functional tight binding (DFTB), and NequIP neural network potential simulations, we investigate the role of NQEs on proton and water transport in bulk water and aqueous electrolytes under confinement in Ti 3 C 2 MXenes. Although overall NQEs are relatively small, especially in bulk, we find that they can alter both quantitative values and qualitative trends on both proton transport and water self-diffusion under confinement relative to classical MD predictions. Therefore, our results suggest the need for NQEs to be considered to simulate aqueous systems under confinement for both qualitative and quantitative accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane

To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. Finally, the results reported here are in reasonable agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Influence of tungsten doping on nonradiative electron–hole recombination in monolayer MoSe 2 with Se vacancies

Two-dimensional transition metal dichalcogenides (TMDs) are receiving significant attention due to their excellent electronic and optoelectronic properties. The material quality is greatly affected by defects that are inevitably generated during material synthesis. Focusing on chalcogenide vacancies, which constitute the most common defect, we use the state-of-the-art simulation methodology developed in our group to demonstrate that W doping of MoSe 2 with Se vacancies reduces charge carrier losses by two mechanisms. First, W doping makes the formation of double Se vacancies unfavorable, while it is favorable in undoped MoSe 2 . Second, if a Se vacancy is present, the charge carrier lifetimes are extended in the W-doped MoSe 2 . Combining ab initio real-time time-dependent density functional theory with nonadiabatic molecular dynamics, the simulations show that the nonradiative carrier losses in the presence of Se vacancies proceed by sub-10 ps electron trapping and relaxation down the manifold of trap states, followed by a 100 ps recombination of trapped electrons with free holes. The electron–vibrational energy exchange is driven by both in-plane and out-of-plane vibrational motions of the MoSe 2 layer. Additionally, the atomistic studies advance our understanding of the influence of defects on charge carrier properties in TMDs and guide improvements of material quality and development of TMD applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular qubits based on photogenerated spin-correlated radical pairs for quantum sensing

Photogenerated spin-correlated radical pairs (SCRPs) in electron donor–bridge–acceptor (D–B–A) molecules can act as molecular qubits and inherently spin qubit pairs. SCRPs can take singlet and triplet spin states, comprising the quantum superposition state. Their synthetic accessibility and well-defined structures, together with their ability to be prepared in an initially pure, entangled spin state and optical addressability, make them one of the promising avenues for advancing quantum information science. Coherence between two spin states and spin selective electron transfer reactions form the foundation of using SCRPs as qubits for sensing. We can exploit the unique sensitivity of the spin dynamics of SCRPs to external magnetic fields for sensing applications including resolution-enhanced imaging, magnetometers, and magnetic switch. Molecular quantum sensors, if realized, can provide new technological developments beyond what is possible with classical counterparts. While the community of spin chemistry has actively investigated magnetic field effects on chemical reactions via SCRPs for several decades, we have not yet fully exploited the synthetic tunability of molecular systems to our advantage. This review offers an introduction to the photogenerated SCRPs-based molecular qubits for quantum sensing, aiming to lay the foundation for researchers new to the field and provide a basic reference for researchers active in the field. We focus on the basic principles necessary to construct molecular qubits based on SCRPs and the examples in quantum sensing explored to date from the perspective of the experimentalist.

Mani, Tomoyasu (ORCID:0000000241255195)↗