Showing 1–50 of 206 results for author: Alonso, D

The Simons Observatory: component separation pipelines for B-modes

Authors: Kevin Wolz, Susanna Azzoni, Carlos Hervías-Caimapo, Josquin Errard, Nicoletta Krachmalnicoff, David Alonso, Benjamin Beringue, Emilie Hertig

Abstract: The upcoming Simons Observatory (SO) Small Aperture Telescopes aim at observing the degree-scale anisotropies of the polarized CMB to constrain the primordial tensor-to-scalar ratio $r$ at the level of $σ(r=0)\lesssim0.003$ to probe models of the very early Universe. We present three complementary $r$ inference pipelines and compare their results on a set of sky simulations that allow us to explor… ▽ More The upcoming Simons Observatory (SO) Small Aperture Telescopes aim at observing the degree-scale anisotropies of the polarized CMB to constrain the primordial tensor-to-scalar ratio $r$ at the level of $σ(r=0)\lesssim0.003$ to probe models of the very early Universe. We present three complementary $r$ inference pipelines and compare their results on a set of sky simulations that allow us to explore a number of Galactic foreground and instrumental noise models, relevant for SO. In most scenarios, the pipelines retrieve consistent and unbiased results. However, several complex foreground scenarios lead to a $>2σ$ bias on $r$ if analyzed with the default versions of these pipelines, highlighting the need for more sophisticated pipeline components that marginalize over foreground residuals. We present two such extensions, using power-spectrum-based and map-based methods, and show that they fully reduce the bias on $r$ to sub-sigma level in all scenarios, and at a moderate cost in terms of $σ(r)$. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 4 pages, 1 figure. Contribution to the 2024 Cosmology session of the 58th Rencontres de Moriond

Tomographic constraints on the production rate of gravitational waves from astrophysical sources

Authors: David Alonso, Mehraveh Nikjoo, Arianna I. Renzini, Emilio Bellini, Pedro G. Ferreira

Abstract: Using an optimal quadratic estimator, we measure the large-scale cross-correlation between maps of the stochastic gravitational-wave intensity, constructed from the first three LIGO-Virgo observing runs, and a suite of tomographic samples of galaxies covering the redshift range $z\lesssim 2$. We do not detect any statistically significant cross-correlation, but the tomographic nature of the data a… ▽ More Using an optimal quadratic estimator, we measure the large-scale cross-correlation between maps of the stochastic gravitational-wave intensity, constructed from the first three LIGO-Virgo observing runs, and a suite of tomographic samples of galaxies covering the redshift range $z\lesssim 2$. We do not detect any statistically significant cross-correlation, but the tomographic nature of the data allows us to place constraints on the (bias-weighted) production rate density of gravitational waves by astrophysical sources as a function of cosmic time. Our constraints range from $\langle b\dotΩ_{\rm GW}\rangle<3.0\times10^{-9}\,{\rm Gyr}^{-1}$ at $z\sim0.06$ to $\langle b\dotΩ_{\rm GW}\rangle<2.7\times10^{-7}\,{\rm Gyr}^{-1}$ at $z\sim1.5$ (95\% C.L.), assuming a frequency spectrum of the form $f^{2/3}$ (corresponding to an astrophysical background of binary mergers), and a reference frequency $f_{\rm ref}=25\,{\rm Hz}$. Although these constraints are $\sim2$ orders of magnitude higher than the expected signal, we show that a detection may be possible with future experiments. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 14 pages, 11 figures

Assessment of Gradient-Based Samplers in Standard Cosmological Likelihoods

Authors: Arrykrishna Mootoovaloo, Jaime Ruiz-Zapatero, Carlos García-García, David Alonso

Abstract: We assess the usefulness of gradient-based samplers, such as the No-U-Turn Sampler (NUTS), by comparison with traditional Metropolis-Hastings algorithms, in tomographic $3 \times 2$ point analyses. Specifically, we use the DES Year 1 data and a simulated future LSST-like survey as representative examples of these studies, containing a significant number of nuisance parameters (20 and 32, respectiv… ▽ More We assess the usefulness of gradient-based samplers, such as the No-U-Turn Sampler (NUTS), by comparison with traditional Metropolis-Hastings algorithms, in tomographic $3 \times 2$ point analyses. Specifically, we use the DES Year 1 data and a simulated future LSST-like survey as representative examples of these studies, containing a significant number of nuisance parameters (20 and 32, respectively) that affect the performance of rejection-based samplers. To do so, we implement a differentiable forward model using JAX-COSMO (Campagne et al. 2023), and we use it to derive parameter constraints from both datasets using the NUTS algorithm as implemented in §4, and the Metropolis-Hastings algorithm as implemented in Cobaya (Lewis 2013). When quantified in terms of the number of effective number of samples taken per likelihood evaluation, we find a relative efficiency gain of $\mathcal{O}(10)$ in favour of NUTS. However, this efficiency is reduced to a factor $\sim 2$ when quantified in terms of computational time, since we find the cost of the gradient computation (needed by NUTS) relative to the likelihood to be $\sim 4.5$ times larger for both experiments. We validate these results making use of analytical multi-variate distributions (a multivariate Gaussian and a Rosenbrock distribution) with increasing dimensionality. Based on these results, we conclude that gradient-based samplers such as NUTS can be leveraged to sample high dimensional parameter spaces in Cosmology, although the efficiency improvement is relatively mild for moderate $(\mathcal{O}(50))$ dimension numbers, typical of tomographic large-scale structure analyses. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 14 pages and 9 figures

LLM+Reasoning+Planning for supporting incomplete user queries in presence of APIs

Authors: Sudhir Agarwal, Anu Sreepathy, David H. Alonso, Prarit Lamba

Abstract: Recent availability of Large Language Models (LLMs) has led to the development of numerous LLM-based approaches aimed at providing natural language interfaces for various end-user tasks. These end-user tasks in turn can typically be accomplished by orchestrating a given set of APIs. In practice, natural language task requests (user queries) are often incomplete, i.e., they may not contain all the… ▽ More Recent availability of Large Language Models (LLMs) has led to the development of numerous LLM-based approaches aimed at providing natural language interfaces for various end-user tasks. These end-user tasks in turn can typically be accomplished by orchestrating a given set of APIs. In practice, natural language task requests (user queries) are often incomplete, i.e., they may not contain all the information required by the APIs. While LLMs excel at natural language processing (NLP) tasks, they frequently hallucinate on missing information or struggle with orchestrating the APIs. The key idea behind our proposed approach is to leverage logical reasoning and classical AI planning along with an LLM for accurately answering user queries including identification and gathering of any missing information in these queries. Our approach uses an LLM and ASP (Answer Set Programming) solver to translate a user query to a representation in Planning Domain Definition Language (PDDL) via an intermediate representation in ASP. We introduce a special API "get_info_api" for gathering missing information. We model all the APIs as PDDL actions in a way that supports dataflow between the APIs. Our approach then uses a classical AI planner to generate an orchestration of API calls (including calls to get_info_api) to answer the user query. Our evaluation results show that our approach significantly outperforms a pure LLM based approach by achieving over 95\% success rate in most cases on a dataset containing complete and incomplete single goal and multi-goal queries where the multi-goal queries may or may not require dataflow among the APIs. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 9 pages main content, 2 pages references, 12 pages appendix, 5 figures

The Simons Observatory: Combining cross-spectral foreground cleaning with multi-tracer $B$-mode delensing for improved constraints on inflation

Authors: Emilie Hertig, Kevin Wolz, Toshiya Namikawa, Antón Baleato Lizancos, Susanna Azzoni, Irene Abril-Cabezas, David Alonso, Carlo Baccigalupi, Erminia Calabrese, Anthony Challinor, Josquin Errard, Giulio Fabbian, Baptiste Jost, Nicoletta Krachmalnicoff, Anto I. Lonappan, Magdy Morshed, Luca Pagano, Blake Sherwin

Abstract: The Simons Observatory (SO), due to start full science operations in early 2025, aims to set tight constraints on inflationary physics by inferring the tensor-to-scalar ratio $r$ from measurements of CMB polarization $B$-modes. Its nominal design targets a precision $σ(r=0) \leq 0.003$ without delensing. Achieving this goal and further reducing uncertainties requires the mitigation of other source… ▽ More The Simons Observatory (SO), due to start full science operations in early 2025, aims to set tight constraints on inflationary physics by inferring the tensor-to-scalar ratio $r$ from measurements of CMB polarization $B$-modes. Its nominal design targets a precision $σ(r=0) \leq 0.003$ without delensing. Achieving this goal and further reducing uncertainties requires the mitigation of other sources of large-scale $B$-modes such as Galactic foregrounds and weak gravitational lensing. We present an analysis pipeline aiming to estimate $r$ by including delensing within a cross-spectral likelihood, and demonstrate it on SO-like simulations. Lensing $B$-modes are synthesised using internal CMB lensing reconstructions as well as Planck-like CIB maps and LSST-like galaxy density maps. This $B$-mode template is then introduced into SO's power-spectrum-based foreground-cleaning algorithm by extending the likelihood function to include all auto- and cross-spectra between the lensing template and the SAT $B$-modes. Within this framework, we demonstrate the equivalence of map-based and cross-spectral delensing and use it to motivate an optimized pixel-weighting scheme for power spectrum estimation. We start by validating our pipeline in the simplistic case of uniform foreground spectral energy distributions (SEDs). In the absence of primordial $B$-modes, $σ(r)$ decreases by 37% as a result of delensing. Tensor modes at the level of $r=0.01$ are successfully detected by our pipeline. Even with more realistic foreground models including spatial variations in the dust and synchrotron spectral properties, we obtain unbiased estimates of $r$ by employing the moment-expansion method. In this case, delensing-related improvements range between 27% and 31%. These results constitute the first realistic assessment of the delensing performance at SO's nominal sensitivity level. (Abridged) △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 24 pages, 9 figures, 3 tables. Submitted to Physical Review D

Relativistic imprints on dispersion measure space distortions

Authors: Shohei Saga, David Alonso

Abstract: We investigate the three-dimensional clustering of sources emitting electromagnetic pulses traveling through cold electron plasma, whose radial distance is inferred from their dispersion measure. As a distance indicator, dispersion measure is systematically affected by inhomogeneities in the electron density along the line of sight and special and general relativistic effects, similar to the case… ▽ More We investigate the three-dimensional clustering of sources emitting electromagnetic pulses traveling through cold electron plasma, whose radial distance is inferred from their dispersion measure. As a distance indicator, dispersion measure is systematically affected by inhomogeneities in the electron density along the line of sight and special and general relativistic effects, similar to the case of redshift surveys. We present analytic expressions for the correlation function of fast radio bursts (FRBs), and for the galaxy-FRB cross-correlation function, in the presence of these dispersion measure-space distortions. We find that the even multipoles of these correlations are primarily dominated by non-local contributions (e.g. the electron density fluctuations integrated along the line of sight), while the dipole also receives a significant contribution from the Doppler effect, one of the major relativistic effects. A large number of FRBs, $\mathcal{O}(10^{5}\sim10^{6})$, expected to be observed in the Square Kilometre Array, would be enough to measure the even multipoles at very high significance, ${\rm S}/{\rm N} \approx 100$, and perhaps to make a first detection of the dipole (${\rm S}/{\rm N} \approx 10$) in the FRB correlation function and FRB-galaxy cross correlation function. This measurement could open a new window to study and test cosmological models. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 16 pages, 10 figures

Cosmic shear with small scales: DES-Y3, KiDS-1000 and HSC-DR1

Authors: Carlos García-García, Matteo Zennaro, Giovanni Aricò, David Alonso, Raul E. Angulo

Abstract: We present a cosmological analysis of the combination of the DES-Y3, KiDS-1000 and HSC-DR1 weak lensing samples under a joint harmonic-space pipeline making use of angular scales down to $\ell_{\rm max}=4500$, corresponding to significantly smaller scales ($δθ\sim2.4'$) than those commonly used in cosmological weak lensing studies. We are able to do so by accurately modelling non-linearities and t… ▽ More We present a cosmological analysis of the combination of the DES-Y3, KiDS-1000 and HSC-DR1 weak lensing samples under a joint harmonic-space pipeline making use of angular scales down to $\ell_{\rm max}=4500$, corresponding to significantly smaller scales ($δθ\sim2.4'$) than those commonly used in cosmological weak lensing studies. We are able to do so by accurately modelling non-linearities and the impact of baryonic effects using Baccoemu. We find $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3}=0.795^{+0.015}_{-0.017}$, in relatively good agreement with CMB constraints from Planck (less than $\sim1.8σ$ tension), although we obtain a low value of $Ω_{\rm m}=0.212^{+0.017}_{-0.032}$, in tension with Planck at the $\sim3σ$ level. We show that this can be recast as an $H_0$ tension if one parametrises the amplitude of fluctuations and matter abundance in terms of variables without hidden dependence on $H_0$. Furthermore, we find that this tension reduces significantly after including a prior on the distance-redshift relationship from BAO data, without worsening the fit. In terms of baryonic effects, we show that failing to model and marginalise over them on scales $\ell\lesssim2000$ does not significantly affect the posterior constraints for DES-Y3 and KiDS-1000, but has a mild effect on deeper samples, such as HSC-DR1. This is in agreement with our ability to only mildly constrain the parameters of the Baryon Correction Model with these data △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 57 pages, 24 figures, comments welcome

Accurate heat currents via reorganised master equation

Authors: Jonas Glatthard, Guillem Aznar-Menargues, José P. Palao, Daniel Alonso, Luis A. Correa

Abstract: The accurate characterisation of energy exchanges between nanoscale quantum systems and their environments is of paramount importance for quantum technologies, and central to quantum thermodynamics. Here, we show that, in order to accurately approximate steady-state heat currents via perturbative master equations, the coupling-induced reorganisation correction to the system's energy must be carefu… ▽ More The accurate characterisation of energy exchanges between nanoscale quantum systems and their environments is of paramount importance for quantum technologies, and central to quantum thermodynamics. Here, we show that, in order to accurately approximate steady-state heat currents via perturbative master equations, the coupling-induced reorganisation correction to the system's energy must be carefully taken into account. Not doing so, may yield sizeable errors, especially at low, or even moderate temperatures. In particular, we show how a 'reorganised master equation' can produce very accurate estimates for the heat currents when the reorganisation energy is weak and one works with environments with a broad spectrum. Notably, such master equation outperforms its 'non-reorganised' counterpart in the calculation of heat currents, at modelling dynamics, and at correctly capturing equilibration. This is so even if both types of equation are derived to the same order of perturbation theory. Most importantly, working with reorganised master equations does not involve additional complications when compared with alternative approaches. Also, invoking the secular approximation to secure thermodynamic consistency does not compromise their precision. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures; comments welcome

The Simons Observatory: impact of bandpass, polarization angle and calibration uncertainties on small-scale power spectrum analysis

Authors: S. Giardiello, M. Gerbino, L. Pagano, D. Alonso, B. Beringue, B. Bolliet, E. Calabrese, G. Coppi, J. Errard, G. Fabbian, I. Harrison, J. C. Hill, H. T. Jense, B. Keating, A. La Posta, M. Lattanzi, A. I. Lonappan, G. Puglisi, C. L. Reichardt, S. M. Simon

Abstract: We study the effects due to mismatches in passbands, polarization angles, and temperature and polarization calibrations in the context of the upcoming cosmic microwave background experiment Simons Observatory (SO). Using the SO multi-frequency likelihood, we estimate the bias and the degradation of constraining power in cosmological and astrophysical foreground parameters assuming different levels… ▽ More We study the effects due to mismatches in passbands, polarization angles, and temperature and polarization calibrations in the context of the upcoming cosmic microwave background experiment Simons Observatory (SO). Using the SO multi-frequency likelihood, we estimate the bias and the degradation of constraining power in cosmological and astrophysical foreground parameters assuming different levels of knowledge of the instrumental effects. We find that incorrect but reasonable assumptions on the values of all the systematics examined here can have important effects in cosmological analyses, hence requiring marginalization approaches at likelihood level. When doing so, we find that the most relevant effect is due to bandpass shifts. When marginalizing over them, the posteriors of parameters describing astrophysical microwave foregrounds (such as radio point sources or dust) get degraded, while cosmological parameters constraints are not significantly affected. Marginalization over polarization angles with up to 0.25$^\circ$ uncertainty causes an irrelevant bias $\lesssim 0.05 σ$ in all parameters. Marginalization over calibration factors in polarization broadens the constraints on the effective number of relativistic degrees of freedom $N_\mathrm{eff}$ by a factor 1.2, interpreted here as a proxy parameter for non standard model physics targeted by high-resolution CMB measurements. △ Less

Submitted 12 March, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 61 pages (27 of which are appendices), 26 figures, 12 tables. (Replacement to update the acknowledgments)

Computing Transiting Exoplanet Parameters with 1D Convolutional Neural Networks

Authors: Santiago Iglesias Álvarez, Enrique Díez Alonso, María Luisa Sánchez Rodríguez, Javier Rodríguez Rodríguez, Saúl Pérez Fernández, Francisco Javier de Cos Juez

Abstract: The transit method allows the detection and characterization of planetary systems by analyzing stellar light curves. Convolutional neural networks appear to offer a viable solution for automating these analyses. In this research, two 1D convolutional neural network models, which work with simulated light curves in which transit-like signals were injected, are presented. One model operates on compl… ▽ More The transit method allows the detection and characterization of planetary systems by analyzing stellar light curves. Convolutional neural networks appear to offer a viable solution for automating these analyses. In this research, two 1D convolutional neural network models, which work with simulated light curves in which transit-like signals were injected, are presented. One model operates on complete light curves and estimates the orbital period, and the other one operates on phase-folded light curves and estimates the semimajor axis of the orbit and the square of the planet-to-star radius ratio. Both models were tested on real data from TESS light curves with confirmed planets to ensure that they are able to work with real data. The results obtained show that 1D CNNs are able to characterize transiting exoplanets from their host star's detrended light curve and, furthermore, reducing both the required time and computational costs compared with the current detection and characterization algorithms. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Growth history and quasar bias evolution at z < 3 from Quaia

Authors: G. Piccirilli, G. Fabbian, D. Alonso, K. Storey-Fisher, J. Carron, A. Lewis, C. García-García

Abstract: We make use of the Gaia-Unwise quasar catalogue, Quaia, to constrain the growth history out to high redshifts from the clustering of quasars and their cross-correlation with maps of the Cosmic Microwave Background (CMB) lensing convergence. Considering three tomographic bins, centered at redshifts $\bar{z}_i = [0.69, 1.59, 2.72]$, we reconstruct the evolution of the amplitude of matter fluctuation… ▽ More We make use of the Gaia-Unwise quasar catalogue, Quaia, to constrain the growth history out to high redshifts from the clustering of quasars and their cross-correlation with maps of the Cosmic Microwave Background (CMB) lensing convergence. Considering three tomographic bins, centered at redshifts $\bar{z}_i = [0.69, 1.59, 2.72]$, we reconstruct the evolution of the amplitude of matter fluctuations $σ_8(z)$ over the last $\sim12$ billion years of cosmic history. In particular, we make one of the highest-redshift measurements of $σ_8$ ($σ_8(z=2.72)=0.22\pm 0.06$), finding it to be in good agreement (at the $\sim1σ$ level) with the value predicted by $Λ$CDM using CMB data from Planck. We also used the data to study the evolution of the linear quasar bias for this sample, finding values similar to those of other quasar samples, although with a less steep evolution at high redshifts. Finally, we study the potential impact of foreground contamination in the CMB lensing maps and, although we find evidence of contamination in cross-correlations at $z\sim1.7$ we are not able to clearly pinpoint its origin as being Galactic or extragalactic. Nevertheless, we determine that the impact of this contamination on our results is negligible. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 28 pages, 11 figures

Out-of-equlibrium inference of stochastic model parameters through population data from generic consumer-resource dynamics

Authors: Jose A. Capitan, David Alonso

Abstract: Consumer-resource dynamics is central in determining biomass transport across ecosystems. The assumptions of mass action, chemostatic conditions and stationarity in stochastic feeding dynamics lead to Holling type II functional responses, whose use is widespread in macroscopic models of population dynamics. However, to be useful for parameter inference, stochastic population models need to be iden… ▽ More Consumer-resource dynamics is central in determining biomass transport across ecosystems. The assumptions of mass action, chemostatic conditions and stationarity in stochastic feeding dynamics lead to Holling type II functional responses, whose use is widespread in macroscopic models of population dynamics. However, to be useful for parameter inference, stochastic population models need to be identifiable, this meaning that model parameters can be uniquely inferred from a large number of model observations. In this contribution we study parameter identifiability in a multi-resource consumer-resource model, for which we can obtain the steady-state and out-of-equilibrium probability distributions of predator's abundances by analytically solving the master equation. Based on these analytical results, we can conduct in silico experiments by tracking the abundance of consumers that are either searching for or handling prey, data then used for maximum likelihood parameter estimation. We show that, when model observations are recorded out of equilibrium, feeding parameters are truly identifiable, whereas if sampling is done at stationarity, only ratios of rates can be inferred from data. We discuss the implications of our results when inferring parameters of general dynamical models. △ Less

Submitted 3 January, 2024; originally announced January 2024.

MSC Class: 92-10 60Gxx 62Fxx (Primary); 05A15 (Secondary)

Improved models for near-Earth asteroids (2100) Ra-Shalom, (3103) Eger, (12711) Tukmit & (161989) Cacus

Authors: Javier Rodríguez Rodríguez, Enrique Díez Alonso, Santiago Iglesias Álvarez, Saúl Pérez Fernández, Javier Licandro, Miguel R. Alarcon, Miquel Serra-Ricart, Noemi Pinilla-Alonso, Susana Fernández Menéndez, Francisco Javier de Cos Juez

Abstract: We present 24 new dense lightcurves of the near-Earth asteroids (3103) Eger, (161989) Cacus, (2100) Ra-Shalom and (12711) Tukmit, obtained with the Instituto Astrofísico Canarias 80 and Telescopio Abierto Remoto 2 telescopes at the Teide Observatory (Tenerife, Spain) during 2021 and 2022, in the framework of projects visible NEAs observations survey and NEO Rapid Observation, Characterization and… ▽ More We present 24 new dense lightcurves of the near-Earth asteroids (3103) Eger, (161989) Cacus, (2100) Ra-Shalom and (12711) Tukmit, obtained with the Instituto Astrofísico Canarias 80 and Telescopio Abierto Remoto 2 telescopes at the Teide Observatory (Tenerife, Spain) during 2021 and 2022, in the framework of projects visible NEAs observations survey and NEO Rapid Observation, Characterization and Key Simulations. The shape models and rotation state parameters ($P$, $λ$, $β$) were computed by applying the lightcurve inversion method to the new data altogether with the archival data. For (3013) Eger and (161989) Cacus, our shape models and rotation state parameters agree with previous works, though they have smaller uncertainties. For (2100) Ra-Shalom, our results also agree with previous studies. Still, we find that a Yarkovsky - O'Keefe - Radzievskii - Paddack acceleration of $\upsilon = (0.223\pm0.237)\times10^{-8}$ rad d$^{-2}$ slightly improves the fit of the lightcurves, suggesting that (2100) Ra-Shalom could be affected by this acceleration. We also present for the first time a shape model for (12711) Tukmit, along with its rotation state parameters ($P=3.484900 \pm 0.000031$ hr, $λ= 27^{\circ}\pm 8^{\circ}$, $β= 9^{\circ} \pm 15^{\circ}$). △ Less

Submitted 15 December, 2023; originally announced December 2023.

One-dimensional Convolutional Neural Networks for Detecting Transiting Exoplanets

Authors: Santiago Iglesias Álvarez, Enrique Díez Alonso, María Luisa Sánchez, Javier Rodríguez Rodríguez, Fernando Sánchez Lasheras, Francisco Javier de Cos Juez

Abstract: The transit method is one of the most relevant exoplanet detection techniques, which consists of detecting periodic eclipses in the light curves of stars. This is not always easy due to the presence of noise in the light curves, which is induced, for example, by the response of a telescope to stellar flux. For this reason, we aimed to develop an artificial neural network model that is able to dete… ▽ More The transit method is one of the most relevant exoplanet detection techniques, which consists of detecting periodic eclipses in the light curves of stars. This is not always easy due to the presence of noise in the light curves, which is induced, for example, by the response of a telescope to stellar flux. For this reason, we aimed to develop an artificial neural network model that is able to detect these transits in light curves obtained from different telescopes and surveys. We created artificial light curves with and without transits to try to mimic those expected for the extended mission of the Kepler telescope (K2) in order to train and validate a 1D convolutional neural network model, which was later tested, obtaining an accuracy of 99.02 % and an estimated error (loss function) of 0.03. These results, among others, helped to confirm that the 1D CNN is a good choice for working with non-phased-folded Mandel and Agol light curves with transits. It also reduces the number of light curves that have to be visually inspected to decide if they present transit-like signals and decreases the time needed for analyzing each (with respect to traditional analysis). △ Less

Submitted 12 December, 2023; originally announced December 2023.

Accuracy requirements on intrinsic alignments for Stage-IV cosmic shear

Authors: Anya Paopiamsap, Natalia Porqueres, David Alonso, Joachim Harnois-Deraps, C. Danielle Leonard

Abstract: In the context of cosmological weak lensing studies, intrinsic alignments (IAs) are one the most complicated astrophysical systematic to model, given the poor understanding of the physical processes that cause them. A number of modelling frameworks for IAs have been proposed in the literature, both purely phenomenological or grounded on a perturbative treatment of symmetry-based arguments. However… ▽ More In the context of cosmological weak lensing studies, intrinsic alignments (IAs) are one the most complicated astrophysical systematic to model, given the poor understanding of the physical processes that cause them. A number of modelling frameworks for IAs have been proposed in the literature, both purely phenomenological or grounded on a perturbative treatment of symmetry-based arguments. However, the accuracy with which any of these approaches is able to describe the impact of IAs on cosmic shear data, particularly on the comparatively small scales ($k\simeq 1\,{\rm Mpc}^{-1}$) to which this observable is sensitive, is not clear. Here we quantify the level of disagreement between the true underlying intrinsic alignments and the theoretical model used to describe them that can be allowed in the context of cosmic shear analyses with future Stage-IV surveys. We consider various models describing this "IA residual", covering both physics-based approaches, as well as completely agnostic prescriptions. The same qualitative results are recovered in all cases explored: for a Stage-IV cosmic shear survey, a mis-modelling of the IA contribution at the $\sim10\%$ level produces shifts of $\lesssim0.5σ$ on the final cosmological parameter constraints. Current and future IA models should therefore aim to achieve this level of accuracy, a prospect that is not unfeasible for models with sufficient flexibility. △ Less

Submitted 8 May, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

A precise symbolic emulator of the linear matter power spectrum

Authors: Deaglan J. Bartlett, Lukas Kammerer, Gabriel Kronberger, Harry Desmond, Pedro G. Ferreira, Benjamin D. Wandelt, Bogdan Burlacu, David Alonso, Matteo Zennaro

Abstract: Computing the matter power spectrum, $P(k)$, as a function of cosmological parameters can be prohibitively slow in cosmological analyses, hence emulating this calculation is desirable. Previous analytic approximations are insufficiently accurate for modern applications, so black-box, uninterpretable emulators are often used. We utilise an efficient genetic programming based symbolic regression fra… ▽ More Computing the matter power spectrum, $P(k)$, as a function of cosmological parameters can be prohibitively slow in cosmological analyses, hence emulating this calculation is desirable. Previous analytic approximations are insufficiently accurate for modern applications, so black-box, uninterpretable emulators are often used. We utilise an efficient genetic programming based symbolic regression framework to explore the space of potential mathematical expressions which can approximate the power spectrum and $σ_8$. We learn the ratio between an existing low-accuracy fitting function for $P(k)$ and that obtained by solving the Boltzmann equations and thus still incorporate the physics which motivated this earlier approximation. We obtain an analytic approximation to the linear power spectrum with a root mean squared fractional error of 0.2% between $k = 9\times10^{-3} - 9 \, h{\rm \, Mpc^{-1}}$ and across a wide range of cosmological parameters, and we provide physical interpretations for various terms in the expression. Our analytic approximation is 950 times faster to evaluate than camb and 36 times faster than the neural network based matter power spectrum emulator BACCO. We also provide a simple analytic approximation for $σ_8$ with a similar accuracy, with a root mean squared fractional error of just 0.1% when evaluated across the same range of cosmologies. This function is easily invertible to obtain $A_{\rm s}$ as a function of $σ_8$ and the other cosmological parameters, if preferred. It is possible to obtain symbolic approximations to a seemingly complex function at a precision required for current and future cosmological analyses without resorting to deep-learning techniques, thus avoiding their black-box nature and large number of parameters. Our emulator will be usable long after the codes on which numerical approximations are built become outdated. △ Less

Submitted 15 April, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 9 pages, 5 figures. Accepted for publication in A&A

Journal ref: A&A 686, A209 (2024)

Characterization in Geant4 of different PET configurations

Authors: M. L. López Toxqui, C. H. Zepeda Fernández, L. F. Rebolledo Herrera, B. De Celis Alonso, E. Moreno Barbosa

Abstract: Positron Emission Tomography (PET) is a Nuclear Medicine technique that creates images that allow the study of metabolic activity and organ function using radiopharmaceuticals. Continuous improvement of scintillation detectors for radiation in PET as well as improvement in electronic detectors (e.g., SiPMs) and signal processing, makes the field of PET a fast and changing environment. If industry… ▽ More Positron Emission Tomography (PET) is a Nuclear Medicine technique that creates images that allow the study of metabolic activity and organ function using radiopharmaceuticals. Continuous improvement of scintillation detectors for radiation in PET as well as improvement in electronic detectors (e.g., SiPMs) and signal processing, makes the field of PET a fast and changing environment. If industry desires to build new systems implementing these technological improvements, it is of its interest to develop modelling strategies that can provide information on how to build them, reducing time and material costs. Bearing this in mind three different PET configurations, were simulated in Geant4, to determine which one presented the best performance according to quality parameters such as spatial resolution (SR), coincident time resolution (CTR) and acceptance value (A). This was done with three different (in size) pairs of LYSO crystals + SiPM detectors. It was found that the 2 Modules system presented worst results than the two Ring detector configurations. Between the Ring configurations the first was marginally better than the second. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Modelling cross-correlations of ultra-high-energy cosmic rays and galaxies

Authors: Federico R. Urban, Stefano Camera, David Alonso

Abstract: The astrophysical engines that power ultra-high-energy cosmic rays (UHECRs) remain to date unknown. Since the propagation horizon of UHECRs is limited to the local, anisotropic Universe, the distribution of UHECR arrival directions should be anisotropic. In this paper we expand the analysis of the potential for the angular, harmonic cross-correlation between UHECRs and galaxies to detect such anis… ▽ More The astrophysical engines that power ultra-high-energy cosmic rays (UHECRs) remain to date unknown. Since the propagation horizon of UHECRs is limited to the local, anisotropic Universe, the distribution of UHECR arrival directions should be anisotropic. In this paper we expand the analysis of the potential for the angular, harmonic cross-correlation between UHECRs and galaxies to detect such anisotropies. We do so by studying simulations performed assuming proton, oxygen and silicon injection models, each simulation containing a number of events comparable to a conservative estimate of currently available datasets, as well as by extending the analytic treatment of the magnetic deflections. Quantitatively, we find that, while the correlations for each given multipole are generally weak, (1) the total harmonic power summed over multipoles is detectable with signal-to-noise ratios well above~5 for both the auto-correlation and the cross-correlation (once optimal weights are applied) in most cases studied here, with peaks of signal-to-noise ratio around between~8 and~10 at the highest energies; (2) if we combine the UHECR auto-correlation and the cross-correlation we are able to reach detection levels of (3σ) and above for individual multipoles at the largest scales, especially for heavy composition. In particular, we predict that the combined-analysis quadrupole could be detected already with existing data. △ Less

Submitted 21 February, 2024; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: Accepted/published version

Journal ref: The Open Journal of Astrophysics 7 (February), 2024

LimberJack.jl: auto-differentiable methods for angular power spectra analyses

Authors: J. Ruiz-Zapatero, D. Alonso, C. García-García, A. Nicola, A. Mootoovaloo, J. M. Sullivan, M. Bonici, P. G. Ferreira

Abstract: We present LimberJack.jl, a fully auto-differentiable code for cosmological analyses of 2 point auto- and cross-correlation measurements from galaxy clustering, CMB lensing and weak lensing data written in Julia. Using Julia's auto-differentiation ecosystem, LimberJack.jl can obtain gradients for its outputs up to an order of magnitude faster than traditional finite difference methods. This makes… ▽ More We present LimberJack.jl, a fully auto-differentiable code for cosmological analyses of 2 point auto- and cross-correlation measurements from galaxy clustering, CMB lensing and weak lensing data written in Julia. Using Julia's auto-differentiation ecosystem, LimberJack.jl can obtain gradients for its outputs up to an order of magnitude faster than traditional finite difference methods. This makes LimberJack.jl greatly synergistic with gradient-based sampling methods, such as Hamiltonian Monte Carlo, capable of efficiently exploring parameter spaces with hundreds of dimensions. We first prove LimberJack.jl's reliability by reanalysing the DES Y1 3$\times$2-point data. We then showcase its capabilities by using a O(100) parameters Gaussian Process to reconstruct the cosmic growth from a combination of DES Y1 galaxy clustering and weak lensing data, eBOSS QSO's, CMB lensing and redshift-space distortions. Our Gaussian process reconstruction of the growth factor is statistically consistent with the $Λ$CDM Planck 2018 prediction at all redshifts. Moreover, we show that the addition of RSD data is extremely beneficial to this type of analysis, reducing the uncertainty in the reconstructed growth factor by $20\%$ on average across redshift. LimberJack.jl is a fully open-source project available on Julia's general repository of packages and GitHub. △ Less

Submitted 15 March, 2024; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: Accepted to OJA, corrected bug displaying wrong Fig. 9

Cosmology from LOFAR Two-metre Sky Survey Data Release 2: Cross-correlation with the cosmic microwave background

Authors: S. J. Nakoneczny, D. Alonso, M. Bilicki, D. J. Schwarz, C. L. Hale, A. Pollo, C. Heneka, P. Tiwari, J. Zheng, M. Brüggen, M. J. Jarvis, T. W. Shimwell

Abstract: We combine the LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) second data release (DR2) catalogue with gravitational lensing maps from the Cosmic Microwave Background (CMB) to place constraints on the bias evolution of LoTSS radio galaxies, and on the amplitude of matter perturbations. We construct a flux-limited catalogue, and analyse its harmonic-space cross-correlation with CMB lensin… ▽ More We combine the LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) second data release (DR2) catalogue with gravitational lensing maps from the Cosmic Microwave Background (CMB) to place constraints on the bias evolution of LoTSS radio galaxies, and on the amplitude of matter perturbations. We construct a flux-limited catalogue, and analyse its harmonic-space cross-correlation with CMB lensing maps from Planck, $C_\ell^{gκ}$, as well as its auto-correlation, $C_\ell^{gg}$. We explore the models describing the redshift evolution of the large-scale radio galaxy bias, discriminating between them through the combination of both $C_\ell^{gκ}$ and $C_\ell^{gg}$. Fixing the bias evolution, we then use these data to place constraints on the amplitude of large scale density fluctuations. We report the significance of the $C_\ell^{gκ}$ signal at a level of $26.6σ$. We determine that a linear bias evolution of the form $b_g(z) = b_{g,D} / D(z)$, where $D(z)$ is the growth rate, is able to provide a good description of the data, and measure $b_{g,D} = 1.41 \pm 0.06$ for a sample flux-limited at $1.5\,{\rm mJy}$, for scales $\ell < 250$ for $C_\ell^{gg}$, and $\ell < 500$ for $C_\ell^{gκ}$. At the sample's median redshift, we obtain $b(z = 0.82) = 2.34 \pm 0.10$. Using $σ_8$ as a free parameter, while keeping other cosmological parameters fixed to the Planck values, we find fluctuations of $σ_8 = 0.75^{+0.05}_{-0.04}$. The result is in agreement with weak lensing surveys, and at $1σ$ difference with Planck CMB constraints. We also attempt to detect the late-time integrated Sachs-Wolfe effect with LOFAR, but with the current sky coverage, the cross-correlation with CMB temperature maps is consistent with zero. Our results are an important step towards constraining cosmology with radio continuum surveys from LOFAR and other future large radio surveys. △ Less

Submitted 15 May, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: The code is available at https://github.com/snakoneczny/cosmo-pipe

Journal ref: A&A, 681, A105 (2024)

Cosmology from LOFAR Two-metre Sky Survey Data Release 2: Angular Clustering of Radio Sources

Authors: C. L. Hale, D. J. Schwarz, P. N. Best, S. J. Nakoneczny, D. Alonso, D. Bacon, L. Böhme, N. Bhardwaj, M. Bilicki, S. Camera, C. S. Heneka, M. Pashapour-Ahmadabadi, P. Tiwari, J. Zheng, K. J. Duncan, M. J. Jarvis, R. Kondapally, M. Magliocchetti, H. J. A. Rottgering, T. W. Shimwell

Abstract: Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $μ$Jy beam$^{-1}$, the LOFAR Two-metre Sky Survey Data Release 2 (LoTSS-DR2) provides the largest low-frequency ($\sim$150 MHz) radio catalogue to date, making it an excellent tool for large-area radio cosmology studies. In this work, we use LoTSS-DR2 sources to investigate the angular two-point correlation function of galaxies wit… ▽ More Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $μ$Jy beam$^{-1}$, the LOFAR Two-metre Sky Survey Data Release 2 (LoTSS-DR2) provides the largest low-frequency ($\sim$150 MHz) radio catalogue to date, making it an excellent tool for large-area radio cosmology studies. In this work, we use LoTSS-DR2 sources to investigate the angular two-point correlation function of galaxies within the survey. We discuss systematics in the data and an improved methodology for generating random catalogues, compared to that used for LoTSS-DR1, before presenting the angular clustering for $\sim$900,000 sources $\geq$$1.5$ mJy and a peak signal-to-noise $\geq$$7.5$ across $\sim$$80\%$ of the observed area. Using the clustering we infer the bias assuming two evolutionary models. When fitting {angular scales of $0.5 \leqθ<5\,°$, using a linear bias model, we find LoTSS-DR2 sources are biased tracers of the underlying matter, with a bias of $b_{C}= 2.14^{+0.22}_{-0.20}$ (assuming constant bias) and $b_{E}(z=0)= 1.79^{+0.15}_{-0.14}$ (for an evolving model, inversely proportional to the growth factor), corresponding to $b_E= 2.81^{+0.24}_{-0.22}$ at the median redshift of our sample, assuming the LoTSS Deep Fields redshift distribution is representative of our data. This reduces to $b_{C}= 2.02^{+0.17}_{-0.16}$ and $b_{E}(z=0)= 1.67^{+0.12}_{-0.12}$ when allowing preferential redshift distributions from the Deep Fields to model our data. Whilst the clustering amplitude is slightly lower than LoTSS-DR1 ($\geq$2 mJy), our study benefits from larger samples and improved redshift estimates. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS. 29 pages, 24 figures

X-ray - cosmic shear cross-correlations: first detection and constraints on baryonic effects

Authors: Tassia Ferreira, David Alonso, Carlos Garcia-Garcia, Nora Elisa Chisari

Abstract: We report a first detection, at very high significance ($25σ$), of the cross-correlation between cosmic shear and the diffuse X-ray background, using data from the Dark Energy Survey and the ROSAT satellite. The X-ray cross-correlation signal is sensitive to the distribution of the surrounding gas in dark matter haloes. This allows us to use our measurements to place constraints on key physical pa… ▽ More We report a first detection, at very high significance ($25σ$), of the cross-correlation between cosmic shear and the diffuse X-ray background, using data from the Dark Energy Survey and the ROSAT satellite. The X-ray cross-correlation signal is sensitive to the distribution of the surrounding gas in dark matter haloes. This allows us to use our measurements to place constraints on key physical parameters that determine the impact of baryonic effects in the matter power spectrum. In particular, we determine the mass of haloes in which feedback has expelled half of their gas content on average to be $\log_{10}(M_c/M_\odot)=13.643^{+0.081}_{-0.12}$, and the polytropic index of the gas to be $Γ= 1.231^{+0.015}_{-0.011}$. This represents a first step in the direct use of X-ray cross-correlations to obtain improved constraints on cosmology and the physics of the intergalactic gas. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: 7 pages, 3 figures, 1 table

Impact of Galactic dust non-Gaussianity on searches for B-modes from inflation

Authors: Irene Abril-Cabezas, Carlos Hervías-Caimapo, Sebastian von Hausegger, Blake D. Sherwin, David Alonso

Abstract: A key challenge in the search for primordial B-modes is the presence of polarized Galactic foregrounds, especially thermal dust emission. Power-spectrum-based analysis methods generally assume the foregrounds to be Gaussian random fields when constructing a likelihood and computing the covariance matrix. In this paper, we investigate how non-Gaussianity in the dust field instead affects CMB and fo… ▽ More A key challenge in the search for primordial B-modes is the presence of polarized Galactic foregrounds, especially thermal dust emission. Power-spectrum-based analysis methods generally assume the foregrounds to be Gaussian random fields when constructing a likelihood and computing the covariance matrix. In this paper, we investigate how non-Gaussianity in the dust field instead affects CMB and foreground parameter inference in the context of inflationary B-mode searches, capturing this effect via modifications to the dust power-spectrum covariance matrix. For upcoming experiments such as the Simons Observatory, we find no dependence of the tensor-to-scalar ratio uncertainty $σ(r)$ on the degree of dust non-Gaussianity or the nature of the dust covariance matrix. We provide an explanation of this result, noting that when frequency decorrelation is negligible, dust in mid-frequency channels is cleaned using high-frequency data in a way that is independent of the spatial statistics of dust. We show that our results hold also for non-zero levels of frequency decorrelation that are compatible with existing data. We find, however, that neglecting the impact of dust non-Gaussianity in the covariance matrix can lead to inaccuracies in goodness-of-fit metrics. Care must thus be taken when using such metrics to test B-mode spectra and models, although we show that any such problems can be mitigated by using only cleaned spectrum combinations when computing goodness-of-fit statistics. △ Less

Submitted 20 December, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 8 figures, accepted version

Journal ref: MNRAS 527 (2024) 5751-5766

Constraints on dark matter and astrophysics from tomographic $γ$-ray cross-correlations

Authors: Anya Paopiamsap, David Alonso, Deaglan J. Bartlett, Maciej Bilicki

Abstract: We study the cross-correlation between maps of the unresolved $γ$-ray background constructed from the 12-year data release of the Fermi Large-Area Telescope, and the overdensity of galaxies in the redshift range $z\lesssim0.4$ as measured by the 2MASS Photometric Redshift survey and the WISE-SuperCOSMOS photometric survey. A signal is detected at the $8-10σ$ level, which we interpret in terms of b… ▽ More We study the cross-correlation between maps of the unresolved $γ$-ray background constructed from the 12-year data release of the Fermi Large-Area Telescope, and the overdensity of galaxies in the redshift range $z\lesssim0.4$ as measured by the 2MASS Photometric Redshift survey and the WISE-SuperCOSMOS photometric survey. A signal is detected at the $8-10σ$ level, which we interpret in terms of both astrophysical $γ$-ray sources, and WIMP dark matter decay and annihilation. The sensitivity achieved allows us to characterise the energy and redshift dependence of the signal, and we show that the latter is incompatible with a pure dark matter origin. We thus use our measurement to place an upper bound on the WIMP decay rate and the annihilation cross-section, finding constraints that are competitive with those found in other analyses. Our analysis is based on the extraction of clean model-independent observables that can then be used to constrain arbitrary astrophysical and particle physics models. In this sense we produce measurements of the $γ$-ray emissivity as a function of redshift and rest-frame energy $ε$, and of a quantity $F(ε)$ encapsulating all WIMP parameters relevant for dark matter decay or annihilation. We make these measurements, together with a full account of their statistical uncertainties, publicly available. △ Less

Submitted 9 May, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Comments: 27 pages, 13 figures. Data and code available in https://github.com/anyabua/FermiX

Galaxy bias in the era of LSST: perturbative bias expansions

Authors: Andrina Nicola, Boryana Hadzhiyska, Nathan Findlay, Carlos García-García, David Alonso, Anže Slosar, Zhiyuan Guo, Nickolas Kokron, Raúl Angulo, Alejandro Aviles, Jonathan Blazek, Jo Dunkley, Bhuvnesh Jain, Marcos Pellejero, James Sullivan, Christopher W. Walter, Matteo Zennaro

Abstract: Upcoming imaging surveys will allow for high signal-to-noise measurements of galaxy clustering at small scales. In this work, we present the results of the LSST bias challenge, the goal of which is to compare the performance of different nonlinear galaxy bias models in the context of LSST Y10 data. Specifically, we compare two perturbative approaches, Lagrangian perturbation theory (LPT) and Euler… ▽ More Upcoming imaging surveys will allow for high signal-to-noise measurements of galaxy clustering at small scales. In this work, we present the results of the LSST bias challenge, the goal of which is to compare the performance of different nonlinear galaxy bias models in the context of LSST Y10 data. Specifically, we compare two perturbative approaches, Lagrangian perturbation theory (LPT) and Eulerian PT (EPT) to two variants of Hybrid Effective Field Theory (HEFT), with our fiducial implementation of these models including terms up to second order in the bias expansion as well as nonlocal bias and deviations from Poissonian stochasticity. We consider different simulated galaxy samples and test the performance of the bias models in a tomographic joint analysis of LSST-Y10-like galaxy clustering, galaxy-galaxy-lensing and cosmic shear. We find both HEFT methods as well as LPT and EPT combined with non-perturbative predictions for the matter power spectrum to yield unbiased constraints on cosmological parameters up to at least a maximal scale of $k_{\mathrm{max}}=0.4 \; \mathrm{Mpc}^{-1}$ for all samples considered, even in the presence of assembly bias. While we find that we can reduce the complexity of the bias model for HEFT without compromising fit accuracy, this is not generally the case for the perturbative models. We find significant detections of non-Poissonian stochasticity in all cases considered, and our analysis shows evidence that small-scale galaxy clustering predominantly improves constraints on galaxy bias rather than cosmological parameters. These results therefore suggest that the systematic uncertainties associated with current nonlinear bias models are likely to be subdominant compared to other sources of error for tomographic analyses of upcoming photometric surveys, which bodes well for future galaxy clustering analyses using these high signal-to-noise data. [abridged] △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: 47 pages, 19 figures, 1 table, to be submitted to JCAP

Quaia, the Gaia-unWISE Quasar Catalog: An All-Sky Spectroscopic Quasar Sample

Authors: Kate Storey-Fisher, David W. Hogg, Hans-Walter Rix, Anna-Christina Eilers, Giulio Fabbian, Michael Blanton, David Alonso

Abstract: We present a new, all-sky quasar catalog, Quaia, that samples the largest comoving volume of any existing spectroscopic quasar sample. The catalog draws on the 6,649,162 quasar candidates identified by the Gaia mission that have redshift estimates from the space observatory's low-resolution BP/RP spectra. This initial sample is highly homogeneous and complete, but has low purity, and 18% of even t… ▽ More We present a new, all-sky quasar catalog, Quaia, that samples the largest comoving volume of any existing spectroscopic quasar sample. The catalog draws on the 6,649,162 quasar candidates identified by the Gaia mission that have redshift estimates from the space observatory's low-resolution BP/RP spectra. This initial sample is highly homogeneous and complete, but has low purity, and 18% of even the bright ($G<20.0$) confirmed quasars have discrepant redshift estimates ($|Δz/(1+z)|>0.2$) compared to those from the Sloan Digital Sky Survey (SDSS). In this work, we combine the Gaia candidates with unWISE infrared data (based on the Wide-field Infrared Survey Explorer survey) to construct a catalog useful for cosmological and astrophysical quasar studies. We apply cuts based on proper motions and Gaia and unWISE colors, reducing the number of contaminants by $\sim$4$\times$. We improve the redshifts by training a $k$-nearest neighbors model on SDSS redshifts, and achieve estimates on the $G<20.0$ sample with only 6% (10%) catastrophic errors with $|Δz/(1+z)|>0.2$ ($0.1$), a reduction of $\sim$3$\times$ ($\sim$2$\times$) compared to the Gaia redshifts. The final catalog has 1,295,502 quasars with $G<20.5$, and 755,850 candidates in an even cleaner $G<20.0$ sample, with accompanying rigorous selection function models. We compare Quaia to existing quasar catalogs, showing that its large effective volume makes it a highly competitive sample for cosmological large-scale structure analyses. The catalog is publicly available at https://zenodo.org/records/10403370. △ Less

Submitted 18 March, 2024; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: Published in the Astrophysical Journal; updated to match published version. Catalog available at https://doi.org/10.5281/zenodo.10403370. Animation available at https://cosmo.nyu.edu/ksf/quaia.mp4

Journal ref: Storey-Fisher (2024) 964, 69

Constraining cosmology with the Gaia-unWISE Quasar Catalog and CMB lensing: structure growth

Authors: David Alonso, Giulio Fabbian, Kate Storey-Fisher, Anna-Christina Eilers, Carlos García-García, David W. Hogg, Hans-Walter Rix

Abstract: We study the angular clustering of Quaia, a Gaia- and unWISE-based catalog of over a million quasars with an exceptionally well-defined selection function. With it, we derive cosmology constraints from the amplitude and growth of structure across cosmic time. We divide the sample into two redshift bins, centered at $z=1.0$ and $z=2.1$, and measure both overdensity auto-correlations and cross-corre… ▽ More We study the angular clustering of Quaia, a Gaia- and unWISE-based catalog of over a million quasars with an exceptionally well-defined selection function. With it, we derive cosmology constraints from the amplitude and growth of structure across cosmic time. We divide the sample into two redshift bins, centered at $z=1.0$ and $z=2.1$, and measure both overdensity auto-correlations and cross-correlations with maps of the Cosmic Microwave Background convergence measured by Planck. From these data, and including a prior from measurements of the baryon acoustic oscillations scale, we place constraints on the amplitude of the matter power spectrum $σ_8=0.766\pm 0.034$, and on the matter density parameter $Ω_m=0.343^{+0.017}_{-0.019}$. These measurements are in reasonable agreement with \planck at the $\sim$ 1.4$σ$ level, and are found to be robust with respect to observational and theoretical uncertainties. We find that our slightly lower value of $σ_8$ is driven by the higher-redshift sample, which favours a low amplitude of matter fluctuations. We present plausible arguments showing that this could be driven by contamination of the CMB lensing map by high-redshift extragalactic foregrounds, which should also affect other cross-correlations with tracers of large-scale structure beyond $z\sim1.5$. Our constraints are competitive with those from state-of-the-art 3$\times$2-point analyses, but arise from a range of scales and redshifts that is highly complementary to those covered by cosmic shear data and most galaxy clustering samples. This, coupled with the unprecedented combination of volume and redshift precision achieved by Quaia allows us to break the usual degeneracy between $Ω_m$ and $σ_8$. △ Less

Submitted 3 July, 2023; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: Submitted to JCAP, corrects links to public data release. The per-bin selection functions used in this analysis, and the complete set of power spectra, covariance, and associated metadata used in our fiducial analysis are publicly available at https://doi.org/10.5281/zenodo.8098635 . The public Quaia catalog can be found at https://doi.org/10.5281/zenodo.8060755

Can we constrain structure growth from galaxy proper motions?

Authors: Iain Duncan, David Alonso, Anže Slosar, Kate Storey-Fisher

Abstract: Galaxy peculiar velocities can be used to trace the growth of structure on cosmological scales. In the radial direction, peculiar velocities cause redshift space distortions, an established cosmological probe, and can be measured individually in the presence of an independent distance indicator. In the transverse direction, peculiar velocities cause proper motions. In this case, however, the prope… ▽ More Galaxy peculiar velocities can be used to trace the growth of structure on cosmological scales. In the radial direction, peculiar velocities cause redshift space distortions, an established cosmological probe, and can be measured individually in the presence of an independent distance indicator. In the transverse direction, peculiar velocities cause proper motions. In this case, however, the proper motions are too small to detect on a galaxy-by-galaxy basis for any realistic experiment in the foreseeable future, but could be detected statistically in cross-correlation with other tracers of the density fluctuations. We forecast the sensitivity for a detection of transverse peculiar velocities through the cross-correlation of a proper motion survey, modelled after existing extragalactic samples measured by Gaia, and an overlaping galaxy survey. In particular, we consider a low-redshift galaxy sample, and a higher-redshift quasar sample. We find that, while the expected cosmological signal is below the expected statistical uncertainties from current data using cross-correlations, the sensitivity can improve fast with future experiments, and the threshold for detection may not be too far away in the future. Quantitatively, we find that the signal-to-noise ratio for detection is in the range $S/N\sim0.3$, with most of the signal concentrated at low redshifts $z\lesssim0.3$. If detected, this signal is sensitive to the product of the expansion and growth rates at late times, and thus would constitute an independent observable, sensitive to both background expansion and large-scale density fluctuations. △ Less

Submitted 1 February, 2024; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 11 pages, 6 figures

The Simons Observatory: Beam characterization for the Small Aperture Telescopes

Authors: Nadia Dachlythra, Adriaan J. Duivenvoorden, Jon E. Gudmundsson, Matthew Hasselfield, Gabriele Coppi, Alexandre E. Adler, David Alonso, Susanna Azzoni, Grace E. Chesmore, Giulio Fabbian, Ken Ganga, Remington G. Gerras, Andrew H. Jaffe, Bradley R. Johnson, Brian Keating, Reijo Keskitalo, Theodore S. Kisner, Nicoletta Krachmalnicoff, Marius Lungu, Frederick Matsuda, Sigurd Naess, Lyman Page, Roberto Puddu, Giuseppe Puglisi, Sara M. Simon , et al. (5 additional authors not shown)

Abstract: We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a map maker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the map maker. We test our reconstruction performance for four different freque… ▽ More We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a map maker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the map maker. We test our reconstruction performance for four different frequency bands against various algorithmic parameters, atmospheric conditions and input beams. We additionally show the reconstruction quality as function of the number of available observations and investigate how different calibration strategies affect the beam uncertainty. For all of the cases considered, we find good agreement between the fitted results and the input beam model within a ~1.5% error for a multipole range l = 30 - 700 and an ~0.5% error for a multipole range l = 50 - 200. We conclude by using a harmonic-domain component separation algorithm to verify that the beam reconstruction errors and biases observed in our analysis do not significantly bias the Simons Observatory r-measurement. △ Less

Submitted 7 May, 2024; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 23 pages, 21 figures, published in ApJ

Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Power Spectra

Authors: Roohi Dalal, Xiangchong Li, Andrina Nicola, Joe Zuntz, Michael A. Strauss, Sunao Sugiyama, Tianqing Zhang, Markus M. Rau, Rachel Mandelbaum, Masahiro Takada, Surhud More, Hironao Miyatake, Arun Kannawadi, Masato Shirasaki, Takanori Taniguchi, Ryuichi Takahashi, Ken Osato, Takashi Hamana, Masamune Oguri, Atsushi J. Nishizawa, Andrés A. Plazas Malagón, Tomomi Sunayama, David Alonso, Anže Slosar, Robert Armstrong , et al. (13 additional authors not shown)

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

Submitted 4 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 35 pages, 18 figures, 6 tables, for coordinated submission to PRD with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

Journal ref: Physical Review D, Volume 108, Issue 12, December 2023, article id.123519

Science with the Einstein Telescope: a comparison of different designs

Authors: Marica Branchesi, Michele Maggiore, David Alonso, Charles Badger, Biswajit Banerjee, Freija Beirnaert, Enis Belgacem, Swetha Bhagwat, Guillaume Boileau, Ssohrab Borhanian, Daniel David Brown, Man Leong Chan, Giulia Cusin, Stefan L. Danilishin, Jerome Degallaix, Valerio De Luca, Arnab Dhani, Tim Dietrich, Ulyana Dupletsa, Stefano Foffa, Gabriele Franciolini, Andreas Freise, Gianluca Gemme, Boris Goncharov, Archisman Ghosh , et al. (51 additional authors not shown)

Abstract: The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogeni… ▽ More The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogenic temperature. Here, we examine the scientific perspectives under possible variations of this reference design. We perform a detailed evaluation of the science case for a single triangular geometry observatory, and we compare it with the results obtained for a network of two L-shaped detectors (either parallel or misaligned) located in Europe, considering different choices of arm-length for both the triangle and the 2L geometries. We also study how the science output changes in the absence of the low-frequency instrument, both for the triangle and the 2L configurations. We examine a broad class of simple `metrics' that quantify the science output, related to compact binary coalescences, multi-messenger astronomy and stochastic backgrounds, and we then examine the impact of different detector designs on a more specific set of scientific objectives. △ Less

Submitted 17 June, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 197 pages, 73 figures. v2: corrections in the part on the sensitivity to stochastic backgrounds. Accepted in JCAP

Report number: ET-0084A-23

First measurement of the nuclear-recoil ionization yield in silicon at 100 eV

Authors: M. F. Albakry, I. Alkhatib, D. Alonso, D. W. P. Amaral, P. An, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, P. S. Barbeau, C. Bathurst, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott , et al. (115 additional authors not shown)

Abstract: We measured the nuclear--recoil ionization yield in silicon with a cryogenic phonon-sensitive gram-scale detector. Neutrons from a mono-energetic beam scatter off of the silicon nuclei at angles corresponding to energy depositions from 4\,keV down to 100\,eV, the lowest energy probed so far. The results show no sign of an ionization production threshold above 100\,eV. These results call for furthe… ▽ More We measured the nuclear--recoil ionization yield in silicon with a cryogenic phonon-sensitive gram-scale detector. Neutrons from a mono-energetic beam scatter off of the silicon nuclei at angles corresponding to energy depositions from 4\,keV down to 100\,eV, the lowest energy probed so far. The results show no sign of an ionization production threshold above 100\,eV. These results call for further investigation of the ionization yield theory and a comprehensive determination of the detector response function at energies below the keV scale. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Journal ref: Physical Review Letters 131.9 (2023): 091801

A Search for Low-mass Dark Matter via Bremsstrahlung Radiation and the Migdal Effect in SuperCDMS

Authors: M. F. Albakry, I. Alkhatib, D. Alonso, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott, J. Cooley, H. Coombes , et al. (108 additional authors not shown)

Abstract: We present a new analysis of previously published of SuperCDMS data using a profile likelihood framework to search for sub-GeV dark matter (DM) particles through two inelastic scattering channels: bremsstrahlung radiation and the Migdal effect. By considering these possible inelastic scattering channels, experimental sensitivity can be extended to DM masses that are undetectable through the DM-nuc… ▽ More We present a new analysis of previously published of SuperCDMS data using a profile likelihood framework to search for sub-GeV dark matter (DM) particles through two inelastic scattering channels: bremsstrahlung radiation and the Migdal effect. By considering these possible inelastic scattering channels, experimental sensitivity can be extended to DM masses that are undetectable through the DM-nucleon elastic scattering channel, given the energy threshold of current experiments. We exclude DM masses down to $220~\textrm{MeV}/c^2$ at $2.7 \times 10^{-30}~\textrm{cm}^2$ via the bremsstrahlung channel. The Migdal channel search provides overall considerably more stringent limits and excludes DM masses down to $30~\textrm{MeV}/c^2$ at $5.0 \times 10^{-30}~\textrm{cm}^2$. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Comments: Submitted to PRD

Report number: 112013

Journal ref: Phys. Rev. D 107, 2023

The Simons Observatory: pipeline comparison and validation for large-scale B-modes

Authors: K. Wolz, S. Azzoni, C. Hervias-Caimapo, J. Errard, N. Krachmalnicoff, D. Alonso, C. Baccigalupi, A. Baleato Lizancos, M. L. Brown, E. Calabrese, J. Chluba, J. Dunkley, G. Fabbian, N. Galitzki, B. Jost, M. Morshed, F. Nati

Abstract: The upcoming Simons Observatory Small Aperture Telescopes aim at achieving a constraint on the primordial tensor-to-scalar ratio $r$ at the level of $σ(r=0)\lesssim0.003$, observing the polarized CMB in the presence of partial sky coverage, cosmic variance, inhomogeneous non-white noise, and Galactic foregrounds. We present three different analysis pipelines able to constrain $r$ given the latest… ▽ More The upcoming Simons Observatory Small Aperture Telescopes aim at achieving a constraint on the primordial tensor-to-scalar ratio $r$ at the level of $σ(r=0)\lesssim0.003$, observing the polarized CMB in the presence of partial sky coverage, cosmic variance, inhomogeneous non-white noise, and Galactic foregrounds. We present three different analysis pipelines able to constrain $r$ given the latest available instrument performance, and compare their predictions on a set of sky simulations that allow us to explore a number of Galactic foreground models and elements of instrumental noise, relevant for the Simons Observatory. The three pipelines employ different combinations of parametric and non-parametric component separation at the map and power spectrum levels, and use B-mode purification to estimate the CMB B-mode power spectrum. We applied them to a common set of simulated realistic frequency maps, and compared and validated them with focus on their ability to extract robust constraints on the tensor-to-scalar ratio $r$. We evaluated their performance in terms of bias and statistical uncertainty on this parameter. In most of the scenarios the three methodologies achieve similar performance. Nevertheless, several simulations with complex foreground signals lead to a $>2σ$ bias on $r$ if analyzed with the default versions of these pipelines, highlighting the need for more sophisticated pipeline components that marginalize over foreground residuals. We show two such extensions, using power-spectrum-based and map-based methods, that are able to fully reduce the bias on $r$ below the statistical uncertainties in all foreground models explored, at a moderate cost in terms of $σ(r)$. △ Less

Submitted 9 July, 2024; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: 20 pages, 15 figures

Journal ref: A&A 686, A16 (2024)

Analytical marginalisation over photometric redshift uncertainties in cosmic shear analyses

Authors: Jaime Ruiz-Zapatero, Boryana Hadzhiyska, David Alonso, Pedro G. Ferreira, Carlos García-García, Arrykrishna Mootoovaloo

Abstract: As the statistical power of imaging surveys grows, it is crucial to account for all systematic uncertainties. This is normally done by constructing a model of these uncertainties and then marginalizing over the additional model parameters. The resulting high dimensionality of the total parameter spaces makes inferring the cosmological parameters significantly more costly using traditional Monte-Ca… ▽ More As the statistical power of imaging surveys grows, it is crucial to account for all systematic uncertainties. This is normally done by constructing a model of these uncertainties and then marginalizing over the additional model parameters. The resulting high dimensionality of the total parameter spaces makes inferring the cosmological parameters significantly more costly using traditional Monte-Carlo sampling methods. A particularly relevant example is the redshift distribution, $p(z)$, of the source samples, which may require tens of parameters to describe fully. However, relatively tight priors can be usually placed on these parameters through calibration of the associated systematics. In this paper we show, quantitatively, that a linearisation of the theoretical prediction with respect to these calibratable systematic parameters allows us to analytically marginalise over these extra parameters, leading to a factor $\sim30$ reduction in the time needed for parameter inference, while accurately recovering the same posterior distributions for the cosmological parameters that would be obtained through a full numerical marginalisation over 160 $p(z)$ parameters. We demonstrate that this is feasible not only with current data and current achievable calibration priors but also for future Stage-IV datasets. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 11 pages, 8 figures, prepared for submission to MNRAS, comments welcome

Cosmology with 6 parameters in the Stage-IV era: efficient marginalisation over nuisance parameters

Authors: Boryana Hadzhiyska, Kevin Wolz, Susanna Azzoni, David Alonso, Carlos García-García, Jaime Ruiz-Zapatero, Anže Slosar

Abstract: The analysis of photometric large-scale structure data is often complicated by the need to account for many observational and astrophysical systematics. The elaborate models needed to describe them often introduce many ``nuisance parameters'', which can be a major inhibitor of an efficient parameter inference. In this paper we introduce an approximate method to analytically marginalise over a larg… ▽ More The analysis of photometric large-scale structure data is often complicated by the need to account for many observational and astrophysical systematics. The elaborate models needed to describe them often introduce many ``nuisance parameters'', which can be a major inhibitor of an efficient parameter inference. In this paper we introduce an approximate method to analytically marginalise over a large number of nuisance parameters based on the Laplace approximation. We discuss the mathematics of the method, its relation to concepts such as volume effects and profile likelihood, and show that it can be further simplified for calibratable systematics by linearising the dependence of the theory on the associated parameters. We quantify the accuracy of this approach by comparing it with traditional sampling methods in the context of existing data from the Dark Energy Survey, as well as futuristic Stage-IV photometric data. The linearised version of the method is able to obtain parameter constraints that are virtually equivalent to those found by exploring the full parameter space for a large number of calibratable nuisance parameters, while reducing the computation time by a factor 3-10. Furthermore, the non-linearised approach is able to analytically marginalise over a large number of parameters, returning constraints that are virtually indistinguishable from the brute-force method in most cases, accurately reproducing both the marginalised uncertainty on cosmological parameters, and the impact of volume effects associated with this marginalisation. We provide simple recipes to diagnose when the approximations made by the method fail and one should thus resort to traditional methods. The gains in sampling efficiency associated with this method enable the joint analysis of multiple surveys, typically hindered by the large number of nuisance parameters needed to describe them. △ Less

Submitted 13 July, 2023; v1 submitted 27 January, 2023; originally announced January 2023.

Comments: 17 pages, 6 figures, 3 tables

Single energy measurement Integral Fluctuation theorem and non-projective measurements

Authors: Daniel Alonso, Antonia Ruiz García

Abstract: We study a Jarzysnki type equality for work in systems that are monitored using non-projective unsharp measurements. The information acquired by the observer from the outcome $f$ of an energy measurement, and the subsequent conditioned normalized state $\hat ρ(t,f)$ evolved up to a final time $t$ are used to define work, as the difference between the final expectation value of the energy and the r… ▽ More We study a Jarzysnki type equality for work in systems that are monitored using non-projective unsharp measurements. The information acquired by the observer from the outcome $f$ of an energy measurement, and the subsequent conditioned normalized state $\hat ρ(t,f)$ evolved up to a final time $t$ are used to define work, as the difference between the final expectation value of the energy and the result $f$ of the measurement. The Jarzynski equality obtained depends on the coherences that the state develops during the process, the characteristics of the meter used to measure the energy, and the noise it induces into the system. We analyze those contributions in some detail to unveil their role. We show that in very particular cases, but not in general, the effect of such noise gives a factor multiplying the result that would be obtained if projective measurements were used instead of non-projective ones. The unsharp character of the measurements used to monitor the energy of the system, which defines the resolution of the meter, leads to different scenarios of interest. In particular, if the distance between neighboring elements in the energy spectrum is much larger than the resolution of the meter, then a similar result to the projective measurement case is obtained, up to a multiplicative factor that depends on the meter. A more subtle situation arises in the opposite case in which measurements may be non-informative, i.e. they may not contribute to update the information about the system. In this case, a correction to the relation obtained in the non-overlapping case appears. We analyze the conditions in which such a correction becomes negligible. We also study the coherences, in terms of the relative entropy of coherence developed by the evolved post-measurement state. We illustrate the results by analyzing a two-level system monitored by a simple meter. △ Less

Submitted 26 December, 2022; originally announced December 2022.

Comments: 14 pages, 3 figures

The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST

Authors: J. Prat, J. Zuntz, Y. Omori, C. Chang, T. Tröster, E. Pedersen, C. García-García, E. Phillips-Longley, J. Sanchez, D. Alonso, X. Fang, E. Gawiser, K. Heitmann, M. Ishak, M. Jarvis, E. Kovacs, P. Larsen, Y. -Y. Mao, L. Medina Varela, M. Paterno, S. D. Vitenti, Z. Zhang, The LSST Dark Energy Science Collaboration

Abstract: We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Sci… ▽ More We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis. △ Less

Submitted 21 April, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 20+11 pages, 10+11 figures. Version accepted in The Open Journal of Astrophysics

The N5K Challenge: Non-Limber Integration for LSST Cosmology

Authors: C. D. Leonard, T. Ferreira, X. Fang, R. Reischke, N. Schoeneberg, T. Tröster, D. Alonso, J. E. Campagne, F. Lanusse, A. Slosar, M. Ishak, the LSST Dark Energy Science Collaboration

Abstract: The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angula… ▽ More The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angular power spectrum of galaxy number counts and cosmic shear data without invoking the so-called 'Limber approximation', in the context of the Rubin Observatory Legacy Survey of Space and Time (LSST). We quantify the performance, in terms of accuracy and speed, of three non-Limber implementations: ${\tt FKEM (CosmoLike)}$, ${\tt Levin}$, and ${\tt matter}$, themselves based on different integration schemes and approximations. We find that in the challenge's fiducial 3x2pt LSST Year 10 scenario, ${\tt FKEM (CosmoLike)}$ produces the fastest run time within the required accuracy by a considerable margin, positioning it favourably for use in Bayesian parameter inference. This method, however, requires further development and testing to extend its use to certain analysis scenarios, particularly those involving a scale-dependent growth rate. For this and other reasons discussed herein, alternative approaches such as ${\tt matter}$ and ${\tt Levin}$ may be necessary for a full exploration of parameter space. We also find that the usual first-order Limber approximation is insufficiently accurate for LSST Year 10 3x2pt analysis on $\ell=200-1000$, whereas invoking the second-order Limber approximation on these scales (with a full non-Limber method at smaller $\ell$) does suffice. △ Less

Submitted 14 February, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 15 pages, 8 figures. Journal-accepted version

Development, manufacturing and testing of small launcher structures from Portugal

Authors: André G. C. Guerra, Daniel Alonso, Catarina Silva, Alexander Costa, Joaquim Rocha, Luis Colaço, Sandra Fortuna, Tiago Pires, Luis Pinheiro, Nuno Carneiro, André João, Gonçalo Araújo, Pedro Meireles, Stephan Schmid

Abstract: During the last decades the industry has seen the number of Earth orbiting satellites rise, mostly due to the need to monitor Earth as well as to establish global communication networks. Nano, micro, and small satellites have been a prime tool for answering these needs, with large and mega constellations planned, leading to a potential launch gap. An effective and commercially appealing solution i… ▽ More During the last decades the industry has seen the number of Earth orbiting satellites rise, mostly due to the need to monitor Earth as well as to establish global communication networks. Nano, micro, and small satellites have been a prime tool for answering these needs, with large and mega constellations planned, leading to a potential launch gap. An effective and commercially appealing solution is the development of small launchers, as these can complement the current available launch opportunity offer, serving a large pool of different types of clients, with a flexible and custom service that large conventional launchers cannot adequately assure. Rocket Factory Augsburg has partnered with CEiiA for the development of several structures for the RFA One rocket. The objective has been the design of solutions that are low-cost, light, and custom-made, applying design and manufacturing concepts as well as technologies from other industries, like the aeronautical and automotive, to the aerospace one. This allows for the implementation of a New Space approach to the launcher segment, while also building a supply chain and a set of solutions that enables the industrialisation of such structures for this and future small launchers. The two main systems under development have been a versatile Kick-Stage, for payload carrying and orbit insertion, and a sturdy Payload Fairing. Even though the use of components off-the-shelf have been widely accepted in the space industry for satellites, these two systems pose different challenges as they must be: highly reliable during the most extreme conditions imposed by the launch, so that they can be considered safe to launch all types of payloads. This paper thus dives deep on the solutions developed in the last few years, presenting also lessons learned during the manufacturing and testing of these structures. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 10 pages, 13 figures, Manuscript presented at the 73rd International Astronautical Congress, IAC 2022, Paris, France, 18 - 22 September 2022

A hybrid map-$C_\ell$ component separation method for primordial CMB $B$-mode searches

Authors: Susanna Azzoni, David Alonso, Maximilian H. Abitbol, Josquin Errard, Nicoletta Krachmalnicoff

Abstract: The observation of the polarised emission from the Cosmic Microwave Background (CMB) from future ground-based and satellite-borne experiments holds the promise of indirectly detecting the elusive signal from primordial tensor fluctuations in the form of large-scale $B$-mode polarisation. Doing so, however, requires an accurate and robust separation of the signal from polarised Galactic foregrounds… ▽ More The observation of the polarised emission from the Cosmic Microwave Background (CMB) from future ground-based and satellite-borne experiments holds the promise of indirectly detecting the elusive signal from primordial tensor fluctuations in the form of large-scale $B$-mode polarisation. Doing so, however, requires an accurate and robust separation of the signal from polarised Galactic foregrounds. We present a component separation method for multi-frequency CMB observations that combines some of the advantages of map-based and power-spectrum-based techniques, and which is direcly applicable to data in the presence of realistic foregrounds and instrumental noise. We demonstrate that the method is able to reduce the contamination from Galactic foregrounds below an equivalent tensor-to-scalar ratio $r_{\rm FG}\lesssim5\times10^{-4}$, as required for next-generation observatories, for a wide range of foreground models with varying degrees of complexity. This bias reduction is associated with a mild $\sim20-30\%$ increase in the final statistical uncertainties, and holds for large sky areas, and for experiments targeting both the reionisation and recombination bumps in the $B$-mode power spectrum. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: 28 pages, 9 figures

Combining cosmic shear data with correlated photo-$z$ uncertainties: constraints from DESY1 and HSC-DR1

Authors: Carlos García-García, David Alonso, Pedro G. Ferreira, Boryana Hadzhiyska, Andrina Nicola, Carles Sánchez, Anže Slosar

Abstract: An accurate calibration of the source redshift distribution $p(z)$ is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often u… ▽ More An accurate calibration of the source redshift distribution $p(z)$ is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often use the same samples for redshift calibration. This introduces a source of statistical and systematic uncertainty that is highly correlated across different weak lensing datasets, and which must be accurately characterised and propagated in order to obtain robust cosmological constraints from their combination. In this paper we introduce a method to quantify and propagate the uncertainties on the source redshift distribution in two different surveys sharing the same calibrating sample. The method is based on an approximate analytical marginalisation of the $p(z)$ statistical uncertainties and the correlated marginalisation of residual systematics. We apply this method to the combined analysis of cosmic shear data from the DESY1 data release and the HSC-DR1 data, using the COSMOS 30-band catalog as a common redshift calibration sample. We find that, although there is significant correlation in the uncertainties on the redshift distributions of both samples, this does not change the final constraints on cosmological parameters significantly. The same is true also for the impact of residual systematic uncertainties from the errors in the COSMOS 30-band photometric redshifts. Additionally, we show that these effects will still be negligible in Stage-IV datasets. Finally, the combination of DESY1 and HSC-DR1 allows us to constrain the ``clumpiness'' parameter to $S_8 = 0.768^{+0.021}_{-0.017}$. This corresponds to a $\sim\sqrt{2}$ improvement in uncertainties with respect to either DES or HSC alone. △ Less

Submitted 21 December, 2022; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 32 pages, 12 figures, comments welcome; v2: minor changes, accepted version in JCAP

Generalised Gillespie Algorithms for Simulations in a Rule-Based Epidemiological Model Framework

Authors: David Alonso, Steffen Bauer, Markus Kirkilionis, Lisa Maria Kreusser, Luca Sbano

Abstract: Rule-based models have been successfully used to represent different aspects of the COVID-19 pandemic, including age, testing, hospitalisation, lockdowns, immunity, infectivity, behaviour, mobility and vaccination of individuals. These rule-based approaches are motivated by chemical reaction rules which are traditionally solved numerically with the standard Gillespie algorithm proposed in the cont… ▽ More Rule-based models have been successfully used to represent different aspects of the COVID-19 pandemic, including age, testing, hospitalisation, lockdowns, immunity, infectivity, behaviour, mobility and vaccination of individuals. These rule-based approaches are motivated by chemical reaction rules which are traditionally solved numerically with the standard Gillespie algorithm proposed in the context of molecular dynamics. When applying reaction system type of approaches to epidemiology, generalisations of the Gillespie algorithm are required due to the time-dependency of the problems. In this article, we present different generalisations of the standard Gillespie algorithm which address discrete subtypes (e.g., incorporating the age structure of the population), time-discrete updates (e.g., incorporating daily imposed change of rates for lockdowns) and deterministic delays (e.g., given waiting time until a specific change in types such as release from isolation occurs). These algorithms are complemented by relevant examples in the context of the COVID-19 pandemic and numerical results. △ Less

Submitted 24 October, 2022; v1 submitted 17 October, 2022; originally announced October 2022.

MSC Class: 92B05; 60G07; ACM Class: J.3; G.3

Constraining the physics of star formation from CIB-cosmic shear cross-correlations

Authors: Baptiste Jego, David Alonso, Carlos García-García, Jaime Ruiz-Zapatero

Abstract: Understanding the physics of star formation is one of the key problems facing modern astrophysics. The Cosmic Infrared Background (CIB), sourced by the emission from all dusty star-forming galaxies since the epoch of reionisation, is a complementary probe to study the star formation history, as well as an important extragalactic foreground for studies of the Cosmic Microwave Background (CMB). Unde… ▽ More Understanding the physics of star formation is one of the key problems facing modern astrophysics. The Cosmic Infrared Background (CIB), sourced by the emission from all dusty star-forming galaxies since the epoch of reionisation, is a complementary probe to study the star formation history, as well as an important extragalactic foreground for studies of the Cosmic Microwave Background (CMB). Understanding the physics of the CIB is therefore of high importance for both cosmology and galaxy formation studies. In this paper, we make high signal-to-noise measurements of the cross-correlation between maps of the CIB from the Planck experiment, and cosmic shear measurements from the Dark Energy Survey and Kilo-Degree Survey. Cosmic shear, sourced mainly by the weak gravitational lensing of photons emitted by background galaxies, is a direct tracer of the matter distribution, and thus we can use its cross-correlation with the CIB to directly test our understanding of the link between the star formation rate (SFR) density and the matter density. We use our measurements to place constraints on a halo-based model of the SFR that parametrises the efficiency with which gas is transformed into stars as a function of halo mass and redshift. These constraints are enhanced by combining our data with model-independent measurements of the bias-weighted SFR density extracted from the tomographic cross-correlation of galaxies and the CIB. We are able to place constraints on the peak efficiency at low redshifts, $η=0.445^{+0.055}_{-0.11}$, and on the halo mass at which this peak efficiency is achieved today $\log_{10}(M_1/M_\odot) = 12.17\pm0.25$. Our constraints are in excellent agreement with direct measurements of the SFR density, as well as other CIB-based studies. △ Less

Submitted 13 January, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: Not too many pages. Comments welcome

The impact of the Universe's expansion rate on constraints on modified growth of structure

Authors: Jaime Ruiz-Zapatero, David Alonso, Pedro G. Ferreira, Carlos Garcia-Garcia

Abstract: In the context of modified gravity, at linear level, the growth of structure in the Universe will be affected by modifications to the Poisson equation and by the background expansion rate of the Universe. It has been shown that these two effects lead to a degeneracy which must be properly accounted for if one is to place reliable constraints on new forces on large scales or, equivalently, modifica… ▽ More In the context of modified gravity, at linear level, the growth of structure in the Universe will be affected by modifications to the Poisson equation and by the background expansion rate of the Universe. It has been shown that these two effects lead to a degeneracy which must be properly accounted for if one is to place reliable constraints on new forces on large scales or, equivalently, modifications to General Relativity. In this paper we show that current constraints are such that assumptions about the background expansion have little impact on constraints on modifications to gravity. We do so by considering the background of a $Λ$ Cold Dark Matter ($Λ$CDM) universe, a universe with a more general equation of state for the dark energy, and finally, a general, model-independent, expansion rate. We use Gaussian Processes to model modifications to Poisson's equation and, in the case of a general expansion rate, to model the redshift dependent Hubble rate. We identify a degeneracy between modifications to Poisson's equation and the background matter density, $Ω_M$, which can only be broken by assuming a model-dependent expansion rate. We show that, with current data, the constraints on modifications to the Poisson equation via measurements of the growth rate range between $10-20\%$ depending on the strength of our assumptions on the Universe's expansion rate. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Comments: Prepared for submission to PRD. Comments welcome

The star formation history in the last 10 billion years from CIB cross-correlations

Authors: Baptiste Jego, Jaime Ruiz-Zapatero, Carlos García-García, Nick Koukoufilippas, David Alonso

Abstract: The Cosmic Infrared Background (CIB) traces the emission of star-forming galaxies throughout all cosmic epochs. Breaking down the contribution from galaxies at different redshifts to the observed CIB maps would allow us to probe the history of star formation. In this paper, we cross-correlate maps of the CIB with galaxy samples covering the range $z\lesssim2$ to measure the bias-weighted star-form… ▽ More The Cosmic Infrared Background (CIB) traces the emission of star-forming galaxies throughout all cosmic epochs. Breaking down the contribution from galaxies at different redshifts to the observed CIB maps would allow us to probe the history of star formation. In this paper, we cross-correlate maps of the CIB with galaxy samples covering the range $z\lesssim2$ to measure the bias-weighted star-formation rate (SFR) density $\langle bρ_{\rm SFR}\rangle$ as a function of time in a model independent way. This quantity is complementary to direct measurements of the SFR density $ρ_{\rm SFR}$, giving a higher weight to more massive haloes, and thus provides additional information to constrain the physical properties of star formation. Using cross-correlations of the CIB with galaxies from the DESI Legacy Survey and the extended Baryon Oscillation Spectroscopic Survey, we obtain high signal-to-noise ratio measurements of $\langle bρ_{\rm SFR}\rangle$, which we then use to place constraints on halo-based models of the star-formation history. We fit halo-based SFR models to our data and compare the recovered $ρ_{\rm SFR}$ with direct measurements of this quantity. We find a qualitatively good agreement between both independent datasets, although the details depend on the specific halo model assumed. This constitutes a useful robustness test for the physical interpretation of the CIB, and reinforces the role of CIB maps as valuable astrophysical probes of the large-scale structure. We report our measurements of $\langle bρ_{\rm SFR}\rangle$ as well as a thorough account of their statistical uncertainties, which can be used to constrain star formation models in combination with other data. △ Less

Submitted 30 June, 2022; originally announced June 2022.

Comments: 17 pages, 11 figures, comments welcome

First measurement of projected phase correlations and large-scale structure constraints

Authors: Felipe Oliveira Franco, Boryana Hadzhiyska, David Alonso

Abstract: Phase correlations are an efficient way to extract astrophysical information that is largely independent from the power spectrum. We develop an estimator for the line correlation function (LCF) of projected fields, given by the correlation between the harmonic-space phases at three equidistant points on a great circle. We make a first, 6.5$σ$ measurement of phase correlations on data from the 2MPZ… ▽ More Phase correlations are an efficient way to extract astrophysical information that is largely independent from the power spectrum. We develop an estimator for the line correlation function (LCF) of projected fields, given by the correlation between the harmonic-space phases at three equidistant points on a great circle. We make a first, 6.5$σ$ measurement of phase correlations on data from the 2MPZ survey. Finally, we show that the LCF can significantly improve constraints on parameters describing the galaxy-halo connection that are typically degenerate using only two-point data. △ Less

Submitted 22 June, 2022; originally announced June 2022.

The CARMENES search for exoplanets around M dwarfs. Rotational variation in activity indicators of Ross 318, YZ CMi, TYC 3529-1437-1, and EV Lac

Authors: P. Schöfer, S. V. Jeffers, A. Reiners, M. Zechmeister, B. Fuhrmeister, M. Lafarga, I. Ribas, A. Quirrenbach, P. J. Amado, J. A. Caballero, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, M. Cortés-Contreras, E. Díez Alonso, S. Dreizler, E. W. Guenther, O. Herbort, E. N. Johnson, A. Kaminski, M. Kürster, D. Montes, J. C. Morales, S. Pedraz, L. Tal-Or

Abstract: The CARMENES instrument is searching for periodic radial-velocity (RV) variations of M dwarfs, which are induced by orbiting planets. However, there are other potential sources of such variations, including rotational modulation caused by stellar activity. We aim to investigate four M dwarfs (Ross 318, YZ CMi, TYC 3529-1437-1, and EV Lac) with different activity levels and spectral sub-types. Our… ▽ More The CARMENES instrument is searching for periodic radial-velocity (RV) variations of M dwarfs, which are induced by orbiting planets. However, there are other potential sources of such variations, including rotational modulation caused by stellar activity. We aim to investigate four M dwarfs (Ross 318, YZ CMi, TYC 3529-1437-1, and EV Lac) with different activity levels and spectral sub-types. Our goal is to compare the periodicities seen in 22 activity indicators and the stellar RVs, and to examine their stability over time. For each star, we calculated GLS periodograms of pseudo-equivalent widths of chromospheric lines, indices of photospheric bands, the differential line width as a measure of the width of the average photospheric absorption line, the RV, the chromatic index that describes the wavelength dependence of the RV, and CCF parameters. We also calculated periodograms for subsets of the data and compared our results to TESS photometry. We find the rotation periods of all four stars to manifest themselves in the RV and photospheric indicators, particularly the TiO~7050 index, whereas the chromospheric lines show clear signals only at lower activity levels. For EV Lac and TYC 3529-1437-1, we find episodes during which indicators vary with the rotation period, and episodes during which they vary with half the rotation period, similarly to photometric light curves. The changing periodicities reflect the evolution of stellar activity features on the stellar surface. We therefore conclude that our results not only emphasise the importance of carefully analysing indicators complementary to the RV in RV surveys, but they also suggest that it is also useful to search for signals in activity indicators in subsets of the dataset, because an activity signal that is present in the RV may not be visible in the activity indicators all the time, in particular for the most active stars. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 8 pages + 9 pages appendix, 5+4 figures, 1+4 tables, accepted for publication in A&A

Journal ref: A&A 663, A68 (2022)

Snowmass 2021 CMB-S4 White Paper

Authors: Kevork Abazajian, Arwa Abdulghafour, Graeme E. Addison, Peter Adshead, Zeeshan Ahmed, Marco Ajello, Daniel Akerib, Steven W. Allen, David Alonso, Marcelo Alvarez, Mustafa A. Amin, Mandana Amiri, Adam Anderson, Behzad Ansarinejad, Melanie Archipley, Kam S. Arnold, Matt Ashby, Han Aung, Carlo Baccigalupi, Carina Baker, Abhishek Bakshi, Debbie Bard, Denis Barkats, Darcy Barron, Peter S. Barry , et al. (331 additional authors not shown)

Abstract: This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan. This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

Propagating spatially-varying multiplicative shear bias to cosmological parameter estimation for stage-IV weak-lensing surveys

Authors: Casey Cragg, Christopher A. J. Duncan, Lance Miller, David Alonso

Abstract: We consider the bias introduced by a spatially-varying multiplicative shear bias (m-bias) on tomographic cosmic shear angular power spectra. To compute the bias in the power spectra, we estimate the mode-coupling matrix associated with an m-bias map using a computationally-efficient pseudo-Cl method. This allows us to consider the effect of the m-bias to high l. We then conduct a Fisher matrix ana… ▽ More We consider the bias introduced by a spatially-varying multiplicative shear bias (m-bias) on tomographic cosmic shear angular power spectra. To compute the bias in the power spectra, we estimate the mode-coupling matrix associated with an m-bias map using a computationally-efficient pseudo-Cl method. This allows us to consider the effect of the m-bias to high l. We then conduct a Fisher matrix analysis to forecast resulting biases in cosmological parameters. For a Euclid-like survey with a spatially-varying m-bias, with zero mean and rms of 0.01, we find that parameter biases reach a maximum of ~10% of the expected statistical error, if multipoles up to l_max = 5000 are included. We conclude that the effect of the spatially-varying m-bias may be a sub-dominant but potentially non-negligible contribution to the error budget in forthcoming weak lensing surveys. We also investigate the dependence of parameter biases on the amplitude and angular scale of spatial variations of the m-bias field, and conclude that requirements should be placed on the rms of spatial variations of the m-bias, in addition to any requirement on the mean value. We find that, for a Euclid-like survey, biases generally exceed ~30% of the statistical error for m-bias rms 0.02 - 0.03 and can exceed the statistical error for rms ~0.04 - 0.05. This allows requirements to be set on the permissible amplitude of spatial variations of the m-bias that will arise due to systematics in forthcoming weak lensing measurements. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 12 pages, 8 figures