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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An assessment of ocean alkalinity enhancement using aqueous hydroxides: kinetics, efficiency, and precipitation thresholds

Abstract. Ocean alkalinity enhancement (OAE) is a promising approach to marine carbon dioxide removal (mCDR) that leverages the large surface area and carbon storage capacity of the oceans to sequester atmospheric CO2 as dissolved bicarbonate (HCO3-). One OAE method involves the conversion of salt in seawater into aqueous alkalinity (NaOH), which is returned to the ocean. The resulting increase in seawater pH and alkalinity causes a shift in dissolved inorganic carbon (DIC) speciation toward carbonate and a decrease in the surface ocean pCO2. The shift in the pCO2 results in enhanced uptake of atmospheric CO2 by the seawater due to gas exchange. In this study, we systematically test the efficiency of CO2 uptake in seawater treated with NaOH at aquarium (15 L) and tank (6000 L) scales to establish operational boundaries for safety and efficiency in advance of scaling up to field experiments. CO2 equilibration occurred on the order of weeks to months, depending on circulation, air forcing, and air bubbling conditions within the test tanks. An increase of ∼0.7–0.9 mol DIC per mol added alkalinity (in the form of NaOH) was observed through analysis of seawater bottle samples and pH sensor data, consistent with the value expected given the values of the carbonate system equilibrium calculations for the range of salinities and temperatures tested. Mineral precipitation occurred when the bulk seawater pH exceeded 10.0 and Ωaragonite exceeded 30.0. This precipitation was dominated by Mg(OH)2 over hours to 1 d before shifting to CaCO3,aragonite precipitation. These data, combined with models of the dilution and advection of alkaline plumes, will allow the estimation of the amount of carbon dioxide removal expected from OAE pilot studies. Future experiments should better approximate field conditions including sediment interactions, biological activity, ocean circulation, air–sea gas exchange rates, and mixing zone dynamics.

Ringham, Mallory C. (ORCID:0000000348020185)↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Uniformly decaying subspaces for error-mitigated quantum computation

Here, we present a general condition to obtain subspaces that decay uniformly in a system governed by the Lindblad master equation and use them to perform error-mitigated quantum computation. The expectation values of dynamics encoded in such subspaces are unbiased estimators of noise-free expectation values. In analogy to the decoherence free subspaces which are left invariant by the action of Lindblad operators, we show that the uniformly decaying subspaces are left invariant (up to orthogonal terms) by the action of the dissipative part of the Lindblad equation. We apply our theory to a system of qubits and qudits undergoing relaxation with varying decay rates and show that such subspaces can be used to eliminate bias up to first-order variations in the decay rates without requiring full knowledge of noise. Since such a bias cannot be corrected through standard symmetry verification, our method can improve error mitigation in dual-rail qubits and, given partial knowledge of noise, can perform better than probabilistic error cancellation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Geometric integration of classical spin dynamics via a mean-field Schrödinger equation

The Landau-Lifshitz equation describes the time evolution of magnetic dipoles and can be derived by taking the classical limit of a quantum mechanical spin Hamiltonian. To take this limit, one constrains the many-body quantum state to a tensor product of coherent states, thereby neglecting entanglement between sites. Expectation values of the quantum spin operators produce the usual classical spin dipoles. One may also consider expectation values of polynomials of the spin operators, leading to quadrupole and higher-order spin moments, which satisfy a dynamical equation of motion that generalizes the Landau-Lifshitz dynamics [Zhang and Batista, Phys. Rev. B 104, 104409 (2021)]. Here we reformulate the dynamics of these N 2 –1 generalized spin components as a mean-field Schrödinger equation on the N-dimensional coherent state. Furthermore, this viewpoint suggests efficient integration methods that respect the local symplectic structure of the classical spin dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Stirred-Reactor Coupon Analysis: An International Round Robin Study

The objective of this task was to determine the precision of the SRCA technique when used to determine the dilute condition corrosion rate. To this end, an interlaboratory round robin study was conducted per the instructions in ASTM Practice E691 to measure the precision with which the SRCA test method can be conducted. Twelve independent labs from eleven different institutions each evaluated four glass compositions in three different conditions. The ASTM procedures recommend at least 6 labs participate in a round robin testing the same 3 materials in the same conditions to determine precision. In this case, 12 labs each performed 12 independent tests. This was only possible thanks to the multi-glass testing capability of the SRCA test. A total of 108 duplicate pairs were used to calculate the repeatability of the tests, with the same glass tested in the same vessel, producing as identical conditions as possible for the replicates. These test results were quite tightly clustered, with a median difference from the average value of the pair of only 2.9%. Based on the calculations outlined in ASTM E177-20 and a measured standard deviation of 4.74%, the intralaboratory repeatability limit (r) was calculated to be within 13.3% of the expected value with 95% confidence level. The reproducibility limit (R) of the test was examined using all 277 data points from the round robin. Because of the differences in dissolution rates due to pH variability and intrinsically for the 12 conditions tested, the reproducibility limits for each condition and overall were calculated from the percent relative residual value for each test. The SRCA test is expected to be reproducible within 64% of the expected value with a 95% confidence level.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Demonstration of 151 Sm analysis using Trinity nuclear debris samples

Here, the long-lived fission product 151 Sm (T 1/2 ≈ 90 y) present in archived dissolved Trinity nuclear test debris samples was analyzed using liquid scintillation counting and the results were compared to the expected values based upon historical analysis results. The measured 151 Sm activities were found to agree with the expected values within experimental error. This demonstration validates the use of this technique in situations where the sample is too old for traditional shorter half-life radioactive fission product signatures such as 99 Mo (T 1/2 = 2.75 d), 97 Zr (T 1/2 = 16.7 h), 147 Nd (T 1/2 = 10.98 d), and 153 Sm (T 1/2 = 46.3 h), or too dilute for determining fissions through stable isotope ratio measurements.

151Sm↗
DOE OSTI2026

Quantum utility in simulating the real-time dynamics of the Fermi–Hubbard model using superconducting quantum computers

The Fermi–Hubbard model is a fundamental model in condensed matter physics that describes strongly correlated electrons. On the other hand, quantum computers are emerging as powerful tools for exploring the complex dynamics of these quantum many-body systems. In this work, we demonstrate the quantum simulation of the one-dimensional Fermi–Hubbard model using IBM's superconducting quantum computers, employing over 100 qubits. We introduce a first-order Trotterization scheme and extend it to an optimized second-order Trotterization for the time evolution in the Fermi–Hubbard model, specifically tailored for the limited qubit connectivity of quantum architectures, such as IBM's platforms. Notably, both Trotterization approaches are scalable and maintain a constant circuit depth at each Trotter step, regardless of the qubit count, enabling us to precisely investigate the relaxation dynamics in the Fermi–Hubbard model by measuring the expectation value of the Néel observable (staggered magnetization) for time-evolved quantum states. Lastly, our successful measurement of expectation values in such large-scale quantum many-body systems, especially at longer time scales with larger entanglement, highlights the quantum utility of superconducting quantum platforms over conventional classical approximation methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Measurement of the CP properties of Higgs boson interactions with $\tau$-leptons with the ATLAS detector

A study of the charge conjugation and parity ($CP$) properties of the interaction between the Higgs boson and $\tau$-leptons is presented. The study is based on a measurement of $CP$-sensitive angular observables defined by the visible decay products of $\tau$-leptons produced in Higgs boson decays. The analysis uses 139 fb -1 of proton–proton collision data recorded at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. Contributions from $\tau$-violating interactions between the Higgs boson and $\tau$-leptons are described by a single mixing angle parameter Φ $\tau$ in the generalised Yukawa interaction. Without constraining the $H$ $\rightarrow$ $\tau$$\tau$ signal strength to its expected value under the Standard Model hypothesis, the mixing angle Φ $\tau$ is measured to be 9° ± 16°, with an expected value of 0° ± 28° at the 68% confidence level. The pure $CP$-odd hypothesis is disfavoured at a level of 3.4 standard deviations. The results are compatible with the predictions for the Higgs boson in the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A technical analog of the cosmological constant problem and a solution thereof

The near vanishing of the cosmological constant is one of the most puzzling open problems in theoretical physics. We consider a system, the so-called framid, that features a technically similar problem. Its stress-energy tensor has a Lorentz-invariant expectation value on the ground state, yet there are no standard, symmetry-based selection rules enforcing this, since the ground state spontaneously breaks boosts. We verify the Lorentz invariance of the expectation value in question with explicit one-loop computations. These, however, yield the expected result only thanks to highly nontrivial cancellations, which are quite mysterious from the low-energy effective theory viewpoint.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On the holographic dual of a symmetry operator at finite temperature

Topological symmetry operators of holographic large 𝑁 CFT 𝐷 ’s are dual to dynamical branes in the gravity dual AdS 𝐷+1 . We use this correspondence to establish a dictionary between thermal expectation values of symmetry operators in the Euclidean CFT 𝐷 and the evaluation of gravitational saddles in the presence of a dynamical brane. Expectation values of 0-form symmetry operators in the CFT 𝐷 are then related to branes wrapped on volume minimizing cycles in the bulk, i.e., the Euclidean continuation of a black hole horizon. We illustrate with some representative examples, including gravity in AdS 3 , duality/triality defects in four-dimensional 𝒩 = 4 super Yang-Mills theory, and the dual of R-symmetry operators probing five-dimensional Bogomol’nyi–Prasad–Sommerfield black holes.

Anomalies↗
DOE OSTI2024

Characterizing non-Markovian and coherent errors in quantum simulation

Quantum simulation of many-body systems, particularly using ultracold atoms and trapped ions, presents a unique form of quantum control—it is a direct implementation of a multi-qubit gate generated by the Hamiltonian. As a consequence, it also faces a unique challenge in terms of benchmarking, because the well-established gate benchmarking techniques are unsuitable for this form of quantum control. Here we show that the symmetries of the target many-body Hamiltonian can be used not only to benchmark but to characterize experimental errors in the quantum simulation. We use our results to develop protocols to characterize these errors, which can be implemented using state-of-the-art technology. We consider two forms of errors: (i) unitary errors arising out of systematic errors in the applied Hamiltonian and (ii) canonical non-Markovian errors arising out of random shot-to-shot fluctuations in the applied Hamiltonian. We show that the dynamics of the expectation value of the target Hamiltonian itself, which is ideally constant in time, can be used to characterize these errors. In the presence of errors, the expectation value of the target Hamiltonian shows a characteristic thermalization dynamics, when it satisfies the operator thermalization hypothesis (OTH). That is, an oscillation in the short time followed by relaxation to a steady-state value in the long time limit. We show that while the steady-state value can be used to characterize the coherent errors, the amplitude of the oscillations can be used to estimate the non-Markovian errors. We prove a sandwich theorem to establish a linear relation between the amplitude of the oscillations and the magnitude of the non-Markovian errors. Moreover, by varying the initial state, we show that the steady state values can be used to completely construct the generator of the coherent errors. Using these results, we develop two experimental protocols to characterize the unitary errors based on these results, one of which requires single-qubit addressing and the other one doesn't. We also develop a protocol to characterize non-Markovian errors. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Variational neural network approach to QFT in the field basis

We present a variational neural network approach for solving quantum field theories in the field basis, focusing on the free Klein-Gordon model formulated in momentum space. While recent studies have explored neural-network-based variational methods for scalar field theory in position space, a systematic benchmark of the analytically solvable Klein-Gordon ground state—particularly in the momentum-space field basis—has been lacking. In this work, we represent the ground-state wavefunctional as a neural network defined on a discretized set of field configurations and train it by minimizing the Hamiltonian expectation value. This framework enables direct comparison to exact analytic results for a range of key observables, including the ground-state energy, two-point correlators, expectation value of the field, and the structure of the learned wavefunctional itself. Our results provide quantitative diagnostics of accuracy and establish a validated foundation for extending neural-network wavefunctional methods to interacting field theories and position-space formulations.

Klein-Gordon model↗
DOE OSTI2024

Path integrals, complex probabilities and the discrete Weyl representation

Abstract A discrete formulation of the real-time path integral as the expectation value of a functional of paths with respect to a complex probability on a sample space of discrete valued paths is explored. The formulation in terms of complex probabilities is motivated by a recent reinterpretation of the real-time path integral as the expectation value of a potential functional with respect to a complex probability distribution on cylinder sets of paths. The discrete formulation in this work is based on a discrete version of the Weyl algebra that can be applied to any observable with a finite number of outcomes. The origin of the complex probability in this work is the completeness relation. In the discrete formulation the complex probability exactly factors into products of conditional probabilities and exact unitarity is maintained at each level of approximation. The approximation of infinite dimensional quantum systems by discrete systems is discussed. The method is illustrated by applying it to scattering theory and quantum field theory. The implications of these applications for quantum computing is discussed.

Physics↗
DOE OSTI2022

The PDF perspective on the tracer-matter connection: Lagrangian bias and non-Poissonian shot noise

ABSTRACT We study the connection of matter density and its tracers from the probability density function (PDF) perspective. One aspect of this connection is the conditional expectation value 〈δtracer|δm〉 when averaging both tracer and matter density over some scale. We present a new way to incorporate a Lagrangian bias expansion of this expectation value into standard frameworks for modelling the PDF of density fluctuations and counts-in-cells statistics. Using N-body simulations and mock galaxy catalogues we confirm the accuracy of this expansion and compare it to the more commonly used Eulerian parametrization. For haloes hosting typical luminous red galaxies, the Lagrangian model provides a significantly better description of 〈δtracer|δm〉 at second order in perturbations. A second aspect of the matter-tracer connection is shot-noise, i.e. the scatter of tracer density around 〈δtracer|δm〉. It is well known that this noise can be significantly non-Poissonian and we validate the performance of a more general, two-parameter shot-noise model for different tracers and simulations. Both parts of our analysis are meant to pave the way for forthcoming applications to survey data.

Friedrich, Oliver↗
DOE OSTI2024

Demonstration of the rodeo algorithm on a quantum computer

The rodeo algorithm is an efficient algorithm for eigenstate preparation and eigenvalue estimation for any observable on a quantum computer. This makes it a promising tool for studying the spectrum and structure of atomic nuclei as well as other fields of quantum many-body physics. The only requirement is that the initial state has sufficient overlap probability with the desired eigenstate. While it is exponentially faster than well-known algorithms such as phase estimation and adiabatic evolution for eigenstate preparation, it has yet to be implemented on an actual quantum device. In this work, we apply the rodeo algorithm to determine the energy levels of a random one-qubit Hamiltonian, resulting in a relative error of 0.08% using mid-circuit measurements on the IBM Q device Casablanca. This surpasses the accuracy of directly-prepared eigenvector expectation values using the same quantum device. We take advantage of the high-accuracy energy determination and use the Hellmann-Feynman theorem to compute eigenvector expectation values for a different random one-qubit observable. For the Hellmann-Feynman calculations, we find a relative error of 0.7%. Here, we conclude by discussing possible future applications of the rodeo algorithm for multi-qubit Hamiltonians.

algorithm↗
DOE OSTI2024

Uranium Single Particle Analysis for Simultaneous Fluorine and Uranium Isotopic Determinations via Laser-Induced Breakdown Spectroscopy/Laser Ablation–Multicollector–Inductively Coupled Plasma–Mass Spectrometry

Uranyl fluoride (UO 2 F 2 ) particles (<20 μm) were subjected to first-of-its-kind analysis via simultaneous laser-induced breakdown spectroscopy (LIBS) and laser ablation multi-collector inductively coupled plasma–mass spectrometry (LA–MC–ICP–MS). Briefly, a nanosecond pulsed high-energy laser was focused onto the sample (particle) surface. In a single laser pulse, the UO 2 F 2 particle was excited/ionized within the microplasma volume, and the emission of light was collected via fiber optics such that emission spectroscopy could be employed for the detection of uranium (U) and fluorine (F). The ablated particle was simultaneously transported into the MC–ICP–MS for high precision isotopic (i.e., 234 U, 235 U, and 238 U) analysis. This method, LIBS/LA–MC–ICP–MS was optimized and employed to rapidly measure 80+ UO 2 F 2 particles, which were subjected to different calcination processes, which results in varying degrees of F loss from the individual particles. In measuring the particles, the average F/U ratios for the populations treated at 100 and 500 °C were 2.78 ± 1.28 and 1.01 ± 0.50, respectively, confirming loss of F through the calcination process. The average 235 U/ 238 U on the particle populations for the 100 and 500 °C were 0.007262 (22) and 0.007231 (23), which was determined to be <0.2% from the expected value. The 234 U/ 238 U ratios on the same particles were 0.000053 (11) and 0.000050 (10) for the 100 and 500 °C, respectively, <10% from the expected value. Notably, each population was analyzed in under 5 min, demonstrating the truly rapid analysis technique presented here in this paper.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

A Risk-Averse Approach for Distribution Grid Expansion Planning

Recent episodes of natural disasters have challenged the resilience of power grids. Adequate distribution grid planning that properly captures the risk aversion of the utility system planner is a key factor to increase the flexibility of distribution networks to circumvent these events. In this paper, we propose a methodology to determine the optimal portfolio of investments in lines and storage devices in order to minimize a convex combination between expected value and CVaR of operational costs, including energy not served, while taking into account the multistage nature of the energy storage management within this context. While the expected value of energy not served has been traditionally employed to tackle routine failures, we also minimize the CVaR of energy not served to address high-impact, low-probability (HILP) events. We illustrate the performance of the proposed methodology with a 54-Bus system test case.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Real classical shadows

Efficiently learning expectation values of a quantum state using classical shadow tomography has become a fundamental task in quantum information theory. In a classical shadows protocol, one measures a state in a chosen basis $\mathcal{W}$ after it has evolved under a unitary transformation randomly sampled from a chosen distribution $\mathcal{U}$. In this work we study the case where $\mathcal{U}$ corresponds to either local or global orthogonal Clifford gates, and $\mathcal{W}$ consists of real-valued vectors. Our results show that for various situations of interest, this ‘real’ classical shadow protocol improves the sample complexity over the standard scheme based on general Clifford unitaries. For example, when one is interested in estimating the expectation values of arbitrary real-valued observables, global orthogonal Cliffords typically decrease the required number of samples by a factor of two. More dramatically, for k-local observables composed only of real-valued Pauli operators, sampling local orthogonal Cliffords leads to a reduction by an exponential-in-k factor in the sample complexity over local unitary Cliffords. Finally, we show that by measuring in a basis containing complex-valued vectors, orthogonal shadows can, in the limit of large system size, exactly reproduce the original unitary shadows protocol.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗