Showing 1–50 of 145 results for author: Szapudi, I

Euclid: ERO -- NISP-only sources and the search for luminous $z=6-8$ galaxies

Authors: J. R. Weaver, S. Taamoli, C. J. R. McPartland, L. Zalesky, N. Allen, S. Toft, D. B. Sanders, H. Atek, R. A. A. Bowler, D. Stern, C. J. Conselice, B. Mobasher, I. Szapudi, P. R. M. Eisenhardt, G. Murphree, I. Valdes, K. Ito, S. Belladitta, P. A. Oesch, S. Serjeant, D. J. Mortlock, N. A. Hatch, M. Kluge, B. Milvang-Jensen, G. Rodighiero , et al. (163 additional authors not shown)

Abstract: This paper presents a search for high redshift galaxies from the Euclid Early Release Observations program "Magnifying Lens." The 1.5 deg$^2$ area covered by the twin Abell lensing cluster fields is comparable in size to the few other deep near-infrared surveys such as COSMOS, and so provides an opportunity to significantly increase known samples of rare UV-bright galaxies at $z\approx6-8$ (… ▽ More This paper presents a search for high redshift galaxies from the Euclid Early Release Observations program "Magnifying Lens." The 1.5 deg$^2$ area covered by the twin Abell lensing cluster fields is comparable in size to the few other deep near-infrared surveys such as COSMOS, and so provides an opportunity to significantly increase known samples of rare UV-bright galaxies at $z\approx6-8$ ($M_{\rm UV}\lesssim-22$). Beyond their still uncertain role in reionisation, these UV-bright galaxies are ideal laboratories from which to study galaxy formation and constrain the bright-end of the UV luminosity function. Of the 501994 sources detected from a combined $Y_{\rm E}$, $J_{\rm E}$, and $H_{\rm E}$ NISP detection image, 168 do not have any appreciable VIS/$I_{\rm E}$ flux. These objects span a range in spectral colours, separated into two classes: 139 extremely red sources; and 29 Lyman-break galaxy candidates. Best-fit redshifts and spectral templates suggest the former is composed of both $z\gtrsim5$ dusty star-forming galaxies and $z\approx1-3$ quiescent systems. The latter is composed of more homogeneous Lyman break galaxies at $z\approx6-8$. In both cases, contamination by L- and T-type dwarfs cannot be ruled out with Euclid images alone. Additional contamination from instrumental persistence is investigated using a novel time series analysis. This work lays the foundation for future searches within the Euclid Deep Fields, where thousands more $z\gtrsim6$ Lyman break systems and extremely red sources will be identified. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 22 pages, 13 figures, paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

The CMB lensing imprint of cosmic voids detected in the WISE-Pan-STARRS luminous red galaxy catalog

Authors: G. Camacho-Ciurana, P. Lee, N. Arsenov, A. Kovács, I. Szapudi, I. Csabai

Abstract: The cross-correlation of cosmic voids with the lensing convergence ($κ$) map of the CMB fluctuations offers a powerful tool to refine our understanding of the dark sector in the consensus cosmological model. Our principal aim is to compare the lensing signature of our galaxy data set with simulations based on the concordance model and characterize the results with an $A_κ$ consistency parameter. I… ▽ More The cross-correlation of cosmic voids with the lensing convergence ($κ$) map of the CMB fluctuations offers a powerful tool to refine our understanding of the dark sector in the consensus cosmological model. Our principal aim is to compare the lensing signature of our galaxy data set with simulations based on the concordance model and characterize the results with an $A_κ$ consistency parameter. In particular, our measurements contribute to the understanding of the "lensing-is-low" tension of the $Λ$CDM model. We selected luminous red galaxies from the WISE-Pan-STARSS data set, allowing an extended 14,200 deg$^2$ sky area, that offers a more precise measurement compared to previous studies. We created 2D and 3D void catalogs to cross-correlate their locations with the Planck lensing map and studied their average imprint signal using a stacking methodology. Applying the same procedure, we also generated a mock catalog from the WebSky simulation for comparison. The 2D void analysis revealed good agreement with the standard cosmological model with $A_κ\approx1.06 \pm 0.08$, i.e. $S/N=13.3$, showing a higher $S/N$ than previous studies using voids detected in the Dark Energy Survey data set. The 3D void analysis exhibited a lower $S/N$ and demonstrated worse agreement with our mock catalog than the 2D voids. These deviations might be attributed to limitations in the mock catalog, such as imperfections in the LRG selection, as well as a potential asymmetry between the North and South patches of the WISE-Pan-STARSS data set in terms of data quality. Overall, we present a significant detection of a CMB lensing signal associated with cosmic voids, largely consistent with the concordance model. Future analyses using even larger data sets also hold great promise of further sharpening these results, given their complementary nature to large-scale structure analyses. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 10 pages, 6 figure, submitted to A&A

Evolution of matter and galaxy clustering in cosmological hydrodynamical simulations

Authors: Jaan Einasto, Gert Hütsi, Lauri-Juhan Liivamägi, Changbom Park, Juhan Kim, Istval Szapudi, Maret Einasto

Abstract: We quantify the evolution of matter and galaxy clustering in cosmological hydrodynamical simulations via correlation and bias functions of matter and galaxies. We use simulations TNG100 and TNG300 with epochs from $z=5$ to $z=0$. We calculate spatial correlation functions of galaxies, $ξ(r)$, for simulated galaxies and dark matter (DM) particles to characterise the evolving cosmic web. We find tha… ▽ More We quantify the evolution of matter and galaxy clustering in cosmological hydrodynamical simulations via correlation and bias functions of matter and galaxies. We use simulations TNG100 and TNG300 with epochs from $z=5$ to $z=0$. We calculate spatial correlation functions of galaxies, $ξ(r)$, for simulated galaxies and dark matter (DM) particles to characterise the evolving cosmic web. We find that bias parameters decrease during the evolution, confirming earlier results. At low and medium luminosities, bias parameters of galaxies, $b_0$, are equal, suggesting that dwarf galaxies reside in the same filamentary web as brighter galaxies. Bias parameters of the lowest luminosity galaxies estimated from CFs are lower relative to CFs of particle density-limited clustered samples of DM. We find that bias parameters $b_0$, estimated from CFs of clustered DM, agree with the expected values from the fraction of particles in the clustered population, $b=1/F_c$. The cosmic web contains filamentary structures of various densities, and fractions of matter in the clustered and the unclustered populations are both less than unity. Thus the CF amplitude of the clustered matter is always higher than for all matter, i.e. bias parameter must be $b>1$. Differences between CFs of galaxies and clustered DM suggest that these functions describe different properties of the cosmic web. △ Less

Submitted 4 June, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 15 pages, 10 figures, accepted in Monthly Notices of Royal Astronomical Society

Introducing the Texas Euclid Survey for Lyman Alpha (TESLA) Survey: Initial Study Correlating Galaxy Properties to Lyman-Alpha Emission

Authors: Oscar A. Chavez Ortiz, Steven L. Finkelstein, Dustin Davis, Gene Leung, Erin Mentuch Cooper, Micaela Bagley, Rebecca Larson, Caitlin M. Casey, Adam P. McCarron, Karl Gebhardt, Yuchen Guo, Chenxu Liu, Isaac Laseter, Jason Rhodes, Ralf Bender, Max Fabricius, Ariel G. Sanchez, Claudia Scarlata, Peter Capak, David Sanders, Istvan Szapudi, Eric Baxter, Conor McPartland, John R. Weaver, Sune Toft , et al. (2 additional authors not shown)

Abstract: We present the Texas Euclid Survey for Lyman-Alpha (TESLA), a spectroscopic survey in the 10 square degree of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyman-alpha emitters (LAEs) correlate with Lyman-alpha emission to understand the escape of Lyman alpha from galaxies at redshifts 2 -- 3.5. We present an analysis of 43 LAEs performed in the N… ▽ More We present the Texas Euclid Survey for Lyman-Alpha (TESLA), a spectroscopic survey in the 10 square degree of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyman-alpha emitters (LAEs) correlate with Lyman-alpha emission to understand the escape of Lyman alpha from galaxies at redshifts 2 -- 3.5. We present an analysis of 43 LAEs performed in the NEP field using early data from the TESLA survey. We use Subaru Hyper Suprime-Cam imaging in the grizy-bands, Spitzer/IRAC channels 1 and 2 from the Hawaii 20 square degree (H20) survey and spectra acquired by the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby-Eberly Telescope. We perform spectral energy distribution (SED) fitting to compute the galaxy properties of 43 LAEs, and study correlations between stellar mass, star formation rate (SFR), and dust, to the Lyman-alpha rest-frame equivalent widths (EW). We uncover marginal (1 sigma significance) correlations between stellar mass and Lyman-alpha EW, and star formation rate (SFR) and Lyman-alpha EW, with a Spearman correlation coefficient of -0.$34_{-.14}^{+.17}$ and -0.$37_{-.14}^{+.16}$ respectively. We show that the Lyman-alpha distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of 150 Angstrom. Once complete the TESLA survey will enable the study of ~ thousands of LAEs to explore correlations between galaxy properties and Lyman-alpha EW. The large sample size will allow the construction of a predictive model for the Lyman-alpha EW as a function of SED-derived galaxy properties, which could be used to improve Lyman-alpha based constraints on reionization. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Complementary Cosmological Simulations

Authors: Gábor Rácz, Alina Kiessling, István Csabai, István Szapudi

Abstract: Cosmic variance limits the accuracy of cosmological N-body simulations, introducing bias in statistics such as the power spectrum, halo mass function, or the cosmic shear. We provide new methods to measure and reduce the effect of cosmic variance in existing and new simulations. We ran pairs of simulations using phase-shifted initial conditions with matching amplitudes. We set the initial amplitud… ▽ More Cosmic variance limits the accuracy of cosmological N-body simulations, introducing bias in statistics such as the power spectrum, halo mass function, or the cosmic shear. We provide new methods to measure and reduce the effect of cosmic variance in existing and new simulations. We ran pairs of simulations using phase-shifted initial conditions with matching amplitudes. We set the initial amplitudes of the Fourier modes to ensure that the average power spectrum of the pair is equal to the cosmic mean power spectrum from linear theory. The average power spectrum of a pair of such simulations remains consistent with the estimated nonlinear spectra of the state-of-the-art methods even at late times. We also show that the effect of cosmic variance on any analysis involving a cosmological simulation can be estimated using the complementary pair of the original simulation. To demonstrate the effectiveness of our novel technique, we simulated a complementary pair of the original Millennium run and quantified the degree to which cosmic variance affected its the power spectrum. The average power spectrum of the original and complementary Millennium simulation was able to directly resolve the baryon acoustic oscillation features. △ Less

Submitted 17 February, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 7 pages, 4 figures, submitted to Astronomy & Astrophysics (A&A)

Journal ref: A&A 672, A59 (2023)

A Machine Learning Approach to Predict Missing Flux Densities in Multi-band Galaxy Surveys

Authors: Nima Chartab, Bahram Mobasher, Asantha Cooray, Shoubaneh Hemmati, Zahra Sattari, Henry C. Ferguson, David B. Sanders, John R. Weaver, Daniel Stern, Henry J. McCracken, Daniel C. Masters, Sune Toft, Peter L. Capak, Iary Davidzon, Mark Dickinson, Jason Rhodes, Andrea Moneti, Olivier Ilbert, Lukas Zalesky, Conor McPartland, Istvan Szapudi, Anton M. Koekemoer, Harry I. Teplitz, Mauro Giavalisco

Abstract: We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with a desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wavebands for measuring the physical properties of galaxies in a Hawaii Two-0 (H20)- and UVISTA-like survey for a sample o… ▽ More We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with a desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wavebands for measuring the physical properties of galaxies in a Hawaii Two-0 (H20)- and UVISTA-like survey for a sample of $i<25$ AB mag galaxies. We find that with available $i$-band fluxes, $r$, $u$, IRAC/$ch2$ and $z$ bands provide most of the information regarding the redshift with importance decreasing from $r$-band to $z$-band. We also find that for the same sample, IRAC/$ch2$, $Y$, $r$ and $u$ bands are the most relevant bands in stellar mass measurements with decreasing order of importance. Investigating the inter-correlation between the bands, we train a model to predict UVISTA observations in near-IR from H20-like observations. We find that magnitudes in $YJH$ bands can be simulated/predicted with an accuracy of $1σ$ mag scatter $\lesssim 0.2$ for galaxies brighter than 24 AB mag in near-IR bands. One should note that these conclusions depend on the selection criteria of the sample. For any new sample of galaxies with a different selection, these results should be remeasured. Our results suggest that in the presence of a limited number of bands, a machine learning model trained over the population of observed galaxies with extensive spectral coverage outperforms template-fitting. Such a machine learning model maximally comprises the information acquired over available extensive surveys and breaks degeneracies in the parameter space of template-fitting inevitable in the presence of a few bands. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: 15 pages, 14 figures, accepted for publication in ApJ

Photometric redshifts for quasars from WISE-PS1-STRM

Authors: Sándor Kunsági-Máté, Róbert Beck, István Szapudi, István Csabai

Abstract: Three-dimensional wide-field galaxy surveys are fundamental for cosmological studies. For higher redshifts (z > 1.0), where galaxies are too faint, quasars still trace the large-scale structure of the Universe. Since available telescope time limits spectroscopic surveys, photometric methods are efficient for estimating redshifts for many quasars. Recently, machine learning methods are increasingly… ▽ More Three-dimensional wide-field galaxy surveys are fundamental for cosmological studies. For higher redshifts (z > 1.0), where galaxies are too faint, quasars still trace the large-scale structure of the Universe. Since available telescope time limits spectroscopic surveys, photometric methods are efficient for estimating redshifts for many quasars. Recently, machine learning methods are increasingly successful for quasar photometric redshifts, however, they hinge on the distribution of the training set. Therefore a rigorous estimation of reliability is critical. We extracted optical and infrared photometric data from the cross-matched catalogue of the WISE All-Sky and PS1 3$π$ DR2 sky surveys. We trained an XGBoost regressor and an artificial neural network on the relation between color indices and spectroscopic redshift. We approximated the effective training set coverage with the K nearest neighbors algorithm. We estimated reliable photometric redshifts of 2,879,298 quasars which overlap with the training set in feature space. We validated the derived redshifts with an independent, clustering-based redshift estimation technique. The final catalog is publicly available. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

Authors: Elcio Abdalla, Guillermo Franco Abellán, Amin Aboubrahim, Adriano Agnello, Ozgur Akarsu, Yashar Akrami, George Alestas, Daniel Aloni, Luca Amendola, Luis A. Anchordoqui, Richard I. Anderson, Nikki Arendse, Marika Asgari, Mario Ballardini, Vernon Barger, Spyros Basilakos, Ronaldo C. Batista, Elia S. Battistelli, Richard Battye, Micol Benetti, David Benisty, Asher Berlin, Paolo de Bernardis, Emanuele Berti, Bohdan Bidenko , et al. (178 additional authors not shown)

Abstract: In this paper we will list a few important goals that need to be addressed in the next decade, also 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 system… ▽ More In this paper we will list a few important goals that need to be addressed in the next decade, also 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 {\it 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 {\it 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.[Abridged] △ Less

Submitted 24 April, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

Journal ref: J. High En. Astrophys. 2204, 002 (2022)

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

Authors: Andrea Moneti, H. J. McCracken, M. Shuntov, O. B. Kauffmann, P. Capak, I. Davidzon, O. Ilbert, C. Scarlata, S. Toft, J. Weaver, R. Chary, J. Cuby, A. L. Faisst, D. C. Masters, C. McPartland, B. Mobasher, D. B. Sanders, R. Scaramella, D. Stern, I. Szapudi, H. Teplitz, L. Zalesky, A. Amara, N. Auricchio, C. Bodendorf , et al. (172 additional authors not shown)

Abstract: We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission tim… ▽ More We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 15 pages with 11 figures, approved by Euclid Consortium Publication Board and submitted to Astronomy and Astrophysics. Data products will become available via the IRSA website once the paper is accepted. This paper is a companion to "COSMOS2020: A panchromatic view of the Universe to z~10 from two complementary catalogs" by John Weaver et al., which is being posted in parallel

Journal ref: A&A 658, A126 (2022)

Indicator Power Spectra: Surgical Excision of Non-linearities and Covariance Matrices for Counts in Cells

Authors: Andrew Repp, István Szapudi

Abstract: We here introduce indicator functions, which identify regions of a given density in order to characterize the density dependence of clustering. After a general introduction to this tool, we show that indicator-function power spectra are biased versions of the linear spectrum on large scales. We provide a calculation from first principles for this bias, we show that it reproduces simulation results… ▽ More We here introduce indicator functions, which identify regions of a given density in order to characterize the density dependence of clustering. After a general introduction to this tool, we show that indicator-function power spectra are biased versions of the linear spectrum on large scales. We provide a calculation from first principles for this bias, we show that it reproduces simulation results, and we provide a simple functional form for the translinear portion of the indicator-function spectra. We also outline two applications: first, these spectra facilitate surgical excision of non-linearity and thus significantly increase the reach of linear theory. Second, indicator-function spectra permit calculation of theoretical covariance matrices for counts-in-cells (CIC), facilitating parameter estimation with complementary CIC methods. △ Less

Submitted 3 November, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: Updated to final published version (relatively minor changes; information added to Figure 2

Evidence for a high-z ISW signal from supervoids in the distribution of eBOSS quasars

Authors: A. Kovács, R. Beck, A. Smith, G. Rácz, I. Csabai, I. Szapudi

Abstract: The late-time integrated Sachs-Wolfe (ISW) imprint of $R\gtrsim 100~h^{-1}{\rm Mpc}$ super-structures is sourced by evolving large-scale potentials due to a dominant dark energy component in the $Λ$CDM model. The aspect that makes the ISW effect distinctly interesting is the repeated observation of stronger-than-expected imprints from supervoids at $z\lesssim0.9$. Here we analyze the un-probed key… ▽ More The late-time integrated Sachs-Wolfe (ISW) imprint of $R\gtrsim 100~h^{-1}{\rm Mpc}$ super-structures is sourced by evolving large-scale potentials due to a dominant dark energy component in the $Λ$CDM model. The aspect that makes the ISW effect distinctly interesting is the repeated observation of stronger-than-expected imprints from supervoids at $z\lesssim0.9$. Here we analyze the un-probed key redshift range $0.8<z<2.2$ where the ISW signal is expected to fade in $Λ$CDM, due to a weakening dark energy component, and eventually become consistent with zero in the matter dominated epoch. On the contrary, alternative cosmological models, proposed to explain the excess low-$z$ ISW signals, predicted a sign-change in the ISW effect at $z\approx1.5$ due to the possible growth of large-scale potentials that is absent in the standard model. To discriminate, we estimated the high-$z$ $Λ$CDM ISW signal using the Millennium XXL mock catalogue, and compared it to our measurements from about 800 supervoids identified in the eBOSS DR16 quasar catalogue. At $0.8<z<1.2$, we found an excess ISW signal with $A_\mathrm{ ISW}\approx3.6\pm2.1$ amplitude. The signal is then consistent with the $Λ$CDM expectation ($A_\mathrm{ ISW}=1$) at $1.2<z<1.5$ where the standard and alternative models predict similar amplitudes. Most interestingly, we also detected an opposite-sign ISW signal at $1.5<z<2.2$ that is in $2.7σ$ tension with the $Λ$CDM prediction. Taken at face value, these moderately significant detections of ISW anomalies suggest an alternative growth rate of structure in low-density environments at $\sim100~h^{-1}{\rm Mpc}$ scales. △ Less

Submitted 18 April, 2022; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 13 pages, 9 figures, accepted by MNRAS after minor corrections to text and some figures

An empirical nonlinear power spectrum overdensity-response

Authors: Gábor Rácz, István Szapudi, István Csabai

Abstract: Context. The overdensity inside a cosmological sub-volume and the tidal fields from its surroundings affect the matter distribution of the region. The resulting difference between the local and global power spectra is characterized by the response function. Aims. Our aim is to provide a new, simple, and accurate formula for the power spectrum overdensity response at highly nonlinear scales based… ▽ More Context. The overdensity inside a cosmological sub-volume and the tidal fields from its surroundings affect the matter distribution of the region. The resulting difference between the local and global power spectra is characterized by the response function. Aims. Our aim is to provide a new, simple, and accurate formula for the power spectrum overdensity response at highly nonlinear scales based on the results of cosmological simulations and paying special attention to the lognormal nature of the density field. Methods. We measured the dark matter power spectrum amplitude as a function of the overdensity ($δ_W$) in $N$-body simulation subsamples. We show that the response follows a power-law form in terms of $(1+δ_W)$, and we provide a new fit in terms of the variance, $σ(L)$, of a sub-volume of size $L$. Results. Our fit has a similar accuracy and a comparable complexity to second-order standard perturbation theory on large scales, but it is also valid for nonlinear (smaller) scales, where perturbation theory needs higher-order terms for a comparable precision. Furthermore, we show that the lognormal nature of the overdensity distribution causes a previously unidentified bias: the power spectrum amplitude for a subsample with an average density is typically underestimated by about $-2σ^2$. Although this bias falls to the sub-percent level above characteristic scales of $200Mpch^{-1}$, taking it into account improves the accuracy of estimating power spectra from zoom-in simulations and smaller high-resolution surveys embedded in larger low-resolution volumes. △ Less

Submitted 11 March, 2022; v1 submitted 30 May, 2021; originally announced May 2021.

Journal ref: A&A 661, A96 (2022)

Constraining Mach's principle with high precision astrometry

Authors: Istvan Szapudi

Abstract: The analyses of high precision astrometric surveys, such as Gaia, implicitly assume a modern version of Mach's Principle: the local inertial frame of our Solar System should be non-rotating in the frame of distant quasars. On the contrary, Einstein's General Relativity allows a rotating universe. Thus, relaxing the assumption of Mach's Principle will allow placing a constraint on a class of rotati… ▽ More The analyses of high precision astrometric surveys, such as Gaia, implicitly assume a modern version of Mach's Principle: the local inertial frame of our Solar System should be non-rotating in the frame of distant quasars. On the contrary, Einstein's General Relativity allows a rotating universe. Thus, relaxing the assumption of Mach's Principle will allow placing a constraint on a class of rotating cosmologies by comparing high precision astrometry of quasars with well-measured solar system orbits. Constraining global rotation will test General Relativity, inflation, and the isotropy of cosmological initial conditions. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 4 pages, Essay written for the Gravity Research Foundation 2021 Awards for Essays on Gravitation

Intensity correlation speckles as a technique for removing Doppler broadening

Authors: R. Merlin, N. Green, I. Szapudi, G. Tarle

Abstract: A method involving intensity correlation measurements is described, which allows for the complete removal of Doppler broadening in the emission of electromagnetic radiation from far-away sources that are inaccessible to conventional Doppler-free measurements. The technique, relying on a correction to g(2) of order N-1, probes the separation between neighboring spectral lines and is also applicable… ▽ More A method involving intensity correlation measurements is described, which allows for the complete removal of Doppler broadening in the emission of electromagnetic radiation from far-away sources that are inaccessible to conventional Doppler-free measurements. The technique, relying on a correction to g(2) of order N-1, probes the separation between neighboring spectral lines and is also applicable to the elimination of broadening due to collisions (N is the number of emitting particles and g(2) is the second-order field correlation function). Possible applications include a determination of cosmological parameters from red shifts of gravitationally-lensed quasars. △ Less

Submitted 20 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. A 103, L041701(2021)

The Variance and Covariance of Counts-in-Cells Probabilities

Authors: Andrew Repp, István Szapudi

Abstract: Counts-in-cells (CIC) measurements contain a wealth of cosmological information yet are seldom used to constrain theories. Although we can predict the shape of the distribution for a given cosmology, to fit a model to the observed CIC probabilities requires the covariance matrix -- both the variance of counts in one probability bin and the covariance between counts in different bins. To date, ther… ▽ More Counts-in-cells (CIC) measurements contain a wealth of cosmological information yet are seldom used to constrain theories. Although we can predict the shape of the distribution for a given cosmology, to fit a model to the observed CIC probabilities requires the covariance matrix -- both the variance of counts in one probability bin and the covariance between counts in different bins. To date, there have been no general expressions for these variances. Here we show that correlations of particular levels, or "slices," of the density field determine the variance and covariance of CIC probabilities. We derive explicit formulae that accurately predict the variance and covariance among subvolumes of a simulated galaxy catalog, opening the door to the use of CIC measurements for cosmological parameter estimation. △ Less

Submitted 30 June, 2020; originally announced July 2020.

Comments: 10 pages, 4 figures; submitted to MNRAS

The anisotropy of the power spectrum in periodic cosmological simulations

Authors: Gábor Rácz, István Szapudi, István Csabai, László Dobos

Abstract: The classical gravitational force on a torus is anisotropic and always lower than Newton's $1/r^2$ law. We demonstrate the effects of periodicity in dark matter only $N$-body simulations of spherical collapse and standard $Λ$CDM initial conditions. Periodic boundary conditions cause an overall negative and anisotropic bias in cosmological simulations of cosmic structure formation. The lower amplit… ▽ More The classical gravitational force on a torus is anisotropic and always lower than Newton's $1/r^2$ law. We demonstrate the effects of periodicity in dark matter only $N$-body simulations of spherical collapse and standard $Λ$CDM initial conditions. Periodic boundary conditions cause an overall negative and anisotropic bias in cosmological simulations of cosmic structure formation. The lower amplitude of power spectra of small periodic simulations are a consequence of the missing large scale modes and the equally important smaller periodic forces. The effect is most significant when the largest mildly non-linear scales are comparable to the linear size of the simulation box, as often is the case for high-resolution hydrodynamical simulations. Spherical collapse morphs into a shape similar to an octahedron. The anisotropic growth distorts the large-scale $Λ$CDM dark matter structures. We introduce the direction-dependent power spectrum invariant under the octahedral group of the simulation volume and show that the results break spherical symmetry. △ Less

Submitted 23 March, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Galaxy Bias and $σ_8$ from Counts in Cells from the SDSS Main Sample

Authors: Andrew Repp, István Szapudi

Abstract: The counts-in-cells (CIC) galaxy probability distribution depends on both the dark matter clustering amplitude $σ_8$ and the galaxy bias $b$. We present a theory for the CIC distribution based on a previous prescription of the underlying dark matter distribution and a linear volume transformation to redshift space. We show that, unlike the power spectrum, the CIC distribution breaks the degeneracy… ▽ More The counts-in-cells (CIC) galaxy probability distribution depends on both the dark matter clustering amplitude $σ_8$ and the galaxy bias $b$. We present a theory for the CIC distribution based on a previous prescription of the underlying dark matter distribution and a linear volume transformation to redshift space. We show that, unlike the power spectrum, the CIC distribution breaks the degeneracy between $σ_8$ and $b$ on scales large enough that both bias and redshift distortions are still linear; thus we obtain a simultaneous fit for both parameters. We first validate the technique on the Millennium Simulation and then apply it to the SDSS Main Galaxy Sample. We find $σ_8 = 0.94^{+.11}_{-.10}$ and $b = 1.36^{+.14}_{-.11}$, consistent with previous complementary results from redshift distortions and from Planck. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Comments: 5 pages, 3 figures; submitted to MNRAS

A common explanation of the Hubble tension and anomalous cold spots in the CMB

Authors: András Kovács, Róbert Beck, István Szapudi, István Csabai, Gábor Rácz, László Dobos

Abstract: The standard cosmological paradigm narrates a reassuring story of a universe currently dominated by an enigmatic dark energy component. Disquietingly, its universal explaining power has recently been challenged by, above all, the $\sim4σ$ tension in the values of the Hubble constant. Another, less studied anomaly is the repeated observation of integrated Sachs-Wolfe imprints $\sim5\times$ stronger… ▽ More The standard cosmological paradigm narrates a reassuring story of a universe currently dominated by an enigmatic dark energy component. Disquietingly, its universal explaining power has recently been challenged by, above all, the $\sim4σ$ tension in the values of the Hubble constant. Another, less studied anomaly is the repeated observation of integrated Sachs-Wolfe imprints $\sim5\times$ stronger than expected in the $Λ$CDM model from R>100 $Mpc/h$ super-structures. Here we show that the inhomogeneous AvERA model of emerging curvature is capable of telling a plausible albeit radically different story that explains both observational anomalies without dark energy. We demonstrate that while stacked imprints of R>100 $Mpc/h$ supervoids in cosmic microwave background temperature maps can discriminate between the AvERA and $Λ$CDM models, their characteristic differences may remain hidden using alternative void definitions and stacking methodologies. Testing the extremes, we then also show that the CMB Cold Spot can plausibly be explained in the AvERA model as an ISW imprint. The coldest spot in the AvERA map is aligned with multiple low-$z$ supervoids with R>100 $Mpc/h$ and central underdensity $δ_{0}\approx-0.3$, resembling the observed large-scale galaxy density field in the Cold Spot area. We hence conclude that the anomalous imprint of supervoids may well be the canary in the coal mine, and existing observational evidence for dark energy should be re-interpreted to further test alternative models. △ Less

Submitted 7 September, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures, accepted by MNRAS

Does Gravity Fall Down? Evidence for Gravitational Wave Deflection Along the Line of Sight to GW 170817

Authors: David Rubin, Istvan Szapudi, Benjamin J. Shappee, Gagandeep S. Anand

Abstract: We present a novel test of general relativity (GR): measuring the geometric component of the time delay due to gravitational lensing. GR predicts that photons and gravitational waves follow the same geodesic paths and thus experience the same geometric time delay. We show that for typical systems, the time delays are tens of seconds, and thus can dominate over astrophysical delays in the timing of… ▽ More We present a novel test of general relativity (GR): measuring the geometric component of the time delay due to gravitational lensing. GR predicts that photons and gravitational waves follow the same geodesic paths and thus experience the same geometric time delay. We show that for typical systems, the time delays are tens of seconds, and thus can dominate over astrophysical delays in the timing of photon emission. For the case of GW 170817, we use a multi-plane lensing code to evaluate the time delay due to four massive halos along the line of sight. From literature mass and distance measurements of these halos, we establish at high confidence (significantly greater than 5 sigma) that the gravitational waves of GW 170817 underwent gravitational deflection to arrive within 1.7 seconds of the photons. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Comments: Under review with ApJL

Empirical Validation of the Ising Galaxy Bias Model

Authors: Andrew Repp, István Szapudi

Abstract: Repp and Szapudi (2019) present a physically-motivated galaxy bias model which remains physical in low-density regions and which also provides a better fit to simulation data than do typical survey-analysis bias models. Given plausible simplifying assumptions, the physics of this model (surprisingly) proves to be analogous to the Ising model of statistical mechanics. In the present work we present… ▽ More Repp and Szapudi (2019) present a physically-motivated galaxy bias model which remains physical in low-density regions and which also provides a better fit to simulation data than do typical survey-analysis bias models. Given plausible simplifying assumptions, the physics of this model (surprisingly) proves to be analogous to the Ising model of statistical mechanics. In the present work we present a method of testing this Ising bias model against empirical galaxy survey data. Using this method, we compare our model (as well as three reference models -- linear, quadratic, and logarithmic) to SDSS, 6dFGS, and COSMOS2015 results, finding that for spectroscopic redshift surveys, the Ising bias model provides a superior fit compared to the reference models. Photometric redshifts, on the other hand, introduce enough error into the radial coordinate that none of the models yields a good fit. A physically meaningful galaxy bias model is necessary for optimal extraction of cosmological information from dense galaxy surveys such as Euclid and WFIRST. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 22 pages, 16 figures; submitted to MNRAS

PS1-STRM: Neural network source classification and photometric redshift catalogue for PS1 $3π$ DR1

Authors: Róbert Beck, István Szapudi, Heather Flewelling, Conrad Holmberg, Eugene Magnier

Abstract: The Pan-STARRS1 (PS1) $3π$ survey is a comprehensive optical imaging survey of three quarters of the sky in the $grizy$ broad-band photometric filters. We present the methodology used in assembling the source classification and photometric redshift (photo-z) catalogue for PS1 $3π$ Data Release 1, titled Pan-STARRS1 Source Types and Redshifts with Machine learning (PS1-STRM). For both main data p… ▽ More The Pan-STARRS1 (PS1) $3π$ survey is a comprehensive optical imaging survey of three quarters of the sky in the $grizy$ broad-band photometric filters. We present the methodology used in assembling the source classification and photometric redshift (photo-z) catalogue for PS1 $3π$ Data Release 1, titled Pan-STARRS1 Source Types and Redshifts with Machine learning (PS1-STRM). For both main data products, we use neural network architectures, trained on a compilation of public spectroscopic measurements that has been cross-matched with PS1 sources. We quantify the parameter space coverage of our training data set, and flag extrapolation using self-organizing maps. We perform a Monte-Carlo sampling of the photometry to estimate photo-z uncertainty. The final catalogue contains $2,902,054,648$ objects. On our validation data set, for non-extrapolated sources, we achieve an overall classification accuracy of $98.1\%$ for galaxies, $97.8\%$ for stars, and $96.6\%$ for quasars. Regarding the galaxy photo-z estimation, we attain an overall bias of $\left<Δz_{\mathrm{norm}}\right>=0.0005$, a standard deviation of $σ(Δz_{\mathrm{norm}})=0.0322$, a median absolute deviation of $\mathrm{MAD}(Δz_{\mathrm{norm}})=0.0161$, and an outlier fraction of $O=1.89\%$. The catalogue will be made available as a high-level science product via the Mikulski Archive for Space Telescopes at https://doi.org/10.17909//t9-rnk7-gr88. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 12 pages, 6 figures. Submitted to MNRAS

Hawaii Two-0: High-redshift galaxy clustering and bias

Authors: Róbert Beck, Conor McPartland, Andrew Repp, David Sanders, István Szapudi

Abstract: We perform an analysis of two-point galaxy clustering and galaxy bias using Subaru Hyper-Suprime Cam (HSC) data taken jointly by the Subaru Strategic Program and the University of Hawaii in the COSMOS field. The depth of the data is similar to the ongoing Hawaii Two-0 (H20) optical galaxy survey, thus the results are indicative of future constraints from tenfold area. We measure the angular auto… ▽ More We perform an analysis of two-point galaxy clustering and galaxy bias using Subaru Hyper-Suprime Cam (HSC) data taken jointly by the Subaru Strategic Program and the University of Hawaii in the COSMOS field. The depth of the data is similar to the ongoing Hawaii Two-0 (H20) optical galaxy survey, thus the results are indicative of future constraints from tenfold area. We measure the angular auto-power spectra of the galaxy overdensity in three redshift bins, defined by dropouts from the g-, r- and i-bands, and compare them to the theoretical expectation from concordance cosmology with linear galaxy bias. We determine the redshift distribution of each bin using a standard template-based photometric redshift method, coupled with a self-organizing map (SOM) to quantify colour space coverage. We also investigate sources of systematic errors to inform the methodology and requirements for Hawaii Two-0. The linear galaxy bias fit results are $b_{\mathrm{gal,g}} = 3.90 \pm 0.33 (\mathrm{stat}) \substack{ +0.64 \\ -0.24 } (\mathrm{sys})$ at redshift $z \simeq 3.7$, $b_{\mathrm{gal,r}} = 8.44 \pm 0.63 (\mathrm{stat}) \substack{ +1.42 \\ -0.72 } (\mathrm{sys})$ at $z \simeq 4.7$, and $b_{\mathrm{gal,i}} = 11.94 \pm 2.24 (\mathrm{stat}) \substack{ +1.82 \\ -1.27 } (\mathrm{sys})$ at $z \simeq 5.9$. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 12 pages, 8 figures. Submitted to MNRAS

An Ising model for galaxy bias

Authors: Andrew Repp, István Szapudi

Abstract: A reliable model of galaxy bias is necessary for interpreting data from future dense galaxy surveys. Conventional bias models are inaccurate, in that they can yield unphysical results ($δ_g < -1$) for voids that might contain half of the available cosmological information. For this reason, we present a physically-motivated bias model based on an analogy with the Ising model. With only two free par… ▽ More A reliable model of galaxy bias is necessary for interpreting data from future dense galaxy surveys. Conventional bias models are inaccurate, in that they can yield unphysical results ($δ_g < -1$) for voids that might contain half of the available cosmological information. For this reason, we present a physically-motivated bias model based on an analogy with the Ising model. With only two free parameters, the model produces sensible results for both high- and low-density regions. We also test the model using a catalog of Millennium Simulation galaxies in cubical survey pixels with side lengths from $2h^{-1}$--$31h^{-1}$Mpc, at redshifts from 0 to 2. We find the Ising model markedly superior to linear and quadratic bias models on scales smaller than $10h^{-1}$Mpc, while those conventional models fare better on scales larger than $30h^{-1}$Mpc. While the largest scale where the Ising model is applicable might vary for a specific galaxy catalog, it should be superior on any scale with a non-negligible fraction of cells devoid of galaxies. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Comments: 13 pages, 8 figures; submitted to MNRAS

A Gravitational Ising Model for the Statistical Bias of Galaxies

Authors: Andrew Repp, István Szapudi

Abstract: Evaluation of gravitational theories by means of cosmological data suffers from the fact that galaxies are biased tracers of dark matter. Current bias models focus primarily on high-density regions, whereas low-density regions carry significant amounts of information relevant to the constraint of dark energy and alternative gravity theories. Thus, proper treatment of both high and low densities is… ▽ More Evaluation of gravitational theories by means of cosmological data suffers from the fact that galaxies are biased tracers of dark matter. Current bias models focus primarily on high-density regions, whereas low-density regions carry significant amounts of information relevant to the constraint of dark energy and alternative gravity theories. Thus, proper treatment of both high and low densities is important for future surveys. Accordingly, we here present an interactionless Ising model for this bias, and we demonstrate that it exhibits a remarkably good fit to both Millennium Simulation and Sloan Digital Sky Survey data, at both density extremes. The quality of the fit indicates that galaxy formation is (to first order) an essentially local process determined by initial conditions. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: Essay written for the Gravity Research Foundation 2019 Awards for Essays on Gravitation; 7 pages, 2 figures

Cosmic voids: a novel probe to shed light on our Universe

Authors: Alice Pisani, Elena Massara, David N. Spergel, David Alonso, Tessa Baker, Yan-Chuan Cai, Marius Cautun, Christopher Davies, Vasiliy Demchenko, Olivier Doré, Andy Goulding, Mélanie Habouzit, Nico Hamaus, Adam Hawken, Christopher M. Hirata, Shirley Ho, Bhuvnesh Jain, Christina D. Kreisch, Federico Marulli, Nelson Padilla, Giorgia Pollina, Martin Sahlén, Ravi K. Sheth, Rachel Somerville, Istvan Szapudi , et al. (4 additional authors not shown)

Abstract: Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust… ▽ More Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White Paper for the Astro2020 Decadal Survey, 5 pages text, 14 pages with references and affiliations

StePS: A Multi-GPU Cosmological N-body Code for Compactified Simulations

Authors: Gábor Rácz, István Szapudi, László Dobos, István Csabai, Alexander S. Szalay

Abstract: We present the multi-GPU realization of the StePS (Stereographically Projected Cosmological Simulations) algorithm with MPI-OpenMP-CUDA hybrid parallelization and nearly ideal scale-out to multiple compute nodes. Our new zoom-in cosmological direct N-body simulation method simulates the infinite universe with unprecedented dynamic range for a given amount of memory and, in contrast to traditional… ▽ More We present the multi-GPU realization of the StePS (Stereographically Projected Cosmological Simulations) algorithm with MPI-OpenMP-CUDA hybrid parallelization and nearly ideal scale-out to multiple compute nodes. Our new zoom-in cosmological direct N-body simulation method simulates the infinite universe with unprecedented dynamic range for a given amount of memory and, in contrast to traditional periodic simulations, its fundamental geometry and topology match observations. By using a spherical geometry instead of periodic boundary conditions, and gradually decreasing the mass resolution with radius, our code is capable of running simulations with a few gigaparsecs in diameter and with a mass resolution of $\sim 10^{9}M_{\odot}$ in the center in four days on three compute nodes with four GTX 1080Ti GPUs in each. The code can also be used to run extremely fast simulations with reasonable resolution for fitting cosmological parameters. These simulations are useful for prediction needs of large surveys. The StePS code is publicly available for the research community. △ Less

Submitted 21 March, 2019; v1 submitted 14 November, 2018; originally announced November 2018.

Predicting the Sufficient-Statistics Power Spectrum for Galaxy Surveys: A Recipe for $P_{A*}(k)$

Authors: Andrew Repp, István Szapudi

Abstract: Future galaxy surveys hope to realize significantly tighter constraints on various cosmological parameters. The higher number densities achieved by these surveys will allow them to probe the smaller scales affected by non-linear clustering. However, in these regimes, the standard power spectrum can extract only a portion of such surveys' cosmological information. In contrast, the alternate statist… ▽ More Future galaxy surveys hope to realize significantly tighter constraints on various cosmological parameters. The higher number densities achieved by these surveys will allow them to probe the smaller scales affected by non-linear clustering. However, in these regimes, the standard power spectrum can extract only a portion of such surveys' cosmological information. In contrast, the alternate statistic $A^*$ has the potential to double these surveys' information return, provided one can predict the $A^*$-power spectrum for a given cosmology. Thus, in this work we provide a prescription for this power spectrum $P_{A^*}(k)$, finding that the prescription is typically accurate to about 5 per cent for near-concordance cosmologies. This prescription will thus allow us to multiply the information gained from surveys such as Euclid and WFIRST. △ Less

Submitted 11 April, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: Minor clarifications and corrections in response to the referee's comments. 16 pages, 9 figures; published in MNRAS

Journal ref: MNRAS (2019), 483 4438R

The integrated Sachs-Wolfe effect in the AvERA cosmology

Authors: Róbert Beck, István Csabai, Gábor Rácz, István Szapudi

Abstract: The recent AvERA cosmological simulation of Rácz et al. (2017) has a $Λ\mathrm{CDM}$-like expansion history and removes the tension between local and Planck (cosmic microwave background) Hubble constants. We contrast the AvERA prediction of the integrated Sachs--Wolfe (ISW) effect with that of $Λ\mathrm{CDM}$. The linear ISW effect is proportional to the derivative of the growth function, thus it… ▽ More The recent AvERA cosmological simulation of Rácz et al. (2017) has a $Λ\mathrm{CDM}$-like expansion history and removes the tension between local and Planck (cosmic microwave background) Hubble constants. We contrast the AvERA prediction of the integrated Sachs--Wolfe (ISW) effect with that of $Λ\mathrm{CDM}$. The linear ISW effect is proportional to the derivative of the growth function, thus it is sensitive to small differences in the expansion histories of the respective models. We create simulated ISW maps tracing the path of light-rays through the Millennium XXL cosmological simulation, and perform theoretical calculations of the ISW power spectrum. AvERA predicts a significantly higher ISW effect than $Λ\mathrm{CDM}$, $A=1.93-5.29$ times larger depending on the $l$ index of the spherical power spectrum, which could be utilized to definitively differentiate the models. We also show that AvERA predicts an opposite-sign ISW effect in the redshift range $z \approx 1.5 - 4.4$, in clear contrast with $Λ\mathrm{CDM}$. Finally, we compare our ISW predictions with previous observations. While at present these cannot distinguish between the two models due to large error bars, and lack of internal consistency suggesting systematics, ISW probes from future surveys will tightly constrain the models. △ Less

Submitted 25 January, 2018; originally announced January 2018.

Comments: 9 pages, 8 figures. Submitted to MNRAS

Compactified Cosmological Simulations of the Infinite Universe

Authors: Gábor Rácz, István Szapudi, István Csabai, László Dobos

Abstract: We present a novel $N$-body simulation method that compactifies the infinite spatial extent of the Universe into a finite sphere with isotropic boundary conditions to follow the evolution of the large-scale structure. Our approach eliminates the need for periodic boundary conditions, a mere numerical convenience which is not supported by observation and which modifies the law of force on large sca… ▽ More We present a novel $N$-body simulation method that compactifies the infinite spatial extent of the Universe into a finite sphere with isotropic boundary conditions to follow the evolution of the large-scale structure. Our approach eliminates the need for periodic boundary conditions, a mere numerical convenience which is not supported by observation and which modifies the law of force on large scales in an unrealistic fashion. We demonstrate that our method outclasses standard simulations executed on workstation-scale hardware in dynamic range, it is balanced in following a comparable number of high and low $k$ modes and, its fundamental geometry and topology match observations. Our approach is also capable of simulating an expanding, infinite universe in static coordinates with Newtonian dynamics. The price of these achievements is that most of the simulated volume has smoothly varying mass and spatial resolution, an approximation that carries different systematics than periodic simulations. Our initial implementation of the method is called StePS which stands for Stereographically Projected Cosmological Simulations. It uses stereographic projection for space compactification and naive $\mathcal{O}(N^2)$ force calculation which is nevertheless faster to arrive at a correlation function of the same quality than any standard (tree or P$^3$M) algorithm with similar spatial and mass resolution. The $N^2$ force calculation is easy to adapt to modern graphics cards, hence our code can function as a high-speed prediction tool for modern large-scale surveys. To learn about the limits of the respective methods, we compare StePS with GADGET-2 \citep{Gadget2_2005MNRAS.364.1105S} running matching initial conditions. △ Less

Submitted 15 February, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

The Bias of the Log Power Spectrum for Discrete Surveys

Authors: Andrew Repp, István Szapudi

Abstract: A primary goal of galaxy surveys is to tighten constraints on cosmological parameters, and the power spectrum $P(k)$ is the standard means of doing so. However, at translinear scales $P(k)$ is blind to much of these surveys' information---information which the log density power spectrum recovers. For discrete fields (such as the galaxy density), $A^*$ denotes the statistic analogous to the log den… ▽ More A primary goal of galaxy surveys is to tighten constraints on cosmological parameters, and the power spectrum $P(k)$ is the standard means of doing so. However, at translinear scales $P(k)$ is blind to much of these surveys' information---information which the log density power spectrum recovers. For discrete fields (such as the galaxy density), $A^*$ denotes the statistic analogous to the log density: $A^*$ is a "sufficient statistic" in that its power spectrum (and mean) capture virtually all of a discrete survey's information. However, the power spectrum of $A^*$ is biased with respect to the corresponding log spectrum for continuous fields, and to use $P_{A^*}(k)$ to constrain the values of cosmological parameters, we require some means of predicting this bias. Here we present a prescription for doing so; for Euclid-like surveys (with cubical cells 16$h^{-1}$ Mpc across) our bias prescription's error is less than 3 per cent. This prediction will facilitate optimal utilization of the information in future galaxy surveys. △ Less

Submitted 27 July, 2018; v1 submitted 2 August, 2017; originally announced August 2017.

Comments: 5 pages, 4 figures; published in MNRASL (here updated to published version)

Journal ref: MNRAS 475, L6-L10 (2018)

Precision Prediction for the Cosmological Density Distribution

Authors: Andrew Repp, István Szapudi

Abstract: The distribution of matter in the universe is, to first order, lognormal. Improving this approximation requires characterization of the third moment (skewness) of the log density field. Thus, using Millennium Simulation phenomenology and building on previous work, we present analytic fits for the mean, variance, and skewness of the log density field $A$. We further show that a Generalized Extreme… ▽ More The distribution of matter in the universe is, to first order, lognormal. Improving this approximation requires characterization of the third moment (skewness) of the log density field. Thus, using Millennium Simulation phenomenology and building on previous work, we present analytic fits for the mean, variance, and skewness of the log density field $A$. We further show that a Generalized Extreme Value (GEV) distribution accurately models $A$; we submit that this GEV behavior is the result of strong intrapixel correlations, without which the smoothed distribution would tend (by the Central Limit Theorem) toward a Gaussian. Our GEV model yields cumulative distribution functions accurate to within 1.7 per cent for near-concordance cosmologies, over a range of redshifts and smoothing scales. △ Less

Submitted 31 July, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 11 pages, 11 figures; published in MNRAS; final revision to post published version

Journal ref: MNRAS 473, 3598-3607 (2018)

Evidence against a supervoid causing the CMB Cold Spot

Authors: Ruari Mackenzie, Tom Shanks, Malcolm N. Bremer, Yan-Chuan Cai, Madusha L. P. Gunawardhana, András Kovács, Peder Norberg, Istvan Szapudi

Abstract: We report the results of the 2dF-VST ATLAS Cold Spot galaxy redshift survey (2CSz) based on imaging from VST ATLAS and spectroscopy from 2dF AAOmega over the core of the CMB Cold Spot. We sparsely surveyed the inner 5$^{\circ}$ radius of the Cold Spot to a limit of $i_{AB} \le 19.2$, sampling $\sim7000$ galaxies at $z<0.4$. We have found voids at $z=$ 0.14, 0.26 and 0.30 but they are interspersed… ▽ More We report the results of the 2dF-VST ATLAS Cold Spot galaxy redshift survey (2CSz) based on imaging from VST ATLAS and spectroscopy from 2dF AAOmega over the core of the CMB Cold Spot. We sparsely surveyed the inner 5$^{\circ}$ radius of the Cold Spot to a limit of $i_{AB} \le 19.2$, sampling $\sim7000$ galaxies at $z<0.4$. We have found voids at $z=$ 0.14, 0.26 and 0.30 but they are interspersed with small over-densities and the scale of these voids is insufficient to explain the Cold Spot through the $Λ$CDM ISW effect. Combining with previous data out to $z\sim1$, we conclude that the CMB Cold Spot could not have been imprinted by a void confined to the inner core of the Cold Spot. Additionally we find that our 'control' field GAMA G23 shows a similarity in its galaxy redshift distribution to the Cold Spot. Since the GAMA G23 line-of-sight shows no evidence of a CMB temperature decrement we conclude that the Cold Spot may have a primordial origin rather than being due to line-of-sight effects. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: Submitted to MNRAS

Density-dependent clustering: I. Pulling back the curtains on motions of the BAO peak

Authors: Mark C. Neyrinck, István Szapudi, Nuala McCullagh, Alex Szalay, Bridget Falck, Jie Wang

Abstract: The most common statistic used to analyze large-scale structure surveys is the correlation function, or power spectrum. Here, we show how `slicing' the correlation function on local density brings sensitivity to interesting non-Gaussian features in the large-scale structure, such as the expansion or contraction of baryon acoustic oscillations (BAO) according to the local density. The sliced correl… ▽ More The most common statistic used to analyze large-scale structure surveys is the correlation function, or power spectrum. Here, we show how `slicing' the correlation function on local density brings sensitivity to interesting non-Gaussian features in the large-scale structure, such as the expansion or contraction of baryon acoustic oscillations (BAO) according to the local density. The sliced correlation function measures the large-scale flows that smear out the BAO, instead of just correcting them as reconstruction algorithms do. Thus, we expect the sliced correlation function to be useful in constraining the growth factor, and modified gravity theories that involve the local density. Out of the studied cases, we find that the run of the BAO peak location with density is best revealed when slicing on a $\sim 40$ Mpc/$h$ filtered density. But slicing on a $\sim100$ Mpc/$h$ filtered density may be most useful in distinguishing between underdense and overdense regions, whose BAO peaks are separated by a substantial $\sim 5$ Mpc/$h$ at $z=0$. We also introduce `curtain plots' showing how local densities drive particle motions toward or away from each other over the course of an $N$-body simulation. △ Less

Submitted 14 May, 2018; v1 submitted 19 October, 2016; originally announced October 2016.

Comments: Matches revision accepted in MNRAS

The lensing and temperature imprints of voids on the Cosmic Microwave Background

Authors: Yan-Chuan Cai, Mark Neyrinck, Qingqing Mao, John A. Peacock, Istvan Szapudi, Andreas A. Berlind

Abstract: We have searched for the signature of cosmic voids in the CMB, in both the Planck temperature and lensing-convergence maps; voids should give decrements in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our analysis on N-body simulations, to avoid a posteriori bias. For the first time, we detect the signature of voids in CMB lensing: the significance is $3.2σ$, close to… ▽ More We have searched for the signature of cosmic voids in the CMB, in both the Planck temperature and lensing-convergence maps; voids should give decrements in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our analysis on N-body simulations, to avoid a posteriori bias. For the first time, we detect the signature of voids in CMB lensing: the significance is $3.2σ$, close to $Λ$CDM in both amplitude and projected density-profile shape. A temperature dip is also seen, at modest significance ($2.3σ$), with amplitude about 6 times the prediction. This temperature signal is induced mostly by voids with radius between 100 and 150 Mpc/h, while the lensing signal is mostly contributed by smaller voids -- as expected; lensing relates directly to density, while ISW depends on gravitational potential. The void abundance in observations and simulations agree, as well. We also repeated the analysis excluding lower-significance voids: no lensing signal is detected, with an upper limit of about twice the $Λ$CDM prediction. But the mean temperature decrement now becomes non-zero at the $3.7σ$ level (similar to that found by Granett et al.), with amplitude about 20 times the prediction. However, the observed dependence of temperature on void size is in poor agreement with simulations, whereas the lensing results are consistent with $Λ$CDM theory. Thus, the overall tension between theory and observations does not favour non-standard theories of gravity, despite the hints of an enhanced amplitude for the ISW effect from voids. △ Less

Submitted 16 December, 2016; v1 submitted 1 September, 2016; originally announced September 2016.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS

Concordance cosmology without dark energy

Authors: Gábor Rácz, László Dobos, Róbert Beck, István Szapudi, István Csabai

Abstract: According to the separate universe conjecture, spherically symmetric sub-regions in an isotropic universe behave like mini-universes with their own cosmological parameters. This is an excellent approximation in both Newtonian and general relativistic theories. We estimate local expansion rates for a large number of such regions, and use a scale parameter calculated from the volume-averaged increme… ▽ More According to the separate universe conjecture, spherically symmetric sub-regions in an isotropic universe behave like mini-universes with their own cosmological parameters. This is an excellent approximation in both Newtonian and general relativistic theories. We estimate local expansion rates for a large number of such regions, and use a scale parameter calculated from the volume-averaged increments of local scale parameters at each time step in an otherwise standard cosmological $N$-body simulation. The particle mass, corresponding to a coarse graining scale, is an adjustable parameter. This mean field approximation neglects tidal forces and boundary effects, but it is the first step towards a non-perturbative statistical estimation of the effect of non-linear evolution of structure on the expansion rate. Using our algorithm, a simulation with an initial $Ω_m=1$ Einstein--de~Sitter setting closely tracks the expansion and structure growth history of the $Λ$CDM cosmology. Due to small but characteristic differences, our model can be distinguished from the $Λ$CDM model by future precision observations. Moreover, our model can resolve the emerging tension between local Hubble constant measurements and the Planck best-fitting cosmology. Further improvements to the simulation are necessary to investigate light propagation and confirm full consistency with cosmic microwave background observations. △ Less

Submitted 12 February, 2017; v1 submitted 29 July, 2016; originally announced July 2016.

Comments: 5 pages, 4 figures, 1 table

Journal ref: Mon Not R Astron Soc Lett 2017 slx026

Precision Prediction of the Log Power Spectrum

Authors: Andrew Repp, István Szapudi

Abstract: At translinear scales, the log power spectrum captures significantly more cosmological information than the standard power spectrum. At high wavenumbers $k$, the Fisher information in the standard power spectrum $P(k)$ fails to increase in proportion to $k$ in part due to correlations between large- and small-scale modes. As a result, $P(k)$ suffers from an information plateau on these translinear… ▽ More At translinear scales, the log power spectrum captures significantly more cosmological information than the standard power spectrum. At high wavenumbers $k$, the Fisher information in the standard power spectrum $P(k)$ fails to increase in proportion to $k$ in part due to correlations between large- and small-scale modes. As a result, $P(k)$ suffers from an information plateau on these translinear scales, so that analysis with the standard power spectrum cannot access the information contained in these small-scale modes. The log power spectrum $P_A(k)$, on the other hand, captures the majority of this otherwise lost information. Until now there has been no means of predicting the amplitude of the log power spectrum apart from cataloging the results of simulations. We here present a cosmology-independent prescription for the log power spectrum; this prescription displays accuracy comparable to that of Smith et al. (2003), over a range of redshifts and smoothing scales, and for wavenumbers up to $1.5h$ Mpc$^{-1}$. △ Less

Submitted 17 December, 2016; v1 submitted 5 July, 2016; originally announced July 2016.

Comments: Updated to published version; 5 pages, 3 figures

Journal ref: MNRAS Letters 2016 464 (1): L21-L25

The Rapid Transient Surveyor

Authors: Christoph Baranec, Jessica R. Lu, Shelley A. Wright, John Tonry, R. Brent Tully, István Szapudi, Marianne Takamiya, Lisa Hunter, Reed Riddle, Shaojie Chen, Mark Chun

Abstract: The Rapid Transient Surveyor (RTS) is a proposed rapid-response, high-cadence adaptive optics (AO) facility for the UH 2.2-m telescope on Maunakea. RTS will uniquely address the need for high-acuity and sensitive near-infrared spectral follow-up observations of tens of thousands of objects in mere months by combining an excellent observing site, unmatched robotic observational efficiency, and an A… ▽ More The Rapid Transient Surveyor (RTS) is a proposed rapid-response, high-cadence adaptive optics (AO) facility for the UH 2.2-m telescope on Maunakea. RTS will uniquely address the need for high-acuity and sensitive near-infrared spectral follow-up observations of tens of thousands of objects in mere months by combining an excellent observing site, unmatched robotic observational efficiency, and an AO system that significantly increases both sensitivity and spatial resolving power. We will initially use RTS to obtain the infrared spectra of ~4,000 Type Ia supernovae identified by the Asteroid Terrestrial-Impact Last Alert System over a two year period that will be crucial to precisely measuring distances and mapping the distribution of dark matter in the z < 0.1 universe. RTS will comprise an upgraded version of the Robo-AO laser AO system and will respond quickly to target-of-opportunity events, minimizing the time between discovery and characterization. RTS will acquire simultaneous-multicolor images with an acuity of 0.07-0.10" across the entire visible spectrum (20% i'-band Strehl in median conditions) and <0.16" in the near infrared, and will detect companions at 0.5" at contrast ratio of ~500. The system will include a high-efficiency prism integral field unit spectrograph: R = 70-140 over a total bandpass of 840-1830 nm with an 8.7" by 6.0" field of view (0.15" spaxels). The AO correction boosts the infrared point-source sensitivity of the spectrograph against the sky background by a factor of seven for faint targets, giving the UH 2.2-m the H-band sensitivity of a 5.7-m telescope without AO. △ Less

Submitted 23 June, 2016; originally announced June 2016.

Comments: 15 pages, 10 figures, 3 tables

Journal ref: Proc. SPIE 9909, Adaptive Optics Systems V, 9909-15 (2016)

What does the N-point function hierarchy of the cosmological matter density field really measure ?

Authors: Julien Carron, István Szapudi

Abstract: The cosmological dark matter field is not completely described by its hierarchy of $N$-point functions, a non-perturbative effect with the consequence that only part of the theory can be probed with the hierarchy. We give here an exact characterization of the joint information of the full set of $N$-point correlators of the lognormal field. The lognormal field is the archetypal example of a field… ▽ More The cosmological dark matter field is not completely described by its hierarchy of $N$-point functions, a non-perturbative effect with the consequence that only part of the theory can be probed with the hierarchy. We give here an exact characterization of the joint information of the full set of $N$-point correlators of the lognormal field. The lognormal field is the archetypal example of a field where this effect occurs, and, at the same time, one of the few tractable and insightful available models to specify fully the statistical properties of the evolved matter density field beyond the perturbative regime. Nonlinear growth in the Universe in that model is set letting the log-density field probability density functional evolve keeping its Gaussian shape, according to the diffusion equation in Euclidean space. We show that the hierarchy probes a different evolution equation, the diffusion equation defined not in Euclidean space but on the compact torus, with uniformity as the long-term solution. The extraction of the hierarchy of correlators can be recast in the form of a nonlinear transformation applied to the field, 'wrapping', undergoing a sharp transition towards complete disorder in the deeply nonlinear regime, where all memory of the initial conditions is lost. △ Less

Submitted 27 April, 2017; v1 submitted 19 August, 2015; originally announced August 2015.

Comments: 4 pages, matches version accepted for publication by MNRAS

The Impact of Non-Gaussianity upon Cosmological Forecasts

Authors: Andrew Repp, István Szapudi, Julien Carron, Melody Wolk

Abstract: The primary science driver for 3D galaxy surveys is their potential to constrain cosmological parameters. Forecasts of these surveys' effectiveness typically assume Gaussian statistics for the underlying matter density, despite the fact that the actual distribution is decidedly non-Gaussian. To quantify the effect of this assumption, we employ an analytic expression for the power spectrum covarian… ▽ More The primary science driver for 3D galaxy surveys is their potential to constrain cosmological parameters. Forecasts of these surveys' effectiveness typically assume Gaussian statistics for the underlying matter density, despite the fact that the actual distribution is decidedly non-Gaussian. To quantify the effect of this assumption, we employ an analytic expression for the power spectrum covariance matrix to calculate the Fisher information for BAO-type model surveys. We find that for typical number densities, at $k_\mathrm{max} = 0.5 h$ Mpc$^{-1}$, Gaussian assumptions significantly overestimate the information on all parameters considered, in some cases by up to an order of magnitude. However, after marginalizing over a six-parameter set, the form of the covariance matrix (dictated by $N$-body simulations) causes the majority of the effect to shift to the "amplitude-like" parameters, leaving the others virtually unaffected. We find that Gaussian assumptions at such wavenumbers can underestimate the dark energy parameter errors by well over 50 per cent, producing dark energy figures of merit almost 3 times too large. Thus, for 3D galaxy surveys probing the non-linear regime, proper consideration of non-Gaussian effects is essential. △ Less

Submitted 24 October, 2015; v1 submitted 30 May, 2015; originally announced June 2015.

Comments: 9 pages, 5 figures, submitted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 454 (4): 3533-3541

Improving constraints on the neutrino mass using sufficient statistics

Authors: M. Wolk, I. Szapudi, J. Bel, C. Carbone, J. Carron

Abstract: We use the "Dark Energy and Massive Neutrino Universe" (DEMNUni) simulations to compare the constraining power of "sufficient statistics" with the standard matter power spectrum on the sum of neutrino masses, $M_ν\equiv \sum m_ν$. In general, the power spectrum, even supplemented with higher moments of the distribution, captures only a fraction of the available cosmological information due to corr… ▽ More We use the "Dark Energy and Massive Neutrino Universe" (DEMNUni) simulations to compare the constraining power of "sufficient statistics" with the standard matter power spectrum on the sum of neutrino masses, $M_ν\equiv \sum m_ν$. In general, the power spectrum, even supplemented with higher moments of the distribution, captures only a fraction of the available cosmological information due to correlations between the Fourier modes. In contrast, the non-linear transform of sufficient statistics, approximated by a logarithmic mapping A=ln(1+δ), was designed to capture all the available cosmological information contained in the matter clustering; in this sense it is an optimal observable. Our analysis takes advantage of the recent analytical model developed by Carron et al. 2014 to estimate both the matter power spectrum and the A-power spectrum covariance matrices. Using a Fisher information approach, we find that using sufficient statistics increases up to 8 times the available information on the total neutrino mass at z=0, thus tightening the constraints by almost a factor of 3 compared to the matter power spectrum. △ Less

Submitted 31 March, 2015; originally announced April 2015.

Comments: 5 pages, 2 figures, 2 tables

On the total cosmological information in galaxy clustering: an analytical approach

Authors: M. Wolk, J. Carron, I. Szapudi

Abstract: Beyond the linear regime of structure formation, part of cosmological information encoded in galaxy clustering becomes inaccessible to the usual power spectrum. "Sufficient statistics", A*, were introduced recently to recapture the lost, and ultimately extract all, cosmological information. We present analytical approximations for the A* and traditional power spectra as well as for their covarianc… ▽ More Beyond the linear regime of structure formation, part of cosmological information encoded in galaxy clustering becomes inaccessible to the usual power spectrum. "Sufficient statistics", A*, were introduced recently to recapture the lost, and ultimately extract all, cosmological information. We present analytical approximations for the A* and traditional power spectra as well as for their covariance matrices in order to calculate analytically their cosmological information content in the context of Fisher information theory. Our approach allows the precise quantitative comparison of the techniques with each other and to the total information in the data, and provides insights into sufficient statistics. In particular, we find that while the A* power spectrum has a similar shape to the usual galaxy power spectrum, its amplitude is strongly modulated by small scale statistics. This effect is mostly responsible for the ability of the A* power spectrum to recapture the information lost for the usual power spectrum. We use our framework to forecast the best achievable cosmological constraints for projected surveys as a function of their galaxy density, and compare the information content of the two power spectra. We find that sufficient statistics extract all cosmological information, resulting in an approximately factor of ~2 gain for dense projected surveys at low redshift. This increase in the effective volume of projected surveys is consistent with previous numerical calculations. △ Less

Submitted 16 March, 2015; originally announced March 2015.

Comments: 9 pages, 7 figures, 1 table

The Subaru FMOS Galaxy Redshift Survey (FastSound). I. Overview of the Survey Targeting on H$α$ Emitters at $z \sim 1.4$

Authors: Motonari Tonegawa, Tomonori Totani, Hiroyuki Okada, Masayuki Akiyama, Gavin Dalton, Karl Glazebrook, Fumihide Iwamuro, Toshinori Maihara, Kouji Ohta, Ikkoh Shimizu, Naruhisa Takato, Naoyuki Tamura, Kiyoto Yabe, Andrew J. Bunker, Jean Coupon, Pedro G. Ferreira, Carlos S. Frenk, Tomotsugu Goto, Chiaki Hikage, Takashi Ishikawa, Takahiko Matsubara, Surhud More, Teppei Okumura, Will J. Percival, Lee R. Spitler , et al. (1 additional authors not shown)

Abstract: FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$α$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$α$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the… ▽ More FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$α$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$α$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$α$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections. △ Less

Submitted 4 May, 2015; v1 submitted 27 February, 2015; originally announced February 2015.

Comments: 12 pages, 9 figures, and 2 tables, accepted for PASJ

On the information content of the matter power spectrum

Authors: Julien Carron, Melody Wolk, István Szapudi

Abstract: We discuss an analytical approximation for the matter power spectrum covariance matrix and its inverse on translinear scales, $k \sim 0.1h - 0.8h/\textrm{Mpc}$ at $z = 0$. We proceed to give an analytical expression for the Fisher information matrix of the nonlinear density field spectrum, and derive implications for its cosmological information content. We find that the spectrum information is ch… ▽ More We discuss an analytical approximation for the matter power spectrum covariance matrix and its inverse on translinear scales, $k \sim 0.1h - 0.8h/\textrm{Mpc}$ at $z = 0$. We proceed to give an analytical expression for the Fisher information matrix of the nonlinear density field spectrum, and derive implications for its cosmological information content. We find that the spectrum information is characterized by a pair of upper bounds, 'plateaux', caused by the trispectrum, and a 'knee' in the presence of white noise. The effective number of Fourier modes, normally growing as a power law, is bounded from above by these plateaux, explaining naturally earlier findings from $N$-body simulations. These plateaux limit best possible measurements of the nonlinear power at the percent level in a $h^{-3}\textrm{Gpc}^3$ volume; the extraction of model parameters from the spectrum is limited explicitly by their degeneracy to the nonlinear amplitude. The value of the first, super-survey (SS) plateau depends on the characteristic survey volume and the large scale power; the second, intra-survey (IS) plateau is set by the small scale power. While both have simple interpretations within the hierarchical \textit{Ansatz}, the SS plateau can be predicted and generalized to still smaller scales within Takada and Hu's spectrum response formalism. Finally, the noise knee is naturally set by the density of tracers. △ Less

Submitted 26 July, 2015; v1 submitted 17 December, 2014; originally announced December 2014.

Comments: 6 pages, 1 fig, matches accepted version

Unveiling the cosmological information beyond linear scales: forecasts for sufficient statistics

Authors: M. Wolk, J. Carron, I. Szapudi

Abstract: Beyond the linear regime, Fourier modes of cosmological random fields become correlated, and the power spectrum of density fluctuations contains only a fraction of the available cosmological information. To unveil this formerly hidden information, the A* non-linear transform was introduced; it is optimized both for the nonlinearities induced by gravity and observational noise. Quantifying the resu… ▽ More Beyond the linear regime, Fourier modes of cosmological random fields become correlated, and the power spectrum of density fluctuations contains only a fraction of the available cosmological information. To unveil this formerly hidden information, the A* non-linear transform was introduced; it is optimized both for the nonlinearities induced by gravity and observational noise. Quantifying the resulting increase of our knowledge of cosmological parameters, we forecast the constraints from the angular power spectrum and that of A* from l ~ 200 to 3000 for upcoming galaxy surveys such as: the Wide-Field Infrared Survey Telescope (WFIRST), the Large Synoptic Survey Telescope (LSST), Euclid, the Hyper Suprime-Cam (HSC) and the Dark Energy Survey (DES). We find that at low redshifts this new data analysis strategy can double the extracted information, effectively doubling the survey area. To test the accuracy of our forecasting and the power of our data analysis methods, we apply the A* transformation to the latest release of the Canada-France-Hawaii-Telescope Legacy Survey (CFHTLS) Wide. While this data set is too sparse to allow for more than modest gains (~1.1-1.2), the realized gain from our method is in excellent agreement with our forecast, thus verifying the robustness of our analysis and prediction pipelines. △ Less

Submitted 22 November, 2014; originally announced November 2014.

Comments: 11, 7 figures, 6 tables

The impact of super-survey modes on cosmological constraints from cosmic shear fields

Authors: Julien Carron, István Szapudi

Abstract: Owing to the mass-sheet degeneracy, cosmic shear maps do not probe directly the Fourier modes of the underlying mass distribution on scales comparable to the survey size and larger. To assess the corresponding effect on attainable cosmological parameter constraints, we quantify the information on super-survey modes in a lognormal model and, when interpreted as nuisance parameters, their degeneraci… ▽ More Owing to the mass-sheet degeneracy, cosmic shear maps do not probe directly the Fourier modes of the underlying mass distribution on scales comparable to the survey size and larger. To assess the corresponding effect on attainable cosmological parameter constraints, we quantify the information on super-survey modes in a lognormal model and, when interpreted as nuisance parameters, their degeneracies to cosmological parameters. Our analytical and numerical calculations clarify the central role of super-sample covariance (SSC) in shaping the statistical power of cosmological observables. Reconstructing the background modes from their non-Gaussian statistical dependence to small scales modes yields the renormalized convergence. This diagonalizes the spectrum covariance matrix, and the information content of the corresponding power spectrum is increased by a factor of two over standard methods. Unfortunately, careful calculation of the Cramer-Rao bound shows that the information recovery can never be made complete, any observable built from shear fields, including optimal sufficient statistics, are subject to severe information loss, typically $80\%$ to $90\%$ below $\ell \sim 3000$ for generic cosmological parameters. The lost information can only be recovered from additional, non-shear based data. Our predictions hold just as well for a tomographic analysis, and/or full sky surveys. △ Less

Submitted 25 November, 2014; v1 submitted 7 August, 2014; originally announced August 2014.

Comments: 11 pages, 4 figures, matches version accepted for publication in MNRAS

Supervoid Origin of the Cold Spot in the Cosmic Microwave Background

Authors: András Kovács, István Szapudi, Benjamin R. Granett, Zsolt Frei, Joseph Silk, Will Burgett, Shaun Cole, Peter W. Draper, Daniel J. Farrow, Nicholas Kaiser, Eugene A. Magnier, Nigel Metcalfe, Jeffrey S. Morgan, Paul Price, John Tonry, Richard Wainscoat

Abstract: We use a WISE-2MASS-Pan-STARRS1 galaxy catalog to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. We obtain photometric redshifts using our multicolor data set to create a tomographic map of the galaxy distribution. The radial density profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic… ▽ More We use a WISE-2MASS-Pan-STARRS1 galaxy catalog to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. We obtain photometric redshifts using our multicolor data set to create a tomographic map of the galaxy distribution. The radial density profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. Our data, combined with an earlier measurement by Granett et al 2010, are consistent with a large $R_{\rm void}=(192 \pm 15)h^{-1} Mpc $ $(2σ)$ supervoid with $δ\simeq -0.13 \pm 0.03$ centered at $z=0.22\pm0.01$. Such a supervoid, constituting a $\sim3.5 σ$ fluctuation in the $ΛCDM$ model, is a plausible cause for the Cold Spot. △ Less

Submitted 6 July, 2014; originally announced July 2014.

Comments: 4 pages, 2 figures, Proceedings of IAU 306 Symposium: Statistical Challenges in 21st Century Cosmology

Optimal observables in galaxy surveys

Authors: Julien Carron, István Szapudi

Abstract: The sufficient statistics of the one-point probability density function of the dark matter density field is worked out using cosmological perturbation theory and tested to the Millennium simulation density field. The logarithmic transformation is recovered for spectral index close to $-1$ as a special case of the family of power transformations. We then discuss how these transforms should be modif… ▽ More The sufficient statistics of the one-point probability density function of the dark matter density field is worked out using cosmological perturbation theory and tested to the Millennium simulation density field. The logarithmic transformation is recovered for spectral index close to $-1$ as a special case of the family of power transformations. We then discuss how these transforms should be modified in the case of noisy tracers of the field and focus on the case of Poisson sampling. This gives us optimal local transformations to apply to galaxy survey data prior the extraction of the spectrum in order to capture most efficiently the information encoded in large scale structures. △ Less

Submitted 1 July, 2014; originally announced July 2014.

Comments: Proceedings of IAU 306 Symposium Statistical Challenges in 21st Century Cosmology

On fast Generation of Cosmological Random Fields

Authors: Julien Carron, Melody Wolk, Istvan Szapudi

Abstract: The statistical translation invariance of cosmological random fields is broken by a finite survey boundary, correlating the observable Fourier modes. Standard methods for generating Gaussian fields either neglect these correlations, or are costly, or both. Here we report on a fast and exact simulation method applicable to a wide class of two-point statistics that requires the simulation of a perio… ▽ More The statistical translation invariance of cosmological random fields is broken by a finite survey boundary, correlating the observable Fourier modes. Standard methods for generating Gaussian fields either neglect these correlations, or are costly, or both. Here we report on a fast and exact simulation method applicable to a wide class of two-point statistics that requires the simulation of a periodic grid of only twice the survey side with fast Fourier transforms. Super-survey modes, dominating the covariance of power spectra beyond linear scales in galaxy surveys and causing the correlation of large and small scales, "beat coupling", or "super-sample" covariance, are precisely accounted for in non-linear transformations of the Gaussian field. As an application, we simulate the CFHTLS $\sim 7^\circ \times 7^\circ$ W1 galaxy density field, modeled as a Poisson sampling of a lognormal density field. We show that our simulations produce power spectra, $A^*$-power spectra, counts-in-cells probability distributions as well as covariances perfectly consistent with the data. In addition, our technique reproduces the information plateau beyond linear scales as observed previously in Sloan Digital Sky Survey galaxy catalogs and in $N$-body simulations. Our method is thus an efficient yet powerful simulation and prediction tool for galaxy survey data and covariances. △ Less

Submitted 28 July, 2014; v1 submitted 23 June, 2014; originally announced June 2014.

Comments: 8 pages, 5 figures, matches version accepted for publication in MNRAS

The Cold Spot in the Cosmic Microwave Background: the Shadow of a Supervoid

Authors: István Szapudi, András Kovács, Benjamin R. Granett, Zsolt Frei, Joseph Silk, Juan Garcia-Bellido, Will Burgett, Shaun Cole, Peter W. Draper, Daniel J. Farrow, Nicholas Kaiser, Eugene A. Magnier, Nigel Metcalfe, Jeffrey S. Morgan, Paul Price, John Tonry, Richard Wainscoat

Abstract: Standard inflationary hot big bang cosmology predicts small fluctuations in the Cosmic Microwave Background (CMB) with isotropic Gaussian statistics. All measurements support the standard theory, except for a few anomalies discovered in the Wilkinson Microwave Anisotropy Probe maps and confirmed recently by the Planck satellite. The Cold Spot is one of the most significant of such anomalies, and t… ▽ More Standard inflationary hot big bang cosmology predicts small fluctuations in the Cosmic Microwave Background (CMB) with isotropic Gaussian statistics. All measurements support the standard theory, except for a few anomalies discovered in the Wilkinson Microwave Anisotropy Probe maps and confirmed recently by the Planck satellite. The Cold Spot is one of the most significant of such anomalies, and the leading explanation of it posits a large void that imprints this extremely cold area via the linear Integrated Sachs-Wolfe (ISW) effect due to the decay of gravitational potentials over cosmic time, or via the Rees-Sciama (RS) effect due to late-time non-linear evolution. Despite several observational campaigns targeting the Cold Spot region, to date no suitably large void was found at higher redshifts $z > 0.3$. Here we report the detection of an $R =(192 \pm 15) h^{-1}Mpc$ size supervoid of depth $δ= -0.13 \pm 0.03$, and centred at redshift $z = 0.22$. This supervoid, possibly the largest ever found, is large enough to significantly affect the CMB via the non-linear RS effect, as shown in our Lemaitre-Tolman-Bondi framework. This discovery presents the first plausible explanation for any of the physical CMB anomalies, and raises the possibility that local large-scale structure could be responsible for other anomalies as well. △ Less

Submitted 15 June, 2014; originally announced June 2014.

Comments: 8 pages, 4 figures, to appear in the proceedings of the Moriond Cosmology Conference 2014