-
Image segmentation for analyzing galaxy-galaxy strong lensing systems
Authors:
Bryan Ostdiek,
Ana Diaz Rivero,
Cora Dvorkin
Abstract:
The goal of this paper is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binn…
▽ More
The goal of this paper is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binned mass range, or neither. Our network is only trained on images with a single smooth lens and either zero or one subhalo near the Einstein ring. On an independent test set with lenses with large ellipticities, quadrupole and octopole moments, and for source apparent magnitudes between 17-25, the area of the main lens is recovered accurately. On average, only 1.3% of the true area is missed and 1.2% of the true area is added to another part of the lens. In addition, subhalos as light as $10^{8.5}M_{\odot}$ can be detected if they lie in bright pixels along the Einstein ring. Furthermore, the model is able to generalize to new contexts it has not been trained on, such as locating multiple subhalos with varying masses or more than one large smooth lens.
△ Less
Submitted 26 January, 2022; v1 submitted 14 September, 2020;
originally announced September 2020.
-
Extracting the Subhalo Mass Function from Strong Lens Images with Image Segmentation
Authors:
Bryan Ostdiek,
Ana Diaz Rivero,
Cora Dvorkin
Abstract:
Detecting substructure within strongly lensed images is a promising route to shed light on the nature of dark matter. However, it is a challenging task, which traditionally requires detailed lens modeling and source reconstruction, taking weeks to analyze each system. We use machine-learning to circumvent the need for lens and source modeling and develop a neural network to both locate subhalos in…
▽ More
Detecting substructure within strongly lensed images is a promising route to shed light on the nature of dark matter. However, it is a challenging task, which traditionally requires detailed lens modeling and source reconstruction, taking weeks to analyze each system. We use machine-learning to circumvent the need for lens and source modeling and develop a neural network to both locate subhalos in an image as well as determine their mass using the technique of image segmentation. The network is trained on images with a single subhalo located near the Einstein ring across a wide range of apparent source magnitudes. The network is then able to resolve subhalos with masses $m\gtrsim 10^{8.5} M_{\odot}$. Training in this way allows the network to learn the gravitational lensing of light, and remarkably, it is then able to detect entire populations of substructure, even for locations further away from the Einstein ring than those used in training. Over a wide range of the apparent source magnitude, the false-positive rate is around three false subhalos per 100 images, coming mostly from the lightest detectable subhalo for that signal-to-noise ratio. With good accuracy and a low false-positive rate, counting the number of pixels assigned to each subhalo class over multiple images allows for a measurement of the subhalo mass function (SMF). When measured over three mass bins from $10^9M_{\odot}$--$10^{10} M_{\odot}$ the SMF slope is recovered with an error of 36% for 50 images, and this improves to 10% for 1000 images with Hubble Space Telescope-like noise.
△ Less
Submitted 14 February, 2022; v1 submitted 14 September, 2020;
originally announced September 2020.
-
Flow-Based Likelihoods for Non-Gaussian Inference
Authors:
Ana Diaz Rivero,
Cora Dvorkin
Abstract:
We investigate the use of data-driven likelihoods to bypass a key assumption made in many scientific analyses, which is that the true likelihood of the data is Gaussian. In particular, we suggest using the optimization targets of flow-based generative models, a class of models that can capture complex distributions by transforming a simple base distribution through layers of nonlinearities. We cal…
▽ More
We investigate the use of data-driven likelihoods to bypass a key assumption made in many scientific analyses, which is that the true likelihood of the data is Gaussian. In particular, we suggest using the optimization targets of flow-based generative models, a class of models that can capture complex distributions by transforming a simple base distribution through layers of nonlinearities. We call these flow-based likelihoods (FBL). We analyze the accuracy and precision of the reconstructed likelihoods on mock Gaussian data, and show that simply gauging the quality of samples drawn from the trained model is not a sufficient indicator that the true likelihood has been learned. We nevertheless demonstrate that the likelihood can be reconstructed to a precision equal to that of sampling error due to a finite sample size. We then apply FBLs to mock weak lensing convergence power spectra, a cosmological observable that is significantly non-Gaussian (NG). We find that the FBL captures the NG signatures in the data extremely well, while other commonly used data-driven likelihoods, such as Gaussian mixture models and independent component analysis, fail to do so. This suggests that works that have found small posterior shifts in NG data with data-driven likelihoods such as these could be underestimating the impact of non-Gaussianity in parameter constraints. By introducing a suite of tests that can capture different levels of NG in the data, we show that the success or failure of traditional data-driven likelihoods can be tied back to the structure of the NG in the data. Unlike other methods, the flexibility of the FBL makes it successful at tackling different types of NG simultaneously. Because of this, and consequently their likely applicability across datasets and domains, we encourage their use for inference when sufficient mock data are available for training.
△ Less
Submitted 6 November, 2020; v1 submitted 10 July, 2020;
originally announced July 2020.
-
Quantifying the Line-of-Sight Halo Contribution to the Dark Matter Convergence Power Spectrum from Strong Gravitational Lenses
Authors:
Atınç Çağan Şengül,
Arthur Tsang,
Ana Diaz Rivero,
Cora Dvorkin,
Hong-Ming Zhu,
Uroš Seljak
Abstract:
Galaxy-galaxy strong gravitational lenses have become a popular probe of dark matter (DM) by providing a window into structure formation on the smallest scales. In particular, the convergence power spectrum of subhalos within lensing galaxies has been suggested as a promising observable to study DM. However, the distances involved in strong-lensing systems are vast, and we expect the relevant volu…
▽ More
Galaxy-galaxy strong gravitational lenses have become a popular probe of dark matter (DM) by providing a window into structure formation on the smallest scales. In particular, the convergence power spectrum of subhalos within lensing galaxies has been suggested as a promising observable to study DM. However, the distances involved in strong-lensing systems are vast, and we expect the relevant volume to contain line-of-sight (LOS) halos that are not associated with the main lens. We develop a formalism to calculate the effect of LOS halos as an effective convergence power spectrum. The multi-lens plane equation couples the angular deflections of consecutive lens planes, but by assuming that the perturbations due to the LOS halos are small, we show that they can be projected onto the main-lens plane as effective subhalos. We test our formalism by simulating lensing systems using the full multi-plane lens equation and find excellent agreement. We show how the relative contribution of LOS halos and subhalos depends on the source and lens redshift, as well as the assumed halo and subhalo mass functions. For a fiducial system with fraction of DM halo mass in substructure $f_{\rm sub}=0.4\%$ for subhalo masses $[10^5-10^8]\rm{M}_{\odot}$, the interloper contribution to the power spectrum is at least several times greater than that of subhalos for source redshifts $z_s\gtrsim0.5$. Furthermore, it is likely that for the SLACS and BELLS lenses the interloper contribution dominates: $f_{\rm sub}\gtrsim2\%$ ($4\%$) is needed for subhalos to dominate in SLACS (BELLS), which is higher than current upper bounds on $f_{\rm sub}$ for our mass range. Since the halo mass function is better understood from first principles, the dominance of interlopers in galaxy-galaxy lenses with high-quality imaging can be seen as a significant advantage when translating this observable into a constraint on DM.
△ Less
Submitted 23 September, 2020; v1 submitted 12 June, 2020;
originally announced June 2020.
-
A Novel CMB Component Separation Method: Hierarchical Generalized Morphological Component Analysis
Authors:
Sebastian Wagner-Carena,
Max Hopkins,
Ana Diaz Rivero,
Cora Dvorkin
Abstract:
We present a novel technique for Cosmic Microwave Background (CMB) foreground subtraction based on the framework of blind source separation. Inspired by previous work incorporating local variation to Generalized Morphological Component Analysis (GMCA), we introduce Hierarchical GMCA (HGMCA), a Bayesian hierarchical graphical model for source separation. We test our method on $N_{\rm side}=256$ sim…
▽ More
We present a novel technique for Cosmic Microwave Background (CMB) foreground subtraction based on the framework of blind source separation. Inspired by previous work incorporating local variation to Generalized Morphological Component Analysis (GMCA), we introduce Hierarchical GMCA (HGMCA), a Bayesian hierarchical graphical model for source separation. We test our method on $N_{\rm side}=256$ simulated sky maps that include dust, synchrotron, free-free and anomalous microwave emission, and show that HGMCA reduces foreground contamination by $25\%$ over GMCA in both the regions included and excluded by the Planck UT78 mask, decreases the error in the measurement of the CMB temperature power spectrum to the $0.02-0.03\%$ level at $\ell>200$ (and $<0.26\%$ for all $\ell$), and reduces correlation to all the foregrounds. We find equivalent or improved performance when compared to state-of-the-art Internal Linear Combination (ILC)-type algorithms on these simulations, suggesting that HGMCA may be a competitive alternative to foreground separation techniques previously applied to observed CMB data. Additionally, we show that our performance does not suffer when we perturb model parameters or alter the CMB realization, which suggests that our algorithm generalizes well beyond our simplified simulations. Our results open a new avenue for constructing CMB maps through Bayesian hierarchical analysis.
△ Less
Submitted 26 April, 2020; v1 submitted 17 October, 2019;
originally announced October 2019.
-
Direct Detection of Dark Matter Substructure in Strong Lens Images with Convolutional Neural Networks
Authors:
Ana Diaz Rivero,
Cora Dvorkin
Abstract:
Strong gravitational lensing is a promising way of uncovering the nature of dark matter, by finding perturbations to images that cannot be well accounted for by modeling the lens galaxy without additional structure, be it subhalos (smaller halos within the smooth lens) or line-of-sight (LOS) halos. We present results attempting to infer the presence of substructure from images without requiring an…
▽ More
Strong gravitational lensing is a promising way of uncovering the nature of dark matter, by finding perturbations to images that cannot be well accounted for by modeling the lens galaxy without additional structure, be it subhalos (smaller halos within the smooth lens) or line-of-sight (LOS) halos. We present results attempting to infer the presence of substructure from images without requiring an intermediate step in which a smooth model has to be subtracted, using a simple convolutional neural network (CNN). We find that the network is only able to infer the presence of subhalos with $>75\%$ accuracy when they have masses of $\geq 5 \times 10^9$M$_{\odot}$ if they lie within the main lens galaxy. Since less massive foreground LOS halos can have the same effect as higher mass subhalos, the CNN can probe lower masses in the halo mass function. The accuracy does not improve significantly if we add a population of less massive subhalos. With the expectation of experiments such as HST and Euclid yielding thousands of high-quality strong lensing images in the next years, having a way of analyzing images quickly to identify candidates that merit further analysis to determine individual subhalo properties while preventing extensive resources being used for images that would yield null detections could be very useful. By understanding the sensitivity as a function of substructure mass, non-detections could be combined with the information from images with substructure to constrain the cold dark matter scenario, in particular if the sensitivity can be pushed to lower masses.
△ Less
Submitted 8 January, 2020; v1 submitted 30 September, 2019;
originally announced October 2019.
-
Dark Matter Science in the Era of LSST
Authors:
Keith Bechtol,
Alex Drlica-Wagner,
Kevork N. Abazajian,
Muntazir Abidi,
Susmita Adhikari,
Yacine Ali-Haïmoud,
James Annis,
Behzad Ansarinejad,
Robert Armstrong,
Jacobo Asorey,
Carlo Baccigalupi,
Arka Banerjee,
Nilanjan Banik,
Charles Bennett,
Florian Beutler,
Simeon Bird,
Simon Birrer,
Rahul Biswas,
Andrea Biviano,
Jonathan Blazek,
Kimberly K. Boddy,
Ana Bonaca,
Julian Borrill,
Sownak Bose,
Jo Bovy
, et al. (155 additional authors not shown)
Abstract:
Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We…
▽ More
Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Neutrino Mass from Cosmology: Probing Physics Beyond the Standard Model
Authors:
Cora Dvorkin,
Martina Gerbino,
David Alonso,
Nicholas Battaglia,
Simeon Bird,
Ana Diaz Rivero,
Andreu Font-Ribera,
George Fuller,
Massimiliano Lattanzi,
Marilena Loverde,
Julian B. Muñoz,
Blake Sherwin,
Anže Slosar,
Francisco Villaescusa-Navarro
Abstract:
Recent advances in cosmic observations have brought us to the verge of discovery of the absolute scale of neutrino masses. Nonzero neutrino masses are known evidence of new physics beyond the Standard Model. Our understanding of the clustering of matter in the presence of massive neutrinos has significantly improved over the past decade, yielding cosmological constraints that are tighter than any…
▽ More
Recent advances in cosmic observations have brought us to the verge of discovery of the absolute scale of neutrino masses. Nonzero neutrino masses are known evidence of new physics beyond the Standard Model. Our understanding of the clustering of matter in the presence of massive neutrinos has significantly improved over the past decade, yielding cosmological constraints that are tighter than any laboratory experiment, and which will improve significantly over the next decade, resulting in a guaranteed detection of the absolute neutrino mass scale.
△ Less
Submitted 8 March, 2019;
originally announced March 2019.
-
Observable Predictions for Massive-Neutrino Cosmologies with Model-Independent Dark Energy
Authors:
Ana Diaz Rivero,
V. Miranda,
Cora Dvorkin
Abstract:
We investigate the bounds on the sum of neutrino masses in a cosmic-acceleration scenario where the equation of state $w(z)$ of dark energy (DE) is constructed in a model-independent way, using a basis of principal components (PCs) that are allowed to cross the phantom barrier $w(z)=-1$. We find that the additional freedom provided to $w(z)$ means the DE can undo changes in the background expansio…
▽ More
We investigate the bounds on the sum of neutrino masses in a cosmic-acceleration scenario where the equation of state $w(z)$ of dark energy (DE) is constructed in a model-independent way, using a basis of principal components (PCs) that are allowed to cross the phantom barrier $w(z)=-1$. We find that the additional freedom provided to $w(z)$ means the DE can undo changes in the background expansion induced by massive neutrinos at low redshifts. This has two significant consequences: (1) it leads to a substantial increase in the upper bound for the sum of the neutrino masses ($M_ν < 0.33 - 0.55$ eV at the 95\% C.L. depending on the data sets and number of PCs included) compared to studies that choose a specific parametrization for $w(z)$; and (2) it causes $\sim1σ$ deviations from $Λ$CDM in the luminosity distance and the Hubble expansion rate at higher redshifts ($z \gtrsim 2$), where the contribution of DE is subdominant and there is little constraining data. The second point consequently means that there are also observable deviations in the shear power spectrum and in the matter power spectrum at low redshift, since the clustering of matter throughout cosmic time depends on the expansion rate. This provides a compelling case to pursue high-$z$ BAO and SN measurements as a way of disentangling the effects of neutrinos and dark energy. Finally, we find that the additional freedom given to the dark energy component has the effect of lowering $S_8$ with respect to $Λ$CDM.
△ Less
Submitted 7 March, 2019;
originally announced March 2019.
-
Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope
Authors:
Alex Drlica-Wagner,
Yao-Yuan Mao,
Susmita Adhikari,
Robert Armstrong,
Arka Banerjee,
Nilanjan Banik,
Keith Bechtol,
Simeon Bird,
Kimberly K. Boddy,
Ana Bonaca,
Jo Bovy,
Matthew R. Buckley,
Esra Bulbul,
Chihway Chang,
George Chapline,
Johann Cohen-Tanugi,
Alessandro Cuoco,
Francis-Yan Cyr-Racine,
William A. Dawson,
Ana Díaz Rivero,
Cora Dvorkin,
Denis Erkal,
Christopher D. Fassnacht,
Juan García-Bellido,
Maurizio Giannotti
, et al. (75 additional authors not shown)
Abstract:
Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d…
▽ More
Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .
△ Less
Submitted 24 April, 2019; v1 submitted 4 February, 2019;
originally announced February 2019.
-
Gravitational Lensing and the Power Spectrum of Dark Matter Substructure: Insights from the ETHOS N-body Simulations
Authors:
Ana Díaz Rivero,
Cora Dvorkin,
Francis-Yan Cyr-Racine,
Jesús Zavala,
Mark Vogelsberger
Abstract:
Strong gravitational lensing has been identified as a promising astrophysical probe to study the particle nature of dark matter. In this paper we present a detailed study of the power spectrum of the projected mass density (convergence) field of substructure in a Milky Way-sized halo. This power spectrum has been suggested as a key observable that can be extracted from strongly lensed images and y…
▽ More
Strong gravitational lensing has been identified as a promising astrophysical probe to study the particle nature of dark matter. In this paper we present a detailed study of the power spectrum of the projected mass density (convergence) field of substructure in a Milky Way-sized halo. This power spectrum has been suggested as a key observable that can be extracted from strongly lensed images and yield important clues about the matter distribution within the lens galaxy. We use two different $N$-body simulations from the ETHOS framework: one with cold dark matter and another with self-interacting dark matter and a cutoff in the initial power spectrum. Despite earlier works that identified $ k \gtrsim 100$ kpc$^{-1}$ as the most promising scales to learn about the particle nature of dark matter we find that even at lower wavenumbers - which are actually within reach of observations in the near future - we can gain important information about dark matter. Comparing the amplitude and slope of the power spectrum on scales $0.1 \lesssim k/$kpc$^{-1} \lesssim 10$ from lenses at different redshifts can help us distinguish between cold dark matter and other exotic dark matter scenarios that alter the abundance and central densities of subhalos. Furthermore, by considering the contribution of different mass bins to the power spectrum we find that subhalos in the mass range $10^7 - 10^8$ M$_{\odot}$ are on average the largest contributors to the power spectrum signal on scales $2 \lesssim k/$kpc$^{-1} \lesssim 15$, despite the numerous subhalos with masses $> 10^8$ M$_{\odot}$ in a typical lens galaxy. Finally, by comparing the power spectra obtained from the subhalo catalogs to those from the particle data in the simulation snapshots we find that the seemingly-too-simple halo model is in fact a fairly good approximation to the much more complex array of substructure in the lens.
△ Less
Submitted 31 August, 2018;
originally announced September 2018.
-
On the Power Spectrum of Dark Matter Substructure in Strong Gravitational Lenses
Authors:
Ana Diaz Rivero,
Francis-Yan Cyr-Racine,
Cora Dvorkin
Abstract:
Studying the smallest self-bound dark matter structure in our Universe can yield important clues about the fundamental particle nature of dark matter. Galaxy-scale strong gravitational lensing provides a unique way to detect and characterize dark matter substructures at cosmological distances from the Milky Way. Within the cold dark matter (CDM) paradigm, the number of low-mass subhalos within len…
▽ More
Studying the smallest self-bound dark matter structure in our Universe can yield important clues about the fundamental particle nature of dark matter. Galaxy-scale strong gravitational lensing provides a unique way to detect and characterize dark matter substructures at cosmological distances from the Milky Way. Within the cold dark matter (CDM) paradigm, the number of low-mass subhalos within lens galaxies is expected to be large, implying that their contribution to the lensing convergence field is approximately Gaussian and could thus be described by their power spectrum. We develop here a general formalism to compute from first principles the substructure convergence power spectrum for different populations of dark matter subhalos. As an example, we apply our framework to two distinct subhalo populations: a truncated Navarro-Frenk-White subhalo population motivated by standard CDM, and a truncated cored subhalo population motivated by self-interacting dark matter (SIDM). We study in detail how the subhalo abundance, mass function, internal density profile, and concentration affect the amplitude and shape of the substructure power spectrum. We determine that the power spectrum is mostly sensitive to a specific combination of the subhalo abundance and moments of the mass function, as well as to the average tidal truncation scale of the largest subhalos included in the analysis. Interestingly, we show that the asymptotic slope of the substructure power spectrum at large wave number reflects the internal density profile of the subhalos. In particular, the SIDM power spectrum exhibits a characteristic steepening at large wave number absent in the CDM power spectrum, opening the possibility of using this observable, if at all measurable, to discern between these two scenarios.
△ Less
Submitted 8 January, 2018; v1 submitted 14 July, 2017;
originally announced July 2017.
-
Probing the small-scale structure in strongly lensed systems via transdimensional inference
Authors:
Tansu Daylan,
Francis-Yan Cyr-Racine,
Ana Diaz Rivero,
Cora Dvorkin,
Douglas P. Finkbeiner
Abstract:
Strong lensing is a sensitive probe of the small-scale density fluctuations in the Universe. We implement a novel approach to modeling strongly lensed systems using probabilistic cataloging, which is a transdimensional, hierarchical, and Bayesian framework to sample from a metamodel (union of models with different dimensionality) consistent with observed photon count maps. Probabilistic cataloging…
▽ More
Strong lensing is a sensitive probe of the small-scale density fluctuations in the Universe. We implement a novel approach to modeling strongly lensed systems using probabilistic cataloging, which is a transdimensional, hierarchical, and Bayesian framework to sample from a metamodel (union of models with different dimensionality) consistent with observed photon count maps. Probabilistic cataloging allows us to robustly characterize modeling covariances within and across lens models with different numbers of subhalos. Unlike traditional cataloging of subhalos, it does not require model subhalos to improve the goodness of fit above the detection threshold. Instead, it allows the exploitation of all information contained in the photon count maps, for instance, when constraining the subhalo mass function. We further show that, by not including these small subhalos in the lens model, fixed-dimensional inference methods can significantly mismodel the data. Using a simulated Hubble Space Telescope (HST) dataset, we show that the subhalo mass function can be probed even when many subhalos in the sample catalogs are individually below the detection threshold and would be absent in a traditional catalog. With the planned Wide Field Infrared Space Telescope (WFIRST), simultaneous probabilistic cataloging of dark subhalos in high-resolution, deep strong lens images has the potential to constrain the subhalo mass function at even lower masses.
△ Less
Submitted 5 March, 2018; v1 submitted 19 June, 2017;
originally announced June 2017.