Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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 379 records · Page 21

Differentiable Predictions for Large Scale Structure with SHAMNet

In simulation-based models of the galaxy-halo connection, theoretical predictions for galaxy clustering and lensing are typically made based on Monte Carlo realizations of a mock universe. In this paper, we use Subhalo Abundance Matching (SHAM) as a toy model to introduce an alternative to stochastic predictions based on mock population, demonstrating how to make simulation-based predictions for clustering and lensing that are both exact and differentiable with respect to the parameters of the model. Conventional implementations of SHAM are based on iterative algorithms such as Richardson-Lucy deconvolution; here we use the JAX library for automatic differentiation to train SHAMNet, a neural network that accurately approximates the stellar-to-halo mass relation (SMHM) defined by abundance matching. In our approach to making differentiable predictions of large-scale structure, we map parameterized PDFs onto each simulated halo, and calculate gradients of summary statistics of the galaxy distribution by using autodiff to propagate the gradients of the SMHM through the statistical estimators used to measure one- and two-point functions. Our techniques are quite general, and we conclude with an overview of how they can be applied in tandem with more complex, higher-dimensional models, creating the capability to make differentiable predictions for the multi-wavelength universe of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria: A DFT cluster expansion study

Magnetic contributions have the potential to significantly influence predicted phase stability within alloy and mineral mixing phase diagrams, yet have been historically challenging to incorporate due to a significant increase to phase space sampling. Here, in this work, we employ a computational protocol that includes spin orientation as an additional configurational component within multi-component cluster expansions between magnetic and non-magnetic metal oxide alloys [calculated using density functional theory (DFT) and the generalized gradient approximation]. This approach was used to determine the effect of magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria from first principles. Two-component cluster expansions of the magnetic components of eskolaite and hematite were first performed showing the ability of this method to properly calculate their respective magnetic properties. Two-component cluster expansions were then performed for non-magnetic Al(III) and ferromagnetic Cr(III) and Fe(III), and phase diagrams were calculated for later comparison. Finally, a non-magnetic Al(III) and “up” and “down” magnetic configurations for anti-ferromagnetic Cr(III) and Fe(III) were performed. Magnetic contributions to the calculated phase diagram for the corundumeskolaite system were shown to be inconsequential but are vital for accurate determination of the corundum-hematite solvus.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Highly Integrated WBG Power Module for Next Generation Plug-In Vehicles

Current inverter technologies used in electric vehicles (EVs) are bulky, heavy, and expensive. The performance of the inverter is also limited by the traditional Si-based semiconductor and packaging technology. These constraints are responsible for extra weight and volume in the vehicle that reduces the driving range along with usable space within the vehicle. Moreover, operating temperature restrictions, and boundaries on the inverter location within the vehicle drastically reduces the inverter performance. Thus an improved inverter design is required for next generation electric vehicles. The objective of this project was to research, develop and demonstrate a highly integrated WBG power module targeting the next generation of plugin hybrid electric vehicles (PHEVs). This power module needed to be functionally and mechanically optimized for General Motors (GM) traction inverter architecture. General Motors was the prime contractor, led the project, and oversaw the design, build, and test of the integrated power module. Virginia Polytechnic Institute and State University (VT) performed modeling and simulation tasks as well as power module design support. Monolith Semiconductor, Inc provided SiC dies and consulted on devices for the inverter. Oak Ridge National Laboratory (ORNL) performed power module technology evaluations, design evaluations of interconnect options, and performed material analysis and cross-sectioning of die activities.

33 ADVANCED PROPULSION SYSTEMS↗

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗

Active Optical Control with Machine Learning: A Proof of Concept for the Vera C. Rubin Observatory

Abstract The Active Optics System of the Vera C. Rubin Observatory (Rubin) uses information provided by four wave front sensors to determine deviations between the reconstructed wave front and the ideal wave front. The observed deviations are used to adjust the control parameters of the optical system to maintain image quality across the 3.°5 field of view. The baseline approach from the project is to obtain amplitudes of the Zernike polynomials describing the distorted wave front from out-of-focus images collected by the wave front sensors. These Zernike amplitudes are related via an “influence matrix” to the control parameters necessary to correct the wave front. In this paper, we use deep-learning methods to extract the control parameters directly from the images captured by the wave front sensors. Our neural net model uses anti-aliasing pooling to boost performance, and a domain-specific loss function to aid learning and generalization. The accuracy of the control parameters derived from our model exceeds Rubin requirements even in the presence of full-moon background levels and mis-centering of reference stars. Although the training process is time consuming, model evaluation requires only a few milliseconds. This low latency should allow for the correction of the optical configuration during the readout and slew interval between successive exposures.

Astronomy & Astrophysics↗

The Stochastic Gause Predator-Prey model: Noise-induced extinctions and invariance

We consider the Gause predator-prey with general bounded or sub‑linear functional responses, – which includes those of Holling types Ⅰ–Ⅳ. – and multiplicative Gaussian noise. In contrast to previous studies, the prey in our model follows logistic dynamics while the predator's population is solely regulated by consumption of the prey. To ensure well-posedeness, we derive explicit Lyapunov‐type criteria ensuring global positivity and moment boundedness of solutions. We find conditions for noise‑induced extinctions, proving that stochasticity can drive either population to collapse even when the deterministic analogue predicts stable coexistence. In the case when the predator becomes extinct, we establish a limiting distribution for the predator's population. Last, for functional responses of Holling type Ⅰ, we provide sufficient conditions on the intensity of the noise for the existence and uniqueness of a stationary distribution.

Gause model↗

Semi-classical Kinetic Theory for Massive Spin-half Fermions with Leading-order Spin Effects

We consider the quantum kinetic-theory description for interacting massive spin-half fermions using the Wigner function formalism. We derive a general kinetic theory description assuming that the spin effects appear at the classical and quantum level. To track the effect of such different contributions we use the semi-classical expansion method to obtain the generalized dynamical equations including spin, analogous to classical Boltzmann equation. This approach can be used to obtain a collision kernel involving local as well as non-local collisions among the microscopic constituent of the system and eventually, a framework of spin hydrodynamics ensuring the conservation of the energy-momentum tensor and total angular momentum tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measuring the Neutron Spin Asymmetry A1n in the Valence Quark Region in Hall C at Jefferson Lab

The quest to understand how the nucleon spin is decomposed into its constituent quark and gluon spin and orbital angular momentum (OAM) components has been at the forefront of nuclear physics for decades. Due to the non-perturbative nature of Quantum Chromodynamics (QCD) - the theory describing how quarks and gluons bind together to form protons and neutrons - making absolute predictions of nucleon spin structure is generally difficult, especially as a function of its quark and gluon longitudinal momentum fraction x. Measurements involving nucleon spin structure serve as a sensitive test for QCD, including ab-initio lattice QCD calculations due to the advent of the quasi-PDF formalism, and various predictions that diverge at large-x. The neutron spin asymmetry An 1 at high?x is a key observable for probing nucleon spin structure. In the valence domain (x > 0.5), sea effects are expected to be negligible, and so the total nucleon spin is considered to be carried by the valence quarks. The valence region can therefore enable us to study the role of quark OAM and other non-perturbative effects of the strong force. An 1 was measured in the deep inelastic scattering region of 0.40 < x < 0.75 and 6 < Q2 < 10 GeV2 in Hall C at Jefferson Lab using a 10.4 GeV longitudinally polarized electron beam, upgraded polarized 3He target, and the High Momentum Spectrometer (HMS) and Super High Momentum Spectrometer (SHMS). E12- 06-110 provides the first precision data in the valence quark region above x = 0.60, and its preliminary results proved consistent with earlier data disqualifying a pQCD model that excluded quark OAM. Combined with previous world proton data, the ratio of the polarized-to-unpolarized up quark momentum distribution (?u + ?u)/(u + u) remained positive at large-x, and the down quark (?d + ?d)/(d + d) remained negative.

Rehfuss, Melanie↗

Pauli energy contribution to the nucleus-nucleus interaction

Background: The Pauli exclusion principle plays a crucial role as a building block of many-body quantal systems comprised of fermions. It also induces a “Pauli repulsion” in the interaction between di-nuclear systems. It has been shown in [Phys. Rev. C 95, 031601(R) (2017)] that the Pauli repulsion widens the nucleus-nucleus potential barrier, thus hindering sub-barrier fusion. Purpose: To investigate the proton and neutron contributions to the Pauli repulsion, both in the bare potential neglecting shape polarization and transfer between the reactants, as well as in the dynamical potential obtained by accounting for such dynamical rearrangements. Methods: As the basis of our study we utilize the Pauli kinetic energy (PKE) obtained by studying the nuclear localization function (NLF). Recently this approach has been generalized to incorporate all of the dynamical and time-odd terms present in the nuclear energy density functional. This approach is employed in the density constrained frozen Hartree-Fock (DCFHF) and in the density constrained time-dependent Hartree-Fock (DC- TDHF) microscopic methods. Results: The PKE spatial distribution shows that a repulsion occurs in the neck between the nuclei when they first touch. Inside the barrier, neutrons can contribute significantly more to the Pauli repulsion in neutron-rich systems. Dynamical effects tend to lower the Pauli repulsion near the barrier. Proton and neutron dynamical contributions to the PKE significantly differ inside the barrier for asymmetric collisions, which is interpreted as an effect of multinucleon transfer. Conclusions: The PKE is shown to make a significant contribution to nuclear interaction potentials. Protons and neutrons can play very different roles in both the bare potential and in the dynamical rearrangement. Further microscopic studies are required to better understand the role of transfer and to investigate the effect of pairing and deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system’s complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Chemistry↗

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗

Higher-order $\beta$-functions in the Standard Model and beyond

In this contribution we consider the recent computation of the gauge coupling $\beta$-function at four loops and the Yukawa matrix $\beta$-function at three loops in the most general, renormalizable and four-dimensional quantum field theory. Furthermore, we discuss ambiguities and divergences arising in Yukawa matrix $\beta$-functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unveiling the Role of Surface Ir-Oxo Species in O 2 Evolution at IrO 2 Electrocatalysts via Embedded Cluster Multireference Calculations

Understanding the mechanisms driving the oxygen evolution reaction (OER) on iridium oxide (IrO 2 )-based catalysts is essential to improving their performance and enabling an actual scale-up of water-splitting photoelectrochemical cells. The mechanistic pathways at IrO 2 interfaces have been extensively investigated computationally using density functional theory (DFT), which predicts a high-energy barrier for the last step of the OER of molecular oxygen detachment and release from the catalyst surface. Nevertheless, surface O 2 over- and under-binding results by standard generalized gradient approximation and hybrid density functionals, respectively, call for further analysis of this crucial step via multireference methods. Aiming at unveiling the nature of such a barrier, we hereby address the formation of O 2 from the most-stable IrO 2 (110) surface with both periodic DFT and an electrostatic embedded cluster approach at the n-electron valence-state perturbation theory. With this multireference approach, we find a value for the aforementioned energy barrier that is much closer to experimental indications than DFT ones. Here, an in-depth analysis of the involved molecular orbitals suggests that the origin of this barrier is related to the breaking of a π interaction between O 2 and Ir surface atom and to a significant additional O 2 interaction with adjacent electrophilic Ir-oxo species, which is present under experimental operating conditions. Besides shedding light on the mechanism of the OER on IrO 2 , these findings point out the importance of multireference methods for dissecting complex reactions at electrocatalytic interfaces and pave the route for further investigations with effective embedding approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Higher-order $β$-functions in the Standard Model and beyond

In this contribution we consider the recent computation of the gauge coupling \beta β -function at four loops and the Yukawa matrix \beta β -function at three loops in the most general, renormalizable and four-dimensional quantum field theory. Furthermore, we discuss ambiguities and divergences arising in Yukawa matrix \beta β -functions.

Herren, Florian↗

Tropical amplitudes for colored Lagrangians

Recently a new formulation for scattering amplitudes in Tr(Φ 3 ) theory has been given based on simple combinatorial ideas in the space of kinematic data. This allows all-loop integrated amplitudes to be expressed as “curve integrals” defined using tropical building blocks — the “headlight functions”. This paper shows how the formulation extends to the amplitudes of more general Lagrangians. We will present a number of different ways of introducing tropical “numerator functions” that allow us to describe general Lagrangian interactions. The simplest family of these “tropical numerators” computes the amplitudes of interesting Lagrangians with infinitely many interactions. We also describe methods for tropically formulating the amplitudes for general Lagrangians. One uses a variant of “Wick contraction” to glue together numerator factors for general interaction vertices. Another uses a natural characterization of polygons on surfaces to give a novel combinatorial description of all possible diagrams associated with arbitrary valence interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparing first-principles density functionals plus corrections for the lattice dynamics of YBa 2 Cu 3 O 6

The enigmatic mechanism underlying unconventional high-temperature superconductivity, especially the role of lattice dynamics, has remained a subject of debate. Theoretical insights have long been hindered due to the lack of an accurate first-principles description of the lattice dynamics of cuprates. Recently, using the r2SCAN meta-generalized gradient approximation (meta-GGA) functional, we have been able to achieve accurate phonon spectra of an insulating cuprate YBa 2 Cu 3 O 6 and discover significant magnetoelastic coupling in experimentally interesting Cu–O bond stretching optical modes [Ning et al., Phys. Rev. B 107, 045126 (2023)]. We extend this work by comparing Perdew–Burke–Ernzerhof and r2SCAN performances with corrections from the on-site Hubbard U and the D4 van der Waals (vdW) methods, aiming at further understanding on both the materials science side and the density functional side. We demonstrate the importance of vdW and self-interaction corrections for accurate first-principles YBa 2 Cu 3 O 6 lattice dynamics. Since r2SCAN by itself partially accounts for these effects, the good performance of r2SCAN is now more fully explained. Finally, in addition, the performances of the Tao–Mo series of meta-GGAs, which are constructed in a different way from the strongly constrained and appropriately normed (SCAN) meta-GGA and its revised version r2SCAN, are also compared and discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Composite anyons on a torus

An adiabatic approach put forward by Greiter and Wilczek interpolates between the integer quantum Hall effects of electrons and composite fermions by varying the statistical flux bound to electrons continuously from zero to an even integer number of flux quanta, such that the intermediate states represent anyons in an external magnetic field with the same “effective” integer filling factor. We consider such anyons on a torus, and construct representative wave functions for their ground as well as excited states. These wave functions involve higher Landau levels in general, but can be explicitly projected into the lowest Landau level for many parameters. We calculate the variational energy gap between the first excited state and ground state and find that it remains open as the statistical phase is varied. Lastly, we obtain from these wave functions, both analytically and numerically, various topological quantities, such as ground-state degeneracy, the Chern number, and the Hall viscosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗