Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Momentum Transfer from the DART Mission Kinetic Impact on Asteroid Dimorphos

The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC on 26 September 2022 as a planetary defence test. DART was the first hypervelocity impact experiment on an asteroid at size and velocity scales relevant to planetary defence, intended to validate kinetic impact as a means of asteroid deflection. Here we report a determination of the momentum transferred to an asteroid by kinetic impact. On the basis of the change in the binary orbit period, we find an instantaneous reduction in Dimorphos’s along-track orbital velocity component of 2.70 ± 0.10 mm s -1 , indicating enhanced momentum transfer due to recoil from ejecta streams produced by the impact. For a Dimorphos bulk density range of 1,500 to 3,300 kg m -3 , we find that the expected value of the momentum enhancement factor, β, ranges between 2.2 and 4.9, depending on the mass of Dimorphos. If Dimorphos and Didymos are assumed to have equal densities of 2,400 kg m -3 , β=3.61$^{+0.19}_{-0.25}$(1σ) . These β values indicate that substantially more momentum was transferred to Dimorphos from the escaping impact ejecta than was incident with DART. Therefore, the DART kinetic impact was highly effective in deflecting the asteroid Dimorphos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Data-driven quantum approximate optimization algorithm for power systems

Abstract Quantum technology provides a ground-breaking methodology to tackle challenging computational issues in power systems. It is especially promising for Distributed Energy Resources (DERs) dominant systems that have been widely developed to promote energy sustainability. In those systems, knowing the maximum sections of power and data delivery is essential for monitoring, operation, and control. However, high computational effort is required. By leveraging quantum resources, Quantum Approximate Optimization Algorithm (QAOA) provides a means to search for these sections efficiently. However, QAOA performance relies heavily on critical parameters, especially for weighted graphs. Here we present a data-driven QAOA, which transfers quasi-optimal parameters between weighted graphs based on the normalized graph density. We verify the strategy with 39,774 expectation value calculations. Without parameter optimization, our data-driven QAOA is comparable with the Goemans-Williamson algorithm. This work advances QAOA and pilots its practical application to power systems in noisy intermediate-scale quantum devices.

Jing, Hang (ORCID:0000000256343489)↗
DOE OSTI2023

Elucidating the mechanism underlying the augmented capacity of MoO 2 as an anode material in Li-ion batteries

Transition-metal oxide anode materials have been observed to possess an intriguing surplus of capacity beyond the expected values based on conversion reaction. However, the mechanisms behind this phenomenon have remained contentious and elusive. Here, this study focuses on synthesized nanosized molybdenum dioxide and its electrochemical performance as an anode material for Li-ion batteries. Our findings reveal a substantial increase in capacity upon cycling, achieving approximately 1688 mA h g -1 , nearly double the theoretical capacity, after 700 cycles at a 1C rate. To elucidate the mechanisms underlying this augmented capacity, a comprehensive analysis employing in situ and ex situ X-ray diffraction, X-ray absorption spectroscopy, scanning electron microscopy, and transmission electron microscopy was conducted at various stages of the Li-ion cell cycling. Our results indicate that no conversion reaction occurs during the initial discharge phase, with Li 2 O and Mo remaining undetected. Instead, Li 0.98 MoO 2 is generated upon lithiation. Further materials characterization employing electron energy loss spectroscopy and energy-dispersive X-ray spectroscopy on the cycled electrode suggests the potential formation of a metallic Li-rich layer at the interface of the Li-ion intercalated phase subsequent to the formation of Li 0.98 MoO 2 , contributing to the surplus Li storage. Moreover, electrochemical impedance spectroscopy coupled with ex situ SEM and TEM analyses reveals that alterations in particle size and morphology, along with changes in the solid electrolyte interphase (SEI) resistance, are instrumental in the capacity variation observed upon cycling.

25 ENERGY STORAGE↗
DOE OSTI2022

Central star formation in double-peak, gas-rich radio galaxies

The respective contributions of gas accretion, galaxy interactions, and mergers to the mass assembly of galaxies, as well as the evolution of their molecular gas and star-formation activity are still not fully understood. In a recent work, a large sample of double-peak (DP) emission-line galaxies have been identified from the SDSS. While the two peaks could represent two kinematic components, they may be linked to the large bulges that their host galaxies tend to have. Star-forming DP galaxies display a central star-formation enhancement and have been discussed as compatible with a sequence of recent minor mergers. In order to probe merger-induced star-formation mechanisms, we conducted observations of the molecular-gas content of 35 star-forming DP galaxies in the upper part of the main sequence (MS) of star formation (SF) with the IRAM 30 m telescope. Including similar galaxies 0.3 dex above the MS and with existing molecular-gas observations from the literature, we finally obtained a sample of 52 such galaxies. We succeeded in fitting the same kinematic parameters to the optical ionised- and molecular-gas emission lines for ten (19%) galaxies. We find a central star-formation enhancement resulting most likely from a galaxy merger or galaxy interaction, which is indicated by an excess of gas extinction found in the centre. This SF is traced by radio continuum emissions at 150 MHz, 1.4 GHz, and 3 GHz, all three of which are linearly correlated in log with the CO luminosity with the same slope. The 52 DP galaxies are found to have a significantly larger amount of molecular gas and longer depletion times, and hence a lower star-formation efficiency, than the expected values at their distance of the MS. The large bulges in these galaxies might be stabilising the gas, hence reducing the SF efficiency. This is consistent with a scenario of minor mergers increasing the mass of bulges and driving gas to the centre. We also excluded the inwards-directed gas migration and central star-formation enhancement as the origin of a bar morphology. Hence, these 52 DP galaxies could be the result of recent minor mergers that funnelled molecular gas towards their centre, triggering SF, but with moderate efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dimension-free path-integral molecular dynamics without preconditioning

Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics—including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods—is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann–Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of “dimension-free” numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi–Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

Korol, Roman (ORCID:0000000193076351)↗
DOE OSTI2020

A light weight regularization for wave function parameter gradients in quantum Monte Carlo

The parameter derivative of the expectation value of the energy, ∂E/∂p, is a key ingredient in variational Monte Carlo (VMC) wave function optimization methods. In some cases, a naïve estimate of this derivative suffers from an infinite variance, which inhibits the efficiency of optimization methods that rely on a stable estimate of the derivative. In this work, we derive a simple regularization of the naïve estimator, which is trivial to implement in existing VMC codes, has finite variance, and a negligible bias, which can be extrapolated to zero bias with no extra cost. We use this estimator to construct an unbiased, finite variance estimation of ∂E/∂p for a multi-Slater–Jastrow trial wave function on the LiH molecule and in the optimization of a multi-Slater–Jastrow trial wave function on the CuO molecule. This regularized estimator is a simple and efficient estimator of ∂E/∂p for VMC optimization techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be used to infer species separation, thermal decoupling, flows, or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to the underlying dependence on the shape and time evolution of the temperature and density profiles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second-order analytical expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second-order expression shows significant improvement over lower-order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling, both first and second-order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with a <10% error, the second-order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

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

Reflectivity and spectrum of relativistic flying plasma mirrors

Flying plasma mirrors induced by intense lasers via laser wake field acceleration scheme have been proposed as a promising way to generate few-cycle EUV or x-ray lasers. In addition, if such a relativistic plasma mirror can accelerate, then it would serve as an analog black hole to investigate the information loss paradox associated with the black hole Hawking evaporation. Among these applications, the reflectivity, which is usually frequency-dependent, would affect the outgoing photon spectrum and, therefore, impact on the analysis of the physics under investigation. In this paper, these two issues are investigated analytically and numerically with one-dimensional particle-in-cell simulations. Based on our simulation results, we propose a new model that provides a better estimate of the reflectivity than those studied previously. Moreover, we found that the peak frequency of the reflected spectrum of a Gaussian incident wave deviates from the expected value, 4γ 2 ω, due to the dependence of reflectivity on the frequency of the incident wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tempering stochastic density functional theory

Here, we introduce a tempering approach with stochastic density functional theory (sDFT), labeled t-sDFT, which reduces the statistical errors in the estimates of observable expectation values. This is achieved by rewriting the electronic density as a sum of a "warm" component complemented by "colder" correction(s). Since the warm component is larger in magnitude but faster to evaluate, we use many more stochastic orbitals for its evaluation than for the smaller-sized colder correction(s). This results in a significant reduction in the statistical fluctuations and systematic deviation compared to sDFT for the same computational effort. We demonstrate the method's performance on large hydrogen-passivated silicon nanocrystals, finding a reduction in the systematic deviation in the energy by more than an order of magnitude, while the systematic deviation in the forces is also quenched. Similarly, the statistical fluctuations are reduced by factors of ≈4-5 for the total energy and ≈1.5-2 for the forces on the atoms. Since the embedding in t-sDFT is fully stochastic, it is possible to combine t-sDFT with other variants of sDFT such as energy-window sDFT and embedded-fragmented sDFT.

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

Simple new methods for deducing lifetimes in recoil distance Doppler-shift measurements

In this work, new approaches for lifetime determination using data from recoil distance Doppler-shift experiments are presented based on the fundamental properties of the functions describing the time evolution of the population of excited nuclear states. To some extent, one of them represents a contraction of the well-known Differential decay-curve method (DDCM) by using the most reliable data point [the maximum of the $n_i(t)$ function describing the population of level $\textit{i}$ in time] and a purely numerical procedure avoiding any fitting of decay curves. The combination with the standard DDCM analysis is promising for improving the reliability and the precision of the results for the lifetimes obtained. The novel part of the approach consists of using a chain of equations at the consecutive maxima of the ni(t) functions, which allow us to precisely determine the ratio of the lifetimes of two consecutive levels and, in the case where one of these lifetimes is known, to determine the unknown one. In addition, a simple integral derivation of the lifetime is presented involving the peak areas measured at different distances, and an application of the first moments (expectation values and centroids in time) of the $n_i(t)$ functions for determining lifetimes is also demonstrated to be useful.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Simultaneous determination of electrical conductivity and thickness of nonmagnetic metallic foils using two-sided multifrequency eddy current techniques

Here, in this paper, we describe a simple method for performing multifrequency eddy current characterization of free-standing uniform-thickness metallic foils using a forked inductive coil arrangement. The method involves measuring the mutual inductance between two coils when a foil is present between the coils, and when it is not present; the ratio of these mutual inductances is compared with an analytical solution, and foil conductivity, thickness, and sheet resistance are simultaneously estimated using numerical inversion and least-squares fitting. This method was used to characterize 34 non-ferrous metallic samples with thicknesses between 50 and 640 μm and with conductivities between 0.8 × 10 7 and 5.8 × 10 7 S/m. The estimated thicknesses from eddy current characterization agreed well with those measured using confocal optical techniques; the two approaches agreed to within 1 μm for samples that were thinner than 200 μm, and to within 0.5% for samples that had a thickness of 200 μm or greater. The estimated conductivities from eddy current characterization were in close agreement with expected values, given knowledge of the materials used. A particular strength of this approach is that the instrumentation needed is broadly available in research and development laboratories and the associated fixturing is easy to manufacture and assemble. A calibration procedure is described that can be used to reduce errors from geometric uncertainties. This calibration requires a sample that has only a known conductivity or thickness; both do not need to be known. The method described herein is likely extensible to conductivities and thickness well outside the ranges measured as part of this work.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Thermoradiative energy conversion in infrared interband cascade InAs/GaSb/AlInSb/GaSb type-II superlattice diodes

Thermoradiative energy conversion presents a means for the direct conversion of thermal energy through radiative transfer to a cold scene. However, much of the study of thermoradiative principles has been based on theory and simulations, with only sparse reports on the experimental demonstration of the concept. This work studies thermoradiative energy conversion in InAs/GaSb/AlInSb/GaSb type-II superlattice cascade devices. The devices exhibit a cutoff wavelength of 3.2 μm at 300 K, corresponding to a bandgap of 0.39 eV. Testing under a temperature-controlled chamber and scene demonstrates a maximum power density of 2.9 mW/m2 at a cell temperature of 121 °C. It is consistent with expected values for device operation limited by Shockley–Read–Hall non-radiative recombination. This result is a significant step in providing an experimental demonstration of thermoradiative energy conversion and a means to characterize cell performance, providing a foundation for further development to achieve practical values for power generation.

Islam, Md Toriqul (ORCID:0000000249871320)↗
DOE OSTI2026

Low-energy ion beam sputtering of polycrystalline diamond

Ion beam milling is a highly versatile and precise technique for the microfabrication of diamond films, enabling controlled material removal at both micrometer and sub-micrometer scales while accommodating complex geometries and surface features. Comprehensive sputtering yield data are essential for process optimization. In this work, we systematically measure the sputter yield of polycrystalline diamond bombarded with $\textrm{N}_2^+$, Ne + , Ar + , Kr + , and Xe + ions over an energy range of 300–900 eV and beam incidence angles of 0°−80°, providing data critical for precision machining applications. Results show that the sputter yield increases with both ion energy and the incidence angle and decreases with increasing ion mass. For $\textrm{N}_2^+$ ions, the sputter yield is largely insensitive to ion energy and exceeds values expected for pure physical sputtering. The concurrent formation of volatile C 2 N 2 molecules during $\textrm{N}_2^+$ ion bombardment suggests that the process is dominated by chemical sputtering. The absence of a peak in the dependence of the sputter yield on the beam incidence angle expected at oblique incidence angles is attributed to the effects of surface roughness.

Materials science↗
DOE OSTI2020

An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices

In this paper, we report a random sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for under- standing the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbor hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbor hopping and 2D lattices. In addition, constraints on real-space con gurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Chiral soft algebras for $\mathcal{N}$ = $2$ gauge theory

Some time ago, Seiberg and Witten solved for moduli spaces of vacua parameterized by scalar vacuum expectation values in $\mathcal{N}$ = $2$ gauge theories. More recently, new vacua associated to soft theorems and asymptotic symmetries have been found. Furthermore, this paper takes some first steps towards a complete picture of the infrared geometry of $\mathcal{N}$ = $2$ gauge theory incorporating both of these infrared structures.

N = 2 gauge theory↗
DOE OSTI2020

Exact sum rules with approximate ground states

Electromagnetic and weak transitions tell us a great deal about the structure of atomic nuclei. Yet modeling transitions can be difficult: it is often easier to compute the ground state, if only as an approximation, than excited states. One alternative is through transition sum rules, in particular the non-energy-weighted and energy-weighted sum rules, which can be computed as expectation values of operators. We investigate by computing sum rules for a variety of nuclei, comparing the numerically exact full configuration-interaction shell model, as a reference, to Hartree-Fock, projected Hartree-Fock, and the nucleon pair approximation. Furthermore, these approximations yield reasonable agreement, which we explain by prior work on the systematics of transition moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Wire, hybrid, and laser-cut X-pinches as Talbot–Lau backlighters for electron density diagnostics

Talbot–Lau x-ray deflectometry (TXD) enables refraction-based imaging for high-energy-density physics experiments, and thus, it has been studied and developed with the goal of diagnosing plasmas relevant to inertial confinement and magnetic liner inertial fusion. X-pinches, known for reliably generating fast (~1 ns), small (~1 µm) x-ray sources, were driven on the compact current driver generator for ablation structure and implosion studies (~200 kA, 150 ns) as a potential backlighter source for TXD. Considering that different X-pinch configurations have characteristic advantages and drawbacks as x-ray generating loads, three distinct copper X-pinch configurations were studied: the wire X-pinch, the hybrid X-pinch, and the laser-cut X-pinch. Furthermore, the Cu K-shell emission from each configuration was characterized and analyzed regarding the specific backlighter requirements for an 8 keV TXD system: spatial and temporal resolution, number of sources, time of emission, spectrum, and reproducibility. Recommendations for future experimental improvements and applications are presented. The electron density of static objects was retrieved from Moiré images obtained through TXD. This allowed to calculate the mass density of static samples within 4% of the expected value for laser-cut X-pinches, which were found to be the optimal X-pinch configuration for TXD due to their high reproducibility, small source size (≤5 µm), short duration (~1 ns), and up to 106 W peak power near 8 keV photon energy. Plasma loads were imaged through TXD for the first-time using laser-cut X-pinch backlighting. Experimental images were compared with simulations from the x-ray wave-front propagation code, demonstrating that TXD can be a powerful x-ray refraction-based diagnostic for dense Z-pinch loads. Future plans for Talbot–Lau interferometry diagnostics in the pulsed-power environment are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗