Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systematic error”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

DESI and other Dark Energy experiments in the era of neutrino mass measurements

Here we present Fisher matrix projections for future cosmological parameter measurements, including neutrino masses, Dark Energy, curvature, modified gravity, the inflationary perturbation spectrum, non-Gaussianity, and dark radiation. We focus on DESI and generally redshift surveys (BOSS, HETDEX, eBOSS, Euclid, and WFIRST), but also include CMB (Planck) and weak gravitational lensing (DES and LSST) constraints. The goal is to present a consistent set of projections, for concrete experiments, which are otherwise scattered throughout many papers and proposals. We include neutrino mass as a free parameter in most projections, as it will inevitably be relevant $-$ DESI and other experiments can measure the sum of neutrino masses to ~ 0.02 eV or better, while the minimum possible sum is 0.06 eV. We note that constraints on Dark Energy are significantly degraded by the presence of neutrino mass uncertainty, especially when using galaxy clustering only as a probe of the BAO distance scale (because this introduces additional uncertainty in the background evolution after the CMB epoch). Using broadband galaxy power becomes relatively more powerful, and bigger gains are achieved by combining lensing survey constraints with redshift survey constraints. Finally, we do not try to be especially innovative, e.g., with complex treatments of potential systematic errors $-$ these projections are intended as a straightforward baseline for comparison to more detailed analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A glimpse at quasar host galaxy far-UV emission using damped Lyα's as natural coronagraphs

In merger-driven models of massive galaxy evolution, the luminous quasar phase is expected to be accompanied by vigorous star formation in quasar host galaxies. In this paper, we use high column density damped Lyα (DLA) systems along quasar sight lines as natural coronagraphs to directly study the far-UV (FUV) radiation from the host galaxies of luminous background quasars. We have stacked the spectra of ∼2000 DLA systems (N {sub H} {sub I} > 10{sup 20.6} cm{sup –2}) with a median absorption redshift (z) = 2.6 selected from quasars observed in the SDSS-III Baryon Oscillation Spectroscopic Survey. We detect residual flux in the dark troughs of the composite DLA spectra. The level of this residual flux significantly exceeds systematic errors in the Sloan Digital Sky Survey fiber sky subtraction; furthermore, the residual flux is strongly correlated with the continuum luminosity of the background quasar, while uncorrelated with DLA column density or metallicity. We conclude that the flux could be associated with the average FUV radiation from the background quasar host galaxies (with medium redshift (z) = 3.1) that is not blocked by the intervening DLA. Assuming that all of the detected flux originates from quasar hosts, for the highest quasar luminosity bin ((L) = 2.5 × 10{sup 13} L {sub ☉}), the host galaxy has an FUV intensity of 1.5 ± 0.2 × 10{sup 40} erg s{sup –1} Å{sup –1}; this corresponds to an unobscured UV star formation rate of 9 M {sub ☉} yr{sup –1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Halo Boundary of Galaxy Clusters in the SDSS

Analytical models and simulations predict a rapid decline in the halo density profile associated with the transition from the “infalling” regime outside the halo to the “collapsed” regime within the halo. Using data from SDSS, we explore evidence for such a feature in the density profiles of galaxy clusters using several different approaches. We first estimate the steepening of the outer galaxy density profile around clusters, finding evidence for truncation of the halo profile. Next, we measure the galaxy density profile around clusters using two sets of galaxies selected on color. We find evidence of an abrupt change in galaxy colors that coincides with the location of the steepening of the density profile. Since galaxies that have completed orbits within the cluster are more likely to be quenched of star formation and thus appear redder, this abrupt change in galaxy color can be associated with the transition from single-stream to multi-stream regimes. We also use a standard model comparison approach to measure evidence for a “splashback”-like feature, but find that this approach is very sensitive to modeling assumptions. Finally, we perform measurements using an independent cluster catalog to test for potential systematic errors associated with cluster selection. We identify several avenues for future work: improved understanding of the small-scale galaxy profile, lensing measurements, identification of proxies for the halo accretion rate, and other tests. With upcoming data from the DES, KiDS, and HSC surveys, we can expect significant improvements in the study of halo boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

How to optimally constrain galaxy assembly bias: supplement projected correlation functions with count-in-cells statistics

Most models for the statistical connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than halo mass, a phenomenon known as galaxy assembly bias. And yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data that are analysed using traditional halo occupation models. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy-halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function w p ( r p ), the galaxy-galaxy lensing signal ΔΣ( r p ), the void probability function VPF( r ), the distributions of counts-in-cylinders P(NCIC), and counts-in-annuli P(NCIA), and the distribution of the ratio of counts in cylinders of different sizes P ( N 2 /N 5 ). We find that despite the frequent use of the combination w p ( r p ) + ΔΣ( r p ) in interpreting galaxy data, the count statistics, P ( N CIC ) and P ( N CIA ), are generally more efficient in constraining galaxy assembly bias when combined with w p ( r p ). Constraints based upon w p ( r p ) and ΔΣ(rp) share common degeneracy directions in the parameter space, while combinations of r p w p ( r p ) with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy-halo connection.

galaxies: evolution↗
DOE OSTI2017

Observation of the Crab Nebula with the HAWC Gamma-Ray Observatory

The Crab Nebula is the brightest TeV gamma-ray source in the sky and has been used for the past 25 years as a reference source in TeV astronomy, for calibration and verification of new TeV instruments. The High Altitude Water Cherenkov Observatory (HAWC), completed in early 2015, has been used to observe the Crab Nebula at high significance across nearly the full spectrum of energies to which HAWC is sensitive. HAWC is unique for its wide field of view, nearly 2 sr at any instant, and its high-energy reach, up to 100 TeV. HAWC’s sensitivity improves with the gamma-ray energy. Above ∼1 TeV the sensitivity is driven by the best background rejection and angular resolution ever achieved for a wide-field ground array. We present a time-integrated analysis of the Crab using 507 live days of HAWC data from 2014 November to 2016 June. The spectrum of the Crab is fit to a function of the form ϕ(E)=ϕ{sub 0}(E/E{sub 0}){sup −α−β⋅ln(E/E{sub 0})}. The data is well fitted with values of α = 2.63 ± 0.03, β = 0.15 ± 0.03, and log{sub 10}(ϕ{sub 0} cm{sup 2} s TeV)=−12.60±0.02 when E {sub 0} is fixed at 7 TeV and the fit applies between 1 and 37 TeV. Study of the systematic errors in this HAWC measurement is discussed and estimated to be ±50% in the photon flux between 1 and 37 TeV. Confirmation of the Crab flux serves to establish the HAWC instrument’s sensitivity for surveys of the sky. The HAWC all-sky survey will be the deepest survey of the northern sky ever conducted in the multi-TeV band.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The Aemulus Project. I. Numerical Simulations for Precision Cosmology

The rapidly growing statistical precision of galaxy surveys has led to a need for ever more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper, we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h(-1) Gpc)(3) simulations with mass resolution and force softening of 3.51 x 10(10) (Omega(m)/0.3) h(-1) M-circle dot and 20 h(-1) kpc, respectively, in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent cosmic microwave background, baryon acoustic oscillation, and Type Ia supernova results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) or to the sample variance of the statistic, whichever is larger, for halos with more than 500 (200) particles, respectively, and scales of r > 200 h(-1) kpc. in real space or k similar to 3 h Mpc(-1). in harmonic space for z <= 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k similar to 4 h Mpc(-1). We provide the halo catalogs and snapshots detailed in this work to the community at. https://AemulusProject.github.io.

large-scale structure of universe↗
DOE OSTI2022

Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies

The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H 0 , the σ 8 – S 8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0 σ tension between the Planck CMB estimate of the Hubble constant H 0 and the SH0ES collaboration measurements. After showing the H 0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ω m , and the amplitude or rate of the growth of structure ( σ 8 , f σ 8 ). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H 0 – S 8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Radio-optical reference frame link using the U.S. Naval observatory astrograph and deep CCD imaging

Between 1997 and 2004 several observing runs were conducted, mainly with the CTIO 0.9 m, to image International Celestial Reference Frame (ICRF) counterparts (mostly QSOs) in order to determine accurate optical positions. Contemporary to these deep CCD images, the same fields were observed with the U.S. Naval Observatory astrograph in the same bandpass. They provide accurate positions on the Hipparcos/Tycho-2 system for stars in the 10-16 mag range used as reference stars for the deep CCD imaging data. Here we present final optical position results of 413 sources based on reference stars obtained by dedicated astrograph observations that were reduced following two different procedures. These optical positions are compared to radio very long baseline interferometry positions. The current optical system is not perfectly aligned to the ICRF radio system with rigid body rotation angles of 3-5 mas (= 3σ level) found between them for all three axes. Furthermore, statistically, the optical-radio position differences are found to exceed the total, combined, known errors in the observations. Systematic errors in the optical reference star positions and physical offsets between the centers of optical and radio emissions are both identified as likely causes. A detrimental, astrophysical, random noise component is postulated to be on about the 10 mas level. If confirmed by future observations, this could severely limit the Gaia to ICRF reference frame alignment accuracy to an error of about 0.5 mas per coordinate axis with the current number of sources envisioned to provide the link. A list of 36 ICRF sources without the detection of an optical counterpart to a limiting magnitude of about R = 22 is provided as well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Characterizing non-Markovian and coherent errors in quantum simulation

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

( Hubble Space Telescope Hubble Space Telescope) proper motion (HSTPROMO) catalogs of galactic globular clusters. IV. Kinematic profiles and average masses of blue straggler stars

We make use of the Hubble Space Telescope proper-motion catalogs derived by Bellini et al. to produce the first radial velocity dispersion profiles σ(R) for blue straggler stars (BSSs) in Galactic globular clusters (GCs), as well as the first dynamical estimates for the average mass of the entire BSS population. We show that BSSs typically have lower velocity dispersions than stars with mass equal to the main-sequence turnoff mass, as one would expect for a more massive population of stars. Since GCs are expected to experience some degree of energy equipartition, we use the relation σ∝M{sup −η}, where η is related to the degree of energy equipartition, along with our velocity dispersion profiles to estimate BSS masses. We estimate η as a function of cluster relaxation from recent Monte Carlo cluster simulations by Bianchini et al. and then derive an average mass ratio M{sub BSS}/M{sub MSTO}=1.50±0.14 and an average mass M{sub BSS}=1.22±0.12 M {sub ⊙} from 598 BSSs across 19 GCs. The final error bars include any systematic errors that are random between different clusters, but not any potential biases inherent to our methodology. Our results are in good agreement with the average mass of M{sub BSS}=1.22±0.06 M {sub ⊙} for the 35 BSSs in Galactic GCs in the literature with properties that have allowed individual mass determination.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Astrometry of OH/IR Stars Using 1612 MHz Hydroxyl Masers. I. Annual Parallaxes of WX Psc and OH138.0+7.2

We report on the measurement of the trigonometric parallaxes of 1612 MHz hydroxyl masers around two asymptotic giant branch stars, WX Psc and OH 138.0+7.2, using the NRAO Very Long Baseline Array with in-beam phase referencing calibration. We obtain a 3σ upper limit of ≤5.3 mas on the parallax of WX Psc, corresponding to a lower limit distance estimate of ≳190 pc. The obtained parallax of OH 138.0+7.2 is 0.52 ± 0.09 mas (±18%), corresponding to a distance of 1.9{sub −0.3}{sup +0.4} kpc, making this the first hydroxyl maser parallax below one milliarcsecond. We also introduce a new method of error analysis for detecting systematic errors in the astrometry. Finally, we compare our trigonometric distances to published phase-lag distances toward these stars and find a good agreement between the two methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A critical assessment of stellar mass measurement methods

This is the second paper in a series aimed at investigating the main sources of uncertainty in measuring the observable parameters in galaxies from their spectral energy distributions (SEDs). In the first paper we presented a detailed account of the photometric redshift measurements and an error analysis of this process. In this paper we perform a comprehensive study of the main sources of random and systematic error in stellar mass estimates for galaxies, and their relative contributions to the associated error budget. Since there is no prior knowledge of the stellar mass of galaxies (unlike their photometric redshifts), we use mock galaxy catalogs with simulated multi-waveband photometry and known redshift, stellar mass, age and extinction for individual galaxies. The multi-waveband photometry for the simulated galaxies were generated in 13 filters spanning from U-band to mid-infrared wavelengths. Given different parameters affecting stellar mass measurement (photometric signal-to-noise ratios (S/N), SED fitting errors and systematic effects), the inherent degeneracies and correlated errors, we formulated different simulated galaxy catalogs to quantify these effects individually. For comparison, we also generated catalogs based on observed photometric data of real galaxies in the Great Observatories Origins Deep Survey-South field, spanning the same passbands. The simulated and observed catalogs were provided to a number of teams within the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey collaboration to estimate the stellar masses for individual galaxies. A total of 11 teams participated, with different combinations of stellar mass measurement codes/methods, population synthesis models, star formation histories, extinction and age. For each simulated galaxy, the differences between the input stellar masses, M{sub input}, and those estimated by each team, M{sub est}, is defined as Δlog(M)≡log(M{sub estimated})−log(M{sub input}), and used to identify the most fundamental parameters affecting stellar mass estimate in galaxies, with the following results. (1) No significant bias in Δ log(M) was found among different codes, with all having comparable scatter (σ(Δlog(M))=0.136 dex). The estimated stellar mass values are seriously affected by low photometric S/N, with the rms scatter increasing for galaxies with H{sub AB}>26 mag; (2) A source of error contributing to the scatter in Δ log(M) is found to be due to photometric uncertainties (0.136 dex) and low resolution in age and extinction grids when generating the SED templates; (3) The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy (σ(Δlog(M))=0.142 dex); (4) The Δ log(M) values are strongly correlated with deviations in age (defined as the difference between the estimated and expected values), with a weaker correlation with extinction; (5) The rms scatter in the estimated stellar masses due to free parameters (after fixing redshifts and initial mass function) are quantified and found to be σ(Δlog(M))=0.110 dex; (6) Using the observed data, we studied the sensitivity of stellar masses to both the population synthesis codes and inclusion of nebular emission lines and found them to affect the stellar mass by 0.2 and 0.3 dex respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Spectroscopic Needs for Imaging Dark Energy Experiments

Ongoing and near-future imaging-based dark energy experiments are critically dependent upon photometric redshifts (a.k.a. photo-z’s): i.e., estimates of the redshifts of objects based only on flux information obtained through broad filters. Higher-quality, lower-scatter photo-z’s will result in smaller random errors on cosmological parameters; while systematic errors in photometric redshift estimates, if not constrained, may dominate all other uncertainties from these experiments. The desired optimization and calibration is dependent upon spectroscopic measurements for secure redshift information; this is the key application of galaxy spectroscopy for imaging-based dark energy experiments. Hence, to achieve their full potential, imaging-based experiments will require large sets of objects with spectroscopically-determined redshifts, for two purposes: Training: Objects with known redshift are needed to map out the relationship between object color and z (or, equivalently, to determine empirically-calibrated templates describing the rest-frame spectra of the full range of galaxies, which may be used to predict the color-z relation). The ultimate goal of training is to minimize each moment of the distribution of differences between photometric redshift estimates and the true redshifts of objects, making the relationship between them as tight as possible. The larger and more complete our “training set” of spectroscopic redshifts is, the smaller the RMS photo-z errors should be, increasing the constraining power of imaging experiments; Requirements: Spectroscopic redshift measurements for ~30,000 objects over >~15 widely-separated regions, each at least ~20 arcmin in diameter, and reaching the faintest objects used in a given experiment, will likely be necessary if photometric redshifts are to be trained and calibrated with conventional techniques. Larger, more complete samples (i.e., with longer exposure times) can improve photo-z algorithms and reduce scatter further, enhancing the science return from planned experiments greatly (increasing the Dark Energy Task Force figure of merit by up to ~50%); Options: This spectroscopy will most efficiently be done by covering as much of the optical and near-infrared spectrum as possible at modestly high spectral resolution (λ/Δλ > ~3000), while maximizing the telescope collecting area, field of view on the sky, and multiplexing of simultaneous spectra. The most efficient instrument for this would likely be either the proposed GMACS/MANIFEST spectrograph for the Giant Magellan Telescope or the OPTIMOS spectrograph for the European Extremely Large Telescope, depending on actual properties when built. The PFS spectrograph at Subaru would be next best and available considerably earlier, c. 2018; the proposed ngCFHT and SSST telescopes would have similar capabilities but start later. Other key options, in order of increasing total time required, are the WFOS spectrograph at TMT, MOONS at the VLT, and DESI at the Mayall 4 m telescope (or the similar 4MOST and WEAVE projects); of these, only DESI, MOONS, and PFS are expected to be available before 2020. Table 2-3 of this white paper summarizes the observation time required at each facility for strawman training samples. To attain secure redshift measurements for a high fraction of targeted objects and cover the full redshift span of future experiments, additional near-infrared spectroscopy will also be required; this is best done from space, particularly with WFIRST-2.4 and JWST; Calibration: The first several moments of redshift distributions (the mean, RMS redshift dispersion, etc.), must be known to high accuracy for cosmological constraints not to be systematics-dominated (equivalently, the moments of the distribution of differences between photometric and true redshifts could be determined instead). The ultimate goal of calibration is to characterize these moments for every subsample used in analyses - i.e., to minimize the uncertainty in their mean redshift, RMS dispersion, etc. – rather than to make the moments themselves small. Calibration may be done with the same spectroscopic dataset used for training if that dataset is extremely high in redshift completeness (i.e., no populations of galaxies to be used in analyses are systematically missed). Accurate photo-z calibration is necessary for all imaging experiments; Requirements: If extremely low levels of systematic incompleteness (<~0.1%) are attained in training samples, the same datasets described above should be sufficient for calibration. However, existing deep spectroscopic surveys have failed to yield secure redshifts for 30–60% of targets, so that would require very large improvements over past experience. This incompleteness would be a limiting factor for training, but catastrophic for calibration. If <~0.1% incompleteness is not attainable, the best known option for calibration of photometric redshifts is to utilize cross-correlation statistics in some form. The most direct method for this uses cross-correlations between positions on the sky of bright objects of known spectroscopic redshift with the sample of objects that we wish to calibrate the redshift distribution for, measured as a function of spectroscopic z. For such a calibration, redshifts of ~100,000 objects over at least several hundred square degrees, spanning the full redshift range of the samples used for dark energy, would be necessary; and Options: The proposed BAO experiment eBOSS would provide sufficient spectroscopy for basic calibrations, particularly for ongoing and near-future imaging experiments. The planned DESI experiment would provide excellent calibration with redundant cross-checks, but will start after the conclusion of some imaging projects. An extension of DESI to the Southern hemisphere would provide the best possible calibration from cross-correlation methods for DES and LSST. We thus anticipate that our two primary needs for spectroscopy – training and calibration of photometric redshifts – will require two separate solutions. For ongoing and future projects to reach their full potential, new spectroscopic samples of faint objects will be needed for training; those new samples may be suitable for calibration, but the latter possibility is uncertain. In contrast, wide-area samples of bright objects are poorly suited for training, but can provide high-precision calibrations via cross-correlation techniques. Additional training/calibration redshifts and/or host galaxy spectroscopy would enhance the use of supernovae and galaxy clusters for cosmology. We also summarize additional work on photometric redshift techniques that will be needed to prepare for data from ongoing and future dark energy experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

$B \to \pi \ell \nu$ and $B_s \to K \ell \nu$ form factors and $|V_{ub}|$ from 2+1-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

We calculate the form factors for B → πℓν and B s → Kℓν decay in dynamical lattice quantum chromodynamics (QCD) using domain-wall light quarks and relativistic b-quarks. We use the (2+1)-flavor gauge-field ensembles generated by the RBC and UKQCD collaborations with the domain-wall fermion action and Iwasaki gauge action. For the b-quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze data at two lattice spacings of a ≈ 0.11, 0.086 fm with unitary pion masses as light as M π ≈ 290 MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) “hard-pion” chiral perturbation theory for heavy-light meson form factors. We provide complete systematic error budgets for the vector and scalar form factors f + (q 2 ) and f 0 (q 2 ) for both B → πℓν and B s → Kℓν at three momenta that span the q 2 range accessible in our numerical simulations. Next we extrapolate these results to q 2 = 0 using a model-independent z-parametrization based on analyticity and unitarity. We present our final results for f + (q 2 ) and f 0 (q 2 )as the coefficients of the series in z and the matrix of correlations between them; this provides a parametrization of the form factors valid over the entire allowed kinematic range. Our results agree with other three-flavor lattice-QCD determinations using staggered light quarks, and have comparable precision, thereby providing important independent cross-checks. Both B → πℓν and B s → Kℓν decays enable determinations of the Cabibbo-Kobayashi-Maskawa matrix element |V ub |. Furthermore, we perform a combined z-fit of our numerical B → πℓν form-factor data with the experimental measurements of the branching fraction from BABAR and Belle leaving the relative normalization as a free parameter; we obtain |V ub | = 3.61(32)×10 -3 , where the error includes statistical and all systematic uncertainties. The same approach can be applied to the decay B s → Kℓν to provide an alternative determination of |V ub | once the process has been measured experimentally. In anticipation of future experimental measurements, we make predictions for B → πℓν and B s → Kℓν differential branching fractions and forward-backward asymmetries in the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Towards a systematic assessment of errors in diffusion Monte Carlo calculations of semiconductors: Case study of zinc selenide and zinc oxide

The fixed node diffusion Monte Carlo (DMC) method has attracted interest in recent years as a way to calculate properties of solid materials with high accuracy. However, the framework for the calculation of properties such as total energies, atomization energies, and excited state energies is not yet fully established. Several outstanding questions remain as to the effect of pseudopotentials, the magnitude of the fixed node error, and the size of supercell finite size effects. Here, we consider in detail the semiconductors ZnSe and ZnO and carry out systematic studies to assess the magnitude of the energy differences arising from controlled and uncontrolled approximations in DMC. The former include time step errors and supercell finite size effects for ground and optically excited states, and the latter include pseudopotentials, the pseudopotential localization approximation, and the fixed node approximation. We find that for these compounds, the errors can be controlled to good precision using modern computational resources and that quantum Monte Carlo calculations using Dirac-Fock pseudopotentials can offer good estimates of both cohesive energy and the gap of these systems. We do however observe differences in calculated optical gaps that arise when different pseudopotentials are used.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

B -meson decay constants from 2 + 1 -flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

We calculate the B-meson decay constants f B , f B s , and their ratio in unquenched lattice QCD using domain-wall light quarks and relativistic b quarks. We use gauge-field ensembles generated by the RBC and UKQCD collaborations using the domain-wall fermion action and Iwasaki gauge action with three flavors of light dynamical quarks. We analyze data at two lattice spacings of a ≈ 0.11, 0.086 fm with unitary pion masses as light as M π ≈ 290 MeV; this enables us to control the extrapolation to the physical light-quark masses and continuum. For the b quarks we use the anisotropic clover action with the relativistic heavy-quark interpretation, such that discretization errors from the heavy-quark action are of the same size as from the light-quark sector. We renormalize the lattice heavy-light axial-vector current using a mostly nonperturbative method in which we compute the bulk of the matching factor nonperturbatively, with a small correction, that is close to unity, in lattice perturbation theory. We also improve the lattice heavy-light current through O(α s a). We extrapolate our results to the physical light-quark masses and continuum using SU(2) heavy-meson chiral perturbation theory, and provide a complete systematic error budget. We obtain f B 0 = 199.5(12.6) MeV, f B + =195.6(14.9) MeV, f B s =235.4(12.2) MeV, f B s /f B 0 =1.197(50), and f B s /f B + =1.223(71), where the errors are statistical and total systematic added in quadrature. Finally, these results are in good agreement with other published results and provide an important independent cross-check of other three-flavor determinations of B-meson decay constants using staggered light quarks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Angular systematics-free cosmological analysis of galaxy clustering in configuration space

ABSTRACT Galaxy redshift surveys are subject to incompleteness and inhomogeneous sampling due to the various constraints inherent to spectroscopic observations. This can introduce systematic errors on the summary statistics of interest, which need to be mitigated in cosmological analysis to achieve high accuracy. Standard practices involve applying weighting schemes based on completeness estimates across the survey footprint, possibly supplemented with additional weighting schemes accounting for density-dependent effects. In this work, we concentrate on pure angular systematics and describe an alternative approach consisting in analysing the galaxy two-point correlation function where angular modes are nulled. By construction, this procedure removes all possible known and unknown sources of angular observational systematics, but also part of the cosmological signal. We use a modified Landy–Szalay estimator for the two-point correlation function that relies on an additional random catalogue where angular positions are randomly drawn from the galaxy catalogue, and provide an analytical model to describe this modified statistic. We test the model by performing an analysis of the full anisotropic clustering in mock catalogues of luminous red and emission-line galaxies at 0.43 < z < 1.1. We find that the model fully accounts for the modified correlation function in redshift space, without introducing new nuisance parameters. The derived cosmological parameters from the analysis of baryon acoustic oscillations and redshift-space distortions display slightly larger statistical uncertainties, mostly for the growth rate of structure parameter fσ8 that exhibits a $50{{\ \rm per\ cent}}$ statistical error increase, but free from angular systematic error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Analytic Chromaticity Formulas for the Muon g -2 experiment at fermilab

The Muon g-2 Experiment (E989) at Fermilab measured the muon anomalous magnetic moment aμ with unprecedented precision of 127 ppb using a storage ring. The experiment pushed the limits in terms of accuracy and systematic errors. Achieving the 78-ppb total systematic uncertainty required precise beam dynamics modeling and corrections for effects on the measured muon precession frequency. Here, we derived the first analytic aberration formulas up to the second order for the muon g-2 storage ring’s combined-function electrostatic quadrupoles (superimposed magnetic dipole and electric quadrupole fields) using an order-by-order perturbation method. From these, we obtained the exact chromaticity formulas for three ring models of different granularity and validated them against numerical calculations using COSY INFINITY, achieving analytic-numerical agreement to ⁠$\mathscr{O}(10^{-10})$. This work resolved discrepancies between previous approximate derivations and provided essential beam dynamics results for Runs 4-6 analyses. We also calculated nonlinear chromaticities up to ninth order. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous muon g-2 experiment at Brookhaven National Laboratory, with the final result announced in June 2025.

aberrations↗