Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equations of motion”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Enhancing the accuracy of XPS calculations: Exploring hybrid basis set schemes for CVS-EOMIP-CCSD calculations

Reliable computational methodologies and basis sets for modeling x-ray spectra are essential for extracting and interpreting electronic and structural information from experimental x-ray spectra. In particular, the trade-off between numerical accuracy and computational cost due to the size of the basis set is a major challenge, since molecular orbitals undergo extreme relaxation in the core-hole state. To gain clarity on the changes in electronic structure induced by the formation of a core-hole, the use of sufficiently flexible basis for expanding the orbitals, particularly for the core region, has been shown to be essential. This work focuses on the refinement of core-hole ionized state calculations using the equation-of-motion coupled cluster family of methods through an extensive analysis on the effectiveness of “hybrid” and mixed basis sets. In this investigation, we utilize the CVS-EOMIP-CCSD method in combination and construct hybrid basis sets piecewise from readily available Dunning’s correlation consistent basis sets in order to calculate x-ray ionization energies (IEs) for a set of small gas phase molecules. Our results provide insights into the impact of basis sets on the CVS-EOMIP-CCSD calculations of K-edge IEs of first-row p-block elements. Furthermore, these insights enable us to understand more about the basis set dependence of the core IEs computed and allow us to establish a protocol for deriving reliable and cost-effective theoretical estimates for computing IEs of small molecules containing such elements.

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

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

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

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

A fractional calculus framework for open quantum dynamics: From Liouville to Lindblad to memory kernels

Open quantum systems exhibit dynamics ranging from unitary evolution to irreversible dissipation. While the Gorini–Kossakowski–Sudarshan–Lindblad equation uniquely characterizes Markovian completely positive and trace-preserving (CPTP) evolution, many physical platforms display non-Markovian features such as algebraic relaxation and coherence backflow. Fractional calculus provides a natural way to model such long-memory behavior through power-law temporal kernels introduced by fractional time derivatives. Here, we develop a unified framework that embeds fractional master equations within the broader hierarchy of open-system formalisms. The fractional equation forms a structured subclass of memory-kernel models, reduces to the Lindblad form at unit order, and, through Bochner–Phillips subordination, admits a CPTP representation as an average over Lindblad semigroups. Its resolvent structure further connects fractional dynamics to established non-Markovian approaches, including Nakajima–Zwanzig kernels and hierarchical equations of motion, providing a compact surrogate for long-memory effects. This formulation positions fractional calculus as a rigorous and practical language for modeling non-Markovian quantum dynamics in chemical physics and physical chemistry, providing a CPTP-preserving, computationally efficient surrogate for structured condensed-phase environments where long-time memory and dissipation play a central role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Degenerate coupled-cluster theory

A size-extensive, converging, black-box, ab initio coupled-cluster (ΔCC) ansatz is introduced that computes the energies and wave functions of states from any degenerate or nondegenerate Slater-determinant references with any numbers of α- and β-spin electrons, any patterns of orbital occupancy, any spin multiplicities, and any spatial symmetries. For a nondegenerate reference, it reduces to the single-reference coupled-cluster ansatz. For a degenerate multireference, it is a natural coupled-cluster extension of degenerate Møller–Plesset perturbation (ΔMP) theory. For ionized and electron-attached references, it is a coupled-cluster Green’s function, although the present theory is convergent toward the full-configuration-interaction limits, while the Feynman–Dyson many-body Green’s function (MBGF) theory generally is not. Its single-excitation instance is a projection Hartree–Fock theory as per the Thouless theorem, which may be useful for core ionizations, high-spin states, and possibly electron affinities. Additionally, a new multireference coupled-cluster theory for a general model space is developed. This quasidegenerate coupled-cluster (QCC) theory is exactly converging, but not black-box, and intended for strong correlation. Determinant-based, general-order algorithms of ΔCC and QCC theories are implemented and compared with configuration-interaction (CI) and equation-of-motion coupled-cluster (EOM-CC) theories through octuple excitations and with ΔMP and MBGF theories up to the nineteenth order. An algebraic, optimal-scaling algorithm of the ΔCC theory is computer-synthesized at the levels of single excitations (ΔCCS) and of single and double excitations (ΔCCSD). As a result, the order of performance is QCC ≈ ΔCC > EOM-CC > CI at the same order or QCC ≈ ΔCC > ΔMP > MBGF at the same cost scaling.

Hirata, So [University of Illinois at Urbana-Champ↗
DOE OSTI2020

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Comparison of state-interaction and spinor-representation calculations of spin-orbit coupling within exact two-component coupled-cluster theories

Here a benchmark study of state-interaction and spinor-representation calculations of spin-orbit coupling using the exact two-component Hamiltonians with atomic mean-field integrals (the X2CAMF schemes) at the equation-of-motion coupled-cluster singles and doubles level is reported. We adopt a version of the X2CAMF scheme with spin-orbit integrals correct to first order for the state-interaction calculations and a version correct to infinite order for the spinor-representation calculations. The differences between the state-interaction calculations with minimum active spaces to account for (quasi-)degeneracy and the spinor-representation calculations thus correspond to higher-order spin-orbit contributions. This state-interaction approach with scalar-relativistic effects accurately included in the reference functions and in the spin-orbit integrals is shown to exhibit robust performance for elements across the periodic table. On the other hand, the more rigorous spinor representation shows more rapid convergence with respect to the number of correlated electrons and is the preferred choice for accurate calculations for heavy elements.

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

From synchrotrons for XFELs: the soft x-ray near-edge spectrum of the ESCA molecule

A predictive understanding of soft x-ray near-edge absorption spectra of small molecules is an enduring theoretical challenge and of current interest for x-ray probes of molecular dynamics. We report the experimental absorption spectrum for the electron spectroscopy for chemical analysis (ESCA) molecule (ethyl trifluoroacetate) near the carbon 1s absorption edge between 285-300 eV. Furthermore, the ESCA molecule with four chemically distinct carbon sites has previously served as a theoretical benchmark for photoelectron spectra and now for photoabsorption spectra. We report a simple edge-specific approach for systematically expanding standard basis sets to properly describe diffuse Rydberg orbitals and the importance of triple excitations in equation-of-motion coupled-cluster calculations of the energy interval between valence and Rydberg excitations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Resonant x-ray emission across the L 3 edge of uranium compounds

A narrow bandwidth x-ray beamline and a multi-crystal von Hamos spectrometer were used to record x-ray absorption and x-ray emission (XES) across the uranium L 3 edge of UO 2 , UO 3 , Cs 2 UO 2 Cl 4 , and Cs 2 UCl 6 . Measurements were made over 17150–17250 eV with an instrumental resolution of ∼2 eV. This resolution allowed the use of Lorentzian peak fits to the L 3 N 4 and L 3 N 5 characteristic x-ray fluorescence lines that have ∼12–13 eV lifetime widths. The fluorescence yields of the four compounds display strong white lines near threshold and multiple scattering features at higher energies. The XES spectra were fit with three Lorentzians, two for the L 3 N 4 and L 3 N 5 lines and one for inelastic scattering by 2p–6d excitations and 4d–6d final states. The intensities of the white lines were largely due to the inelastic scattering component. To support the measurements, relativistic equation-of-motion coupled-cluster and restricted active space configuration interaction calculations were performed for uranium core-excitation energies in Cs 2 UO 2 Cl 4 and Cs 2 UCl 6 . These calculations with rigorous treatments of relativistic effects are shown to provide accurate uranium L 3 -edge binding energies, L 3 N 5 emission energies, as well as the energy losses in the inelastic scattering process. The measured and calculated results show variations between compounds with U(IV) oxidation states that contain 6d 0 5f 2 electrons in their ground configurations (UO 2 and Cs 2 UCl 6 ) compared with U(VI) compounds (UO 3 and Cs 2 UO 2 Cl 4 ) with empty 5f and 6d configurations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Unified theory of spin and charge in a ferromagnet

Here, we derive a unified theory of spin and charge degrees of freedom in a ferromagnet. The spin-transfer torque and spin electromotive force are examined from the coarse-grained perspective of collective coordinates. The resulting equations of motion reflect a balance of conservative, gyroscopic (Berry-phase), and dissipative forces. We then expand the space of collective coordinates by adding the electric charge. The adiabatic spin-transfer torque and spin electromotive force (emf) turn out to be a gyroscopic force; their nonadiabatic counterparts are a dissipative force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Spin–orbit torque in a three-fold-symmetric bilayer and its effect on magnetization dynamics

Field-free switching of perpendicular magnetization has previously been observed in an epitaxial L1 1 -ordered CoPt/CuPt bilayer and attributed to spin–orbit torque (SOT) arising from the crystallographic 3m point group of the interface. Using a first-principles nonequilibrium Green's function formalism combined with the Anderson disorder model, we calculate the angular dependence of the SOT in a CoPt/CuPt bilayer and find that the magnitude of the 3m SOT is about 20% of the conventional dampinglike SOT. We further study the magnetization dynamics in perpendicularly magnetized films in the presence of 3m SOT and Dzyaloshinskii–Moriya interaction, using the equations of motion for domain wall dynamics and micromagnetic simulations. For systems where strong interfacial DMI results in the Néel character of domain walls, we find that a very large current density is required to achieve deterministic switching, because reorientation of the magnetization inside the domain wall is necessary to induce the switching asymmetry. For thicker films with relatively weak interfacial DMI and the Bloch character of domain walls, deterministic switching is possible at much smaller current densities, which may explain the recent experimental findings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A classical, non-singular, bouncing universe

We present a model for a classical, non-singular bouncing cosmology without violation of the null energy condition (NEC). The field content is General Relativity plus a real scalar field with a canonical kinetic term and only renormalizable, polynomial-type self-interactions for the scalar field in the Jordan frame. The universe begins vacuum-energy dominated and is contracting at t=-∞. We consider a closed universe with a positive spatial curvature, which is responsible for the universe bouncing without any NEC violation. An Rφ 2 coupling between the Ricci scalar and the scalar field drives the scalar field from the initial false vacuum to the true vacuum during the bounce. The model is sub-Planckian throughout its evolution and every dimensionful parameter is below the effective-field-theory scale M P , so we expect no ghost-type or tachyonic instabilities. This model solves the horizon problem and extends co-moving particle geodesics to past infinity, resulting in a geodesically complete universe without singularities. We solve the Friedman equations and the scalar-field equation of motion numerically, and analytically under certain approximations.

Astronomy & Astrophysics↗
DOE OSTI2021

Ultralocality and slow contraction

Abstract We study the detailed process by which slow contraction smooths and flattens the universe using an improved numerical relativity code that accepts initial conditions with non-perturbativedeviations from homogeneity and isotropy along two independent spatial directions.Contrary to common descriptions of the early universe,we find that the geometry first rapidly converges to an inhomogeneous, spatially-curved and anisotropic ultralocal state in which all spatialgradient contributions to the equations of motion decrease as an exponential in time to negligible values. This is followed by a second stage in which the geometry converges to a homogeneous, spatially flat and isotropic spacetime. In particular, the decay appears to follow the same history whether the entire spacetime or only parts of it are smoothed by the end of slow contraction.

Astronomy & Astrophysics↗
DOE OSTI2023

Observational constraints on warm natural inflation

Warm natural inflation is studied for the case of the original cosine potential. The radiation bath during inflation induces a dissipation (friction) rate in the equation of motion for the inflation field, which can potentially reduce the field excursion needed for an observationally viable period of inflation. Here we examine if the dissipation thus provides a mechanism to avoid the large decay constant f ≳ M pl of cold cosine natural inflation. Whereas temperature independent dissipation has previously been shown to alleviate the need for a trans-Planckian decay constant f, we illustrate here the difficulties of accommodating a significantly sub-Planckian decay constant (f < 10 -1 M pl ) for the case of a temperature dependent dissipation rate in the form Γ ∝ T c , with c = {1,3}. Such dissipation rates represent physically well-motivated constructions in the literature. For each model, we map its location in the r-n s plane and compare with cosmic microwave background (CMB) data. For c = 1 (c = 3), we find that agreement with CMB data requires that dissipation be in the weak (moderate) regime and that the minimum allowed value of the decay constant in the potential is f min = 0.3 (0.8) M pl , respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Determination of the characteristic magnetic pre-sheath length at divertor surfaces using micro-engineered targets on DiMES at DIII-D

The magnetic pre-sheath (MPS) width, L MPS , is a critical parameter to define the sheath potential, which controls the ion trajectory of low-Z species (D, T, He, and C), as well as the prompt re-deposition of high-Z species. In this work, to determine L MPS , we fabricated micro-trenches (30×30×4 µm) via focused ion beam (FIB) milling on a silicon surface and exposed them to L-mode deuterium plasmas in DIII-D via the Divertor Material Evaluation System (DiMES) removable sample exposure probe. The areal distribution of impurity depositions, mainly consisting of carbon, was measured by energy-dispersive X-ray spectroscopy (EDS) to reveal the deuterium ion shadowing effect on the trench floors. The carbon deposition profiles showed that the erosion was maximized for the azimuthal direction of φ = -40° (referenced to the toroidal magnetic field direction) as well as the polar angle of θ = 80°. A Monte Carlo equation-of-motion model, based on a collisionless MPS, was used to calculate the azimuthal and polar deuterium ion angle distributions (IADs) for a range of L MPS = k × ρ i , where ρ i is the ion gyro radius and k = 0.5-4. Then, gross erosion profiles were calculated by a Monte Carlo micro-patterning and roughness (MPR) code for ion sputtering using as input the calculated azimuthal and polar IADs for each value of k . Good agreement with the experimental C deposition profiles was obtained for the case k = 2.5-3.5. This result is consistent with a previous kinetic modeling prediction of k ~ 3, as well as previous analytical investigations that predicted the L MPS to be several ion gyro radii. A validation of theoretical sheath models supports its applicability to ITER and pilot plant divertors to successfully predict plasma-materials interactions.

Erosion↗
DOE OSTI2020

Far From Equilibrium Hydrodynamics and the Beam Energy Scan

The existence of hydrodynamic attractors in rapidly expanding relativistic systems has shed light on the success of relativistic hydrodynamics in describing heavy-ion collisions at zero chemical potential. As the search for the QCD critical point continues, it is important to investigate how out of equilibrium effects influence the trajectories on the QCD phase diagram. In this proceedings, we study a Bjorken expanding hydrodynamic system based on DMNR equations of motion with initial out of equilibrium effects and finite chemical potential in a system with a critical point. We find that the initial conditions are not unique for a specific freeze-out point, but rather the system can evolve to the same final state freeze-out point with a wide range of initial baryon chemical potential, µ B . For the same initial energy density and baryon density, depending on how far out of equilibrium the system begins, the initial µ B can vary by Δµ B ~ 350 MeV. Furthermore, our results indicate that knowledge of the out-of-equilibrium effects in the initial state provide vital information that influences the search for the QCD critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Excited-state downfolding using ground-state formalisms

Downfolding coupled cluster (CC) techniques are powerful tools for reducing the dimensionality of many-body quantum problems. This work investigates how ground-state downfolding formalisms can target excited states using non-Aufbau reference determinants, paving the way for applications of quantum computing in excited-state chemistry. This study focuses on doubly excited states for which canonical equation-of-motion CC approaches struggle to describe unless one includes higher-than-double excitations. The downfolding technique results in state-specific effective Hamiltonians that, when diagonalized in their respective active spaces, provide ground- and excited-state total energies (and therefore excitation energies) comparable to high-level CC methods. The performance of this procedure is examined with doubly excited states of H 2 , Methylene, Formaldehyde, and Nitroxyl.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Newton versus the machine: solving the chaotic three-body problem using deep neural networks

ABSTRACT Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system’s chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.

Breen, Philip G.↗