Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variance to Mean”

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 145 records · Page 8

Probing triaxial deformation of atomic nuclei in high-energy heavy ion collisions

Most atomic nuclei are deformed with a quadrupole shape described by its overall strength β 2 and triaxiality γ. The deformation can be accessed in high-energy heavy ion collisions by measuring the collective flow response of the produced quark-gluon plasma to the eccentricity ε 2 and the density gradient d ⊥ in the initial state. Using an analytical estimate and a Glauber model, I show that the variances $\langleε^{2}_{2}\rangle$ or $\langle$(δd ⊥ /d ⊥ ) 2 $\rangle$ and skewnesses $\langleε^{2}_{2}δd_{⊥}/d_{⊥}\rangle$ or $\langle(δd_{⊥}/d_{⊥})^{3}\rangle$ have a simple analytical form of a'+b'β$^{2}_{2}$ and a'+[b'+c' cos(3γ)]β$^{3}_{2}$, respectively. From these, I constructed several normalized skewnesses to isolate the γ dependence from that of β 2 and show that the correlations between a normalized skewness and a variance can constrain simultaneously β 2 and γ. Assuming a linear relation with elliptic flow v 2 and mean-transverse momentum [p T ] of final-state particles, v 2 ∝ ε 2 and δ [p T ] / [p T ] ∝ δd ⊥ /d ⊥ , similar conclusions are also expected for the variances and skewnesses of v 2 and [p T ], i.e., a + bβ$^{2}_{2}$ for $\langle v^{2}_{2}\rangle$ and $\langle(δ[p_{T}] / [p_{T}])^{2}\rangle$ and a + [b + c cos (3γ)] β$^{3}_{2}$ for $\langle v^{2}_{2}δ [p_{T}] / [p_{T}]\rangle$ or $\langle(δ[p_{T}] / [p_{T}])^{3}\rangle$. My findings motivate a dedicated system scan of high-energy heavy ion collisions at RHIC and LHC to measure triaxiality of atomic nuclei: one first determines the coefficients b and c by collisions of isobaric near prolate nuclei, cos (3γ) ≈ 1, and near oblate nuclei, cos (3γ) ≈ –1 , with known β 2 values, followed by collisions of other species of interest with similar mass number. The (β 2 , γ) values for this species can be inferred directly from the measured variance and skewness observables from these collisions. Finally, the results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine learning and process-based modeling of spatiotemporal changes in active layer thickness across Alaska

Permafrost degradation poses a growing threat to infrastructure stability and ecosystem resilience in the rapidly warming Arctic. We investigated the spatiotemporal dynamics of active layer thickness (ALT) across Alaska by integrating field observations, environmental datasets, a physically based Stefan model, and machine learning (ML) techniques. Using weather projections from the Coupled Model Intercomparison Project Phase 6 under two Shared Socioeconomic Pathways (SSP 2-4.5 and SSP 5-8.5), we assessed ALT sensitivity to projected future weather conditions. The random forest (RF) model outperformed the Stefan approach in predicting ALT on the training dataset (R² = 0.84 vs. 0.53) but demonstrated lower generalizability on the test dataset (R² = 0.24 vs. 0.54). The root mean square error (RMSE) for the RF model for training and testing ranged from 14 to 22 cm, compared to 17 and 18 cm for the Stefan model. Variable importance analysis revealed that mean annual temperature and slope angle were the strongest predictors of ALT, accounting for 19% and 18% of the variance, respectively, followed by sediment transport index (14%) and stream power index (11%). Comparative analysis of baseline ALT predictions showed the Stefan model tended to project a thicker active layer (mean ± SD: 65 ± 16 cm), compared to the RF model (mean ± SD: 59 ± 8.8) cm). Both models indicated a latitudinal gradient in ALT, with shallower depths at higher latitudes. Projected ALT increases by 2100 were estimated at 3.3 ± 2.2 cm under SSP 2-4.5 and 5.9 ± 4.0 cm under SSP 5-8.5 for the ML model, whereas the Stefan model projected substantially larger increases of 13 ± 2.6 cm (SSP 2-4.5) and 28 ± 4.4 cm (SSP5-8.5). Spatial analysis showed the greatest ALT increases in northern Alaska, with relatively smaller changes in southern regions. These findings highlight the complex, multifactorial nature of ALT dynamics and the value of hybrid modeling approaches. As rising temperatures accelerate permafrost thaw, changes in ALT can disrupt ecosystems, damage infrastructures, and enhance the release of stored soil carbon, highlighting the urgent need for improved predictive capabilities to inform adaptation strategies in the Arctic.

Climate sciences↗

Enabling Robust Exoplanet Atmospheric Retrievals with Gaussian Processes

Atmospheric retrievals are essential tools for interpreting exoplanet transmission and eclipse spectra, enabling quantitative constraints on the chemical composition, aerosol properties, and thermal structure of planetary atmospheres. The James Webb Space Telescope (JWST) offers unprecedented spectral precision, resolution, and wavelength coverage, unlocking transformative insights into the formation, evolution, climate, and potential habitability of planetary systems. However, this opportunity is accompanied by challenges: modeling assumptions and unaccounted-for noise or signal sources can bias retrieval outcomes and their interpretation. To address these limitations, we introduce a Gaussian process (GP)-aided atmospheric retrieval framework that flexibly accounts for unmodeled features and correlated noise in exoplanet spectra. We validate this method on synthetic JWST observations, and show that GP-aided retrievals reduce bias in inferred abundances and better capture model–data mismatches than traditional approaches. We also introduce the concept of mean squared error to quantify the trade-off between bias and variance, arguing that this metric more accurately reflects retrieval performance than bias alone. We then reanalyze the NIRISS/SOSS JWST transmission spectrum of WASP-96 b, finding that GP-aided retrievals yield broader constraints on CO 2 and H 2 O, possibly alleviating tension between previous retrieval results and equilibrium predictions. Our GP framework provides precise and accurate constraints while highlighting regions where models fail to explain the data. As JWST matures and future facilities come online, a deeper understanding of the limitations of both data and models will be essential, and GP-enabled retrievals like the one presented here offer a principled path forward.

Rotman, Yoav [Arizona State Univ., Tempe, AZ (Unit↗

A temporal analysis and response to nitrate availability of 3D root system architecture in diverse pennycress (Thlaspi arvense L.) accessions

Introduction Roots have a central role in plant resource capture and are the interface between the plant and the soil that affect multiple ecosystem processes. Field pennycress ( Thlaspi arvense L.) is a diploid annual cover crop species that has potential utility for reducing soil erosion and nutrient losses; and has rich seeds (30-35% oil) amenable to biofuel production and as a protein animal feed. The objective of this research was to (1) precisely characterize root system architecture and development, (2) understand plastic responses of pennycress roots to nitrate nutrition, (3) and determine genotypic variance available in root development and nitrate plasticity. Methods Using a root imaging and analysis pipeline, the 4D architecture of the pennycress root system was characterized under four nitrate regimes, ranging from zero to high nitrate concentrations. These measurements were taken at four time points (days 5, 9, 13, and 17 after sowing). Results Significant nitrate condition response and genotype interactions were identified for many root traits, with the greatest impact observed on lateral root traits. In trace nitrate conditions, a greater lateral root count, length, density, and a steeper lateral root angle was observed compared to high nitrate conditions. Additionally, genotype-by-nitrate condition interaction was observed for root width, width:depth ratio, mean lateral root length, and lateral root density. Discussion These findings illustrate root trait variance among pennycress accessions. These traits could serve as targets for breeding programs aimed at developing improved cover crops that are responsive to nitrate, leading to enhanced productivity, resilience, and ecosystem service.

59 BASIC BIOLOGICAL SCIENCES↗

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗

Statistical uncertainties of the v n { 2 k } harmonics from Q cumulants

Analytic formulas to calculate statistical uncertainties of v n { 2 k } harmonics extracted from the Q cumulants are presented. The Q cumulants are multivariate polynomial functions of the weighted means of 2m-particle azimuthal correlations, $\langle\langle 2 m \rangle\rangle$. Variances and covariances of the $\langle\langle 2 m \rangle\rangle$ are included in the analytic formulas of the uncertainties that can be calculated simultaneously with the calculations of the v n { 2 k } harmonics. The calculations are performed using a simple toy model, which roughly simulates elliptic flow azimuthal anisotropy with magnitudes around 0.05. The results are compared with the results obtained by the many data subsets, and by the bootstrapping method. The first one is a common way of estimation of the statistical uncertainties of the v n { 2 k } harmonics in a real experiment. In order to increase precision in the measurement of the v n { 2 k } harmonics, a large number of 15 000 subsets and the same number of the resampling in the bootstrap method is used. Unlike the other ways of the analytic calculation of the v n { 2 k } statistical uncertainties, our proposal that includes the use of squared weights in the calculation of both the variances and covariances, gives the best agreement with the results obtained from the subsets and bootstrap method. Additionally, a recurrence equation between Q cumulants of any order is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Using ground-based lidar data to investigate the water–vapor budget in the daytime atmospheric boundary layer

The moisture advection term in the water–vapor budget equation is investigated with a combination of a vertically-staring water–vapor lidar and Doppler lidar systems. These instruments make it possible to get the mean profile of moisture tendency and the latent heat flux (LHF) divergence. We use data of the Land–Atmosphere Feedback Experiment (LAFE) at the Atmospheric Radiation Measurement (ARM) Program’s Southern Great Plains (SGP) site, Oklahoma, USA, collected on 30 August 2017 between 15 and 24 UTC, which corresponds to 09 to 18 LT. The lidars provide turbulence resolving profiles of moisture and vertical wind fluctuations. The LHF profile is derived from the covariance of these moisture and vertical wind fluctuations. The mean boundary layer height z i is determined from the peak of the moisture variance. The results demonstrate that the combination of two remote sensing instruments can be applied for determining the dominant water–vapor budget terms, namely moisture tendency, latent heat flux divergence and moisture advection.

Advection↗

A New Coupled Biogeochemical Modeling Approach Provides Accurate Predictions of Methane and Carbon Dioxide Fluxes Across Diverse Tidal Wetlands

Abstract Tidal wetlands provide valuable ecosystem services, including storing large amounts of carbon. However, the net exchanges of carbon dioxide (CO 2 ) and methane (CH 4 ) in tidal wetlands are highly uncertain. While several biogeochemical models can operate in tidal wetlands, they have yet to be parameterized and validated against high‐frequency, ecosystem‐scale CO 2 and CH 4 flux measurements across diverse sites. We paired the Cohort Marsh Equilibrium Model (CMEM) with a version of the PEPRMT model called PEPRMT‐Tidal, which considers the effects of water table height, sulfate, and nitrate availability on CO 2 and CH 4 emissions. Using a model‐data fusion approach, we parameterized the model with three sites and validated it with two independent sites, with representation from the three marine coasts of North America. Gross primary productivity (GPP) and ecosystem respiration (R eco ) modules explained, on average, 73% of the variation in CO 2 exchange with low model error (normalized root mean square error (nRMSE) <1). The CH 4 module also explained the majority of variance in CH 4 emissions in validation sites ( R 2 = 0.54; nRMSE = 1.15). The PEPRMT‐Tidal‐CMEM model coupling is a key advance toward constraining estimates of greenhouse gas emissions across diverse North American tidal wetlands. Further analyses of model error and case studies during changing salinity conditions guide future modeling efforts regarding four main processes: (a) the influence of salinity and nitrate on GPP, (b) the influence of laterally transported dissolved inorganic C on R eco , (c) heterogeneous sulfate availability and methylotrophic methanogenesis impacts on surface CH 4 emissions, and (d) CH 4 responses to non‐periodic changes in salinity.

54 ENVIRONMENTAL SCIENCES↗

Comptonization by reconnection plasmoids in black hole coronae I: Magnetically dominated pair plasma

We perform 2D particle-in-cell simulations of reconnection in magnetically dominated electron–positron plasmas subject to strong Compton cooling. We vary the magnetization σ $\gg$ 1, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting black holes, which produce hard X-rays through Comptonization of seed soft photons. We find that the particle energy spectrum is dominated by a peak at mildly relativistic energies, which results from bulk motions of cooled plasmoids. The peak has a quasi-Maxwellian shape with an effective temperature of ~100 keV, which depends only weakly on the flow magnetization and the strength of radiative cooling. The mean bulk energy of the reconnected plasma is roughly independent of σ, whereas the variance is larger for higher magnetizations. The spectra also display a high-energy tail, which receives ~25 percent of the dissipated reconnection power for σ = 10 and ~40 percent for σ = 40. We complement our particle-in-cell studies with a Monte Carlo simulation of the transfer of seed soft photons through the reconnection layer, and find the escaping X-ray spectrum. The simulation demonstrates that Comptonization is dominated by the bulk motions in the chain of Compton-cooled plasmoids and, for σ ~ 10, yields a spectrum consistent with the typical hard state of accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

Fast-neutron-induced fission of Pu 240 and Pu 242

Herein we report the measurement of the total kinetic energy (TKE) release in the fast neutron induced fission of 240 Pu and 242 Pu. The results are compared to the predictions of the GEF model, the CGMF model, and the model of Denisov and Sedykh as well as previous exptl. work on these reactions. Our absolute measurements of the TKE release are in good agreement with the previous measurements of Nethaway et al. for the interaction of 14.8 MeV neutrons with 240 Pu [Phys. Rev. C16, 1907 (1977)] and of Winkelmann and Aumann for the interaction of 15 MeV neutrons with 242 Pu [Phys. Rev. C30, 934 (1984)]. The general trends of the measured TKE values agree with phenomenol. models but the variances of the TKE distributions are significantly less than predicted by various models. The mean postneutron emission TKE release decreases nonlinearly with increasing neutron energy and can be represented as TKE(MeV) = 175.8 ± 0.3 - (2.4 ± 0.8) lo g 10 E n - (1.4 ± 0.4) lo g 10 $E^{2}_{n}$ for 240 Pu and TKE(MeV) = 177.1 ± 0.3 - (1.2 ± 0.9) lo g 10 E n - (1.8 ± 0.5)lo g 10 $E_{n}^{2}$ for 242 Pu.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Genetic correlations between switchgrass performance in sward conditions and surrogate measures

Switchgrass (Panicum virgatum L.) breeders use spaced-seedling plots for evaluations of biomass yield, but prior studies are ambiguous about the relationship between performance in spaced plots and commercial seeded sward performance. This study used a split-plot incomplete block design to evaluate 52 half-sib switchgrass families across two years to determine the correlations among four different plot types (including seeded swards) and four individual plant traits. Plots (1.62 m 2 ) that contained 3, 6, and 18 plants had mean yields of 5.92, 8.04, and 11.06 Mg ha –1 , respectively, while seeded-sward plots had a mean of 5.37 Mg ha –1 . In a mixed-model analysis including family-, year-, and plot-type interactions with a spatial-spline correction, the largest proportion of family variance was accounted for by the main family effect (48%), while the plot type-by-family effect was only 17% of family variance. These results indicate that spaced-plant plots are reasonable surrogates for biomass performance in sward conditions. Yield reliability was 0.50 across all plot types with a mean of 0.46 for individual plot types. Surrogates for sward yield varied due to genetic background of the families. Yield of families derived from the Liberty cultivar (n = 32) had genetic correlations with winter survivorship scores (0.66). Diverse lowland ecotype families (n = 20) had genetic correlations with plant height and heading date. This study highlighted a major difference in performance among the plot types, specifically that seeded swards have depressed yield due to high plant density. This should encourage further research for the mechanism behind this difference in performance.

59 BASIC BIOLOGICAL SCIENCES↗

Denoising Autoencoder for Reconstructing Sensor Observation Data and Predicting Evapotranspiration: Noisy and Missing Values Repair and Uncertainty Quantification

Abstract Machine learning (ML) methods applied in scientific research often deal with interrelated features in high‐dimensional data. Reducing data noise and redundancy is needed to increase prediction accuracy and efficiency especially when dealing with data from field sensors. We explored an unsupervised learning method, the denoising autoencoder (DAE), to extract the underlying data structure from noisy raw data in the context of predicting hydrologic quantities from multiple field sensors. These sensors have intrinsic instrumental noise and occasional malfunctions that cause missing values. Our DAE neural network reconstructed meteorological sensor data containing noise and missing values to predict evapotranspiration in a mountainous watershed. The DAE reconstructed the sensor variables with a mean coefficient of determination value of 0.77 across 15 dimensions representing individual sensors. It reduced variance and bias uncertainties compared to a classical autoencoder model. The reconstruction quality varied across dimensions depending on their cross‐correlation and alignment with the underlying data structure. Uncertainties arising from the model structure were overall higher than those resulting from data corruption. We attached the DAE structure to a downstream ET‐prediction neural network in three formats and achieved reasonably accurate ET predictions . The use of the DAE notably reduced variance uncertainty in ET prediction. However, excessive variance reduction may be accompanied by an increase in bias due to the intrinsic bias‐variance tradeoff. Our method of evaluating and reducing uncertainties in aggregated data from different sources can be used to improve predictive models, process understanding, and uncertainty quantification for better water resource management. Plain Language Summary We present a machine learning method, namely the denoising autoencoder, which reduces the effects of data noise and missing values typically present in scientific data sets collected through sensor measurements. This method selects the most relevant information from noisy raw data collected by the instruments and fills in missing values. To demonstrate the effectiveness of our method, we applied it to predict evapotranspiration, a hydrologic variable that represents the water moved from the land surface to the atmosphere through a combination of evaporation and plant water use (transpiration). We also used a random sampling technique (the Monte Carlo method) to compare the uncertainty in the predictions when using the raw and noisy data versus the reconstructed data. The denoising process produced more accurate predictions of evapotranspiration with less uncertainty. Improved predictions of evapotranspiration can lead to a better understanding and accounting of water budgets. This ML approach is broadly suitable for a wide variety of applications that involve noisy sensor data with missing values. Key Points We used a denoising autoencoder (DAE) neural network to reduce noise in meteorological and soil sensor observations by on average We used Monte Carlo sampling to estimate the bias and variance of all model outputs, including uncertainty sources from data and the model We attached the DAE component to a downstream neural network to predict ET with the variance reduced by , compared to that without the DAE

denoising autoencoder↗

GALIC: hybrid multi-qubitwise pauli grouping for quantum computing measurement

Abstract Observable estimation is a core primitive in NISQ-era algorithms targeting quantum chemistry applications. To reduce the state preparation overhead required for accurate estimation, recent works have proposed various simultaneous measurement schemes to lower estimator variance. Two primary grouping schemes have been proposed: full commutativity (FC) and qubit-wise commutativity (QWC), with no compelling means of interpolation. In this work we propose a generalized framework for designing and analyzing context-aware hybrid FC/QWC commutativity relations. We use our framework to propose a noise-and-connectivity aware grouping strategy: Generalized backend-Aware pauLI Commutation (GALIC). We demonstrate how GALIC interpolates between FC and QWC, maintaining estimator accuracy in Hamiltonian estimation while lowering variance by an average of 20% compared to QWC. We also explore the design space of near-term quantum devices using the GALIC framework, specifically comparing device noise levels and connectivity. We find that error suppression has a more than 10 × larger impact on device-aware estimator variance than qubit connectivity with even larger correlation differences in estimator biases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extreme metrics from large ensembles: investigating the effects of ensemble size on their estimates

Abstract. We consider the problem of estimating the ensemble sizes required to characterize the forced component and the internal variability of a number of extreme metrics. While we exploit existing large ensembles, our perspective is that of a modeling center wanting to estimate a priori such sizes on the basis of an existing small ensemble (we assume the availability of only five members here). We therefore ask if such a small-size ensemble is sufficient to estimate accurately the population variance (i.e., the ensemble internal variability) and then apply a well-established formula that quantifies the expected error in the estimation of the population mean (i.e., the forced component) as a function of the sample size n, here taken to mean the ensemble size. We find that indeed we can anticipate errors in the estimation of the forced component for temperature and precipitation extremes as a function of n by plugging into the formula an estimate of the population variance derived on the basis of five members. For a range of spatial and temporal scales, forcing levels (we use simulations under Representative Concentration Pathway 8.5) and two models considered here as our proof of concept, it appears that an ensemble size of 20 or 25 members can provide estimates of the forced component for the extreme metrics considered that remain within small absolute and percentage errors. Additional members beyond 20 or 25 add only marginal precision to the estimate, and this remains true when statistical inference through extreme value analysis is used. We then ask about the ensemble size required to estimate the ensemble variance (a measure of internal variability) along the length of the simulation and – importantly – about the ensemble size required to detect significant changes in such variance along the simulation with increased external forcings. Using the F test, we find that estimates on the basis of only 5 or 10 ensemble members accurately represent the full ensemble variance even when the analysis is conducted at the grid-point scale. The detection of changes in the variance when comparing different times along the simulation, especially for the precipitation-based metrics, requires larger sizes but not larger than 15 or 20 members. While we recognize that there will always exist applications and metric definitions requiring larger statistical power and therefore ensemble sizes, our results suggest that for a wide range of analysis targets and scales an effective estimate of both forced component and internal variability can be achieved with sizes below 30 members. This invites consideration of the possibility of exploring additional sources of uncertainty, such as physics parameter settings, when designing ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

A Standardized Analysis Process Using Digital Image Correlation to Calculate In Situ Cladding Strain from Modified Burst Tests for Fuel Performance Code Validation

Historical data collection on nuclear fuel cladding materials has focused on generating a statistically significant amount of data to assess the material and its failure behavior. Furthermore, data generated to support material model and failure criteria development were previously posttest evaluations, so a large number of tests was required to gain new understanding. A way to expedite this process is to develop techniques capable of generating large, high-fidelity data sets from a single test with lower uncertainty or quantified uncertainty. One such example of this approach is Oak Ridge National Laboratory’s use of modified burst tests (MBTs) to analyze the mechanical behavior and failure conditions of cladding during a simulated reactivity-initiated accident (RIA). Each test incorporates digital image correlation (DIC) analysis techniques that are used to assess the accumulated strain in situ, as well as eventual cladding failure. This work has been fruitful in defining strain-to-failure conditions for materials like silicon carbide (SiC) fiber–reinforced/SiC matrix composite tubes (SiC/SiC), iron-chromium-aluminum (FeCrAl) alloy tubes, and chromium-coated Zircaloy-4 tubes. However, there are numerous DIC software available, including open-source and proprietary software. The different DIC software use various algorithms to process images and calculate displacement values. Using these different software and algorithms can lead to varying results, and perhaps larger-than-expected uncertainties. In the present study, previously published MBT data encompassing a variety of test conditions were reanalyzed with two different DIC software to assess the variance in the calculated strain results. The data consisted of SiC/SiC, FeCrAl, and chromium-coated Zircaloy-4 tubes. Plots of the calculated strains during the transient revealed good agreement between the two DIC software. The average root-mean-square errors between the two software was 0.20% strain, which is slightly larger than a previously reported error value for these tests. In conclusion, this variance in results is low enough that this analysis method can be used for code validation.

Reactivity-initiated accident↗

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗

Toward quantum Monte Carlo forces on heavier ions: Scaling properties

Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Z eff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Z eff , increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Zeff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Z e ff 6.5 ( 3 ) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Z e ff &#x2212; 2 , insensitive to model assumptions or the use of the “space warp” variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Z eff .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An increase in marine heatwaves without significant changes in surface ocean temperature variability

Marine heatwaves (MHWs)—extremely warm, persistent sea surface temperature (SST) anomalies causing substantial ecological and economic consequences—have increased worldwide in recent decades. Concurrent increases in global temperatures suggest that climate change impacted MHW occurrences, beyond random changes arising from natural internal variability. Moreover, the long-term SST warming trend was not constant but instead had more rapid warming in recent decades. Here we show that this nonlinear trend can—on its own—appear to increase SST variance and hence MHW frequency. Using a Linear Inverse Model to separate climate change contributions to SST means and internal variability, both in observations and CMIP6 historical simulations, we find that most MHW increases resulted from regional mean climate trends that alone increased the probability of SSTs exceeding a MHW threshold. Our results suggest the need to carefully attribute global warming-induced changes in climate extremes, which may not always reflect underlying changes in variability.

54 ENVIRONMENTAL SCIENCES↗