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

The Simons Observatory: gain, bandpass and polarization-angle calibration requirements for B-mode searches

We quantify the calibration requirements for systematic uncertainties for next-generation ground-based observatories targeting the large-angle B-mode polarization of the Cosmic Microwave Background, with a focus on the Simons Observatory (SO). We explore uncertainties on gain calibration, bandpass center frequencies, and polarization angles, including the frequency variation of the latter across the bandpass. We find that gain calibration and bandpass center frequencies must be known to percent levels or less to avoid biases on the tensor-to-scalar ratio r on the order of Δ r∼10-3, in line with previous findings. Polarization angles must be calibrated to the level of a few tenths of a degree, while their frequency variation between the edges of the band must be known to O(10) degrees. Given the tightness of these calibration requirements, we explore the level to which residual uncertainties on these systematics would affect the final constraints on r if included in the data model and marginalized over. We find that the additional parameter freedom does not degrade the final constraints on r significantly, broadening the error bar by O(10%) at most. We validate these results by reanalyzing the latest publicly available data from the collaboration within an extended parameter space covering both cosmological, foreground and systematic parameters. Finally, our results are discussed in light of the instrument design and calibration studies carried out within SO.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Evaluation of advanced automatic PET segmentation methods using nonspherical thin-wall inserts

Purpose: The use of positron emission tomography (PET) within radiotherapy treatment planning requires the availability of reliable and accurate segmentation tools. PET automatic segmentation (PET-AS) methods have been recommended for the delineation of tumors, but there is still a lack of thorough validation and cross-comparison of such methods using clinically relevant data. In particular, studies validating PET segmentation tools mainly use phantoms with thick plastic walls inserts of simple spherical geometry and have not specifically investigated the effect of the target object geometry on the delineation accuracy. Our work therefore aimed at generating clinically realistic data using nonspherical thin-wall plastic inserts, for the evaluation and comparison of a set of eight promising PET-AS approaches. Methods: Sixteen nonspherical inserts were manufactured with a plastic wall of 0.18 mm and scanned within a custom plastic phantom. These included ellipsoids and toroids derived with different volumes, as well as tubes, pear- and drop-shaped inserts with different aspect ratios. A set of six spheres of volumes ranging from 0.5 to 102 ml was used for a baseline study. A selection of eight PET-AS methods, written in house, was applied to the images obtained. The methods represented promising segmentation approaches such as adaptive iterative thresholding, region-growing, clustering and gradient-based schemes. The delineation accuracy was measured in terms of overlap with the computed tomography reference contour, using the dice similarity coefficient (DSC), and error in dimensions. Results: The delineation accuracy was lower for nonspherical inserts than for spheres of the same volume in 88% cases. Slice-by-slice gradient-based methods, showed particularly lower DSC for tori (DSC < 0.5), caused by a failure to recover the object geometry. The region-growing method reached high levels of accuracy for most inserts (DSC > 0.76 except for tori) but showed the largest errors in the recovery of pears and drops dimensions (higher than 10% and 30% of the true length, respectively). Large errors were visible for one of the gradient-based contouring methods when delineating drop-shaped inserts. Low DSC due to systematic underestimation of the volumes was observed for our fuzzy clustering method when using nonspherical inserts. The adaptive iterative thresholding method produced the highest DSC score on our nonspherical dataset (DSC > 0.83, except for tori) and showed robustness to the insert geometry. Conclusions: This study investigated the accuracy of eight PET-AS methods for the delineation of objects with a range of nonspherical geometries. The authors’ results confirmed the robustness of some segmentation approaches, but also showed the weaknesses of some of the other methods implemented, which were not observed with spherical inserts. This work therefore highlights the importance of using a variety of thin-wall inserts with complex geometries for the validation of PET-AS methods and provided useful information for further development of the methods tested.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Design and performance of a variable gap system for thermal conductivity measurements of high temperature, corrosive, and reactive fluids

We report high-temperature fluids such as molten salts, liquid metals, and gasses are being proposed for many advanced energy systems including thermal energy storage devices, concentrating solar plants, and advanced nuclear reactor designs. However, the chemical behavior and thermophysical properties of many of these fluids have not been well characterized, which hinders the design, modeling, safety analysis, and deployment of these systems. Thermal conductivity is a property that is especially limited by existing measurement capabilities, which are subject to errors caused by convection, material interaction, radiative heat transfer, and instrument degradation. Therefore, there is a lack of standard, systematic measurement techniques for high-temperature, reactive, and corrosive fluids. In this work, the development of a variable gap thermal conductivity measurement system is detailed. The system is designed to measure the thermal conductivity of highly corrosive and reactive fluids, and survive operation between 100 °C and 800 °C. The effects of convection are minimized by limiting the thickness of the specimen to thin sizes (<0.3 mm). Corrections for radiative heat transfer were included in the working equations to consider specimens with varying optical properties. The design, construction, instrumentation, operating principles, and data analysis techniques are discussed in detail. The system was tested up to 500 °C using helium gas and molten KNO 3 -NaNO 3 to verify the measurement technique and determine the sources of error. At 300 and 400 °C KNO 3 -NaNO 3 , results showed maximum relative error of 6% when compared to results in the literature. The helium results were within 13% of those in the literature at 300 and 400 °C. Higher errors were observed at 500 °C for both fluids, and the sources of these errors are discussed.

42 ENGINEERING↗
DOE OSTI2016

Proton radius from electron scattering data

Background: The proton charge radius extracted from recent muonic hydrogen Lamb shift measurements is significantly smaller than that extracted from atomic hydrogen and electron scattering measurements. The discrepancy has become known as the proton radius puzzle. Purpose: In an attempt to understand the discrepancy, we review high-precision electron scattering results from Mainz, Jefferson Lab, Saskatoon and Stanford. Methods: We make use of stepwise regression techniques using the F-test as well as the Akaike information criterion to systematically determine the predictive variables to use for a given set and range of electron scattering data as well as to provide multivariate error estimates. Results: Starting with the precision, low four-momentum transfer (Q 2 ) data from Mainz (1980) and Saskatoon (1974), we find that a stepwise regression of the Maclaurin series using the F-test as well as the Akaike information criterion justify using a linear extrapolation which yields a value for the proton radius that is consistent with the result obtained from muonic hydrogen measurements. Applying the same Maclaurin series and statistical criteria to the 2014 Rosenbluth results on GE from Mainz, we again find that the stepwise regression tends to favor a radius consistent with the muonic hydrogen radius but produces results that are extremely sensitive to the range of data included in the fit. Making use of the high-Q 2 data on G E to select functions which extrapolate to high Q 2 , we find that a Pad´e (N = M = 1) statistical model works remarkably well, as does a dipole function with a 0.84 fm radius, G E (Q 2 ) = (1 + Q 2 /0.66 GeV 2 ) -2 . Conclusions: Rigorous applications of stepwise regression techniques and multivariate error estimates result in the extraction of a proton charge radius that is consistent with the muonic hydrogen result of 0.84 fm; either from linear extrapolation of the extreme low-Q 2 data or by use of the Pad´e approximant for extrapolation using a larger range of data. Thus, based on a purely statistical analysis of electron scattering data, we conclude that the electron scattering result and the muonic hydrogen result are consistent. Lastly, it is the atomic hydrogen results that are the outliers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Elemental ratio measurements of organic compounds using aerosol mass spectrometry: characterization, improved calibration, and implications

Elemental compositions of organic aerosol (OA) particles provide useful constraints on OA sources, chemical evolution, and effects. The Aerodyne high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) is widely used to measure OA elemental composition. This study evaluates AMS measurements of atomic oxygen-to-carbon (O : C), hydrogen-to-carbon (H : C), and organic mass-to-organic carbon (OM : OC) ratios, and of carbon oxidation state ( OS C ) for a vastly expanded laboratory data set of multifunctional oxidized OA standards. For the expanded standard data set, the method introduced by Aiken et al. (2008), which uses experimentally measured ion intensities at all ions to determine elemental ratios (referred to here as "Aiken-Explicit"), reproduces known O : C and H : C ratio values within 20% (average absolute value of relative errors) and 12%, respectively. The more commonly used method, which uses empirically estimated H 2 O + and CO + ion intensities to avoid gas phase air interferences at these ions (referred to here as "Aiken-Ambient"), reproduces O : C and H : C of multifunctional oxidized species within 28 and 14% of known values. The values from the latter method are systematically biased low, however, with larger biases observed for alcohols and simple diacids. A detailed examination of the H 2 O + , CO + , and CO 2 + fragments in the high-resolution mass spectra of the standard compounds indicates that the Aiken-Ambient method underestimates the CO + and especially H 2 O + produced from many oxidized species. Combined AMS–vacuum ultraviolet (VUV) ionization measurements indicate that these ions are produced by dehydration and decarboxylation on the AMS vaporizer (usually operated at 600 °C). Thermal decomposition is observed to be efficient at vaporizer temperatures down to 200 °C. These results are used together to develop an "Improved-Ambient" elemental analysis method for AMS spectra measured in air. The Improved-Ambient method uses specific ion fragments as markers to correct for molecular functionality-dependent systematic biases and reproduces known O : C (H : C) ratios of individual oxidized standards within 28% (13%) of the known molecular values. The error in Improved-Ambient O : C (H : C) values is smaller for theoretical standard mixtures of the oxidized organic standards, which are more representative of the complex mix of species present in ambient OA. For ambient OA, the Improved-Ambient method produces O : C (H : C) values that are 27% (11%) larger than previously published Aiken-Ambient values; a corresponding increase of 9% is observed for OM : OC values. These results imply that ambient OA has a higher relative oxygen content than previously estimated. The OS C values calculated for ambient OA by the two methods agree well, however (average relative difference of 0.06 OS C units). This indicates that OS C is a more robust metric of oxidation than O : C, likely since OS C is not affected by hydration or dehydration, either in the atmosphere or during analysis.

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

OBSERVATION AND SPECTRAL MEASUREMENTS OF THE CRAB NEBULA WITH MILAGRO

The Crab Nebula was detected with the Milagro experiment at a statistical significance of 17 standard deviations over the lifetime of the experiment. The experiment was sensitive to approximately 100 GeV-100 TeV gamma-ray air showers by observing the particle footprint reaching the ground. The fraction of detectors recording signals from photons at the ground is a suitable proxy for the energy of the primary particle and has been used to measure the photon energy spectrum of the Crab Nebula between {approx}1 and {approx}100 TeV. The TeV emission is believed to be caused by inverse-Compton upscattering of ambient photons by an energetic electron population. The location of a TeV steepening or cutoff in the energy spectrum reveals important details about the underlying electron population. We describe the experiment and the technique for distinguishing gamma-ray events from the much more-abundant hadronic events. We describe the calculation of the significance of the excess from the Crab and how the energy spectrum is fitted. The differential photon energy spectrum, including the statistical errors from the fit, obtained using a simple power-law hypothesis for data between 2005 September and 2008 March is (6.5 {+-} 0.4) Multiplication-Sign 10{sup -14}(E/10 TeV){sup -3.1{+-}0.1}(cm{sup 2} s TeV ){sup -1} between {approx}1 TeV and {approx}100 TeV. Allowing for a possible exponential cutoff, the photon energy spectrum is fitted as (2.5{sup +0.7}{sub -0.4}) Multiplication-Sign 10{sup -12}(E/3 TeV){sup -2.5{+-}0.4}exp (- E/32{sup +39}{sub -18} TeV) (cm{sup 2} s TeV){sup -1}. The results are subject to an {approx}30% systematic uncertainty in the overall flux and an {approx}0.1 systematic uncertainty in the power-law indices quoted. Uncertainty in the overall energy scale has been absorbed into these errors. Fixing the spectral index to values that have been measured below 1 TeV by IACT experiments (2.4-2.6), the fit to the Milagro data suggests that Crab exhibits a spectral steepening or cutoff between about 20-40 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Constraining the baryon-dark matter relative velocity with the large-scale 3-point correlation function of the SDSS BOSS DR12 CMASS galaxies

We search for a galaxy clustering bias due to a modulation of galaxy number with the baryon-dark matter relative velocity resulting from recombination-era physics. We find no detected signal and place the constraint bv <0.01 on the relative velocity bias for the CMASS galaxies. This bias is an important potential systematic of Baryon Acoustic Oscillation (BAO) method measurements of the cosmic distance scale using the 2-point clustering. Our limit on the relative velocity bias indicates a systematic shift of no more than 0.3% rms in the distance scale inferred from the BAO feature in the BOSS 2-point clustering, well below the 1% statistical error of this measurement. In conclusion, this constraint is the most stringent currently available and has important implications for the ability of upcoming large-scale structure surveys such as DESI to self-protect against the relative velocity as a possible systematic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Systematically improvable multiscale solver for correlated electron systems

The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe a framework for a general multiscale method based on embedding a self-energy of a strongly correlated subsystem into a self-energy generated by a method able to treat large weakly correlated systems approximately. As an example, we present the embedding of an exact diagonalization self-energy into a self-energy generated from self-consistent second-order perturbation theory. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the two-dimensional Hubbard model, as a benchmark, we illustrate that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the method in multiple regimes of interaction strengths and dopings of the model. Here, the general embedding framework we present avoids difficulties such as double counting corrections, frequency-dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantification of modeling uncertainty in the Rayleigh damping model

Understanding and accurately characterizing energy dissipation mechanisms in civil structures during earthquakes is an important element of seismic assessment and design. The most commonly used model is attributed to Rayleigh. This paper proposes a systematic approach to quantify the uncertainty associated with Rayleigh's damping model. Bayesian calibration with embedded model error is employed to treat the coefficients of the Rayleigh model as random variables using modal damping ratios. Through a numerical example, we illustrate how this approach works and how the calibrated model can address modeling uncertainty associated with the Rayleigh damping model.

42 ENGINEERING↗
DOE OSTI2023

Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT

Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.

58 GEOSCIENCES↗
DOE OSTI2015

CORONAL MAGNETOGRAPHY OF A SIMULATED SOLAR ACTIVE REGION FROM MICROWAVE IMAGING SPECTROPOLARIMETRY

We have simulated the Expanded Owens Valley Solar Array (EOVSA) radio images generated at multiple frequencies from a model solar active region, embedded in a realistic solar disk model, and explored the resulting data cube for different spectral analysis schemes to evaluate the potential for realizing one of EOVSA’s most important scientific goals—coronal magnetography. In this paper, we focus on modeling the gyroresonance and free–free emission from an on-disk solar active region model with realistic complexities in electron density, temperature and magnetic field distribution. We compare the magnetic field parameters extrapolated from the image data cube along each line of sight after folding through the EOVSA instrumental profile with the original (unfolded) parameters used in the model. We find that even the most easily automated, image-based analysis approach (Level-0) provides reasonable quantitative results, although they are affected by systematic effects due to finite sampling in the Fourier (UV) plane. Finally, we note the potential for errors due to misidentified harmonics of the gyrofrequency, and discuss the prospects for applying a more sophisticated spectrally based analysis scheme (Level-1) to resolve the issue in cases where improved UV coverage and spatial resolution are available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Using angular two-point correlations to self-calibrate the photometric redshift distributions of DECaLS DR9

Calibrating the redshift distributions of photometric galaxy samples is essential in weak lensing studies. The self-calibration method combines angular auto- and cross-correlations between galaxies in multiple photometric redshift (photo-z) bins to reconstruct the scattering rates matrix between redshift bins. In this paper, we test a recently proposed self-calibration algorithm using the DECaLS Data Release 9 and investigate to what extent the scattering rates are determined. We first mitigate the spurious angular correlations due to imaging systematics by a machine learning based method. We then improve the algorithm for χ 2 minimization and error estimation. Finally, we solve for the scattering matrices, carry out a series of consistency tests, and find reasonable agreements: (1) finer photo-z bins return a high-resolution scattering matrix, and it is broadly consistent with the low-resolution matrix from wider bins; (2) the scattering matrix from the Northern Galactic Cap is almost identical to that from Southern Galactic Cap; (3) the scattering matrices are in reasonable agreement with those constructed from the power spectrum and the weighted spectroscopic subsample. We also evaluate the impact of cosmic magnification. Although it changes little the diagonal elements of the scattering matrix, it affects the off-diagonals significantly. The scattering matrix also shows some dependence on scale cut of input correlations, which may be related to a known numerical degeneracy between certain scattering pairs. In conclusion, this work demonstrates the feasibility of the self-calibration method in real data and provides a practical alternative to calibrate the redshift distributions of photometric samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Direct Measurement of the Nickel Spectrum in Cosmic Rays in the Energy Range from 8.8 GeV/$\mathcal{n}$ to 240 GeV/$\mathcal{n}$ with CALET on the International Space Station

The relative abundance of cosmic ray nickel nuclei with respect to iron is by far larger than for all other transiron elements; therefore it provides a favorable opportunity for a low background measurement of its spectrum. Since nickel, as well as iron, is one of the most stable nuclei, the nickel energy spectrum and its relative abundance with respect to iron provide important information to estimate the abundances at the cosmic ray source and to model the Galactic propagation of heavy nuclei. However, only a few direct measurements of cosmic-ray nickel at energy larger than ~ 3 GeV/n are available at present in the literature, and they are affected by strong limitations in both energy reach and statistics. In this Letter, we present a measurement of the differential energy spectrum of nickel in the energy range from 8.8 to 240 GeV/n , carried out with unprecedented precision by the Calorimetric Electron Telescope (CALET) in operation on the International Space Station since 2015. The CALET instrument can identify individual nuclear species via a measurement of their electric charge with a dynamic range extending far beyond iron (up to atomic number Z = 40 ). The particle’s energy is measured by a homogeneous calorimeter (1.2 proton interaction lengths, 27 radiation lengths) preceded by a thin imaging section (3 radiation lengths) providing tracking and energy sampling. This Letter follows our previous measurement of the iron spectrum [O. Adriani et al. (CALET Collaboration), Phys. Rev. Lett. 126, 241101 (2021).], and it extends our investigation on the energy dependence of the spectral index of heavy elements. It reports the analysis of nickel data collected from November 2015 to May 2021 and a detailed assessment of the systematic uncertainties. In the region from 20 to 240 GeV/n our present data are compatible within the errors with a single power law with spectral index -2.51 ± 0.07 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Single Pass Collider Memo: Gradient Perturbations of the SLC arc

As the beam passes through the arcs, the gradient it encounters at each magnet differs from the design value. This deviation may be in part random and in part systematic. In this note we make estimates of the effects to be expected from both kinds of errors.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

A 2.4% DETERMINATION OF THE LOCAL VALUE OF THE HUBBLE CONSTANT

We use the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) to reduce the uncertainty in the local value of the Hubble constant from 3.3% to 2.4%. The bulk of this improvement comes from new near-infrared (NIR) observations of Cepheid variables in 11 host galaxies of recent type Ia supernovae (SNe Ia), more than doubling the sample of reliable SNe Ia having a Cepheid-calibrated distance to a total of 19; these in turn leverage the magnitude-redshift relation based on ∼300 SNe Ia at z < 0.15. All 19 hosts as well as the megamaser system NGC 4258 have been observed with WFC3 in the optical and NIR, thus nullifying cross-instrument zeropoint errors in the relative distance estimates from Cepheids. Other noteworthy improvements include a 33% reduction in the systematic uncertainty in the maser distance to NGC 4258, a larger sample of Cepheids in the Large Magellanic Cloud (LMC), a more robust distance to the LMC based on late-type detached eclipsing binaries (DEBs), HST observations of Cepheids in M31, and new HST -based trigonometric parallaxes for Milky Way (MW) Cepheids. We consider four geometric distance calibrations of Cepheids: (i) megamasers in NGC 4258, (ii) 8 DEBs in the LMC, (iii) 15 MW Cepheids with parallaxes measured with HST /FGS, HST /WFC3 spatial scanning and/or Hipparcos , and (iv) 2 DEBs in M31. The Hubble constant from each is 72.25 ± 2.51, 72.04 ± 2.67, 76.18 ± 2.37, and 74.50 ± 3.27 km s{sup 1} Mpc{sup 1}, respectively. Our best estimate of H {sub 0} = 73.24 ± 1.74 km s{sup 1} Mpc{sup 1} combines the anchors NGC 4258, MW, and LMC, yielding a 2.4% determination (all quoted uncertainties include fully propagated statistical and systematic components). This value is 3.4 σ higher than 66.93 ± 0.62 km s{sup 1} Mpc{sup 1} predicted by ΛCDM with 3 neutrino flavors having a mass of 0.06 eV and the new Planck data, but the discrepancy reduces to 2.1 σ relative to the prediction of 69.3 ± 0.7 km s{sup 1} Mpc{sup 1} based on the comparably precise combination of WMAP +ACT+SPT+BAO observations, suggesting that systematic uncertainties in CMB radiation measurements may play a role in the tension. If we take the conflict between Planck high-redshift measurements and our local determination of H {sub 0} at face value, one plausible explanation could involve an additional source of dark radiation in the early universe in the range of Δ N {sub eff} ≈ 0.4–1. We anticipate further significant improvements in H {sub 0} from upcoming parallax measurements of long-period MW Cepheids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Joint analysis of galaxy-galaxy lensing and galaxy clustering: Methodology and forecasts for Dark Energy Survey

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large scale structure. Our analysis will be carried out on data from the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the tangential shears of background galaxies induced by these foreground lenses. We develop a practical approach to modeling the assumptions and systematic effects affecting small scale lensing, which provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement errors, we study how external priors on different subsets of these parameters affect our growth constraints. Degeneracies within the HOD model, as well as between the HOD and the growth function, are identified as the dominant source of complication, with other systematic effects sub-dominant. The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of the galaxy sample that we employ. Finally, we conclude that DES data will provide powerful constraints on the evolution of structure growth in the universe, conservatively/optimistically constraining the growth function to 7.9%/4.8% with its first-year data that covered over 1000 square degrees, and to 3.9%/2.3% with its full five-year data that will survey 5000 square degrees, including both statistical and systematic uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Nucleon momentum fraction, helicity and transversity from 2+1-flavor lattice QCD

A detailed analysis of the systematic uncertainties in the calculation of the isovector momentum fraction, ( x ) u - d , helicity moment, ( x ) Δ u - Δ d , and the transversity moment, ( x ) δu - δd , of the nucleon is presented using high-statistics data on seven ensembles of gauge configurations generated by the JLab/W&M/LANL/MIT collaborations using 2 + 1-flavors of dynamical Wilson-clover quarks. The much higher statistics have facilitated better control over all systematics compared to previous lattice calculations. The least understood systematic — excited-state contamination — is quantified by studying the variation of the results as a function of different estimates of the mass gap of the first excited state, obtained from two- and three-point correlation functions, and as a function of the pion mass M π . The final results are obtained using a simultaneous fit in the lattice spacing a , pion mass M π and the finite volume parameter M π L keeping leading order corrections. The data show no significant dependence on the lattice spacing and some evidence for finite-volume corrections. Our final results, in the $ \overline{\mathrm{MS}} $ scheme at 2 GeV, are ( x ) u - d = 0 . 155(17)(20), ( x ) Δ u - Δ d = 0 . 183(14)(20) and ( x ) δu - δd = 0 . 220(18)(20), where the first error is the overall analysis uncertainty assuming excited-state contributions have been removed, and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Volume-preserving algorithm for secular relativistic dynamics of charged particles

Secular dynamics of relativistic charged particles has theoretical significance and a wide range of applications. However, conventional algorithms are not applicable to this problem due to the coherent accumulation of numerical errors. To overcome this difficulty, we develop a volume-preserving algorithm (VPA) with long-term accuracy and conservativeness via a systematic splitting method. Applied to the simulation of runaway electrons with a time-span over 10 magnitudes, the VPA generates accurate results and enables the discovery of new physics for secular runaway dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗