Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical calculations”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2020

Reactants, products, and transition states of elementary chemical reactions based on quantum chemistry

Reaction times, activation energies, branching ratios, yields, and many other quantitative attributes are important for precise organic syntheses and generating detailed reaction mechanisms. Often, it would be useful to be able to classify proposed reactions as fast or slow. However, quantitative chemical reaction data, especially for atom-mapped reactions, are difficult to find in existing databases. Therefore, we used automated potential energy surface exploration to generate 12,000 organic reactions involving H, C, N, and O atoms calculated at the ω B97X-D3/def2-TZVP quantum chemistry level. We report the results of geometry optimizations and frequency calculations for reactants, products, and transition states of all reactions. Additionally, we extracted atom-mapped reaction SMILES, activation energies, and enthalpies of reaction. We believe that this data will accelerate progress in automated methods for organic synthesis and reaction mechanism generation—for example, by enabling the development of novel machine learning models for quantitative reaction prediction.

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

SA-GAT-SR: self-adaptable graph attention networks with symbolic regression for high-fidelity material property prediction

Recent advances in machine learning have demonstrated an enormous utility of deep learning approaches, particularly Graph Neural Networks (GNNs) for materials science. These methods have emerged as powerful tools for high-throughput prediction of material properties, offering a compelling enhancement and alternative to traditional first-principles calculations. While the community has predominantly focused on developing increasingly complex and universal models to enhance predictive accuracy, such approaches often lack physical interpretability and insights into materials behavior. Here, we introduce a novel computational paradigm—Self-Adaptable Graph Attention Networks integrated with Symbolic Regression (SA-GAT-SR)—that synergistically combines the predictive capability of GNNs with the interpretative power of symbolic regression. Our framework employs a self-adaptable encoding algorithm that automatically identifies and adjust attention weights so as to screen critical features from an expansive 180-dimensional feature space while maintaining O(n) computational scaling. The integrated SR module subsequently distills these features into compact analytical expressions that explicitly reveal quantum-mechanically meaningful relationships, achieving 23 × acceleration compared to conventional SR implementations that heavily rely on first-principle calculations-derived features as input. This work suggests a new framework in computational materials science, bridging the gap between predictive accuracy and physical interpretability, offering valuable physical insights into material behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Large ordered moment with strong easy-plane anisotropy and vortex-domain pattern in the kagome ferromagnet Fe 3 Sn

We report the magnetic anisotropy of kagome bilayer ferromagnet Fe 3 Sn probed by the bulk magnetometry and magnetic force microscopy (MFM) on high-quality single crystals. The dependence of magnetization on the orientation of the external magnetic field reveals strong easyplane magnetocrystalline anisotropy and anisotropy of the saturation magnetization. The leading magnetocrystalline anisotropy constant shows a monotonous increase from $\mathrm{K_1$≈ –1.0 \times 10^6 J/m^3}$ at 300 K to $\mathrm{–1.3 \times 10^6 J/m^3}$ at 2 K. Our ab initio electronic structure calculations yield the value of total magnetic moment of 7.1$μ_Β$/$f.u.$ and a magnetocrystalline anisotropy energy density of –0.57 $\mathrm{meV/}$$f.u.$ $\mathrm{(–1.62 \times 10^6J/m^3)}$ both being in reasonable agreement with the experimental values. The MFM imaging reveals micrometer-scale magnetic vortices with weakly pinned cores that vanish at the saturation field of ~3T applied perpendicular to the kagome plane. The observed vortex-domain structure is well reproduced by the micromagnetic simulations, using the experimentally determined value of the anisotropy and exchange stiffness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Full Counting Statistics of Charge in Chaotic Many-Body Quantum Systems

We investigate the full counting statistics of charge transport in U(1)-symmetric random unitary circuits. We consider an initial mixed state prepared with a chemical potential imbalance between the left and right halves of the system and study the fluctuations of the charge transferred across the central bond in typical circuits. Using an effective replica statistical mechanics model and a mapping onto an emergent classical stochastic process valid at large on-site Hilbert space dimension, we show that charge transfer fluctuations approach those of the symmetric exclusion process at long times, with subleading t –1/2 quantum corrections. Here, we discuss our results in the context of fluctuating hydrodynamics and macroscopic fluctuation theory of classical nonequilibrium systems and check our predictions against direct matrix-product state calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels

Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

High reactivity of mesoporous CeO 2 to dissociate chemical warfare agent sarin

We perform spectroscopic surface characterization on mesoporous ceria (CeO 2 ) under exposure of sarin gas, also known as GB. We find mesoporous CeO 2 displays extraordinary reactivity towards GB dissociation under high vacuum. The spectroscopic results combined with density functional theory (DFT) calculations suggest that the cleavage of the P–F and P–OR bonds are both possible pathways for GB to dissociate on the CeO 2 surfaces. The modeling reveals differences in dissociation pathways of GB on the crystallographic (111)- and (110)-surfaces of ceria. Here we also report some dissociation of GB on mesoporous CeO 2 under ambient conditions. The finding that mesoporous CeO 2 can degrade GB under mild conditions is encouraging and implies the strong potential of ceria-based materials to be applied in combating chemical warfare agents related to the GB series.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1972

Quantum statistical mechanics of dense partially ionized hydrogen

The theory of dense hydrogen plasmas beginning with the two component quantum grand partition function is reviewed. It is shown that ionization equilibrium and molecular dissociation equilibrium can be treated in the same manner with proper consideration of all two-body states. A quantum perturbation expansion is used to give an accurate calculation of the equation of state of the gas for any degree of dissociation and ionization. The statistical mechanical calculation of the plasma equation of state is intended for stellar interiors. The general approach is extended to the calculation of the equation of state of the outer layers of large planets.

Dewitt, H. E.↗
DOE OSTI2022

The LPM effect in sequential bremsstrahlung: 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ corrections

An important question concerning in-medium high-energy parton showers in a quark-gluon plasma or other QCD medium is whether consecutive splittings of the partons in a given shower can be treated as quantum mechanically independent, or whether the formation times for two consecutive splittings instead have significant overlap. Various previous calculations of the effect of overlapping formation times have either (i) restricted attention to a soft bremsstrahlung limit, or else (ii) used the large-N c limit (where N c = 3 is the number of quark colors). In this paper, we make a first study of the accuracy of the large-N c limit used by those calculations of overlap effects that avoid a soft bremsstrahlung approximation. Specifically, we calculate the 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ correction to previous N c = ∞ results for overlap g → gg → ggg of two consecutive gluon splittings g → gg. At order 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ , there is interesting and non-trivial color dynamics that must be accounted for during the overlap of the formation times.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

From PAH to fullerenes: Closing a nano bowl with a nano lid

Quantum chemical DLPNO-CCSD(T)/cc-pVDZ//B3LYP/6–31G(d) calculations have been performed to unravel the formation mechanism of buckminsterfullerene (C 60 ) via the reaction of the C 40 nano bowl (C 40 H 10 ) and corannulene (C 20 H 10 ). The generated potential energy surfaces and molecular properties were further utilized to evaluate equilibrium constants and rate constants of elementary chemical reactions involved in the C 60 synthesis. The overall C 40 H 10 + C 20 H 10 → C 60 + 10H 2 reaction is shown to be highly exoergic and hence thermodynamically favorable at temperatures above 700 K and a kinetically feasible pathway under high-temperature conditions was revealed. This route involves hydrogen atom abstraction steps to activate closed-shell reactants/intermediates alternating with cyclodehydrogenation reactions eventually zipping the edges between the reacting C 40 and C 20 units by consecutive closures of six- and five-membered rings between them. The reaction is driven/catalyzed by hydrogen atoms which carry out the hydrogen abstractions from C 60 H x and are later recovered on the cyclodehydrogenation stages. The initial association reaction, C 40 H 9 + C 20 H 10 → C 60 H 18 + H, linking the two units by a covalent C-C bond appears to represent the kinetic bottleneck of the whole multistep process. This reaction is fast at low temperatures but slows down and has its equilibrium shifted toward the reactants at high temperatures. Here, the feasibility of the proposed mechanism corroborates the hypothesis that C 60 can be produced through a series of elementary reactions typical for the PAH growth.

Equilibrium constant↗
NASA NTRS1978

A semiclassical approach to collisional ionization with application to the Ar-He system

A semiclassical approach based on the propagation of classical trajectories on potential surfaces analytically continued into the complex plane, together with a discretization procedure, has been developed for the problem of collisional ionization. Based on Franck-Condon considerations the formalism is reduced to that of the two-state approximation. Preionization loss and tunnelling beyond turning points have not been considered. Calculated partial ionization cross sections for the Ar-He system show good agreement with a fully quantum mechanical treatment.

Lam, K.-S.↗
NASA NTRS2013

Experimental Investigation of Charging Properties of Interstellar Type Silica Dust Grains by Secondary Electron Emissions

The dust charging by electron impact is an important dust charging processes in astrophysical and planetary environments. Incident low energy electrons are reflected or stick to the grains charging the dust grains negatively. At sufficiently high energies electrons penetrate the grains, leading to excitation and emission of electrons referred to as secondary electron emission (SEE). Available classical theoretical models for calculations of SEE yields are generally applicable for neutral, planar, or bulk surfaces. These models, however, are not valid for calculations of the electron impact charging properties of electrostatically charged micron/submicron-size dust grains in astrophysical environments. Rigorous quantum mechanical models are not yet available, and the SEE yields have to be determined experimentally for development of more accurate models for charging of individual dust grains. At the present time, very limited experimental data are available for charging of individual micron-size dust grains, particularly for low energy electron impact. The experimental results on individual, positively charged, micron-size lunar dust grains levitated carried out by us in a unique facility at NASA-MSFC, based on an electrodynamic balance, indicate that the SEE by electron impact is a complex process. The electron impact may lead to charging or discharging of dust grains depending upon the grain size, surface potential, electron energy, electron flux, grain composition, and configuration (Abbas et al, 2010, 2012). In this paper, we discuss SEE charging properties of individual micron-size silica microspheres that are believed to be analogs of a class of interstellar dust grains. The measurements indicate charging of the 0.2m silica particles when exposed to 25 eV electron beams and discharging when exposed to higher energy electron beams. Relatively large size silica particles (5.2-6.82m) generally discharge to lower equilibrium potentials at both electron energies. These measurements conducted on silica microspheres are qualitatively similar in nature to our previous SEE measurements on lunar Apollo missions dust samples.

Tankosic, D.↗
DOE OSTI2021

Real-time quantum calculations of phase shifts using wave packet time delays

We present a method to extract the phase shift of a scattering process using the real-time evolution in the early and intermediate stages of the collision in order to estimate the time delay of a wave packet. This procedure is convenient when using noisy quantum computers for which the asymptotic out-state behavior is unreachable. We demonstrate that the challenging Fourier transforms involved in the state preparation and measurements can be implemented in 1+1 dimensions with current trapped ion devices and IBM quantum computers. We compare quantum computations of the time delays obtained in the one-particle quantum mechanics limit and the scalable quantum field theory formulation with accurate numerical results. We discuss the finite volume effects in the Wigner formula connecting time delays to phase shifts. The results reported involve two- and four-qubit calculations, and we discuss the possibility of larger scale computations in the near future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations

DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. Here, the flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.

36 MATERIALS SCIENCE↗
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 OSTI2021

Logarithmic criticality in transverse thermoelectric conductivity of the ferromagnetic topological semimetal CoMnSb

Here we report the results of our experimental studies on the magnetic, transport, and thermoelectric properties of the ferromagnetic metal CoMnSb. A sizable anomalous Hall conductivity $σ_{yx}$ and transverse thermoelectric conductivity $α_{yx}$ are found experimentally and comparable in size to the values estimated from first-principles calculation. Our experiment further reveals that CoMnSb exhibits -T log T critical behavior in $α_{yx}$ (T), deviating from Fermi liquid behavior $α_{yx}$ ~ T over a decade of temperature between 10 and 400 K, similar to ferromagnetic Weyl and nodal-line semimetals. Our theoretical calculation for CoMnSb also predicts the - T log T behavior when the Fermi energy locates near the Weyl nodes in momentum space.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electric-field-induced crossover of polarization reversal mechanisms in Al 1–x Sc x N ferroelectrics

Scandium-doped aluminum nitride, Al 1–x Sc x N, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al–Sc–N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al 1–x Sc x N to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms—inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Finally, understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗