Engineering Papers⌕ Search

Engineering topics

Padmanabhan, Nikhil

Publications and source records attributed to Padmanabhan, Nikhil.

At least 19 records

Mitigating the noise of DESI mocks using analytic control variates

In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological N-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

First detection of the BAO signal from early DESI data

We present the first detection of the baryon acoustic oscillations (BAOs) signal obtained using unblinded data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). From a selected sample of 261 291 luminous red galaxies spanning the redshift interval 0.4 < z < 1.1 and covering 1651 square degrees with a 57.9 per cent completeness level, we report a ∼5σ level BAO detection and the measurement of the BAO location at a precision of 1.7 per cent. Using a bright galaxy sample of 109 523 galaxies in the redshift range 0.1 < z < 0.5, over 3677 square degrees with a 50.0 per cent completeness, we also detect the BAO feature at ∼3σ significance with a 2.6 per cent precision. These first BAO measurements represent an important milestone, acting as a quality control on the optimal performance of the complex robotically actuated, fibre-fed DESI spectrograph, as well as an early validation of the DESI spectroscopic pipeline and data management system. Based on these first promising results, we forecast that DESI is on target to achieve a high-significance BAO detection at sub-per cent precision with the completed 5-yr survey data, meeting the top-level science requirements on BAO measurements. This exquisite level of precision will set new standards in cosmology and confirm DESI as the most competitive BAO experiment for the remainder of this decade.

79 ASTRONOMY AND ASTROPHYSICS↗

The High Latitude Spectroscopic Survey on the Nancy Grace Roman Space Telescope

The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0–1.93 μm, R = 435–865) and the 0.28 deg 2 wide-field camera. We present a reference HLSS that maps 2000 deg 2 and achieves an emission-line flux limit of 10 -16 erg s -1 cm -2 at 6.5σ, requiring ~0.6 yr of observing time. We summarize the flowdown of the Roman science objectives to the science and technical requirements of the HLSS. We construct a mock redshift survey over the full HLSS volume by applying a semianalytic galaxy formation model to a cosmological N-body simulation and use this mock survey to create pixel-level simulations of 4 deg 2 of HLSS grism spectroscopy. We find that the reference HLSS would measure ~10 million H α galaxy redshifts that densely map large-scale structure at z = 1–2 and 2 million [O III] galaxy redshifts that sparsely map structures at z = 2–3. We forecast the performance of this survey for measurements of the cosmic expansion history with baryon acoustic oscillations and the growth of large-scale structure with redshift-space distortions. We also study possible deviations from the reference design and find that a deep HLSS at f line > 7 × 10 -17 erg s -1 cm -2 over 4000 deg 2 (requiring ~1.5 yr of observing time) provides the most compelling stand-alone constraints on dark energy from Roman alone. This provides a useful reference for future optimizations. The reference survey, simulated data sets, and forecasts presented here will inform community decisions on the final scope and design of the Roman HLSS.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal 1D Ly α forest power spectrum estimation – I. DESI-lite spectra

ABSTRACT The 1D Ly α forest flux power spectrum P1D is sensitive to scales smaller than a typical galaxy survey, and hence ties to the intergalactic medium’s thermal state, suppression from neutrino masses, and new dark matter models. It has emerged as a competitive framework to study new physics, but also has come with various challenges and systematic errors in analysis. In this work, we revisit the optimal quadratic estimator for P1D, which is robust against the relevant problems such as pixel masking, time evolution within spectrum, and quasar continuum errors. We further improve the estimator by introducing a fiducial power spectrum, which enables us to extract more information by alleviating the discreteness of band powers. We meticulously apply our method to synthetic Dark Energy Spectroscopic Instrument (DESI) spectra and demonstrate how the estimator overcomes each challenge. We further apply an optimization scheme that approximates the Fisher matrix to three elements per row and reduces computation time by 60 per cent. We show that we can achieve per cent precision in P1D with 5-yr DESI data in the absence of systematics and provide forecasts for different spectral qualities.

79 ASTRONOMY AND ASTROPHYSICS↗

Astro2020 APC White Paper: The MegaMapper: a z > 2 Spectroscopic Instrument for the Study of Inflation and Dark Energy

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2¡z¡5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Ruler To Measure The Universe : Probing Dark Energy with Baryon Acoustic Oscillations

Understanding the nature of the dark energy that is responsible for roughly 70% of the energy content of the Universe and is causing the accelerated expansion of the Universe is one of the foremost problems in cosmology today. This grant supported the development of the Baryon Acoustic Oscillation (BAO) technique to measure the expansion rate of the Universe and thereby constrain the properties of the dark energy. The primary achievement over the course of the grant was the development of the reconstruction technique for BAO which enhances the constraining power of this method. At the start of this grant, this reconstruction technique, although theoretically studied, had not been applied to data. The work supported by this grant helped make reconstruction a standard technique in the analysis of BAO data. This grant also supported analyses using the Baryon Oscillation Spectroscopic Survey (BOSS, supported by DOE), as well as theoretical work in preparation for the next generations of surveys, most notably the Dark Energy Spectroscopic Instrument (DESI).

79 ASTRONOMY AND ASTROPHYSICS↗

Anatomy of cosmic tidal reconstruction

Abstract 21-cm intensity surveys aim to map neutral hydrogen atoms in the universe through hyper-fine emission. Unfortunately, long-wavelength (low-wavenumber) radial modes are highly contaminated by smooth astrophysical foregrounds that are six orders of magnitude brighter than the cosmological signal. This contamination also leaks into higher radial and angular wavenumber modes and forms a foreground wedge. Cosmic tidal reconstruction aims to extract the large-scale signal from anisotropic features in the local small-scale power spectrum through non-linear tidal interactions; losing small-scale modes to foreground wedge will impair its performance. In this paper, we review tidal interaction theory and estimator construction, and derive the theoretical expressions for the reconstructed spectra. We show the reconstruction is robust against peculiar velocities. Removing low line-of-sight k modes, we demonstrate cross-correlation coefficient r is greater than 0.7 on large scales (k ≲ 0.1 h Mpc−1) even with a cut-off value $k^c_{\Vert }=0.1$h Mpc−1. Discarding wedge modes yields 0.3 ≲ r ≲ 0.5 and completely removes the dependency on $k^c_{\Vert }$. Our theoretical predictions agree with these numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Ensemble photometric redshifts

Upcoming imaging surveys, such as LSST, will provide an unprecedented view of the Universe, but with limited resolution along the line-of-sight. Common ways to increase resolution in the third dimension, and reduce misclassifications, include observing a wider wavelength range and/or combining the broad-band imaging with higher spectral resolution data. The challenge with these approaches is matching the depth of these ancillary data with the original imaging survey. However, while a full 3D map is required for some science, there are many situations where only the statistical distribution of objects (dN/dz) in the line-of-sight direction is needed. In such situations, there is no need to measure the fluxes of individual objects in all of the surveys. Rather a stacking procedure can be used to perform an `ensemble photo-z'. We show how a shallow, higher spectral resolution survey can be used to measure dN/dz for stacks of galaxies which coincide in a deeper, lower resolution survey. The galaxies in the deeper survey do not even need to appear individually in the shallow survey. We give a toy model example to illustrate tradeoffs and considerations for applying this method. This approach will allow deep imaging surveys to leverage the high resolution of spectroscopic and narrow/medium band surveys underway, even when the latter do not have the same reach to high redshift.

galaxies:distances and redshifts↗

The extended Baryon Oscillation Spectroscopic Survey: testing a new approach to measure the evolution of the structure growth

We report that the extended Baryon Oscillation Spectroscopic Survey (eBOSS) is one of the first of a new generation of galaxy redshift surveys that will cover a large range in redshift with sufficient resolution to measure the baryon acoustic oscillations signal. For surveys covering a large redshift range we can no longer ignore cosmological evolution, meaning that either the redshift shells analysed have to be significantly narrower than the survey, or we have to allow for the averaging over evolving quantities. Both of these have the potential to remove signal: analysing small volumes increases the size of the Fourier window function, reducing the large-scale information, while averaging over evolving quantities can, if not performed carefully, remove differential information. It will be important to measure cosmological evolution from these surveys to explore and discriminate between models. We apply a method to optimally extract this differential information to mock catalogues designed to mimic the eBOSS quasar sample. By applying a set of weights to extract redshift-space distortion measurements as a function of redshift, we demonstrate an analysis that does not invoke the problems discussed above. Finally, we show that our estimator gives unbiased constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

The clustering of theSDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: measuring the anisotropic baryon acoustic oscillations with redshift weights

We present an anisotropic analysis of baryon acoustic oscillation (BAO) signal from the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample. The sample consists of 147 000 quasars distributed over a redshift range of 0.8 < z < 2.2. We apply the redshift weights technique to the clustering of quasars in this sample and achieve a 4.6 per cent measurement of the angular distance measurement DM at z = 2.2 and Hubble parameter H at z = 0.8. We parametrize the distance-redshift relation, relative to a fiducial model, as a Taylor series. The coefficients of this expansion are used to reconstruct the distance-redshift relation and obtain distance and Hubble parameter measurements at all redshifts within the redshift range of the sample. Reporting the result at two characteristic redshifts, we determine DM(z = 1) = 3405 ± 305 (rd/rd, fid) Mpc, H(z=1) = 120.7± 7.3 (r_d,fid / r_d) km s^{-1} Mpc^{-1} and DM(z = 2) = 5325 ± 249 (rd/rd, fid) Mpc, H(z=2) = 189.9± 32.9 (r_d,fid / r_d) km s^{-1} Mpc^{-1}. These measurements are highly correlated. We assess the outlook of BAO analysis from the final quasar sample by testing the method on a set of mocks that mimic the noise level in the final sample. We demonstrate on these mocks that redshift weighting shrinks the measurement error by over 25 per cent on average. We conclude redshift weighting can bring us closer to the cosmological goal of the final quasar sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Matched filtering with interferometric 21 cm experiments

A new generation of interferometric instruments is emerging, which aims to use intensity mapping of redshifted 21 cm radiation to measure the large-scale structure of the Universe at z ≃ 1–2 over wide areas of the sky. While these instruments typically have limited angular resolution, they cover huge volumes and thus can be used to provide large samples of rare objects. We study how well such instruments could find spatially extended large-scale structures, such as cosmic voids, using a matched filter formalism. Such a formalism allows us to work in Fourier space, the natural space for interferometers, and to study the impact of finite u - v coverage, noise and foregrounds on our ability to recover voids. We find that in the absence of foregrounds, such instruments would provide enormous catalogs of voids, with high completeness, but that control of foregrounds is key to realizing this goal.

79 ASTRONOMY AND ASTROPHYSICS↗

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample

Here we present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg 2 and volume of 18.7 Gpc 3 , divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51 and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method, in combination with a cosmic microwave background prior on the sound horizon scale, after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DMH from the Alcock–Paczynski (AP) effect and the growth of structure, quantified by fσ8(z), from redshift-space distortions (RSD). We combine individual measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method; in particular, the AP measurement from sub-BAO scales sharpens constraints from post-reconstruction BAOs by breaking degeneracy between DM and H. Combined with Planck 2016 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature ΩK = 0.0003 ± 0.0026 and a dark energy equation-of-state parameter w = -1.01 ± 0.06, in strong affirmation of the spatially flat cold dark matter (CDM) model with a cosmological constant (ΛCDM). Our RSD measurements of fσ8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3 ± 1.0 km s -1 Mpc -1 even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8 ± 1.2km s -1 Mpc -1 . Assuming flat ΛCDM, we find Ωm = 0.310 ± 0.005 and H0 = 67.6 ± 0.5 km s -1 Mpc -1 , and we find a 95 per cent upper limit of 0.16 eV c -2 on the neutrino mass sum.

79 ASTRONOMY AND ASTROPHYSICS↗

Mitigating the impact of the DESI fiber assignment on galaxy clustering

In this paper, we present a simple strategy to mitigate the impact of an incomplete spectroscopic redshift galaxy sample as a result of fiber assignment and survey tiling. The method has been designed for the Dark Energy Spectroscopic Instrument (DESI) galaxy survey but may have applications beyond this. We propose a modification to the usual correlation function that nulls the almost purely angular modes affected by survey incompleteness due to fiber assignment. Predictions of this modified statistic can be calculated given a model of the two point correlation function. The new statistic can be computed with a slight modification to the data catalogues input to the standard correlation function code and does not incur any additional computational time. Lastly, we show that the spherically averaged baryon acoustic oscillation signal is not biased by the new statistic.

79 ASTRONOMY AND ASTROPHYSICS↗

The Clustering of Galaxies in the Completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmic Flows and Cosmic Web from Luminous Red Galaxies

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the Sloan Digital Sky Survey-III Baryon Oscillations Spectroscopic Survey Data Release 12 CMASS galaxy clustering catalogue. We rely on a given Λ cold dark matter cosmology, a mesh resolution in the range of 6-10 h -1 Mpc, and a lognormal-Poisson model with a redshift-dependent non-linear bias. The bias parameters are derived from the data and a general renormalized perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the ARGO code, to iteratively reconstruct the darkmatter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions. Our tests relying on accurate N-body-based mock galaxy catalogues show unbiased real space power spectra of the non-linear density field up to k ~ 0.2 h Mpc -1 , and vanishing quadrupoles down to r ~ 20 h -1 Mpc. We also demonstrate that the non-linear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of r ~ 0.68 including about 10 per cent of satellite galaxies with virial motions (about r = 0.75 without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to k ~ 0.2 h Mpc -1 . This work will be especially useful to improve, for example, baryon acoustic oscillation reconstructions, kinematic Sunyaev- Zeldovich, integrated Sachs-Wolfe measurements or environmental studies.

Ata, Metin↗

High Energy Physics Exascale Requirements Review. An Office of Science review sponsored jointly by Advanced Scientific Computing Research and High Energy Physics, June 10-12, 2015, Bethesda, Maryland

The U.S. Department of Energy (DOE) Office of Science (SC) Offices of High Energy Physics (HEP) and Advanced Scientific Computing Research (ASCR) convened a programmatic Exascale Requirements Review on June 10–12, 2015, in Bethesda, Maryland. This report summarizes the findings, results, and recommendations derived from that meeting. The high-level findings and observations are as follows. Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude — and in some cases greater — than that available currently. The growth rate of data produced by simulations is overwhelming the current ability of both facilities and researchers to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. Data rates and volumes from experimental facilities are also straining the current HEP infrastructure in its ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. A close integration of high-performance computing (HPC) simulation and data analysis will greatly aid in interpreting the results of HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. Long-range planning between HEP and ASCR will be required to meet HEP’s research needs. To best use ASCR HPC resources, the experimental HEP program needs (1) an established, long-term plan for access to ASCR computational and data resources, (2) the ability to map workflows to HPC resources, (3) the ability for ASCR facilities to accommodate workflows run by collaborations potentially comprising thousands of individual members, (4) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, (5) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.

97 MATHEMATICS AND COMPUTING↗