Deep Learning (DL) has shown remarkable results in solving inverse problems in various domains. In particular, the Tikhonet approach is very powerful to deconvolve optical astronomical images (Sureau ...et al. 2020). Yet, this approach only uses the \(\ell_2\) loss, which does not guarantee the preservation of physical information (e.g. flux and shape) of the object reconstructed in the image. In Nammour et al. (2021), a new loss function was proposed in the framework of sparse deconvolution, which better preserves the shape of galaxies and reduces the pixel error. In this paper, we extend Tikhonet to take into account this shape constraint, and apply our new DL method, called ShapeNet, to optical and radio-interferometry simulated data set. The originality of the paper relies on i) the shape constraint we use in the neural network framework, ii) the application of deep learning to radio-interferometry image deconvolution for the first time, and iii) the generation of a simulated radio data set that we make available for the community. A range of examples illustrates the results.
We introduce a novel approach to reconstruct dark matter mass maps from weak gravitational lensing measurements. The cornerstone of the proposed method lies in a new modelling of the matter density ...field in the Universe as a mixture of two components:(1) a sparsity-based component that captures the non-Gaussian structure of the field, such as peaks or halos at different spatial scales; and (2) a Gaussian random field, which is known to well represent the linear characteristics of the field.Methods. We propose an algorithm called MCALens which jointly estimates these two components. MCAlens is based on an alternating minimization incorporating both sparse recovery and a proximal iterative Wiener filtering. Experimental results on simulated data show that the proposed method exhibits improved estimation accuracy compared to state-of-the-art mass map reconstruction methods.
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.
We present reconstructed convergence maps, \textit{mass maps}, from the Dark Energy Survey (DES) third year (Y3) weak gravitational lensing data set. The mass maps are weighted projections of the ...density field (primarily dark matter) in the foreground of the observed galaxies. We use four reconstruction methods, each is a \textit{maximum a posteriori} estimate with a different model for the prior probability of the map: Kaiser-Squires, null B-mode prior, Gaussian prior, and a sparsity prior. All methods are implemented on the celestial sphere to accommodate the large sky coverage of the DES Y3 data. We compare the methods using realistic \(\Lambda\)CDM simulations with mock data that are closely matched to the DES Y3 data. We quantify the performance of the methods at the map level and then apply the reconstruction methods to the DES Y3 data, performing tests for systematic error effects. The maps are compared with optical foreground cosmic-web structures and are used to evaluate the lensing signal from cosmic-void profiles. The recovered dark matter map covers the largest sky fraction of any galaxy weak lensing map to date.
Upcoming spectroscopic galaxy surveys are extremely promising to help in addressing the major challenges of cosmology, in particular in understanding the nature of the dark universe. The strength of ...these surveys comes from their unprecedented depth and width. Optimal extraction of their three-dimensional information is of utmost importance to best constrain the properties of the dark universe. Although there is theoretical motivation and novel tools to explore these surveys using the 3D spherical Fourier-Bessel (SFB) power spectrum of galaxy number counts \(C_\ell(k,k^\prime)\), most survey optimisations and forecasts are based on the tomographic spherical harmonics power spectrum \(C^{(ij)}_\ell\). We performed a new investigation of the information that can be extracted from the tomographic and 3D SFB techniques by comparing the forecast cosmological parameter constraints obtained from a Fisher analysis in the context of planned stage IV wide-field galaxy surveys. The comparison was made possible by careful and coherent treatment of non-linear scales in the two analyses. Nuisance parameters related to a scale- and redshift-dependent galaxy bias were also included for the first time in the computation of both the 3D SFB and tomographic power spectra. Tomographic and 3D SFB methods can recover similar constraints in the absence of systematics. However, constraints from the 3D SFB analysis are less sensitive to unavoidable systematics stemming from a redshift- and scale-dependent galaxy bias. Even for surveys that are optimised with tomography in mind, a 3D SFB analysis is more powerful. In addition, for survey optimisation, the figure of merit for the 3D SFB method increases more rapidly with redshift, especially at higher redshifts, suggesting that the 3D SFB method should be preferred for designing and analysing future wide-field spectroscopic surveys.
We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with ...realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of \(0.4\,\rm{deg}^2\) as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of \(22.5\,\rm{mag}\,\rm{arcsec}^{-2}\), and \(24.9\,\rm{mag}\,\rm{arcsec}^{-2}\) for the Euclid Deep Survey. This corresponds to approximately \(250\) million galaxies at the end of the mission and a \(50\,\%\) complete sample for stellar masses above \(10^{10.6}\,\rm{M}_\odot\) (resp. \(10^{9.6}\,\rm{M}_\odot\)) at a redshift \(z\sim0.5\) for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.
The primordial power spectrum describes the initial perturbations that seeded the large-scale structure we observe today. It provides an indirect probe of inflation or other structure-formation ...mechanisms. In this letter, we recover the primordial power spectrum from the Planck PR1 dataset, using our recently published algorithm PRISM. PRISM is a sparsity-based inversion method, that aims at recovering features in the primordial power spectrum from the empirical power spectrum of the cosmic microwave background (CMB). This ill-posed inverse problem is regularised using a sparsity prior on features in the primordial power spectrum in a wavelet dictionary. Although this non-parametric method does not assume a strong prior on the shape of the primordial power spectrum, it is able to recover both its general shape and localised features. As a results, this approach presents a reliable way of detecting deviations from the currently favoured scale-invariant spectrum. We applied PRISM to 100 simulated Planck data to investigate its performance on Planck-like data. We also tested the algorithm's ability to recover a small localised feature at \(k \sim 0.125\) Mpc\(^{-1}\), which caused a large dip at \(\ell \sim 1800\) in the angular power spectrum. We then applied PRISM to the Planck PR1 power spectrum to recover the primordial power spectrum. We find no significant departures from the fiducial Planck PR1 near scale-invariant primordial power spectrum with \(A_s=2.215\times10^{-9}\) and \(n_s = 0.9624\).
Detection of supernovae and, more generally, of transient events in large surveys can provide numerous false detections.In the case of a deferred processing of survey images, this implies ...reconstructing complete light curves for all detections, requiring sizable processing time and resources.Optimizing the detection of transient events is thus an important issue for both present and future surveys.We present here the optimization done in the SuperNova Legacy Survey (SNLS) for the 5-year data deferred photometric analysis. In this analysis, detections are derived from stacks of subtracted images with one stack per lunation.The 3-year analysis provided 300,000 detections dominated by signals of bright objects that were not perfectly subtracted.Allowing these artifacts to be detected leads not only to a waste of resources but also to possible signal coordinate contamination.We developed a subtracted image stack treatment to reduce the number of non SN-like events using morphological component analysis.This technique exploits the morphological diversity of objects to be detected to extract the signal of interest.At the level of our subtraction stacks, SN-like events are rather circular objects while most spurious detections exhibit different shapes.A two-step procedure was necessary to have a proper evaluation of the noise in the subtracted image stacks and thus a reliable signal extraction.We also set up a new detection strategy to obtain coordinates with good resolution for the extracted signal.SNIa MC generated images were used to study detection efficiency and coordinate resolution.When tested on SNLS 3 data this procedure decreases the number of detections by a factor of two, while losing only 10% of SN-like events, almost all faint.MC results show that SNIa detection efficiency is equivalent to that of the original method for bright events, while the coordinate resolution is improved.
The primordial power spectrum is an indirect probe of inflation or other structure-formation mechanisms. We introduce a new method, named \textbf{PRISM}, to estimate this spectrum from the empirical ...cosmic microwave background (CMB) power spectrum. This is a sparsity-based inversion method, which leverages a sparsity prior on features in the primordial spectrum in a wavelet dictionary to regularise the inverse problem. This non-parametric approach is able to reconstruct the global shape as well as localised features of spectrum accurately and proves to be robust for detecting deviations from the currently favoured scale-invariant spectrum. We investigate the strength of this method on a set of WMAP nine-year simulated data for three types of primordial spectra and then process the WMAP nine-year data as well as the Planck PR1 data. We find no significant departures from a near scale-invariant spectrum.