Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Born 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.

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

At least 127 records · Page 7

Convective-storm environments in subtropical South America from high-frequency soundings during RELAMPAGO-CACTI

During the RELAMPAGO-CACTI field experiments in 2018-19, an unprecedented number of balloon-borne soundings were collected in Argentina. Radiosondes were launched from both fixed and mobile platforms, yielding 2712 soundings during the period 15 October 2018-30 April 2019. Approximately 20% of these soundings were collected by highly mobile platforms, strategically positioned for each intensive observing period, and launching approximately once per hour. The combination of fixed and mobile soundings capture both the overall conditions characterizing the RELAMPAGO-CACTI campaign, as well as the detailed evolution of environments supporting the initiation and upscale growth of deep convective storms, including some that produced hazardous hail and heavy rainfall. Episodes of frequent convection were characterized by sufficient quantities of moisture and instability for deep convection, along with deep-layer vertical wind shear supportive of organized or rotating storms. Eleven soundings showed most-unstable convective available potential energy (MUCAPE) exceeding 6000 J kg-1, comparable to the extreme instability observed in other parts of the world with intense deep convection. Parameters used to diagnose severe-storm potential showed that conditions were often favorable for supercells and severe hail, but not for tornadoes, primarily owing to insufficient low-level wind shear. High-frequency soundings also revealed the structure and evolution of the boundary layer leading up to convection initiation, convectively generated cold pools, the South American Low-Level Jet (SALLJ), and elevated nocturnal convection. This sounding dataset will enable improved understanding and prediction of convective storms and their surroundings in subtropical South America, as well as comparisons with other heavily studied regions such as the central United States that have not previously been possible.

Schumacher, Russ S.↗

Direct computation of magnetic surfaces in Boozer coordinates and coil optimization for quasisymmetry

We propose a new method to compute magnetic surfaces that are parametrized in Boozer coordinates for vacuum magnetic fields. We also propose a measure for quasisymmetry on the computed surfaces and use it to design coils that generate a magnetic field that is quasisymmetric on those surfaces. The rotational transform of the field and complexity measures for the coils are also controlled in the design problem. Using an adjoint approach, we are able to obtain analytic derivatives for this optimization problem, yielding an efficient gradient-based algorithm. Starting from an initial coil set that presents nested magnetic surfaces for a large fraction of the volume, our method converges rapidly to coil systems generating fields with excellent quasisymmetry and low particle losses. In particular for low complexity coils, we are able to significantly improve the performance compared with coils obtained from the standard two-stage approach, e.g. reduce losses of fusion-produced alpha particles born at half-radius from $17.7\,\%$ to $6.6\,\%$ . We also demonstrate 16-coil configurations with alpha loss ${<}1\,\%$ and neoclassical transport magnitude $\epsilon _{\text {eff}}^{3/2}$ less than approximately $5\times 10^{-9}$ .

Physics↗

Observation of e⁺e⁻ → pp$\bar{p}\bar{n}\pi^-$ + c. c.

Using data taken at 29 center-of-mass energies between 4.16 and 4.70 GeV with the BESIII detector at the Beijing Electron Positron Collider corresponding to a total integrated luminosity of approximately 18.8 fb -1 , the process e⁺e⁻ → pp$\overline{p}\overline{n}\pi^-$ + c. c. is observed for the first time with a statistical significance of 11.5σ. The average Born cross sections in the energy ranges of (4.160, 4.380) GeV, (4.400, 4.600) GeV and (4.610, 4.700) GeV are measured to be (21.5 ± 5.7 ± 1.2) fb, (46.3 ± 10.6 ± 2.5) fb and (59.0 ± 9.4 ± 3.2) fb, respectively, where the first uncertainties are statistical and the second are systematic. The line shapes of the $\overline{p}\overline{n}$ and ppπ⁻ invariant mass spectra are consistent with phase space distributions, indicating that no hexaquark or di-baryon state is observed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nonharmonic contributions to the high-temperature phonon thermodynamics of Cr

Here, phonon densities of states (DOSs) of body-centered cubic chromium were measured by time-of-flight inelastic neutron scattering at temperatures up to 1493 K. Density functional theory calculations with both quasiharmonic (QH) and anharmonic (AH) methods were performed at temperatures above the Néel temperature. Features in the phonon DOSs decrease in energy (soften) substantially with temperature. A Born–von Kármán analysis using fits to the experimental DOSs reveals a softening of almost 17% of the high-transverse phonon branch between 330 and 1493 K. The low-transverse branch changes by approximately half this amount. The AH calculations capture the observed behavior of the two transverse phonon branches, but the QH calculations give some inverted trends. Vibrational entropies from phonons and electrons are obtained, and their sum is in excellent agreement with the entropy of chromium obtained by calorimetry, indicating that above 330 K, no explicit temperature-dependent magnetic contributions are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Strong-field induced fragmentation and isomerization of toluene probed by ultrafast femtosecond electron diffraction and mass spectrometry

We investigate the fragmentation and isomerization of toluene molecules induced by strong-field ionization with a femtosecond near-infrared laser pulse. Momentum-resolved coincidence time-of-flight ion mass spectrometry is used to determine the relative yield of different ionic products and fragmentation channels as a function of laser intensity. Ultrafast electron diffraction is used to capture the structure of the ions formed on a picosecond time scale by comparing the diffraction signal with theoretical predictions. Through the combination of the two measurements and theory, we are able to determine the main fragmentation channels and to distinguish between ions with identical mass but different structures. In addition, our diffraction measurements show that the independent atom model, which is widely used to analyze electron diffraction patterns, is not a good approximation for diffraction from ions. Here, we show that the diffraction data is in very good agreement with ab initio scattering calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Complete lowest order radiative corrections in semi-inclusive scattering of polarized particles

For this work, the lowest order radiative corrections to the cross section and asymmetries measured in experiments on semi-inclusive deep inelastic scattering of polarized particles were calculated. Both exact and leading log expressions were presented and discussed for the total correction that include the contributions from the processes of i) real photon emission with semi-inclusive processes, ii) loop diagrams, and iii) real photon emission with exclusive processes. Radiative corrections to the Sivers and Collins asymmetries in $\pi^+$ electroproduction were studied numerically within the kinematical conditions of modern experimental environments at Jefferson Laboratory (JLab). The Wandzura-Wilczek approximation for the semi-inclusive structure functions and MAID2007 parameterization for the six amplitudes of exclusive processes were used in numeric analyses. The results show that i) radiative effects can generate a correction comparable to the size of Sivers and Collins asymmetries at the Born level, ii) there is good agreement between the exact and leading-order corrections, iii) external functions (that is, other than the Sivers and Collins functions in the respective asymmetries) can generate a contribution to the radiative correction up to 20%, and iv) there exists a strong dependence of the radiative correction on the models for semi-inclusive and exclusive structure functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Oxygen–Chlorine Chemisorption Scaling for Seawater Electrolysis on Transition Metals: The Role of Redox

To clarify what controls species oxidation selectivity in seawater electrolysis, density functional theory (DFT) is used to identify chemisorption enthalpy trends and scaling relations for the simplest relevant adsorbates (O, Cl, and H) on relevant surfaces of 3d transition metals, as well as Pd and Pt, in face-centered-cubic and, if different, their ground-state crystal structures. Approximations are tested for electron exchange-correlation (XC) and van der Waals interactions to assess their ability to reproduce experimental adsorption enthalpies of H and O on Pt(111). The vdW-uncorrected generalized gradient approximation to XC of Perdew, Burke, and Ernzerhof (PBE) agrees most closely with experiments. Using DFT-PBE thereafter, it is determined that the O chemisorption enthalpy on this wide range of transition-metal surfaces is proportional to the sum of first and second atomic ionization energies, akin to a Born–Haber cycle for a redox reaction, indicating that metal redox activity controls O chemisorption strength. Then it is shown that the O and Cl chemisorption enthalpies are strongly correlated, suggesting that the transition metals considered will oxidize unselectively water and Cl – . This strong correlation appears also for crystal reduction potentials of binary oxides and chlorides, indicating a fundamental challenge for future seawater electrode materials design.

08 HYDROGEN↗

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗

Phonon thermodynamics and elastic behavior of GaAs at high temperatures and pressures

The phonons of wurtzite and zinc blende GaAs were calculated at simultaneously elevated temperature and pressure, and elastic constants were calculated as functions of pressure. Pressure caused instabilities of shorter-wavelength transverse acoustic modes in both wurtzite and zinc blende GaAs, causing them to fall to zero at 18 and 20GPa, respectively. Here, the Born stability criteria, which depend on elastic constants and only long wavelength phonons, therefore overestimated the pressure needed to induce instability at 0 K. At elevated temperatures, explicit anharmonicity pushes the onset of instability to higher pressures in both wurtzite and zincblende GaAs. Phonon linewidth and densities of states data showed that the quasiharmonic approximation failed to account for temperature-induced phonon frequency shifts, and the quasiharmonic approximation became less reliable at elevated pressure. In general, the number of three-phonon processes increased with pressure, thereby increasing the temperature-driven broadening of phonon spectral lineshapes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson’s theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs–Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born–Oppenheimer molecular dynamics performed at a single temperature. In conclusion, this advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Engineered thermalization and cooling of quantum many-body systems

We develop a scheme for engineering genuine thermal states in analog quantum simulation platforms by coupling local degrees of freedom to driven, dissipative ancilla pseudospins. We demonstrate the scheme in a many-body quantum spin lattice simulation setting. A Born-Markov master equation describing the dynamics of the many-body system is developed, and we show that if the ancilla energies are periodically modulated, with a carefully chosen hierarchy of timescales, one can effectively thermalize the many-body system. Through analysis of the time-dependent dynamical generator, we determine the conditions under which the true thermal state is an approximate dynamical fixed point for general system Hamiltonians. Finally, we evaluate the thermalization protocol through numerical simulation and discuss prospects for implementation on current quantum simulation hardware.

36 MATERIALS SCIENCE↗

Electrical conductivity of a warm neutron star crust in magnetic fields: Neutron-drip regime

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at nonzero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve protoneutron stars born in supernova explosions, binary neutron star mergers, and accreting neutron stars. The temperature-density range studied covers the transition from a semidegenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (i) the screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (ii) the correlations of the ionic component in a one-component plasma, and (iii) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropy introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust, which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviations of relaxation time and components of the conductivity tensor from the average are below ≤25% except close to crust-core transition, where nonspherical nuclear structures are expected. Finally, our results can be used in dissipative magnetohydrodynamics simulations of warm compact stars.

Physics↗

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Simulation-based Method for Correcting Mode Coupling in CMB Angular Power Spectra

Modern cosmic microwave background (CMB) analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-C ℓ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, F ℓ , that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, ${J}_{{\ell }{\ell }^{\prime} }$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the Spider balloon-borne CMB experiment.

79 ASTRONOMY AND ASTROPHYSICS↗

Primordial black holes from bubble collisions during a first-order phase transition

We study the possibility of production of primordial black holes (PBHs) from bubble collisions during a first-order phase transition. While typical colliding bubbles are small and irrelevant for PBH production, we find that those that can produce PBHs must have a macroscopically thick fluid shell and have been born much before the typical nucleation time. Particularly large uncertainties arise from an exponential sensitivity of the nucleation rate on the required duration of bubble growth which depends on the details of the collisions and the evolution of the spacetime metric toward the end of the phase transition. We introduce a few parameters to be obtained from future numerical simulation to represent those unknowns, and estimate the PBH abundance in an Abelian Higgs benchmark model and show that it can be significant. We predict an approximately monochromatic PBH mass spectrum, and find regions in the parameter space where the PBHs can constitute entire dark matter or even overclose the Universe. Our result thus shows that models with a first-order phase transition can be constrained by overabundant PBHs or null results of other PBH searches. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Precise stellarator quasi-symmetry can be achieved with electromagnetic coils

Magnetic fields with quasi-symmetry are known to provide good confinement of charged particles and plasmas, but the extent to which quasi-symmetry can be achieved in practice has remained an open question. Recent work [M. Landreman and E. Paul, Phys. Rev. Lett. 128, 035001, 2022] reports the discovery of toroidal magnetic fields that are quasi-symmetric to orders-of-magnitude higher precision than previously known fields. We show that these fields can be accurately produced using electromagnetic coils of only moderate engineering complexity, that is, coils that have low curvature and that are sufficiently separated from each other. Our results demonstrate that these new quasi-symmetric fields are relevant for applications requiring the confinement of energetic charged particles for long time scales, such as nuclear fusion. The coils’ length plays an important role for how well the quasi-symmetric fields can be approximated. For the longest coil set considered and a mean field strength of 1 T, the departure from quasi-symmetry is of the order of Earth’s magnetic field. Additionally, we find that magnetic surfaces extend far outside the plasma boundary used by Landreman and Paul, providing confinement far from the core. Simulations confirm that the magnetic fields generated by the new coils confine particles with high kinetic energy substantially longer than previously known coil configurations. In particular, when scaled to a reactor, the best found configuration loses only 0.04% of energetic particles born at midradius when following guiding center trajectories for 200 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗