Multi-armed bandit with sub-exponential rewards
Not provided.
SEARCH · Search NASA
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.
Not provided.
We present measurements of the cosmic microwave background (CMB) lensing potential using the final Planck 2018 temperature and polarization data. Using polarization maps filtered to account for the noise anisotropy, we increase the significance of the detection of lensing in the polarization maps from 5σ to 9σ. Combined with temperature, lensing is detected at 40σ. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles 8 ≤ L ≤ 400 (extending the range to lower L compared to 2015), which we use to constrain cosmological parameters. We find good consistency between lensing constraints and the results from the Planck CMB power spectra within the ΛCDM model. Combined with baryon density and other weak priors, the lensing analysis alone constrains σ 8 Ωm 0.25 = 0.589 ± 0.020 (1σ errors). Also combining with baryon acoustic oscillation data, we find tight individual parameter constraints, σ8 = 0.811 ± 0.019, H0 = 67.9-1.3+1.2 km s -1 Mpc -1 , and Ω m = 0.303 -0.018 +0.016 . Combining with Planck CMB power spectrum data, we measure σ 8 to better than 1% precision, finding σ 8 = 0.811 ± 0.006. CMB lensing reconstruction data are complementary to galaxy lensing data at lower redshift, having a different degeneracy direction in σ 8 - Ω m space; we find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using the Planck cosmic infrared background (CIB) maps as an additional tracer of high-redshift matter, we make a combined Planck-only estimate of the lensing potential over 60% of the sky with considerably more small-scale signal. We additionally demonstrate delensing of the Planck power spectra using the joint and individual lensing potential estimates, detecting a maximum removal of 40% of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance.
We study the asymmetric primordial fluctuations in a model of inflation in which translational invariance is broken by a domain wall. We calculate the corrections to the power spectrum of curvature perturbations; they are anisotropic and contain dipole, quadrupole, and higher multipoles with non-trivial scale-dependent amplitudes. Inspired by observations of these multipole asymmetries in terms of two-point correlations and variance in real space, we demonstrate that this model can explain the observed anomalous power asymmetry of the cosmic microwave background (CMB) sky, including its characteristic feature that the dipole dominates over higher multipoles. We test the viability of the model and place approximate constraints on its parameters by using observational values of dipole, quadrupole, and octopole amplitudes of the asymmetry measured by a local-variance estimator. We find that a configuration of the model in which the CMB sphere does not intersect the domain wall during inflation provides a good fit to the data. We further derive analytic expressions for the corrections to the CMB temperature covariance matrix, or angular power spectra, which can be used in future statistical analysis of the model in spherical harmonic space.
Chemical evolution models are powerful tools for interpreting stellar abundance surveys and understanding galaxy evolution. However, their predictions depend heavily on the treatment of inflow, outflow, star formation efficiency (SFE), the stellar initial mass function, the SN Ia delay time distribution, stellar yields, and stellar population mixing. Using flexCE, a flexible one-zone chemical evolution code, we investigate the effects of and trade-offs between parameters. Two critical parameters are SFE and the outflow mass-loading parameter, which shift the knee in [O/Fe]–[Fe/H] and the equilibrium abundances that the simulations asymptotically approach, respectively. One-zone models with simple star formation histories follow narrow tracks in [O/Fe]–[Fe/H] unlike the observed bimodality (separate high- α and low- α sequences) in this plane. A mix of one-zone models with inflow timescale and outflow mass-loading parameter variations, motivated by the inside-out galaxy formation scenario with radial mixing, reproduces the two sequences better than a one-zone model with two infall epochs. We present [X/Fe]–[Fe/H] tracks for 20 elements assuming three different supernova yield models and find some significant discrepancies with solar neighborhood observations, especially for elements with strongly metallicity-dependent yields. We apply principal component abundance analysis to the simulations and existing data to reveal the main correlations among abundances and quantify their contributions to variation in abundance space. For the stellar population mixing scenario, the abundances of α -elements and elements with metallicity-dependent yields dominate the first and second principal components, respectively, and collectively explain 99% of the variance in the model. flexCE is a python package available at https://github.com/bretthandrews/flexCE.
This work was carried out as part of the Heliostat Consortium (HelioCon) Field Deployment subtask with the aim to develop a reduced order model framework for correlating wind speed and pointing deviation of a heliostat facet. There are only sparse field measurements of heliostat pointing deviations and accompanying wind conditions published in the literature. Heliostat test standards, such as IEC 62862-4-3, propose a suite of tests including laser pointing repeatability at wind speeds below 4 m/s, and provide technical requirements for heliostat slope and tracking deviations in coarse average wind speed bins of 4 m/s, 6 m/s, and 8 m/s. In addressing the gap of the variation of heliostat pointing deviation with wind speed, field measurements of laser pointing on a grid target and wind conditions were analyzed in this study at the Third-Party Metrology Platform at the National Laboratory of the Rockies (NLR) Flatirons Campus. Horizontal pointing deviations were found to follow a logarithmic relationship with peak wind speed, whereas vertical pointing deviations follow an exponential relationship with peak wind speed. Both horizontal and vertical pointing deviations also follow a second order polynomial relationship, as expected from the proportionality of elastic loads and deformations with the square of wind speed. The results indicate that heliostat facet pointing deviations in the vertical direction increase at a faster rate than in the horizontal direction with increasing wind speed over the tested range, however these are dependent on the heliostat structural design. Next steps are recommended for additional field measurements to confirm a linear relationship of pointing deviation with applied moment on a heliostat facet, and to distinguish between gravity-induced and wind-induced pointing deviations at different elevation angles. The derived correlations in the preliminary analysis in this report serve as a case study for heliostat developers and plant operators to estimate the wind-induced pointing deviations and their variation with peak gust wind speed. Next steps in future work would recommend higher resolution and longer duration datasets for different elevation angles and wind directions to reduce uncertainties and variance of collected laser beam spot data and their correlations with bin-averaged wind speed.
One of the main goals of investigations using present and future giant extensive air shower (EAS) arrays is the mass composition of ultra-high energy cosmic rays (UHECRs). A new approach to the problem is presented, combining the analysis of arrival directions with the statistical test of the paired EAS samples. One of the ideas of the method is to search for possible correlations between UHECR masses and their separate sources; for instance, if there are two sources in different areas of the celestial sphere injecting different nuclei, but the fluxes are comparable so that arrival directions are isotropic, then the aim is to reveal a difference in the mass composition of cosmic-ray fluxes. The method is based on a non-parametric statistical test-the Wilcoxon signed-rank routine-which does not depend on the populations fitting any parameterized distributions. Two particular algorithms are proposed: first, using measurements of the depth of the EAS maximum position in the atmosphere; and second, relying on the age variance of air showers initiated by different primary particles. The formulated method is applied to the Yakutsk array data, in order to demonstrate the possibility of searching for a difference in average mass composition between the two UHECR sets, arriving particularly from the supergalactic plane and a complementary region.
The resonance mechanism in the initial of wind-wave generation proposed by Phillips is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.
The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterizing remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way’s accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution lambda-cold dark matter magnetohydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum.
First-principles molecular dynamics (FPMD) simulations based on density functional theory are becoming increasingly popular for the description of liquids. In view of the high computational cost of these simulations, the choice of an appropriate equilibration protocol is critical. In this study, we assess two methods of estimation of equilibration times using a large dataset of first-principles molecular dynamics simulations of water. The Gelman-Rubin potential scale reduction factor [A. Gelman and D. B. Rubin, Stat. Sci. 7, 457 (1992)] and the marginal standard error rule heuristic proposed by White [Simulation 69, 323 (1997)] are evaluated on a set of 32 independent 64-molecule simulations of 58 ps each, amounting to a combined cumulative time of 1.85 ns. The availability of multiple independent simulations also allows for an estimation of the variance of averaged quantities, both within MD runs and between runs. We analyze atomic trajectories, focusing on correlations of the Kohn-Sham energy, pair correlation functions, number of hydrogen bonds, and diffusion coefficient. The observed variability across samples provides a measure of the uncertainty associated with these quantities, thus facilitating meaningful comparisons of different approximations used in the simulations. We find that the computed diffusion coefficient and average number of hydrogen bonds are affected by a significant uncertainty in spite of the large size of the dataset used. Lastly, a comparison with classical simulations using the TIP4P/2005 model confirms that the variability of the diffusivity is also observed after long equilibration times. Complete atomic trajectories and simulation output files are available online for further analysis.
Abstract Much of the stochastic groundwater literature is devoted to the analysis of flow and transport in Gaussian or multi‐Gaussian log hydraulic conductivity (or transmissivity) fields, ( x being a position vector), characterized by one or (less frequently) a multiplicity of spatial correlation scales. Yet Y and many other variables and their (spatial or temporal) increments, Δ Y , are known to be generally non‐Gaussian. One common manifestation of non‐Gaussianity is that whereas frequency distributions of Y often exhibit mild peaks and light tails, those of increments are generally symmetric with peaks that grow sharper, and tails that become heavier, as separation scale or lag between pairs of Y values decreases. A statistical model that captures these disparate, scale‐dependent distributions of Y and in a unified and consistent manner has been recently proposed by us. This new “generalized sub‐Gaussian (GSG)” model has the form where is (generally, but not necessarily) a multiscale Gaussian random field and is a nonnegative subordinator independent of G . The purpose of this paper is to explore analytically, in an elementary manner, lead‐order effects that non‐Gaussian heterogeneity described by the GSG model have on the stochastic description of flow and transport. Recognizing that perturbation expansion of hydraulic conductivity diverges when Y is sub‐Gaussian, we render the expansion convergent by truncating Y 's domain of definition. We then demonstrate theoretically and illustrate by way of numerical examples that, as the domain of truncation expands, (a) the variance of truncated Y (denoted by Y t ) approaches that of Y and (b) the pdf (and thereby moments) of Y t increments approach those of Y increments and, as a consequence, the variogram of approaches that of . This in turn guarantees that perturbing to second order in (the standard deviation of Y t ) yields results which approach those we obtain upon perturbing to second order in even as the corresponding series diverges. Our analysis is rendered mathematically tractable by considering mean‐uniform steady state flow in an unbounded, two‐dimensional domain of mildly heterogeneous Y with a single‐scale function G having an isotropic exponential covariance. Results consist of expressions for (a) lead‐order autocovariance and cross‐covariance functions of hydraulic head, velocity, and advective particle displacement and (b) analogues of preasymptotic as well as asymptotic Fickian dispersion coefficients. We compare these theoretically and graphically with corresponding expressions developed in the literature for Gaussian Y . We find the former to differ from the latter by a factor k = ( denoting ensemble expectation) and the GSG covariance of longitudinal velocity to contain an additional nugget term depending on this same factor. In the limit as Y becomes Gaussian, k reduces to one and the nugget term drops out.
Abstract The current study examined the associations between peer victimization, body image shame, self‐criticism, self‐reassurance, depressive symptoms and eating psychopathology in 609 female adolescents. Correlational analyses showed that being the victim of bullying was positively associated with body image shame, self‐criticism, with low self‐reassurance, depressive symptoms and eating psychopathology. A path analysis indicated that victimization experiences were associated with increased depressive symptoms partially through increased levels of body image shame, and a severe form of self‐criticism – hated self. Body image shame and hated‐self self‐criticism fully mediated the association between victimization experiences eating psychopathology. The tested model accounted for a total of 51% of depressive symptoms variance and for 52% of eating psychopathology variance. These findings may have important intervention and prevention implications, by suggesting that bullying experiences fuel body image shame and consequent self‐directed hostility and anger, which, in turn, predict increased depressive symptomatology and eating psychopathology in female adolescents.
Forthcoming radio continuum surveys will cover large volumes of the observable Universe and will reach to high redshifts, making them potentially powerful probes of dark energy, modified gravity and non-Gaussianity. We consider the continuum surveys with LOFAR, WSRT and ASKAP, and examples of continuum surveys with the SKA. We extend recent work on these surveys by including redshift space distortions and lensing convergence in the radio source auto-correlation. In addition we compute the general relativistic (GR) corrections to the angular power spectrum. These GR corrections to the standard Newtonian analysis of the power spectrum become significant on scales near and beyond the Hubble scale at each redshift. We find that the GR corrections are at most percent-level in LOFAR, WODAN and EMU surveys, but they can produce O(10%) changes for high enough sensitivity SKA continuum surveys. The signal is however dominated by cosmic variance, and multiple-tracer techniques will be needed to overcome this problem. The GR corrections are suppressed in continuum surveys because of the integration over redshift — we expect that GR corrections will be enhanced for future SKA HI surveys in which the source redshifts will be known. We also provide predictions for the angular power spectra in the case where the primordial perturbations have local non-Gaussianity. We find that non-Gaussianity dominates over GR corrections, and rises above cosmic variance when f{sub NL}∼>5 for SKA continuum surveys.
A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.
The Nuclear Fuels Storage and Transportation Planning Project and the Used Fuel Disposition Campaign, both under the US Department of Energy (DOE) Office of Nuclear Energy, require the ability to predict k{sub eff}, heat generation and temperature distributions, and external dose rates from any combination of commercial spent nuclear fuel assemblies loaded into any type of licensed cask. Work in the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory has begun to develop the UNF-ST[and]DARDS system (Used Nuclear Fuel Storage, Transportation and Disposal Analysis Resource and Data System) that uses many of the nuclear analysis codes from the SCALE package, a comprehensive database of used nuclear fuel inventories from utilities, and a library of model template files to produce complete input files for the various types of UNF computational analyses. UNF-ST[and]DARDS uses the template files and input parameters that are applicable to a safety analysis to assemble and execute a complete input file, collect the output files and relevant information for subsequent calculations, and illustrate calculation results as a function of decay time and other relevant analysis parameters. This paper will focus on the shielding aspects and the ability to predict external dose rates from as-loaded casks. Realistic dose rate values based on actual canister fuel contents may be used in developing appropriate plans, procedures, and methods to keep personnel radiation exposures as low as reasonably achievable during cask repackaging and transportation operations, as well as in performing total system analyses and transportation risk evaluations. The MAVRIC shielding sequence of the SCALE package was created for the US Nuclear Regulatory Commission for analyzing dose rates outside of storage, transfer, or transportation casks loaded with spent nuclear fuel. MAVRIC contains automated variance reduction capabilities based on the CADIS and FW-CADIS methods that are effective for optimizing deep penetration Monte Carlo shielding problems for achieving low relative uncertainties in short computing times. MAVRIC uses the Denovo discrete ordinates code to compute importance maps and biased source terms for the Monaco fixed-source Monte Carlo code, which can use multi-group or continuous-energy cross section data. ORIGEN (Oak Ridge Isotope Generation and Depletion) is the SCALE computer code for calculating time-dependent isotopic concentrations during irradiation and decay for more than 2,200 nuclides that can be produced by nuclear fuel irradiation and activation. A series of SCALE enhancements was implemented to enable automatic decay heat and radiation source term generation and shielding safety analysis of as-loaded casks. The enhancements to SCALE made under this project will be a part of the SCALE 6.2 release expected in fall 2014. (authors)
Cosmic polarization rotation, which may result from parity-violating new physics or the presence of primordial magnetic fields, converts E-mode polarization of the cosmic microwave background (CMB) into B-mode polarization. Anisotropic cosmic polarization rotation leads to statistical anisotropy in CMB polarization and can be reconstructed with quadratic estimator techniques similar to those designed for gravitational lensing of the CMB. At the sensitivity of upcoming CMB surveys, lensing-induced B-mode polarization will act as a limiting factor in the search for anisotropic cosmic polarization rotation, meaning that an analysis which incorporates some form of delensing will be required to improve constraints on the effect with future surveys. In this study we extend the ResUNet-CMB convolutional neural network to reconstruct anisotropic cosmic polarization rotation in the presence of gravitational lensing and patchy reionization, and we show that the network simultaneously reconstructs all three effects with variance that is lower than that from the standard quadratic estimator nearly matching the performance of an iterative reconstruction method.
By applying the Wannier-based extended Kugel-Khomskii model, we carry out first-principles calculations and electronic structure analysis to understand the spin-phonon coupling effect in transition-metal perovskites. In this work, we demonstrate the successful application of our approach to SrMnO 3 and BiFeO 3 . We show that both the electron orbitals under crystal-field splitting and the electronic configuration should be taken into account in order to understand the large variances of spin-phonon coupling effects among various phonon modes as well as in different materials.
It is crucial that the U.S. electrical grid be reliable, due to its critical role in the economy and our national defense/security. Many common tasks and industrial operations rely on proper electrical power distribution. In designing new bus systems or evaluating existing ones, the type of configuration employed affects the level of risk, reliability, maintenance, and cost involved. The present work utilizes Systems Analysis Programs for Hands-on Integrated Reliability Evaluations (SAPHIRE) to perform a probabilistic risk assessment (PRA)-based sensitivity study on substation reliability in regard to different bus configurations. This study examines single bus, main and transfer, breaker and a half, double bus/double breaker, and ring bus configurations applied to busses with different numbers of lines. The case variances, methods, and rankings presented herein can help guide bus configurations for electrical grid design/analysis in the future.
In the field of nuclear safeguards, passive neutron multiplicity counting (PNMC) is a method typically employed in non-destructive assay (NDA) of special nuclear material (SNM) for nonproliferation, verification and accountability purposes. PNMC is generally performed using a well-type thermal neutron counter and relies on the detection of correlated pairs or higher order multiplets of neutrons emitted by an assayed item. To assay SNM, a set of parameters for a given well-counter is required to link the measured multiplicity rates to the assayed item properties. Detection efficiency, die-away time, gate utilization factors (tightly connected to die-away time) as well as optimum gate width setting are among the key parameters. These parameters along with the underlying model assumptions directly affect the accuracy of the SNM assay. In this paper we examine the role of gate utilization factors and the single exponential die-away time assumption and their impact on the measurements for a range of plutonium materials. In addition, we examine the importance of item-optimized coincidence gate width setting as opposed to using a universal gate width value. Finally, the traditional PNMC based on multiplicity shift register electronics is extended to Feynman-type analysis and application of this approach to Pu mass assay is demonstrated.