Source author record

Tom Charnock

Tom Charnock appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.IM hep-th Machine Learning

Catalog footprint

What is connected

7works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Catalog-free modeling of galaxy types in deep images: Massive dimensional reduction with neural networks

Current models of galaxy evolution are constrained by the analysis of catalogs containing the flux and size of galaxies extracted from multiband deep fields carrying inevitable observational and extraction-related biases which can be highly correlated. In practice, taking all of these effects simultaneously into account is difficult, and derived models are inevitably biased. To address this issue, we use robust likelihood-free methods for the inference of luminosity function parameters, made possible via massive compression of multiband images using artificial neural networks. This technique makes the use of catalogs unnecessary when comparing observed and simulated multiband deep fields and constraining model parameters. A forward modeling approach generates galaxies of multiple types depending on luminosity function parameters and paints them on photometric multiband deep fields including both the instrumental and observational characteristics. The simulated and the observed images present the same selection effects and can therefore be properly compared. We train a fully-convolutional neural network to extract the most model-parameter-sensitive summary statistics out of these realistic simulations, shrinking down the dimensionality of the summary space. Finally, using the trained network to compress both observed and simulated deep fields, the model parameter values are constrained through Population Monte Carlo likelihood-free inference. Using synthetic photometric multiband deep fields similar to the CFHTLS and D1/D2 deep fields and massively compressing them through the convolutional neural network, we demonstrate the robustness, accuracy and consistency of this new catalog-free inference method. We are able to constrain the parameters of luminosity functions of different types of galaxies and our results are fully compatible with the classic catalog extraction approaches.

preprint2022arXiv

The Cosmic Graph: Optimal Information Extraction from Large-Scale Structure using Catalogues

We present an implicit likelihood approach to quantifying cosmological information over discrete catalogue data, assembled as graphs. To do so, we explore cosmological parameter constraints using mock dark matter halo catalogues. We employ Information Maximising Neural Networks (IMNNs) to quantify Fisher information extraction as a function of graph representation. We a) demonstrate the high sensitivity of modular graph structure to the underlying cosmology in the noise-free limit, b) show that graph neural network summaries automatically combine mass and clustering information through comparisons to traditional statistics, c) demonstrate that networks can still extract information when catalogues are subject to noisy survey cuts, and d) illustrate how nonlinear IMNN summaries can be used as asymptotically optimal compressed statistics for Bayesian simulation-based inference. We reduce the area of joint $Ω_m, σ_8$ parameter constraints with small ($\sim$100 object) halo catalogues by a factor of 42 over the two-point correlation function, and demonstrate that the networks automatically combine mass and clustering information. This work utilises a new IMNN implementation over graph data in Jax, which can take advantage of either numerical or auto-differentiability. We also show that graph IMNNs successfully compress simulations away from the fiducial model at which the network is fitted, indicating a promising alternative to n-point statistics in catalogue simulation-based analyses.

preprint2020arXiv

Detecting outliers in astronomical images with deep generative networks

With the advent of future big-data surveys, automated tools for unsupervised discovery are becoming ever more necessary. In this work, we explore the ability of deep generative networks for detecting outliers in astronomical imaging datasets. The main advantage of such generative models is that they are able to learn complex representations directly from the pixel space. Therefore, these methods enable us to look for subtle morphological deviations which are typically missed by more traditional moment-based approaches. We use a generative model to learn a representation of expected data defined by the training set and then look for deviations from the learned representation by looking for the best reconstruction of a given object. In this first proof-of-concept work, we apply our method to two different test cases. We first show that from a set of simulated galaxies, we are able to detect $\sim90\%$ of merging galaxies if we train our network only with a sample of isolated ones. We then explore how the presented approach can be used to compare observations and hydrodynamic simulations by identifying observed galaxies not well represented in the models.

preprint2020arXiv

Super-resolution emulator of cosmological simulations using deep physical models

We present an extension of our recently developed Wasserstein optimized model to emulate accurate high-resolution features from computationally cheaper low-resolution cosmological simulations. Our deep physical modelling technique relies on restricted neural networks to perform a mapping of the distribution of the low-resolution cosmic density field to the space of the high-resolution small-scale structures. We constrain our network using a single triplet of high-resolution initial conditions and the corresponding low- and high-resolution evolved dark matter simulations from the Quijote suite of simulations. We exploit the information content of the high-resolution initial conditions as a well constructed prior distribution from which the network emulates the small-scale structures. Once fitted, our physical model yields emulated high-resolution simulations at low computational cost, while also providing some insights about how the large-scale modes affect the small-scale structure in real space.

preprint2019arXiv

Neural physical engines for inferring the halo mass distribution function

An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques.

preprint2016arXiv

CMB constraints on cosmic strings and superstrings

We present the first complete Markov chain Monte Carlo analysis of cosmological models with evolving cosmic (super)string networks, using the unconnected segment model in the unequal-time correlator formalism. For ordinary cosmic string networks, we derive joint constraints on Lambda cold dark matter (CDM) and string network parameters, namely the string tension Gmu, the loop-chopping efficiency c_r and the string wiggliness α. For cosmic superstrings, we obtain joint constraints on the fundamental string tension Gmu_F, the string coupling g_s, the self-interaction coefficient c_s, and the volume of compact extra dimensions w. This constitutes the most comprehensive CMB analysis of LambdaCDM cosmology + strings to date. For ordinary cosmic string networks our updated constraint on the string tension is, in relativistic units, Gmu<1.1x10^-7, while for cosmic superstrings our constraint on the fundamental string tension is Gmu_F<2.8x10^-8, both obtained using Planck2015 temperature and polarisation data.

preprint2014arXiv

Tension between the power spectrum of density perturbations measured on large and small scales

There is a tension between measurements of the amplitude of the power spectrum of density perturbations inferred using the Cosmic Microwave Background (CMB) and directly measured by Large-Scale Structure (LSS) on smaller scales. We show that this tension exists, and is robust, for a range of LSS indicators including clusters, lensing and redshift space distortions and using CMB data from either $Planck$ or WMAP+SPT/ACT. One obvious way to try to reconcile this is the inclusion of a massive neutrino which could be either active or sterile. Using $Planck$ and a combination of all the LSS data we find that (i) for an active neutrino $\sum m_ν= (0.357\pm0.099)\,{\rm eV}$ and (ii) for a sterile neutrino $m_{\rm sterile}^{\rm eff}= (0.67\pm0.18)\,{\rm eV}$ and $ΔN_{\rm eff}= 0.32\pm0.20$. This is, however, at the expense of a degraded fit to $Planck$ temperature data, and we quantify the residual tension at $2.5σ$ and $1.6 σ$ for massive and sterile neutrinos respectively. We also consider alternative explanations including a lower redshift for reionization that would be in conflict with polarisation measurements made by WMAP and $ad$-$hoc$ modifications to primordial power spectrum.