SEARCH · Engineering Papers

Results for “quantum memories”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2024

First-principles investigation of near-field energy transfer between localized quantum emitters in solids

We present a predictive and general approach to investigate near-field energy transfer processes between localized defects in semiconductors, which couples first-principles electronic structure calculations and a nonrelativistic quantum electrodynamics description of photons in the weak-coupling regime. The approach is general and can be readily applied to investigate broad classes of defects in solids. We apply our approach to investigate an exemplar point defect in an oxide, the F center in MgO, and we show that the energy transfer from a magnetic source, e.g., a rare-earth impurity, to the vacancy can lead to spin nonconserving long-lived excitations that are dominant processes in the near field, at distances relevant to the design of photonic devices and ultrahigh dense memories. We also define a descriptor for coherent energy transfer to predict geometrical configurations of emitters to enable long-lived excitations, that are useful to design optical memories in semiconductor and insulators. Published by the American Physical Society 2024

Chattaraj, Swarnabha (ORCID:0000000329333581)↗

DOE OSTI2020

Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction

We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N 3 ) operations and O(N 2 ) memory, where N is the number of basis functions. Furthermore, this limits the overall computational scaling to O(N 4 ) operations and O(N 2 ) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Quasiprobabilistic Readout Correction of Midcircuit Measurements for Adaptive Feedback via Measurement Randomized Compiling

Quantum measurements are a fundamental component of quantum computing. However, on present-day quantum computers, measurements can be more error prone than quantum gates and are susceptible to nonunital errors as well as nonlocal correlations due to measurement crosstalk. While readout errors can be mitigated in postprocessing, this is inefficient in the number of qubits due to a combinatorially large number of possible states that need to be characterized. In this work, we show that measurement errors can be tailored into a simple stochastic error model using randomized compiling, enabling the efficient mitigation of readout errors via quasiprobability distributions reconstructed from the measurement of a single preparation state in an exponentially large confusion matrix. We demonstrate the scalability and power of this approach by correcting readout errors without matrix inversion on a large number of different preparation states applied to a register of eight superconducting transmon qubits. Moreover, we show that this method can be extended to midcircuit measurements used for active feedback via quasiprobabilistic error cancellation, and we demonstrate the correction of measurement errors on an ancilla qubit used to detect and actively correct bit-flip errors on an entangled memory qubit. Our approach enables the correction of readout errors on large numbers of qubits and offers a strategy for correcting readout errors in adaptive circuits in which the results of midcircuit measurements are used to perform conditional operations on nonlocal qubits in real time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2005

Simulation of Laser Cooling and Trapping in Engineering Applications

An advanced computer code is undergoing development for numerically simulating laser cooling and trapping of large numbers of atoms. The code is expected to be useful in practical engineering applications and to contribute to understanding of the roles that light, atomic collisions, background pressure, and numbers of particles play in experiments using laser-cooled and -trapped atoms. The code is based on semiclassical theories of the forces exerted on atoms by magnetic and optical fields. Whereas computer codes developed previously for the same purpose account for only a few physical mechanisms, this code incorporates many more physical mechanisms (including atomic collisions, sub-Doppler cooling mechanisms, Stark and Zeeman energy shifts, gravitation, and evanescent-wave phenomena) that affect laser-matter interactions and the cooling of atoms to submillikelvin temperatures. Moreover, whereas the prior codes can simulate the interactions of at most a few atoms with a resonant light field, the number of atoms that can be included in a simulation by the present code is limited only by computer memory. Hence, the present code represents more nearly completely the complex physics involved when using laser-cooled and -trapped atoms in engineering applications. Another advantage that the code incorporates is the possibility to analyze the interaction between cold atoms of different atomic number. Some properties that cold atoms of different atomic species have, like cross sections and the particular excited states they can occupy when interacting with each other and light fields, play important roles not yet completely understood in the new experiments that are under way in laboratories worldwide to form ultracold molecules. Other research efforts use cold atoms as holders of quantum information, and more recent developments in cavity quantum electrodynamics also use ultracold atoms to explore and expand new information-technology ideas. These experiments give a hint on the wide range of applications and technology developments that can be tackled using cold atoms and light fields. From more precise atomic clocks and gravity sensors to the development of quantum computers, there will be a need to completely understand the whole ensemble of physical mechanisms that play a role in the development of such technologies. The code also permits the study of the dynamic and steady-state operations of technologies that use cold atoms. The physical characteristics of lasers and fields can be time-controlled to give a realistic simulation of the processes involved such that the design process can determine the best control features to use. It is expected that with the features incorporated into the code it will become a tool for the useful application of ultracold atoms in engineering applications. Currently, the software is being used for the analysis and understanding of simple experiments using cold atoms, and for the design of a modular compact source of cold atoms to be used in future research and development projects. The results so far indicate that the code is a useful design instrument that shows good agreement with experimental measurements (see figure), and a Windows-based user-friendly interface is also under development.

Ramirez-Serrano, Jaime↗

DOE OSTI2024

Giant resistance switch in twisted transition metal dichalcogenide tunnel junctions

Resistance switching in multilayer structures are typically based on materials possessing ferroic orders. Here, in this work, we predict an extremely large resistance switching based on the relative spin–orbit splitting in twisted transition metal dichalcogenide (TMD) monolayers tunnel junctions. Because of the valence band spin splitting which depends on the valley index in the Brillouin zone, the perpendicular electronic transport through the junction depends on the relative reciprocal space overlap of the spin-dependent Fermi surfaces of both layers, which can be tuned by twisting one layer. Our quantum transport calculations reveal a switching resistance larger than 10 6 % when the relative alignment of TMDs goes from 0° to 60° and when the angle is kept fixed at 60° and the Fermi level is varied. By creating vacancies, we evaluate how inter-valley scattering affects the efficiency and find that the resistance switching remains large (10 4 %) for typical values of vacancy concentration. Not only should this resistance switching be observed at room temperature due to the large spin splitting, but our results also show how twist angle engineering and control of van der Waals heterostructures could be used for next-generation memory and electronic applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Slow Quasiparticle Dynamics and Anyonic Statistics in a Fractional Quantum Hall Fabry-Pérot Interferometer

Anyons are two-dimensional particles with fractional exchange statistics that emerge as elementary excitations of fractional quantum Hall phases. Experimentally, their exchange statistics can be measured in the edge-state Fabry-Pérot interferometer, wherein the presence of 𝑁 𝑞⁢𝑝 localized anyons contributes a phase 𝑁 𝑞⁢𝑝⁢ 𝜃 𝑎 to the interference pattern where 𝜃 𝑎 is twice the exchange phase. Here we report the observation of large, hysteretic phase jumps in a monolayer graphene Fabry-Pérot interferometer at 𝜈 = 1/3. When the filling factor is increased from 𝜈 < 1/3 toward the center of the plateau, we observe phase slips with magnitude Δ⁢𝜃 ≈ 2⁢𝜋/3, consistent with the addition of individual quasiparticles to the interferometer bulk. These phase slips occur as instantaneous jumps in the interference signal, with intervals between the jumps indicating quasiparticle equilibration times exceeding 20 min. We use this long timescale to investigate the effect of changes in interferometer area 𝐴 𝐼 and 𝑁 𝑞⁢𝑝 independently at fixed magnetic field, revealing a striking memory effect in the phase slip magnitude. In particular, as the 𝜈 =1/3 plateau is approached from higher filling, we observed phase slips with Δ⁢𝜃 significantly larger than 2⁢𝜋/3 over the same range of gate voltage where quantized jumps are seen for increasing 𝜈. We discuss this asymmetry in terms of bulk-edge coupling of quasiparticles localized near the edge or in the bulk, and argue that this effect can be qualitatively reconciled with theoretical expectations for strongly interacting quasiparticles in the presence of weak disorder and strongly nonequilibrium charge dynamics. Besides providing a replication of interferometric measurements sensitive to 𝜃 𝑎 , our results highlight the key role played by charge dynamics on signatures of the anyon phase, and demonstrate that fractional quasiparticles can be indefinitely localized in nonequilibrium configurations.

anyons↗

NASA NTRS2000

Far Ultraviolet Imaging from the Image Spacecraft: Wideband FUV Imaging - 2

The Far Ultraviolet Wideband Imaging Camera (WIC) complements the magnetospheric images taken by the IMAGE satellite instruments with simultaneous global maps of the terrestrial aurora. Thus, a primary requirement of WIC is to image the total intensity of the aurora in wavelength regions most representative of the aurora] source and least contaminated by dayglow, have sufficient field of view to cover the entire polar region from spacecraft apogee and have resolution that is Sufficient to resolve auroras on a scale of 1 to 2 latitude degrees, The instrument is sensitive in the spectral region from 140- 190 nm. The WIC is mounted on the rotating, IMAGE spacecraft viewing radially outward and has a field of view of 17 deg in the direction parallel to the spacecraft spin axis. Its field of view is 30 deg in the direction perpendicular to the spin axis, although only a 17 deg x 17 deg image of the Earth is recorded. The optics was an all-reflective, inverted Cassegrain Burch camera using concentric optics with a small convex primary and a large concave secondary mirror. The mirrors were coated by a special multi-layer coating, which has low reflectivity in the visible and near UV region, The detector consists of a MCP-Intensified CCD. The MCP is curved to accommodate the focal surface of the concentric optics. Tile phosphor of the image intensifier is deposited on a concave fiberoptic window, which is then Coupled to the CCD with a fiberoptic taper. The camera head operates in a fast frame transfer mode with the CCD being read approximately 30 full frames (512 by 256 pixel) per second with an exposure time of 0.033 s. The image motion (file to the satellite spin is minimal during such a short exposure. Each image is electronically distortion corrected using the look up table scheme. An offset is added to each memory address that is proportional to the image shift due to satellite rotation, and the charge signal is digitally summed in memory. On orbit, approximately 300 frames will be added to produce one WIC image in memory. The advantage of the electronic motion compensation and distortion correction is that it is extremely flexible, permitting several kinds of corrections including motions parallel and perpendicular to the predicted axis of rotation. File instrument was calibrated by applying ultraviolet light through a vacuum monochromator and measuring the absolute responsivity of the instrument. To obtain the data for the distortion look up table the camera was turned through various angles and the input angles corresponding to a pixel matrix were recorded. It was found that the spectral response peaked at 150 nm and fell off in either direction. The equivalent aperture of the camera, including mirror reflectivities and effective photocathode quantum efficiency, is about 0.04 sq cm. Thus, a 100 Rayleigh LBH aurora is expected to produce 23 equivalent counts per pixel per 10 s exposure at the peak of instrument response.

Mende, S. B.↗

DOE OSTI2023

Strong Axiality in a Dysprosium(III) Bis(borolide) Complex Leads to Magnetic Blocking at 65 K

Substituted dysprosocenium complexes of the type [Dy(Cp R ) 2 ] + exhibit slow magnetic relaxation at cryogenic temperatures and have emerged as top-performing single-molecule magnets. The remarkable properties of these compounds derive in part from the strong axial ligand field afforded by the cyclopentadiene anions, and the design of analogous compounds with even stronger ligand fields is one promising route toward identifying new single-molecule magnets that retain a magnetic memory at even higher temperatures. Here, we report the synthesis and characterization of a dysprosium bis(borolide) compound, [K(18-crown-6)][Dy(BC 4 Ph 5 ) 2 ] (1), featuring the dysprosocenate anion [Dy(BC 4 Ph 5 ) 2 ] - with a pseudoaxial coordination environment afforded by two dianionic pentaphenyl borolide ligands. Variable-field magnetization data reveal open magnetic hysteresis up to 66 K, establishing 1 as a top-performing single-molecule magnet among its dysprosocenium analogues. Ac magnetic susceptibility data indicate that 1 relaxes via an Orbach mechanism above ~80 K with U eff = 1500(100) cm -1 and τ 0 = 10 -12.0(9) s, whereas Raman relaxation and quantum tunneling of the magnetization dominate at lower temperatures. Compound 1 exhibits a 100 s blocking temperature of 65 K, among the highest reported for dysprosium-based single-molecule magnets. Ab initio spin dynamics calculations support the experimental U eff and τ 0 values and enable a quantitative comparison of the relaxation dynamics of 1 and two representative dysprosocenium cations, yielding additional insights into the impact of the crystal field splitting and vibronic coupling on the observed relaxation behavior. Importantly, compound 1 represents a step toward the development of alternatives to substituted dysprosocenium single-molecule magnets with increased axiality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

System and method for controlling an energy beam of an additive manufacturing system

An additive manufacturing system includes an energy beam generator that generates an energy beam to melt and fuse raw materials for a part and a computing device that controls operation of the energy beam generator. The computing device includes a memory element that stores or access a three-dimensional model of the part, and a processing element that receives at least a portion of the three-dimensional model and controls a parameter of the energy beam generator according to a gradient function so as to apply a variable amount of heat to the raw material that forms an overhang on the part.

42 ENGINEERING↗

DOE OSTI2020

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Magnetic Blocking in Fluoflavine Radical-Bridged Dilanthanide Complexes

Magnetic exchange coupling is difficult to foster in polynuclear lanthanide (Ln) complexes and poorly understood. While coupling Ln ions through closed-shell ligands is inherently weak due to the contracted 4f orbitals, placing open-shell ligands instead has proven to promote orders of magnitude stronger coupling, giving rise to single-molecule magnets (SMMs) innate to real magnetic memory effect in the case of the anisotropic Ln ions. Notably, the impact of radical bridges with differing oxidation states on magnetic blocking remains unexplored due to lack of Ln SMMs with radicals in two distinct charge states. Herein, the first dilanthanide complexes (Ln = Gd, Dy) containing fluoflavine (flv) bridges, [(Cp* 2 Ln) 2 (μ-flv z )]X, (where X = [Al(OC{CF 3 } 3 ) 4 ] − (z = 1−•), 1-Ln; X = 0 (z = 2−), 2-Ln; X =[K(crypt-222)] + (z = 3−•), 3-Ln) are reported. 1-Ln and 3-Ln, comprising the flv 1−• and flv 3−• radical bridges, were investigated via single-crystal X-ray diffraction (SCXRD), ultraviolet−visible (UV−vis) spectroscopy, Superconducting Quantum Interference Device (SQUID) magnetometry, high-field electron paramagnetic resonance (HF-EPR) spectroscopy and broken-symmetry density functional theory (BS-DFT) calculations. 1-Dy and 3-Dy constitute the first SMMs innate to radicals in two differing oxidation states. 1-Dy exhibits a spin-reversal barrier U eff of 28.36 cm −1 and open magnetic hysteresis loops below 3 K. By contrast, 3-Dy displays a much higher U eff of 143(2) cm −1 and open hysteresis loops until 9.5 K, representing a record for dilanthanide SMMs containing an organic radical bridge. The boost in SMM properties in 3-Dy is attributed to spin-phonon coupling and improved frontier orbital structure. This study paves the way for advanced design strategies of polynuclear Ln SMMs.

Benner, Florian [Michigan State Univ., East Lansin↗

DOE OSTI2022

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Experiments on the single-mode Richtmyer–Meshkov instability with reshock at high energy densities

Here, the hydrodynamic instability growth of a reshocked single-mode interface between high energy density fluids is studied. A laser-driven shock wave is used to drive an initially solid, sinusoidal interface between a dense plastic (1.43 g/cc) and a light foam (≈ 0.110 g/cc). After the interface has grown to a nonlinear state where the amplitude is of order of the wavelength, it is reshocked. The reshock compresses the nonlinear perturbation, which then grows at about twice the rate. While the pre-reshock growth rate is sensitive to the initial amplitude and wavelength of the perturbation, the post-reshock growth rate is comparatively insensitive to the initial condition. Qualitatively, we observe that the perturbations are less coherent after reshock, consistent with the idea that having a reshock accelerates the transition to turbulence. We find that some memory of the initial condition remains, even after reshock at late time: it appears if the initial perturbations have large enough wavelengths, and the flow structure of size comparable to the initial wavelength persists through reshock. Our results agree with design simulations and are consistent with the phenomenology of reshock studies in conventional gaseous shock tubes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Spin–Lattice Relaxation Decoherence Suppression in Vanishing Orbital Angular Momentum Qubits

Multifrequency electron paramagnetic resonance spectroscopy on oriented single crystals of magnetically dilute Gd(III) ions in Gd 0.004 Y 0.996 (trensal) is used to determine the Hamiltonian parameters of the ground 8 S 7/2 term and its phase memory time, Tm, characterizing its coherent spin dynamics. Here, the vanishing orbital angular momentum of the 8 S 7/2 term makes it relatively insensitive to spin–lattice relaxation mediated by magnetoelastic coupling and leads to a T m of 12 μs at 3 K, which is not limited by spin–lattice relaxation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Velocity-Space Hybridization of Direct Simulation Monte Carlo and a Quasi-Particle Boltzmann Solver

This paper presents a new method for modeling rarefied gas flows based on hybridization of direct simulation Monte Carlo (DSMC) and discrete velocity method (DVM)-based quasi-particle representations of the velocity distribution function. It is aimed at improving the resolution of the tails of the distribution function (compared with DSMC) and computational efficiency (compared with DVM). Details of the method, such as the collision algorithm and the particle merging scheme, are discussed. The hybrid approach is applied to the study of noise in a Maxwellian distribution, computation of electron-impact ionization rate coefficient, as well as numerical simulation of a supersonic Couette flow. The hybrid-based solver is compared with pure DSMC and DVM approaches in terms of accuracy, computational speed, and memory use. It is shown that such a hybrid approach can provide a lower computational cost than a pure DVM approach, while being able to retain accuracy in modeling high-velocity tails of the distribution function. For problems where trace species have a significant impact on the flow physics, the proposed method is shown to be capable of providing better computational efficiency and accuracy compared with standard fixed-weight DSMC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗