-
Lens Modeling of STRIDES Strongly Lensed Quasars using Neural Posterior Estimation
Authors:
Sydney Erickson,
Sebastian Wagner-Carena,
Phil Marshall,
Martin Millon,
Simon Birrer,
Aaron Roodman,
Thomas Schmidt,
Tommaso Treu,
Stefan Schuldt,
Anowar Shajib,
Padma Venkatraman,
The LSST Dark Energy Science Collaboration
Abstract:
Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory's Legacy Survey of Space and Time (LSST) and the Euclid Wide Survey, we re…
▽ More
Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory's Legacy Survey of Space and Time (LSST) and the Euclid Wide Survey, we require fast and standardizable modeling techniques. To address this need, we apply neural posterior estimation (NPE) for modeling galaxy-scale strongly lensed quasars from the Strong Lensing Insights into the Dark Energy Survey (STRIDES) sample. NPE brings two advantages: speed and the ability to implicitly marginalize over nuisance parameters. We extend this method by employing sequential NPE to increase precision of mass model posteriors. We then fold individual lens models into a hierarchical Bayesian inference to recover the population distribution of lens mass parameters, accounting for out-of-distribution shift. After verifying our method using simulated analogs of the STRIDES lens sample, we apply our method to 14 Hubble Space Telescope single-filter observations. We find the population mean of the power-law elliptical mass distribution slope, $γ_{\text{lens}}$, to be $\mathcal{M}_{γ_{\text{lens}}}=2.13 \pm 0.06$. Our result represents the first population-level constraint for these systems. This population-level inference from fully automated modeling is an important stepping stone towards cosmological inference with large samples of strongly lensed quasars.
△ Less
Submitted 13 October, 2024;
originally announced October 2024.
-
Simulation-Based Inference Benchmark for LSST Weak Lensing Cosmology
Authors:
Justine Zeghal,
Denise Lanzieri,
François Lanusse,
Alexandre Boucaud,
Gilles Louppe,
Eric Aubourg,
Adrian E. Bayer,
The LSST Dark Energy Science Collaboration
Abstract:
Standard cosmological analysis, which relies on two-point statistics, fails to extract the full information of the data. This limits our ability to constrain with precision cosmological parameters. Thus, recent years have seen a paradigm shift from analytical likelihood-based to simulation-based inference. However, such methods require a large number of costly simulations. We focus on full-field i…
▽ More
Standard cosmological analysis, which relies on two-point statistics, fails to extract the full information of the data. This limits our ability to constrain with precision cosmological parameters. Thus, recent years have seen a paradigm shift from analytical likelihood-based to simulation-based inference. However, such methods require a large number of costly simulations. We focus on full-field inference, considered the optimal form of inference. Our objective is to benchmark several ways of conducting full-field inference to gain insight into the number of simulations required for each method. We make a distinction between explicit and implicit full-field inference. Moreover, as it is crucial for explicit full-field inference to use a differentiable forward model, we aim to discuss the advantages of having this property for the implicit approach. We use the sbi_lens package which provides a fast and differentiable log-normal forward model. This forward model enables us to compare explicit and implicit full-field inference with and without gradient. The former is achieved by sampling the forward model through the No U-Turns sampler. The latter starts by compressing the data into sufficient statistics and uses the Neural Likelihood Estimation algorithm and the one augmented with gradient. We perform a full-field analysis on LSST Y10 like weak lensing simulated mass maps. We show that explicit and implicit full-field inference yield consistent constraints. Explicit inference requires 630 000 simulations with our particular sampler corresponding to 400 independent samples. Implicit inference requires a maximum of 101 000 simulations split into 100 000 simulations to build sufficient statistics (this number is not fine tuned) and 1 000 simulations to perform inference. Additionally, we show that our way of exploiting the gradients does not significantly help implicit inference.
△ Less
Submitted 26 September, 2024;
originally announced September 2024.
-
Impact of Large-Scale Structure Systematics on Cosmological Parameter Estimation
Authors:
Humna Awan,
Eric Gawiser,
Javier Sanchez,
Ignacio Sevilla-Noarbe,
the LSST Dark Energy Science Collaboration
Abstract:
Large near-future galaxy surveys offer sufficient statistical power to make our cosmology analyses data-driven, limited primarily by systematic errors. Understanding the impact of systematics is therefore critical. We perform an end-to-end analysis to investigate the impact of some of the systematics that affect large-scale structure studies by doing an inference analysis using simulated density m…
▽ More
Large near-future galaxy surveys offer sufficient statistical power to make our cosmology analyses data-driven, limited primarily by systematic errors. Understanding the impact of systematics is therefore critical. We perform an end-to-end analysis to investigate the impact of some of the systematics that affect large-scale structure studies by doing an inference analysis using simulated density maps with various systematics; these include systematics caused by photometric redshifts (photo-$z$s), Galactic dust, structure induced by the telescope observing strategy and observing conditions, and incomplete covariance matrices. Specifically, we consider the impacts of incorrect photo-$z$ distributions (photometric biases, scatter, outliers; spectroscopic calibration biases), dust map resolution, incorrect dust law, selecting none or only some contaminant templates for deprojection, and using a diagonal covariance matrix instead of a full one. We quantify the biases induced by these systematics on cosmological parameter estimation using tomographic galaxy angular power spectra, with a focus on identifying whether the maximum plausible level of each systematic has an adverse impact on the estimation of key cosmological parameters from a galaxy clustering analysis with Rubin Observatory Legacy Survey of Space and Time (LSST). We find photo-$z$ systematics to be the most pressing out of the systematics investigated, with spectroscopic calibration biases leading to the greatest adverse impact while helpfully being flagged by a high $χ^2$ value for the best fit model. Larger-than-expected photo-$z$ scatter, on the other hand, has a significant impact without necessarily indicating a poor fit. In contrast, in the analysis framework used in this work, biases from observational systematics and incomplete covariance matrices are comfortably subdominant.
△ Less
Submitted 21 September, 2024;
originally announced September 2024.
-
Dark Energy Survey: 2.1% measurement of the Baryon Acoustic Oscillation scale from the final dataset
Authors:
Juan Mena-Fernández,
Dark Energy Survey Collaboration
Abstract:
Here, we present the angular diameter distance measurement obtained from the measurement of the Baryonic Acoustic Oscillation (BAO) feature using the completed Dark Energy Survey (DES) data, summarizing the main results of [Phys. Rev. D 110, 063514] and [Phys. Rev. D 110, 063515]. We use a galaxy sample optimized for BAO science in the redshift range 0.6 < z < 1.2, with an effective redshift of…
▽ More
Here, we present the angular diameter distance measurement obtained from the measurement of the Baryonic Acoustic Oscillation (BAO) feature using the completed Dark Energy Survey (DES) data, summarizing the main results of [Phys. Rev. D 110, 063514] and [Phys. Rev. D 110, 063515]. We use a galaxy sample optimized for BAO science in the redshift range 0.6 < z < 1.2, with an effective redshift of $z_{\rm eff}$ = 0.85. Our consensus measurement constrains the ratio of the angular distance to the sound horizon scale to $D_M(z_{\rm eff})/r_d$ = 19.51 $\pm$ 0.41. This measurement is found to be 2.13$σ$ below the angular BAO scale predicted by Planck. To date, it represents the most precise measurement from purely photometric data, and the most precise from any Stage-III experiment at such high redshift. The analysis was performed blinded to the BAO position and is shown to be robust against analysis choices, data removal, redshift calibrations and observational systematics.
△ Less
Submitted 13 September, 2024;
originally announced September 2024.
-
Selective Dynamical Imaging of Interferometric Data
Authors:
Joseph Farah,
Peter Galison,
Kazunori Akiyama,
Katherine L. Bouman,
Geoffrey C. Bower,
Andrew Chael,
Antonio Fuentes,
José L. Gómez,
Mareki Honma,
Michael D. Johnson,
Yutaro Kofuji,
Daniel P. Marrone,
Kotaro Moriyama,
Ramesh Narayan,
Dominic W. Pesce,
Paul Tiede,
Maciek Wielgus,
Guang-Yao Zhao,
The Event Horizon Telescope Collaboration
Abstract:
Recent developments in very long baseline interferometry (VLBI) have made it possible for the Event Horizon Telescope (EHT) to resolve the innermost accretion flows of the largest supermassive black holes on the sky. The sparse nature of the EHT's $(u, v)$-coverage presents a challenge when attempting to resolve highly time-variable sources. We demonstrate that the changing (u, v)-coverage of the…
▽ More
Recent developments in very long baseline interferometry (VLBI) have made it possible for the Event Horizon Telescope (EHT) to resolve the innermost accretion flows of the largest supermassive black holes on the sky. The sparse nature of the EHT's $(u, v)$-coverage presents a challenge when attempting to resolve highly time-variable sources. We demonstrate that the changing (u, v)-coverage of the EHT can contain regions of time over the course of a single observation that facilitate dynamical imaging. These optimal time regions typically have projected baseline distributions that are approximately angularly isotropic and radially homogeneous. We derive a metric of coverage quality based on baseline isotropy and density that is capable of ranking array configurations by their ability to produce accurate dynamical reconstructions. We compare this metric to existing metrics in the literature and investigate their utility by performing dynamical reconstructions on synthetic data from simulated EHT observations of sources with simple orbital variability. We then use these results to make recommendations for imaging the 2017 EHT Sgr A* data set.
△ Less
Submitted 12 September, 2024;
originally announced September 2024.
-
The Blending ToolKit: A simulation framework for evaluation of galaxy detection and deblending
Authors:
Ismael Mendoza,
Andrii Torchylo,
Thomas Sainrat,
Axel Guinot,
Alexandre Boucaud,
Maxime Paillassa,
Camille Avestruz,
Prakruth Adari,
Eric Aubourg,
Biswajit Biswas,
James Buchanan,
Patricia Burchat,
Cyrille Doux,
Remy Joseph,
Sowmya Kamath,
Alex I. Malz,
Grant Merz,
Hironao Miyatake,
Cécile Roucelle,
Tianqing Zhang,
the LSST Dark Energy Science Collaboration
Abstract:
We present an open source Python library for simulating overlapping (i.e., blended) images of galaxies and performing self-consistent comparisons of detection and deblending algorithms based on a suite of metrics. The package, named Blending Toolkit (BTK), serves as a modular, flexible, easy-to-install, and simple-to-use interface for exploring and analyzing systematic effects related to blended g…
▽ More
We present an open source Python library for simulating overlapping (i.e., blended) images of galaxies and performing self-consistent comparisons of detection and deblending algorithms based on a suite of metrics. The package, named Blending Toolkit (BTK), serves as a modular, flexible, easy-to-install, and simple-to-use interface for exploring and analyzing systematic effects related to blended galaxies in cosmological surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time (LSST). BTK has three main components: (1) a set of modules that perform fast image simulations of blended galaxies, using the open source image simulation package GalSim; (2) a module that standardizes the inputs and outputs of existing deblending algorithms; (3) a library of deblending metrics commonly defined in the galaxy deblending literature. In combination, these modules allow researchers to explore the impacts of galaxy blending in cosmological surveys. Additionally, BTK provides researchers who are developing a new deblending algorithm a framework to evaluate algorithm performance and make principled comparisons with existing deblenders. BTK includes a suite of tutorials and comprehensive documentation. The source code is publicly available on GitHub at https://github.com/LSSTDESC/BlendingToolKit.
△ Less
Submitted 26 September, 2024; v1 submitted 10 September, 2024;
originally announced September 2024.
-
YOLO-CL cluster detection in the Rubin/LSST DC2 simulation
Authors:
Kirill Grishin,
Simona Mei,
Stephane Ilic,
Michel Aguena,
Dominique Boutigny,
Marie Paturel,
the LSST Dark Energy Science Collaboration
Abstract:
LSST will provide galaxy cluster catalogs up to z$\sim$1 that can be used to constrain cosmological models once their selection function is well-understood. We have applied the deep convolutional network YOLO for CLuster detection (YOLO-CL) to LSST simulations from the Dark Energy Science Collaboration Data Challenge 2 (DC2), and characterized the LSST YOLO-CL cluster selection function. We have t…
▽ More
LSST will provide galaxy cluster catalogs up to z$\sim$1 that can be used to constrain cosmological models once their selection function is well-understood. We have applied the deep convolutional network YOLO for CLuster detection (YOLO-CL) to LSST simulations from the Dark Energy Science Collaboration Data Challenge 2 (DC2), and characterized the LSST YOLO-CL cluster selection function. We have trained and validated the network on images from a hybrid sample of (1) clusters observed in the Sloan Digital Sky Survey and detected with the red-sequence Matched-filter Probabilistic Percolation, and (2) simulated DC2 dark matter haloes with masses $M_{200c} > 10^{14} M_{\odot}$. We quantify the completeness and purity of the YOLO-CL cluster catalog with respect to DC2 haloes with $M_{200c} > 10^{14} M_{\odot}$. The YOLO-CL cluster catalog is 100% and 94% complete for halo mass $M_{200c} > 10^{14.6} M_{\odot}$ at $0.2<z<0.8$, and $M_{200c} > 10^{14} M_{\odot}$ and redshift $z \lesssim 1$, respectively, with only 6% false positive detections. All the false positive detections are dark matter haloes with $ 10^{13.4} M_{\odot} \lesssim M_{200c} \lesssim 10^{14} M_{\odot}$. The YOLO-CL selection function is almost flat with respect to the halo mass at $0.2 \lesssim z \lesssim 0.9$. The overall performance of YOLO-CL is comparable or better than other cluster detection methods used for current and future optical and infrared surveys. YOLO-CL shows better completeness for low mass clusters when compared to current detections in surveys using the Sunyaev Zel'dovich effect, and detects clusters at higher redshifts than X-ray-based catalogs. The strong advantage of YOLO-CL over traditional galaxy cluster detection techniques is that it works directly on images and does not require photometric and photometric redshift catalogs, nor does it need to mask stellar sources and artifacts.
△ Less
Submitted 5 September, 2024;
originally announced September 2024.
-
Impact of survey spatial variability on galaxy redshift distributions and the cosmological $3\times2$-point statistics for the Rubin Legacy Survey of Space and Time (LSST)
Authors:
Qianjun Hang,
Benjamin Joachimi,
Eric Charles,
John Franklin Crenshaw,
Patricia Larsen,
Alex I. Malz,
Sam Schmidt,
Ziang Yan,
Tianqing Zhang,
the LSST Dark Energy Science Collaboration
Abstract:
We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry dataset degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-$z$ estimator, BPZ, and a…
▽ More
We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry dataset degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-$z$ estimator, BPZ, and a machine learning method, FlexZBoost. We select the Gold sample, defined as $i<25.3$ for 10 year LSST data, with an adjusted magnitude cut for each year and divide it into five tomographic redshift bins for the weak lensing lens and source samples. We quantify the change in the number of objects, mean redshift, and width of each tomographic bin as a function of the coadd $i$-band depth for 1-year (Y1), 3-year (Y3), and 5-year (Y5) data. In particular, Y3 and Y5 have large non-uniformity due to the rolling cadence of LSST, hence provide a worst-case scenario of the impact from non-uniformity. We find that these quantities typically increase with depth, and the variation can be $10-40\%$ at extreme depth values. Based on these results and using Y3 as an example, we propagate the variable depth effect to the weak lensing $3\times2$pt data vector in harmonic space. We find that galaxy clustering is most susceptible to variable depth, causing significant deviations at large scales if not corrected for, due to the depth-dependent number density variations. For galaxy-shear and shear-shear power spectra, we find little impact given the expected LSST Y3 noise.
△ Less
Submitted 4 September, 2024;
originally announced September 2024.
-
Modeling the 3-point correlation function of projected scalar fields on the sphere
Authors:
Abraham Arvizu,
Alejandro Aviles,
Juan Carlos Hidalgo,
Eladio Moreno,
Gustavo Niz,
Mario A. Rodriguez-Meza,
Sofía Samario,
The LSST Dark Energy Science Collaboration
Abstract:
One of the main obstacles for the signal extraction of the three point correlation function using photometric surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST), will be the prohibitive computation time required for dealing with a vast quantity of sources. Brute force algorithms, which naively scales as $\mathcal{O}(N^3)$ with the number of objects, can be further improv…
▽ More
One of the main obstacles for the signal extraction of the three point correlation function using photometric surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST), will be the prohibitive computation time required for dealing with a vast quantity of sources. Brute force algorithms, which naively scales as $\mathcal{O}(N^3)$ with the number of objects, can be further improved with tree methods but not enough to deal with large scale correlations of Rubin's data. However, a harmonic basis decomposition of these higher order statistics reduces the time dramatically, to scale as a two-point correlation function with the number of objects, so that the signal can be extracted in a reasonable amount of time. In this work, we aim to develop the framework to use these expansions within the Limber approximation for scalar (or spin-0) fields, such as galaxy counts, weak lensing convergence or aperture masses. We develop an estimator to extract the signal from catalogs and different phenomenological and theoretical models for its description. The latter includes halo model and standard perturbation theory, to which we add a simple effective field theory prescription based on the short range of non-locality of cosmic fields, significantly improving the agreement with simulated data. In parallel to the modeling of the signal, we develop a code that can efficiently calculate three points correlations of more than 200 million data points (a full sky simulation with Nside=4096) in $\sim$40 minutes on a single high-performance computing node, enabling a feasible analysis for the upcoming LSST data.
△ Less
Submitted 29 August, 2024;
originally announced August 2024.
-
MADNESS Deblender: Maximum A posteriori with Deep NEural networks for Source Separation
Authors:
Biswajit Biswas,
Eric Aubourg,
Alexandre Boucaud,
Axel Guinot,
Junpeng Lao,
Cécile Roucelle,
the LSST Dark Energy Science Collaboration
Abstract:
Due to the unprecedented depth of the upcoming ground-based Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory, approximately two-thirds of the galaxies are likely to be affected by blending - the overlap of physically separated galaxies in images. Thus, extracting reliable shapes and photometry from individual objects will be limited by our ability to correct blending and con…
▽ More
Due to the unprecedented depth of the upcoming ground-based Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory, approximately two-thirds of the galaxies are likely to be affected by blending - the overlap of physically separated galaxies in images. Thus, extracting reliable shapes and photometry from individual objects will be limited by our ability to correct blending and control any residual systematic effect. Deblending algorithms tackle this issue by reconstructing the isolated components from a blended scene, but the most commonly used algorithms often fail to model complex realistic galaxy morphologies.
As part of an effort to address this major challenge, we present MADNESS, which takes a data-driven approach and combines pixel-level multi-band information to learn complex priors for obtaining the maximum a posteriori solution of deblending. MADNESS is based on deep neural network architectures such as variational auto-encoders and normalizing flows. The variational auto-encoder reduces the high-dimensional pixel space into a lower-dimensional space, while the normalizing flow models a data-driven prior in this latent space.
Using a simulated test dataset with galaxy models for a 10-year LSST survey and a galaxy density ranging from 48 to 80 galaxies per arcmin2 we characterize the aperture-photometry g-r color, structural similarity index, and pixel cosine similarity of the galaxies reconstructed by MADNESS. We compare our results against state-of-the-art deblenders including scarlet. With the r-band of LSST as an example, we show that MADNESS performs better than in all the metrics. For instance, the average absolute value of relative flux residual in the r-band for MADNESS is approximately 29% lower than that of scarlet. The code is publicly available on GitHub.
△ Less
Submitted 27 August, 2024;
originally announced August 2024.
-
Analytical Weak Lensing Shear Inference for Precision Cosmology
Authors:
Xiangchong Li,
Rachel Mandelbaum,
The LSST Dark Energy Science Collaboration
Abstract:
Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the 'renoising' no…
▽ More
Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the 'renoising' noise bias correction approach introduced by Metacalibration. We implement this analytical noise bias correction within the AnaCal shear estimation framework and propose several enhancements to the noise bias correction algorithm. We evaluate the improved AnaCal using simulations designed to replicate Rubin LSST imaging data. These simulations feature semi-realistic galaxies and stars, complete with representative distributions of magnitudes and Galactic spatial density. We conduct tests under various observational challenges, including cosmic rays, defective CCD columns, bright star saturation, bleed trails, and spatially variable point spread functions. Our results indicate a multiplicative bias in weak lensing shear recovery of less than a few tenths of a percent, meeting LSST DESC requirements without requiring calibration from external image simulations. Additionally, our algorithm achieves rapid processing, handling one galaxy in less than a millisecond.
△ Less
Submitted 12 August, 2024;
originally announced August 2024.
-
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.
-
H.E.S.S. observations of the 2021 periastron passage of PSR B1259-63/LS 2883
Authors:
H. E. S. S. Collaboration,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
H. Ashkar,
M. Backes,
V. Barbosa Martins,
R. Batzofin,
Y. Becherini,
D. Berge,
K. Bernlöhr,
M. Böttcher,
C. Boisson,
J. Bolmont,
M. de Bony de Lavergne,
J. Borowska,
M. Bouyahiaoui,
R. Brose,
A. Brown,
F. Brun,
B. Bruno,
T. Bulik,
C. Burger-Scheidlin,
S. Caroff,
S. Casanova
, et al. (119 additional authors not shown)
Abstract:
PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ day…
▽ More
PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ days to $t_p+127$ days around the system's 2021 periastron passage. We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve -- although sparsely sampled in this time period -- we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index $Γ= 2.65 \pm 0.04_{\text{stat}}$ $\pm 0.04_{\text{sys}}$, a value consistent with the previous H.E.S.S. observations of the source. We report spectral variability with a difference of $ΔΓ= 0.56 ~\pm~ 0.18_{\text{stat}}$ $~\pm~0.10_{\text{sys}}$ at 95% c.l., between sub-periods of the 2021 dataset. We also find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after $t_p+25$ days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from $\sim t_p-23$ days to $\sim t_p+126$ days. This suggests that the GeV and TeV emission originate from different electron populations.
△ Less
Submitted 26 June, 2024;
originally announced June 2024.
-
Joint Modelling of Astrophysical Systematics for Cosmology with LSST
Authors:
Nikolina Šarčević,
C. Danielle Leonard,
Markus M. Rau,
LSST Dark Energy Science Collaboration
Abstract:
We present a novel framework for jointly modelling the weak lensing source galaxy redshift distribution and the intrinsic alignment of galaxies via a shared luminosity function (LF). Considering this framework within the context of a Rubin Observatory's Legacy Survey of Space and Time (LSST) Year 1 and Year 10 cosmic shear analysis, we first demonstrate the substantial impact of the LF on both sou…
▽ More
We present a novel framework for jointly modelling the weak lensing source galaxy redshift distribution and the intrinsic alignment of galaxies via a shared luminosity function (LF). Considering this framework within the context of a Rubin Observatory's Legacy Survey of Space and Time (LSST) Year 1 and Year 10 cosmic shear analysis, we first demonstrate the substantial impact of the LF on both source galaxy redshift distributions and the intrinsic alignment contamination. We establish how the individual parameters of a Schechter LF model impact the redshift distribution of a magnitude-limited sample, and we demonstrate the effect of marginalising over the LF parameters as incorporated in the intrinsic alignment modelling of a standard cosmic shear analysis set-up. We forecast the impact of our joint modelling framework on cosmological parameter constraints. Our preliminary results are promising, indicating that this framework can yield cosmological constraints consistent with those expected from standard analyses, enhanced by the flexibility of not fixing LF parameters. We plan to further validate these findings with comprehensive Markov chain Monte Carlo simulations to robustly quantify bias avoidance and underscore the framework's efficacy. Taking advantage of our forecasting results and the parameter degeneracies, we identify the specific impact of the shape of the LF of source galaxies on the cosmic shear data vector. We also discuss the potential of this method in providing a way to model generic selection functions in redshift distribution estimation, as well as its possibilities for extension to a 3x2pt analysis, particularly with respect to incorporating galaxy bias in this luminosity-function-based framework. Although we consider the context of LSST cosmic shear in this work, the proposed joint modelling framework is generically applicable to weak lensing surveys.
△ Less
Submitted 12 June, 2024; v1 submitted 5 June, 2024;
originally announced June 2024.
-
Probabilistic Forward Modeling of Galaxy Catalogs with Normalizing Flows
Authors:
John Franklin Crenshaw,
J. Bryce Kalmbach,
Alexander Gagliano,
Ziang Yan,
Andrew J. Connolly,
Alex I. Malz,
Samuel J. Schmidt,
The LSST Dark Energy Science Collaboration
Abstract:
Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo-z) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo-z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce P…
▽ More
Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo-z) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo-z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of "true" redshift posteriors that can be used for photo-z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo-z estimation. We discuss how PZFlow will be used to validate the photo-z estimation pipeline of the Dark Energy Science Collaboration (DESC), and the wider applicability of PZFlow for statistical modeling of any tabular data.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Astrometric Redshifts of Supernovae
Authors:
Jaemyoung Jason Lee,
Masao Sako,
Richard Kessler,
Alex I. Malz,
The LSST Dark Energy Science Collaboration
Abstract:
Differential Chromatic Refraction (DCR) is caused by the wavelength dependence of our atmosphere's refractive index, which shifts the apparent positions of stars and galaxies and distorts their shapes depending on their spectral energy distributions (SEDs). While this effect is typically mitigated and corrected for in imaging observations, we investigate how DCR can instead be used to our advantag…
▽ More
Differential Chromatic Refraction (DCR) is caused by the wavelength dependence of our atmosphere's refractive index, which shifts the apparent positions of stars and galaxies and distorts their shapes depending on their spectral energy distributions (SEDs). While this effect is typically mitigated and corrected for in imaging observations, we investigate how DCR can instead be used to our advantage to infer the redshifts of supernovae from multi-band, time-series imaging data. We simulate Type Ia supernovae (SNe Ia) in the proposed Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Deep Drilling Field (DDF), and evaluate astrometric redshifts. We find that the redshift accuracy improves dramatically with the statistical quality of the astrometric measurements as well as with the accuracy of the astrometric solution. For a conservative choice of a 5-mas systematic uncertainty floor, we find that our redshift estimation is accurate at $z < 0.6$. We then combine our astrometric redshifts with both host galaxy photometric redshifts and supernovae photometric (light-curve) redshifts and show that this considerably improves the overall redshift estimates. These astrometric redshifts will be valuable especially since Rubin will discover a vast number of supernovae for which we will not be able to obtain spectroscopic redshifts.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Broadband Multi-wavelength Properties of M87 during the 2018 EHT Campaign including a Very High Energy Flaring Episode
Authors:
J. C. Algaba,
M. Balokovic,
S. Chandra,
W. Y. Cheong,
Y. Z. Cui,
F. D'Ammando,
A. D. Falcone,
N. M. Ford,
M. Giroletti,
C. Goddi,
M. A. Gurwell,
K. Hada,
D. Haggard,
S. Jorstad,
A. Kaur,
T. Kawashima,
S. Kerby,
J. Y. Kim,
M. Kino,
E. V. Kravchenko,
S. S. Lee,
R. S. Lu,
S. Markoff,
J. Michail,
J. Neilsen
, et al. (14 additional authors not shown)
Abstract:
The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physi…
▽ More
The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling.
△ Less
Submitted 9 September, 2024; v1 submitted 24 April, 2024;
originally announced April 2024.
-
A new approach for deducing rms proton radii from charge-changing reactions of neutron-rich nuclei and the reaction-target dependence
Authors:
J. -C. Zhang,
B. -H. Sun,
I. Tanihata,
R. Kanungo,
C. Scheidenberger,
S. Terashima,
Feng Wang,
F. Ameil,
J. Atkinson,
Y. Ayyad,
S. Bagchi,
D. Cortina-Gil,
I. Dillmann,
A. Estrade,
A. Evdokimov,
F. Farinon,
H. Geissel,
G. Guastalla,
R. Janik,
S. Kaur,
R. Knobel,
J. Kurcewicz,
Yu. A. Litvinov,
M. Marta,
M. Mostazo
, et al. (13 additional authors not shown)
Abstract:
We report the charge-changing cross sections ($σ_{\text{cc}}$) of 24 $p$-shell nuclides on both hydrogen and carbon at about 900$A$ MeV, of which $^{8,9}$Li, $^{10\textendash12}$Be, $^{10,14,15}$B, $^{14,15,17\textendash22}$N and $^{16}$O on hydrogen and $^{8,9}$Li on carbon are for the first time. Benefiting from the data set,we found a new and robust relationship between the scaling factor of th…
▽ More
We report the charge-changing cross sections ($σ_{\text{cc}}$) of 24 $p$-shell nuclides on both hydrogen and carbon at about 900$A$ MeV, of which $^{8,9}$Li, $^{10\textendash12}$Be, $^{10,14,15}$B, $^{14,15,17\textendash22}$N and $^{16}$O on hydrogen and $^{8,9}$Li on carbon are for the first time. Benefiting from the data set,we found a new and robust relationship between the scaling factor of the Glauber model calculations and the separation energies of the nuclei of interest on both targets.This allows us to deduce proton radii ($R_p$) for the first time from the cross sections on hydrogen. Nearly identical $R_p$ values are deduced from both target data for the neutron-rich carbon isotopes, however, the $R_p$ from the hydrogen target is systematically smaller in the neutron-rich nitrogen isotopes.This calls for further experimental and theoretical investigations.
△ Less
Submitted 31 March, 2024;
originally announced April 2024.
-
Machine Learning LSST 3x2pt analyses -- forecasting the impact of systematics on cosmological constraints using neural networks
Authors:
Supranta S. Boruah,
Tim Eifler,
Vivian Miranda,
Elyas Farah,
Jay Motka,
Elisabeth Krause,
Xiao Fang,
Paul Rogozenski,
The LSST Dark Energy Science Collaboration
Abstract:
Validating modeling choices through simulated analyses and quantifying the impact of different systematic effects will form a major computational bottleneck in the preparation for 3$\times$2 analysis with Stage-IV surveys such as Vera Rubin Observatory's Legacy Survey of Space and Time (LSST). We can significantly reduce the computational requirements by using machine learning based emulators, whi…
▽ More
Validating modeling choices through simulated analyses and quantifying the impact of different systematic effects will form a major computational bottleneck in the preparation for 3$\times$2 analysis with Stage-IV surveys such as Vera Rubin Observatory's Legacy Survey of Space and Time (LSST). We can significantly reduce the computational requirements by using machine learning based emulators, which allow us to run fast inference while maintaining the full realism of the data analysis pipeline. In this paper, we use such an emulator to run simulated 3$\times$2 (cosmic shear, galaxy-galaxy lensing, and galaxy clustering) analyses for mock LSST-Y1/Y3/Y6/Y10 surveys and study the impact of various systematic effects (galaxy bias, intrinsic alignment, baryonic physics, shear calibration and photo-$z$ uncertainties). Closely following the DESC Science Requirement Document (with several updates) our main findings are: {\it a)} The largest contribution to the `systematic error budget' of LSST 3$\times$2 analysis comes from galaxy bias uncertainties, while the contribution of baryonic and shear calibration uncertainties are significantly less important. {\it b)} Tighter constraints on intrinsic alignment and photo-$z$ parameters can improve cosmological constraints noticeably, which illustrates synergies of LSST and spectroscopic surveys. {\it c)} The scale cuts adopted in the DESC SRD may be too conservative and pushing to smaller scales can increase cosmological information significantly. {\it d)} We investigate the impact of photo-$z$ outliers on 3$\times$2 pt analysis and find that we need to determine the outlier fraction to within $5-10\%$ accuracy to ensure robust cosmological analysis. We caution that these findings depend on analysis choices (parameterizations, priors, scale cuts) and can change for different settings.
△ Less
Submitted 18 March, 2024;
originally announced March 2024.
-
Exploring Eco-Friendly Gas Mixtures for Resistive Plate Chambers: A Comprehensive Study on Performance and Aging
Authors:
The RPC ECOGas@GIF++ collaboration,
:,
L. Quaglia,
M. Abbrescia,
G. Aielli,
R. Aly,
M. C. Arena,
M. Barroso,
L. Benussi,
S. Bianco,
D. Boscherini,
F. Bordon,
A. Bruni,
S. Buontempo,
M. Busato,
P. Camarri,
R. Cardarelli,
L. Congedo,
D. De Jesus Damiao,
M. De Serio,
A. Di Ciaccio,
L. Di Stante,
P. Dupieux,
J. Eysermans,
A. Ferretti
, et al. (35 additional authors not shown)
Abstract:
Resistive Plate Chambers (RPCs) are gaseous detectors widely used in high energy physics experiments, operating with a gas mixture primarily containing Tetrafluoroethane (C$_{2}$H$_{2}$F$_{4}$), commonly known as R-134a, which has a global warming potential (GWP) of 1430. To comply with European regulations and explore environmentally friendly alternatives, the RPC EcoGas@GIF++ collaboration, invo…
▽ More
Resistive Plate Chambers (RPCs) are gaseous detectors widely used in high energy physics experiments, operating with a gas mixture primarily containing Tetrafluoroethane (C$_{2}$H$_{2}$F$_{4}$), commonly known as R-134a, which has a global warming potential (GWP) of 1430. To comply with European regulations and explore environmentally friendly alternatives, the RPC EcoGas@GIF++ collaboration, involving ALICE, ATLAS, CMS, LHCb/SHiP, and EP-DT communities, has undertaken intensive R\&D efforts to explore new gas mixtures for RPC technology.
A leading alternative under investigation is HFO1234ze, boasting a low GWP of 6 and demonstrating reasonable performance compared to R-134a. Over the past few years, RPC detectors with slightly different characteristics and electronics have been studied using HFO and CO$_{2}$-based gas mixtures at the CERN Gamma Irradiation Facility. An aging test campaign was launched in August 2022, and during the latest test beam in July 2023, all detector systems underwent evaluation. This contribution will report the results of the aging studies and the performance evaluations of the detectors with and without irradiation.
△ Less
Submitted 29 February, 2024;
originally announced February 2024.
-
Recovered SN Ia rate from simulated LSST images
Authors:
Vincenzo Petrecca,
Maria Teresa Botticella,
Enrico Cappellaro,
Laura Greggio,
Bruno Sánchez,
Anais Möller,
Masao Sako,
Melissa Graham,
Maurizio Paolillo,
Federica Bianco,
the LSST Dark Energy Science Collaboration
Abstract:
The Legacy Survey of Space and Time (LSST) will revolutionize Time Domain Astronomy by detecting millions of transients. In particular, it is expected to increment the number of type Ia supernovae (SNIa) of a factor of 100 compared to existing samples up to z~1.2. Such a high number of events will dramatically reduce statistical uncertainties in the analysis of SNIa properties and rates. However,…
▽ More
The Legacy Survey of Space and Time (LSST) will revolutionize Time Domain Astronomy by detecting millions of transients. In particular, it is expected to increment the number of type Ia supernovae (SNIa) of a factor of 100 compared to existing samples up to z~1.2. Such a high number of events will dramatically reduce statistical uncertainties in the analysis of SNIa properties and rates. However, the impact of all other sources of uncertainty on the measurement must still be evaluated. The comprehension and reduction of such uncertainties will be fundamental both for cosmology and stellar evolution studies, as measuring the SNIa rate can put constraints on the evolutionary scenarios of different SNIa progenitors. We use simulated data from the DESC Data Challenge 2 (DC2) and LSST Data Preview 0 (DP0) to measure the SNIa rate on a 15 deg2 region of the Wide-Fast-Deep area. We select a sample of SN candidates detected on difference images, associate them to the host galaxy, and retrieve their photometric redshifts (z-phot). Then, we test different light curves classification methods, with and without redshift priors. We discuss how the distribution in redshift measured for the SN candidates changes according to the selected host galaxy and redshift estimate. We measure the SNIa rate analyzing the impact of uncertainties due to z-phot, host galaxy association and classification on the distribution in redshift of the starting sample. We found a 17% average lost fraction of SNIa with respect to the simulated sample. As 10% of the bias is due to the uncertainty on the z-phot alone (which also affects classification when used as a prior), it results to be the major source of uncertainty. We discuss possible reduction of the errors in the measurement of the SNIa rate, including synergies with other surveys, which may help using the rate to discriminate different progenitor models.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Improving Photometric Redshift Estimates with Training Sample Augmentation
Authors:
Irene Moskowitz,
Eric Gawiser,
John Franklin Crenshaw,
Brett H. Andrews,
Alex I. Malz,
Samuel Schmidt,
The LSST Dark Energy Science Collaboration
Abstract:
Large imaging surveys will rely on photometric redshifts (photo-z's), which are typically estimated through machine learning methods. Currently planned spectroscopic surveys will not be deep enough to produce a representative training sample for LSST, so we seek methods to improve the photo-z estimates that arise from non-representative training samples. Spectroscopic training samples for photo-z'…
▽ More
Large imaging surveys will rely on photometric redshifts (photo-z's), which are typically estimated through machine learning methods. Currently planned spectroscopic surveys will not be deep enough to produce a representative training sample for LSST, so we seek methods to improve the photo-z estimates that arise from non-representative training samples. Spectroscopic training samples for photo-z's are biased towards redder, brighter galaxies, which also tend to be at lower redshift than the typical galaxy observed by LSST, leading to poor photo-z estimates with outlier fractions nearly 4 times larger than for a representative training sample. In this paper, we apply the concept of training sample augmentation, where we augment simulated non-representative training samples with simulated galaxies possessing otherwise unrepresented features. When we select simulated galaxies with (g-z) color, i-band magnitude and redshift outside the range of the original training sample, we are able to reduce the outlier fraction of the photo-z estimates for simulated LSST data by nearly 50% and the normalized median absolute deviation (NMAD) by 56%. When compared to a fully representative training sample, augmentation can recover nearly 70% of the degradation in the outlier fraction and 80% of the degradation in NMAD. Training sample augmentation is a simple and effective way to improve training samples for photo-z's without requiring additional spectroscopic samples.
△ Less
Submitted 14 May, 2024; v1 submitted 23 February, 2024;
originally announced February 2024.
-
Curvature in the very-high energy gamma-ray spectrum of M87
Authors:
H. E. S. S. Collaboration,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
H. Ashkar,
M. Backes,
V. Barbosa Martins,
R. Batzofin,
Y. Becherini,
D. Berge,
K. Bernlöhr,
M. Böttcher,
C. Boisson,
J. Bolmont,
M. de Bony de Lavergne,
F. Bradascio,
R. Brose,
F. Brun,
B. Bruno,
T. Bulik C. Burger-Scheidlin,
T. Bylund,
S. Casanova,
R. Cecil,
J. Celic,
M. Cerruti
, et al. (110 additional authors not shown)
Abstract:
The radio galaxy M87 is a variable very-high energy (VHE) gamma-ray source, exhibiting three major flares reported in 2005, 2008, and 2010. Despite extensive studies, the origin of the VHE gamma-ray emission is yet to be understood. In this study, we investigate the VHE gamma-ray spectrum of M87 during states of high gamma-ray activity, utilizing 20.2$\,$ hours the H.E.S.S. observations. Our findi…
▽ More
The radio galaxy M87 is a variable very-high energy (VHE) gamma-ray source, exhibiting three major flares reported in 2005, 2008, and 2010. Despite extensive studies, the origin of the VHE gamma-ray emission is yet to be understood. In this study, we investigate the VHE gamma-ray spectrum of M87 during states of high gamma-ray activity, utilizing 20.2$\,$ hours the H.E.S.S. observations. Our findings indicate a preference for a curved spectrum, characterized by a log-parabola model with extra-galactic background light (EBL) model above 0.3$\,$TeV at the 4$σ$ level, compared to a power-law spectrum with EBL. We investigate the degeneracy between the absorption feature and the EBL normalization and derive upper limits on EBL models mainly sensitive in the wavelength range 12.4$\,$$μ$m - 40$\,$$μ$m.
△ Less
Submitted 25 April, 2024; v1 submitted 20 February, 2024;
originally announced February 2024.
-
TeV flaring activity of the AGN PKS 0625-354 in November 2018
Authors:
H. E. S. S. Collaboration,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
H. Ashkar,
M. Backes,
A. Baktash,
V. Barbosa Martins,
J. Barnard,
R. Batzofin,
Y. Becherini,
D. Berge,
K. Bernlöhr,
B. Bi,
M. Böttcher,
C. Boisson,
J. Bolmont,
M. de Bony de Lavergne,
J. Borowska,
F. Bradascio,
M. Breuhaus,
R. Brose,
A. Brown,
F. Brun,
B. Bruno
, et al. (117 additional authors not shown)
Abstract:
Most $γ$-ray detected active galactic nuclei are blazars with one of their relativistic jets pointing towards the Earth. Only a few objects belong to the class of radio galaxies or misaligned blazars. Here, we investigate the nature of the object PKS 0625-354, its $γ$-ray flux and spectral variability and its broad-band spectral emission with observations from H.E.S.S., Fermi-LAT, Swift-XRT, and U…
▽ More
Most $γ$-ray detected active galactic nuclei are blazars with one of their relativistic jets pointing towards the Earth. Only a few objects belong to the class of radio galaxies or misaligned blazars. Here, we investigate the nature of the object PKS 0625-354, its $γ$-ray flux and spectral variability and its broad-band spectral emission with observations from H.E.S.S., Fermi-LAT, Swift-XRT, and UVOT taken in November 2018. The H.E.S.S. light curve above 200 GeV shows an outburst in the first night of observations followed by a declining flux with a halving time scale of 5.9h. The $γγ$-opacity constrains the upper limit of the angle between the jet and the line of sight to $\sim10^\circ$. The broad-band spectral energy distribution shows two humps and can be well fitted with a single-zone synchrotron self Compton emission model. We conclude that PKS 0625-354, as an object showing clear features of both blazars and radio galaxies, can be classified as an intermediate active galactic nuclei. Multi-wavelength studies of such intermediate objects exhibiting features of both blazars and radio galaxies are sparse but crucial for the understanding of the broad-band emission of $γ$-ray detected active galactic nuclei in general.
△ Less
Submitted 13 January, 2024;
originally announced January 2024.
-
Detecting strongly-lensed type Ia supernovae with LSST
Authors:
Nikki Arendse,
Suhail Dhawan,
Ana Sagués Carracedo,
Hiranya V. Peiris,
Ariel Goobar,
Radek Wojtak,
Catarina Alves,
Rahul Biswas,
Simon Huber,
Simon Birrer,
The LSST Dark Energy Science Collaboration
Abstract:
Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (…
▽ More
Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (SNIa) in LSST by quantifying the expected annual number of detections, the impact of stellar microlensing, follow-up feasibility, and how to best separate lensed and unlensed SNIa. We simulate SNIa lensed by galaxies, using the current LSST baseline v3.0 cadence, and find an expected number of 44 lensed SNIa detections per year. Microlensing effects by stars in the lens galaxy are predicted to lower the lensed SNIa detections by $\sim 8 \%$. The lensed events can be separated from the unlensed ones by jointly considering their colours and peak magnitudes. We define a `gold sample' of $\sim 10$ lensed SNIa per year with time delay $> 10$ days, $> 5$ detections before light-curve peak, and sufficiently bright ($m_i < 22.5$ mag) for follow-up observations. In three years of LSST operations, such a sample is expected to yield a $1.5\%$ measurement of the Hubble constant.
△ Less
Submitted 16 May, 2024; v1 submitted 7 December, 2023;
originally announced December 2023.
-
Polarimetric Geometric Modeling for mm-VLBI Observations of Black Holes
Authors:
Freek Roelofs,
Michael D. Johnson,
Andrew Chael,
Michael Janssen,
Maciek Wielgus,
Avery E. Broderick,
The Event Horizon Telescope Collaboration
Abstract:
The Event Horizon Telescope (EHT) is a millimeter very long baseline interferometry (VLBI) array that has imaged the apparent shadows of the supermassive black holes M87* and Sagittarius A*. Polarimetric data from these observations contain a wealth of information on the black hole and accretion flow properties. In this work, we develop polarimetric geometric modeling methods for mm-VLBI data, foc…
▽ More
The Event Horizon Telescope (EHT) is a millimeter very long baseline interferometry (VLBI) array that has imaged the apparent shadows of the supermassive black holes M87* and Sagittarius A*. Polarimetric data from these observations contain a wealth of information on the black hole and accretion flow properties. In this work, we develop polarimetric geometric modeling methods for mm-VLBI data, focusing on approaches that fit data products with differing degrees of invariance to broad classes of calibration errors. We establish a fitting procedure using a polarimetric "m-ring" model to approximate the image structure near a black hole. By fitting this model to synthetic EHT data from general relativistic magnetohydrodynamic models, we show that the linear and circular polarization structure can be successfully approximated with relatively few model parameters. We then fit this model to EHT observations of M87* taken in 2017. In total intensity and linear polarization, the m-ring fits are consistent with previous results from imaging methods. In circular polarization, the m-ring fits indicate the presence of event-horizon-scale circular polarization structure, with a persistent dipolar asymmetry and orientation across several days. The same structure was recovered independently of observing band, used data products, and model assumptions. Despite this broad agreement, imaging methods do not produce similarly consistent results. Our circular polarization results, which imposed additional assumptions on the source structure, should thus be interpreted with some caution. Polarimetric geometric modeling provides a useful and powerful method to constrain the properties of horizon-scale polarized emission, particularly for sparse arrays like the EHT.
△ Less
Submitted 17 November, 2023;
originally announced November 2023.
-
The impact of environment on size: Galaxies are 50% smaller in the Fornax Cluster compared to the field
Authors:
Nushkia Chamba,
Matthew Hayes,
The LSST Dark Energy Science Collaboration
Abstract:
Size is a fundamental parameter for measuring the growth of galaxies and the role of the environment on their evolution. However, the conventional size definitions used for this purpose are often biased and miss the diffuse, outermost signatures of galaxy growth, including star formation and gas accretion. This issue is addressed by examining low surface brightness truncations or galaxy ``edges''…
▽ More
Size is a fundamental parameter for measuring the growth of galaxies and the role of the environment on their evolution. However, the conventional size definitions used for this purpose are often biased and miss the diffuse, outermost signatures of galaxy growth, including star formation and gas accretion. This issue is addressed by examining low surface brightness truncations or galaxy ``edges'' as a physically motivated tracer of size based on star formation thresholds. Our total sample consists of $\sim900$ galaxies with stellar masses ranging from $10^5 M_{\odot} < M_{\star} < 10^{11} M_{\odot}$. This sample of nearby cluster, group satellite and nearly isolated field galaxies was compiled using multi-band imaging from the Fornax Deep Survey, deep IAC Stripe 82 and Dark Energy Camera Legacy Surveys. We find that the edge radii scale as $R_{\rm edge} \propto M_{\star}^{0.42}$ with a very small intrinsic scatter ($\sim 0.07$ dex). The scatter is driven by the morphology and environment of galaxies. In both the cluster and field, early-type dwarfs are systematically smaller by $\sim20\%$ than the late-types. However, compared to the field galaxies in the Fornax cluster are the most impacted. At a fixed stellar mass, edges in the cluster can be found at $\sim$ 50\% smaller radii and the average stellar surface density at the edges is a factor of two higher $\sim 1\,M_{\odot}$/pc$^2$. Our findings support the rapid removal of loosely bound neutral hydrogen in hot, crowded environments which truncates galaxies outside-in earlier, preventing the formation of more extended sizes and lower density edges. Our results highlight the importance of deep imaging surveys to study the low surface brightness imprints of the large scale structure and environment on galaxy evolution.
△ Less
Submitted 30 May, 2024; v1 submitted 16 November, 2023;
originally announced November 2023.
-
An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations
Authors:
Nicholas Van Alfen,
Duncan Campbell,
Jonathan Blazek,
C. Danielle Leonard,
Francois Lanusse,
Andrew Hearin,
Rachel Mandelbaum,
The LSST Dark Energy Science Collaboration
Abstract:
We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce n…
▽ More
We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the IllustrisTNG suite. We create an HOD mock with TNG-like correlations using our method, constraining the associated IA model parameters, with the $χ^2_{\rm dof}$ between our model's correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation--position and orientation--orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to $0.1~h^{-1}{\rm Mpc}$. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.
△ Less
Submitted 11 September, 2024; v1 submitted 13 November, 2023;
originally announced November 2023.
-
Generation of realistic input parameters for simulating atmospheric point-spread functions at astronomical observatories
Authors:
Claire-Alice Hébert,
Joshua E. Meyers,
My H. Do,
Patricia R. Burchat,
the LSST Dark Energy Science Collaboration
Abstract:
High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear -- correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all…
▽ More
High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear -- correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all other sources of correlated image distortions must be modeled or removed. One such source is the correlated blurring of images due to optical turbulence in the atmosphere, which dominates the point-spread function (PSF) for ground-based instruments. In this work, we leverage data from weather forecasting models to produce wind speeds and directions, and turbulence parameters, that are realistically correlated with altitude. To study the resulting correlations in the size and shape of the PSF, we generate simulated images of the PSF across a ~10 square-degree field of view -- the size of the camera focal plane for the Vera C. Rubin Observatory in Chile -- using weather data and historical seeing for a geographic location near the Observatory. We make quantitative predictions for two-point correlation functions (2PCF) that are used in analyses of cosmic shear. We observe a strong anisotropy in the two-dimensional 2PCF, which is expected based on observations in real images, and study the dependence of the orientation of the anisotropy on dominant wind directions near the ground and at higher altitudes.
The code repository for producing the correlated weather parameters for input to simulations (psf-weather-station) is public at https://github.com/LSSTDESC/psf-weather-station.
△ Less
Submitted 21 March, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
Discovery of a Radiation Component from the Vela Pulsar Reaching 20 Teraelectronvolts
Authors:
The H. E. S. S. Collaboration,
:,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
H. Ashkar,
M. Backes,
V. Barbosa Martins,
R. Batzofin,
Y. Becherini,
D. Berge,
K. Bernlöhr,
B. Bi,
M. Böttcher,
C. Boisson,
J. Bolmont,
M. de Bony de Lavergne,
J. Borowska,
F. Bradascio,
M. Breuhaus,
R. Brose,
F. Brun,
B. Bruno,
T. Bulik,
C. Burger-Scheidlin
, et al. (157 additional authors not shown)
Abstract:
Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a fe…
▽ More
Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a few gigaelectronvolt (GeV). Using the H.E.S.S. array of Cherenkov telescopes, we discovered a novel radiation component emerging beyond this generic GeV cutoff in the Vela pulsar's broadband spectrum. The extension of gamma-ray pulsation energies up to at least 20 teraelectronvolts (TeV) shows that Vela pulsar can accelerate particles to Lorentz factors higher than $4\times10^7$. This is an order of magnitude larger than in the case of the Crab pulsar, the only other pulsar detected in the TeV energy range. Our results challenge the state-of-the-art models for high-energy emission of pulsars while providing a new probe, i.e. the energetic multi-TeV component, for constraining the acceleration and emission processes in their extreme energy limit.
△ Less
Submitted 9 October, 2023;
originally announced October 2023.
-
Cosmology with imaging galaxy surveys: the impact of evolving galaxy bias and magnification
Authors:
Shivam Pandey,
Carles Sánchez,
Bhuvnesh Jain,
The LSST Dark Energy Science Collaboration
Abstract:
The joint analysis of galaxy clustering and galaxy-shear cross-correlations (galaxy-galaxy lensing) in imaging surveys constitutes one of the main avenues to obtain cosmological information. Analyses from Stage III surveys have assumed constant galaxy bias and weak lensing magnification coefficients within redshift bins. We use simple models of bias and magnification evolution to show that plausib…
▽ More
The joint analysis of galaxy clustering and galaxy-shear cross-correlations (galaxy-galaxy lensing) in imaging surveys constitutes one of the main avenues to obtain cosmological information. Analyses from Stage III surveys have assumed constant galaxy bias and weak lensing magnification coefficients within redshift bins. We use simple models of bias and magnification evolution to show that plausible levels of evolution impact cosmological parameter inference significantly for Stage IV surveys and possibly for current datasets as well. The distance weighting factors in galaxy-galaxy lensing cause it to have a different effective mean redshift within a lens bin than galaxy clustering. If the evolution is not included carefully in the theory modeling, this produces an apparent de-correlation of clustering and lensing even on large scales and leads to biased cosmological constraints. The level of bias is sensitive to the redshift and the width of the lens redshift distributions: for magnitude-limited-like samples it is comparable to the statistical precision even for the full DES survey, and exceeds it for Rubin-LSST. We find that ignoring the evolution can lead to approximately $1σ$ biases in the cosmological constraints in future Rubin-LSST analyses. We describe how these effects can be mitigated from the modeling and data sides, in particular by using cross-correlations to constrain the evolution parameters. We also demonstrate that the combination of galaxy clustering and CMB lensing does not suffer from this problem, due to the different distance weighting.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
Increasing the rate capability for the cryogenic stopping cell of the FRS Ion Catcher
Authors:
J. W. Zhao,
D. Amanbayev,
T. Dickel,
I. Miskun,
W. R. Plass,
N. Tortorelli,
S. Ayet San Andres,
Soenke Beck,
J. Bergmann,
Z. Brencic,
P. Constantin,
H. Geissel,
F. Greiner,
L. Groef,
C. Hornung,
N. Kuzminzuk,
G. Kripko-Koncz,
I. Mardor,
I. Pohjalainen,
C. Scheidenberger,
P. G. Thirolf,
S. Bagchi,
E. Haettner,
E. Kazantseva,
D. Kostyleva
, et al. (23 additional authors not shown)
Abstract:
At the FRS Ion Catcher (FRS-IC), projectile and fission fragments are produced at relativistic energies, separated in-flight, energy-bunched, slowed down, and thermalized in the ultra-pure helium gas-filled cryogenic stopping cell (CSC). Thermalized nuclei are extracted from the CSC using a combination of DC and RF electric fields and gas flow. This CSC also serves as the prototype CSC for the Sup…
▽ More
At the FRS Ion Catcher (FRS-IC), projectile and fission fragments are produced at relativistic energies, separated in-flight, energy-bunched, slowed down, and thermalized in the ultra-pure helium gas-filled cryogenic stopping cell (CSC). Thermalized nuclei are extracted from the CSC using a combination of DC and RF electric fields and gas flow. This CSC also serves as the prototype CSC for the Super-FRS, where exotic nuclei will be produced at unprecedented rates making it possible to go towards the extremes of the nuclear chart. Therefore, it is essential to efficiently extract thermalized exotic nuclei from the CSC under high beam rate conditions, in order to use the rare exotic nuclei which come as cocktail beams. The extraction efficiency dependence on the intensity of the impinging beam into the CSC was studied with a primary beam of 238U and its fragments. Tests were done with two different versions of the DC electrode structure inside the cryogenic chamber, the standard 1 m long and a short 0.5 m long DC electrode. In contrast to the rate capability of 10^4 ions/s with the long DC electrode, results show no extraction efficiency loss up to the rate of 2x10^5 ions/s with the new short DC electrode. This order of magnitude increase of the rate capability paves the way for new experiments at the FRS-IC, including exotic nuclei studies with in-cell multi-nucleon transfer reactions. The results further validate the design concept of the CSC for the Super-FRS, which was developed to effectively manage beams of even higher intensities.
△ Less
Submitted 4 August, 2023;
originally announced August 2023.
-
Slitless spectrophotometry with forward modelling: principles and application to atmospheric transmission measurement
Authors:
Jérémy Neveu,
Vincent Brémaud,
Pierre Antilogus,
Florent Barret,
Sébastien Bongard,
Yannick Copin,
Sylvie Dagoret-Campagne,
Claire Juramy,
Laurent Le-Guillou,
Marc Moniez,
Eduardo Sepulveda,
The LSST Dark Energy Science Collaboration
Abstract:
In the next decade, many optical surveys will aim to tackle the question of dark energy nature, measuring its equation of state parameter at the permil level. This requires trusting the photometric calibration of the survey with a precision never reached so far, controlling many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or on av…
▽ More
In the next decade, many optical surveys will aim to tackle the question of dark energy nature, measuring its equation of state parameter at the permil level. This requires trusting the photometric calibration of the survey with a precision never reached so far, controlling many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or on average for each season or for the full survey, can help reach the permil precision for magnitudes. This work aims at proving the ability to use slitless spectroscopy for standard star spectrophotometry and its use to monitor on-site atmospheric transmission as needed, for example, by the Vera C. Rubin Observatory Legacy Survey of Space and Time supernova cosmology program. We fully deal with the case of a disperser in the filter wheel, which is the configuration chosen in the Rubin Auxiliary Telescope. The theoretical basis of slitless spectrophotometry is at the heart of our forward model approach to extract spectroscopic information from slitless data. We developed a publicly available software called Spectractor (https://github.com/LSSTDESC/Spectractor) that implements each ingredient of the model and finally performs a fit of a spectrogram model directly on image data to get the spectrum. We show on simulations that our model allows us to understand the structure of spectrophotometric exposures. We also demonstrate its use on real data, solving specific issues and illustrating how our procedure allows the improvement of the model describing the data. Finally, we discuss how this approach can be used to directly extract atmospheric transmission parameters from data and thus provide the base for on-site atmosphere monitoring. We show the efficiency of the procedure on simulations and test it on the limited data set available.
△ Less
Submitted 10 June, 2024; v1 submitted 10 July, 2023;
originally announced July 2023.
-
Fat cosmic ray tracks in charge-coupled devices
Authors:
Theodore A. Grosson,
Andrei Nomerotski,
The LSST Dark Energy Science Collaboration
Abstract:
Cosmic rays are particles from the upper atmosphere which often leave bright spots and trails in images from telescope CCDs. We investigate so-called ``fat" cosmic rays seen in images from Vera C. Rubin Observatory and the Subaru Telescope. These tracks are much wider and brighter than typical cosmic ray tracks, and therefore are more capable of obscuring data in science images. By understanding t…
▽ More
Cosmic rays are particles from the upper atmosphere which often leave bright spots and trails in images from telescope CCDs. We investigate so-called ``fat" cosmic rays seen in images from Vera C. Rubin Observatory and the Subaru Telescope. These tracks are much wider and brighter than typical cosmic ray tracks, and therefore are more capable of obscuring data in science images. By understanding the origins of these tracks, we can better ensure that they do not interfere with on-sky data. We compare the properties of these tracks to simulated and theoretical models in order to identify both the particles causing these tracks as well as the reason for their excess spread. We propose that the origin of these tracks is cosmic ray protons, which deposit much greater charge in the CCDs than typical cosmic rays due to their lower velocities. The generated charges then repel each other while drifting through the detector, resulting in a track which is much wider than typical tracks.
△ Less
Submitted 19 June, 2023;
originally announced June 2023.
-
Mean range bunching of exotic nuclei produced by in-flight fragmentation and fission -- Stopped-beam experiments with increased efficiency
Authors:
Timo Dickel,
Christine Hornung,
Daler Amanbayev,
Samuel Ayet San Andres,
Soenke Beck,
Julian Bergmann,
Hans Geissel,
Juergen Gerl,
Magdalena Gorska,
Lizzy Groef,
Emma Haettner,
Jan-Paul Hucka,
Daria A. Kostyleva,
Gabriella Kripko-Koncz,
Ali Mollaebrahimi,
Ivan Mukha,
Stephane Pietri,
Wolfgang R. Plaß,
Zsolt Podolyak,
Sivaji Purushothaman,
Moritz Pascal Reiter,
Heidi Roesch,
Christoph Scheidenberger,
Yoshiki K. Tanaka,
Helmut Weick
, et al. (2 additional authors not shown)
Abstract:
The novel technique of mean range bunching has been developed and applied at the projectile fragment separator FRS at GSI in four experiments of the FAIR phase-0 experimental program. Using a variable degrader system at the final focal plane of the FRS, the ranges of the different nuclides can be aligned, allowing to efficiently implant a large number of different nuclides simultaneously in a gas-…
▽ More
The novel technique of mean range bunching has been developed and applied at the projectile fragment separator FRS at GSI in four experiments of the FAIR phase-0 experimental program. Using a variable degrader system at the final focal plane of the FRS, the ranges of the different nuclides can be aligned, allowing to efficiently implant a large number of different nuclides simultaneously in a gas-filled stopping cell or an implantation detector. Stopping and studying a cocktail beam overcomes the present limitations of stopped-beam experiments. The conceptual idea of mean range bunching is described and illustrated using simulations. In a single setting of the FRS, 37 different nuclides were stopped in the cryogenic stopping cell and were measured in a single setting broadband mass measurement with the multiple-reflection time-of-flight mass spectrometer of the FRS Ion Catcher.
△ Less
Submitted 30 May, 2023;
originally announced June 2023.
-
Study of a pulsar wind nebula candidate around the intermediate-age pulsar PSR J1413-6205 with H.E.S.S
Authors:
P. Chambéry,
Y. Gallant,
A. Jardin-Blicq,
M. Lemoine-Goumard,
V. Marandon,
A. Sinha,
M. Tsirou,
H. E. S. S. Collaboration
Abstract:
Very-high-energy $γ$-ray emission provides constraints on the morphology and the physics mechanisms involved in the evolution of pulsar wind nebulae (PWNe). In the Galactic plane, around $312 ^{\circ}$ of Galactic longitude, a promising region two-degree wide containing five powerful pulsars may offer a new insight on the transition between TeV-emitting PWNe and pulsar halos. Their rotational ener…
▽ More
Very-high-energy $γ$-ray emission provides constraints on the morphology and the physics mechanisms involved in the evolution of pulsar wind nebulae (PWNe). In the Galactic plane, around $312 ^{\circ}$ of Galactic longitude, a promising region two-degree wide containing five powerful pulsars may offer a new insight on the transition between TeV-emitting PWNe and pulsar halos. Their rotational energies range from $10^{35}$ to $10^{37}$ erg s$^{-1}$ for ages between 13.6 and 62.8 kyr. Extended emission is detected with H.E.S.S. (High Energy Stereoscopic System) in their vicinity, notably around the pulsar PSR J1413-6205.
We processed 124 hours of H.E.S.S observations with an algorithm improving background fitting for the study of extended sources. We applied a three-dimensional likelihood analysis technique to model the different sources in the region using a configuration that optimizes the collection area at the highest energies.
This contribution focuses on the detection of a new extended source around PSR J1413-6205 over 5$σ$ with a hard spectrum. Preliminary results on this source show a radius of $0.12 ^{\circ}$ $\pm$ $0.01 ^{\circ}_{\rm stat}$, an index of 2.06 $\pm$ 0.20$_{\rm stat}$ and a lower limit on a cut-off energy of 17 TeV, at a 90% confidence level. The detected emission is consistent with previous PWN models.
△ Less
Submitted 4 June, 2023;
originally announced June 2023.
-
Are classification metrics good proxies for SN Ia cosmological constraining power?
Authors:
Alex I. Malz,
Mi Dai,
Kara A. Ponder,
Emille E. O. Ishida,
Santiago Gonzalez-Gaitain,
Rupesh Durgesh,
Alberto Krone-Martins,
Rafael S. de Souza,
Noble Kennamer,
Sreevarsha Sreejith,
Lluis Galbany,
The LSST Dark Energy Science Collaboration,
The Cosmostatistics Initiative
Abstract:
Context: When selecting a classifier to use for a supernova Ia (SN Ia) cosmological analysis, it is common to make decisions based on metrics of classification performance, i.e. contamination within the photometrically classified SN Ia sample, rather than a measure of cosmological constraining power. If the former is an appropriate proxy for the latter, this practice would save those designing an…
▽ More
Context: When selecting a classifier to use for a supernova Ia (SN Ia) cosmological analysis, it is common to make decisions based on metrics of classification performance, i.e. contamination within the photometrically classified SN Ia sample, rather than a measure of cosmological constraining power. If the former is an appropriate proxy for the latter, this practice would save those designing an analysis pipeline from the computational expense of a full cosmology forecast. Aims: This study tests the assumption that classification metrics are an appropriate proxy for cosmology metrics. Methods: We emulate photometric SN Ia cosmology samples with controlled contamination rates of individual contaminant classes and evaluate each of them under a set of classification metrics. We then derive cosmological parameter constraints from all samples under two common analysis approaches and quantify the impact of contamination by each contaminant class on the resulting cosmological parameter estimates. Results: We observe that cosmology metrics are sensitive to both the contamination rate and the class of the contaminating population, whereas the classification metrics are insensitive to the latter. Conclusions: We therefore discourage exclusive reliance on classification-based metrics for cosmological analysis design decisions, e.g. classifier choice, and instead recommend optimizing using a metric of cosmological parameter constraining power.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
Constraining the cosmic-ray pressure in the inner Virgo Cluster using H.E.S.S. observations of M 87
Authors:
H. E. S. S. Collaboration,
:,
F. Aharonian,
F. Ait Benkhali,
C. Arcaro,
J. Aschersleben,
M. Backes,
V. Barbosa Martins,
R. Batzofin,
Y. Becherini,
D. Berge,
K. Bernlöhr,
B. Bi,
M. Böttcher,
C. Boisson,
J. Bolmont,
J. Borowska,
F. Bradascio,
M. Breuhaus,
R. Brose,
F. Brun,
B. Bruno,
T. Bulik,
C. Burger-Scheidlin,
T. Bylund
, et al. (139 additional authors not shown)
Abstract:
The origin of the gamma-ray emission from M87 is currently a matter of debate. This work aims to localize the VHE (100 GeV-100 TeV) gamma-ray emission from M87 and probe a potential extended hadronic emission component in the inner Virgo Cluster. The search for a steady and extended gamma-ray signal around M87 can constrain the cosmic-ray energy density and the pressure exerted by the cosmic rays…
▽ More
The origin of the gamma-ray emission from M87 is currently a matter of debate. This work aims to localize the VHE (100 GeV-100 TeV) gamma-ray emission from M87 and probe a potential extended hadronic emission component in the inner Virgo Cluster. The search for a steady and extended gamma-ray signal around M87 can constrain the cosmic-ray energy density and the pressure exerted by the cosmic rays onto the intra-cluster medium, and allow us to investigate the role of the cosmic rays in the active galactic nucleus feedback as a heating mechanism in the Virgo Cluster. H.E.S.S. telescopes are sensitive to VHE gamma rays and have been utilized to observe M87 since 2004. We utilized a Bayesian block analysis to identify M87 emission states with H.E.S.S. observations from 2004 until 2021, dividing them into low, intermediate, and high states. Because of the causality argument, an extended ($\gtrsim$kpc) signal is allowed only in steady emission states. Hence, we fitted the morphology of the 120h low state data and found no significant gamma-ray extension. Therefore, we derived for the low state an upper limit of 58"(corresponding to $\approx$4.6kpc) in the extension of a single-component morphological model described by a rotationally symmetric 2D Gaussian model at 99.7% confidence level. Our results exclude the radio lobes ($\approx$30 kpc) as the principal component of the VHE gamma-ray emission from the low state of M87. The gamma-ray emission is compatible with a single emission region at the radio core of M87. These results, with the help of two multiple-component models, constrain the maximum cosmic-ray to thermal pressure ratio $X_{CR,max.}$$\lesssim$$0.32$ and the total energy in cosmic-ray protons (CRp) to $U_{CR}$$\lesssim$5$\times10^{58}$ erg in the inner 20kpc of the Virgo Cluster for an assumed CRp power-law distribution in momentum with spectral index $α_{p}$=2.1.
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
Forecasting the power of Higher Order Weak Lensing Statistics with automatically differentiable simulations
Authors:
Denise Lanzieri,
François Lanusse,
Chirag Modi,
Benjamin Horowitz,
Joachim Harnois-Déraps,
Jean-Luc Starck,
The LSST Dark Energy Science Collaboration
Abstract:
We present the Differentiable Lensing Lightcone (DLL), a fully differentiable physical model designed for being used as a forward model in Bayesian inference algorithms requiring access to derivatives of lensing observables with respect to cosmological parameters. We extend the public FlowPM N-body code, a particle-mesh N-body solver, simulating lensing lightcones and implementing the Born approxi…
▽ More
We present the Differentiable Lensing Lightcone (DLL), a fully differentiable physical model designed for being used as a forward model in Bayesian inference algorithms requiring access to derivatives of lensing observables with respect to cosmological parameters. We extend the public FlowPM N-body code, a particle-mesh N-body solver, simulating lensing lightcones and implementing the Born approximation in the Tensorflow framework. Furthermore, DLL is aimed at achieving high accuracy with low computational costs. As such, it integrates a novel Hybrid Physical-Neural parameterisation able to compensate for the small-scale approximations resulting from particle-mesh schemes for cosmological N-body simulations. We validate our simulations in an LSST setting against high-resolution $κ$TNG simulations by comparing both the lensing angular power spectrum and multiscale peak counts. We demonstrate an ability to recover lensing $C_\ell$ up to a 10% accuracy at $\ell=1000$ for sources at redshift 1, with as few as $\sim 0.6$ particles per Mpc/h. As a first use case, we use this tool to investigate the relative constraining power of the angular power spectrum and peak counts statistic in an LSST setting. Such comparisons are typically very costly as they require a large number of simulations, and do not scale well with the increasing number of cosmological parameters. As opposed to forecasts based on finite differences, these statistics can be analytically differentiated with respect to cosmology, or any systematics included in the simulations at the same computational cost of the forward simulation. We find that the peak counts outperform the power spectrum on the cold dark matter parameter $Ω_c$, on the amplitude of density fluctuations $σ_8$, and on the amplitude of the intrinsic alignment signal $A_{IA}$.
△ Less
Submitted 12 May, 2023;
originally announced May 2023.
-
Detection of extended gamma-ray emission around the Geminga pulsar with H.E.S.S
Authors:
H. E. S. S. Collaboration,
:,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
H. Ashkar,
M. Backes,
V. Barbosa Martins,
R. Batzofin,
Y. Becherini,
D. Berge,
K. Bernlöhr,
B. Bi,
M. Böttcher,
C. Boisson,
J. Bolmont,
J. Borowska,
M. Bouyahiaoui,
F. Bradascio,
R. Brose,
F. Brun,
B. Bruno,
T. Bulik,
C. Burger Scheidlin,
F. Cangemi
, et al. (143 additional authors not shown)
Abstract:
Geminga is an enigmatic radio-quiet gamma-ray pulsar located at a mere 250 pc distance from Earth. Extended very-high-energy gamma-ray emission around the pulsar was discovered by Milagro and later confirmed by HAWC, which are both water Cherenkov detector-based experiments. However, evidence for the Geminga pulsar wind nebula in gamma rays has long evaded detection by imaging atmospheric Cherenko…
▽ More
Geminga is an enigmatic radio-quiet gamma-ray pulsar located at a mere 250 pc distance from Earth. Extended very-high-energy gamma-ray emission around the pulsar was discovered by Milagro and later confirmed by HAWC, which are both water Cherenkov detector-based experiments. However, evidence for the Geminga pulsar wind nebula in gamma rays has long evaded detection by imaging atmospheric Cherenkov telescopes (IACTs) despite targeted observations. The detection of gamma-ray emission on angular scales > 2 deg poses a considerable challenge for the background estimation in IACT data analysis. With recent developments in understanding the complementary background estimation techniques of water Cherenkov and atmospheric Cherenkov instruments, the H.E.S.S. IACT array can now confirm the detection of highly extended gamma-ray emission around the Geminga pulsar with a radius of at least 3 deg in the energy range 0.5-40 TeV. We find no indications for statistically significant asymmetries or energy-dependent morphology. A flux normalisation of $(2.8\pm0.7)\times10^{-12}$ cm$^{-2}$s$^{-1}$TeV$^{-1}$ at 1 TeV is obtained within a 1 deg radius region around the pulsar. To investigate the particle transport within the halo of energetic leptons around the pulsar, we fitted an electron diffusion model to the data. The normalisation of the diffusion coefficient obtained of $D_0 = 7.6^{+1.5}_{-1.2} \times 10^{27}$ cm$^2$s$^{-1}$, at an electron energy of 100 TeV, is compatible with values previously reported for the pulsar halo around Geminga, which is considerably below the Galactic average.
△ Less
Submitted 5 April, 2023;
originally announced April 2023.
-
Search for the evaporation of primordial black holes with H.E.S.S
Authors:
H. E. S. S. collaboration,
:,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
M. Boettcher,
M. Backes,
V. Barbosa Martins,
R. Batzo,
Y. Becherini,
D. Berge,
B. Bi,
C. Boisson,
J. Bolmont,
M. de Bony de Lavergne,
J. Borowska,
F. Bradascio,
R. Brose,
F. Brun,
B. Bruno,
T. Bulik,
C. Burger-Scheidlin,
S. Caro,
S. Casanova,
J. Celic
, et al. (124 additional authors not shown)
Abstract:
Primordial Black Holes (PBHs) are hypothetical black holes predicted to have been formed from density fluctuations in the early Universe. PBHs with an initial mass around $10^{14}-10^{15}$g are expected to end their evaporation at present times in a burst of particles and very-high-energy (VHE) gamma rays. Those gamma rays may be detectable by the High Energy Stereoscopic System (H.E.S.S.), an arr…
▽ More
Primordial Black Holes (PBHs) are hypothetical black holes predicted to have been formed from density fluctuations in the early Universe. PBHs with an initial mass around $10^{14}-10^{15}$g are expected to end their evaporation at present times in a burst of particles and very-high-energy (VHE) gamma rays. Those gamma rays may be detectable by the High Energy Stereoscopic System (H.E.S.S.), an array of imaging atmospheric Cherenkov telescopes. This paper reports on the search for evaporation bursts of VHE gamma rays with H.E.S.S., ranging from 10 to 120 seconds, as expected from the final stage of PBH evaporation and using a total of 4816 hours of observations. The most constraining upper limit on the burst rate of local PBHs is $2000$ pc$^{-3}$ yr$^{-1}$ for a burst interval of 120 seconds, at the 95\% confidence level. The implication of these measurements for PBH dark matter are also discussed.
△ Less
Submitted 22 March, 2023;
originally announced March 2023.
-
H.E.S.S. follow-up observations of GRB221009A
Authors:
H. E. S. S. Collaboration,
:,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
H. Ashkar,
M. Backes,
A. Baktash,
V. Barbosa Martins,
R. Batzofin,
Y. Becherini,
D. Berge,
K. Bernlöhr,
B. Bi,
M. Böttcher,
C. Boisson,
J. Bolmont,
M. de Bony de Lavergne,
J. Borowska,
M. Bouyahiaoui,
F. Bradascio,
M. Breuhaus,
R. Brose,
F. Brun,
B. Bruno
, et al. (138 additional authors not shown)
Abstract:
GRB221009A is the brightest gamma-ray burst ever detected. To probe the very-high-energy (VHE, $>$\!100 GeV) emission, the High Energy Stereoscopic System (H.E.S.S.) began observations 53 hours after the triggering event, when the brightness of the moonlight no longer precluded observations. We derive differential and integral upper limits using H.E.S.S. data from the third, fourth, and ninth nigh…
▽ More
GRB221009A is the brightest gamma-ray burst ever detected. To probe the very-high-energy (VHE, $>$\!100 GeV) emission, the High Energy Stereoscopic System (H.E.S.S.) began observations 53 hours after the triggering event, when the brightness of the moonlight no longer precluded observations. We derive differential and integral upper limits using H.E.S.S. data from the third, fourth, and ninth nights after the initial GRB detection, after applying atmospheric corrections. The combined observations yield an integral energy flux upper limit of $Φ_\mathrm{UL}^{95\%} = 9.7 \times 10^{-12}~\mathrm{erg\,cm^{-2}\,s^{-1}}$ above $E_\mathrm{thr} = 650$ GeV. The constraints derived from the H.E.S.S. observations complement the available multiwavelength data. The radio to X-ray data are consistent with synchrotron emission from a single electron population, with the peak in the SED occurring above the X-ray band. Compared to the VHE-bright GRB190829A, the upper limits for GRB221009A imply a smaller gamma-ray to X-ray flux ratio in the afterglow. Even in the absence of a detection, the H.E.S.S. upper limits thus contribute to the multiwavelength picture of GRB221009A, effectively ruling out an IC dominated scenario.
△ Less
Submitted 18 March, 2023;
originally announced March 2023.
-
Metadetection Weak Lensing for the Vera C. Rubin Observatory
Authors:
Erin S. Sheldon,
Matthew R. Becker,
Michael Jarvis,
Robert Armstrong,
The LSST Dark Energy Science Collaboration
Abstract:
Forthcoming astronomical imaging surveys will use weak gravitational lensing shear as a primary probe to study dark energy, with accuracy requirements at the 0.1% level. We present an implementation of the Metadetection shear measurement algorithm for use with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This new code works with the data products produced by the LSST Scien…
▽ More
Forthcoming astronomical imaging surveys will use weak gravitational lensing shear as a primary probe to study dark energy, with accuracy requirements at the 0.1% level. We present an implementation of the Metadetection shear measurement algorithm for use with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This new code works with the data products produced by the LSST Science Pipelines, and uses the pipeline algorithms when possible. We tested the code using a new set of simulations designed to mimic LSST imaging data. The simulated images contained semi-realistic galaxies, stars with representative distributions of magnitudes and galactic spatial density, cosmic rays, bad CCD columns and spatially variable point spread functions. Bright stars were saturated and simulated ``bleed trails'' were drawn. Problem areas were interpolated, and the images were coadded into small cells, excluding images not fully covering the cell to guarantee a continuous point spread function. In all our tests the measured shear was accurate within the LSST requirements.
△ Less
Submitted 11 May, 2023; v1 submitted 6 March, 2023;
originally announced March 2023.
-
Search for Two-neutrino Double-Beta Decay of $^{136}\rm Xe$ to the $0^+_1$ excited state of $^{136}\rm Ba$ with the Complete EXO-200 Dataset
Authors:
EXO-200 Collaboration,
:,
S. Al Kharusi,
G. Anton,
I. Badhrees,
P. S. Barbeau,
D. Beck,
V. Belov,
T. Bhatta,
M. Breidenbach,
T. Brunner,
G. F. Cao,
W. R. Cen,
C. Chambers,
B. Cleveland,
M. Coon,
A. Craycraft,
T. Daniels,
L. Darroch,
S. J. Daugherty,
J. Davis,
S. Delaquis,
A. Der Mesrobian-Kabakian,
R. DeVoe,
J. Dilling
, et al. (83 additional authors not shown)
Abstract:
A new search for two-neutrino double-beta ($2νββ$) decay of $^{136}\rm Xe$ to the $0^+_1$ excited state of $^{136}\rm Ba$ is performed with the full EXO-200 dataset. A deep learning-based convolutional neural network is used to discriminate signal from background events. Signal detection efficiency is increased relative to previous searches by EXO-200 by more than a factor of two. With the additio…
▽ More
A new search for two-neutrino double-beta ($2νββ$) decay of $^{136}\rm Xe$ to the $0^+_1$ excited state of $^{136}\rm Ba$ is performed with the full EXO-200 dataset. A deep learning-based convolutional neural network is used to discriminate signal from background events. Signal detection efficiency is increased relative to previous searches by EXO-200 by more than a factor of two. With the addition of the Phase II dataset taken with an upgraded detector, the median 90$\%$ confidence level half-life sensitivity of $2νββ$ decay to the $0^+_1$ state of $^{136}\rm Ba$ is $2.9 \times 10^{24}~\rm yr$ using a total $^{136}\rm Xe$ exposure of $234.1~\rm kg~yr$. No statistically significant evidence for $2νββ$ decay to the $0^+_1$ state is observed, leading to a lower limit of $T^{2ν}_{1/2}(0^+ \rightarrow 0^+_1) > 1.4\times10^{24}~\rm yr$ at 90$\%$ confidence level, improved by 70$\%$ relative to the current world's best constraint.
△ Less
Submitted 16 October, 2023; v1 submitted 2 March, 2023;
originally announced March 2023.
-
HESS J1809$-$193: a halo of escaped electrons around a pulsar wind nebula?
Authors:
H. E. S. S. Collaboration,
:,
F. Aharonian,
F. Ait Benkhali,
J. Aschersleben,
H. Ashkar,
M. Backes,
V. Barbosa Martins,
R. Batzofin,
Y. Becherini,
D. Berge,
M. Böttcher,
C. Boisson,
J. Bolmont,
J. Borowska,
M. Bouyahiaoui,
F. Bradascio,
M. Breuhaus,
R. Brose,
F. Brun,
B. Bruno,
T. Bulik,
C. Burger-Scheidlin,
T. Bylund,
S. Caroff
, et al. (130 additional authors not shown)
Abstract:
Context. HESS J1809$-$193 is an unassociated very-high-energy $γ$-ray source located on the Galactic plane. While it has been connected to the nebula of the energetic pulsar PSR J1809$-$1917, supernova remnants and molecular clouds present in the vicinity also constitute possible associations. Recently, the detection of $γ$-ray emission up to energies of $\sim$100 TeV with the HAWC observatory has…
▽ More
Context. HESS J1809$-$193 is an unassociated very-high-energy $γ$-ray source located on the Galactic plane. While it has been connected to the nebula of the energetic pulsar PSR J1809$-$1917, supernova remnants and molecular clouds present in the vicinity also constitute possible associations. Recently, the detection of $γ$-ray emission up to energies of $\sim$100 TeV with the HAWC observatory has led to renewed interest in HESS J1809$-$193.
Aims. We aim to understand the origin of the $γ$-ray emission of HESS J1809$-$193.
Methods. We analysed 93.2 h of data taken on HESS J1809$-$193 above 0.27 TeV with the High Energy Stereoscopic System (H.E.S.S.), using a multi-component, three-dimensional likelihood analysis. In addition, we provide a new analysis of 12.5 yr of Fermi-LAT data above 1 GeV within the region of HESS J1809$-$193. The obtained results are interpreted in a time-dependent modelling framework.
Results. For the first time, we were able to resolve the emission detected with H.E.S.S. into two components: an extended component that exhibits a spectral cut-off at $\sim$13 TeV, and a compact component that is located close to PSR J1809$-$1917 and shows no clear spectral cut-off. The Fermi-LAT analysis also revealed extended $γ$-ray emission, on scales similar to that of the extended H.E.S.S. component.
Conclusions. Our modelling indicates that based on its spectrum and spatial extent, the extended H.E.S.S. component is likely caused by inverse Compton emission from old electrons that form a halo around the pulsar wind nebula. The compact component could be connected to either the pulsar wind nebula or the supernova remnant and molecular clouds. Due to its comparatively steep spectrum, modelling the Fermi-LAT emission together with the H.E.S.S. components is not straightforward. (abridged)
△ Less
Submitted 27 February, 2023;
originally announced February 2023.
-
Measurement of telescope transmission using a Collimated Beam Projector
Authors:
Nicholas Mondrik,
Michael Coughlin,
Marc Betoule,
Sébastien Bongard,
Joseph P. Rice,
Ping-Shine Shaw,
Christopher W. Stubbs,
John T. Woodward,
LSST Dark Energy Science Collaboration
Abstract:
With the increasingly large number of type Ia supernova being detected by current-generation survey telescopes, and even more expected with the upcoming Rubin Observatory Legacy Survey of Space and Time, the precision of cosmological measurements will become limited by systematic uncertainties in flux calibration rather than statistical noise. One major source of systematic error in determining SN…
▽ More
With the increasingly large number of type Ia supernova being detected by current-generation survey telescopes, and even more expected with the upcoming Rubin Observatory Legacy Survey of Space and Time, the precision of cosmological measurements will become limited by systematic uncertainties in flux calibration rather than statistical noise. One major source of systematic error in determining SNe Ia color evolution (needed for distance estimation) is uncertainty in telescope transmission, both within and between surveys. We introduce here the Collimated Beam Projector (CBP), which is meant to measure a telescope transmission with collimated light. The collimated beam more closely mimics a stellar wavefront as compared to flat-field based instruments, allowing for more precise handling of systematic errors such as those from ghosting and filter angle-of-incidence dependence. As a proof of concept, we present CBP measurements of the StarDICE prototype telescope, achieving a standard (1 sigma) uncertainty of 3 % on average over the full wavelength range measured with a single beam illumination.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
The Event Horizon Telescope Image of the Quasar NRAO 530
Authors:
Svetlana Jorstad,
Maciek Wielgus,
Rocco Lico,
Sara Issaoun,
Avery E. Broderick,
Dominic W. Pesce,
Jun Liu,
Guang-Yao Zhao,
Thomas P. Krichbaum,
Lindy Blackburn,
Chi-Kwan Chan,
Michael Janssen,
Venkatessh Ramakrishnan,
Kazunori Akiyama,
Antxon Alberdi,
Juan Carlos Algaba,
Katherine L. Bouman,
Ilje Cho,
Antonio Fuentes,
Jose L. Gomez,
Mark Gurwell,
Michael D. Johnson,
Jae-Young Kim,
Ru-Sen Lu,
Ivan Marti-Vidal
, et al. (5 additional authors not shown)
Abstract:
We report on the observations of the quasar NRAO 530 with the Event Horizon Telescope (EHT) on 2017 April 5-7, when NRAO 530 was used as a calibrator for the EHT observations of Sagittarius A*. At z=0.902 this is the most distant object imaged by the EHT so far. We reconstruct the first images of the source at 230 GHz, at an unprecedented angular resolution of $\sim$ 20 $μ$as, both in total intens…
▽ More
We report on the observations of the quasar NRAO 530 with the Event Horizon Telescope (EHT) on 2017 April 5-7, when NRAO 530 was used as a calibrator for the EHT observations of Sagittarius A*. At z=0.902 this is the most distant object imaged by the EHT so far. We reconstruct the first images of the source at 230 GHz, at an unprecedented angular resolution of $\sim$ 20 $μ$as, both in total intensity and in linear polarization. We do not detect source variability, allowing us to represent the whole data set with static images. The images reveal a bright feature located on the southern end of the jet, which we associate with the core. The feature is linearly polarized, with a fractional polarization of $\sim$5-8% and has a sub-structure consisting of two components. Their observed brightness temperature suggests that the energy density of the jet is dominated by the magnetic field. The jet extends over 60 $μ$as along a position angle PA$\sim -$28$^\circ$. It includes two features with orthogonal directions of polarization (electric vector position angle, EVPA), parallel and perpendicular to the jet axis, consistent with a helical structure of the magnetic field in the jet. The outermost feature has a particularly high degree of linear polarization, suggestive of a nearly uniform magnetic field. Future EHT observations will probe the variability of the jet structure on $μ$as scales, while simultaneous multi-wavelength monitoring will provide insight into the high energy emission origin.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST
Authors:
J. Prat,
J. Zuntz,
Y. Omori,
C. Chang,
T. Tröster,
E. Pedersen,
C. García-García,
E. Phillips-Longley,
J. Sanchez,
D. Alonso,
X. Fang,
E. Gawiser,
K. Heitmann,
M. Ishak,
M. Jarvis,
E. Kovacs,
P. Larsen,
Y. -Y. Mao,
L. Medina Varela,
M. Paterno,
S. D. Vitenti,
Z. Zhang,
The LSST Dark Energy Science Collaboration
Abstract:
We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Sci…
▽ More
We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis.
△ Less
Submitted 21 April, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Improved Tomographic Binning of 3x2pt Lens Samples: Neural Network Classifiers and Optimal Bin Assignments
Authors:
Irene Moskowitz,
Eric Gawiser,
Abby Bault,
Adam Broussard,
Jeffrey A. Newman,
Joe Zuntz,
the LSST Dark Energy Science Collaboration
Abstract:
Large imaging surveys, such as the Legacy Survey of Space and Time, rely on photometric redshifts and tomographic binning for 3x2pt analyses that combine galaxy clustering and weak lensing. In this paper, we propose a method for optimizing the tomographic binning choice for the lens sample of galaxies. We divide the CosmoDC2 and Buzzard simulated galaxy catalogs into a training set and an applicat…
▽ More
Large imaging surveys, such as the Legacy Survey of Space and Time, rely on photometric redshifts and tomographic binning for 3x2pt analyses that combine galaxy clustering and weak lensing. In this paper, we propose a method for optimizing the tomographic binning choice for the lens sample of galaxies. We divide the CosmoDC2 and Buzzard simulated galaxy catalogs into a training set and an application set, where the training set is nonrepresentative in a realistic way, and then estimate photometric redshifts for the application sets. The galaxies are sorted into redshift bins covering equal intervals of redshift or comoving distance, or with an equal number of galaxies in each bin, and we consider a generalized extension of these approaches. We find that bins of equal comoving distance produce the highest dark energy figure of merit of the initial binning choices, but that the choice of bin edges can be further optimized. We then train a neural network classifier to identify galaxies that are either highly likely to have accurate photometric redshift estimates or highly likely to be sorted into the correct redshift bin. The neural network classifier is used to remove poor redshift estimates from the sample, and the results are compared to the case when none of the sample is removed. We find that the neural network classifiers are able to improve the figure of merit by ~13% and are able to recover ~25% of the loss in the figure of merit that occurs when a nonrepresentative training sample is used.
△ Less
Submitted 14 June, 2023; v1 submitted 13 December, 2022;
originally announced December 2022.
-
The N5K Challenge: Non-Limber Integration for LSST Cosmology
Authors:
C. D. Leonard,
T. Ferreira,
X. Fang,
R. Reischke,
N. Schoeneberg,
T. Tröster,
D. Alonso,
J. E. Campagne,
F. Lanusse,
A. Slosar,
M. Ishak,
the LSST Dark Energy Science Collaboration
Abstract:
The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angula…
▽ More
The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angular power spectrum of galaxy number counts and cosmic shear data without invoking the so-called 'Limber approximation', in the context of the Rubin Observatory Legacy Survey of Space and Time (LSST). We quantify the performance, in terms of accuracy and speed, of three non-Limber implementations: ${\tt FKEM (CosmoLike)}$, ${\tt Levin}$, and ${\tt matter}$, themselves based on different integration schemes and approximations. We find that in the challenge's fiducial 3x2pt LSST Year 10 scenario, ${\tt FKEM (CosmoLike)}$ produces the fastest run time within the required accuracy by a considerable margin, positioning it favourably for use in Bayesian parameter inference. This method, however, requires further development and testing to extend its use to certain analysis scenarios, particularly those involving a scale-dependent growth rate. For this and other reasons discussed herein, alternative approaches such as ${\tt matter}$ and ${\tt Levin}$ may be necessary for a full exploration of parameter space. We also find that the usual first-order Limber approximation is insufficiently accurate for LSST Year 10 3x2pt analysis on $\ell=200-1000$, whereas invoking the second-order Limber approximation on these scales (with a full non-Limber method at smaller $\ell$) does suffice.
△ Less
Submitted 14 February, 2023; v1 submitted 8 December, 2022;
originally announced December 2022.