Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Window”

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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Performance demonstration of vacuum microwave components critical for the operation of the ITER low-field side reflectometer

Final design studies in preparation for manufacturing have been performed for functional components of the vacuum portion of the ITER Low-Field Side Reflectometer (LFSR). These components consist of an antenna array, electron cyclotron heating (ECH) protection mirrors, phase calibration mirrors, and vacuum windows. Evaluation of these components was conducted at the LFSR test facility and DIII-D. The antenna array consists of six corrugated-waveguide antennas for simultaneous profile, fluctuation, and Doppler measurements. A diffraction grating, incorporated into the plasma-facing miter bend, provides protection of sensitive components from stray ECH at 170 GHz. For in situ phase calibration of the LFSR profile reflectometer, an embossed mirror is incorporated into the adjacent miter bend. Measurements of the radiated beam profile indicate that these components have a small, acceptable effect on mode conversion and beam quality. Baseline transmission characteristics of the dual-disk vacuum window are obtained and are used to guide ongoing developments. Preliminary simulations indicate that a surface-relief structure on the window surfaces can greatly improve transmission. The workability of real-time phase measurements was demonstrated on the DIII-D profile reflectometer. The new automated real-time analysis agrees well with the standard post-processing routine.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX. I. Sheath equivalent dielectric layer for modeling the RF sheath

The pulsed linear plasma device Prototype Material Plasma Exposure eXperiment (Proto-MPEX) uses a radio frequency (RF) helicon antenna with an aluminum nitride ceramic window for plasma production. The RF sheath created under the helicon antenna is sufficient to cause ion impact energies to be greater than the sputtering threshold of the AlN helicon window material and for impurities to be created. In this study, we investigate the RF sheath on the inner diameter of the helicon window and its impact on the impurity production rates in Proto-MPEX. Three models—a 3D COMSOL finite element RF model of the Proto-MPEX helicon region, a rectified DC sheath potential model, and the 3D Global Impurity TRansport code—are coupled together to study impurity production and transportation. This novel method of impurity generation and transport modeling spans length scales ranging from the sheath (millimeters or less) up to the full device (meters) and can be applied to other radio frequency sources and antennas in a wide range of plasma physics studies, including basic plasmas, low-temperature processing plasmas, plasma thrusters, and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurement of Photoneutrons from Depleted Uranium and Comparison Study Using MCNPX-PoliMi

High-energy photon interrogation is a nondestructive technique that is used to detect special nuclear materials and characterize nuclear waste. The development of such systems is complex and requires Monte Carlo simulations to optimize system performance. Monte Carlo simulations rely on various scattering, absorption, and photonuclear cross-section data. While the scattering and absorption cross-section data have been extensively studied and validated with experiments, the results obtained from photonuclear simulations are often found to underpredict measured results, indicating uncertainties in the photonuclear cross sections themselves. Thus, there is a need for new measured results that can be used to quantify underpredictions in simulations using photonuclear cross-section data. In the present work, we interrogated depleted uranium with a 9-MV electron linac and detected photoneutrons with trans-stilbene organic scintillators. The measurement of photoneutrons with organic scintillators is challenging due to the presence of the intense photon flux, which causes issues such as pulse pile-up, detector saturation, and poor signal-to-background ratio. To mitigate these challenges, we used iron and polyethylene shielding of varying thicknesses around the depleted uranium target and a neural network–based digital pulse processing algorithm to recover neutron and photon information from piled-up events. Our goal was to compare the measured photoneutron count rate with the simulated rate obtained using the MCNPX-PoliMi transport code. For a light output window of 0.28 to 2.67 MeVee (1.66- to 6.85-MeV proton recoil energy), we found that the simulated count rate obtained using the ENDF/B-VII photonuclear cross-section library underpredicts the measured rate by 32.8% 3.2%. Additionally, we compared the simulated and measured photoneutron light output distributions. For the least thicknesses of shielding, the simulation was found to underpredict measurements in the 0.70- to 2.67-MeVee light output window. For the greatest thicknesses of shielding, the simulation was found to underpredict the measurement across the entire light output window of 0.28 to 2.67 MeVee.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2013

THE RADIAL VELOCITY DETECTION OF EARTH-MASS PLANETS IN THE PRESENCE OF ACTIVITY NOISE: THE CASE OF {alpha} CENTAURI Bb

We present an analysis of the publicly available HARPS radial velocity (RV) measurements for {alpha} Cen B, a star hosting an Earth-mass planet candidate in a 3.24 day orbit. The goal is to devise robust ways of extracting low-amplitude RV signals of low-mass planets in the presence of activity noise. Two approaches were used to remove the stellar activity signal which dominates the RV variations: (1) Fourier component analysis (pre-whitening), and (2) local trend filtering (LTF) of the activity using short time windows of the data. The Fourier procedure results in a signal at P = 3.236 days and K = 0.42 m s{sup -1}, which is consistent with the presence of an Earth-mass planet, but the false alarm probability for this signal is rather high at a few percent. The LTF results in no significant detection of the planet signal, although it is possible to detect a marginal planet signal with this method using a different choice of time windows and fitting functions. However, even in this case the significance of the 3.24 day signal depends on the details of how a time window containing only 10% of the data is filtered. Both methods should have detected the presence of {alpha} Cen Bb at a higher significance than is actually seen. We also investigated the influence of random noise with a standard deviation comparable to the HARPS data and sampled in the same way. The distribution of the noise peaks in the period range 2.8-3.3 days has a maximum of Almost-Equal-To 3.2 days and amplitudes approximately one-half of the K-amplitude for the planet. The presence of the activity signal may boost the velocity amplitude of these signals to values comparable to the planet. It may be premature to attribute the 3.24 day RV variations to an Earth-mass planet. A better understanding of the noise characteristics in the RV data as well as more measurements with better sampling will be needed to confirm this exoplanet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A BROADBAND RADIO STUDY OF THE AVERAGE PROFILE AND GIANT PULSES FROM PSR B1821-24A

We present the results of a wideband (720–2400 MHz) study of PSR B1821–24A (J1824–2452A, M28A), an energetic millisecond pulsar (MSP) visible in radio, X-rays and γ-rays. In radio, the pulsar has a complex average profile that spans ≳85% of the spin period and exhibits strong evolution with observing frequency. For the first time we measure phase-resolved polarization properties and spectral indices of radio emission throughout almost all of the on-pulse window. We synthesize our findings with high-energy information to compare M28A to other known γ-ray MSPs and to speculate that M28A’s radio emission originates in multiple regions within its magnetosphere (i.e., both in the slot or outer gaps near the light cylinder and at lower altitudes above the polar cap). M28A is one of a handful of pulsars that are known to emit giant radio pulses (GRPs)—short, bright radio pulses of unknown nature. We report a drop in the linear polarization of the average profile in both windows of GRP generation and also a “W”-shaped absorption feature (resembling a double notch), partly overlapping with one of the GRP windows. The GRPs themselves have broadband spectra consisting of multiple patches with Δν/ν ∼ 0.07. Although our time resolution was not sufficient to resolve the GRP structure on the μs scale, we argue that GRPs from this pulsar most closely resemble the GRPs from the main pulse of the Crab pulsar, which consist of a series of narrowband nanoshots.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

3D vacuum magnetic field modelling of the ITER ELM control coil during standard operating scenarios

In-vessel, non-axisymmetric, control coils have proven to be an important option for mitigating and suppressing edge-localized modes (ELMs) in high performance operating regimes on a growing number of tokamaks. Additionally, an in-vessel non-axisymmetric ELM control coil is being considered in the ITER baseline design. For this paper, in preparing for the initial operation of this coil set, a comprehensive study was carried out to characterize the linear superposition of the 3D vacuum magnetic field, produced by the ELM coil, on a series of equilibria representing nine standard ITER operating scenarios. Here, the spatial phase angle of toroidally distributed currents, specified with a cosine waveform, in the upper and lower rows of the ITER ELM coil (IEC) set is varied in 2° steps while holding the current in the equatorial row of coils constant. The peak current in each of the three toroidal rows of window-frame coils making up the IEC is scanned between 5 kAt and 90 kAt in 5 kAt steps and the width of the edge region covered by overlapping vacuum field magnetic islands is calculated. This width is compared to a vacuum field ELM suppression correlation criterion found in DIII-D. A minimum coil current satisfying the DIII-D criterion, along with an associated set of phase angles, is identified for each ITER operating scenario. These currents range from 20 kAt to 75 kAt depending on the operating scenario being used and the toroidal mode number (n) of the cosine waveform. Comparisons between the scaling of the divertor footprint area in cases with n = 3 perturbation fields versus those with n = 4 show significant advantages when using n = 3. In addition, it is found that the DIII-D correlation criterion can be satisfied in the event that various combinations of individual IEC window-frame coils need to be turned off due to malfunctioning components located inside the vacuum vessel. Finally, details of these results for both the full set of 27 window-frame coils and various reduced sets, using either n = 3 and n = 4 perturbation fields, are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Regarding the optimization of O1-mode ECRH and the feasibility of EBW startup on NSTX-U

Recently published scenarios for fully non-inductive startup and operation on the National Spherical Torus eXperiment Upgrade (NSTX-U) [Menard J et al 2012 Nucl. Fusion 52 083015] show Electron Cyclotron Resonance Heating (ECRH) as an important component in preparing a target plasma for efficient High Harmonic Fast Wave and Neutral Beam heating. The modelling of the propagation and absorption of EC waves in the evolving plasma is required to define the most effective window of operation, and to optimize the launcher geometry for maximal heating and current drive during this window. Here in this paper, we extend a previous optimization of O1-mode ECRH on NSTX-U to account for the full time-dependent performance of the ECRH using simulations performed with TRANSP. We find that the evolution of the density profile has a prominent role in the optimization by defining the time window of operation, which in certain cases may be a more important metric to compare launcher performance than the average power absorption. This feature cannot be captured by analysis on static profiles, and should be accounted for when optimizing ECRH on any device that operates near the cutoff density. Additionally, the utility of the electron Bernstein wave (EBW) in driving current and generating closed flux surfaces in the early startup phase has been demonstrated on a number of devices. Using standalone GENRAY simulations, we find that efficient EBW current drive is possible on NSTX-U if the injection angle is shifted below the midplane and aimed towards the top half of the vacuum vessel. However, collisional damping of the EBW is projected to be significant, in some cases accounting for up to 97\% of the absorbed EBW power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

PISCES-RF: a liquid-cooled high-power steady-state helicon plasma device

Radio-frequency (RF) driven helicon plasma sources can produce relatively high-density plasmas (n > 1019 m -3 ) at relatively moderate powers (<2 kW) in argon. However, to produce similar high-density plasmas for fusion relevant gases such as hydrogen (H), deuterium (D) and helium (He), much higher RF powers are needed. For very high RF powers, thermal issues of the RF-transparent dielectric window, used in the RF source design, limit the plasma operation timescales. To mitigate this constraint, we have designed, built and tested a novel helicon plasma source assembly with a fully liquid-cooled RF-transparent window which allows steady state operations at high power (up to 20 kW) and successfully produces high-density plasma with both argon and H. Deionized (DI) water, flowing between two concentric dielectric RF windows, is used as the coolant. We show that a full azimuthal blanket of DI water does not prevent high-density plasma production. From calorimetry on the DI water, we measure the net heat removed by the coolant at steady state conditions. Using infra-red imaging, we calculate the constant plasma heat deposition and measure the final steady state temperature distribution patterns on the inner surface of the ceramic layer. The heat deposition pattern follows the helical shape of the antenna. We also show the consistency between the heat absorbed by the DI water, as measured by calorimetry, and the total heat due to the combined effect of the plasma heating and the absorbed RF. These results are being used to answer critical engineering questions for the 200 kW RF device materials plasma exposure experiment being designed at the Oak Ridge National Laboratory as a next generation plasma material interaction device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

GravitoMagnetic force in modified Newtonian dynamics

We introduce the Gauge Vector-Tensor (GVT) theory by extending the AQUAL's approach to the GravitoElectroMagnetism (GEM) approximation of gravity. GVT is a generally covariant theory of gravity composed of a pseudo Riemannian metric and two U(1) gauge connections that reproduces MOND in the limit of very weak gravitational fields while remains consistent with the Einstein-Hilbert gravity in the limit of strong and Newtonian gravitational fields. GVT also provides a simple framework to study the GEM approximation to gravity. We illustrate that the gravitomagnetic force at the edge of a galaxy can be in accord with either GVT or ΛCDM but not both. We also study the physics of the GVT theory around the gravitational saddle point of the Sun and Jupiter system. We notice that the conclusive refusal of the GVT theory demands measuring either both of the gravitoelectric and gravitomagnetic fields inside the Sun-Jupiter MOND window, or the gravitoelectric field inside two different solar GVT MOND windows. The GVT theory, however, will be favored by observing an anomaly in the gravitoelectric field inside a single MOND window.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Disconnected covariance of 2-point functions in large-scale structure

Optimal analyses using the 2-point functions of large-scale structure probes require accurate covariance matrices. A covariance matrix of the 2-point function comprises the disconnected part and the connected part. While the connected covariance only becomes important on small scales, the disconnected covariance is dominant on large scales, where the survey window has a significant impact. In this work, we develop an analytical method to compute the disconnected covariance, accounting for the window effect. Derived under the flat-sky approximation, our formalism is applicable to wide surveys by swapping in the curved-sky window functions. Our method works for both the power spectrum and the correlation function, and applies to the covariances of various probes including the multipoles and the wedges of 3D clustering, the angular and the projected statistics of clustering and shear, as well as the cross covariances between different probes. We verify the analytic covariance against the sample covariance from the galaxy mock simulations in two test cases: (1) the power spectrum multipole covariance, and (2) the joint covariance of the projected correlation function and the correlation function multipoles. Our method achieve good agreement with the mocks, while at a negligible computational cost. Unlike mocks, our analytic covariance is free of sampling noise, which often leads to numerical problems and the need to inflate the errors. In addition, our method can use the best-fit power spectrum as input, in contrast to the standard procedure of using a fiducial model that may deviate significantly from the truth. We also show that a naive diagonal power spectrum covariance underestimates the signal-to-noise ratio compared to our analytic covariance. The code that accompanies this paper is available at https://github.com/eelregit/covdisc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Search For Electron-Antineutrinos Associated With 2 Gravitational-Wave Events GW150914, GW151012, GW151226, 3 GW170104, GW170608, GW170814, and GW170817 at Daya Bay

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using the data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of ±10 s, ±500 s, and ±1000 s relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of (1.13 - 2.44)×1011 cm-2 at 5 MeV to 8.0×107cm-2at 100 MeV for the three timewindows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be (5.4 - 7.0)×109 cm-2 for the three time windows.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Achieving a robust grassy-ELM operation regime in CFETR

We have identified a robust grassy-edge localized mode (ELM) operation regime for future tokamak reactors. The regime exists within a pedestal top electron collisionality (ν*) window at high global poloidal beta (β p ). The existence of an upper ν* limit for grassy-ELMs is consistent with results previously reported in experiments (Oyama et al 2010 Nucl. Fusion 50 064014), while the existence of a lower ν* limit has not been reported previously. Using EPED and BOUT++, a theoretical model that quantitatively explains the physics of the grassy-ELMs within the window, which distinguishes them from the small mixed-ELMs at lower ν*, is presented for the first time. A peeling-ballooning stability boundary is obtained by scanning the operating density space. The change in density corresponds to a change in ν* that affects the pedestal bootstrap current. High β p leads to a strong Shafranov shift, which affects the flux surface averaged pressure drive. The two effects combine to create a peeling-dominated window in intermediate ν* buffered by ballooning-dominated regimes. Only the peeling-dominated regime shows a cyclic behavior in the perturbed pressure during the nonlinear simulation of an ELM crash, reminiscent of grassy-ELM dynamics. Similarly, the energy released across the separatrix is demonstrated to be significantly smaller. The quick recovery of the ELM crash is explainable by the rapid rise of a low n kink-peeling instability when the pedestal current I ped exceeds a threshold at high β p . It minimizes the excursion beyond marginal stability and is absent in the ballooning-dominated regime. Comparison with recent experiments over a range of β p and ν* strongly supports the physical picture proposed by the modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A baseline for ensemble-based, time-resolved inflow reconstruction for a single turbine using large-eddy simulations and latent diffusion models

We are interested in reconstructing winds flowing through a turbine on a second-by-second basis over a 10 min window. Previously, we developed a machine learning algorithm that takes in a snapshot of wind speed measurements and generates ensembles of three-dimensional wind field estimates. Here, we use these estimates as initial conditions in large-eddy simulations and reconstruct atmospheric and turbine response dynamic quantities in a synthetic field campaign. In doing so, we establish a baseline for model validation that future time-aware data assimilation techniques will be compared to. In turbine-free case studies, ground truth wind speeds consistently fall within our estimated wind speed distribution for the first 100 s after the simulation start. In simulations with turbines, the wind estimates show a small bias of 0.10 m s 1 and good correlation of 0.80 during the first 100 s. During this window, our estimates of the Blade 1 bending moment and generator power typically span the ground truth, with the estimate of the former performing better overall. In summary, this approach shows promise as a stand-alone technique for reconstructing real-world inflow and turbine dynamics in 1-2 min windows and as a foundation for future time-aware data assimilation techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Frequency domain analysis of errors in cross-correlations of ambient seismic noise

We analyse random errors (variances) in cross-correlations of ambient seismic noise in the frequency domain, which differ from previous time domain methods. Extending previous theoretical results on ensemble averaged cross-spectrum, we estimate confidence interval of stacked cross-spectrum of finite amount of data at each frequency using non-overlapping windows with fixed length. The extended theory also connects amplitude and phase variances with the variance of each complex spectrum value. Analysis of synthetic stationary ambient noise is used to estimate the confidence interval of stacked cross-spectrum obtained with different length of noise data corresponding to different number of evenly spaced windows of the same duration. This method allows estimating Signal/Noise Ratio (SNR) of noise cross-correlation in the frequency domain, without specifying filter bandwidth or signal/noise windows that are needed for time domain SNR estimations. Based on synthetic ambient noise data, we also compare the probability distributions, causal part amplitude and SNR of stacked cross-spectrum function using one-bit normalization or pre-whitening with those obtained without these pre-processing steps. Natural continuous noise records contain both ambient noise and small earthquakes that are inseparable from the noise with the existing pre-processing steps. Using probability distributions of random cross-spectrum values based on the theoretical results provides an effective way to exclude such small earthquakes, and additional data segments (outliers) contaminated by signals of different statistics (e.g. rain, cultural noise), from continuous noise waveforms. This technique is applied to constrain values and uncertainties of amplitude and phase velocity of stacked noise cross-spectrum at different frequencies, using data from southern California at both regional scale (~35 km) and dense linear array (~20 m) across the plate-boundary faults. A block bootstrap resampling method is used to account for temporal correlation of noise cross-spectrum at low frequencies (0.05–0.2 Hz) near the ocean microseismic peaks.

58 GEOSCIENCES↗
DOE OSTI2018

Umbrella sampling: a powerful method to sample tails of distributions

We present the umbrella sampling (US) technique and show that it can be used to sample extremely low-probability areas of the posterior distribution that may be required in statistical analyses of data. In this approach sampling of the target likelihood is split into sampling of multiple biased likelihoods confined within individual umbrella windows. We show that the US algorithm is efficient and highly parallel and that it can be easily used with other existing Markov Chain Monte Carlo (MCMC) samplers. The method allows the user to capitalize on their intuition and define umbrella windows and increase sampling accuracy along specific directions in the parameter space. Alternatively, one can define umbrella windows using an approach similar to parallel tempering. We provide a public code that implements US as a standalone PYTHON package. We present a number of tests illustrating the power of the US method in sampling low-probability areas of the posterior and show that this ability allows a considerably more robust sampling of multimodal distributions compared to standard sampling methods. We also present an application of the method in a real world example of deriving cosmological constraints using the supernova type Ia data. We show that US can sample the posterior accurately down to the approximate to 15 sigma credible region in the Omega(m) - Omega(Lambda) plane, while for the same computational effort the affine-invariant MCMC sampling implemented in the emcee code samples the posterior reliably only to approximate to 3 sigma.

cosmology: cosmological parameters↗
DOE OSTI2020

Density weighted angular redshift fluctuations: a new cosmological observable

We propose the use of angular fluctuations in the galaxy redshift field as a new way to extract cosmological information in the Universe. This new probe $\delta z (\hat{\mathbf {n}})$ consists of the statistics of sky maps built by projecting redshifts under a Gaussian window of width σz centred upon a redshift zobs, and weighted by the galaxy density field. We compute the angular power spectrum of the $\delta z (\hat{\mathbf {n}})$ field in both numerical simulations and in linear perturbation theory. From these, we find that the $\delta z (\hat{\mathbf {n}})$ field (i) is sensitive to the underlying density and peculiar velocity fields; (ii) is highly correlated, at the $\gtrsim 60\, {{\ \rm per\ cent}}$ level, to the line-of-sight projected peculiar velocity field; (iii) for narrow windows (σz < 0.03), it is almost completely uncorrelated to the projected galaxy angular density field under the same redshift window; and (iv) it is largely unaffected by multiplicative and additive systematic errors on the observed number of galaxies that are redshift-independent over ~σz. We conclude that $\delta z (\hat{\mathbf {n}})$ is a simple and robust tomographic measure of the cosmic density and velocity fields, complementary to angular clustering, that will contribute to more complete exploitations of current and upcoming galaxy redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Hadronic Vacuum Polarization for the Muon 𝑔 − 2 from Lattice QCD: Long-Distance and Full Light-Quark Connected Contribution

We present results for the dominant light-quark connected contribution to the long-distance window of the hadronic vacuum polarization (HVP) contribution to the muon 𝑔 − 2 from lattice quantum chromodynamics. Specifically, with a new determination of the lattice scale on MILC’s physical-mass HISQ ensembles, using the Ω− baryon mass, we obtain a result of 𝑎$^{𝑙⁢𝑙,LD}_{𝜇}$⁡(conn) = 400.2⁢(2.3) stat ⁢(3.7) syst ⁢[4.3] total × 10 −10 . Summing this result with our recent determinations of the light-quark connected contributions to the short- and intermediate-distance windows, we obtain a subpercent precision determination of the light-quark-connected contribution to HVP of 𝑎$^{𝑙⁢𝑙}_{𝜇}$⁡(conn) = 655.2⁢(2.3) stat ⁢(3.9) syst ⁢[4.5]total×10 −10 . Finally, as a consistency check, we verify that an independent analysis of the full contribution is in agreement with the sum of individual windows. We discuss our future plans for improvements of our HVP calculations to meet the target precision of the Fermilab 𝑔 − 2 experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗