Showing 1–50 of 276 results for author: Melchior, P

The little coadd that could: Estimating shear from coadded images

Authors: Robert Armstrong, Erin Sheldon, Eric Huff, Jim Bosch, Eli Rykoff, Rachel Mandelbaum, Arun Kannawadi, Peter Melchior, Robert Lupton, Matthew R. Becker, Yusra Al-Sayyed, The LSST Dark Energy Science Collaboration

Abstract: Upcoming wide field surveys will have many overlapping epochs of the same region of sky. The conventional wisdom is that in order to reduce the errors sufficiently for systematics-limited measurements, like weak lensing, we must do simultaneous fitting of all the epochs. Using current algorithms this will require a significant amount of computing time and effort. In this paper, we revisit the pote… ▽ More Upcoming wide field surveys will have many overlapping epochs of the same region of sky. The conventional wisdom is that in order to reduce the errors sufficiently for systematics-limited measurements, like weak lensing, we must do simultaneous fitting of all the epochs. Using current algorithms this will require a significant amount of computing time and effort. In this paper, we revisit the potential of using coadds for shear measurements. We show on a set of image simulations that the multiplicative shear bias can be constrained below the 0.1% level on coadds, which is sufficient for future lensing surveys. We see no significant differences between simultaneous fitting and coadded approaches for two independent shear codes: Metacalibration and BFD. One caveat of our approach is the assumption of a principled coadd, i.e. the PSF is mathematically well-defined for all the input images. This requires us to reject CCD images that do not fully cover the coadd region. We estimate that the number of epochs that must be rejected for a survey like LSST is on the order of 20%, resulting in a small loss in depth of less than 0.1 magnitudes. We also put forward a cell-based coaddition scheme that meets the above requirements for unbiased weak lensing shear estimation in the context of LSST. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Constraining Protoplanetary Disk Winds from Forbidden Line Profiles with Simulation-based Inference

Authors: Ahmad Nemer, ChangHoon Hahn, Jiaxuan Li, Peter Melchior, Jeremy Goodman

Abstract: Protoplanetary disks are the sites of vigorous hydrodynamic processes, such as accretion and outflows, and ultimately establish the conditions for the formation of planets. The properties of disk outflows are often inferred through analysis of forbidden emission lines. These lines contain multiple overlapping components, tracing different emission regions with different processes that excite them:… ▽ More Protoplanetary disks are the sites of vigorous hydrodynamic processes, such as accretion and outflows, and ultimately establish the conditions for the formation of planets. The properties of disk outflows are often inferred through analysis of forbidden emission lines. These lines contain multiple overlapping components, tracing different emission regions with different processes that excite them: a high-velocity component (tracing a jet), a broad low-velocity component (tracing inner disk wind), and a narrow low-velocity component (tracing outer disk wind). They are also heavily contaminated by background spectral features. All of these challenges call into question the traditional approach of fitting Gaussian components to the line profiles, and cloud the physical interpretation of those components. We introduce a novel statistical technique to analyze emission lines in protoplanetary disks. Simulation-Based Inference is a computationally efficient machine learning technique that produces posterior distributions of the parameters (e.g. magnetic field, radiation sources, geometry) of a representative wind model when given a spectrum, without any prior assumption about line shapes (e.g.symmetry). In this pathfinder study, we demonstrate that this technique indeed accurately recovers the parameters from simulated spectra without noise and background. A following work will deal with the analysis of observed spectra. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Score-matching neural networks for improved multi-band source separation

Authors: Matt L. Sampson, Peter Melchior, Charlotte Ward, Sufia Birmingham

Abstract: We present the implementation of a score-matching neural network that represents a data-driven prior for non-parametric galaxy morphologies. The gradients of this prior can be included in the optimization routine of the recently developed multi-band modeling framework Scarlet2, a redesign of the Scarlet method currently employed as deblender in the pipelines of the HyperSuprimeCam survey and the R… ▽ More We present the implementation of a score-matching neural network that represents a data-driven prior for non-parametric galaxy morphologies. The gradients of this prior can be included in the optimization routine of the recently developed multi-band modeling framework Scarlet2, a redesign of the Scarlet method currently employed as deblender in the pipelines of the HyperSuprimeCam survey and the Rubin Observatory. The addition of the prior avoids the requirement of nondifferentiable constraints, which can lead to convergence failures we discovered in Scarlet. We present the architecture and training details of our score-matching neural network and show with simulated Rubin-like observations that Scarlet2 outperforms Scarlet in accuracy of total flux and morphology estimates, while maintaining excellent performance for colors. We also demonstrate significant improvements in the robustness to inaccurate initializations. Scarlet2 is written in python, extendend by JAX and equinox, and is fully GPU compatible. The implementation and data package of the score model are publicly available at https://github.com/pmelchior/scarlet2. △ Less

Submitted 14 January, 2024; originally announced January 2024.

Comments: 17 pages, 12 figures, submitted to Astronomy and Computing

Joint Cosmic Density Reconstruction from Photometric and Spectroscopic Samples

Authors: Benjamin Horowitz, Peter Melchior

Abstract: We reconstruct the dark matter density field from spatially overlapping spectroscopic and photometric redshift catalogs through a forward modelling approach. Instead of directly inferring the underlying density field, we find the best fitting initial Gaussian fluctuations that will evolve into the observed cosmic volume. To account for the substantial uncertainty of photometric redshifts we employ… ▽ More We reconstruct the dark matter density field from spatially overlapping spectroscopic and photometric redshift catalogs through a forward modelling approach. Instead of directly inferring the underlying density field, we find the best fitting initial Gaussian fluctuations that will evolve into the observed cosmic volume. To account for the substantial uncertainty of photometric redshifts we employ a differentiable continuous Poisson process. In the context of the upcoming Prime Focus Spectrograph (PFS), we find improvements in cosmic structure classification equivalent to 50-100\% more spectroscopic targets by combining relatively sparse spectroscopic with dense photometric samples. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 7 pages, 6 figures

AESTRA: Deep Learning for Precise Radial Velocity Estimation in the Presence of Stellar Activity

Authors: Yan Liang, Joshua N. Winn, Peter Melchior

Abstract: Stellar activity interferes with precise radial velocity measurements and limits our ability to detect and characterize planets, particularly Earth-like planets. We introduce \aestra (Auto-Encoding STellar Radial-velocity and Activity), a deep learning method for precise radial velocity measurements. It combines a spectrum auto-encoder, which learns to create realistic models of the star's rest-fr… ▽ More Stellar activity interferes with precise radial velocity measurements and limits our ability to detect and characterize planets, particularly Earth-like planets. We introduce \aestra (Auto-Encoding STellar Radial-velocity and Activity), a deep learning method for precise radial velocity measurements. It combines a spectrum auto-encoder, which learns to create realistic models of the star's rest-frame spectrum, and a radial-velocity estimator, which learns to identify true Doppler shifts in the presence of spurious shifts due to line-profile variations. Being self-supervised, \aestra does not need "ground truth" radial velocities for training, making it applicable to exoplanet host stars for which the truth is unknown. In tests involving 1,000 simulated spectra, \aestra can detect planetary signals as low as 0.1 m/s even in the presence of 3 m/s of activity-induced noise and 0.3 m/s of photon noise per spectrum. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 14 pages, 12 figures. Accepted for publication in The Astronomical Journal

Multiscale Feature Attribution for Outliers

Authors: Jeff Shen, Peter Melchior

Abstract: Machine learning techniques can automatically identify outliers in massive datasets, much faster and more reproducible than human inspection ever could. But finding such outliers immediately leads to the question: which features render this input anomalous? We propose a new feature attribution method, Inverse Multiscale Occlusion, that is specifically designed for outliers, for which we have littl… ▽ More Machine learning techniques can automatically identify outliers in massive datasets, much faster and more reproducible than human inspection ever could. But finding such outliers immediately leads to the question: which features render this input anomalous? We propose a new feature attribution method, Inverse Multiscale Occlusion, that is specifically designed for outliers, for which we have little knowledge of the type of features we want to identify and expect that the model performance is questionable because anomalous test data likely exceed the limits of the training data. We demonstrate our method on outliers detected in galaxy spectra from the Dark Energy Survey Instrument and find its results to be much more interpretable than alternative attribution approaches. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 6 pages, 2 figures, accepted to NeurIPS 2023 Workshop on Machine Learning and the Physical Sciences. Code available at https://github.com/al-jshen/imo

Constructing Impactful Machine Learning Research for Astronomy: Best Practices for Researchers and Reviewers

Authors: D. Huppenkothen, M. Ntampaka, M. Ho, M. Fouesneau, B. Nord, J. E. G. Peek, M. Walmsley, J. F. Wu, C. Avestruz, T. Buck, M. Brescia, D. P. Finkbeiner, A. D. Goulding, T. Kacprzak, P. Melchior, M. Pasquato, N. Ramachandra, Y. -S. Ting, G. van de Ven, S. Villar, V. A. Villar, E. Zinger

Abstract: Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best pr… ▽ More Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best practices, challenges, and drawbacks, which, at present, are often reported on incompletely in the astrophysical literature. With this paper, we aim to provide a primer to the astronomical community, including authors, reviewers, and editors, on how to implement machine learning models and report their results in a way that ensures the accuracy of the results, reproducibility of the findings, and usefulness of the method. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 14 pages, 3 figures; submitted to the Bulletin of the American Astronomical Society

Cosmology with Galaxy Photometry Alone

Authors: ChangHoon Hahn, Francisco Villaescusa-Navarro, Peter Melchior, Romain Teyssier

Abstract: We present the first cosmological constraints using only the observed photometry of galaxies. Villaescusa-Navarro et al. (2022; arXiv:2201.02202) recently demonstrated that the internal physical properties of a single simulated galaxy contain a significant amount of cosmological information. These physical properties, however, cannot be directly measured from observations. In this work, we present… ▽ More We present the first cosmological constraints using only the observed photometry of galaxies. Villaescusa-Navarro et al. (2022; arXiv:2201.02202) recently demonstrated that the internal physical properties of a single simulated galaxy contain a significant amount of cosmological information. These physical properties, however, cannot be directly measured from observations. In this work, we present how we can go beyond theoretical demonstrations to infer cosmological constraints from actual galaxy observables (e.g. optical photometry) using neural density estimation and the CAMELS suite of hydrodynamical simulations. We find that the cosmological information in the photometry of a single galaxy is limited. However, we combine the constraining power of photometry from many galaxies using hierarchical population inference and place significant cosmological constraints. With the observed photometry of $\sim$20,000 NASA-Sloan Atlas galaxies, we constrain $Ω_m = 0.323^{+0.075}_{-0.095}$ and $σ_8 = 0.799^{+0.088}_{-0.085}$. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 15 pages, 7 figures, submitted to ApJL, comments welcome

PopSED: Population-Level Inference for Galaxy Properties from Broadband Photometry with Neural Density Estimation

Authors: Jiaxuan Li, Peter Melchior, ChangHoon Hahn, Song Huang

Abstract: We present PopSED, a framework for the population-level inference of galaxy properties from photometric data. Unlike the traditional approach of first analyzing individual galaxies and then combining the results to determine the physical properties of the entire galaxy population, we directly make the population distribution the inference objective. We train normalizing flows to approximate the po… ▽ More We present PopSED, a framework for the population-level inference of galaxy properties from photometric data. Unlike the traditional approach of first analyzing individual galaxies and then combining the results to determine the physical properties of the entire galaxy population, we directly make the population distribution the inference objective. We train normalizing flows to approximate the population distribution by minimizing the Wasserstein distance between the synthetic photometry of the galaxy population and the observed data. We validate our method using mock observations and apply it to galaxies from the GAMA survey. PopSED reliably recovers the redshift and stellar mass distribution of $10^{5}$ galaxies using broadband photometry within $<1$ GPU hr, being $10^{5-6}$ times faster than the traditional spectral energy distribution modeling method. From the population posterior, we also recover the star-forming main sequence for GAMA galaxies at $z<0.1$. With the unprecedented number of galaxies in upcoming surveys, our method offers an efficient tool for studying galaxy evolution and deriving redshift distributions for cosmological analyses. △ Less

Submitted 19 December, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 13 pages, 4 figures. AJ accepted

Outlier Detection in the DESI Bright Galaxy Survey

Authors: Yan Liang, Peter Melchior, ChangHoon Hahn, Jeff Shen, Andy Goulding, Charlotte Ward

Abstract: We present an unsupervised search for outliers in the Bright Galaxy Survey (BGS) dataset from the DESI Early Data Release. This analysis utilizes an autoencoder to compress galaxy spectra into a compact, redshift-invariant latent space, and a normalizing flow to identify low-probability objects. The most prominent outliers show distinctive spectral features such as irregular or double-peaked emiss… ▽ More We present an unsupervised search for outliers in the Bright Galaxy Survey (BGS) dataset from the DESI Early Data Release. This analysis utilizes an autoencoder to compress galaxy spectra into a compact, redshift-invariant latent space, and a normalizing flow to identify low-probability objects. The most prominent outliers show distinctive spectral features such as irregular or double-peaked emission lines, or originate from galaxy mergers, blended sources, and rare quasar types, including one previously unknown Broad Absorption Line system. A significant portion of the BGS outliers are stars spectroscopically misclassified as galaxies. By building our own star model trained on spectra from the DESI Milky Way Survey, we have determined that the misclassification likely stems from the Principle Component Analysis of stars in the DESI pipeline. To aid follow-up studies, we make the full probability catalog of all BGS objects and our pre-trained models publicly available. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures

Spotting Hallucinations in Inverse Problems with Data-Driven Priors

Authors: Matt L. Sampson, Peter Melchior

Abstract: Hallucinations are an inescapable consequence of solving inverse problems with deep neural networks. The expressiveness of recent generative models is the reason why they can yield results far superior to conventional regularizers; it can also lead to realistic-looking but incorrect features, potentially undermining the trust in important aspects of the reconstruction. We present a practical and c… ▽ More Hallucinations are an inescapable consequence of solving inverse problems with deep neural networks. The expressiveness of recent generative models is the reason why they can yield results far superior to conventional regularizers; it can also lead to realistic-looking but incorrect features, potentially undermining the trust in important aspects of the reconstruction. We present a practical and computationally efficient method to determine, which regions in the solutions of inverse problems with data-driven priors are prone to hallucinations. By computing the diagonal elements of the Fisher information matrix of the likelihood and the data-driven prior separately, we can flag regions where the information is prior-dominated. Our diagnostic can directly be compared to the reconstructed solutions and enables users to decide if measurements in such regions are robust for their application. Our method scales linearly with the number of parameters and is thus applicable in high-dimensional settings, allowing it to be rolled out broadly for the large-volume data products of future wide-field surveys. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 7 pages 3 figures, Accepted at the ICML 2023 Workshop on Machine Learning for Astrophysics

Cool Cores in Clusters of Galaxies in the Dark Energy Survey

Authors: K. Graham, J. O'Donnell, M. M. Silverstein, O. Eiger, T. E. Jeltema, D. L. Hollowood, D. Cross, S. Everett, P. Giles, J. Jobel, D. Laubner, A. McDaniel, A. K. Romer, A. Swart, M. Aguena, S. Allam, O. Alves, D. Brooks, M. Carrasco Kind, J. Carretero, M. Costanzi, L. N. da Costa, M. E. S. Pereira, J. De Vicente, S. Desai , et al. (32 additional authors not shown)

Abstract: We search for the presence of cool cores in optically-selected galaxy clusters from the Dark Energy Survey (DES) and investigate their prevalence as a function of redshift and cluster richness. Clusters were selected from the redMaPPer analysis of three years of DES observations that have archival Chandra X-ray observations, giving a sample of 99 clusters with a redshift range of… ▽ More We search for the presence of cool cores in optically-selected galaxy clusters from the Dark Energy Survey (DES) and investigate their prevalence as a function of redshift and cluster richness. Clusters were selected from the redMaPPer analysis of three years of DES observations that have archival Chandra X-ray observations, giving a sample of 99 clusters with a redshift range of $0.11 < z < 0.87$ and a richness range of $25 < λ< 207$. Using the X-ray data, the core temperature was compared to the outer temperature to identify clusters where the core temperature is a factor of 0.7 or less than the outer temperature. We found a cool core fraction of approximately 20% with no significant trend in the cool core fraction with either redshift or richness. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: shortened version accepted to RNAAS

The Intrinsic Alignment of Red Galaxies in DES Y1 redMaPPer Galaxy Clusters

Authors: C. Zhou, A. Tong, M. A. Troxel, J. Blazek, C. Lin, D. Bacon, L. Bleem, A. Carnero Rosell, C. Chang, M. Costanzi, J. DeRose, J. P. Dietrich, A. Drlica-Wagner, D. Gruen, R. A. Gruendl, B. Hoyle, M. Jarvis, N. MacCrann, B. Mawdsley, T. McClintock, P. Melchior, J. Prat, A. Pujol, E. Rozo, E. S. Rykoff , et al. (57 additional authors not shown)

Abstract: Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We meas… ▽ More Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We measure this intrinsic alignment in Dark Energy Survey (DES) Year 1 redMaPPer clusters. We find evidence of a non-zero mean radial alignment of galaxies within clusters between redshift 0.1-0.7. We find a significant systematic in the measured ellipticities of cluster satellite galaxies that we attribute to the central galaxy flux and other intracluster light. We attempt to correct this signal, and fit a simple model for intrinsic alignment amplitude ($A_{\textrm{IA}}$) to the measurement, finding $A_{\textrm{IA}}=0.15\pm 0.04$, when excluding data near the edge of the cluster. We find a significantly stronger alignment of the central galaxy with the cluster dark matter halo at low redshift and with higher richness and central galaxy absolute magnitude (proxies for cluster mass). This is an important demonstration of the ability of large photometric data sets like DES to provide direct constraints on the intrinsic alignment of galaxies within clusters. These measurements can inform improvements to small-scale modeling and simulation of the intrinsic alignment of galaxies to help improve the separation of the intrinsic alignment signal in weak lensing studies. △ Less

Submitted 5 September, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 13 figures. Accepted to MNRAS

Rates and properties of type Ia supernovae in galaxy clusters within the Dark Energy Survey

Authors: M. Toy, P. Wiseman, M. Sullivan, C. Frohmaier, O. Graur, A. Palmese, B. Popovic, T. M. Davis, L. Galbany, L. Kelsey, C. Lidman, D. Scolnic, S. Allam, S. Desai, T. M. C. Abbott, M. Aguena, O. Alves, J. Annis, D. Bacon, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero , et al. (37 additional authors not shown)

Abstract: We identify 66 photometrically classified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) that have occurred within red-sequence selected galaxy clusters. We compare light-curve and host galaxy properties of the cluster SNe to 1024 DES SNe Ia located in field galaxies, the largest comparison of two such samples at high redshift (z > 0.1). We find that cluster SN light curves decline… ▽ More We identify 66 photometrically classified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) that have occurred within red-sequence selected galaxy clusters. We compare light-curve and host galaxy properties of the cluster SNe to 1024 DES SNe Ia located in field galaxies, the largest comparison of two such samples at high redshift (z > 0.1). We find that cluster SN light curves decline faster than those in the field (97.7 per cent confidence). However, when limiting these samples to host galaxies of similar colour and mass, there is no significant difference in the SN light curve properties. Motivated by previous detections of a higher-normalised SN Ia delay time distribution in galaxy clusters, we measure the intrinsic rate of SNe Ia in cluster and field environments. We find the average ratio of the SN Ia rate per galaxy between high mass ($10\leq\log\mathrm{(M_{*}/M_{\odot})} \leq 11.25$) cluster and field galaxies to be $0.594 \pm0.068$. This difference is mass-dependent, with the ratio declining with increasing mass, which suggests that the stellar populations in cluster hosts are older than those in field hosts. We show that the mass-normalised rate (or SNe per unit mass) in massive-passive galaxies is consistent between cluster and field environments. Additionally, both of these rates are consistent with rates previously measured in clusters at similar redshifts. We conclude that in massive-passive galaxies, which are the dominant hosts of cluster SNe, the cluster DTD is comparable to the field. △ Less

Submitted 28 September, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Designing an Optimal Kilonova Search using DECam for Gravitational Wave Events

Authors: C. R. Bom, J. Annis, A. Garcia, A. Palmese, N. Sherman, M. Soares-Santos, L. Santana-Silva, R. Morgan, K. Bechtol, T. Davis, H. T. Diehl, S. S. Allam, T. G. Bachmann, B. M. O. Fraga, J. Garcıa-Bellido, M. S. S. Gill, K. Herner, C. D. Kilpatrick, M. Makler, F. Olivares E., M. E. S. Pereira, J. Pineda, A. Santos, D. L. Tucker, M. P. Wiesner , et al. (45 additional authors not shown)

Abstract: We address the problem of optimally identifying all kilonovae detected via gravitational wave emission in the upcoming LIGO/Virgo/KAGRA Collaboration observing run, O4, which is expected to be sensitive to a factor of $\sim 7$ more Binary Neutron Stars alerts than previously. Electromagnetic follow-up of all but the brightest of these new events will require $>1$ meter telescopes, for which limite… ▽ More We address the problem of optimally identifying all kilonovae detected via gravitational wave emission in the upcoming LIGO/Virgo/KAGRA Collaboration observing run, O4, which is expected to be sensitive to a factor of $\sim 7$ more Binary Neutron Stars alerts than previously. Electromagnetic follow-up of all but the brightest of these new events will require $>1$ meter telescopes, for which limited time is available. We present an optimized observing strategy for the Dark Energy Camera during O4. We base our study on simulations of gravitational wave events expected for O4 and wide-prior kilonova simulations. We derive the detectabilities of events for realistic observing conditions. We optimize our strategy for confirming a kilonova while minimizing telescope time. For a wide range of kilonova parameters, corresponding to a fainter kilonova compared to GW170817/AT2017gfo we find that, with this optimal strategy, the discovery probability for electromagnetic counterparts with the Dark Energy Camera is $\sim 80\%$ at the nominal binary neutron star gravitational wave detection limit for the next LVK observing run (190 Mpc), which corresponds to a $\sim 30\%$ improvement compared to the strategy adopted during the previous observing run. For more distant events ($\sim 330$ Mpc), we reach a $\sim 60\%$ probability of detection, a factor of $\sim 2$ increase. For a brighter kilonova model dominated by the blue component that reproduces the observations of GW170817/AT2017gfo, we find that we can reach $\sim 90\%$ probability of detection out to 330 Mpc, representing an increase of $\sim 20 \%$, while also reducing the total telescope time required to follow-up events by $\sim 20\%$. △ Less

Submitted 1 November, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Comments: 26 pages, 11 figures, accepted by ApJ

Report number: DES-2022-0714, FERMILAB-PUB-23-048-PPD

Autoencoding Galaxy Spectra II: Redshift Invariance and Outlier Detection

Authors: Yan Liang, Peter Melchior, Sicong Lu, Andy Goulding, Charlotte Ward

Abstract: We present an unsupervised outlier detection method for galaxy spectra based on the spectrum autoencoder architecture spender, which reliably captures spectral features and provides highly realistic reconstructions for SDSS galaxy spectra. We interpret the sample density in the autoencoder latent space as a probability distribution, and identify outliers as low-probability objects with a normalizi… ▽ More We present an unsupervised outlier detection method for galaxy spectra based on the spectrum autoencoder architecture spender, which reliably captures spectral features and provides highly realistic reconstructions for SDSS galaxy spectra. We interpret the sample density in the autoencoder latent space as a probability distribution, and identify outliers as low-probability objects with a normalizing flow. However, we found that the latent-space position is not, as expected from the architecture, redshift invariant, which introduces stochasticity into the latent space and the outlier detection method. We solve this problem by adding two novel loss terms during training, which explicitly link latent-space distances to data-space distances, preserving locality in the autoencoding process. Minimizing the additional losses leads to a redshift-invariant, non-degenerate latent space distribution with clear separations between common and anomalous data. We inspect the spectra with the lowest probability and find them to include blends with foreground stars, extremely reddened galaxies, galaxy pairs and triples, and stars that are misclassified as galaxies. We release the newly trained spender model and the latent-space probability for the entire SDSS-I galaxy sample to aid further investigations. △ Less

Submitted 5 February, 2023; originally announced February 2023.

Comments: 12 pages, 7 figures, submitted to AJ. Code available at https://github.com/pmelchior/spender

The Dark Energy Survey Year 3 and eBOSS: constraining galaxy intrinsic alignments across luminosity and colour space

Authors: S. Samuroff, R. Mandelbaum, J. Blazek, A. Campos, N. MacCrann, G. Zacharegkas, A. Amon, J. Prat, S. Singh, J. Elvin-Poole, A. J. Ross, A. Alarcon, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, C. Chang, R. Chen, A. Choi, M. Crocce, C. Davis, J. DeRose , et al. (82 additional authors not shown)

Abstract: We present direct constraints on galaxy intrinsic alignments using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift $z\sim0.2-1.0$, luminous red galaxies (LRGs) from eBOSS a… ▽ More We present direct constraints on galaxy intrinsic alignments using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift $z\sim0.2-1.0$, luminous red galaxies (LRGs) from eBOSS at $z\sim0.8$, and also a SDSS-III BOSS CMASS sample at $z\sim0.5$. We measure two point intrinsic alignment correlations, which we fit using a model that includes lensing, magnification and photometric redshift error. Fitting on scales $6<r_{\rm p} < 70$ Mpc$/h$, we make a detection of intrinsic alignments in each sample, at $5σ-22σ$ (assuming a simple one parameter model for IAs). Using these red samples, we measure the IA-luminosity relation. Our results are statistically consistent with previous results, but offer a significant improvement in constraining power, particularly at low luminosity. With this improved precision, we see detectable dependence on colour between broadly defined red samples. It is likely that a more sophisticated approach than a binary red/blue split, which jointly considers colour and luminosity dependence in the IA signal, will be needed in future. We also compare the various signal components at the best fitting point in parameter space for each sample, and find that magnification and lensing contribute $\sim2-18\%$ of the total signal. As precision continues to improve, it will certainly be necessary to account for these effects in future direct IA measurements. Finally, we make equivalent measurements on a sample of Emission Line Galaxies (ELGs) from eBOSS at $z\sim 0.8$. We report a null detection, constraining the IA amplitude (assuming the nonlinear alignment model) to be $A_1=0.07^{+0.32}_{-0.42}$ ($|A_1|<0.78$ at $95\%$ CL). △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Submitted to MNRAS; 22 pages, 11 figures

Report number: FERMILAB-PUB-22-914-PPD

Plausible Adversarial Attacks on Direct Parameter Inference Models in Astrophysics

Authors: Benjamin Horowitz, Peter Melchior

Abstract: In this abstract we explore the possibility of introducing biases in physical parameter inference models from adversarial-type attacks. In particular, we inject small amplitude systematics into inputs to a mixture density networks tasked with inferring cosmological parameters from observed data. The systematics are constructed analogously to white-box adversarial attacks. We find that the analysis… ▽ More In this abstract we explore the possibility of introducing biases in physical parameter inference models from adversarial-type attacks. In particular, we inject small amplitude systematics into inputs to a mixture density networks tasked with inferring cosmological parameters from observed data. The systematics are constructed analogously to white-box adversarial attacks. We find that the analysis network can be tricked into spurious detection of new physics in cases where standard cosmological estimators would be insensitive. This calls into question the robustness of such networks and their utility for reliably detecting new physics. △ Less

Submitted 27 November, 2022; originally announced November 2022.

Comments: Accepted submission to Machine Learning and the Physical Sciences workshop, NeurIPS 2022

A Sample of Dust Attenuation Laws for DES Supernova Host Galaxies

Authors: J. Duarte, S. González-Gaitán, A. Mourao, A. Paulino-Afonso, P. Guilherme-Garcia, J. Aguas, L. Galbany, L. Kelsey, D. Scolnic, M. Sullivan, D. Brout, A. Palmese, P. Wiseman, A. Pieres, A. A. Plazas Malagón, A. Carnero Rosell, C. To, D. Gruen, D. Bacon, D. Brooks, D. L. Burke, D. W. Gerdes, D. J. James, D. L. Hollowood, D. Friedel , et al. (36 additional authors not shown)

Abstract: Type Ia supernovae (SNe Ia) are useful distance indicators in cosmology, provided their luminosity is standardized by applying empirical corrections based on light-curve properties. One factor behind these corrections is dust extinction, accounted for in the color-luminosity relation of the standardization. This relation is usually assumed to be universal, which could potentially introduce systema… ▽ More Type Ia supernovae (SNe Ia) are useful distance indicators in cosmology, provided their luminosity is standardized by applying empirical corrections based on light-curve properties. One factor behind these corrections is dust extinction, accounted for in the color-luminosity relation of the standardization. This relation is usually assumed to be universal, which could potentially introduce systematics into the standardization. The ``mass-step'' observed for SNe Ia Hubble residuals has been suggested as one such systematic. We seek to obtain a completer view of dust attenuation properties for a sample of 162 SN Ia host galaxies and to probe their link to the ``mass-step''. We infer attenuation laws towards hosts from both global and local (4 kpc) Dark Energy Survey photometry and Composite Stellar Population model fits. We recover a optical depth/attenuation slope relation, best explained by differing star/dust geometry for different galaxy orientations, which is significantly different from the optical depth/extinction slope relation observed directly for SNe. We obtain a large variation of attenuation slopes and confirm these change with host properties, like stellar mass and age, meaning a universal SN Ia correction should ideally not be assumed. Analyzing the cosmological standardization, we find evidence for a ``mass-step'' and a two dimensional ``dust-step'', both more pronounced for red SNe. Although comparable, the two steps are found no to be completely analogous. We conclude that host galaxy dust data cannot fully account for the ``mass-step'', using either an alternative SN standardization with extinction proxied by host attenuation or a ``dust-step'' approach. △ Less

Submitted 19 December, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 20 pages, 10 figues, 9 tables. Supplementary material included (10 pages). Accepted for publication on A&A

Report number: DES-2022-0694; FERMILAB-PUB-22-760-PPD

Journal ref: A&A 680, A56 (2023)

Autoencoding Galaxy Spectra I: Architecture

Authors: Peter Melchior, Yan Liang, ChangHoon Hahn, Andy Goulding

Abstract: We introduce the neural network architecture SPENDER as a core differentiable building block for analyzing, representing, and creating galaxy spectra. It combines a convolutional encoder, which pays attention to up to 256 spectral features and compresses them into a low-dimensional latent space, with a decoder that generates a restframe representation, whose spectral range and resolution exceeds t… ▽ More We introduce the neural network architecture SPENDER as a core differentiable building block for analyzing, representing, and creating galaxy spectra. It combines a convolutional encoder, which pays attention to up to 256 spectral features and compresses them into a low-dimensional latent space, with a decoder that generates a restframe representation, whose spectral range and resolution exceeds that of the observing instrument. The decoder is followed by explicit redshift, resampling, and convolution transformations to match the observations. The architecture takes galaxy spectra at arbitrary redshifts and is robust to glitches like residuals of the skyline subtraction, so that spectra from a large survey can be ingested directly without additional preprocessing. We demonstrate the performance of SPENDER by training on the entire spectroscopic galaxy sample of SDSS-II; show its ability to create highly accurate reconstructions with substantially reduced noise; perform deconvolution and oversampling for a super-resolution model that resolves the [OII] doublet; introduce a novel method to interpret attention weights as proxies for important spectral features; and infer the main degrees of freedom represented in the latent space. We conclude with a discussion of future improvements and applications. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 17 pages, 5 figures, submitted to AJ. Code available at https://github.com/pmelchior/spender

Deblending Galaxies with Generative Adversarial Networks

Authors: Shoubaneh Hemmati, Eric Huff, Hooshang Nayyeri, Agnès Ferté, Peter Melchior, Bahram Mobasher, Jason Rhodes, Abtin Shahidi, Harry Teplitz

Abstract: Deep generative models including generative adversarial networks (GANs) are powerful unsupervised tools in learning the distributions of data sets. Building a simple GAN architecture in PyTorch and training on the CANDELS data set, we generate galaxy images with the Hubble Space Telescope resolution starting from a noise vector. We proceed by modifying the GAN architecture to improve the Subaru Hy… ▽ More Deep generative models including generative adversarial networks (GANs) are powerful unsupervised tools in learning the distributions of data sets. Building a simple GAN architecture in PyTorch and training on the CANDELS data set, we generate galaxy images with the Hubble Space Telescope resolution starting from a noise vector. We proceed by modifying the GAN architecture to improve the Subaru Hyper Suprime-Cam ground-based images by increasing their resolution to the HST resolution. We use the super resolution GAN on a large sample of blended galaxies which we create using CANDELS cutouts. In our simulated blend sample, $\sim 20 \%$ would unrecognizably be blended even in the HST resolution cutouts. In the HSC-like cutouts this fraction rises to $\sim 90\%$. With our modified GAN we can lower this value to $\sim 50\%$. We quantify the blending fraction in the high, low and GAN resolutions over the whole manifold of angular separation, flux ratios, sizes and redshift difference between the two blended objects. The two peaks found by the GAN deblender result in ten times improvement in the photometry measurement of the blended objects. Modifying the architecture of the GAN, we also train a Multi-wavelength GAN with seven band optical+NIR HST cutouts. This multi-wavelength GAN improves the fraction of detected blends by another $\sim 10\%$ compared to the single-band GAN. This is most beneficial to the current and future precision cosmology experiments (e.g., LSST, SPHEREx, Euclid, Roman), specifically those relying on weak gravitational lensing, where blending is a major source of systematic error. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 12 pages, 9 figures, accepted for publication in the Astrophysical Journal

Beyond Ultra-diffuse Galaxies. I. Mass--Size Outliers among the Satellites of Milky Way Analogs

Authors: Jiaxuan Li, Jenny E. Greene, Johnny P. Greco, Song Huang, Peter Melchior, Rachael Beaton, Kirsten Casey, Shany Danieli, Andy Goulding, Remy Joseph, Erin Kado-Fong, Ji Hoon Kim, Lauren A. MacArthur

Abstract: Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru's Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and mode… ▽ More Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru's Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and modeling them. As a case study, we surveyed 922 Milky Way analogs in the nearby Universe ($0.01 < z < 0.04$) and build a large sample of satellite galaxies that are outliers in the mass-size relation. These ``ultra-puffy'' galaxies (UPGs), defined to be $1.5σ$ above the average mass-size relation, represent the tail of the satellite size distribution. We find that each MW analog hosts $N_{\rm UPG} = 0.31\pm 0.05$ ultra-puffy galaxies on average, which is consistent with but slightly lower than the observed abundance at this halo mass in the Local Volume. We also construct a sample of ultra-diffuse galaxies (UDGs) in MW analogs and find an abundance of $N_{\rm UDG} = 0.44\pm0.05$ per host. With literature results, we confirm that the UDG abundance scales with the host halo mass following a sublinear power law. We argue that our definition for ultra-puffy galaxies, which is based on the mass-size relation, is more physically-motivated than the common definition of ultra-diffuse galaxies, which depends on surface brightness and size cuts and thus yields different surface mass density cuts for quenched and star-forming galaxies. △ Less

Submitted 11 September, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, ApJ accepted. Also see the companion paper arXiv:2302.14108

$\texttt{Mangrove}$: Learning Galaxy Properties from Merger Trees

Authors: Christian Kragh Jespersen, Miles Cranmer, Peter Melchior, Shirley Ho, Rachel S. Somerville, Austen Gabrielpillai

Abstract: Efficiently mapping baryonic properties onto dark matter is a major challenge in astrophysics. Although semi-analytic models (SAMs) and hydrodynamical simulations have made impressive advances in reproducing galaxy observables across cosmologically significant volumes, these methods still require significant computation times, representing a barrier to many applications. Graph Neural Networks (GNN… ▽ More Efficiently mapping baryonic properties onto dark matter is a major challenge in astrophysics. Although semi-analytic models (SAMs) and hydrodynamical simulations have made impressive advances in reproducing galaxy observables across cosmologically significant volumes, these methods still require significant computation times, representing a barrier to many applications. Graph Neural Networks (GNNs) have recently proven to be the natural choice for learning physical relations. Among the most inherently graph-like structures found in astrophysics are the dark matter merger trees that encode the evolution of dark matter halos. In this paper we introduce a new, graph-based emulator framework, $\texttt{Mangrove}$, and show that it emulates the galactic stellar mass, cold gas mass and metallicity, instantaneous and time-averaged star formation rate, and black hole mass -- as predicted by a SAM -- with root mean squared error up to two times lower than other methods across a $(75 Mpc/h)^3$ simulation box in 40 seconds, 4 orders of magnitude faster than the SAM. We show that $\texttt{Mangrove}$ allows for quantification of the dependence of galaxy properties on merger history. We compare our results to the current state of the art in the field and show significant improvements for all target properties. $\texttt{Mangrove}$ is publicly available. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: 15 pages, 9 figures, 3 tables, 10 pages of Appendices. Accepted for publication in ApJ

Identification of Galaxy-Galaxy Strong Lens Candidates in the DECam Local Volume Exploration Survey Using Machine Learning

Authors: E. A. Zaborowski, A. Drlica-Wagner, F. Ashmead, J. F. Wu, R. Morgan, C. R. Bom, A. J. Shajib, S. Birrer, W. Cerny, L. Buckley-Geer, B. Mutlu-Pakdil, P. S. Ferguson, K. Glazebrook, S. J. Gonzalez Lozano, Y. Gordon, M. Martinez, V. Manwadkar, J. O'Donnell, J. Poh, A. Riley, J. D. Sakowska, L. Santana-Silva, B. X. Santiago, D. Sluse, C. Y. Tan , et al. (66 additional authors not shown)

Abstract: We perform a search for galaxy-galaxy strong lens systems using a convolutional neural network (CNN) applied to imaging data from the first public data release of the DECam Local Volume Exploration Survey (DELVE), which contains $\sim 520$ million astronomical sources covering $\sim 4,000$ $\mathrm{deg}^2$ of the southern sky to a $5σ$ point-source depth of $g=24.3$, $r=23.9$, $i=23.3$, and… ▽ More We perform a search for galaxy-galaxy strong lens systems using a convolutional neural network (CNN) applied to imaging data from the first public data release of the DECam Local Volume Exploration Survey (DELVE), which contains $\sim 520$ million astronomical sources covering $\sim 4,000$ $\mathrm{deg}^2$ of the southern sky to a $5σ$ point-source depth of $g=24.3$, $r=23.9$, $i=23.3$, and $z=22.8$ mag. Following the methodology of similar searches using DECam data, we apply color and magnitude cuts to select a catalog of $\sim 11$ million extended astronomical sources. After scoring with our CNN, the highest scoring 50,000 images were visually inspected and assigned a score on a scale from 0 (definitely not a lens) to 3 (very probable lens). We present a list of 581 strong lens candidates, 562 of which are previously unreported. We categorize our candidates using their human-assigned scores, resulting in 55 Grade A candidates, 149 Grade B candidates, and 377 Grade C candidates. We additionally highlight eight potential quadruply lensed quasars from this sample. Due to the location of our search footprint in the northern Galactic cap ($b > 10$ deg) and southern celestial hemisphere (${\rm Dec.}<0$ deg), our candidate list has little overlap with other existing ground-based searches. Where our search footprint does overlap with other searches, we find a significant number of high-quality candidates which were previously unidentified, indicating a degree of orthogonality in our methodology. We report properties of our candidates including apparent magnitude and Einstein radius estimated from the image separation. △ Less

Submitted 25 August, 2023; v1 submitted 19 October, 2022; originally announced October 2022.

Comments: 24 pages; published version (ApJ)

Journal ref: ApJ 954 68 (2023)

Mapping gas around massive galaxies: cross-correlation of DES Y3 galaxies and Compton-$y$-maps from SPT and Planck

Authors: J. Sánchez, Y. Omori, C. Chang, L. E. Bleem, T. Crawford, A. Drlica-Wagner, S. Raghunathan, G. Zacharegkas, T. M. C. Abbott, M. Aguena, A. Alarcon, S. Allam, O. Alves, A. Amon, S. Avila, E. Baxter, K. Bechtol, B. A. Benson, G. M. Bernstein, E. Bertin, S. Bocquet, D. Brooks, D. L. Burke, A. Campos, J. E. Carlstrom , et al. (102 additional authors not shown)

Abstract: We cross-correlate positions of galaxies measured in data from the first three years of the Dark Energy Survey with Compton-$y$-maps generated using data from the South Pole Telescope (SPT) and the {\it Planck} mission. We model this cross-correlation measurement together with the galaxy auto-correlation to constrain the distribution of gas in the Universe. We measure the hydrostatic mass bias or,… ▽ More We cross-correlate positions of galaxies measured in data from the first three years of the Dark Energy Survey with Compton-$y$-maps generated using data from the South Pole Telescope (SPT) and the {\it Planck} mission. We model this cross-correlation measurement together with the galaxy auto-correlation to constrain the distribution of gas in the Universe. We measure the hydrostatic mass bias or, equivalently, the mean halo bias-weighted electron pressure $\langle b_{h}P_{e}\rangle$, using large-scale information. We find $\langle b_{h}P_{e}\rangle$ to be $[0.16^{+0.03}_{-0.04},0.28^{+0.04}_{-0.05},0.45^{+0.06}_{-0.10},0.54^{+0.08}_{-0.07},0.61^{+0.08}_{-0.06},0.63^{+0.07}_{-0.08}]$ meV cm$^{-3}$ at redshifts $z \sim [0.30, 0.46, 0.62,0.77, 0.89, 0.97]$. These values are consistent with previous work where measurements exist in the redshift range. We also constrain the mean gas profile using small-scale information, enabled by the high-resolution of the SPT data. We compare our measurements to different parametrized profiles based on the cosmo-OWLS hydrodynamical simulations. We find that our data are consistent with the simulation that assumes an AGN heating temperature of $10^{8.5}$K but are incompatible with the model that assumes an AGN heating temperature of $10^{8.0}$K. These comparisons indicate that the data prefer a higher value of electron pressure than the simulations within $r_{500c}$ of the galaxies' halos. △ Less

Submitted 18 October, 2022; v1 submitted 16 October, 2022; originally announced October 2022.

Comments: 20 pages, 9 figures. Submitted to MNRAS

Report number: FERMILAB-PUB-22-603-PPD

The PSZ-MCMF catalogue of Planck clusters over the DES region

Authors: D. Hernández-Lang, M. Klein, J. J. Mohr, S. Grandis, J. -B. Melin, P. Tarrío, M. Arnaud, G. W. Pratt, T. M. C. Abbott, M. Aguena, O. Alves, F. Andrade-Oliveira, D. Bacon, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, F. J. Castander, M. Costanzi, L. N. da Costa, M. E. S. Pereira, S. Desai, H. T. Diehl , et al. (36 additional authors not shown)

Abstract: We present the first systematic follow-up of Planck Sunyaev-Zeldovich effect (SZE) selected candidates down to signal-to-noise (S/N) of 3 over the 5000 deg$^2$ covered by the Dark Energy Survey. Using the MCMF cluster confirmation algorithm, we identify optical counterparts, determine photometric redshifts and richnesses and assign a parameter, $f_{\rm cont}$, that reflects the probability that ea… ▽ More We present the first systematic follow-up of Planck Sunyaev-Zeldovich effect (SZE) selected candidates down to signal-to-noise (S/N) of 3 over the 5000 deg$^2$ covered by the Dark Energy Survey. Using the MCMF cluster confirmation algorithm, we identify optical counterparts, determine photometric redshifts and richnesses and assign a parameter, $f_{\rm cont}$, that reflects the probability that each SZE-optical pairing represents a random superposition of physically unassociated systems rather than a real cluster. The new PSZ-MCMF cluster catalogue consists of 853 MCMF confirmed clusters and has a purity of 90%. We present the properties of subsamples of the PSZ-MCMF catalogue that have purities ranging from 90% to 97.5%, depending on the adopted $f_{\rm cont}$ threshold. Halo mass estimates $M_{500}$, redshifts, richnesses, and optical centers are presented for all PSZ-MCMF clusters. The PSZ-MCMF catalogue adds 589 previously unknown Planck identified clusters over the DES footprint and provides redshifts for an additional 50 previously published Planck selected clusters with S/N>4.5. Using the subsample with spectroscopic redshifts, we demonstrate excellent cluster photo-$z$ performance with an RMS scatter in $Δz/(1+z)$ of 0.47%. Our MCMF based analysis allows us to infer the contamination fraction of the initial S/N>3 Planck selected candidate list, which is ~50%. We present a method of estimating the completeness of the PSZ-MCMF cluster sample. In comparison to the previously published Planck cluster catalogues. this new S/N>3 MCMF confirmed cluster catalogue populates the lower mass regime at all redshifts and includes clusters up to z$\sim$1.3. △ Less

Submitted 25 August, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 20 pages, 5 Appendices, 17 figures. The newest submission matches the accepted MNRAS version

Mapping Variations of Redshift Distributions with Probability Integral Transforms

Authors: J. Myles, D. Gruen, A. Amon, A. Alarcon, J. DeRose, S. Everett, S. Dodelson, G. M. Bernstein, A. Campos, I. Harrison, N. MacCrann, J. McCullough, M. Raveri, C. Sánchez, M. A. Troxel, B. Yin, T. M. C. Abbott, S. Allam, O. Alves, F. Andrade-Oliveira, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind , et al. (44 additional authors not shown)

Abstract: We present a method for mapping variations between probability distribution functions and apply this method within the context of measuring galaxy redshift distributions from imaging survey data. This method, which we name PITPZ for the probability integral transformations it relies on, uses a difference in curves between distribution functions in an ensemble as a transformation to apply to anothe… ▽ More We present a method for mapping variations between probability distribution functions and apply this method within the context of measuring galaxy redshift distributions from imaging survey data. This method, which we name PITPZ for the probability integral transformations it relies on, uses a difference in curves between distribution functions in an ensemble as a transformation to apply to another distribution function, thus transferring the variation in the ensemble to the latter distribution function. This procedure is broadly applicable to the problem of uncertainty propagation. In the context of redshift distributions, for example, the uncertainty contribution due to certain effects can be studied effectively only in simulations, thus necessitating a transfer of variation measured in simulations to the redshift distributions measured from data. We illustrate the use of PITPZ by using the method to propagate photometric calibration uncertainty to redshift distributions of the Dark Energy Survey Year 3 weak lensing source galaxies. For this test case, we find that PITPZ yields a lensing amplitude uncertainty estimate due to photometric calibration error within 1 per cent of the truth, compared to as much as a 30 per cent underestimate when using traditional methods. △ Less

Submitted 4 February, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Dark Energy Survey Year 3 Results: Redshift Calibration of the MagLim Lens Sample from the combination of SOMPZ and clustering and its impact on Cosmology

Authors: G. Giannini, A. Alarcon, M. Gatti, A. Porredon, M. Crocce, G. M. Bernstein, R. Cawthon, C. Sánchez, C. Doux, J. Elvin-Poole, M. Raveri, J. Myles, A. Amon, S. Allam, O. Alves, F. Andrade-Oliveira, E. Baxter, K. Bechtol, M. R. Becker, J. Blazek, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, A. Choi , et al. (89 additional authors not shown)

Abstract: We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift ca… ▽ More We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift calibration of the lens sample. We describe in detail the methodology, we validate it on simulations and discuss the main effects dominating our error budget. The new calibration is in fair agreement with the fiducial DES Y3 redshift distributions calibration, with only mild differences ($<3σ$) in the means and widths of the distributions. We study the impact of this new calibration on cosmological constraints, analysing DES Y3 galaxy clustering and galaxy-galaxy lensing measurements, assuming a $Λ$CDM cosmology. We obtain $Ω_{\rm m} = 0.30\pm 0.04$, $σ_8 = 0.81\pm 0.07 $ and $S_8 = 0.81\pm 0.04$, which implies a $\sim 0.4σ$ shift in the $Ω_{\rm}-S_8$ plane compared to the fiducial DES Y3 results, highlighting the importance of the redshift calibration of the lens sample in multi-probe cosmological analyses. △ Less

Submitted 18 October, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Journal ref: MNRAS, stad2945 (2023)

Characterising the Intracluster Light over the Redshift Range $0.2 < z < 0.8$ in the DES-ACT Overlap

Authors: Jesse B. Golden-Marx, Y. Zhang, R. L. C. Ogando, S. Allam, D. L. Tucker, C. J. Miller, M. Hilton, B. Mutlu-Pakdil, T. M. C. Abbott, M. Aguena, O. Alves, F. Andrade-Oliveira, J. Annis, D. Bacon, E. Bertin, S. Bocquet, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, F. J. Castander, C. Conselice, M. Costanzi, L. N. da Costa, M. E. S. Pereira , et al. (44 additional authors not shown)

Abstract: We characterise the properties and evolution of Bright Central Galaxies (BCGs) and the surrounding intracluster light (ICL) in galaxy clusters identified in overlapping regions of the Dark Energy Survey and Atacama Cosmology Telescope Survey (DES-ACT), covering the redshift range $0.20<z<0.80$. Using this sample, we measure no change in the ICL's stellar content (between 50-300\,kpc) over this red… ▽ More We characterise the properties and evolution of Bright Central Galaxies (BCGs) and the surrounding intracluster light (ICL) in galaxy clusters identified in overlapping regions of the Dark Energy Survey and Atacama Cosmology Telescope Survey (DES-ACT), covering the redshift range $0.20<z<0.80$. Using this sample, we measure no change in the ICL's stellar content (between 50-300\,kpc) over this redshift range in clusters with log$_{10}(M_{\rm 200m,SZ}$/M$_{\odot})>$14.4. We also measure the stellar mass - halo mass (SMHM) relation for the BCG+ICL system and find that the slope, $β$, which characterises the dependence of $M_{\rm 200m,SZ}$ on the BCG+ICL stellar mass, increases with radius. The outskirts are more strongly correlated with the halo than the core, which supports that the BCG+ICL system follows a two-phase growth, where recent growth ($z<2$) occurs beyond the BCG's core. Additionally, we compare our observed SMHM relation results to the IllustrisTNG 300-1 cosmological hydrodynamic simulations and find moderate qualitative agreement in the amount of diffuse light. However, the SMHM relation's slope is steeper in TNG300-1 and the intrinsic scatter is lower, likely from the absence of projection effects in TNG300-1. Additionally, we find that the ICL exhibits a colour gradient such that the outskirts are bluer than the core. Moreover, for the lower halo mass clusters (log$_{10}(M_{\rm 200m,SZ}$/M$_{\odot})<$14.59 ), we detect a modest change in the colour gradient's slope with lookback time, which combined with the absence of stellar mass growth may suggest that lower mass clusters have been involved in growth via tidal stripping more recently than their higher mass counterparts. △ Less

Submitted 14 September, 2022; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 19 pages, 12 Figures, 3 Tables. Submitted to MNRAS on 9/7/2022 -- Metadata typo corrected

Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck III: Combined cosmological constraints

Authors: T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, J. Annis, B. Ansarinejad, S. Avila, D. Bacon, E. J. Baxter, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, E. Bertin, J. Blazek, L. E. Bleem, S. Bocquet, D. Brooks, E. Buckley-Geer, D. L. Burke, H. Camacho, A. Campos, J. E. Carlstrom , et al. (146 additional authors not shown)

Abstract: We present cosmological constraints from the analysis of two-point correlation functions between galaxy positions and galaxy lensing measured in Dark Energy Survey (DES) Year 3 data and measurements of cosmic microwave background (CMB) lensing from the South Pole Telescope (SPT) and Planck. When jointly analyzing the DES-only two-point functions and the DES cross-correlations with SPT+Planck CMB l… ▽ More We present cosmological constraints from the analysis of two-point correlation functions between galaxy positions and galaxy lensing measured in Dark Energy Survey (DES) Year 3 data and measurements of cosmic microwave background (CMB) lensing from the South Pole Telescope (SPT) and Planck. When jointly analyzing the DES-only two-point functions and the DES cross-correlations with SPT+Planck CMB lensing, we find $Ω_{\rm m} = 0.344\pm 0.030$ and $S_8 \equiv σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.773\pm 0.016$, assuming $Λ$CDM. When additionally combining with measurements of the CMB lensing autospectrum, we find $Ω_{\rm m} = 0.306^{+0.018}_{-0.021}$ and $S_8 = 0.792\pm 0.012$. The high signal-to-noise of the CMB lensing cross-correlations enables several powerful consistency tests of these results, including comparisons with constraints derived from cross-correlations only, and comparisons designed to test the robustness of the galaxy lensing and clustering measurements from DES. Applying these tests to our measurements, we find no evidence of significant biases in the baseline cosmological constraints from the DES-only analyses or from the joint analyses with CMB lensing cross-correlations. However, the CMB lensing cross-correlations suggest possible problems with the correlation function measurements using alternative lens galaxy samples, in particular the redMaGiC galaxies and high-redshift MagLim galaxies, consistent with the findings of previous studies. We use the CMB lensing cross-correlations to identify directions for further investigating these problems. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 20 pages, 15 figures

Report number: FERMILAB-PUB-22-475-PPD

Constraining the Baryonic Feedback with Cosmic Shear Using the DES Year-3 Small-Scale Measurements

Authors: A. Chen, G. Aricò, D. Huterer, R. Angulo, N. Weaverdyck, O. Friedrich, L. F. Secco, C. Hernández-Monteagudo, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, J. Blazek, A. Brandao-Souza, S. L. Bridle, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, C. Chang , et al. (117 additional authors not shown)

Abstract: We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package \texttt{Baccoemu} to accelerate the evaluation of the baryonic nonlinear matter power spectrum. We design our ana… ▽ More We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package \texttt{Baccoemu} to accelerate the evaluation of the baryonic nonlinear matter power spectrum. We design our analysis pipeline to focus on the constraints of the baryonic suppression effects, utilizing the implication given by a principal component analysis on the Fisher forecasts. Our constraint on the baryonic effects can then be used to better model and ameliorate the effects of baryons in producing cosmological constraints from the next generation large-scale structure surveys. We detect the baryonic suppression on the cosmic shear measurements with a $\sim 2 σ$ significance. The characteristic halo mass for which half of the gas is ejected by baryonic feedback is constrained to be $M_c > 10^{13.2} h^{-1} M_{\odot}$ (95\% C.L.). The best-fit baryonic suppression is $\sim 5\%$ at $k=1.0 {\rm Mpc}\ h^{-1}$ and $\sim 15\%$ at $k=5.0 {\rm Mpc} \ h^{-1}$. Our findings are robust with respect to the assumptions about the cosmological parameters, specifics of the baryonic model, and intrinsic alignments. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Comments: 20 pages, 10 figures. DES Collaboration, Year-3 analysis

Lightweight starshade position sensing with convolutional neural networks and simulation-based inference

Authors: Andrew Chen, Anthony Harness, Peter Melchior

Abstract: Starshades are a leading technology to enable the direct detection and spectroscopic characterization of Earth-like exoplanets. To keep the starshade and telescope aligned over large separations, reliable sensing of the peak of the diffracted light of the occluded star is required. Current techniques rely on image matching or model fitting, both of which put substantial computational burdens on re… ▽ More Starshades are a leading technology to enable the direct detection and spectroscopic characterization of Earth-like exoplanets. To keep the starshade and telescope aligned over large separations, reliable sensing of the peak of the diffracted light of the occluded star is required. Current techniques rely on image matching or model fitting, both of which put substantial computational burdens on resource-limited spacecraft computers. We present a lightweight image processing method based on a convolutional neural network paired with a simulation-based inference technique to estimate the position of the spot of Arago and its uncertainty. The method achieves an accuracy of a few centimeters across the entire pupil plane, while only requiring 1.6 MB in stored data structures and 5.3 MFLOPs (million floating point operations) per image at test time. By deploying our method at the Princeton Starshade Testbed, we demonstrate that the neural network can be trained on simulated images and used on real images, and that it can successfully be integrated in the control system for closed-loop formation flying. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: submitted to JATIS

Snowmass2021 Cosmic Frontier White Paper: Enabling Flagship Dark Energy Experiments to Reach their Full Potential

Authors: Jonathan A. Blazek, Doug Clowe, Thomas E. Collett, Ian P. Dell'Antonio, Mark Dickinson, Lluís Galbany, Eric Gawiser, Katrin Heitmann, Renée Hložek, Mustapha Ishak, Saurabh W. Jha, Alex G. Kim, C. Danielle Leonard, Anja von der Linden, Michelle Lochner, Rachel Mandelbaum, Peter Melchior, Joel Meyers, Jeffrey A. Newman, Peter Nugent, Saul Perlmutter, Daniel J. Perrefort, Javier Sánchez, Samuel J. Schmidt, Sukhdeep Singh , et al. (3 additional authors not shown)

Abstract: A new generation of powerful dark energy experiments will open new vistas for cosmology in the next decade. However, these projects cannot reach their utmost potential without data from other telescopes. This white paper focuses in particular on the compelling benefits of ground-based spectroscopic and photometric observations to complement the Vera C. Rubin Observatory, as well as smaller program… ▽ More A new generation of powerful dark energy experiments will open new vistas for cosmology in the next decade. However, these projects cannot reach their utmost potential without data from other telescopes. This white paper focuses in particular on the compelling benefits of ground-based spectroscopic and photometric observations to complement the Vera C. Rubin Observatory, as well as smaller programs in aid of a DESI-2 experiment and CMB-S4. These additional data sets will both improve dark energy constraints from these flagship projects beyond what would possible on their own and open completely new windows into fundamental physics. For example, additional photometry and single-object spectroscopy will provide necessary follow-up information for supernova and strong lensing cosmology, while highly-multiplexed spectroscopy both from smaller facilities over wide fields and from larger facilities over narrower regions of sky will yield more accurate photometric redshift estimates for weak lensing and galaxy clustering measurements from the Rubin Observatory, provide critical spectroscopic host galaxy redshifts for supernova Hubble diagrams, provide improved understanding of limiting astrophysical systematic effects, and enable new measurements that probe the nature of gravity. A common thread is that access to complementary data from a range of telescopes/instruments would have a substantial impact on the rate of advance of dark energy science in the coming years. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: White paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

The DECam Local Volume Exploration Survey Data Release 2

Authors: A. Drlica-Wagner, P. S. Ferguson, M. Adamów, M. Aguena, F. Andrade-Oliveira, D. Bacon, K. Bechtol, E. F. Bell, E. Bertin, P. Bilaji, S. Bocquet, C. R. Bom, D. Brooks, D. L. Burke, J. A. Carballo-Bello, J. L. Carlin, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, F. J. Castander, W. Cerny, C. Chang, Y. Choi, C. Conselice, M. Costanzi , et al. (99 additional authors not shown)

Abstract: We present the second public data release (DR2) from the DECam Local Volume Exploration survey (DELVE). DELVE DR2 combines new DECam observations with archival DECam data from the Dark Energy Survey, the DECam Legacy Survey, and other DECam community programs. DELVE DR2 consists of ~160,000 exposures that cover >21,000 deg^2 of the high Galactic latitude (|b| > 10 deg) sky in four broadband optica… ▽ More We present the second public data release (DR2) from the DECam Local Volume Exploration survey (DELVE). DELVE DR2 combines new DECam observations with archival DECam data from the Dark Energy Survey, the DECam Legacy Survey, and other DECam community programs. DELVE DR2 consists of ~160,000 exposures that cover >21,000 deg^2 of the high Galactic latitude (|b| > 10 deg) sky in four broadband optical/near-infrared filters (g, r, i, z). DELVE DR2 provides point-source and automatic aperture photometry for ~2.5 billion astronomical sources with a median 5σ point-source depth of g=24.3, r=23.9, i=23.5, and z=22.8 mag. A region of ~17,000 deg^2 has been imaged in all four filters, providing four-band photometric measurements for ~618 million astronomical sources. DELVE DR2 covers more than four times the area of the previous DELVE data release and contains roughly five times as many astronomical objects. DELVE DR2 is publicly available via the NOIRLab Astro Data Lab science platform. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 23 pages, 10 figures, 5 tables; to be submitted to AAS Journals; public data release at https://datalab.noirlab.edu/delve/. arXiv admin note: substantial text overlap with arXiv:2103.07476

Report number: FERMILAB-PUB-22-209-LDRD-PPD

Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck II: Cross-correlation measurements and cosmological constraints

Authors: C. Chang, Y. Omori, E. J. Baxter, C. Doux, A. Choi, S. Pandey, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, K. Bechtol, M. R. Becker, G. M. Bernstein, F. Bianchini, J. Blazek, L. E. Bleem, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, R. Chen, J. Cordero, T. M. Crawford, M. Crocce , et al. (141 additional authors not shown)

Abstract: Cross-correlations of galaxy positions and galaxy shears with maps of gravitational lensing of the cosmic microwave background (CMB) are sensitive to the distribution of large-scale structure in the Universe. Such cross-correlations are also expected to be immune to some of the systematic effects that complicate correlation measurements internal to galaxy surveys. We present measurements and model… ▽ More Cross-correlations of galaxy positions and galaxy shears with maps of gravitational lensing of the cosmic microwave background (CMB) are sensitive to the distribution of large-scale structure in the Universe. Such cross-correlations are also expected to be immune to some of the systematic effects that complicate correlation measurements internal to galaxy surveys. We present measurements and modeling of the cross-correlations between galaxy positions and galaxy lensing measured in the first three years of data from the Dark Energy Survey with CMB lensing maps derived from a combination of data from the 2500 deg$^2$ SPT-SZ survey conducted with the South Pole Telescope and full-sky data from the Planck satellite. The CMB lensing maps used in this analysis have been constructed in a way that minimizes biases from the thermal Sunyaev Zel'dovich effect, making them well suited for cross-correlation studies. The total signal-to-noise of the cross-correlation measurements is 23.9 (25.7) when using a choice of angular scales optimized for a linear (nonlinear) galaxy bias model. We use the cross-correlation measurements to obtain constraints on cosmological parameters. For our fiducial galaxy sample, which consist of four bins of magnitude-selected galaxies, we find constraints of $Ω_{m} = 0.272^{+0.032}_{-0.052}$ and $S_{8} \equiv σ_8 \sqrt{Ω_{m}/0.3}= 0.736^{+0.032}_{-0.028}$ ($Ω_{m} = 0.245^{+0.026}_{-0.044}$ and $S_{8} = 0.734^{+0.035}_{-0.028}$) when assuming linear (nonlinear) galaxy bias in our modeling. Considering only the cross-correlation of galaxy shear with CMB lensing, we find $Ω_{m} = 0.270^{+0.043}_{-0.061}$ and $S_{8} = 0.740^{+0.034}_{-0.029}$. Our constraints on $S_8$ are consistent with recent cosmic shear measurements, but lower than the values preferred by primary CMB measurements from Planck. △ Less

Submitted 31 March, 2022; v1 submitted 23 March, 2022; originally announced March 2022.

Comments: 25 pages, 19 figures, submitted to PRD

Report number: FERMILAB-PUB-22-098-PPD

Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck I: Construction of CMB Lensing Maps and Modeling Choices

Authors: Y. Omori, E. J. Baxter, C. Chang, O. Friedrich, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, K. Bechtol, M. R. Becker, G. M. Bernstein, J. Blazek, L. E. Bleem, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, R. Chen, A. Choi, J. Cordero, T. M. Crawford, M. Crocce, C. Davis, J. DeRose , et al. (138 additional authors not shown)

Abstract: Joint analyses of cross-correlations between measurements of galaxy positions, galaxy lensing, and lensing of the cosmic microwave background (CMB) offer powerful constraints on the large-scale structure of the Universe. In a forthcoming analysis, we will present cosmological constraints from the analysis of such cross-correlations measured using Year 3 data from the Dark Energy Survey (DES), and… ▽ More Joint analyses of cross-correlations between measurements of galaxy positions, galaxy lensing, and lensing of the cosmic microwave background (CMB) offer powerful constraints on the large-scale structure of the Universe. In a forthcoming analysis, we will present cosmological constraints from the analysis of such cross-correlations measured using Year 3 data from the Dark Energy Survey (DES), and CMB data from the South Pole Telescope (SPT) and Planck. Here we present two key ingredients of this analysis: (1) an improved CMB lensing map in the SPT-SZ survey footprint, and (2) the analysis methodology that will be used to extract cosmological information from the cross-correlation measurements. Relative to previous lensing maps made from the same CMB observations, we have implemented techniques to remove contamination from the thermal Sunyaev Zel'dovich effect, enabling the extraction of cosmological information from smaller angular scales of the cross-correlation measurements than in previous analyses with DES Year 1 data. We describe our model for the cross-correlations between these maps and DES data, and validate our modeling choices to demonstrate the robustness of our analysis. We then forecast the expected cosmological constraints from the galaxy survey-CMB lensing auto and cross-correlations. We find that the galaxy-CMB lensing and galaxy shear-CMB lensing correlations will on their own provide a constraint on $S_8=σ_8 \sqrt{Ω_{\rm m}/0.3}$ at the few percent level, providing a powerful consistency check for the DES-only constraints. We explore scenarios where external priors on shear calibration are removed, finding that the joint analysis of CMB lensing cross-correlations can provide constraints on the shear calibration amplitude at the 5 to 10% level. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 30 pages, 20 figures, To be submitted to PRD

Report number: FERMILAB-PUB-22-194-PPD

Dark Energy Survey Year 3 results: imprints of cosmic voids and superclusters in the Planck CMB lensing map

Authors: A. Kovács, P. Vielzeuf, I. Ferrero, P. Fosalba, U. Demirbozan, R. Miquel, C. Chang, N. Hamaus, G. Pollina, K. Bechtol, M. Becker, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, M. Crocce, A. Drlica-Wagner, J. Elvin-Poole, M. Gatti, G. Giannini, R. A. Gruendl, A. Porredon, A. J. Ross, E. S. Rykoff, I. Sevilla-Noarbe, E. Sheldon , et al. (60 additional authors not shown)

Abstract: The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set ($\sim$5,000 deg$^{2}$), expanding on previous measurements that used Y1 catalogues ($\sim$1,300 deg$^{2}$). Given the increased statistical power… ▽ More The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set ($\sim$5,000 deg$^{2}$), expanding on previous measurements that used Y1 catalogues ($\sim$1,300 deg$^{2}$). Given the increased statistical power compared to Y1 data, we report a $6.6σ$ detection of negative CMB convergence ($κ$) imprints using approximately 3,600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the $Λ$CDM model (parameters $Ω_{\rm m} = 0.25$, $σ_8 = 0.8$), and the discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_κ=κ_{\rm DES}/κ_{\rm MICE}$ to a simulation-based template with fixed shape and found a moderate $2σ$ deviation in the signal with $A_κ\approx0.79\pm0.12$. We also examined the WebSky simulation that is based on a Planck 2018 $Λ$CDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the $8.4σ$ level; again with a lower-than-expected $A_κ=0.84\pm0.10$ amplitude. The combination of voids and superclusters yields a $10.3σ$ detection with an $A_κ=0.82\pm0.08$ constraint on the CMB lensing amplitude, thus the overall signal is $2.3σ$ weaker than expected from MICE. △ Less

Submitted 14 July, 2022; v1 submitted 21 March, 2022; originally announced March 2022.

Comments: 14 pages, 8 figures, accepted by MNRAS after minor corrections

Accelerated Bayesian SED Modeling using Amortized Neural Posterior Estimation

Authors: ChangHoon Hahn, Peter Melchior

Abstract: State-of-the-art spectral energy distribution (SED) analyses use a Bayesian framework to infer the physical properties of galaxies from observed photometry or spectra. They require sampling from a high-dimensional space of SED model parameters and take $>10-100$ CPU hours per galaxy, which renders them practically infeasible for analyzing the $billions$ of galaxies that will be observed by upcomin… ▽ More State-of-the-art spectral energy distribution (SED) analyses use a Bayesian framework to infer the physical properties of galaxies from observed photometry or spectra. They require sampling from a high-dimensional space of SED model parameters and take $>10-100$ CPU hours per galaxy, which renders them practically infeasible for analyzing the $billions$ of galaxies that will be observed by upcoming galaxy surveys ($e.g.$ DESI, PFS, Rubin, Webb, and Roman). In this work, we present an alternative scalable approach to rigorous Bayesian inference using Amortized Neural Posterior Estimation (ANPE). ANPE is a simulation-based inference method that employs neural networks to estimate the posterior probability distribution over the full range of observations. Once trained, it requires no additional model evaluations to estimate the posterior. We present, and publicly release, ${\rm SED}{flow}$, an ANPE method to produce posteriors of the recent Hahn et al. (2022) SED model from optical photometry. ${\rm SED}{flow}$ takes ${\sim}1$ $second~per~galaxy$ to obtain the posterior distributions of 12 model parameters, all of which are in excellent agreement with traditional Markov Chain Monte Carlo sampling results. We also apply ${\rm SED}{flow}$ to 33,884 galaxies in the NASA-Sloan Atlas and publicly release their posteriors: see https://changhoonhahn.github.io/SEDflow. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 21 pages, 5 figures; submitted to ApJ; code available at https://changhoonhahn.github.io/SEDflow

Robust sampling for weak lensing and clustering analyses with the Dark Energy Survey

Authors: P. Lemos, N. Weaverdyck, R. P. Rollins, J. Muir, A. Ferté, A. R. Liddle, A. Campos, D. Huterer, M. Raveri, J. Zuntz, E. Di Valentino, X. Fang, W. G. Hartley, M. Aguena, S. Allam, J. Annis, E. Bertin, S. Bocquet, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, F. J. Castander, A. Choi , et al. (46 additional authors not shown)

Abstract: Recent cosmological analyses rely on the ability to accurately sample from high-dimensional posterior distributions. A variety of algorithms have been applied in the field, but justification of the particular sampler choice and settings is often lacking. Here we investigate three such samplers to motivate and validate the algorithm and settings used for the Dark Energy Survey (DES) analyses of the… ▽ More Recent cosmological analyses rely on the ability to accurately sample from high-dimensional posterior distributions. A variety of algorithms have been applied in the field, but justification of the particular sampler choice and settings is often lacking. Here we investigate three such samplers to motivate and validate the algorithm and settings used for the Dark Energy Survey (DES) analyses of the first 3 years (Y3) of data from combined measurements of weak lensing and galaxy clustering. We employ the full DES Year 1 likelihood alongside a much faster approximate likelihood, which enables us to assess the outcomes from each sampler choice and demonstrate the robustness of our full results. We find that the ellipsoidal nested sampling algorithm $\texttt{MultiNest}$ reports inconsistent estimates of the Bayesian evidence and somewhat narrower parameter credible intervals than the sliced nested sampling implemented in $\texttt{PolyChord}$. We compare the findings from $\texttt{MultiNest}$ and $\texttt{PolyChord}$ with parameter inference from the Metropolis-Hastings algorithm, finding good agreement. We determine that $\texttt{PolyChord}$ provides a good balance of speed and robustness, and recommend different settings for testing purposes and final chains for analyses with DES Y3 data. Our methodology can readily be reproduced to obtain suitable sampler settings for future surveys. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: 16 pages, 11 figures

The Evolution of AGN Activity in Brightest Cluster Galaxies

Authors: T. Somboonpanyakul, M. McDonald, A. Noble, M. Aguena, S. Allam, A. Amon, F. Andrade-Oliveira, D. Bacon, M. B. Bayliss, E. Bertin, S. Bhargava, D. Brooks, E. Buckley-Geer, D. L. Burke, M. Calzadilla, R. Canning, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, M. Costanzi L. N. da Costa, M. E. S. Pereira J. De Vicente P. Doel P. Eisenhardt S. Everett A. E. Evrard, I. Ferrero, B. Flaugher, B. Floyd, J. García-Bellido , et al. (51 additional authors not shown)

Abstract: We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations on the full 2500 deg^2 South Pole Telescope (SPT)-SZ cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-infrared color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-host… ▽ More We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations on the full 2500 deg^2 South Pole Telescope (SPT)-SZ cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-infrared color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-hosting BCG fraction, which is defined as the fraction of BCGs hosting bright central AGNs over all possible BCGs. Assuming {\bf an evolving} single-burst stellar population model, we find statistically significant evidence (>99.9%) for a mid-IR excess at high redshift compared to low redshift, suggesting that the fraction of AGN-hosting BCGs increases with redshift over the range of 0 < z < 1.3. The best-fit redshift trend of the AGN-hosting BCG fraction has the form (1+z)^(4.1+/-1.0). These results are consistent with previous studies in galaxy clusters as well as field galaxies. One way to explain this result is that member galaxies at high redshift tend to have more cold gas. While BCGs in nearby galaxy clusters grow mostly by dry mergers with cluster members, leading to no increase in AGN activity, BCGs at high redshift could primarily merge with gas-rich satellites, providing fuel for feeding AGNs. If this observed increase in AGN activity is linked to gas-rich mergers, rather than ICM cooling, we would expect to see an increase in scatter in the P_cav vs L_cool relation at z > 1. Lastly, this work confirms that the runaway cooling phase, as predicted by the classical cooling flow model, in the Phoenix cluster is extremely rare and most BCGs have low (relative to Eddington) black hole accretion rates. △ Less

Submitted 9 February, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 13 pages, 9 figures, Accepted for publication in ApJ

Rubin-Euclid Derived Data Products: Initial Recommendations

Authors: Leanne P. Guy, Jean-Charles Cuillandre, Etienne Bachelet, Manda Banerji, Franz E. Bauer, Thomas Collett, Christopher J. Conselice, Siegfried Eggl, Annette Ferguson, Adriano Fontana, Catherine Heymans, Isobel M. Hook, Éric Aubourg, Hervé Aussel, James Bosch, Benoit Carry, Henk Hoekstra, Konrad Kuijken, Francois Lanusse, Peter Melchior, Joseph Mohr, Michele Moresco, Reiko Nakajima, Stéphane Paltani, Michael Troxel , et al. (95 additional authors not shown)

Abstract: This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum… ▽ More This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology. △ Less

Submitted 13 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

The DES view of the Eridanus supervoid and the CMB Cold Spot

Authors: A. Kovács, N. Jeffrey, M. Gatti, C. Chang, L. Whiteway, N. Hamaus, O. Lahav, G. Pollina, D. Bacon, T. Kacprzak, B. Mawdsley, S. Nadathur, D. Zeurcher, J. García-Bellido, A. Alarcon, A. Amon, K. Bechtol, G. M. Bernstein, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, R. Chen, A. Choi, J. Cordero , et al. (97 additional authors not shown)

Abstract: The Cold Spot is a puzzling large-scale feature in the Cosmic Microwave Background temperature maps and its origin has been subject to active debate. As an important foreground structure at low redshift, the Eridanus supervoid was recently detected, but it was subsequently determined that, assuming the standard $Λ$CDM model, only about 10-20$\%$ of the observed temperature depression can be accoun… ▽ More The Cold Spot is a puzzling large-scale feature in the Cosmic Microwave Background temperature maps and its origin has been subject to active debate. As an important foreground structure at low redshift, the Eridanus supervoid was recently detected, but it was subsequently determined that, assuming the standard $Λ$CDM model, only about 10-20$\%$ of the observed temperature depression can be accounted for via its Integrated Sachs-Wolfe imprint. However, $R\gtrsim100~h^{-1}\mathrm{Mpc}$ supervoids elsewhere in the sky have shown ISW imprints $A_{\mathrm{ISW}}\approx5.2\pm1.6$ times stronger than expected from $Λ$CDM ($A_{\mathrm{ISW}}=1$), which warrants further inspection. Using the Year-3 redMaGiC catalogue of luminous red galaxies from the Dark Energy Survey, here we confirm the detection of the Eridanus supervoid as a significant under-density in the Cold Spot's direction at $z<0.2$. We also show, with $\mathrm{S/N}\gtrsim5$ significance, that the Eridanus supervoid appears as the most prominent large-scale under-density in the dark matter mass maps that we reconstructed from DES Year-3 gravitational lensing data. While we report no significant anomalies, an interesting aspect is that the amplitude of the lensing signal from the Eridanus supervoid at the Cold Spot centre is about $30\%$ lower than expected from similar peaks found in N-body simulations based on the standard $Λ$CDM model with parameters $Ω_{\rm m} = 0.279$ and $σ_8 = 0.82$. Overall, our results confirm the causal relation between these individually rare structures in the cosmic web and in the CMB, motivating more detailed future surveys in the Cold Spot region. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: accepted for publication by MNRAS, 14 pages, 10 figures

Lensing Without Borders. I. A Blind Comparison of the Amplitude of Galaxy-Galaxy Lensing Between Independent Imaging Surveys

Authors: A. Leauthaud, A. Amon, S. Singh, D. Gruen, J. U. Lange, S. Huang, N. C. Robertson, T. N. Varga, Y. Luo, C. Heymans, H. Hildebrandt, C. Blake, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Annis, E. Bertin, S. Bhargava, J. Blazek, S. L. Bridle, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero , et al. (82 additional authors not shown)

Abstract: Lensing Without Borders is a cross-survey collaboration created to assess the consistency of galaxy-galaxy lensing signals ($ΔΣ$) across different data-sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of $ΔΣ$ using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported… ▽ More Lensing Without Borders is a cross-survey collaboration created to assess the consistency of galaxy-galaxy lensing signals ($ΔΣ$) across different data-sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of $ΔΣ$ using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported systematic errors which agree to better than 2.3$σ$ in four lens bins and three radial ranges. For lenses with $z_{\rm L}>0.43$ and considering statistical errors, we detect a 3-4$σ$ correlation between lensing amplitude and survey depth. This correlation could arise from the increasing impact at higher redshift of unrecognised galaxy blends on shear calibration and imperfections in photometric redshift calibration. At $z_{\rm L}>0.54$ amplitudes may additionally correlate with foreground stellar density. The amplitude of these trends is within survey-defined systematic error budgets which are designed to include known shear and redshift calibration uncertainty. Using a fully empirical and conservative method, we do not find evidence for large unknown systematics. Systematic errors greater than 15% (25%) ruled out in three lens bins at 68% (95%) confidence at $z<0.54$. Differences with respect to predictions based on clustering are observed to be at the 20-30% level. Our results therefore suggest that lensing systematics alone are unlikely to fully explain the "lensing is low" effect at $z<0.54$. This analysis demonstrates the power of cross-survey comparisons and provides a promising path for identifying and reducing systematics in future lensing analyses. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: 41 page, 20 figures

Cosmic Shear in Harmonic Space from the Dark Energy Survey Year 1 Data: Compatibility with Configuration Space Results

Authors: H. Camacho, F. Andrade-Oliveira, A. Troja, R. Rosenfeld, L. Faga, R. Gomes, C. Doux, X. Fang, M. Lima, V. Miranda, T. F. Eifler, O. Friedrich, M. Gatti, G. M. Bernstein, J. Blazek, S. L. Bridle, A. Choi, C. Davis, J. DeRose, E. Gaztanaga, D. Gruen, W. G. Hartley, B. Hoyle, M. Jarvis, N. MacCrann , et al. (74 additional authors not shown)

Abstract: We perform a cosmic shear analysis in harmonic space using the first year of data collected by the Dark Energy Survey (DES-Y1). We measure the cosmic weak lensing shear power spectra using the Metacalibration catalogue and perform a likelihood analysis within the framework of CosmoSIS. We set scale cuts based on baryonic effects contamination and model redshift and shear calibration uncertainties… ▽ More We perform a cosmic shear analysis in harmonic space using the first year of data collected by the Dark Energy Survey (DES-Y1). We measure the cosmic weak lensing shear power spectra using the Metacalibration catalogue and perform a likelihood analysis within the framework of CosmoSIS. We set scale cuts based on baryonic effects contamination and model redshift and shear calibration uncertainties as well as intrinsic alignments. We adopt as fiducial covariance matrix an analytical computation accounting for the mask geometry in the Gaussian term, including non-Gaussian contributions. A suite of 1200 lognormal simulations is used to validate the harmonic space pipeline and the covariance matrix. We perform a series of stress tests to gauge the robustness of the harmonic space analysis. Finally, we use the DES-Y1 pipeline in configuration space to perform a similar likelihood analysis and compare both results, demonstrating their compatibility in estimating the cosmological parameters $S_8$, $σ_8$ and $Ω_m$. The methods implemented and validated in this paper will allow us to perform a consistent harmonic space analysis in the upcoming DES data. △ Less

Submitted 10 October, 2022; v1 submitted 13 November, 2021; originally announced November 2021.

Comments: 16 pages, 11 figures. This version matches the published one

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 516, Issue 4, November 2022, Page 5799

Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps

Authors: M. Gatti, B. Jain, C. Chang, M. Raveri, D. Zürcher, L. Secco, L. Whiteway, N. Jeffrey, C. Doux, T. Kacprzak, D. Bacon, P. Fosalba, A. Alarcon, A. Amon, K. Bechtol, M. Becker, G. Bernstein, J. Blazek, A. Campos, A. Choi, C. Davis, J. Derose, S. Dodelson, F. Elsner, J. Elvin-Poole , et al. (85 additional authors not shown)

Abstract: We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing… ▽ More We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analysed in the context of the $Λ$CDM model, varying 5 cosmological parameters and 19 nuisance parameters modelling astrophysical and measurement systematics. Our modelling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7\% measurement of the amplitude of fluctuations parameter $S_8\equiv σ_8 (Ω_m/0.3)^{0.5} = 0.784\pm 0.013$. The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the $Λ$CDM model. The addition of the third moment improves the constraints on $S_8$ and $Ω_{\rm m}$ by $\sim$15\% and $\sim$25\% compared to an analysis that only uses second moments. We compare our results with {\it Planck} constraints from the Cosmic Microwave Background (CMB), finding a $2.2$ \textendash $2.8σ$ tension in the full parameter space, depending on the combination of moments considered. The third moment independently is in $2.8σ$ tension with {\it Planck}, and thus provides a cross-check on analyses of 2-point correlations. △ Less

Submitted 9 September, 2022; v1 submitted 19 October, 2021; originally announced October 2021.

Comments: 27 pages, 20 figures, accepted for publication in PRD

Dark Energy Survey Year 3 results: Cosmology with peaks using an emulator approach

Authors: D. Zürcher, J. Fluri, R. Sgier, T. Kacprzak, M. Gatti, C. Doux, L. Whiteway, A. Refregier, C. Chang, N. Jeffrey, B. Jain, P. Lemos, D. Bacon, A. Alarcon, A. Amon, K. Bechtol, M. Becker, G. Bernstein, A. Campos, R. Chen, A. Choi, C. Davis, J. Derose, S. Dodelson, F. Elsner , et al. (97 additional authors not shown)

Abstract: We constrain the matter density $Ω_{\mathrm{m}}$ and the amplitude of density fluctuations $σ_8$ within the $Λ$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps create… ▽ More We constrain the matter density $Ω_{\mathrm{m}}$ and the amplitude of density fluctuations $σ_8$ within the $Λ$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv~σ_8\sqrt{Ω_{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68% confidence limits), with a precision of 1.8%, an improvement of ~38% compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in $S_8$ is recorded. We find a mild tension of $1.5 \thinspace σ$ between our study and the results from Planck 2018, with our analysis yielding a lower $S_8$. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment $A_{\mathrm{IA}}$ and $S_8$, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data. △ Less

Submitted 21 October, 2021; v1 submitted 19 October, 2021; originally announced October 2021.

The DES Bright Arcs Survey: Candidate Strongly Lensed Galaxy Systems from the Dark Energy Survey 5,000 Sq. Deg. Footprint

Authors: J. H. O'Donnell, R. D. Wilkinson, H. T. Diehl, C. Aros-Bunster, K. Bechtol, S. Birrer, E. J. Buckley-Geer, A. Carnero Rosell, M. Carrasco Kind, L. N. da Costa, S. J. Gonzalez Lozano, R. A. Gruendl, M. Hilton, H. Lin, K. A. Lindgren, J. Martin, A. Pieres, E. S. Rykoff, I. Sevilla-Noarbe, E. Sheldon, C. Sifón, D. L. Tucker, B. Yanny, T. M. C. Abbott, M. Aguena , et al. (57 additional authors not shown)

Abstract: We report the combined results of eight searches for strong gravitational lens systems in the full 5,000 sq. deg. of Dark Energy Survey (DES) observations. The observations accumulated by the end of the third observing season fully covered the DES footprint in 5 filters (grizY), with an $i-$band limiting magnitude (at $10σ$) of 23.44. In four searches, a list of potential candidates was identified… ▽ More We report the combined results of eight searches for strong gravitational lens systems in the full 5,000 sq. deg. of Dark Energy Survey (DES) observations. The observations accumulated by the end of the third observing season fully covered the DES footprint in 5 filters (grizY), with an $i-$band limiting magnitude (at $10σ$) of 23.44. In four searches, a list of potential candidates was identified using a color and magnitude selection from the object catalogs created from the first three observing seasons. Three other searches were conducted at the locations of previously identified galaxy clusters. Cutout images of potential candidates were then visually scanned using an object viewer. An additional set of candidates came from a data-quality check of a subset of the color-coadd "tiles" created from the full DES six-season data set. A short list of the most promising strong lens candidates was then numerically ranked according to whether or not we judged them to be bona fide strong gravitational lens systems. These searches discovered a diverse set of 247 strong lens candidate systems, of which 81 are identified for the first time. We provide the coordinates, magnitudes, and photometric properties of the lens and source objects, and an estimate of the Einstein radius for 81 new systems and 166 previously reported. This catalog will be of use for selecting interesting systems for detailed follow-up, studies of galaxy cluster and group mass profiles, as well as a training/validation set for automated strong lens searches. △ Less

Submitted 3 January, 2022; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 38 pages, 17 figures, 4 tables, accepted by ApJS

Graph Neural Network-based Resource Allocation Strategies for Multi-Object Spectroscopy

Authors: Tianshu Wang, Peter Melchior

Abstract: Resource allocation problems are often approached with linear programming techniques. But many concrete allocation problems in the experimental and observational sciences cannot or should not be expressed in the form of linear objective functions. Even if the objective is linear, its parameters may not be known beforehand because they depend on the results of the experiment for which the allocatio… ▽ More Resource allocation problems are often approached with linear programming techniques. But many concrete allocation problems in the experimental and observational sciences cannot or should not be expressed in the form of linear objective functions. Even if the objective is linear, its parameters may not be known beforehand because they depend on the results of the experiment for which the allocation is to be determined. To address these challenges, we present a bipartite Graph Neural Network architecture for trainable resource allocation strategies. Items of value and constraints form the two sets of graph nodes, which are connected by edges corresponding to possible allocations. The GNN is trained on simulations or past problem occurrences to maximize any user-supplied, scientifically motivated objective function, augmented by an infeasibility penalty. The amount of feasibility violation can be tuned in relation to any available slack in the system. We apply this method to optimize the astronomical target selection strategy for the highly multiplexed Subaru Prime Focus Spectrograph instrument, where it shows superior results to direct gradient descent optimization and extends the capabilities of the currently employed solver which uses linear objective functions. The development of this method enables fast adjustment and deployment of allocation strategies, statistical analyses of allocation patterns, and fully differentiable, science-driven solutions for resource allocation problems. △ Less

Submitted 29 September, 2021; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: The GNN code used in this paper is available at https://github.com/tianshu-wang/PFS-GNN-bipartite

Cross-correlation of DES Y3 lensing and ACT/${\it Planck}$ thermal Sunyaev Zel'dovich Effect II: Modeling and constraints on halo pressure profiles

Authors: S. Pandey, M. Gatti, E. Baxter, J. C. Hill, X. Fang, C. Doux, G. Giannini, M. Raveri, J. DeRose, H. Huang, E. Moser, N. Battaglia, A. Alarcon, A. Amon, M. Becker, A. Campos, C. Chang, R. Chen, A. Choi, K. Eckert, J. Elvin-Poole, S. Everett, A. Ferte, I. Harrison, N. Maccrann , et al. (100 additional authors not shown)

Abstract: Hot, ionized gas leaves an imprint on the cosmic microwave background via the thermal Sunyaev Zel'dovich (tSZ) effect. The cross-correlation of gravitational lensing (which traces the projected mass) with the tSZ effect (which traces the projected gas pressure) is a powerful probe of the thermal state of ionized baryons throughout the Universe, and is sensitive to effects such as baryonic feedback… ▽ More Hot, ionized gas leaves an imprint on the cosmic microwave background via the thermal Sunyaev Zel'dovich (tSZ) effect. The cross-correlation of gravitational lensing (which traces the projected mass) with the tSZ effect (which traces the projected gas pressure) is a powerful probe of the thermal state of ionized baryons throughout the Universe, and is sensitive to effects such as baryonic feedback. In a companion paper (Gatti et al. 2021), we present tomographic measurements and validation tests of the cross-correlation between galaxy shear measurements from the first three years of observations of the Dark Energy Survey, and tSZ measurements from a combination of Atacama Cosmology Telescope and ${\it Planck}$ observations. In this work, we use the same measurements to constrain models for the pressure profiles of halos across a wide range of halo mass and redshift. We find evidence for reduced pressure in low mass halos, consistent with predictions for the effects of feedback from active galactic nuclei. We infer the hydrostatic mass bias ($B \equiv M_{500c}/M_{\rm SZ}$) from our measurements, finding $B = 1.8\pm0.1$ when adopting the ${\it Planck}$-preferred cosmological parameters. We additionally find that our measurements are consistent with a non-zero redshift evolution of $B$, with the correct sign and sufficient magnitude to explain the mass bias necessary to reconcile cluster count measurements with the ${\it Planck}$-preferred cosmology. Our analysis introduces a model for the impact of intrinsic alignments (IA) of galaxy shapes on the shear-tSZ correlation. We show that IA can have a significant impact on these correlations at current noise levels. △ Less

Submitted 24 November, 2022; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: 22 pages, 13 figures. Comments welcome

Cross-correlation of DES Y3 lensing and ACT/${\it Planck}$ thermal Sunyaev Zel'dovich Effect I: Measurements, systematics tests, and feedback model constraints

Authors: M. Gatti, S. Pandey, E. Baxter, J. C. Hill, E. Moser, M. Raveri, X. Fang, J. DeRose, G. Giannini, C. Doux, H. Huang, N. Battaglia, A. Alarcon, A. Amon, M. Becker, A. Campos, C. Chang, R. Chen, A. Choi, K. Eckert, J. Elvin-Poole, S. Everett, A. Ferte, I. Harrison, N. Maccrann , et al. (104 additional authors not shown)

Abstract: We present a tomographic measurement of the cross-correlation between thermal Sunyaev-Zeldovich (tSZ) maps from ${\it Planck}$ and the Atacama Cosmology Telescope (ACT) and weak galaxy lensing shears measured during the first three years of observations of the Dark Energy Survey (DES Y3). This correlation is sensitive to the thermal energy in baryons over a wide redshift range, and is therefore a… ▽ More We present a tomographic measurement of the cross-correlation between thermal Sunyaev-Zeldovich (tSZ) maps from ${\it Planck}$ and the Atacama Cosmology Telescope (ACT) and weak galaxy lensing shears measured during the first three years of observations of the Dark Energy Survey (DES Y3). This correlation is sensitive to the thermal energy in baryons over a wide redshift range, and is therefore a powerful probe of astrophysical feedback. We detect the correlation at a statistical significance of $21σ$, the highest significance to date. We examine the tSZ maps for potential contaminants, including cosmic infrared background (CIB) and radio sources, finding that CIB has a substantial impact on our measurements and must be taken into account in our analysis. We use the cross-correlation measurements to test different feedback models. In particular, we model the tSZ using several different pressure profile models calibrated against hydrodynamical simulations. Our analysis marginalises over redshift uncertainties, shear calibration biases, and intrinsic alignment effects. We also marginalise over $Ω_{\rm m}$ and $σ_8$ using ${\it Planck}$ or DES priors. We find that the data prefers the model with a low amplitude of the pressure profile at small scales, compatible with a scenario with strong AGN feedback and ejection of gas from the inner part of the halos. When using a more flexible model for the shear profile, constraints are weaker, and the data cannot discriminate between different baryonic prescriptions. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: submitted to PRD