Showing 1–50 of 458 results for author: Collaboration, S

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.

Comments: accepted to A&A

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.

Comments: 19 pages, 16 figures, 3 tables; v2: added LSST DESC to arxiv author list

Initial measurement of reactor antineutrino oscillation at SNO+

Authors: SNO+ Collaboration, :, A. Allega, M. R. Anderson, S. Andringa, M. Askins, D. J. Auty, A. Bacon, J. Baker, F. Barão, N. Barros, R. Bayes, E. W. Beier, T. S. Bezerra, A. Bialek, S. D. Biller, E. Blucher, E. Caden, E. J. Callaghan, M. Chen, S. Cheng, B. Cleveland, D. Cookman, J. Corning, M. A. Cox , et al. (96 additional authors not shown)

Abstract: The SNO+ collaboration reports its first spectral analysis of long-baseline reactor antineutrino oscillation using 114 tonne-years of data. Fitting the neutrino oscillation probability to the observed energy spectrum yields constraints on the neutrino mass-squared difference $Δm^2_{21}$. In the ranges allowed by previous measurements, the best-fit $Δm^2_{21}$ is (8.85$^{+1.10}_{-1.33}$) $\times$ 1… ▽ More The SNO+ collaboration reports its first spectral analysis of long-baseline reactor antineutrino oscillation using 114 tonne-years of data. Fitting the neutrino oscillation probability to the observed energy spectrum yields constraints on the neutrino mass-squared difference $Δm^2_{21}$. In the ranges allowed by previous measurements, the best-fit $Δm^2_{21}$ is (8.85$^{+1.10}_{-1.33}$) $\times$ 10$^{-5}$ eV$^2$. This measurement is continuing in the next phases of SNO+ and is expected to surpass the present global precision on $Δm^2_{21}$ with about three years of data. △ Less

Submitted 30 May, 2024; originally announced May 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.

Comments: 19 pages, 13 figures, submitted to AJ

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.

Comments: 28 pages, 24 figures, submitted to The Astrophysical Journal

Broadband Multi-wavelength Properties of M87 during the 2018 EHT Campaign including a Very High Energy Flaring Episode

Authors: The Event Horizon Telescope - Multi-wavelength science working group, The Event Horizon Telescope Collaboration, The Fermi Large Area Telescope Collaboration, H. E. S. S. Collaboration, MAGIC Collaboration, VERITAS Collaboration, EAVN Collaboration

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 24 April, 2024; originally announced April 2024.

Comments: 46 pages, 23 figures, submitted to Astronomy & Astrophysics

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.

Comments: 5 pages, 3 figures

Journal ref: Science Bulletin 69,1647(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.

Comments: 16 pages, 10 Figures, To be submitted to PRD. Comments welcome

The classical T Tauri star CI Tau observed with SPIRou: magnetospheric accretion and planetary formation

Authors: J. -F. Donati, B. Finociety, P. I. Cristofari, S. H. P. Alencar, C. Moutou, X. Delfosse, P Fouqué, L. Arnold, C. Baruteau, Á. Kóspál, F. Ménard, A. Carmona, K. Grankin, M. Takami, E. Artigau, R. Doyon, G. Hébrard, the SLS collaboration

Abstract: We report new observations of the classical T~Tauri star CI~Tau with the SPIRou near-infrared spectropolarimeter and velocimeter at the Canada-France-Hawaii Telescope (CFHT) in late 2019, 2020 and 2022, complemented with observations obtained with the ESPaDOnS optical spectropolarimeter at CFHT in late 2020. From our SPIRou and ESPaDOnS spectra, to which we applied Least-Squares Deconvolution, we… ▽ More We report new observations of the classical T~Tauri star CI~Tau with the SPIRou near-infrared spectropolarimeter and velocimeter at the Canada-France-Hawaii Telescope (CFHT) in late 2019, 2020 and 2022, complemented with observations obtained with the ESPaDOnS optical spectropolarimeter at CFHT in late 2020. From our SPIRou and ESPaDOnS spectra, to which we applied Least-Squares Deconvolution, we infer longitudinal fields clearly modulated with the 9-d rotation period of CI~Tau. Using Zeeman-Doppler imaging, we reconstruct the large-scale magnetic topology, first from SPIRou data only in all three seasons, then from our 2020 SPIRou and ESPaDOnS data simultaneously. We find that CI~Tau hosts a mainly axisymmetric poloidal field, with a 1~kG dipole slightly tilted to the rotation axis and dark spots close to the pole that coincide with the footpoints of accretion funnels linking the star to the inner disc. Our results also suggest that CI~Tau accretes mass from the disc in a stable fashion. We further find that radial velocities (RV) derived from atomic and CO lines in SPIRou spectra are both rotationally modulated, but with a much lower amplitude than that expected from the putative candidate planet CI~Tau~b. We confirm the presence of a RV signal at a period of 23.86~d reported in a separate analysis, but detect it clearly in CO lines only and not in atomic lines, suggesting that it likely traces a non-axisymmetric structure in the inner disc of CI~Tau rather than a massive close-in planet. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: MNRAS, in press (22 pages, 11 figures, 6 tables)

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.

Comments: Accepted for publication in A&A

Journal ref: A&A 686, A11 (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.

Comments: 11 pages, 4 figures, published in ApJ Letters

Journal ref: ApJL 967 L6 (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.

Comments: 10 pages, 7 figures. Accepted for publication in A&A. Corresponding authors: Victor Barbosa Martins, Rahul Cecil, Iryna Lypova, Manuel Meyer, Perri Zilberman. Supplementary material: https://zenodo.org/records/10781524

Journal ref: A&A, 685, A96 (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.

Comments: 9 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

Solar neutrino measurements using the full data period of Super-Kamiokande-IV

Authors: Super-Kamiokande Collaboration, :, K. Abe, C. Bronner, Y. Hayato, K. Hiraide, K. Hosokawa, K. Ieki, M. Ikeda, S. Imaizumi, K. Iyogi, J. Kameda, Y. Kanemura, R. Kaneshima, Y. Kashiwagi, Y. Kataoka, Y. Kato, Y. Kishimoto, S. Miki, S. Mine, M. Miura, T. Mochizuki, S. Moriyama, Y. Nagao, M. Nakahata , et al. (305 additional authors not shown)

Abstract: An analysis of solar neutrino data from the fourth phase of Super-Kamiokande~(SK-IV) from October 2008 to May 2018 is performed and the results are presented. The observation time of the data set of SK-IV corresponds to $2970$~days and the total live time for all four phases is $5805$~days. For more precise solar neutrino measurements, several improvements are applied in this analysis: lowering th… ▽ More An analysis of solar neutrino data from the fourth phase of Super-Kamiokande~(SK-IV) from October 2008 to May 2018 is performed and the results are presented. The observation time of the data set of SK-IV corresponds to $2970$~days and the total live time for all four phases is $5805$~days. For more precise solar neutrino measurements, several improvements are applied in this analysis: lowering the data acquisition threshold in May 2015, further reduction of the spallation background using neutron clustering events, precise energy reconstruction considering the time variation of the PMT gain. The observed number of solar neutrino events in $3.49$--$19.49$ MeV electron kinetic energy region during SK-IV is $65,443^{+390}_{-388}\,(\mathrm{stat.})\pm 925\,(\mathrm{syst.})$ events. Corresponding $\mathrm{^{8}B}$ solar neutrino flux is $(2.314 \pm 0.014\, \rm{(stat.)} \pm 0.040 \, \rm{(syst.)}) \times 10^{6}~\mathrm{cm^{-2}\,s^{-1}}$, assuming a pure electron-neutrino flavor component without neutrino oscillations. The flux combined with all SK phases up to SK-IV is $(2.336 \pm 0.011\, \rm{(stat.)} \pm 0.043 \, \rm{(syst.)}) \times 10^{6}~\mathrm{cm^{-2}\,s^{-1}}$. Based on the neutrino oscillation analysis from all solar experiments, including the SK $5805$~days data set, the best-fit neutrino oscillation parameters are $\rm{sin^{2} θ_{12,\,solar}} = 0.306 \pm 0.013 $ and $Δm^{2}_{21,\,\mathrm{solar}} = (6.10^{+ 0.95}_{-0.81}) \times 10^{-5}~\rm{eV}^{2}$, with a deviation of about 1.5$σ$ from the $Δm^{2}_{21}$ parameter obtained by KamLAND. The best-fit neutrino oscillation parameters obtained from all solar experiments and KamLAND are $\sin^{2} θ_{12,\,\mathrm{global}} = 0.307 \pm 0.012 $ and $Δm^{2}_{21,\,\mathrm{global}} = (7.50^{+ 0.19}_{-0.18}) \times 10^{-5}~\rm{eV}^{2}$. △ Less

Submitted 20 February, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 47 pages, 61 figures

Journal ref: Phys. Rev. D 109, 092001 (2024)

Constraining nucleon effective masses with flow and stopping observables from the S$π$RIT experiment

Authors: C. Y. Tsang, M. Kurata-Nishimura, M. B. Tsang, W. G. Lynch, Y. X. Zhang, J. Barney, J. Estee, G. Jhang, R. Wang, M. Kaneko, J. W. Lee, T. Isobe, T. Murakami, D. S. Ahn, L. Atar, T. Aumann, H. Baba, K. Boretzky, J. Brzychczyk, G. Cerizza, N. Chiga, N. Fukuda, I. Gasparic, B. Hong, A. Horvat , et al. (30 additional authors not shown)

Abstract: Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow ($v_1$), elliptic flow ($v_2$) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the S$π$RIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultane… ▽ More Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow ($v_1$), elliptic flow ($v_2$) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the S$π$RIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultaneously within the Improved Quantum Molecular Dynamics-Skyrme (ImQMD-Sky) transport code to obtain a multivariate correlated constraint. The varied parameters include symmetry energy, $S_0$, and slope of the symmetry energy, $L$, at saturation density, isoscalar effective mass, $m_{s}^*/m_{N}$, isovector effective mass, $m_{v}^{*}/m_{N}$ and the in-medium cross-section enhancement factor $η$. We find that the flow and VarXZ observables are sensitive to the splitting of proton and neutron effective masses and the in-medium cross-section. Comparisons of ImQMD-Sky predictions to the S$π$RIT data suggest a narrow range of preferred values for $m_{s}^*/m_{N}$, $m_{v}^{*}/m_{N}$ and $η$. △ Less

Submitted 8 December, 2023; originally announced December 2023.

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.

Comments: Accepted for publication in MNRAS; 15 pages, 13 figures

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.

Comments: 21 pages, 15 figures, 3 tables, accepted for publication in A&A. Example galaxies in Figs. 2, 6 and 7. Key results in Figs. 8, 9, 13 and 14

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 3 June, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 19 pages, 13 figures, 3 tables, accepted for publication by The Open Journal of Astrophysics, code available at https://github.com/astropy/halotools

Journal ref: Open Journal of Astrophysics, 7 (2024)

Atmospheric neutrino oscillation analysis with neutron tagging and an expanded fiducial volume in Super-Kamiokande I-V

Authors: Super-Kamiokande Collaboration, :, T. Wester, K. Abe, C. Bronner, Y. Hayato, K. Hiraide, K. Hosokawa, K. Ieki, M. Ikeda, J. Kameda, Y. Kanemura, R. Kaneshima, Y. Kashiwagi, Y. Kataoka, S. Miki, S. Mine, M. Miura, S. Moriyama, Y. Nakano, M. Nakahata, S. Nakayama, Y. Noguchi, K. Sato, H. Sekiya , et al. (212 additional authors not shown)

Abstract: We present a measurement of neutrino oscillation parameters with the Super-Kamiokande detector using atmospheric neutrinos from the complete pure-water SK I-V (April 1996-July 2020) data set, including events from an expanded fiducial volume. The data set corresponds to 6511.3 live days and an exposure of 484.2 kiloton-years. Measurements of the neutrino oscillation parameters $Δm^2_{32}$,… ▽ More We present a measurement of neutrino oscillation parameters with the Super-Kamiokande detector using atmospheric neutrinos from the complete pure-water SK I-V (April 1996-July 2020) data set, including events from an expanded fiducial volume. The data set corresponds to 6511.3 live days and an exposure of 484.2 kiloton-years. Measurements of the neutrino oscillation parameters $Δm^2_{32}$, $\sin^2θ_{23}$, $\sin^2 θ_{13}$, $δ_{CP}$, and the preference for the neutrino mass ordering are presented with atmospheric neutrino data alone, and with constraints on $\sin^2 θ_{13}$ from reactor neutrino experiments. Our analysis including constraints on $\sin^2 θ_{13}$ favors the normal mass ordering at the 92.3% level. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 24 pages, 18 figures

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.

Comments: 14 pages, 7 figures

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.

Comments: 38 pages, 6 figures. This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, Nat Astron (2023), and is available online at https://doi.org/10.1038/s41550-023-02052-3

Measurement of the muon flux at the SND@LHC experiment

Authors: SND@LHC Collaboration

Abstract: The Scattering and Neutrino Detector at the LHC (\SND) started taking data at the beginning of Run 3 of the LHC. The experiment is designed to perform measurements with neutrinos produced in proton-proton collisions at the LHC in an energy range between 100GeV and 1 TeV. It covers a previously unexplored pseudo-rapidity range of $7.2<η<8.4$. The detector is located 480 m downstream of the ATLAS in… ▽ More The Scattering and Neutrino Detector at the LHC (\SND) started taking data at the beginning of Run 3 of the LHC. The experiment is designed to perform measurements with neutrinos produced in proton-proton collisions at the LHC in an energy range between 100GeV and 1 TeV. It covers a previously unexplored pseudo-rapidity range of $7.2<η<8.4$. The detector is located 480 m downstream of the ATLAS interaction point in the TI18 tunnel. It comprises a veto system, a target consisting of tungsten plates interleaved with nuclear emulsion and scintillating fiber (SciFi) trackers, followed by a muon detector (UpStream, US and DownStream, DS). In this article we report the measurement of the muon flux in three subdetectors: the emulsion, the SciFi trackers and the DownStream Muon detector. The muon flux per integrated luminosity through an 18$\times$18 cm$^{2}$ area in the emulsion is $1.5 \pm 0.1(\textrm{stat}) \times 10^4\,\textrm{fb/cm}^{2}$. The muon flux per integrated luminosity through a 31$\times$31 cm$^{2}$ area in the centre of the SciFi is $2.06\pm0.01(\textrm{stat})\pm0.12(\textrm{sys}) \times 10^{4} \textrm{fb/cm}^{2}$. The muon flux per integrated luminosity through a 52$\times$52 cm$^{2}$ area in the centre of the downstream muon system is $2.35\pm0.01(\textrm{stat})\pm0.10(\textrm{sys}) \times 10^{4}\,\textrm{fb/cm}^{2}$. The total relative uncertainty of the measurements by the electronic detectors is 6 $\%$ for the SciFi and 4 $\%$ for the DS measurement. The Monte Carlo simulation prediction of these fluxes is 20-25 $\%$ lower than the measured values. △ Less

Submitted 16 November, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 11 pages, 14 figures

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.

Comments: 12 pages, 11 figures

Star formation in CALIFA survey perturbed galaxies.III. Stellar & ionized-gas kinematic distributions

Authors: Abdías Morales-Vargas, Juan Pablo Torres-Papaqui, Fernando Fabián Rosales-Ortega, Marcel Chow-Martínez, René Alberto Ortega-Minakata, Aitor Carlos Robleto-Orús, the CALIFA survey Collaboration

Abstract: We obtain the kinematic distributions of stars (synthetic model line absorption) and ionized gas (H$α$ line emission) for star-forming regions residing in CALIFA survey tidally perturbed (perturbed) and non-tidally perturbed (control) galaxies. We set the uncertainties of the velocity dispersion by measuring the statistical variability of the datasets themselves. Using these adopted uncertainties… ▽ More We obtain the kinematic distributions of stars (synthetic model line absorption) and ionized gas (H$α$ line emission) for star-forming regions residing in CALIFA survey tidally perturbed (perturbed) and non-tidally perturbed (control) galaxies. We set the uncertainties of the velocity dispersion by measuring the statistical variability of the datasets themselves. Using these adopted uncertainties and considering the sensitivity of the grating device, we establish thresholds of reliability that allow us to select reliable velocity dispersions. From this selection, we pair the star-forming spaxels between control and perturbed galaxies at the closest shifts in velocity (de-redshifting). We compare their respective distributions of velocity dispersion. In perturbed galaxies, median velocity dispersions for the stellar and gaseous components are minimally higher and equal, respectively, than those in control galaxies. The spread in velocity dispersion and the velocity shift - velocity dispersion space agree with this result. Unlike the well-known trend in strongly interacting systems, the stellar and ionized-gas motions are not disturbed by the influence of close companions. For the gaseous component, this result is due to the poor statistical variability of its data, a consequence of the tightness in velocity dispersion derived from high spectral line intensities. This analysis concludes the series, which previously showed star-forming regions in galaxies with close companions undergoing more prominent gas inflows, resulting in differences in their star formation and consequent metal content. △ Less

Submitted 17 September, 2023; originally announced September 2023.

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 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.

Magnetic fields & rotation periods of M dwarfs from SPIRou spectra

Authors: J. -F. Donati, L. T. Lehmann, P. I. Cristofari, P. Fouqué, C. Moutou, P. Charpentier, M. Ould-Elhkim, A. Carmona, X. Delfosse, E. Artigau, S. H. P. Alencar, C. Cadieux, L. Arnold, P. Petit, J. Morin, T. Forveille, R. Cloutier, R. Doyon, G. Hébrard, the SLS collaboration

Abstract: We present near-infrared spectropolarimetric observations of a sample of 43 weakly- to moderately-active M dwarfs, carried with SPIRou at the Canada-France-Hawaii Telescope in the framework of the SPIRou Legacy Survey from early 2019 to mid 2022. We use the 6700 circularly polarised spectra collected for this sample to investigate the longitudinal magnetic field and its temporal variations for all… ▽ More We present near-infrared spectropolarimetric observations of a sample of 43 weakly- to moderately-active M dwarfs, carried with SPIRou at the Canada-France-Hawaii Telescope in the framework of the SPIRou Legacy Survey from early 2019 to mid 2022. We use the 6700 circularly polarised spectra collected for this sample to investigate the longitudinal magnetic field and its temporal variations for all sample stars, from which we diagnose, through quasi-periodic Gaussian process regression, the periodic modulation and longer-term fluctuations of the longitudinal field. We detect the large-scale field for 40 of our 43 sample stars, and infer a reliable or tentative rotation period for 38 of them, using a Bayesian framework to diagnose the confidence level at which each rotation period is detected. We find rotation periods ranging from 14 to over 60d for the early-M dwarfs, and from 70 to 200d for most mid- and late-M dwarfs (potentially up to 430d for one of them). We also find that the strength of the detected large-scale fields does not decrease with increasing period or Rossby number for the slowly rotating dwarfs of our sample as it does for higher-mass, more active stars, suggesting that these magnetic fields may be generated through a different dynamo regime than those of more rapidly rotating stars. We also show that the large-scale fields of most sample stars evolve on long timescales, with some of them globally switching sign as stars progress on their putative magnetic cycles. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: MNRAS, in press (25 pages, 15 figures, 3 tables)

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.

Comments: 30 pages, 36 figures, published version in A&A

Journal ref: A&A, 684, A21 (2024)

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.

Comments: 13 pages, 7 figures. Accepted to JATIS

Journal ref: J. Astron. Telesc. Instrum. Syst. 9(2), 028006 (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.

Journal ref: NIM B 541 (2023) 275-278

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.

Comments: 6 pages, 5 figures, 2 tables, 7th Heidelberg International Symposium on High-Energy Gamma-Ray Astronomy (Gamma2022), pre-published in Proceedings of science, Speaker P. Chambéry

Journal ref: PoS, 417, 2023, 135

GECKOS: Turning galaxy evolution on its side with deep observations of edge-on galaxies

Authors: J. van de Sande, A. Fraser-McKelvie, D. B. Fisher, M. Martig, M. R. Hayden, the GECKOS Survey collaboration

Abstract: We present GECKOS (Generalising Edge-on galaxies and their Chemical bimodalities, Kinematics, and Outflows out to Solar environments), a new ESO VLT/MUSE large program. The main aim of GECKOS is to reveal the variation in key physical processes of disk formation by connecting Galactic Archaeology with integral field spectroscopic observations of nearby galaxies. Edge-on galaxies are ideal for this… ▽ More We present GECKOS (Generalising Edge-on galaxies and their Chemical bimodalities, Kinematics, and Outflows out to Solar environments), a new ESO VLT/MUSE large program. The main aim of GECKOS is to reveal the variation in key physical processes of disk formation by connecting Galactic Archaeology with integral field spectroscopic observations of nearby galaxies. Edge-on galaxies are ideal for this task: they allow us to disentangle the assembly history imprinted in thick disks and provide the greatest insights into outflows. The GECKOS sample of 35 nearby edge-on disk galaxies is designed to trace the assembly histories and properties of galaxies across a large range of star formation rates, bulge-to-total ratios, and boxy and non-boxy bulges. GECKOS will deliver spatially resolved measurements of stellar abundances, ages, and kinematics, as well as ionised gas metallicities, ionisation parameters, pressure, and inflow and outflow kinematics; all key parameters for building a complete chemodynamical picture of disk galaxies. With these data, we aim to extend Galactic analysis methods to the wider galaxy population, reaping the benefits of detailed Milky Way studies, while probing the diverse mechanisms of galaxy evolution. △ Less

Submitted 31 May, 2023; originally announced June 2023.

Comments: 6 pages, 2 figures, Proceedings of IAUS 377, eds. F. Tabatabaei, B. Barbuy, and Y. Ting

DES Y3 + KiDS-1000: Consistent cosmology combining cosmic shear surveys

Authors: Dark Energy Survey, Kilo-Degree Survey Collaboration, :, T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, M. Asgari, S. Avila, D. Bacon, K. Bechtol, M. R. Becker, G. M. Bernstein, E. Bertin, M. Bilicki, J. Blazek, S. Bocquet, D. Brooks, P. Burger, D. L. Burke, H. Camacho, A. Campos, A. Carnero Rosell , et al. (138 additional authors not shown)

Abstract: We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of… ▽ More We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of $0.790^{+0.018}_{-0.014}$. The mean marginal is lower than the maximum a posteriori estimate, $S_8=0.801$, owing to skewness in the marginal distribution and projection effects in the multi-dimensional parameter space. Our results are consistent with $S_8$ constraints from observations of the cosmic microwave background by Planck, with agreement at the $1.7σ$ level. We use a Hybrid analysis pipeline, defined from a mock survey study quantifying the impact of the different analysis choices originally adopted by each survey team. We review intrinsic alignment models, baryon feedback mitigation strategies, priors, samplers and models of the non-linear matter power spectrum. △ Less

Submitted 19 October, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 21 figures, 15 tables, accepted Open Journal of Astrophysics. Download the chains from https://des.ncsa.illinois.edu/releases/y3a2/Y3key-joint-des-kids or create your own chains with CosmoSIS using https://github.com/joezuntz/cosmosis-standard-library/blob/main/examples/des-y3_and_kids-1000.ini Watch the core team discuss this analysis at https://cosmologytalks.com/2023/05/26/des-kids

Report number: FERMILAB-PUB-23-267-PPD

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.

Comments: 9 pages, 6 figures; submitted to A&A

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.

Comments: 15 pages, 7 figures. Accepted for publication in A&A. Corresponding authors: Victor Barbosa Martins, Stefan Ohm, Cornelia Arcaro, Natalia Żywucka, Mathieu de Naurois

Journal ref: A&A 675, A138 (2023)

Observation of collider muon neutrinos with the SND@LHC experiment

Authors: The SND@LHC Collaboration

Abstract: We report the direct observation of muon neutrino interactions with the SND@LHC detector at the Large Hadron Collider. A data set of proton-proton collisions at $\sqrt{s} = 13.6\,$TeV collected by SND@LHC in 2022 is used, corresponding to an integrated luminosity of 36.8$\,\rm{fb}^{-1}$. The search is based on information from the active electronic components of the SND@LHC detector, which covers… ▽ More We report the direct observation of muon neutrino interactions with the SND@LHC detector at the Large Hadron Collider. A data set of proton-proton collisions at $\sqrt{s} = 13.6\,$TeV collected by SND@LHC in 2022 is used, corresponding to an integrated luminosity of 36.8$\,\rm{fb}^{-1}$. The search is based on information from the active electronic components of the SND@LHC detector, which covers the pseudo-rapidity region of $7.2 < η< 8.4$, inaccessible to the other experiments at the collider. Muon neutrino candidates are identified through their charged-current interaction topology, with a track propagating through the entire length of the muon detector. After selection cuts, 8 $ν_μ$ interaction candidate events remain with an estimated background of 0.076 events, yielding a significance of seven standard deviations for the observed $ν_μ$ signal. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: Submitted to Phys. Rev. Lett

Report number: CERN-EP-2023-092

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.

Comments: Submitted to A&A, 18 pages, 14 figures, comments are welcome

Journal ref: A&A 679, A61 (2023)

The magnetic field and multiple planets of the young dwarf AU~Mic

Authors: J. -F. Donati, P. I. Cristofari, B. Finociety, B. Klein, C. Moutou, E. Gaidos, C. Cadieux, E. Artigau, A. C. M. Correia, G. Boué, N. J. Cook, A. Carmona, L. T. Lehmann, J. Bouvier, E. Martioli, J. Morin, P. Fouqué, X. Delfosse, R. Royon, G. Hébrard, S. H. P. Alencar, J. Laskar, L. Arnold, P. Petit, A. Kospal , et al. (3 additional authors not shown)

Abstract: In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures o… ▽ More In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures of spectral lines. We find that both are modulated with the stellar rotation period (4.86 d), and evolve on a timescale of months under differential rotation and intrinsic variability. The small-scale field, estimated from the broadening of spectral lines, reaches $2.61\pm0.05$ kG. The large-scale field, inferred with Zeeman-Doppler imaging from Least-Squares Deconvolved profiles of circularly-polarized and unpolarized spectral lines, is mostly poloidal and axisymmetric, with an average intensity of $550\pm30$ G. We also find that surface differential rotation, as derived from the large-scale field, is $\simeq$30% weaker than that of the Sun. We detect the radial velocity (RV) signatures of transiting planets b and c, although dwarfed by activity, and put an upper limit on that of candidate planet d, putatively causing the transit-timing variations of b and c. We also report the detection of the RV signature of a new candidate planet (e) orbiting further out with a period of $33.39\pm0.10$ d, i.e., near the 4:1 resonance with b. The RV signature of e is detected at 6.5$σ$ while those of b and c show up at $\simeq$4$σ$, yielding masses of $10.2^{+3.9}_{-2.7}$ and $14.2^{+4.8}_{-3.5}$ Earth masses for b and c, and a minimum mass of $35.2^{+6.7}_{-5.4}$ Earth masses for e. △ Less

Submitted 24 April, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: MNRAS, in press (20 pages and 12 figures + 9 pages of supplementary material)

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.

Comments: 16 pages, 15 figures, 7 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 673, A148 (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.

Comments: Accepted for publication in JCAP, corresponding authors: F.Brun, J-F. Glicenstein, V. Marandon, T. Tavernier

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.

Comments: 10 pages, 4 figures. Accepted for publication in APJL. Corresponding authors: J. Damascene Mbarubucyeye, H. Ashkar, S. J. Zhu, B. Reville, F. Schüssler

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.

Comments: 10 pages, 5 figures; minor text updates and additional references

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.

Comments: 14 pages, 10 figures. Accepted for publication in A&A. Corresponding authors: Vikas Joshi, Lars Mohrmann

Journal ref: A&A 672, A103 (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.

Measurement of the cosmogenic neutron yield in Super-Kamiokande with gadolinium loaded water

Authors: Super-Kamiokande Collaboration, :, M. Shinoki, K. Abe, Y. Hayato, K. Hiraide, K. Hosokawa, K. Ieki, M. Ikeda, J. Kameda, Y. Kanemura, R. Kaneshima, Y. Kashiwagi, Y. Kataoka, S. Miki, S. Mine, M. Miura, S. Moriyama, Y. Nakano, M. Nakahata, S. Nakayama, Y. Noguchi, K. Okamoto, K. Sato, H. Sekiya , et al. (217 additional authors not shown)

Abstract: Cosmic-ray muons that enter the Super-Kamiokande detector cause hadronic showers due to spallation in water, producing neutrons and radioactive isotopes. Those are a major background source for studies of MeV-scale neutrinos and searches for rare events. Since 2020, gadolinium was introduced in the ultra-pure water in the Super-Kamiokande detector to improve the detection efficiency of neutrons. I… ▽ More Cosmic-ray muons that enter the Super-Kamiokande detector cause hadronic showers due to spallation in water, producing neutrons and radioactive isotopes. Those are a major background source for studies of MeV-scale neutrinos and searches for rare events. Since 2020, gadolinium was introduced in the ultra-pure water in the Super-Kamiokande detector to improve the detection efficiency of neutrons. In this study, the cosmogenic neutron yield was measured using data acquired during the period after the gadolinium loading. The yield was found to be $(2.76 \pm 0.02\,\mathrm{(stat.) \pm 0.19\,\mathrm{(syst.)}}) \times 10^{-4}\,μ^{-1} \mathrm{g^{-1} cm^{2}}$ at 259 GeV of average muon energy at the Super-Kamiokande detector. △ Less

Submitted 25 October, 2023; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: 10 pages, 10 figures, 3 tables

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.

Comments: 20+11 pages, 10+11 figures. Version accepted in The Open Journal of Astrophysics

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.

Comments: 18 pages, 14 figures

Journal ref: ApJ 950 49 (2023)

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.

Comments: 15 pages, 8 figures. Journal-accepted version

Fringing Analysis and Simulation for the Vera C. Rubin Observatory's Legacy Survey of Space and Time

Authors: Zhiyuan Guo, Chris W. Walter, Craig Lage, Robert H. Lupton, The LSST Dark Energy Science Collaboration

Abstract: The presence of fringing in astronomical CCD images will impact photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for Charge-Coupled Devices (CCDs) already implemented on the Rubin Observatory LSST Camera's focal plane. After making physical measurements… ▽ More The presence of fringing in astronomical CCD images will impact photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for Charge-Coupled Devices (CCDs) already implemented on the Rubin Observatory LSST Camera's focal plane. After making physical measurements and knowing the compositions, we have developed a model for the e2v CCDs. We present a method to fit for the internal height variation of the epoxy layer within the sensors based on fringing measurements in a laboratory setting. This method is generic enough that it can be easily modified to work for other CCDs. Using the derived fringing model, we successfully reproduce comparable fringing amplitudes that match the observed levels in images taken by existing telescopes with different optical designs. This model is then used to forecast the expected level of fringing in a single LSST y-band sky background exposure with Rubin telescope optics in the presence of a realistic time varying sky spectrum. The predicted fringing amplitude in LSST images ranges from $0.04\%$ to $0.2\%$ depending on the location of a CCD on the focal plane. We find that the predicted variation in surface brightness caused by fringing in LSST y-band skybackground images is about $0.6\ μ\rm{Jy}\ \rm{arcsec}^{-2}$, which is 40 times larger than the current measurement error. We conclude that it is necessary to include fringing correction in the Rubin's LSST image processing pipeline. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 17 pages, 16 figures, submitted to PASP