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

Toward alcohol synthesis from CO hydrogenation on Cu(111)-supported MoS 2 – predictions from DFT+KMC

In the quest for cheap and efficient catalysts for alcohol synthesis from syngas, a material of interest is single-layer MoS2 owing to its low cost, abundancy, and flexible structure. Because of the inertness of its basal plane, however, it is essential to find ways that make it catalytically active. Herein, by means of density functional theory based calculations of reaction pathways and activation energy barriers and accompanying kinetic Monte Carlo simulations, we show that while S vacancy row structures activate the MoS 2 basal plane, further enhancement of chemical activity and selectivity can be achieved by interfacing the MoS 2 layer with a metallic support. When defect-laden MoS2 is grown on Cu(111), there is not only an increase in the active region (surface area of active sites) but also charge transfer from Cu to MoS 2 , resulting in a shift of the Fermi level such that the frontier states (d orbitals of the exposed Mo atoms) appear close to it, making the MoS 2 /Cu(111) system ready for catalytic activity. Finally, our calculated thermodynamics of reaction pathways lead to the conclusion that the Cu(111) substrate promotes both methanol and ethanol as the products, while kinetic Monte Carlo simulations suggest a high selectivity toward the formation of ethanol.

2D materials↗
DOE OSTI2020

Effect of finite ion mass on relativistic self-induced transparency of plasma layers with a sharp boundary

We consider the effect of ion mobility on the threshold of relativistic self-induced transparency under irradiation of thin plasma layers by circularly polarised laser radiation. An analytical model of the motion of ions during the removal of electrons from the layer surface by the ponderomotive force of laser radiation is constructed. The model is used to analyse the motion of probe electrons in the resulting electromagnetic field. It is shown that the higher the ion mobility and the longer the laser pulse, the more stable the plasma layer to longitudinal perturbations of a finite value and the higher the threshold of self-induced transparency. These conclusions are verified by one-dimensional numerical simulation of a complete system of kinetic equations for plasma and Maxwell’s equations. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

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

Understanding how chemical structure affects ignition-delay-time $\phi$-sensitivity

$\phi$-sensitivity is the change in ignition delay time (IDT) with respect to the fuel-to-air equivalence ratio ($\phi$). High $\phi$-sensitivity is a desirable fuel property for applications in advanced compression ignition and multi-mode engine designs. Understanding how $\phi$-sensitivity depends on chemical structure is essential for selecting promising biofuels from the ever-growing list of proposed candidates. Here, we investigate the effect of chemical structure on $\phi$-sensitivity with experiment, simulation, and theory. Experimental Advanced Fuel Ignition Delay Analyzer (AFIDA) measurements for 2,4-dimethylpentane and diisopropyl ether provide evidence that branching and functional groups strongly impact $\phi$-sensitivity. Further insights into this dependence are obtained with 0-D kinetic simulations with existing mechanisms for n-pentane, diethyl ether, 3-pentanone, n-heptane, 2-methylhexane, 2,4-dimethylpentane, and 2,2,3-trimethylbutane. Quantum mechanical (QM) G4 calculations of low-temperature reactions help explain the observed experimental and simulation trends. Specifically, these QM calculations provide theoretical estimates of the ketohydroperoxide (KHP) dissociation rates, the HO 2 formation rates from peroxy radical (ROO), and the “cross-over” temperatures, i.e., the temperature at which ROO dissociation is favored compared to hydroperoxyl radical (QOOH) formation. Each of these reaction rates is compared to the n-alkane reference point to determine the impact of branching and different functional groups. Although kinetic mechanisms typically assume that KHP dissociation rates are invariant of chemical environment, our QM results suggest that this rate can span a range of roughly two orders of magnitude. We also discuss the importance of including the peroxy-hydroperoxy (OO-OOH) hydrogen transfer reaction for branched ethers. Finally, the insights gained assist in proposing a highly $\phi$-sensitive compound, namely, isopropyl propyl ether.

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

Energetics and kinetics of various cyano radical hydrogen abstractions

The cyano radical (CN) is an abundant, open-shell molecule found in a variety of environments, including the atmosphere, the interstellar medium and combustion processes. In these environments, it often reacts with small, closed-shell molecules via hydrogen abstraction. Both carbon and nitrogen atoms of the cyano radical are reactive sites, however the carbon is more reactive with reaction barrier heights generally between 2–15 kcal mol -1 lower than those of the analogous nitrogen. The CN + HX → HCN/HNC + X, with X = H, CH 3 , NH 2 , OH, F, SiH 3 , PH 2 , SH, Cl, C 2 H, CN reactions have been studied at a high-level of theory, including CCSD(T)-F12a. Finally, kinetics were obtained over the 100–1000 K temperature range, showing excellent agreement with those rate constants that have been determined experimentally.

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

Kinklike structures in an arcsin real scalar dynamics

Highlights: • Kinklike analyses in a real scalar theory with an arcsin dynamic. • First-order formalism implementation and deformation procedure in an arcsin scalar dynamics. In this paper, we analyze kink-like analytical solutions in a real scalar theory with an arcsin dynamics inspired by the arcsin electrodynamics presented in Kruglov (2015). This analysis is done by means of the first-order formalism. This formalism provides a framework where the equations of motion can be simplified by preserving the linear stability of the theory. In this work, the deformation procedure is implemented with the aim of finding exact solutions in systems with generalized dynamics. Along the paper, we explore how the first-order formalism is implemented in the arcsin kinetics and how such a term influences the kink-like solutions. As a part of the result of our paper, we show that the kink-like solutions are similar to the ones obtained in the standard scalar kinetic theory. We also show that the extra parameter, that controls the non-linearities of the model, plays an essential role in the energy densities and stability potentials. These quantities vary according to this parameter. The goals here are to show how the first-order framework is implemented in this arcsin scenario and to present the analytical kink-like solutions that can be found by means of the first-order framework and deformation method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2023

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

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

Timescales of Quantum Equilibration, Dissipation and Fluctuation in Nuclear Collisions

Understanding the dynamics of equilibration processes in quantum systems as well as their interplay with dissipation and fluctuation is a major challenge in quantum many-body theory. The timescales of such processes are investigated in collisions of atomic nuclei using fully microscopic approaches. Findings from time-dependent Hartree-Fock and time-dependent random-phase approximation calculations are compared for 13 systems over a broad range of energies. The timescale for full mass equilibration ( ~ 2 × 10 - 20 s ) is found to be much larger than timescales for neutron-to-proton equilibration, kinetic energy, and angular momentum dissipations which are on the order of 10 - 21 s . Fluctuations of mass numbers in the fragments and correlations between their neutron and proton numbers build up within only a few 10 - 21 s . This suggests that dissipation is basically not impacted by mass equilibration, but is mostly driven by the exchange of nucleons between the fragments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2023

Mixed quantum/classical theory for rotational energy exchange in symmetric-top-rotor + linear-rotor collisions and a case study of the ND 3 + D 2 system

The extension of mixed quantum/classical theory (MQCT) to describe collisional energy transfer is developed for a symmetric-top-rotor + linear-rotor system and is applied to ND 3 + D 2 . State-to-state transition cross sections are computed in a broad energy range for all possible processes: when both ND 3 and D 2 molecules are excited or both are quenched, when one is excited while the other is quenched and vice versa, when the ND 3 state changes its parity while D 2 is excited or quenched, and when ND 3 is excited or quenched while D 2 remains in the same state, ground or excited. In all these processes the results of MQCT are found to approximately satisfy the principle of microscopic reversibility. For a set of sixteen state-to-state transitions available from the literature for a collision energy of 800 cm –1 the values of cross sections predicted by MQCT are within 8% of accurate full-quantum results. A useful time-dependent insight is obtained by monitoring the evolution of state populations along MQCT trajectories. It is shown that, if before the collision, D 2 is in its ground state, the excitation of ND 3 rotational states proceeds through a two-step mechanism in which the kinetic energy of molecule–molecule collision is first used to excite D 2 and only then is transferred to the excited rotational states of ND 3 . Finally, it is found that both potential coupling and Coriolis coupling play important roles in ND 3 + D 2 collisions.

inelastic cross-sections↗
DOE OSTI2020

Theoretical H + O 3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O 3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. Additionally, the theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O 3 → OD + O 2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

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

Catalytic resonance theory for the kinetic signatures of multiple wavelength photocatalysis

Multi-step surface chemistry promoted with one or two wavelengths of light exhibited distinct kinetic signatures indicative of the light-adsorbate interactions and the number of photon-sensitive elementary steps. In this work, kinetic Monte Carlo simulations identified the kinetic response of general surface mechanisms to variation in the per-site photon flux of one or two wavelengths of incident light. Photocatalytic rates were described via a non-dimensional photon flux, identifying multiple kinetic regimes with unique degrees of rate control. Under photon-controlled kinetic conditions, maximum quantum yield and turnover frequency were obtained under constant illumination, while pulsing of one or more light sources was shown to exhibit slower rates and less efficient photocatalytic promotion due to the inherent dynamic nature of light, which comprises a stream of photons. Simulations provided distinct kinetic signatures in catalytic rates and Arrhenius kinetics for specific light-surface-adsorbate interactions corresponding to photocatalytic promotion of specific steps in surface chemistry.

catalytic resonance↗
DOE OSTI2025

Kinetic Isotope Effect in the Radical Step-Growth Termination of On-Surface Synthesized Graphene Nanoribbons

Bottom-up on-surface synthesis has emerged as a versatile tool to access and finely tune the electronic structure of nanographenes. The controlled generation of reactive intermediates catalyzed and stabilized by a supporting substrate has enabled the design, assembly, and characterization of a wide range of exotic tailor-made quantum materials. Even under these tightly controlled conditions, the growth of extended structures remains limited by termination processes and undesired side reactions. Here, we identify an H atom transfer as one principal contributor to the radical step growth termination that limits the on-surface growth of N = 7 armchair graphene nanoribbons (7-AGNRs) on Au(111) surfaces. Analysis of 7-AGNR lengths grown from protiated and deuterated molecular precursors reveals a primary kinetic isotope effect of KIE ∼ 1.4. First-principles density functional theory calculations suggest that a concerted H atom transfer mechanism that involves the breaking of a C–H/D bond in the transition state is associated with radical chain termination.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
DOE OSTI2022

Electrochemical Hydrogenation of CO on Cu(100): Insights from Accurate Multiconfigurational Wavefunction Methods

Copper (Cu) remains the most efficacious electrocatalyst for electrochemical CO 2 reduction (CO 2 R). Its activity and selectivity are highly facet-dependent. We recently examined the commonly proposed rate-limiting CO hydrogenation step on Cu(111) via embedded correlated wavefunction (ECW) theory and demonstrated that only this higher-level theory yields predictions consistent with potential-dependent experimental kinetics. Here, to understand the differing activities of Cu(111) and Cu(100) in catalyzing CO 2 R, we explore CO hydrogenation on Cu(100) using ECW theory. We predict that the preferred pathway involves the reduction of adsorbed CO (*CO) to *COH via proton-coupled electron transfer (PCET) at working potentials, although *CHO also may form with a kinetically accessible but higher barrier. In contrast, our earlier work on Cu(111) concluded that *COH and *CHO formation via PCET are equally feasible. Furthermore, this work illustrates one possible origin of the facet dependence of CO 2 R mechanisms and products on Cu electrodes and sheds light on how the selectivity of CO 2 R electrocatalysts can be controlled by the surface morphology.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Order-by-order uncertainties of nucleon-nucleon Wolfenstein amplitudes in chiral effective field theory

Quantum mechanical invariance principles dictate the most general operator structure that can be present in the nucleon-nucleon (NN) interaction. Five independent operators appear in the on-shell NN amplitude together with five corresponding coefficient functions. The usual choice for these coefficient functions is known as the NN Wolfenstein amplitudes. We analyze the order-by-order convergence of each of the five NN Wolfenstein amplitudes predicted by a semilocal coordinate space potential implementation of chiral effective field theory (𝜒⁢EFT). We do this at laboratory kinetic energies between 25 and 200 MeV for both neutron-proton and proton-proton scattering. Our analysis uses the Gaussian-process methods developed by the BUQEYE Collaboration to describe the contributions of each 𝜒⁢EFT order, and so yields truncation uncertainties for each Wolfenstein amplitude that are correlated across scattering angles. We combine information on the size of different orders in the EFT to infer the 𝜒⁢EFT breakdown scale for each amplitude, finding, on average, Λ 𝑏 between 750 and 800 MeV. Furthermore, with this choice of Λ 𝑏 , the EFT truncation uncertainties cover both higher-order results and empirical Wolfenstein amplitudes well for all orders other than the leading order.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A thermodynamic analysis of hypervelocity impacts on metals

Hypervelocity impacts on solid surfaces generate and propagate shock waves and create extreme density, pressure, and temperature conditions. Understanding hypervelocity impacts requires a multi-physics approach at the intersection of fluid dynamics, solid mechanics, and quantum theory, which has the capability of spanning multiple different dominating physics regimes. These strong shock conditions found in hypervelocity impacts are studied in the context of high-speed collisions between spacecraft and meteoroids and debris. We study the problem of hypervelocity impacts on metals with the lens of understanding the processes that lead to the creation of an impact generated plasma. In order to characterize the impacting particle as well as its threat to spacecraft from the subsequent radiating plasma an analysis was performed that highlights how the kinetic energy of the impactor is deposited into a target, how the shock wave and shocked material in the target evolve, and how these effects can lead towards the creation of an impact plasma through heating. The analysis suggests that at low velocities the majority of the energy from a hypervelocity impact is stored in reversible elastic energy and does not contribute to strong heating. At higher impact velocities, nucleus and electron thermal components dominate, increasing the energy available for ionization. Overall, this work provides the complete set of state variable initial conditions for any study of hypervelocity impact plasma expansion and its effects.

42 ENGINEERING↗