The new generation of deep photometric surveys requires unprecedentedly precise shape and photometry measurements of billions of galaxies to achieve their main science goals. At such depths, one major limiting factor is the blending of galaxies due to line-of-sight projection, with an expected fraction of blended galaxies of up to 50%. Current deblending approaches are in most cases either too slow or not accurate enough to reach the level of requirements. This work explores the use of deep neural networks to estimate the photometry of blended pairs of galaxies in monochrome space images, similar to the ones that will be delivered by the Euclid space telescope. Using a clean sample of isolated galaxies from the CANDELS survey, we artificially blend them and train two different network models to recover the photometry of the two galaxies. We show that our approach can recover the original photometry of the galaxies before being blended with $\sim$7% accuracy without any human intervention and without any assumption on the galaxy shape. This represents an improvement of at least a factor of 4 compared to the classical SExtractor approach. We also show that forcing the network to simultaneously estimate a binary segmentation map results in a slightly improved photometry. All data products and codes will be made public to ease the comparison with other approaches on a common data set.

Open clusters (OCs) are popular tracers of the structure and evolutionary history of the Galactic disk. The OC population is often considered to be complete within 1.8 kpc of the Sun. The recent Gaia Data Release 2 (DR2) allows the latter claim to be challenged. We perform a systematic search for new OCs in the direction of Perseus using precise and accurate astrometry from Gaia DR2. We implement a coarse-to-fine search method. First, we exploit spatial proximity using a fast density-aware partitioning of the sky via a k-d tree in the spatial domain of Galactic coordinates, (l, b). Secondly, we employ a Gaussian mixture model in the proper motion space to quickly tag fields around OC candidates. Thirdly, we apply an unsupervised membership assignment method, UPMASK, to scrutinise the candidates. We visually inspect colour-magnitude diagrams to validate the detected objects. Finally, we perform a diagnostic to quantify the significance of each identified overdensity in proper motion and in parallax space We report the discovery of 41 new stellar clusters. This represents an increment of at least 20% of the previously known OC population in this volume of the Milky Way. We also report on the clear identification of NGC 886, an object previously considered an asterism. This letter challenges the previous claim of a near-complete sample of open clusters up to 1.8 kpc. Our results reveal that this claim requires revision, and a complete census of nearby open clusters is yet to be found.

Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogues. Here we present cosmoabc, a Python ABC sampler featuring a Population Monte Carlo (PMC) variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code is very flexible and can be easily coupled to an external simulator, while allowing to incorporate arbitrary distance and prior functions. As an example of practical application, we coupled cosmoabc with the numcosmo library and demonstrate how it can be used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function. cosmoabc is published under the GPLv3 license on PyPI and GitHub and documentation is available at http://goo.gl/SmB8EX

We explore the dynamics and evolution of the Universe at early and late times, focusing on both dark energy and extended gravity models and their astrophysical and cosmological consequences. Modified theories of gravity not only provide an alternative explanation for the recent expansion history of the universe, but they also offer a paradigm fundamentally distinct from the simplest dark energy models of cosmic acceleration. In this review, we perform a detailed theoretical and phenomenological analysis of different modified gravity models and investigate their consistency. We also consider the cosmological implications of well motivated physical models of the early universe with a particular emphasis on inflation and topological defects. Astrophysical and cosmological tests over a wide range of scales, from the solar system to the observable horizon, severely restrict the allowed models of the Universe. Here, we review several observational probes -- including gravitational lensing, galaxy clusters, cosmic microwave background temperature and polarization, supernova and baryon acoustic oscillations measurements -- and their relevance in constraining our cosmological description of the Universe.

The existence of multiple subclasses of type Ia supernovae (SNeIa) has been the subject of great debate in the last decade. One major challenge inevitably met when trying to infer the existence of one or more subclasses is the time consuming, and subjective, process of subclass definition. In this work, we show how machine learning tools facilitate identification of subtypes of SNeIa through the establishment of a hierarchical group structure in the continuous space of spectral diversity formed by these objects. Using Deep Learning, we were capable of performing such identification in a 4 dimensional feature space (+1 for time evolution), while the standard Principal Component Analysis barely achieves similar results using 15 principal components. This is evidence that the progenitor system and the explosion mechanism can be described by a small number of initial physical parameters. As a proof of concept, we show that our results are in close agreement with a previously suggested classification scheme and that our proposed method can grasp the main spectral features behind the definition of such subtypes. This allows the confirmation of the velocity of lines as a first order effect in the determination of SNIa subtypes, followed by 91bg-like events. Given the expected data deluge in the forthcoming years, our proposed approach is essential to allow a quick and statistically coherent identification of SNeIa subtypes (and outliers). All tools used in this work were made publicly available in the Python package Dimensionality Reduction And Clustering for Unsupervised Learning in Astronomy (DRACULA) and can be found within COINtoolbox (https://github.com/COINtoolbox/DRACULA).

Galaxy clusters are a valuable source of cosmological information. Their formation and evolution depends on the underlying cosmology and on the statistical nature of the primordial density fluctuations. In this work we investigate the impact of primordial non-gaussianities (PNG) on the scaling properties of galaxy clusters. We performed a series of cosmological hydrodynamic N-body simulations featuring adiabatic gas physics and different levels of non-Gaussian initial conditions within the $\Lambda$CDM framework. We focus on the T-M, S-M, Y-M and Yx-M scalings relating the total cluster mass with temperature, entropy and SZ cluster integrated pressure that reflect the thermodynamical state of the intra-cluster medium. Our results show that PNG have an impact on cluster scalings laws. The mass power-law indexes of the scalings are almost unaffected by the existence of PNG but the amplitude and redshift evolution of their normalizations are clearly affected. The effect is stronger for the evolution of the Y-M and Yx-M normalizations, which change by as much as 22% and 16% when $f_{NL}$ varies from -500 to 500, respectively. These results are consistent with the view that positive/negative $f_{NL}$ affect cluster profiles due to an increase/decrease of cluster concentrations. At low values of $f_{NL}$, as suggested by present Planck constraints on a scale invariant $f_{NL}$, the impact on the scalings normalizations is only a few percent, which is small when compared with the effect of additional gas physics and other cosmological effects such as dark energy. However if $f_{NL}$ is in fact a scale dependent parameter, PNG may have larger positive/negative amplitudes at clusters scales and therefore our results suggest that PNG should be taken into account when galaxy cluster data is used to infer cosmological parameters or to asses the constraining power of future cluster surveys.

The redshift dependence of the abundance of galaxy clusters is very sensitive to the statistical properties of primordial density perturbations. It can thus be used to probe small deviations from Gaussian initial conditions. Such deviations constitute a very important signature of many inflationary scenarios, and are thus expected to provide crucial information on physical processes which took place in the very early Universe. We have determined the biases which may be introduced in the estimation of cosmological parameters by wrongly assuming the absence of primordial non-Gaussianities. Although we find that the estimation of the present-day dark energy density using cluster counts is not very sensitive to the non-Gaussian properties of the density field, we show that the biases can be considerably larger in the estimation of the dark energy equation of state parameter $w$ and of the amplitude of the primordial density perturbations. Our results suggest that a significant level of non-Gaussianity at cluster scales may be able to reconcile the constraint on the amplitude of the primordial perturbations obtained using galaxy cluster number counts from the Planck Sunyaev-Zeldovich Catalog with that obtained from the primary Cosmic Microwave Background anisotropies measured by the Planck satellite.

The evolution with time of the abundance of galaxy clusters is very sensitive to the statistical properties of the primordial density perturbations. It can thus be used to probe small deviations from Gaussianity in the initial conditions. The characterization of such deviations would help distinguish between different inflationary scenarios, and provide us with information on physical processes which took place in the early Universe. We have found that when the information contained in the galaxy cluster counts is used to reconstruct the dark energy equation of state as a function of redshift, assuming erroneously that no primordial non-Gaussianities exist, an apparent evolution with time in the effective dark energy equation of state arises,characterized by the appearance of a clear discontinuity.

Current constraints on the dark energy equation of state parameter, $w$, are expected to be improved by more than one order of magnitude in the next decade. If $|w-1| \gsim 0.01$ around the present time, but the dark energy dynamics is sufficiently slow, it is possible that future constraints will rule out a cosmological constant while being consistent with a time-independent equation of state parameter. In this paper, we show that although models with such behavior can be constructed, they do require significant fine-tuning. Therefore, if the observed acceleration of the Universe is induced by a dark energy component, then finding $w \neq -1$ would, on its own, constitute very strong evidence for a dynamical dark energy equation of state.