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

Synthesis, Stability, and Magnetic Properties of Antiperovskite Co 3 PdN

Experimental synthesis and characterization of theoretically predicted compounds are important steps in the materials discovery pipeline. Here, we report on the synthesis of Co 3 PdN, which was recently predicted to be a stable magnetic antiperovskite. The Co 3 PdN thin films were grown by reactive sputtering and were confirmed to form in an antiperovskite crystal structure. The thermal stability of the compound is demonstrated up to 600 K by in situ X-ray diffraction, though the phase persists at slightly higher temperatures (700 K) in an air-free magnetometer. Both ab initio calculations and magnetization measurements find Co 3 PdN to be ferromagnetic with an experimentally determined Curie temperature of T C = 560 ± 5 K. The saturation magnetization of 1.2 μ B /Co found in the experiment is slightly lower than the 1.7 μ B /Co value expected by theory. A narrow magnetic hysteresis loop with a coercive field of 100 Oe at low temperature suggests that Co 3 PdN might be useful in electronic applications requiring fast switching of the magnetization vector. While prior prediction of Co 3 PdN showed a gapped electronic band structure for each spin channel, we show that this was due to incomplete sampling of Brillouin zone paths and that band crossings exist along R-X|M and X|M-R paths. The metallic nature of Co 3 PdN is further confirmed by temperature-dependent transport measurements, which also show a considerable anomalous Hall effect. Altogether, this work represents an appreciable step toward understanding the synthesis, structure, stability, and properties of a new magnetic material.

14 SOLAR ENERGY↗
DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

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

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗
DOE OSTI2020

Quadrupole topological photonic crystals

Quadrupole topological phases, exhibiting protected boundary states that are themselves topological insulators of lower dimensions, have recently been of great interest. Extensions of these ideas from current tight binding models to continuum theories for realistic materials require the identification of quantized invariants describing the bulk quadrupole order. Here we identify the analog of quadrupole order in Maxwell’s equations for a gyromagnetic photonic crystal (PhC) through a double-band-inversion process. The quadrupole moment is quantized by the simultaneous presence of crystalline symmetry and broken time-reversal symmetry, which is confirmed using three independent methods: analysis of symmetry eigenvalues, numerical calculations of the nested Wannier bands and the expectation value of the quadrupole operator. Furthermore, we reveal the boundary manifestations of quadrupole phases as quantized edge polarizations and fractional corner charges. The latter are the consequence of a filling anomaly of energy bands as first predicted in electronic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Lie algebraic theory of barren plateaus for deep parameterized quantum circuits

Variational quantum computing schemes train a loss function by sending an initial state through a parametrized quantum circuit, and measuring the expectation value of some operator. Despite their promise, the trainability of these algorithms is hindered by barren plateaus (BPs) induced by the expressiveness of the circuit, the entanglement of the input data, the locality of the observable, or the presence of noise. Up to this point, these sources of BPs have been regarded as independent. In this work, we present a general Lie algebraic theory that provides an exact expression for the variance of the loss function of sufficiently deep parametrized quantum circuits, even in the presence of certain noise models. Our results allow us to understand under one framework all aforementioned sources of BPs. This theoretical leap resolves a standing conjecture about a connection between loss concentration and the dimension of the Lie algebra of the circuit’s generators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Short-depth QAOA circuits and quantum annealing on higher-order ising models

Abstract We present a direct comparison between QAOA (Quantum Alternating Operator Ansatz), and QA (Quantum Annealing) on 127 qubit problem instances. QAOA with p = 1, 2 rounds is executed on the 127 qubit heavy-hex graph gate-model quantum computer ibm_washington, using on-device grid-searches for angle finding, and QA is executed on two Pegasus-chip D-Wave quantum annealers. The problems are random Ising models whose connectivity matches heavy-hex graphs and the Pegasus graph connectivity, and optionally include hardware-compatible cubic terms ( Z Z Z terms). The QAOA circuits are heavily optimized and of extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of the heavy-hex lattice. QAOA and QA are both compared against simulated annealing and the optimal solutions are computed exactly using CPLEX. The noiseless mean QAOA expectation values for p = 1, 2 are computed using classical light-cone based simulations. We find QA outperforms QAOA on the evaluated devices.

127 qubits↗
DOE OSTI2024

Group-theoretic error mitigation enabled by classical shadows and symmetries

Abstract Estimating expectation values is a key subroutine in quantum algorithms. Near-term implementations face two major challenges: a limited number of samples required to learn a large collection of observables, and the accumulation of errors in devices without quantum error correction. To address these challenges simultaneously, we develop a quantum error-mitigation strategy called symmetry-adjusted classical shadows , by adjusting classical-shadow tomography according to how symmetries are corrupted by device errors. As a concrete example, we highlight global U(1) symmetry, which manifests in fermions as particle number and in spins as total magnetization, and illustrate their group-theoretic unification with respective classical-shadow protocols. We establish rigorous sampling bounds under readout errors obeying minimal assumptions, and perform numerical experiments with a more comprehensive model of gate-level errors derived from existing quantum processors. Our results reveal symmetry-adjusted classical shadows as a low-cost strategy to mitigate errors from noisy quantum experiments in the ubiquitous presence of symmetry.

Zhao, Andrew (ORCID:0000000202990277)↗
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↗