Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum kinetic theory”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

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

Superdiffusion from Emergent Classical Solitons in Quantum Spin Chains

Finite-temperature spin transport in the quantum Heisenberg spin chain is known to be superdiffusive, and has been conjectured to lie in the Kardar-Parisi-Zhang (KPZ) universality class. Using a kinetic theory of transport, we compute the KPZ coupling strength for the Heisenberg chain as a function of temperature, directly from microscopics; the results agree well with density-matrix renormalization group simulations. Here, we establish a rigorous quantum-classical correspondence between the “giant quasiparticles” that govern superdiffusion and solitons in the classical continuous Landau-Lifshitz ferromagnet. We conclude that KPZ universality has the same origin in classical and quantum integrable isotropic magnets: a finite-temperature gas of low-energy classical solitons.

1-dimensional spin chains↗
DOE OSTI2025

The kinetics of the reactions of ground-state NH with N 2 O and implications for ammonia combustion

Quantum chemistry and canonical transition state theory have been applied to derive rate constants k for reactions of 3 Σ − imidogen with nitrous oxide. 3 channels were quantified. The fastest leads to NNH+NO and may be summarized as k = 2.8 × 10 7 T 1.68 exp.(−23.1 kcal mol −1 /RT) cm 3 mole −1 s −1 . Two much slower pathways lead to N 2 + 3 HNO and spin-forbidden 1 HNO. Modeling of literature data, from an NH 3 /O 2 flame and jet-stirred reactor measurements on dilute NH 3 /N 2 O mixtures, with the new kinetic information suggests that these channels are too slow to consume N 2 O significantly, and confirms that a widely used proposed rate constant for NH+N 2 O overestimates the reactivity.

Marshall, Paul [University of North Texas, Denton,↗
DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Liouvillian dynamics of the open Schwinger model: String breaking and kinetic dissipation in a thermal medium

Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an antifermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. Here in this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system’s von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Conformal maps of viscous electron flow in the Gurzhi crossover

Herein we investigate the impact of geometric constriction on the viscous flow of electron liquid through quantum point contacts. We provide analysis on the electric potential distribution given the setup of a slit configuration and use the method of conformal mapping to obtain analytical results. The potential profile can be tested and contrasted experimentally with the scanning tunneling potentiometry technique. We discuss intricate physics that underlies the Gurzhi effect, i.e., the enhancement of conductivity in the viscous flow, and compare results for different boundary conditions. In addition, we calculate the temperature dependence of the momentum relaxation time as a result of impurity assisted quasiballistic interference effects and discuss various correlational corrections that lead to the violation of Matthiessen's rule in the hydrodynamic regime. We caution that spatially inhomogeneous profiles of current in the Gurzhi crossover between Ohmic and Stokes flows might also appear in the nonhydrodynamic limit where nonlocality plays an important role. This conclusion is corroborated by calculation of dispersive conductivity in the weakly impure limit.

2-dimensional systems↗
DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Collisional stopping power of ions in warm dense matter

Here, a model for the collisional stopping of ions on free electrons in warm dense matter is developed and explored. It is based on plasma kinetic theory, but with modifications to address the warm dense matter regime. Specifically, it uses the Boltzmann-Uehling-Uhlenbeck kinetic equation to incorporate effects of Fermi degeneracy of electrons. The cross section is computed from quantum scattering of electrons and ions occurring via the potential of mean force derived from an average atom model, which incorporates effects of strong Coulomb correlations. Predictions from this model show comparable accuracy to results from time-dependent density functional theory calculations for deuterium near solid density and a temperature of several electronvolts, at a fraction of the computational cost. Further, the model captures the transition of a plasma from the classical limit to the degenerate limit, including qualitative behaviors of solid state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mixed stochastic-deterministic time-dependent density functional theory: application to stopping power of warm dense carbon

Warm dense matter (WDM) describes an intermediate phase, between condensed matter and classical plasmas, found in natural and man-made systems. In a laboratory setting, WDM is often created dynamically. It is typically laser or pulse-power generated and can be difficult to characterize experimentally. Measuring the energy loss of high energy ions, caused by a WDM target, is both a promising diagnostic and of fundamental importance to inertial confinement fusion research. However, electron coupling, degeneracy, and quantum effects limit the accuracy of easily calculable kinetic models for stopping power, while high temperatures make the traditional tools of condensed matter, e.g. time-dependent density functional theory (TD-DFT), often intractable. In this study, we have developed a mixed stochastic-deterministic approach to TD-DFT which provides more efficient computation while maintaining the required precision for model discrimination. Recently, this approach showed significant improvement compared to models when compared to experimental energy loss measurements in WDM carbon. Here, we describe this approach and demonstrate its application to warm dense carbon stopping across a range of projectile velocities. We compare direct stopping-power calculation to approaches based on combining homogeneous electron gas response with bound electrons, with parameters extracted from our TD-DFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Theoretical and experimental study on the O( 3 P) + 2,5-dimethylfuran reaction in the gas phase

In this work we report a joint experimental and computational study on the 2,5-dimethylfuran oxidation reaction in the gas phase initiated by atomic oxygen O( 3 P). The experiments have been performed by using vacuum-ultraviolet synchrotron radiation at the Advanced Light Source (ALS) of the Lawrence Berkeley National Laboratory (LBNL), at a temperature of 550 K and a pressure of 8 Torr. The experimental data were supported by quantum-chemical calculations along with a kinetic model, also taking into account the possible involvement of different magnetic states, performed in the framework of the RRKM theory. Propyne, acetaldehyde, methylglyoxal, dimethylglyoxal, 3-penten-2-one, 2,5-dimethylfuran-3(2 H )-one, and 1,2-diacetyl ethylene have been identified as the main primary products arising under the conditions of the experiment. Our computational model suggests that these species can be formed at the concentration and branching ratio experimentally observed only in the presence of a non-negligible fraction of non-thermalized intermediates.

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

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2024

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

Second-order spectral line shift comparisons

The second-order spectral line width formulae from the projection operator and kinetic theory methods were recently compared. It was shown that a systematic expansion of the projection operator width expression including initial correlations formally agrees with the second-order kinetic theory result. It is now shown that the second-order dynamic shifts are also formally the same. The static shifts, however, differ due to an ad hoc treatment of electron-electron correlations in the projection operator method. The approximation is necessary in order to screen the radiator-electron interactions. The differences, however, are expected to be small. Finally, the results suggest using the rigorous and more compact second-order width and shift expressions from the kinetic theory method as the starting point for spectral line shape calculations. At line center, however, the projection operator second-order expression for the width and shift simplifies and reduces to the kinetic theory result.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nanocrystal Assemblies: Current Advances and Open Problems

Here we explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.

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

Transport coefficients of quasiparticle models within a new relaxation time approximation of the Boltzmann equation

We investigate the transport properties of a kinetic theory model that is tuned to describe the thermodynamic properties of QCD at zero chemical potential using a new formulation of the relaxation time approximation. In contrast to previous approaches, the latter is constructed to preserve the fundamental properties of the collision term of the Boltzmann equation for any energy dependence of the relaxation time. A novel choice of matching conditions is implemented to ensure that the background mean-field depends only on the temperature even when the system is out of equilibrium. We provide a consistent analysis of how the transport coefficients of relativistic Navier-Stokes theory vary with the energy dependence of the relaxation time. We also show that the entropy production of this theory is consistent with the second law of thermodynamics and verify that it is independent of the matching conditions employed. We used this fact to calculate the matching independent combination of transport coefficients.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the Finite Nuclear Effect and Gaussian Basis Sets for Four-Component Dirac Hartree-Fock Calculations

In this perspective we describe an overview of Dirac Hartree-Fock theory, which is the foundation of relativistic quantum chemistry methods vital in the molecular modeling of heavy element chemistry and physics. A special emphasis was placed upon the physical connection between the finite nuclear model and the kinetic-balance condition between the large and small components. Additionally, we review the consequences of basis set choices for Dirac Hartree-Fock calculations in the restricted kinetic balance formalism, and show that the basis set limit does not have a hard lower bound due to the nature of the Hamiltonian, referred to in the literature as prolapse. Standard Gaussian basis set contraction schemes can break down in Dirac Hartree-Fock calculations, and special care must be taken for the description of core orbitals near the nucleus. The phenomena described in this perspective are vital for theoreticians and practitioners to understand when selecting a basis set for heavy element quantum chemical calculations using Dirac Hartree-Fock theory.

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

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗