SEARCH · Engineering Papers

Results for “Adiabatic approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗

DOE OSTI2024

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Thermodynamically consistent versions of approximations used in modelling moist air

Abstract Some existing approaches to modelling the thermodynamics of moist air make approximations that break thermodynamic consistency , such that the resulting thermodynamics does not obey the first and second laws or has other inconsistencies. Recently, an approach to avoid such inconsistency has been suggested: the use of thermodynamic potentials in terms of their natural variables , from which all thermodynamic quantities and relationships (equations of state) are derived. In this article, we develop this approach for unapproximated moist‐air thermodynamics and two widely used approximations: the constant‐ approximation and the dry heat capacities approximation. The (consistent) constant‐ approximation is particularly attractive because it leads to, with the appropriate choice of thermodynamic variable, adiabatic dynamics that depend only on total mass and are independent of the breakdown between water forms. Additionally, a wide variety of material from different sources in the literature on thermodynamics in atmospheric modelling is brought together. It is hoped that this article provides a comprehensive reference for the use of thermodynamic potentials in atmospheric modelling, especially for the three systems considered here.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Dark radiation isocurvature from cosmological phase transitions

Abstract Cosmological first order phase transitions are typically associated with physics beyond the Standard Model, and thus of great theoretical and observational interest. Models of phase transitions where the energy is mostly converted to dark radiation can be constrained through limits on the dark radiation energy density (parameterized by ΔN eff ). However, the current constraint (ΔN eff < 0.3) assumes the perturbations are adiabatic. We point out that a broad class of non-thermal first order phase transitions that start during inflation but do not complete until after reheating leave a distinct imprint in the scalar field from bubble nucleation. Dark radiation inherits the perturbation from the scalar field when the phase transition completes, leading to large-scale isocurvature that would be observable in the CMB. We perform a detailed calculation of the isocurvature power spectrum and derive constraints on ΔN eff based on CMB+BAO data. For a reheating temperature ofT rh and a nucleation temperatureT * , the constraint is approximately ΔN eff ≲ 10 -5 (T * /T rh ) -4 , which can be much stronger than the adiabatic result. We also point out that since perturbations of dark radiation have a non-Gaussian origin, searches for non-Gaussianity in the CMB could place a stringent bound on ΔN eff as well.

Astronomy & Astrophysics↗

DOE OSTI2025

Counterdiabatic Driving with Performance Guarantees

Counterdiabatic (CD) driving has the potential to speed up adiabatic quantum state preparation by suppressing unwanted excitations. However, existing approaches either require intractable classical computations or are based on approximations that do not have performance guarantees. We propose and analyze a nonvariational, system-agnostic CD expansion method and analytically show that it converges exponentially quickly in the expansion order. In finite systems, the required resources scale inversely with the spectral gap, which we argue is asymptotically optimal. To extend our method to the thermodynamic limit and suppress errors stemming from high-frequency transitions, we leverage finite-time adiabatic protocols. In particular, we show that a time determined by the quantum speed limit is sufficient to prepare the desired ground state, without the need to optimize the adiabatic trajectory. Numerical tests of our method on the quantum Ising chain show that our method can outperform state-of-the-art variational CD approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Multi-seasonal measurements of the ground-level atmospheric ice-nucleating particle abundance on the North Slope of Alaska

Atmospheric ice-nucleating particles (INPs) are an important subset of aerosol particles that are responsible for the heterogeneous formation of ice crystals. INPs modulate the arctic cloud phase (liquid vs. ice), resulting in implications for radiative feedbacks. The number of arctic INP studies investigating specific INP episodes or sources increased recently. However, existing studies are based on short-duration field data, and long-term datasets are lacking. Continuous, long-term measurements are key to determining the abundance and variability of ambient arctic INPs and constraining aerosol–cloud interactions, e.g., to verify and/or improve simulations of mixed-phase clouds. Here, we present a new long-duration INP dataset from the Arctic: 2 years of predominantly immersion-mode INP concentrations (n INP ) measured continuously at the National Oceanic and Atmospheric Administration's Barrow Atmospheric Baseline Observatory (BRW) on the North Slope of Alaska. A portable ice nucleation experiment chamber (PINE-03), which simulates adiabatic expansion cooling, was used to directly measure the ground-level INP abundance with an approximately 12 min time resolution from October 2021 to December 2023. We document PINE-03 n INP measurements as well as estimated ice nucleation active surface site density (n s ) over a wide range of heterogeneous freezing temperatures from −16 to −31 °C from which we introduce new season-specific parameterizations suitable for modeling mixed-phase clouds. Collocated aerosol and meteorological data were analyzed to assess the correlation between ambient n INP , air mass origin region, and meteorological variability. Our findings suggest (1) very high freezing efficiency of INPs across the measured temperatures (n s ≈ 2×10 8 –10 10 m −2 for −16 to −31 °C), which is a factor of 10–1000 times greater efficiency as compared to that found in the previous mid-latitude INP measurements in fall using the same instrument; (2) surprisingly high n INP (≥ 1 L −1 at −25 °C) for the examined temperatures throughout the year that were not measured by PINE-03 at other sites; and (3) high n INP in spring, possibly related to arctic haze episodes. Relatively low concentrations of aerosol surface area and contrasting high-INP concentrations at BRW relative to mid-latitude sites are the possible reasons for the observed high freezing efficiency.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Dual-RF Phase Control for Adiabatic Energy Ramping in the EIC Rapid Cycling Synchrotron

The Electron-Ion Collider (EIC) Rapid Cycling Synchrotron (RCS) accelerates electrons from 750 MeV to 10 GeV over approximately 64.7 ms. While the baseline design simultaneously ramps both RF voltage and phase, this study evaluates an alternative where RF cavities are divided into two groups, each operating at a fixed 6.5 MV, with the effective voltage synthesized through phase control. In the symmetric scheme (Option A), both groups vary symmetrically about the synchronous phase to provide exact matching of energy gain and synchrotron tune throughout the ramp. Conversely, the asymmetric scheme (Option B) pins one group at the final operating phase while the other varies to supply the required energy gain. The feasibility of Option B is determined by the adiabaticity of the resulting synchrotron tune evolution, with ramp profiles, phase trajectories, and adiabaticity metrics analyzed for both implementations.

Gamage, B. [Thomas Jefferson National Accelerator ↗

DOE OSTI2023

General relativistic pulsations of ultra-massive ZZ Ceti stars

ABSTRACT Ultra-massive white dwarf stars are currently being discovered at a considerable rate, thanks to surveys such as the Gaia space mission. These dense and compact stellar remnants likely play a major role in Type Ia supernova explosions. It is possible to probe the interiors of ultra-massive white dwarfs through asteroseismology. In the case of the most massive white dwarfs, general relativity could affect their structure and pulsations substantially. In this work, we present results of relativistic pulsation calculations employing relativistic ultra-massive ONe-core white dwarf models with hydrogen-rich atmospheres and masses ranging from 1.29 to $1.369 \ \mathrm{M}_{\odot }$ with the aim of assessing the impact of general relativity on the adiabatic gravity (g)-mode period spectrum of very high mass ZZ Ceti stars. Employing the relativistic Cowling approximation for the pulsation analysis, we find that the critical buoyancy (Brunt–Väisälä) and acoustic (Lamb) frequencies are larger for the relativistic case, compared to the Newtonian case, due to the relativistic white dwarf models having smaller radii and higher gravities for a fixed stellar mass. In addition, the g-mode periods are shorter in the relativistic case than those in the Newtonian computations, with relative differences of up to ∼$50$ per cent for the highest mass models ($1.369 \ \mathrm{M}_{\odot }$) and for effective temperatures typical of the ZZ Ceti instability strip. Hence, the effects of general relativity on the structure, evolution, and pulsations of white dwarfs with masses larger than ∼$1.29 \ \mathrm{M}_{\odot }$ cannot be ignored in the asteroseismological analysis of ultra-massive ZZ Ceti stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Electron Drift Kinetics in Space and in the Lab [Slides]

In a magnetized plasma, dynamics are fundamentally split between the directions parallel and perpendicular to the magnetic field. Generally, gyration about the magnetic field will occur on smaller spatial and time scales than the bulk dynamics of the system. From Hamiltonian mechanics, a periodic motion on a time scale much shorter than the other system dynamics leads to an adiabatic invariant. The adiabatic invariant for this motion is the magnetic moment $\mu=\frac{mv⟂^2}{2B}$, which will be approximately conserved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. Here, these methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Iterative quantum optimization of spin glass problems with rapidly oscillating transverse fields

In this work, we introduce a new iterative quantum algorithm, called Iterative Symphonic Tunneling for Satisfiability problems (IST-SAT), which solves quantum spin glass optimization problems using high-frequency oscillating transverse fields. IST-SAT operates as a sequence of iterations, in which bitstrings returned from one iteration are used to set spin-dependent phases in oscillating transverse fields in the next iteration. Over several iterations, the novel mechanism of the algorithm steers the system toward the problem ground state. We benchmark IST-SAT on sets of hard MAX-3-XORSAT problem instances with exact state vector simulation, and report polynomial speedups over Trotterized adiabatic quantum computation and the best known semi-greedy classical algorithm. When IST-SAT is seeded with a sufficiently good initial approximation, the algorithm converges to exact solution(s) in a polynomial number of iterations. Our numerical results identify a critical Hamming radius, or quality of initial approximation, where the time-to-solution crosses from exponential to polynomial scaling in problem size. This work proposes IST-SAT a new quantum algorithm, which improves upon solutions obtained from initial classical or quantum optimization algorithms. The steering mechanism we introduce through IST-SAT presents a new path toward achieving quantum advantage in optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Efficient state preparation for the Schwinger model with a theta term

We present a comparison of different quantum state preparation algorithms and their overall efficiency for the Schwinger model with a theta term. While adiabatic state preparation is proved to be effective, in practice it leads to large gate counts to prepare the ground state. The quantum approximate optimization algorithm (QAOA) provides excellent results while keeping the counts small by design, at the cost of an expensive classical minimization process. We introduce a “blocked” modification of the Schwinger Hamiltonian to be used in the QAOA that further decreases the length of the algorithms as the size of the problem is increased. The rodeo algorithm (RA) provides a powerful tool to efficiently prepare any eigenstate of the Hamiltonian, as long as its overlap with the initial guess is large enough. We obtain the best results when combining the blocked QAOA ansatz and the RA, as this provides an excellent initial state with a relatively short algorithm without the need to perform any classical steps for large problem sizes. Published by the American Physical Society 2025

Bazavov, Alexei (ORCID:0000000321411901)↗

DOE OSTI2025

ExINP NSA Ice Nucleating Particle Concentrations

This data set comprises cumulative ambient ice-nucleating particle (INP) concentrations measured at the National Oceanic and Atmospheric Administration's (NOAA’s) Barrow Atmospheric Baseline Observatory (71.3230° N, 156.6114° W, “BRW” hereafter), next to the Atmospheric Radiation Measurement (ARM) North Slope of Alaska (NSA) site and ~ 6 km northeast of the town of Utqiaġvik. Our INP abundance data were generated using a combination of online instrument, the Portable Ice Nucleation Experiment chamber ver. 3 (PINE-03), and an offline cold stage, the West Texas Cryogenic Refrigerator Applied to Freezing Test system (WT-CRAFT). Our online INP data are all from the Examining the Ice-Nucleating Particles from NSA (ExINP-NSA) campaign conducted from October 19, 2021 to May 24, 2024. The offline INP concentration analysis was performed at West Texas A&M University for aerosol particle samples collected on polycarbonate filters (with 0.2-micron diameter pores). The PINE-03 measurements, as well as sampling activities for offline INP measurements, were conducted using the BRW site. An inset laminar sampling stack was mounted to the instrument platform, allowing PINE-03 to intake particle-laden air. For most of the campaign period, the semi-autonomous PINE-03 chamber was remotely controlled from West Texas A&M University using the LabView interface through the BeyondTrust remote-access console. PINE-03 was set to conduct an immersion freezing expansion experiment (i.e., simulated adiabatic cooling along with RHw at or above 100%). Except during the scheduled maintenance periods, the time resolution of each expansion experiment was approximately 12 minutes. PINE-03 continuously measured INP concentrations during the entire campaign without any substantial breaks. For most of the campaign period, PINE scanned its set-point vessel air temperatures from -14 °C to -31 °C and back to -14 °C about every 120 minutes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

ExINP NSA Ice Nucleating Particle Concentrations

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Nonadiabatic strong-field photoionization revisited: Detailed analysis

We develop a theoretical framework to interpret “bump-like” enhancements in the energy distribution of electrons emitted from cesium atoms subjected to intense infrared laser fields, as observed in a recent experiment. These emissions, detected at energies approaching twice the ponderomotive energy (2⁢𝑈 p ), indicate a mechanism where electrons are nonadiabatically ejected with substantial inward velocity following ionization. The bump-like features arise from returning trajectories that scatter off the residual ion. We interpret these experimental results classically using a periodically driven version of the Bohr model and quantum mechanically through Keldysh's adiabaticity parameter (𝛾 K ) and Reiss's continuum-state intensity parameter 𝒵 in a Coulomb-corrected version of the strong-field approximation (SFA) that incorporates “complex” electronic trajectories. Notably, our theoretical analysis shows that these emissions are independent of the specific short-range form of the atomic potential.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Spin-crossover complexes: Self-interaction correction vs density correction

Complexes containing a transition metal atom with a 3d 4 –3d 7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi–Löwdin self-interaction correction. Considering that Hartree–Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree–Fock method and the Perdew–Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. Here, we find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Investigation of plasmon relaxation mechanisms using nonadiabatic molecular dynamics

Hot carriers generated from the decay of plasmon excitation can be harvested to drive a wide range of physical or chemical processes. However, their generation efficiency is limited by the concomitant phonon-induced relaxation processes by which the energy in excited carriers is transformed into heat. However, simulations of dynamics of nanoscale clusters are challenging due to the computational complexity involved. Here, in this paper, we adopt our newly developed Trajectory Surface Hopping (TSH) nonadiabatic molecular dynamics algorithm to simulate plasmon relaxation in Au 20 clusters, taking the atomistic details into account. The electronic properties are treated within the Linear Response Time-Dependent Tight-binding Density Functional Theory (LR-TDDFTB) framework. The relaxation of plasmon due to coupling to phonon modes in Au 20 beyond the Born–Oppenheimer approximation is described by the TSH algorithm. The numerically efficient LR-TDDFTB method allows us to address a dense manifold of excited states to ensure the inclusion of plasmon excitation. Starting from the photoexcited plasmon states in Au 20 cluster, we find that the time constant for relaxation from plasmon excited states to the lowest excited states is about 2.7 ps, mainly resulting from a stepwise decay process caused by low-frequency phonons of the Au 20 cluster. Furthermore, our simulations show that the lifetime of the phonon-induced plasmon dephasing process is ~10.4 fs and that such a swift process can be attributed to the strong nonadiabatic effect in small clusters. Our simulations demonstrate a detailed description of the dynamic processes in nanoclusters, including plasmon excitation, hot carrier generation from plasmon excitation dephasing, and the subsequent phonon-induced relaxation process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗