Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “galaxy structure”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Intrinsic alignment as an RSD contaminant in the DESI survey

ABSTRACT We measure the tidal alignment of the major axes of luminous red galaxies (LRGs) from the Legacy Imaging Survey and use it to infer the artificial redshift-space distortion signature that will arise from an orientation-dependent, surface-brightness selection in the Dark Energy Spectroscopic Instrument (DESI) survey. Using photometric redshifts to downweight the shape–density correlations due to weak lensing, we measure the intrinsic tidal alignment of LRGs. Separately, we estimate the net polarization of LRG orientations from DESI’s fibre-magnitude target selection to be of order 10−2 along the line of sight. Using these measurements and a linear tidal model, we forecast a 0.5 per cent fractional decrease on the quadrupole of the two-point correlation function for projected separations of 40–80 h−1 Mpc. We also use a halo catalogue from the Abacussummit cosmological simulation suite to reproduce this false quadrupole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The imprint of cosmic voids from the DESI Legacy Survey DR9 Luminous Red Galaxies in the Planck 2018 lensing map through spectroscopically calibrated mocks

The cross-correlation of cosmic voids with the lensing convergence (κ) map of the Cosmic Microwave Background (CMB) fluctuations provides a powerful tool to refine our understanding of the current cosmological model. However, several studies have reported a moderate tension (up to ∼2σ) between the lensing imprint of cosmic voids on the observed CMB and the ΛCDM signal predicted by simulations. To address this “lensing-is-low” tension and to obtain new, precise measurements of the signal, we exploit the large DESI Legacy Survey Luminous Red Galaxy (LRG) data set, covering approximately 19 500 deg 2 of the sky and including about 10 million LRGs at z < 1.05. Our ΛCDM template was created using the Buzzard mocks, which we specifically calibrated to match the clustering properties of the observed galaxy sample by exploiting more than one million DESI spectra. We identified our catalogs of 3D voids in the range 0.35 < z < 0.95 and cross-correlated them through a stacking methodology, dividing the sample into bins according to the redshift and λv values of the voids. For the full void sample, we report a 14σ detection of the lensing signal, with A κ = 1.016 ± 0.054, which increases to 17σ when considering the void-in-void (A κ = 0.944 ± 0.064) and the void-in-cloud (A κ = 0.975 ± 0.060) populations individually, the highest detection significance for studies of this kind. We observe a full agreement between observations and ΛCDM mocks across all redshift bins, sky regions, and void populations considered. In addition to these findings, our analysis highlights the importance of accurately matching sparseness and redshift error distributions between mocks and observations, as well as the role of λ v in enhancing the signal-to-noise ratio through void population discrimination.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Cosmology from cosmic shear power spectra

Here, we measure weak lensing cosmic shear power spectra from the 3-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers 416 deg 2 of the northern sky, with a mean i-band seeing of 0.59 arcsec and an effective galaxy number density of 15 arcmin —2 within our adopted redshift range. With an i-band magnitude limit of 24.5 mag, and four tomographic redshift bins spanning 0.3 ≤ z ph ≤ 1.5 based on photometric redshifts, we obtain a high-significance measurement of the cosmic shear power spectra, with a signal-to-noise ratio of approximately 26.4 in the multipole range 300 < ℓ < 1800. The accuracy of our power spectrum measurement is tested against realistic mock shear catalogs, and we use these catalogs to get a reliable measurement of the covariance of the power spectrum measurements. We use a robust blinding procedure to avoid confirmation bias, and model various uncertainties and sources of bias in our analysis, including point spread function systematics, redshift distribution uncertainties, the intrinsic alignment of galaxies and the modeling of the matter power spectrum. For a flat Λ CDM model, we find S 8 ≡ σ 8 (Ω m /0.3) 0.5 = ${0.776}_{—0.033}^{+0.032}$, which is in excellent agreement with the constraints from the other HSC Year 3 cosmology analyses, as well as those from a number of other cosmic shear experiments. This result implies a ~2σ-level tension with the Planck 2018 cosmology. We study the effect that various systematic errors and modeling choices could have on this value, and find that they can shift the best-fit value of S 8 by no more than ~ 0.5σ, indicating that our result is robust to such systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Probing quasar viewing angle with the variability structure function

Given the anisotropic emission from quasar accretion discs, their viewing angle affects estimates of the quasar luminosity, black hole mass, and Eddington ratio. Discs appear overluminous when viewed pole-on and underluminous when viewed at high inclination. In radio-quiet quasars, the viewing angle is usually unknown, although spectroscopic indicators have been proposed. Here, we use a recently discovered universality in the variability structure function (SF) of quasar light curves (LCs), where all quasars show the same SF when clocks run in units of orbital time-scale. As an offset from the mean relation can be caused by incorrect orbital time-scales and thus incorrect luminosities, we correlate these offsets with suggested inclination indicators. We derive SFs from NASA/Asteroid Terrestrial-impact Last Alert System (ATLAS) LCs spanning ~6 yr of observation, using a sample of 183 luminous quasars with measured H β lines as well as 753 quasars with C $\tiny{IV}$ and Mg $\tiny{II}$ lines. Starting from the proposed orientation indicators, we expect quasars with narrower H β lines and with more blueshifted C $\tiny{IV}$ lines to be viewed more pole-on and thus appear overluminous. In contrast, our SF analysis finds that presumed pole-on discs appear underluminous, consistently for both line indicators. We discuss possible explanations for the behaviour of quasars with highly blueshifted C $\tiny{IV}$ lines irrespective of inclination angle, including dusty outflows that might render the accretion disc underluminous and flatter disc temperature profiles with longer orbital time-scales than in thin-disc models but reach no satisfying conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Blinding multiprobe cosmological experiments

ABSTRACT The goal of blinding is to hide an experiment’s critical results – here the inferred cosmological parameters – until all decisions affecting its analysis have been finalized. This is especially important in the current era of precision cosmology, when the results of any new experiment are closely scrutinized for consistency or tension with previous results. In analyses that combine multiple observational probes, like the combination of galaxy clustering and weak lensing in the Dark Energy Survey (DES), it is challenging to blind the results while retaining the ability to check for (in)consistency between different parts of the data. We propose a simple new blinding transformation, which works by modifying the summary statistics that are input to parameter estimation, such as two-point correlation functions. The transformation shifts the measured statistics to new values that are consistent with (blindly) shifted cosmological parameters while preserving internal (in)consistency. We apply the blinding transformation to simulated data for the projected DES Year 3 galaxy clustering and weak lensing analysis, demonstrating that practical blinding is achieved without significant perturbation of internal-consistency checks, as measured here by degradation of the χ2 between the data and best-fitting model. Our blinding method’s performance is expected to improve as experiments evolve to higher precision and accuracy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

CODEX clusters

Context: Large area catalogs of galaxy clusters constructed from ROSAT All-Sky Survey provide the basis for our knowledge of the population of clusters thanks to long-term multiwavelength efforts to follow up observations of these clusters. Aims: The advent of large area photometric surveys superseding previous, in-depth all-sky data allows us to revisit the construction of X-ray cluster catalogs, extending the study to lower cluster masses and higher redshifts and providing modeling of the selection function. Methods: We performed a wavelet detection of X-ray sources and made extensive simulations of the detection of clusters in the RASS data. We assigned an optical richness to each of the 24 788 detected X-ray sources in the 10 382 square degrees of the Baryon Oscillation Spectroscopic Survey area using red sequence cluster finder redMaPPer version 5.2 run on Sloan Digital Sky Survey photometry. We named this survey COnstrain Dark Energy with X-ray (CODEX) clusters. Results: We show that there is no obvious separation of sources on galaxy clusters and active galactic nuclei (AGN) based on the distribution of systems on their richness. This is a combination of an increasing number of galaxy groups and their selection via the identification of X-ray sources either by chance or by groups hosting an AGN. To clean the sample, we use a cut on the optical richness at the level corresponding to the 10% completeness of the survey and include it in the modeling of the cluster selection function. We present the X-ray catalog extending to a redshift of 0.6. Conclusions: The CODEX suvey is the first large area X-ray selected catalog of northern clusters reaching fluxes of 10 –13 ergs s –1 cm –2 . We provide modeling of the sample selection and discuss the redshift evolution of the high end of the X-ray luminosity function (XLF). Our results on z < 0.3 XLF agree with previous studies, while we provide new constraints on the 0.3 < z < 0.6 XLF. We find a lack of strong redshift evolution of the XLF, provide exact modeling of the effect of low number statistics and AGN contamination, and present the resulting constraints on the flat ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Probing parity violation with the four-point correlation function of BOSS galaxies

Parity-violating physics in the early universe can leave detectable traces in late-time observables. While vector- and tensor-type parity violation can be observed in the B -modes of the cosmic microwave background, scalar-type signatures are visible only in the four-point correlation function (4PCF) and beyond. This work presents a blind test for parity violation in the 4PCF of the baryon oscillation spectroscopic survey (BOSS) CMASS sample, considering galaxy separations in the range [20; 160]h -1 Mpc The parity-odd 4PCF contains no contributions from standard Λ CDM physics and can be efficiently measured using recently developed estimators. Data are analyzed using both a nonparametric rank test (comparing the BOSS 4PCFs to those of realistic simulations) and a compressed χ 2 analysis, with the former avoiding the assumption of a Gaussian likelihood. These find similar results, with the rank test giving a detection probability of 99.6% (2.9 σ). This provides significant evidence for parity violation, from either cosmological sources or systematics. We perform a number of systematic tests: although these do not reveal any observational artifacts, we cannot exclude the possibility that our detection is caused by the simulations not faithfully representing the statistical properties of the BOSS data. Our measurements can be used to constrain physical models of parity violation. As an example, we consider a coupling between the inflaton and a U (1) gauge field and place bounds on the latter’s energy density, which are several orders of magnitude stronger than those previously reported. Upcoming probes such as DESI and Euclid will reveal whether our detection of parity violation is due to new physics, and strengthen the bounds on a variety of models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Minkowski Functionals of SDSS-III BOSS: Hints of Possible Anisotropy in the Density Field?

We present measurements of the Minkowski functionals extracted from the SDSS-III BOSS catalogs. After defining the Minkowski functionals, we describe how an unbiased reconstruction of these statistics can be obtained from a field with masked regions and survey boundaries, validating our methodology with Gaussian random fields and mock galaxy snapshot data. From the BOSS galaxy data, we generate a set of four density fields in three dimensions corresponding to the northern and southern skies of LOWZ and CMASS catalogs, smoothing over large scales (Gaussian smoothing scale of 35 Mpc) such that the field is perturbatively non-Gaussian. We extract the Minkowski functionals from each data set separately, and measure their shapes and amplitudes by fitting a Hermite polynomial expansion. For the shape parameter of the Minkowski functional curves a0, that is related to the bispectrum of the field, we find that the LOWZ-South data presents a systematically lower value of a 0 = –0.080 ± 0.040 than its northern sky counterpart a 0 = 0.032 ± 0.024. Although the significance of this discrepancy is not very high, it potentially indicates some systematics in the data or that the matter density field exhibits anisotropy at low redshift. By assuming a standard isotropic flat ΛCDM cosmology, the amplitudes of the Minkowski functionals from the combination of northern and southern sky data give the constraints Ω c h 2 n s = 0.110 ± 0.006 and 0.111 ± 0.008 for CMASS and LOWZ, respectively, which is in agreement with the Planck ΛCDM best-fit Ω c h 2 n s = 0.116 ± 0.001.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Probing circular polarization and magnetic field structure in active galactic nuclei

The composition and magnetic field morphology of relativistic jets can be studied using circular polarization (CP) measurement. Recent three-dimensional relativistic magnetohydrodynamic (3D RMHD) simulations coupled with radiative transfer (RT) calculations make strong predictions about the level (and morphology) of the jet’s CP emission. These simulations show that the sign of CP and the electric vector position angle (EVPA) are both sensitive to the jet’s magnetic field morphology within the radio core. We probe this theory by exploring whether the jet’s radio core EVPA orientation is consistent with the observed sign of the core CP in deep full-track polarimetric observations. Based on a selection of sources from earlier MOJAVE observations, we aim to probe the nature of linear polarization (LP) and CP in the innermost regions of jets from a small sample of nine blazars. This sample includes sources that have exhibited: (i) positive CP; (ii) negative CP; or (iii) positive and negative CP simultaneously in the radio core region. By coupling deep polarimetric observations of a carefully selected sample of blazars with state-of-the-art RMHD and RT calculations, we hope to gain a deeper understanding of the physics of blazar jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Constraining gravity with a new precision 𝐸 𝐺 estimator using Planck + SDSS BOSS data

The 𝐸 𝐺 statistic is a discriminating probe of gravity developed to test the prediction of general relativity (GR) for the relation between gravitational potential and clustering on the largest scales in the observable Universe. We present a novel high-precision estimator for the 𝐸 𝐺 statistic using CMB lensing and galaxy clustering correlations that carefully matches the effective redshifts across the different measurement components to minimize corrections. A suite of detailed tests is performed to characterize the estimator’s accuracy, its sensitivity to assumptions and analysis choices, and the non-Gaussianity of the estimator’s uncertainty is characterized. After finalization of the estimator, it is applied to Planck CMB lensing and SDSS CMASS and LOWZ galaxy data. We report the first harmonic space measurement of 𝐸 𝐺 using the LOWZ sample and CMB lensing and also updated constraints using the final CMASS sample and the latest Planck CMB lensing map. We find $\hat{𝐸}$$^{Planck+CMASS}_{𝐺}$ = 0.3⁢6$^{+0.06}_{−0.05}$⁢(68.27%) and $\hat{𝐸}$$^{Planck+LOWZ}_{𝐺}$ = 0.4⁢0$^{+0.11}_{−0.09}$⁢(68.27%), with additional subdominant systematic error budget estimates of 2% and 3%, respectively. Using Ω m,0 constraints from Planck and SDSS BAO observations, Λ⁢CDM-GR predicts 𝐸$^{GR}_ {𝐺}$⁡(𝑧 =0.555) = 0.401 ± 0.005 and 𝐸$^{GR}_{𝐺}$⁡(𝑧 =0.316) = 0.452 ± 0.005 at the effective redshifts of the CMASS and LOWZ based measurements. We report the measurement to be in good statistical agreement with the Λ⁢CDM-GR prediction and report that the measurement is also consistent with the more general GR prediction of scale independence for 𝐸 𝐺 . Furthermore, this work provides a carefully constructed and calibrated statistic with which 𝐸 𝐺 measurements can be confidently and accurately obtained with upcoming survey data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellar kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within 1.2σ, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a ~4σ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within ~0.6σ<, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Information content of higher order galaxy correlation functions

The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance–redshift relationship. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence, all the cosmological information can be extracted from their two-point correlation function. Subsequent non-linear evolution under gravity, as well as halo and then galaxy formation, generates higher order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field’s power spectrum is equal to that of all the higher order functions of the final, non-linear field. This claim implies that reconstruction of the initial density field from the non-linear field renders analysis of higher order correlation functions of the latter superfluous. We show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the two and three-point correlation functions can, in some cases, exceed those offered by the initial power spectrum. We provide a mathematical justification for this claim and demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot-noise, taking modes up to k max = 0.2 h Mpc −1 , using the bispectrum alone offers a factor of 2 reduction in the variance on the cosmic distance scale relative to that available from the linear power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Merger Response of Halo Anisotropy Properties

Anisotropy properties—halo spin, shape, position offset, velocity offset, and orientation—are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of halos. These properties reflect the dynamical state of halos, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy properties in response to merger activity using the IllustrisTNG simulations. We find that the response trajectories of the anisotropy properties significantly deviate from secular evolution. These trajectories have the same qualitative features and timescales across a wide range of merger and host properties. We propose explanations for the behavior of these properties and connect their evolution to the relevant stages of merger dynamics. We measure the relevant dynamical timescales. We also explore the dependence of the strength of the response on time of merger, merger ratio, and mass of the main halo. These results provide insight into the physics of halo mergers and their effects on the statistical behavior of halo properties. This study paves the way toward a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The H I absorption zoo: JVLA extension to z ~ 0.4

Here, we present an H I 21 cm absorption study of a sample of 26 radio-loud active galactic nuclei (AGN) at 0.25 < z < 0.4 carried out with the Karl G. Jansky Very Large Array. Our aim was to study the rate of incidence of H I in different classes of radio AGN, the morphology and kinematics of the detected H I, and the nature of the interaction between the H I and the radio source at these redshifts. Our sample consists of 14 sources with sizes of up to tens of kpc and 12 compact sources (< a few kpc) in the radio-power range 10 25.7 W Hz –1 –10 26.5 W Hz –1 . We detect H I in five sources, corresponding to a detection rate of ~19%. Within the error bars, this agrees with the detection rate found at lower redshifts. We find that the rest-frame UV luminosities of most of the sources in the sample, including all the detections, are below the proposed threshold above which the H I is supposed to have been ionised. An analysis of the optical emission-line spectra of the sources shows that despite their high radio powers, about one-third of the sample, including two detections, are low-ionisation sources. The radio continuum emission from the sources detected in H I is unresolved at ~5 to 10 kpc scales in our observations, but shows extended structure on parsec scales. We analysed the H I 21 cm absorption spectra of the detections to understand the morphology and kinematics of H I. The absorption profiles are mostly complex with widths between the nulls ranging from ~60 km s –1 to 700 km s –1 . These detections also exhibit remarkably high H I column densities in the range ~10 21 cm –2 –10 22 cm –2 for T spin = 100 K and unit covering factor. Our modelling of the H I 21 cm absorption profiles suggests that in two sources the gas appears to be disturbed, and in three cases, including one with disturbed H I, the majority of the absorption is consistent with it arising from an H I disc. Despite the high radio power of our sources, we do not detect fast outflows. However, the optical emission lines in these detections show the presence of significantly disturbed gas in the nuclear regions in the form of very wide and highly blueshifted emission-line components. Since some of our detections are also low-ionisation AGN, it is quite possible that this disturbance is caused by the radio jets. Overall, our findings point towards a continuation of the low-z trends in the H I detection rate and the incidence of H I in radio AGN up to z ~ 0.4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The kinematic properties of Milky Way stellar halo populations

ABSTRACT In the Gaia era, stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern data sets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ‘action diamond’) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially anisotropic Gaia–Enceladus merger remnant. Assuming a 1:1 ratio of Gaia–Enceladus stars to more isotropic halo, we find our selection achieves a sample purity of 86 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic data sets can noticeably impact the E−Lz distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Chemical Cartography with APOGEE: Mapping Disk Populations with a 2-process Model and Residual Abundances

We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] versus [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core-collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals Δ[X/H] from this two-parameter fit. The rms residuals range from ~0.01–0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow for the comparison of populations controlled for differences in metallicity and [α/Fe]. Relative to the main disk (R = 3–13 kpc), we find nearly identical abundance patterns in the outer disk (R = 15–17 kpc), 0.05–0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4–1 dex) of multiple elements in ω Cen. The residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Scaling pair count to next galaxy surveys

ABSTRACT Counting pairs of galaxies or stars according to their distance is at the core of real-space correlation analyses performed in astrophysics and cosmology. Upcoming galaxy surveys (LSST, Euclid) will measure properties of billions of galaxies challenging our ability to perform such counting in a minute-scale time relevant for the usage of simulations. The problem is only limited by efficient access to the data, hence belongs to the big data category. We use the popular Apache Spark framework to address it and design an efficient high-throughput algorithm to deal with hundreds of millions to billions of input data. To optimize it, we revisit the question of non-hierarchical sphere pixelization based on cube symmetries and develop a new one dubbed the ‘Similar Radius Sphere Pixelization’ (SARSPix) with very close to square pixels. It provides the most adapted indexing over the sphere for all distance-related computations. Using LSST-like fast simulations, we compute autocorrelation functions on tomographic bins containing between a hundred million to one billion data points. In each case, we achieve the construction of a standard pair-distance histogram in about 2 min, using a simple algorithm that is shown to scale, over a moderate number of nodes (16–64). This illustrates the potential of this new techniques in the field of astronomy where data access is becoming the main bottleneck. They can be easily adapted to other use-cases as nearest-neighbours search, catalogue cross-match or cluster finding. The software is publicly available from https://github.com/astrolabsoftware/SparkCorr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Imaging systematics and clustering of DESI main targets

ABSTRACT We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG), and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits $r_0 = 7.78 \pm 0.26\, h^{-1}$Mpc, γ = 1.98 ± 0.02 for LRGs and $r_0 = 5.45 \pm 0.1\, h^{-1}$Mpc, γ = 1.54 ± 0.01 for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale-dependent biases in the weakly non-linear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer dN/dz and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments.

79 ASTRONOMY AND ASTROPHYSICS↗