Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Do Arctic mixed-phase clouds sometimes dissipate due to insufficient aerosol? Evidence from comparisons between observations and idealized simulations

Abstract. Mixed-phase clouds are ubiquitous in the Arctic. These clouds can persist for days and dissipate in a matter of hours. It is sometimes unknown what causes this sudden dissipation, but aerosol–cloud interactions may be involved. Arctic aerosol concentrations can be low enough to affect cloud formation and structure, and it has been hypothesized that, in some instances, concentrations can drop below some critical value needed to maintain a cloud. We use observations from a Department of Energy ARM site on the northern slope of Alaska at Oliktok Point (OLI), the Arctic Summer Cloud Ocean Study (ASCOS) field campaign in the high Arctic Ocean, and the Integrated Characterisation of Energy, Clouds, Atmospheric state, and Precipitation at Summit – Aerosol Cloud Experiment (ICECAPS-ACE) project at the NSF (National Science Foundation) Summit Station in Greenland (SMT) to identify one case per site where Arctic boundary layer clouds dissipated coincidentally with a decrease in surface aerosol concentrations. These cases are used to initialize idealized large eddy simulations (LESs) in which aerosol concentrations are held constant until, at a specified time, all aerosols are removed instantaneously – effectively creating an extreme case of aerosol-limited dissipation which represents the fastest a cloud could possibly dissipate via this process. These LESs are compared against the observed data to determine whether cases could, potentially, be dissipating due to insufficient aerosol. The OLI case's observed liquid water path (LWP) dissipated faster than its simulation, indicating that other processes are likely the primary drivers of the dissipation. The ASCOS and SMT observed LWP dissipated at similar rates to their respective simulations, suggesting that aerosol-limited dissipation may be occurring in these instances. We also find that the microphysical response to this extreme aerosol forcing depends greatly on the specific case being simulated. Cases with drizzling liquid layers are simulated to dissipate by accelerating precipitation when aerosol is removed while the case with a non-drizzling liquid layer dissipates quickly, possibly glaciating via the Wegener–Bergeron–Findeisen (WBF) process. The non-drizzling case is also more sensitive to ice-nucleating particle (INP) concentrations than the drizzling cases. Overall, the simulations suggest that aerosol-limited cloud dissipation in the Arctic is plausible and that there are at least two microphysical pathways by which aerosol-limited dissipation can occur.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global Dust Cycle and Direct Radiative Effect in E3SM Version 1: Impact of Increasing Model Resolution

Quantification of dust aerosols in Earth System Models (ESMs) has important implications for water cycle and biogeochemistry studies. This study examines the global life cycle and direct radiative effects (DREs) of dust in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1), and the impact of increasing model resolution both horizontally and vertically. The default 1° E3SMv1 captures the spatial and temporal variability in the observed dust aerosol optical depth (DAOD) reasonably well, but overpredicts dust absorption in the shortwave (SW). Simulations underestimate the dust vertical and long-range transport, compared with the satellite dust extinction profiles. After updating dust refractive indices and correcting for a bias in partitioning size-segregated emissions, both SW cooling and longwave (LW) warming of dust simulated by E3SMv1 are increased and agree better with other recent studies. The estimated net dust DRE of –0.42 Wm –2 represents a stronger cooling effect than the observationally based estimate –0.2 Wm –2 (–0.48 to +0.2), due to a smaller LW warming. Constrained by a global mean DAOD, model sensitivity studies of increasing horizontal and vertical resolution show strong influences on the simulated global dust burden and lifetime primarily through the change of dust dry deposition rate; there are also remarkable differences in simulated spatial distributions of DAOD, DRE, and deposition fluxes. Thus, constraining the global DAOD is insufficient for accurate representation of dust climate effects, especially in transitioning to higher- or variable-resolution ESMs. Better observational constraints of dust vertical profiles, dry deposition, size, and LW properties are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A new process-based and scale-aware desert dust emission scheme for global climate models – Part I: Description and evaluation against inverse modeling emissions

Desert dust accounts for most of the atmosphere's aerosol burden by mass and produces numerous important impacts on the Earth system. However, current global climate models (GCMs) and land-surface models (LSMs) struggle to accurately represent key dust emission processes, in part because of inadequate representations of soil particle sizes that affect the dust emission threshold, surface roughness elements that absorb wind momentum, and boundary-layer characteristics that control wind fluctuations. Furthermore, because dust emission is driven by small-scale (~ 1 km or smaller) processes, simulating the global cycle of desert dust in GCMs with coarse horizontal resolutions (~ 100 km) presents a fundamental challenge. This representation problem is exacerbated by dust emission fluxes scaling nonlinearly with wind speed above a threshold wind speed that is sensitive to land-surface characteristics. Here, we address these fundamental problems underlying the simulation of dust emissions in GCMs and LSMs by developing improved descriptions of (1) the effect of soil texture on the dust emission threshold, (2) the effects of nonerodible roughness elements (both rocks and green vegetation) on the surface wind stress, and (3) the effects of boundary-layer turbulence on driving intermittent dust emissions. We then use the resulting revised dust emission parameterization to simulate global dust emissions in a standalone model forced by reanalysis meteorology and land-surface fields. We further propose (4) a simple methodology to rescale lower-resolution dust emission simulations to match the spatial variability of higher-resolution emission simulations in GCMs. The resulting dust emission simulation shows substantially improved agreement against regional dust emissions observationally constrained by inverse modeling. We thus find that our revised dust emission parameterization can substantially improve dust emission simulations in GCMs and LSMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Boosting H I -Galaxy Cross-Clustering Signal through Higher-Order Cross-Correlations

After reionization, neutral hydrogen (${\rm H\, \small {I}}$) traces the large-scale structure (LSS) of the Universe, enabling ${\rm H\, \small {I}}$ intensity mapping (IM) to capture the LSS in 3D and constrain key cosmological parameters. We present a new framework utilizing higher-order cross-correlations to study ${\rm H\, \small {I}}$ clustering around galaxies, tested using real-space data from the IllustrisTNG300 simulation. This approach computes the joint distributions of k-nearest neighbor (kNN) optical galaxies and the ${\rm H\, \small {I}}$ brightness temperature field smoothed at relevant scales (the kNN-field framework), providing sensitivity to all higher-order cross-correlations, unlike two-point statistics. To simulate ${\rm H\, \small {I}}$ data from actual surveys, we add random thermal noise and apply a simple foreground cleaning model, filtering out Fourier modes of the brightness temperature field with k ∥ < k min,∥ . Under current levels of thermal noise and foreground cleaning, typical of a Canadian Hydrogen Intensity Mapping Experiment (CHIME)-like survey, the ${\rm H\, \small {I}}$-galaxy cross-correlation signal in our simulations, using the kNN-field framework, is detectable at >30σ across r = [3, 12] h –1 Mpc. In contrast, the detectability of the standard two-point correlation function (2PCF) over the same scales depends strongly on the foreground filter: a sharp k ∥ filter can spuriously boost detection to 8σ due to position-space ringing, whereas a less sharp filter yields no detection. Nonetheless, we conclude that kNN-field cross-correlations are robustly detectable across a broad range of foreground filtering and thermal noise conditions, suggesting their potential for enhanced constraining power over 2PCFs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Enhancing the sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment to low energy signals

Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Enhancing the sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment to low energy signals

Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$.

Akerib, DS↗
DOE OSTI2026

Titanium alloy response sensitivity to variations in spectral reconstructions of National Ignition Facility xenon line-emission x-ray sources

Thermomechanical shock experiments on the National Ignition Facility (NIF) aim to study high strain rate dynamic material response. In such experiments, the NIF laser is used to generate high fluence x-ray emission sources, which irradiate material samples of interest. Under sufficiently high x-ray energy deposition, thermomechanical impulses are generated in the materials. While it is known that the characteristics of x-ray generated impulses vary as a function of incident x-ray spectra, it remains unclear how spectral assumptions and uncertainties in NIF spectral reconstructions affect our interpretation of impulsive loading. Here, in this paper, we simulate the response of a standard titanium alloy baseline sample to synthetic analytically derived and measured NIF xenon line-emission x-ray sources with a radiation hydrodynamics code. We vary the source spectral characteristics based on different source reconstruction techniques to understand the resulting variation in baseline sample response and compare the simulated response with experimental results. We find that the response is highly sensitive to assumptions made about the spectral contents and that knowledge of spectral uncertainties bounds our understanding of the resulting material response. The results of this effort help to extend our ability to use baseline material samples to extract quantitative properties from x-ray experiments on the NIF.

Alloys↗
DOE OSTI2022

Increasing sensitivity of dryland vegetation greenness to precipitation due to rising atmospheric CO 2

Water availability plays a critical role in shaping terrestrial ecosystems, particularly in low- and mid-latitude regions. The sensitivity of vegetation growth to precipitation strongly regulates global vegetation dynamics and their responses to drought, yet sensitivity changes in response to climate change remain poorly understood. Here we use long-term satellite observations combined with a dynamic statistical learning approach to examine changes in the sensitivity of vegetation greenness to precipitation over the past four decades. We observe a robust increase in precipitation sensitivity (0.624% yr –1 ) for drylands, and a decrease (–0.618% yr –1 ) for wet regions. Using model simulations, we show that the contrasting trends between dry and wet regions are caused by elevated atmospheric CO 2 (eCO 2 ). eCO 2 universally decreases the precipitation sensitivity by reducing leaf-level transpiration, particularly in wet regions. However, in drylands, this leaf-level transpiration reduction is overridden at the canopy scale by a large proportional increase in leaf area. The increased sensitivity for global drylands implies a potential decrease in ecosystem stability and greater impacts of droughts in these vulnerable ecosystems under continued global change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of Urbanization on Convection, Lightning, and Precipitation Over the Houston Metropolitan Area: Case Study Simulation From the TRACER Campaign

This study investigates the effects of urbanization, specifically land use change and anthropogenic emissions (AE), on convection, lightning, and surface precipitation for a case of summertime sea‐breeze convection observed over the Houston metropolitan area. The unique capabilities of the NASA‐Unified Weather Research and Forecasting model allows us to conduct a series of sensitivity experiments with complex configurations, in particular including multi‐year land model spin‐up simulations, treatment of aerosols and their precursors, and explicit cloud charging and lightning. The simulation results show that urban land use primarily alters the temporal evolution of convection, lightning, and surface precipitation, leading to late afternoon thunderstorm development. The decrease in latent heat flux from the land surface caused by urbanization weakens convection in the early afternoon, while a condition suitable for convection development is maintained in the late afternoon due to less stabilization of the lower troposphere by the weaker convection development and high sensible heat flux from the surface. On the other hand, anthropogenic aerosols directly enhance convection, lightning, and surface precipitation by increasing convective updrafts due to the aerosol‐induced convective invigoration. The combined effects of urban land use and AE lead to even stronger thunderstorms in the late afternoon, mostly consistent with observations. These results indicate that urbanization increases the probability of late afternoon thunderstorms over the Houston area during the summer season. Advanced weather forecasting models that incorporate these urbanization effects might support sustainable urban planning to better mitigate the impacts of urbanization on local weather and public safety.

Iguchi, T. [Univ. of Maryland, College Park, MD (U↗
DOE OSTI2025

The DESI-Lensing Mock Challenge: large-scale cosmological analysis of 3x2-pt statistics

The current generation of large galaxy surveys will test the cosmological model by combining multiple types of observational probes. Realising the statistical promise of these new datasets requires rigorous attention to all aspects of analysis including cosmological measurements, modelling, covariance and parameter likelihood. In this paper we present the results of an end-to-end simulation study designed to test the analysis pipeline for the combination of the Dark Energy Spectroscopic Instrument (DESI) Year 1 galaxy redshift dataset and separate weak gravitational lensing information from the Kilo-Degree Survey, Dark Energy Survey and Hyper-Suprime-Cam Survey. Our analysis employs the 3x2-pt correlation functions including cosmic shear and galaxy-galaxy lensing, together with the projected correlation function of the spectroscopic DESI lenses. We build realistic simulations of these datasets including galaxy halo occupation distributions, photometric redshift errors, weights, multiplicative shear calibration biases and magnification. We calculate the analytical covariance of these correlation functions including the Gaussian, noise and super-sample contributions, and show that our covariance determination agrees with estimates based on the ensemble of simulations. We use a Bayesian inference platform to demonstrate that we can recover the fiducial cosmological parameters of the simulation within the statistical error margin of the experiment, investigating the sensitivity to scale cuts. This study is the first in a sequence of papers in which we present and validate the large-scale 3x2-pt cosmological analysis of DESI-Y1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

EDGE2D-EIRENE and ERO2.0 predictions of nitrogen molecular break-up and transport in the divertor of JET low-confinement mode plasmas

EDGE2D-EIRENE simulations of nitrogen-seeded partially detached JET L-mode plasmas show that the divertor N I to N V radiation distributions are highly sensitive to the upstream electron density, and less sensitive to changes in the assumed cross-field particle diffusivity for nitrogen ions. The EDGE2D-EIRENE simulations reproduce the peak intensities of N I to N V as measured by vertically viewing divertor filterscopes to within 50 % for a narrow range of the upstream electron density within the experimental uncertainties, while the predicted profiles are narrower than the measured ones. Including nitrogen atoms only in the ERO2.0 simulations implies lower N III and N IV peak intensities by one and two thirds, respectively, compared to EDGE2D-EIRENE. If nitrogen is assumed to recycle exclusively as molecules instead of atoms, ERO2.0 predicts that the N III and N IV intensities in the divertor increase by up to a factor of two and that the time-averaged, volume-integrated number of N 5+ to N 7+ ions in the plasma also increases by approximately a factor of 2.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sensitivity Study of Type-I X-ray Bursts to Nuclear Reaction Rates

Type-I X-ray bursts are frequently occurring thermonuclear runaways on the surface of an accreting neutron star. The different nuclear reactions that empower the cataclysmic event play a key role in accurately comparing theoretical models and observations. We investigate the effect of the uncertainties in the nuclear reactions using a ONEZONE model for a set of different compositions and accretion rates that are within the range of the standard observed burst sources. A combination of a full-reaction network and a semi-analytic model is used to obtain the conditions at the time of X-ray burst ignition via simulating the settling process of the accreted material. We then evaluate the sensitivity of the X-ray burst model by varying the proton and alpha-induced reaction rates in JINA REACLIBV2.2 within representative nuclear physics uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Deficient precipitation sensitivity to Sahel land surface forcings among CMIP5 models

Abstract The overall performance of the simulated seasonal precipitation response to local terrestrial forcings, namely vegetation abundance and soil moisture, in the Sahel among the Coupled Model Intercomparison Project Phase Five (CMIP5) Earth System Models (ESMs) is systematically investigated and compared with its observational counterpart using a multivariate statistical method. The observed seasonal precipitation response is evaluated against a large ensemble of observational, reanalysis, and satellite data sets to provide quantification of uncertainties. The behaviour of models with and without a Dynamic Global Vegetation Model (DGVM) component is also explored, along with the mechanisms responsible for terrestrial feedback on rainfall. In general, the CMIP5 models can reasonably capture the seasonal evolution of Sahel precipitation and soil moisture, albeit with wet biases during the pre‐monsoon period and dry biases during the peak monsoon period. The non‐DGVM ESMs simulate comparable leaf area indices (LAIs) with observations, while DGVM‐enabled ESMs simulate too much year‐round LAI. The variance of precipitation that is attributed to oceanic forcings in CMIP5 is comparable with observations; however, the variance of precipitation that is attributed to terrestrial forcings is smaller in CMIP5 models than observed, especially for non‐DGVM ESMs. CMIP5 models, especially those without DGVMs, undervalue precipitation's observed response strength to soil moisture anomalies. In both observations and CMIP5 models, none of the atmospheric variables show significant responses to direct vegetation forcing, except for the response in transpiration. Although vegetation has minimal direct effect on the atmospheric state, it can affect the atmosphere by modifying soil moisture and transpiration rate indirectly, which helps explain the more realistic simulation of rainfall in DGVM‐enabled ESMs than non‐DGVM ESMs. Coupling of an ESM to a DGVM is critical in generating reasonable land–atmosphere feedback and examining future ecological and climatic changes over the Sahel.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Automated detector simulation and reconstruction parametrization using machine learning

Rapidly applying the effects of detector response to physics objects (e.g. electrons, muons, showers of particles) is essential in high energy physics. Presently available tools for the transformation from truth-level physics objects to reconstructed detector-level physics objects involve manually defining resolution functions. These resolution functions are typically derived in bins of variables that are correlated with the resolution (e.g. pseudorapidity and transverse momentum). This process is time consuming, requires manual updates when detector conditions change, and can miss important correlations. Machine learning offers a way to automate the process of building these truth-to-reconstructed object transformations and can capture complex correlation for any given set of input variables. Such machine learning algorithms, with sufficient optimization, could have a wide range of applications: improving phenomenological studies by using a better detector representation, allowing for more efficient production of Geant4 simulation by only simulating events within an interesting part of phase space, and studies on future experimental sensitivity to new physics.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Projected WIMP sensitivity of the LUX-ZEPLIN dark matter experiment

LUX-ZEPLIN (LZ) is a next generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with Weakly Interacting Massive Particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above $1.6 \times 10^{-48}$ cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$ discovery potential is projected reaching cross sections below the existing and projected exclusion limits of similar experiments that are currently operating. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of $2.7 \times 10^{-43}$ cm$^{2}$ ($8.1 \times 10^{-42}$ cm$^{2}$) for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗
DOE OSTI2021

Dakota A Multilevel Parallel Object-Oriented Framework for Design Optimization Parameter Estimation Uncertainty Quantification and Sensitivity Analysis (V.6.14) (Theory Manual)

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota's iterative analysis capabilities.

97 MATHEMATICS AND COMPUTING↗