-
ELEPHANT: ExtragaLactic alErt Pipeline for Hostless AstroNomical Transients
Authors:
P. J. Pessi,
R. Durgesh,
L. Nakazono,
E. E. Hayes,
R. A. P. Oliveira,
E. E. O. Ishida,
A. Moitinho,
A. Krone-Martins,
B. Moews,
R. S. de Souza,
R. Beck,
M. A. Kuhn,
K. Nowak,
S. Vaughan
Abstract:
Context. Transient astronomical events that exhibit no discernible association with a host galaxy are commonly referred to as hostless. These rare phenomena are associated with extremely energetic events, and they can offer unique insights into the properties and evolution of stars and galaxies. However, the sheer number of transients captured by contemporary high-cadence astronomical surveys rend…
▽ More
Context. Transient astronomical events that exhibit no discernible association with a host galaxy are commonly referred to as hostless. These rare phenomena are associated with extremely energetic events, and they can offer unique insights into the properties and evolution of stars and galaxies. However, the sheer number of transients captured by contemporary high-cadence astronomical surveys renders the manual identification of all potential hostless transients impractical. Therefore, creating a systematic identification tool is crucial for studying these elusive events. Aims. We present the ExtragaLactic alErt Pipeline for Hostless AstroNomical Transients (ELEPHANT), a framework for filtering hostless transients in astronomical data streams. Methods. We used Fink to access all the ZTF alerts produced between January/2022 and December/2023, selecting only those associated with extragalactic transients. We then processed the associated stamps using a sequence of image analysis techniques to retrieve hostless candidates. Results. We find that less than 2% of all analyzed transients are potentially hostless. Among them, approximately 10% have a spectroscopic class reported on TNS, with Type Ia supernova being the most common class, followed by SLSN. Among the hostless candidates retrieved by our pipeline, there was SN 2018ibb, which has been proposed to be a PISN candidate; and SN 2022ann, one of only five known SNe Icn. When no class is reported on TNS, the dominant classes are QSO and SN candidates, the former obtained from SIMBAD and the latter inferred using the Fink ML classifier. Conclusions. ELEPHANT represents an effective strategy to filter extragalactic events within large and complex astronomical alert streams. There are many applications for which this pipeline will be useful, ranging from transient selection for follow-up to studies of transient environments.
△ Less
Submitted 28 April, 2024;
originally announced April 2024.
-
Physics-informed neural networks in the recreation of hydrodynamic simulations from dark matter
Authors:
Zhenyu Dai,
Ben Moews,
Ricardo Vilalta,
Romeel Dave
Abstract:
Physics-informed neural networks have emerged as a coherent framework for building predictive models that combine statistical patterns with domain knowledge. The underlying notion is to enrich the optimization loss function with known relationships to constrain the space of possible solutions. Hydrodynamic simulations are a core constituent of modern cosmology, while the required computations are…
▽ More
Physics-informed neural networks have emerged as a coherent framework for building predictive models that combine statistical patterns with domain knowledge. The underlying notion is to enrich the optimization loss function with known relationships to constrain the space of possible solutions. Hydrodynamic simulations are a core constituent of modern cosmology, while the required computations are both expensive and time-consuming. At the same time, the comparatively fast simulation of dark matter requires fewer resources, which has led to the emergence of machine learning algorithms for baryon inpainting as an active area of research; here, recreating the scatter found in hydrodynamic simulations is an ongoing challenge. This paper presents the first application of physics-informed neural networks to baryon inpainting by combining advances in neural network architectures with physical constraints, injecting theory on baryon conversion efficiency into the model loss function. We also introduce a punitive prediction comparison based on the Kullback-Leibler divergence, which enforces scatter reproduction. By simultaneously extracting the complete set of baryonic properties for the Simba suite of cosmological simulations, our results demonstrate improved accuracy of baryonic predictions based on dark matter halo properties, successful recovery of the fundamental metallicity relation, and retrieve scatter that traces the target simulation's distribution.
△ Less
Submitted 19 October, 2023; v1 submitted 24 March, 2023;
originally announced March 2023.
-
Photometric Redshift Uncertainties in Weak Gravitational Lensing Shear Analysis: Models and Marginalization
Authors:
Tianqing Zhang,
Markus Michael Rau,
Rachel Mandelbaum,
Xiangchong Li,
Ben Moews
Abstract:
Recovering credible cosmological parameter constraints in a weak lensing shear analysis requires an accurate model that can be used to marginalize over nuisance parameters describing potential sources of systematic uncertainty, such as the uncertainties on the sample redshift distribution $n(z)$. Due to the challenge of running Markov Chain Monte-Carlo (MCMC) in the high dimensional parameter spac…
▽ More
Recovering credible cosmological parameter constraints in a weak lensing shear analysis requires an accurate model that can be used to marginalize over nuisance parameters describing potential sources of systematic uncertainty, such as the uncertainties on the sample redshift distribution $n(z)$. Due to the challenge of running Markov Chain Monte-Carlo (MCMC) in the high dimensional parameter spaces in which the $n(z)$ uncertainties may be parameterized, it is common practice to simplify the $n(z)$ parameterization or combine MCMC chains that each have a fixed $n(z)$ resampled from the $n(z)$ uncertainties. In this work, we propose a statistically-principled Bayesian resampling approach for marginalizing over the $n(z)$ uncertainty using multiple MCMC chains. We self-consistently compare the new method to existing ones from the literature in the context of a forecasted cosmic shear analysis for the HSC three-year shape catalog, and find that these methods recover similar cosmological parameter constraints, implying that using the most computationally efficient of the approaches is appropriate. However, we find that for datasets with the constraining power of the full HSC survey dataset (and, by implication, those upcoming surveys with even tighter constraints), the choice of method for marginalizing over $n(z)$ uncertainty among the several methods from the literature may significantly impact the statistical uncertainties on cosmological parameters, and a careful model selection is needed to ensure credible parameter intervals.
△ Less
Submitted 18 April, 2023; v1 submitted 21 June, 2022;
originally announced June 2022.
-
Hybrid analytic and machine-learned baryonic property insertion into galactic dark matter haloes
Authors:
Ben Moews,
Romeel Davé,
Sourav Mitra,
Sultan Hassan,
Weiguang Cui
Abstract:
While cosmological dark matter-only simulations relying solely on gravitational effects are comparably fast to compute, baryonic properties in simulated galaxies require complex hydrodynamic simulations that are computationally costly to run. We explore the merging of an extended version of the equilibrium model, an analytic formalism describing the evolution of the stellar, gas, and metal content…
▽ More
While cosmological dark matter-only simulations relying solely on gravitational effects are comparably fast to compute, baryonic properties in simulated galaxies require complex hydrodynamic simulations that are computationally costly to run. We explore the merging of an extended version of the equilibrium model, an analytic formalism describing the evolution of the stellar, gas, and metal content of galaxies, into a machine learning framework. In doing so, we are able to recover more properties than the analytic formalism alone can provide, creating a high-speed hydrodynamic simulation emulator that populates galactic dark matter haloes in N-body simulations with baryonic properties. While there exists a trade-off between the reached accuracy and the speed advantage this approach offers, our results outperform an approach using only machine learning for a subset of baryonic properties. We demonstrate that this novel hybrid system enables the fast completion of dark matter-only information by mimicking the properties of a full hydrodynamic suite to a reasonable degree, and discuss the advantages and disadvantages of hybrid versus machine learning-only frameworks. In doing so, we offer an acceleration of commonly deployed simulations in cosmology.
△ Less
Submitted 15 November, 2022; v1 submitted 10 December, 2020;
originally announced December 2020.
-
Ridges in the Dark Energy Survey for cosmic trough identification
Authors:
Ben Moews,
Morgan A. Schmitz,
Andrew J. Lawler,
Joe Zuntz,
Alex I. Malz,
Rafael S. de Souza,
Ricardo Vilalta,
Alberto Krone-Martins,
Emille E. O. Ishida
Abstract:
Cosmic voids and their corresponding redshift-projected mass densities, known as troughs, play an important role in our attempt to model the large-scale structure of the Universe. Understanding these structures enables us to compare the standard model with alternative cosmologies, constrain the dark energy equation of state, and distinguish between different gravitational theories. In this paper,…
▽ More
Cosmic voids and their corresponding redshift-projected mass densities, known as troughs, play an important role in our attempt to model the large-scale structure of the Universe. Understanding these structures enables us to compare the standard model with alternative cosmologies, constrain the dark energy equation of state, and distinguish between different gravitational theories. In this paper, we extend the subspace-constrained mean shift algorithm, a recently introduced method to estimate density ridges, and apply it to 2D weak lensing mass density maps from the Dark Energy Survey Y1 data release to identify curvilinear filamentary structures. We compare the obtained ridges with previous approaches to extract trough structure in the same data, and apply curvelets as an alternative wavelet-based method to constrain densities. We then invoke the Wasserstein distance between noisy and noiseless simulations to validate the denoising capabilities of our method. Our results demonstrate the viability of ridge estimation as a precursor for denoising weak lensing observables to recover the large-scale structure, paving the way for a more versatile and effective search for troughs.
△ Less
Submitted 14 November, 2022; v1 submitted 18 May, 2020;
originally announced May 2020.
-
On the road to percent accuracy II: calibration of the non-linear matter power spectrum for arbitrary cosmologies
Authors:
Benjamin Giblin,
Matteo Cataneo,
Ben Moews,
Catherine Heymans
Abstract:
We introduce an emulator approach to predict the non-linear matter power spectrum for broad classes of beyond-$Λ$CDM cosmologies, using only a suite of $Λ$CDM $N$-body simulations. By including a range of suitably modified initial conditions in the simulations, and rescaling the resulting emulator predictions with analytical `halo model reactions', accurate non-linear matter power spectra for gene…
▽ More
We introduce an emulator approach to predict the non-linear matter power spectrum for broad classes of beyond-$Λ$CDM cosmologies, using only a suite of $Λ$CDM $N$-body simulations. By including a range of suitably modified initial conditions in the simulations, and rescaling the resulting emulator predictions with analytical `halo model reactions', accurate non-linear matter power spectra for general extensions to the standard $Λ$CDM model can be calculated. We optimise the emulator design by substituting the simulation suite with non-linear predictions from the standard {\sc halofit} tool. We review the performance of the emulator for artificially generated departures from the standard cosmology as well as for theoretically motivated models, such as $f (R)$ gravity and massive neutrinos. For the majority of cosmologies we have tested, the emulator can reproduce the matter power spectrum with errors $\lesssim 1\%$ deep into the highly non-linear regime. This work demonstrates that with a well-designed suite of $Λ$CDM simulations, extensions to the standard cosmological model can be tested in the non-linear regime without any reliance on expensive beyond-$Λ$CDM simulations.
△ Less
Submitted 13 September, 2019; v1 submitted 6 June, 2019;
originally announced June 2019.
-
Gaussbock: Fast parallel-iterative cosmological parameter estimation with Bayesian nonparametrics
Authors:
Ben Moews,
Joe Zuntz
Abstract:
We present and apply Gaussbock, a new embarrassingly parallel iterative algorithm for cosmological parameter estimation designed for an era of cheap parallel computing resources. Gaussbock uses Bayesian nonparametrics and truncated importance sampling to accurately draw samples from posterior distributions with an orders-of-magnitude speed-up in wall time over alternative methods. Contemporary pro…
▽ More
We present and apply Gaussbock, a new embarrassingly parallel iterative algorithm for cosmological parameter estimation designed for an era of cheap parallel computing resources. Gaussbock uses Bayesian nonparametrics and truncated importance sampling to accurately draw samples from posterior distributions with an orders-of-magnitude speed-up in wall time over alternative methods. Contemporary problems in this area often suffer from both increased computational costs due to high-dimensional parameter spaces and consequent excessive time requirements, as well as the need for fine tuning of proposal distributions or sampling parameters. Gaussbock is designed specifically with these issues in mind. We explore and validate the performance and convergence of the algorithm on a fast approximation to the Dark Energy Survey Year 1 (DES Y1) posterior, finding reasonable scaling behavior with the number of parameters. We then test on the full DES Y1 posterior using large-scale supercomputing facilities, and recover reasonable agreement with previous chains, although the algorithm can underestimate the tails of poorly-constrained parameters. Additionally, we discuss and demonstrate how Gaussbock recovers complex posterior shapes very well at lower dimensions, but faces challenges to perform well on such distributions in higher dimensions. In addition, we provide the community with a user-friendly software tool for accelerated cosmological parameter estimation based on the methodology described in this paper.
△ Less
Submitted 15 May, 2020; v1 submitted 23 May, 2019;
originally announced May 2019.
-
Photometry of high-redshift blended galaxies using deep learning
Authors:
Alexandre Boucaud,
Marc Huertas-Company,
Caroline Heneka,
Emille E. O. Ishida,
Nima Sedaghat,
Rafael S. de Souza,
Ben Moews,
Hervé Dole,
Marco Castellano,
Emiliano Merlin,
Valerio Roscani,
Andrea Tramacere,
Madhura Killedar,
Arlindo M. M. Trindade
Abstract:
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 slo…
▽ More
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.
△ Less
Submitted 3 May, 2019;
originally announced May 2019.
-
Stress testing the dark energy equation of state imprint on supernova data
Authors:
Ben Moews,
Rafael S. de Souza,
Emille E. O. Ishida,
Alex I. Malz,
Caroline Heneka,
Ricardo Vilalta,
Joe Zuntz
Abstract:
This work determines the degree to which a standard Lambda-CDM analysis based on type Ia supernovae can identify deviations from a cosmological constant in the form of a redshift-dependent dark energy equation of state w(z). We introduce and apply a novel random curve generator to simulate instances of w(z) from constraint families with increasing distinction from a cosmological constant. After pr…
▽ More
This work determines the degree to which a standard Lambda-CDM analysis based on type Ia supernovae can identify deviations from a cosmological constant in the form of a redshift-dependent dark energy equation of state w(z). We introduce and apply a novel random curve generator to simulate instances of w(z) from constraint families with increasing distinction from a cosmological constant. After producing a series of mock catalogs of binned type Ia supernovae corresponding to each w(z) curve, we perform a standard Lambda-CDM analysis to estimate the corresponding posterior densities of the absolute magnitude of type Ia supernovae, the present-day matter density, and the equation of state parameter. Using the Kullback-Leibler divergence between posterior densities as a difference measure, we demonstrate that a standard type Ia supernova cosmology analysis has limited sensitivity to extensive redshift dependencies of the dark energy equation of state. In addition, we report that larger redshift-dependent departures from a cosmological constant do not necessarily manifest easier-detectable incompatibilities with the Lambda-CDM model. Our results suggest that physics beyond the standard model may simply be hidden in plain sight.
△ Less
Submitted 5 July, 2019; v1 submitted 23 December, 2018;
originally announced December 2018.
-
Forging new worlds: high-resolution synthetic galaxies with chained generative adversarial networks
Authors:
Levi Fussell,
Ben Moews
Abstract:
Astronomy of the 21st century increasingly finds itself with extreme quantities of data. This growth in data is ripe for modern technologies such as deep image processing, which has the potential to allow astronomers to automatically identify, classify, segment and deblend various astronomical objects. In this paper, we explore the use of chained generative adversarial networks (GANs), a class of…
▽ More
Astronomy of the 21st century increasingly finds itself with extreme quantities of data. This growth in data is ripe for modern technologies such as deep image processing, which has the potential to allow astronomers to automatically identify, classify, segment and deblend various astronomical objects. In this paper, we explore the use of chained generative adversarial networks (GANs), a class of generative models that learn mappings from latent spaces to data distributions by modelling the joint distribution of the data, to produce physically realistic galaxy images as one use case of such models. In cosmology, such datasets can aid in the calibration of shape measurements for weak lensing by augmenting data with synthetic images. By measuring the distributions of multiple physical properties, we show that images generated with our approach closely follow the distributions of real galaxies, further establishing state-of-the-art GAN architectures as a valuable tool for modern-day astronomy.
△ Less
Submitted 16 March, 2019; v1 submitted 7 November, 2018;
originally announced November 2018.
-
Gaia DR2 unravels incompleteness of nearby cluster population: New open clusters in the direction of Perseus
Authors:
T. Cantat-Gaudin,
A. Krone-Martins,
N. Sedaghat,
A. Farahi,
R. S. de Souza,
R. Skalidis,
A. I. Malz,
S. Macêdo,
B. Moews,
C. Jordi,
A. Moitinho,
A. Castro-Ginard,
E. E. O. Ishida,
C. Heneka,
A. Boucaud,
A. M. M. Trindade
Abstract:
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…
▽ More
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.
△ Less
Submitted 21 March, 2019; v1 submitted 12 October, 2018;
originally announced October 2018.