Search

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., Comment: 15 pages, 8 figures. Journal-accepted version

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

PAON4 is an L-band (1250-1500 MHz) small interferometer operating in transit mode deployed at the Nan\c{c}ay observatory in France, designed as a prototype instrument for Intensity Mapping. It features four 5~meter diameter dishes in a compact triangular configuration, with a total geometric collecting area of $\sim75 \mathrm{m^2}$, and equipped with dual polarization receivers. A total of 36 visibilities are computed from the 8 independent RF signals by the software correlator over the full 250~MHz RF band. The array operates in transit mode, with the dishes pointed toward a fixed declination, while the sky drifts across the instrument. Sky maps for each frequency channel are then reconstructed by combining the time-dependent visibilities from the different baselines observed at different declinations. This paper presents an overview of the PAON4 instrument design and goals, as a prototype for dish arrays to map the Large Scale Structure in radio, using intensity mapping of the atomic hydrogen $21~\mathrm{cm}$ line. We operated PAON4 over several years and use data from observations in different periods to assess the array performance. We present preliminary analysis of a large fraction of this data and discuss crucial issues for this type of instrument, such as the calibration strategy, instrument response stability, and noise behaviour., Comment: 17 pages, 16 figures, accepted to MNRAS

Apache Spark is a Big Data framework for working on large distributed datasets. Although widely used in the industry, it remains rather limited in the academic community or often restricted to software engineers. The goal of this paper is to show with practical uses-cases that the technology is mature enough to be used without excessive programming skills by astronomers or cosmologists in order to perform standard analyses over large datasets, as those originating from future galaxy surveys. To demonstrate it, we start from a realistic simulation corresponding to 10 years of LSST data taking (6 billions of galaxies). Then, we design, optimize and benchmark a set of Spark python algorithms in order to perform standard operations as adding photometric redshift errors, measuring the selection function or computing power spectra over tomographic bins. Most of the commands execute on the full 110 GB dataset within tens of seconds and can therefore be performed interactively in order to design full-scale cosmological analyses. A jupyter notebook summarizing the analysis is available at https://github.com/astrolabsoftware/1807.03078.

In the near future, cosmology will enter the wide and deep galaxy survey area allowing high-precision studies of the large scale structure of the universe in three dimensions. To test cosmological models and determine their parameters accurately, it is natural to confront data with exact theoretical expectations expressed in the observational parameter space (angles and redshift). The data-driven galaxy number count fluctuations on redshift shells, can be used to build correlation functions $C(\theta; z_1, z_2)$ on and between shells which can probe the baryonic acoustic oscillations, the distance-redshift distortions as well as gravitational lensing and other relativistic effects. Transforming the model to the data space usually requires the computation of the angular power spectrum $C_\ell(z_1, z_2)$ but this appears as an artificial and inefficient step plagued by apodization issues. In this article we show that it is not necessary and present a compact expression for $C(\theta; z_1, z_2)$ that includes directly the leading density and redshift space distortions terms from the full linear theory. It can be evaluated using a fast integration method based on Clenshaw-Curtis quadrature and Chebyshev polynomial series. This new method to compute the correlation functions without any Limber approximation, allows us to produce and discuss maps of the correlation function directly in the observable space and is a significant step towards disentangling the data from the tested models.

The statistical distribution of galaxies is a powerful probe to constrain cosmological models and gravity. In particular the matter power spectrum $P(k)$ brings information about the cosmological distance evolution and the galaxy clustering together. However the building of $P(k)$ from galaxy catalogues needs a cosmological model to convert angles on the sky and redshifts into distances, which leads to difficulties when comparing data with predicted $P(k)$ from other cosmological models, and for photometric surveys like LSST. The angular power spectrum $C_\ell(z_1,z_2)$ between two bins located at redshift $z_1$ and $z_2$ contains the same information than the matter power spectrum, is free from any cosmological assumption, but the prediction of $C_\ell(z_1,z_2)$ from $P(k)$ is a costly computation when performed exactly. The Angpow software aims at computing quickly and accurately the auto ($z_1=z_2$) and cross ($z_1 \neq z_2$) angular power spectra between redshift bins. We describe the developed algorithm, based on developments on the Chebyshev polynomial basis and on the Clenshaw-Curtis quadrature method. We validate the results with other codes, and benchmark the performance. Angpow is flexible and can handle any user defined power spectra, transfer functions, and redshift selection windows. The code is fast enough to be embedded inside programs exploring large cosmological parameter spaces through the $C_\ell(z_1,z_2)$ comparison with data. We emphasize that the Limber's approximation, often used to fasten the computation, gives wrong $C_\ell$ values for cross-correlations., Comment: Published in Astronomy & Astrophysics

We have observed regions of three galaxy clusters at z$\sim$ [0.06, 0.09] (Abell85, Abell1205, Abell2440), as well as calibration sources with the Nancay radiotelescope (NRT) to search for 21 cm emission and fully characterize the FPGA based BAORadio digital backend. The total observation time of few hours per source have been distributed over few months, from March 2011 to January 2012, due to scheduling constraints of the NRT, which is a transit telescope. Data have been acquired in parallel with the NRT standard correlator (ACRT) back-end, as well as with the BAORadio data acquisition system. The latter enables wide band instantaneous observation of the [1250, 1500]MHz frequency range, as well as the use of powerful RFI mitigation methods thanks to its fine time sampling. A number of questions related to instrument stability, data processing and calibration are discussed. We have obtained the radiometer curves over the integration time range [0.01,10 000] seconds and we show that sensitivities of few mJy over most of the wide frequency band can be reached with the NRT. It is clearly shown that in blind line search, which is the context of HI intensity mapping for Baryon Acoustic Oscillations, the new acquisition system and processing pipeline outperforms the standard one. We report a positive detection of 21 cm emission at 3 sigma-level from galaxies in the outer region of Abell85 at 1352 MHz (14 400 km/s) corresponding to a line strength of 0.8 Jy km/s. We observe also an excess power around 1318 MHz (21 600 km/s), although at lower statistical significance, compatible with emission from Abell1205 galaxies. Detected radio line emissions have been cross matched with optical catalogs and we have derived hydrogen mass estimates., Comment: 18 pages, 17 figures

Large Scale Structures (LSS) in the universe can be traced using the neutral atomic hydrogen HI through its 21cm emission. Such a 3D matter distribution map can be used to test the Cosmological model and to constrain the Dark Energy properties or its equation of state. A novel approach, called intensity mapping can be used to map the HI distribution, using radio interferometers with large instantaneous field of view and waveband. In this paper, we study the sensitivity of different radio interferometer configurations, or multi-beam instruments for the observation of large scale structures and BAO oscillations in 21cm and we discuss the problem of foreground removal. For each configuration, we determine instrument response by computing the (u,v) or Fourier angular frequency plane coverage using visibilities. The (u,v) plane response is the noise power spectrum, hence the instrument sensitivity for LSS P(k) measurement. We describe also a simple foreground subtraction method to separate LSS 21 cm signal from the foreground due to the galactic synchrotron and radio sources emission. We have computed the noise power spectrum for different instrument configuration as well as the extracted LSS power spectrum, after separation of 21cm-LSS signal from the foregrounds. We have also obtained the uncertainties on the Dark Energy parameters for an optimized 21 cm BAO survey. We show that a radio instrument with few hundred simultaneous beams and a collecting area of ~10000 m^2 will be able to detect BAO signal at redshift z ~ 1 and will be competitive with optical surveys., Comment: 16 pages, 18 figures submitted to Astronomy & Astrophysics (A&A)

3D mapping of matter distribution in the universe through the 21 cm radio emission of atomic hydrogen is a complementary approach to optical surveys for the study of the Large Scale Structures, in particular for measuring the BAO (Baryon Acoustic Oscillation) scale up to redshifts z <~ 3 and constrain dark energy. We propose to carry such a survey through a novel method, called intensity mapping, without detecting individual galaxies radio emission. This method requires a wide band instrument, 100 MHz or larger, and multiple beams, while a rather modest angular resolution of 10 arcmin would be sufficient. The instrument would have a few thousand square meters of collecting area and few hundreds of simultaneous beams. These constraints could be fulfilled with a dense array of receivers in interferometric mode, or a phased array at the focal plane of a large antenna.

PARISROC is a complete read out chip, in AMS SiGe 0.35 micron technology [1], for photomultipliers array. It allows triggerless acquisition for next generation neutrino experiments and it belongs to an R&D program funded by the French national agency for research (ANR) called PMm2: Innovative electronics for photodetectors array used in High Energy Physics and Astroparticles [2] (ref.ANR-06-BLAN- 0186). The ASIC integrates 16 independent and auto triggered channels with variable gain and provides charge and time measurement by a Wilkinson ADC and a 24-bit Counter. The charge measurement should be performed from 1 up to 300 photo-electrons (p.e.) with a good linearity. The time measurement allowed to a coarse time with a 24-bit counter at 10 MHz and a fine time on a 100ns ramp to achieve a resolution of 1 ns. The ASIC sends out only the relevant data through network cables to the central data storage. This paper describes the front-end electronics ASIC called PARISROC., Comment: Topical Workshop on Electronics for Particle Physics (TWEPP-09, Paris, France, Sept. 2009)

The feasibility of a next generation neutrino observatory in Europe is being considered within the LAGUNA design study. To accommodate giant neutrino detectors and shield them from cosmic rays, a new very large underground infrastructure is required. Seven potential candidate sites in different parts of Europe and at several distances from CERN are being studied: Boulby (UK), Canfranc (Spain), Fr\'ejus (France/Italy), Pyh\"asalmi (Finland), Polkowice-Sieroszowice (Poland), Slanic (Romania) and Umbria (Italy). The design study aims at the comprehensive and coordinated technical assessment of each site, at a coherent cost estimation, and at a prioritization of the sites within the summer 2010., Comment: 5 pages, contribution to the Workshop "European Strategy for Future Neutrino Physics", CERN, Oct. 2009

PARISROC is a complete read out chip, in AMS SiGe 0.35 !m technology, for photomultipliers array. It allows triggerless acquisition for next generation neutrino experiments and it belongs to an R&D program funded by the French national agency for research (ANR) called PMm2: ?Innovative electronics for photodetectors array used in High Energy Physics and Astroparticles? (ref.ANR-06-BLAN-0186). The ASIC (Application Specific Integrated Circuit) integrates 16 independent and auto triggered channels with variable gain and provides charge and time measurement by a Wilkinson ADC (Analog to Digital Converter) and a 24-bit Counter. The charge measurement should be performed from 1 up to 300 photo- electrons (p.e.) with a good linearity. The time measurement allowed to a coarse time with a 24-bit counter at 10 MHz and a fine time on a 100ns ramp to achieve a resolution of 1 ns. The ASIC sends out only the relevant data through network cables to the central data storage. This paper describes the front-end electronics ASIC called PARISROC., Comment: IEEE Nuclear Science Symposium an Medical Imaging Conference (2009 NSS/MIC)

The next generation of proton decay and neutrino experiments, the post-SuperKamiokande detectors as those that will take place in megaton size water tanks, will require very large surfaces of photodetection and a large volume of data. Even with large hemispherical photomultiplier tubes, the expected number of channels should reach hundreds of thousands. A funded R&D program to implement a solution is presented here. The very large surface of photodetection is segmented in macro pixels made of 16 hemispherical (12 inches) photomultiplier tubes connected to an autonomous front-end which works on a triggerless data acquisition mode. The expected data transmission rate is 5 Mb/s per cable, which can be achieved with existing techniques. This architecture allows to reduce considerably the cost and facilitate the industrialization. This document presents the simulations and measurements which define the requirements for the photomultipliers and the electronics. A proto-type of front-end electronics was successfully tested with 16 photomultiplier tubes supplied by a single high voltage, validating the built-in gain adjustment and the calibration principle. The first tests and calculations on the photomultiplier glass led to the study of a new package optimized for a 10 bar pressure in order to sustain the high underwater pressure., Comment: 1 pdf file, 4 pages, 4 figures, NDIP08, submitted to Nucl. Instr. and Meth. Phys. Res. A

This report summarises the conclusions from the detector group of the International Scoping Study of a future Neutrino Factory and Super-Beam neutrino facility. The baseline detector options for each possible neutrino beam are defined as follows: 1. A very massive (Megaton) water Cherenkov detector is the baseline option for a sub-GeV Beta Beam and Super Beam facility. 2. There are a number of possibilities for either a Beta Beam or Super Beam (SB) medium energy facility between 1-5 GeV. These include a totally active scintillating detector (TASD), a liquid argon TPC or a water Cherenkov detector. 3. A 100 kton magnetized iron neutrino detector (MIND) is the baseline to detect the wrong sign muon final states (golden channel) at a high energy (20-50 GeV) neutrino factory from muon decay. A 10 kton hybrid neutrino magnetic emulsion cloud chamber detector for wrong sign tau detection (silver channel) is a possible complement to MIND, if one needs to resolve degeneracies that appear in the $\delta$-$\theta_{13}$ parameter space., Comment: Detector report of the International Scoping Study of a future Neutrino Factory and Super-Beam facility, 86 pages, 49 figures

This document reports on a series of experimental and theoretical studies conducted to assess the astro-particle physics potential of three future large-scale particle detectors proposed in Europe as next generation underground observatories. The proposed apparatus employ three different and, to some extent, complementary detection techniques: GLACIER (liquid Argon TPC), LENA (liquid scintillator) and MEMPHYS (\WC), based on the use of large mass of liquids as active detection media. The results of these studies are presented along with a critical discussion of the performance attainable by the three proposed approaches coupled to existing or planned underground laboratories, in relation to open and outstanding physics issues such as the search for matter instability, the detection of astrophysical- and geo-neutrinos and to the possible use of these detectors in future high-intensity neutrino beams., Comment: 50 pages, 26 figures

A water \v{C}erenkov detector project, of megaton scale, to be installed in the Fr\'ejus underground site and dedicated to nucleon decay, neutrinos from supernovae, solar and atmospheric neutrinos, as well as neutrinos from a super-beam and/or a beta-beam coming from CERN, is presented and compared with competitor projects in Japan and in the USA. The performances of the European project are discussed, including the possibility to measure the mixing angle $\theta_{13}$ and the CP-violating phase $\delta$., Comment: 1+33 pages, 14 figures, Expression of Interest of MEMPHYS project

We consider the physics potential of CERN based neutrino oscillation experiments consisting of a Beta Beam (BB) and a Super Beam (SPL) sending neutrinos to MEMPHYS, a 440 kt water \v{C}erenkov detector at Frejus, at a distance of 130 km from CERN. The $\theta_{13}$ discovery reach and the sensitivity to CP violation are investigated, including a detailed discussion of parameter degeneracies and systematical errors. For SPL sensitivities similar to the ones of the phase II of the T2K experiment (T2HK) are obtained, whereas the BB may reach significantly better sensitivities, depending on the achieved number of total ion decays. The results for the CERN-MEMPHYS experiments are less affected by systematical uncertainties than T2HK. We point out that by a combination of data from BB and SPL a measurement with antineutrinos is not necessary and hence the same physics results can be obtained within about half of the measurement time compared to one single experiment. Furthermore, it is shown how including data from atmospheric neutrinos in the MEMPHYS detector allows to resolve parameter degeneracies and, in particular, provides sensitivity to the neutrino mass hierarchy and the octant of $\theta_{23}$., Comment: 32 pages, 17 figures, minor improvements on the text wrt to v2, version to appear in JHEP

An optimization of the CERN-SPL beam line has been performed which leads to better sensitivities to the $\theta\_{13}$ mixing angle and to the $\delta\_{CP}$ violating phase than those advocated considering baseline scenario., Comment: to be published in Nuclear Physics B (Proceedings Supplement)

The poster presented the OPERA experiment and the SPL-Fr\'{e}jus neutrino super beam project. OPERA (1) is an experiment willing to see tau neutrino appearance. The detector is described in section 2 and the expected results are given. Section 3 is devoted to the SPL neutrino super beam (2), which search for $\theta\_{13}$. An optimisation of the SPL energy is proposed.