Search

With GW170817 being the only multimessenger gravitational wave (GW) event with an associated kilonova (KN) detected so far, there exists a pressing need for realistic estimation of the GW localization uncertainties and rates, as well as optimization of available telescope time to enable the detection of new KNe. For this purpose, we simulate GW events assuming a data-driven, GW-motivated distribution of binary parameters for the LIGO/Virgo/KAGRA (LVK) fourth and fifth observing runs (O4 and O5). For our particular case of estimating KN detection rates in the O4 and O5 runs, we map the binary neutron star (BNS) and neutron star-black hole (NSBH) properties to the optical light curves arising from r-process nucleosynthesis in the ejecta. We use the simulated population of KNe to generate follow-up observing plans, with the primary goal of optimizing detection with the Gravitational Wave Multi-Messenger Astronomy DECam Survey (GW-MMADS). We explore KN detectability as a function of mass, distance, and spin of the binaries. Out of the mergers that produce a KN in our simulations, we expect detectable KNe for DECam-like instruments at a per-year rate of: $1-16$ ($0-1$) for BNS (NSBH) in O4, and $16-165$ ($1-12$) for BNS (NSBH) in O5, using a fiducial exposure time and conditional on the uncertainty on the equation of state (EOS) and volumetric rates of the mergers. Taking into account observability constraints, our scheduler covers the location of the KN $\sim 30-38\%$ of the times for our fiducial EOS. We provide the depths needed to detect a significant fraction of our simulated mergers for the astronomical community to use in their follow-up campaigns., Comment: 18 pages, 10 figures, 10 tables

The Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory is planned to begin in the Fall of 2025. The LSST survey cadence has been designed via a community-driven process regulated by the Survey Cadence Optimization Committee (SCOC), which recommended up to 3% of the observing time to carry out Target of Opportunity (ToO) observations. Experts from the scientific community, Rubin Observatory personnel, and members of the SCOC were brought together to deliver a recommendation for the implementation of the ToO program during a workshop held in March 2024. Four main science cases were identified: gravitational wave multi-messenger astronomy, high energy neutrinos, Galactic supernovae, and small potentially hazardous asteroids possible impactors. Additional science cases were identified and briefly addressed in the documents, including lensed or poorly localized gamma-ray bursts and twilight discoveries. Trigger prioritization, automated response, and detailed strategies were discussed for each science case. This document represents the outcome of the Rubin ToO 2024 workshop, with additional contributions from members of the Rubin Science Collaborations. The implementation of the selection criteria and strategies presented in this document has been endorsed in the SCOC Phase 3 Recommendations document (PSTN-056). Although the ToO program is still to be finalized, this document serves as a baseline plan for ToO observations with the Rubin Observatory.

Modeling the intrinsic alignment (IA) of galaxies poses a challenge to weak lensing analyses. The Dark Energy Survey is expected to be less impacted by IA when limited to blue, star-forming galaxies. The cosmological parameter constraints from this blue cosmic shear sample are stable to IA model choice, unlike passive galaxies in the full DES Y3 sample, the goodness-of-fit is improved and the $\Omega_{m}$ and $S_8$ better agree with the cosmic microwave background. Mitigating IA with sample selection, instead of flexible model choices, can reduce uncertainty in $S_8$ by a factor of 1.5., Comment: Data access available at https://jamiemccullough.github.io/data/blueshear/

We present the results of our multi-wavelength (X-ray to radio) follow-up campaign of the Einstein Probe transient EP240408a. The initial 10 s trigger displayed bright soft X-ray (0.5-4 keV) radiation with peak luminosity $L_\textrm{X} \gtrsim 10^{49}$ ($10^{50}$) erg s$^{-1}$ for an assumed redshift z>0.5 (2.0). The Neil Gehrels Swift Observatory and Neutron star Interior Composition ExploreR discovered a fading X-ray counterpart lasting for $\sim$5 d (observer frame), which showed a long-lived (~4 d) plateau-like emission ($t^{-0.5}$) before a sharp powerlaw decline ($t^{-7}$). The plateau emission was in excess of $L_\textrm{X} \gtrsim 10^{46}$ ($10^{47}$) erg s$^{-1}$ at z>0.5 (2.0). Deep optical and radio observations resulted in non-detections of the transient. Our observations with Gemini South revealed a faint potential host galaxy ($r \approx 24$ AB mag) near the edge of the X-ray localization. The faint candidate host, and lack of other potential hosts ($r \gtrsim 26$ AB mag; $J \gtrsim 23$ AB mag), implies a higher redshift origin (z>0.5), which produces extreme X-ray properties that are inconsistent with many known extragalactic transient classes. In particular, the lack of a bright gamma-ray counterpart, with the isotropic-equivalent energy ($10 - 10,000$ keV) constrained by GECam and Konus-Wind to $E_{\gamma,\textrm{iso}} \lesssim 4\times10^{51}$ ($6\times10^{52}$) erg at z>0.5 (2.0), conflicts with known gamma-ray bursts (GRBs) of similar X-ray luminosities. We therefore favor a jetted tidal disruption event (TDE) as the progenitor of EP240408a at z>1.0, possibly caused by the disruption of a white dwarf by an intermediate mass black hole. The alternative is that EP240408a may represent a new, previously unknown class of transient., Comment: Submitted; 33 pages; 9 figures

Following the wealth of new results enabled by multimessenger observations of the binary neutron star (BNS) merger GW170817, the next goal is increasing the number of detections of electromagnetic (EM) counterparts to gravitational wave (GW) events. Here, we study the detectability of the prompt emission and afterglows produced by the relativistic jets launched by BNS mergers that will be detected by LIGO-Virgo-KAGRA (LVK) during their fifth observing run (O5), and by next generation (XG) GW detectors (Einstein Telescope and Cosmic Explorer). We quantify the impact of various BNS merger and jet afterglow parameters on the likelihood of detection for a wide range of telescopes, focusing on the impact of the observer's viewing angle $\theta_\textrm{v}$ and the jet's core half-opening angle $\theta_\textrm{c}$. We find that during the LVK O5 run, the James Webb Space Telescope (JWST) is expected to reach the largest number of afterglow detections, up to $32\%$ of BNS mergers with $27\%$ detectable on timescales $>10$ d. Overall the detection of $20-30\%$ of events is possible across the EM spectrum with the available instruments, including the Chandra X-ray Observatory, the Nancy Grace Roman Space Telescope, the Square Kilometer Array (SKA), and the Very Large Array (VLA). In the XG era, afterglows for $20-25\%$ of events will be detectable with AXIS, Athena, Lynx, UVEX, JWST and the next generation VLA (ngVLA), reaching beyond redshift $z \gtrsim 1$. We also find that the majority of detected afterglows are expected to accompany prompt gamma-ray emission detectable by the Neil Gehrels Swift Observatory, the Space-based multi-band astronomical Variable Objects Monitor (SVOM), and the Transient High-Energy Sky and Early Universe Surveyor (THESEUS), with at least $\sim 6\%$ of afterglows being orphan., Comment: Submitted to Physical Review D, 36 pages, 13 figures

Gravitational lensing magnification of Type Ia supernovae (SNe Ia) allows information to be obtained about the distribution of matter on small scales. In this paper, we derive limits on the fraction $\alpha$ of the total matter density in compact objects (which comprise stars, stellar remnants, small stellar groupings and primordial black holes) of mass $M > 0.03 M_{\odot}$ over cosmological distances. Using 1,532 SNe Ia from the Dark Energy Survey Year 5 sample (DES-SN5YR) combined with a Bayesian prior for the absolute magnitude $M$, we obtain $\alpha < 0.12$ at the 95\% confidence level after marginalisation over cosmological parameters, lensing due to large-scale structure, and intrinsic non-Gaussianity. Similar results are obtained using priors from the cosmic microwave background, baryon acoustic oscillations and galaxy weak lensing, indicating our results do not depend on the background cosmology. We argue our constraints are likely to be conservative (in the sense of the values we quote being higher than the truth), but discuss scenarios in which they could be weakened by systematics of the order of $\Delta \alpha \sim 0.04$, Comment: Accepted by MNRAS

The detection of the gravitational-wave event GW230529, presumably a neutron star-black hole (NSBH) merger, by the LIGO-Virgo-KAGRA (LVK) Collaboration is an exciting discovery for multimessenger astronomy. The black hole (BH) has a high probability of falling within the ''mass gap'' between the peaks of the neutron star (NS) and the BH mass distributions. Because of the low primary mass, the binary is more likely to produce an electromagnetic counterpart than previously detected NSBH mergers. We investigate the possible kilonova (KN) emission from GW230529, and find that if it was an NSBH, there is a $\sim$ 2-41% probability (depending on the assumed equation of state) that GW230925 produced a KN with magnitude peaking at $\sim 1-2$ day post merger at $g \lesssim 23.5$, $i<23$. Hence, it could have been detected by ground-based telescopes. If it was a binary neutron star (BNS) merger, we find $\sim$ 0-12% probability that it produced a KN. Motivated by these numbers, we simulated a broader population of mgNSBH mergers that may be detected in O4, and we obtained a 9-21% chance of producing a KN, which would be detectable with $g\lesssim 25$ and $ i \lesssim 24$, typically fainter than what is expected from GW230529. Based on these findings, DECam-like instruments may be able to detect up to 80% of future mgNSBH KNe, thus up to $\sim1$ multimessenger mgNSBH per year may be discoverable at the current level of sensitivity (O4)., Comment: 13 pages, 4 figures, 2 tables

We present the detection and analysis of GRB 231115A, a candidate extragalactic magnetar giant flare (MGF) observed by Fermi/GBM and localized by INTEGRAL to the starburst galaxy M82. This burst exhibits distinctive temporal and spectral characteristics that align with known MGFs, including a short duration and a high peak energy. Gamma-ray analyses reveal significant insights into this burst, supporting conclusions already established in the literature: our time-resolved spectral studies provide further evidence that GRB 231115A is indeed a MGF. Significance calculations also suggest a robust association with M82, further supported by a high Bayes factor that minimizes the probability of chance alignment with a neutron star merger. Despite extensive follow-up efforts, no contemporaneous gravitational wave or radio emissions were detected. The lack of radio emission sets stringent upper limits on possible radio luminosity. Constraints from our analysis show no fast radio bursts (FRBs) associated with two MGFs. X-ray observations conducted post-burst by Swift/XRT and XMM/Newton provided additional data, though no persistent counterparts were identified. Our study underscores the importance of coordinated multi-wavelength follow-up and highlights the potential of MGFs to enhance our understanding of short GRBs and magnetar activities in the cosmos. Current MGF identification and follow-up implementation are insufficient for detecting expected counterparts; however, improvements in these areas may allow for the recovery of follow-up signals with existing instruments. Future advancements in observational technologies and methodologies will be crucial in furthering these studies.

Using a sample of 2800 galaxy clusters identified in the Dark Energy Survey across the redshift range $0.20 < z < 0.60$, we characterize the hierarchical assembly of Bright Central Galaxies (BCGs) and the surrounding intracluster light (ICL). To quantify hierarchical formation we use the stellar mass - halo mass (SMHM) relation for the BCG+ICL system and incorporate the magnitude gap (M14), the difference in brightness between the BCG (measured within 30kpc) and 4th brightest cluster member galaxy within 0.5 $R_{200,c}$. The inclusion of M14, which traces BCG hierarchical growth, increases the slope and decreases the intrinsic scatter in the SMHM relation, highlighting that it is a latent variable within the BCG+ICL SMHM relation. Moreover, the correlation with M14 decreases at large radii from the BCG's centre. However, the stellar light within the BCG+ICL transition region (30kpc - 80kpc) most strongly correlates with the dark matter halo mass and has a statistically significant correlation with M14. As the light in the transition region and M14 are independent measurements, the transition region may grow as a result of the BCG's hierarchical two-phase formation. Additionally, as M14 and ICL result from hierarchical growth, we use a stacked sample and find that clusters with large M14 values are characterized by larger ICL and BCG+ICL fractions, which illustrates that the merger processes that build the BCG stellar mass also grow the ICL. Furthermore, this may suggest that M14 combined with the ICL fraction can be used as a method to identify dynamically relaxed clusters., Comment: 16 pages, 5 Figures, submitted to MNRAS on 8/30/2024

Using 1533 type Ia supernovae (SNe Ia) from the five-year sample of the Dark Energy Survey (DES), we investigate the effects of projected galactocentric separation between the SNe and their host galaxies on their light curves and standardization. We show, for the first time, that the difference in SN Ia post-standardization brightnesses between high and low-mass hosts reduces from $0.078\pm0.011$ mag in the full sample to $0.036 \pm 0.018$ mag for SNe Ia located in the outer regions of their host galaxies, while increasing to $0.100 \pm 0.014$ mag for SNe in the inner regions. In these inner regions, the step can be reduced (but not removed) using a model where the $R_V$ of dust along the line-of-sight to the SN changes as a function of galaxy properties. To explain the remaining difference, we use the distributions of the SN Ia stretch parameter to test whether the inferred age of SN progenitors are more varied in the inner regions of galaxies. We find that the proportion of high-stretch SNe Ia in red (older) environments is more prominent in outer regions and that the outer regions stretch distributions are overall more homogeneous compared to inner regions, but conclude that this effect cannot explain the reduction in significance of any Hubble residual step in outer regions. We conclude that the standardized distances of SNe Ia located in the outer regions of galaxies are less affected by their global host galaxy properties than those in the inner regions., Comment: 17 pages, 13 figures

The present state of cosmology is facing a crisis where there is a fundamental disagreement in measurements of the Hubble constant ($H_{0}$), with significant tension between the early and late universe methods. Type Ia supernovae (SNe Ia) are important to measuring $H_{0}$ through the astronomical distance ladder. However, there remains potential to better standardise SN Ia light curves by using known dependencies on host galaxy properties after the standard light curve width and colour corrections have been applied to the peak SN Ia luminosities. To explore this, we use the 5-year photometrically identified SNe Ia sample obtained by the Dark Energy Survey, along with host galaxy spectra obtained by the Australian Dark Energy Survey. Using host galaxy spectroscopy, we find a significant trend with the equivalent width (EW) of the [OII] $\lambda\lambda$ 3727, 29 doublet, a proxy for specific star formation rate, and Hubble residuals. We find that the correlation with [OII] EW is a powerful alternative to the commonly used mass step after initial light curve corrections. We applied our [OII] EW correction to a sample of 20 SN Ia hosted by calibrator galaxies observed using WiFeS, and examined the impact on both the SN Ia absolute magnitude and $H_{0}$. We then explored different [OII] EW corrections and found $H_{0}$ values ranging between $72.80$ to $73.28~\mathrm{km} \mathrm{s}^{-1} \mathrm{Mpc}^{-1}$. Notably, even after using an additional [OII] EW correction, the impact of host galaxy properties in standardising SNe Ia appears limited in reducing the current tension ($\sim$5$\sigma$) with the Cosmic Microwave Background result for $H_{0}$., Comment: 16 pages, 13 figures. Submitting to MNRAS

We present galaxy-galaxy lensing measurements using a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data as lenses. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. These dark-matter-dominated objects are intriguing due to potentially unusual formation channels that lead to their diffuse stellar component. Given the faintness of LSBGs, using standard observational techniques to characterize their total masses proves challenging. Weak gravitational lensing, which is less sensitive to the stellar component of galaxies, could be a promising avenue to estimate the masses of LSBGs. Our LSBG sample consists of 23,790 galaxies separated into red and blue color types at $g-i\ge 0.60$ and $g-i< 0.60$, respectively. Combined with the DES Y3 shear catalog, we measure the tangential shear around these LSBGs and find signal-to-noise ratios of 6.67 for the red sample, 2.17 for the blue sample, and 5.30 for the full sample. We use the clustering redshifts method to obtain redshift distributions for the red and blue LSBG samples. Assuming all red LSBGs are satellites, we fit a simple model to the measurements and estimate the host halo mass of these LSBGs to be $\log(M_{\rm host}/M_{\odot}) = 12.98 ^{+0.10}_{-0.11}$. We place a 95% upper bound on the subhalo mass at $\log(M_{\rm sub}/M_{\odot})<11.51$. By contrast, we assume the blue LSBGs are centrals, and place a 95% upper bound on the halo mass at $\log(M_\mathrm{host}/M_\odot) < 11.84$. We find that the stellar-to-halo mass ratio of the LSBG samples is consistent with that of the general galaxy population. This work illustrates the viability of using weak gravitational lensing to constrain the halo masses of LSBGs., Comment: 20 pages, 14 figures

Cosmological analyses with Type Ia Supernovae (SNe Ia) have traditionally been reliant on spectroscopy for both classifying the type of supernova and obtaining reliable redshifts to measure the distance-redshift relation. While obtaining a host-galaxy spectroscopic redshift for most SNe is feasible for small-area transient surveys, it will be too resource intensive for upcoming large-area surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time, which will observe on the order of millions of SNe. Here we use data from the Dark Energy Survey (DES) to address this problem with photometric redshifts (photo-z) inferred directly from the SN light-curve in combination with Gaussian and full p(z) priors from host-galaxy photo-z estimates. Using the DES 5-year photometrically-classified SN sample, we consider several photo-z algorithms as host-galaxy photo-z priors, including the Self-Organizing Map redshifts (SOMPZ), Bayesian Photometric Redshifts (BPZ), and Directional-Neighbourhood Fitting (DNF) redshift estimates employed in the DES 3x2 point analyses. With detailed catalog-level simulations of the DES 5-year sample, we find that the simulated w can be recovered within $\pm$0.02 when using SN+SOMPZ or DNF prior photo-z, smaller than the average statistical uncertainty for these samples of 0.03. With data, we obtain biases in w consistent with simulations within ~1$\sigma$ for three of the five photo-z variants. We further evaluate how photo-z systematics interplay with photometric classification and find classification introduces a subdominant systematic component. This work lays the foundation for next-generation fully photometric SNe Ia cosmological analyses., Comment: 19 pages, 9 figures. Submitting to MNRAS, comments welcome

We carried out long-term monitoring of the LIGO/Virgo/KAGRA binary black hole (BBH) merger candidate S230922g in search of electromagnetic emission from the interaction of the merger remnant with an embedding active galactic nuclei (AGN) accretion disk. Using a dataset primarily composed of wide-field imaging from the Dark Energy Camera (DECam) and supplemented by additional photometric and spectroscopic resources, we searched ~ 70% of the sky area probability for transient phenomena, and discovered 6 counterpart candidates. One especially promising candidate - AT 2023aagj - exhibited temporally varying asymmetric components in spectral broad line regions, a feature potentially indicative of an off-center event such as a BBH merger. This represents the first live search and multiwavelength, photometric, and spectroscopic monitoring of a GW BBH optical counterpart candidate in the disk of an AGN., Comment: 25 pages, 9 figures. Comments welcome

We analyze the confident binary black hole (BBH) detections from the third Gravitational-Wave Transient Catalog (GWTC-3) with an alternative mass population model in order to capture features in the mass distribution beyond the Powerlaw + Peak model. We find that the peak of a second power law characterizes the $\sim 30-35~ M_\odot$ bump, such that the data marginally prefers a mixture of two power laws for the mass distribution of binary components over a Powerlaw + Peak model with a Bayes Factor $\log_{10}\mathcal{B}$ of 0.1. This result may imply that the $\sim 30-35~ M_\odot$ feature represents the onset of a second population of BBH mergers (e.g. from a dynamical formation channel) rather than a specific mass feature over a broader distribution. When an additional Gaussian bump is allowed within our power law mixture model, we find a new feature in the BH mass spectrum at $\sim65-70~M_\odot$. This new feature may be consistent with hierarchical mergers, and constitute $\sim2\%$ of the BBH population. This model also recovers a maximum mass of $58^{+30}_{-14}~M_\odot$ for the second power law, consistent with the onset of a pair-instability supernova mass gap., Comment: Comments are welcomed

With the number of supernovae observed expected to drastically increase thanks to large-scale surveys like the Dark Energy Spectroscopic Instrument (DESI), it is necessary that the tools we use to classify these objects keep up with this increase. We previously created Supernova Tagging and Classification (STag) to address this problem by employing machine learning techniques alongside logistic regression in order to assign 'tags' to spectra based on spectral features. STag II is a continuation of this work, which now makes use of model supernova spectra combined with real DESI spectra in order to train STag to better deal with realistic data. We also make use of the rlap score as a trustworthiness cut, making for a more robust and accurate supernova classifier than before., Comment: 27 pages, 7 figures. The STag code is available at https://github.com/wdavison909/STag

We present the joint tomographic analysis of galaxy-galaxy lensing and galaxy clustering in harmonic space, using galaxy catalogues from the first three years of observations by the Dark Energy Survey (DES Y3). We utilise the redMaGiC and MagLim catalogues as lens galaxies and the METACALIBRATION catalogue as source galaxies. The measurements of angular power spectra are performed using the pseudo-$C_\ell$ method, and our theoretical modelling follows the fiducial analyses performed by DES Y3 in configuration space, accounting for galaxy bias, intrinsic alignments, magnification bias, shear magnification bias and photometric redshift uncertainties. We explore different approaches for scale cuts based on non-linear galaxy bias and baryonic effects contamination. Our fiducial covariance matrix is computed analytically, accounting for mask geometry in the Gaussian term, and including non-Gaussian contributions and super-sample covariance terms. To validate our harmonic space pipelines and covariance matrix, we used a suite of 1800 log-normal simulations. We also perform a series of stress tests to gauge the robustness of our harmonic space analysis. In the $\Lambda$CDM model, the clustering amplitude $S_8 =\sigma_8(\Omega_m/0.3)^{0.5}$ is constrained to $S_8 = 0.704\pm 0.029$ and $S_8 = 0.753\pm 0.024$ ($68\%$ C.L.) for the redMaGiC and MagLim catalogues, respectively. For the $w$CDM, the dark energy equation of state is constrained to $w = -1.28 \pm 0.29$ and $w = -1.26^{+0.34}_{-0.27}$, for redMaGiC and MagLim catalogues, respectively. These results are compatible with the corresponding DES Y3 results in configuration space and pave the way for harmonic space analyses using the DES Y6 data., Comment: To be submitted to MNRAS

The residuals of the distance moduli of Type Ia supernovae (SN Ia) relative to a Hubble diagram fit contain information about the inhomogeneity of the universe, due to weak lensing magnification by foreground matter. By correlating the residuals of the Dark Energy Survey Year 5 SN Ia sample (DES-SN5YR) with extra-galactic foregrounds from the DES Y3 Gold catalog, we detect the presence of lensing at $6.0 \sigma$ significance. This is the first detection with a significance level above $5\sigma$. Constraints on the effective mass-to-light ratios and radial profiles of dark-matter haloes surrounding individual galaxies are also obtained. We show that the scatter of SNe Ia around the Hubble diagram is reduced by modifying the standardisation of the distance moduli to include an easily calculable de-lensing (i.e., environmental) term. We use the de-lensed distance moduli to recompute cosmological parameters derived from SN Ia, finding in Flat $w$CDM a difference of $\Delta \Omega_{\rm M} = +0.036$ and $\Delta w = -0.056$ compared to the unmodified distance moduli, a change of $\sim 0.3\sigma$. We argue that our modelling of SN Ia lensing will lower systematics on future surveys with higher statistical power. We use the observed dispersion of lensing in DES-SN5YR to constrain $\sigma_8$, but caution that the fit is sensitive to uncertainties at small scales. Nevertheless, our detection of SN Ia lensing opens a new pathway to study matter inhomogeneity that complements galaxy-galaxy lensing surveys and has unrelated systematics., Comment: Submitted to MNRAS

We present $griz$ photometric light curves for the full 5 years of the Dark Energy Survey Supernova program (DES-SN), obtained with both forced Point Spread Function (PSF) photometry on Difference Images (DIFFIMG) performed during survey operations, and Scene Modelling Photometry (SMP) on search images processed after the survey. This release contains $31,636$ DIFFIMG and $19,706$ high-quality SMP light curves, the latter of which contains $1635$ photometrically-classified supernovae that pass cosmology quality cuts. This sample spans the largest redshift ($z$) range ever covered by a single SN survey ($0.1

Type Ia Supernovae (SNe Ia) are a critical tool in measuring the accelerating expansion of the universe. Recent efforts to improve these standard candles have focused on incorporating the effects of dust on distance measurements with SNe Ia. In this paper, we use the state-of-the-art Dark Energy Survey 5 year sample to evaluate two different families of dust models: empirical extinction models derived from SNe Ia data, and physical attenuation models from the spectra of galaxies. Among the SNe Ia-derived models, we find that a logistic function of the total-to-selective extinction RV best recreates the correlations between supernova distance measurements and host galaxy properties, though an additional 0.02 magnitudes of grey scatter are needed to fully explain the scatter in SNIa brightness in all cases. These empirically-derived extinction distributions are highly incompatible with the physical attenuation models from galactic spectral measurements. From these results, we conclude that SNe Ia must either preferentially select extreme ends of galactic dust distributions, or that the characterisation of dust along the SNe Ia line-of-sight is incompatible with that of galactic dust distributions.

We present a precise measurement of cosmological time dilation using the light curves of 1504 type Ia supernovae from the Dark Energy Survey spanning a redshift range $0.1\lesssim z\lesssim 1.2$. We find that the width of supernova light curves is proportional to $(1+z)$, as expected for time dilation due to the expansion of the Universe. Assuming type Ia supernovae light curves are emitted with a consistent duration $\Delta t_{\rm em}$, and parameterising the observed duration as $\Delta t_{\rm obs}=\Delta t_{\rm em}(1+z)^b$, we fit for the form of time dilation using two methods. Firstly, we find that a power of $b \approx 1$ minimises the flux scatter in stacked subsamples of light curves across different redshifts. Secondly, we fit each target supernova to a stacked light curve (stacking all supernovae with observed bandpasses matching that of the target light curve) and find $b=1.003\pm0.005$ (stat) $\pm\,0.010$ (sys). Thanks to the large number of supernovae and large redshift-range of the sample, this analysis gives the most precise measurement of cosmological time dilation to date, ruling out any non-time-dilating cosmological models at very high significance., Comment: 14 pages, 13 figures. Updated in response to reviewer feedback. Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

We measure the current expansion rate of the Universe, Hubble's constant $H_0$, by calibrating the absolute magnitudes of supernovae to distances measured by Baryon Acoustic Oscillations. This `inverse distance ladder' technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning $0.01\lt z \lt1.13$ anchored to the recent Baryon Acoustic Oscillation measurements from DESI spanning $0.30 \lt z_{\mathrm{eff}} \lt 2.33$. To trace cosmology to $z=0$, we use the third-, fourth- and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher-redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find $H_0=67.19^{+0.66}_{-0.64}\mathrm{~km} \mathrm{~s}^{-1} \mathrm{~Mpc}^{-1}$, in agreement with the value found by Planck without the need to assume Flat-$\Lambda$CDM. However the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

Reverberation mapping is the leading technique used to measure direct black hole masses outside of the local Universe. Additionally, reverberation measurements calibrate secondary mass-scaling relations used to estimate single-epoch virial black hole masses. The Australian Dark Energy Survey (OzDES) conducted one of the first multi-object reverberation mapping surveys, monitoring 735 AGN up to $z\sim4$, over 6 years. The limited temporal coverage of the OzDES data has hindered recovery of individual measurements for some classes of sources, particularly those with shorter reverberation lags or lags that fall within campaign season gaps. To alleviate this limitation, we perform a stacking analysis of the cross-correlation functions of sources with similar intrinsic properties to recover average composite reverberation lags. This analysis leads to the recovery of average lags in each redshift-luminosity bin across our sample. We present the average lags recovered for the H$\beta$, Mg II and C IV samples, as well as multi-line measurements for redshift bins where two lines are accessible. The stacking analysis is consistent with the Radius-Luminosity relations for each line. Our results for the H$\beta$ sample demonstrate that stacking has the potential to improve upon constraints on the $R-L$ relation, which have been derived only from individual source measurements until now., Comment: 20 pages, 15 figures. Accepted by MNRAS

We present the MOST Hosts survey (Multi-Object Spectroscopy of Transient Hosts). The survey is planned to run throughout the five years of operation of the Dark Energy Spectroscopic Instrument (DESI) and will generate a spectroscopic catalog of the hosts of most transients observed to date, in particular all the supernovae observed by most public, untargeted, wide-field, optical surveys (PTF/iPTF, SDSS II, ZTF, DECAT, DESIRT). Scientific questions for which the MOST Hosts survey will be useful include Type Ia supernova cosmology, fundamental plane and peculiar velocity measurements, and the understanding of the correlations between transients and their host galaxy properties. Here, we present the first release of the MOST Hosts survey: 21,931 hosts of 20,235 transients. These numbers represent 36% of the final MOST Hosts sample, consisting of 60,212 potential host galaxies of 38,603 transients (a transient can be assigned multiple potential hosts). Of these galaxies, 40% do not appear in the DESI primary target list and therefore require a specific program like MOST Hosts. Of all the transients in the MOST Hosts list, only 26.7% have existing classifications, and so the survey will provide redshifts (and luminosities) for nearly 30,000 transients. A preliminary Hubble diagram and a transient luminosity-duration diagram are shown as examples of future potential uses of the MOST Hosts survey. The survey will also provide a training sample of spectroscopically observed transients for photometry-only classifiers, as we enter an era when most newly observed transients will lack spectroscopic classification. The MOST Hosts DESI survey data will be released through the Wiserep platform on a rolling cadence and updated to match the DESI releases. Dates of future releases and updates are available through the https://mosthosts.desi.lbl.gov website., Comment: Submitted to ApJS

We present a new constraint on the Hubble constant ($H_0$) from the standard dark siren method using a sample of $5$ well-covered gravitational waves (GW) alerts reported during the first part of the fourth LIGO/Virgo/KAGRA observing run and with $3$ updated standard dark sirens from third observation run in combination the previous constraints from the first three runs. Our methodology relies on the galaxy catalogue method alone. We use a Deep Learning method to derive the full probability density estimation of photometric redshifts using the Legacy Survey catalogues. We add the constraints from well localized Binary Black Hole mergers to the sample of standard dark sirens analysed in our previous work. We combine the $H_0$ posterior for $5$ new standard sirens with other $10$ previous events (using the most recent available data for the 5 novel events and updated 3 previous posteriors from O3), finding $H_0 = 70.4^{+13.6}_{-11.7}~{\rm km~s^{-1}~Mpc^{-1}}$ (68\% Confidence interval) with the catalogue method only. This result represents an improvement of $\sim 23\%$ comparing the new $15$ dark siren constrain with the previous $10$ dark siren constraint and a reduction in uncertainty of $\sim 40\%$ from the combination of $15$ dark and bright sirens compared with the GW170817 bright siren alone. The combination of dark and bright siren GW170817 with recent jet constraints yields $H_0$ of $68.0^{+4.4}_{-3.8}~{\rm km~s^{-1}~Mpc^{-1}}$, a $\sim 6\%$ precision from Standard Sirens, reducing the previous constraint uncertainty by $\sim 10\%$ ., Comment: 18 pages, 11 figures and 2 tables. Article accepted for publication in MNRAS

Recent work by Moroianu et al. (2022) has suggested that the binary neutron star (BNS) merger GW190425 might have a potential fast radio burst (FRB) counterpart association, FRB 20190425A, at the 2.8$\sigma$ level of confidence with a likely host galaxy association, namely UGC10667. The authors argue that the observations are consistent with a long-lived hypermassive neutron star (HMNS) that formed promptly after the BNS merger and was stable for approximately 2.5 hours before promptly collapsing into a black hole. Recently, Bhardwaj et al. (2023) conclusively associated FRB 20190425A with UGC10667, potentially providing a direct host galaxy candidate for GW190425. In this work, we examine the multi-messenger association based on the space-time localization overlaps between GW190425 and the FRB host galaxy UGC10667 and find that the odds for a coincident association are $\mathcal{O}(5)$. We validate this estimate by using a Gaussian Process (GP) density estimator. Assuming that the association is indeed real, we then perform Bayesian parameter estimation on GW190425 assuming that the BNS event took place in UGC10667. We find that the viewing angle of GW190425 excludes an on-axis system at $p(\theta_v>30^o)\approx99.99$\%, highly favouring an off-axis system similar to GRB 170817A. We also find a slightly higher source frame total mass for the binary, namely, $m_{\rm{total}} = 3.42^{+0.34}_{-0.11} M_{\odot}$, leading to an increase on the probability of prompt collapse into a black hole and therefore disfavors the long-lived HMNS formation scenario. Given our findings, we conclude that the association between GW190425 and FRB 20190425A is disfavoured by current state-of-the-art gravitational-wave analyses., Comment: Comments are welcomed

We measure the stacked lensing signal in the direction of galaxy clusters in the Dark Energy Survey Year 3 (DES Y3) redMaPPer sample, using cosmic microwave background (CMB) temperature data from SPT-3G, the third-generation CMB camera on the South Pole Telescope (SPT). We estimate the lensing signal using temperature maps constructed from the initial 2 years of data from the SPT-3G 'Main' survey, covering 1500 deg$^2$ of the Southern sky. We then use this signal as a proxy for the mean cluster mass of the DES sample. In this work, we employ three versions of the redMaPPer catalogue: a Flux-Limited sample containing 8865 clusters, a Volume-Limited sample with 5391 clusters, and a Volume&Redshift-Limited sample with 4450 clusters. For the three samples, we find the mean cluster masses to be ${M}_{200{\rm{m}}}=1.66\pm0.13$ [stat.]$\pm0.03$ [sys.], $1.97\pm0.18$ [stat.]$\pm0.05$ [sys.], and $2.11\pm0.20$ [stat.]$\pm0.05$ [sys.]$\times{10}^{14}\ {\rm{M}}_{\odot }$, respectively. This is a factor of $\sim2$ improvement relative to the precision of measurements with previous generations of SPT surveys and the most constraining cluster mass measurements using CMB cluster lensing to date. Overall, we find no significant tensions between our results and masses given by redMaPPer mass-richness scaling relations of previous works, which were calibrated using CMB cluster lensing, optical weak lensing, and velocity dispersion measurements from various combinations of DES, SDSS and Planck data. We then divide our sample into 3 redshift and 3 richness bins, finding no significant tensions with optical weak-lensing calibrated masses in these bins. We forecast a $5.7\%$ constraint on the mean cluster mass of the DES Y3 sample with the complete SPT-3G surveys when using both temperature and polarization data and including an additional $\sim1400$ deg$^2$ of observations from the 'Extended' SPT-3G survey., Comment: 23 pages, 9 figures, accepted for publication in JCAP. Minor changes and corrections have been made relative to v1

Current and future Type Ia Supernova (SN Ia) surveys will need to adopt new approaches to classifying SNe and obtaining their redshifts without spectra if they wish to reach their full potential. We present here a novel approach that uses only photometry to identify SNe Ia in the 5-year Dark Energy Survey (DES) dataset using the SuperNNova classifier. Our approach, which does not rely on any information from the SN host-galaxy, recovers SNe Ia that might otherwise be lost due to a lack of an identifiable host. We select 2,298 high-quality SNe Ia from the DES 5-year dataset an almost complete sample of detected SNe Ia. More than 700 of these have no spectroscopic host redshift and are potentially new SNIa compared to the DES-SN5YR cosmology analysis. To analyse these SNe Ia, we derive their redshifts and properties using only their light-curves with a modified version of the SALT2 light-curve fitter. Compared to other DES SN Ia samples with spectroscopic redshifts, our new sample has in average higher redshift, bluer and broader light-curves, and fainter host-galaxies. Future surveys such as LSST will also face an additional challenge, the scarcity of spectroscopic resources for follow-up. When applying our novel method to DES data, we reduce the need for follow-up by a factor of four and three for host-galaxy and live SN respectively compared to earlier approaches. Our novel method thus leads to better optimisation of spectroscopic resources for follow-up., Comment: Accepted MNRAS

High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a timescale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these timescales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. We find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be $R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0} (\mathrm{systematic})\times10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$. The mass-normalised rate is $R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times10^{-17}~\mathrm{M_\odot^{-1}}~\mathrm{yr}^{-1}$ and the volumetric rate is $R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10 to 40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included., Comment: 29 pages, 14 figures. Accepted by MNRAS

In this paper we present and validate the galaxy sample used for the analysis of the baryon acoustic oscillation (BAO) signal in the Dark Energy Survey (DES) Y6 data. The definition is based on a color and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.6, while ensuring a high-quality photo-$z$ determination. The optimization is performed using a Fisher forecast algorithm, finding the optimal $i$-magnitude cut to be given by $i$<19.64+2.894$z_{\rm ph}$. For the optimal sample, we forecast an increase in precision in the BAO measurement of $\sim$25% with respect to the Y3 analysis. Our BAO sample has a total of 15,937,556 galaxies in the redshift range 0.6<$z_{\rm ph}$<1.2, and its angular mask covers 4,273.42 deg${}^2$ to a depth of $i$=22.5. We validate its redshift distributions with three different methods: directional neighborhood fitting algorithm (DNF), which is our primary photo-$z$ estimation; direct calibration with spectroscopic redshifts from VIPERS; and clustering redshift using SDSS galaxies. The fiducial redshift distribution is a combination of these three techniques performed by modifying the mean and width of the DNF distributions to match those of VIPERS and clustering redshift. In this paper we also describe the methodology used to mitigate the effect of observational systematics, which is analogous to the one used in the Y3 analysis. This paper is one of the two dedicated to the analysis of the BAO signal in DES Y6. In its companion paper, we present the angular diameter distance constraints obtained through the fitting to the BAO scale., Comment: 23 pages, 10 figures. Submitted to PRD

We present the angular diameter distance measurement obtained with the Baryonic Acoustic Oscillation feature from galaxy clustering in the completed Dark Energy Survey, consisting of six years (Y6) of observations. We use the Y6 BAO galaxy sample, optimized for BAO science in the redshift range 0.6<$z$<1.2, with an effective redshift at $z_{\rm eff}$=0.85 and split into six tomographic bins. The sample has nearly 16 million galaxies over 4,273 square degrees. Our consensus measurement constrains the ratio of the angular distance to sound horizon scale to $D_M(z_{\rm eff})/r_d$ = 19.51$\pm$0.41 (at 68.3% confidence interval), resulting from comparing the BAO position in our data to that predicted by Planck $\Lambda$CDM via the BAO shift parameter $\alpha=(D_M/r_d)/(D_M/r_d)_{\rm Planck}$. To achieve this, the BAO shift is measured with three different methods, Angular Correlation Function (ACF), Angular Power Spectrum (APS), and Projected Correlation Function (PCF) obtaining $\alpha=$ 0.952$\pm$0.023, 0.962$\pm$0.022, and 0.955$\pm$0.020, respectively, which we combine to $\alpha=$ 0.957$\pm$0.020, including systematic errors. When compared with the $\Lambda$CDM model that best fits Planck data, this measurement is found to be 4.3% and 2.1$\sigma$ below the angular BAO scale predicted. To date, it represents the most precise angular BAO measurement at $z$>0.75 from any survey and the most precise measurement at any redshift from photometric surveys. The analysis was performed blinded to the BAO position and it is shown to be robust against analysis choices, data removal, redshift calibrations and observational systematics., Comment: Submitted to PRD, 39 pages, 12 figures

Cosmological analyses using galaxy clusters in optical/NIR photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for {\it all} galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89\% membership accuracy with a mild dependency on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass--cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities ($\mu_{\star}$) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation. These results indicate the potential of Copacabana and $\mu_{\star}$ to be used in cosmological analyses of optically selected clusters in the future., Comment: 14 pages, 10 figures, submitted to MNRAS

We present the full Hubble diagram of photometrically-classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey supernova program (DES-SN). DES-SN discovered more than 20,000 SN candidates and obtained spectroscopic redshifts of 7,000 host galaxies. Based on the light-curve quality, we select 1635 photometrically-identified SNe Ia with spectroscopic redshift 0.10$< z <$1.13, which is the largest sample of supernovae from any single survey and increases the number of known $z>0.5$ supernovae by a factor of five. In a companion paper, we present cosmological results of the DES-SN sample combined with 194 spectroscopically-classified SNe Ia at low redshift as an anchor for cosmological fits. Here we present extensive modeling of this combined sample and validate the entire analysis pipeline used to derive distances. We show that the statistical and systematic uncertainties on cosmological parameters are $\sigma_{\Omega_M,{\rm stat+sys}}^{\Lambda{\rm CDM}}=$0.017 in a flat $\Lambda$CDM model, and $(\sigma_{\Omega_M},\sigma_w)_{\rm stat+sys}^{w{\rm CDM}}=$(0.082, 0.152) in a flat $w$CDM model. Combining the DES SN data with the highly complementary CMB measurements by Planck Collaboration (2020) reduces uncertainties on cosmological parameters by a factor of 4. In all cases, statistical uncertainties dominate over systematics. We show that uncertainties due to photometric classification make up less than 10% of the total systematic uncertainty budget. This result sets the stage for the next generation of SN cosmology surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time., Comment: 39 pages, 19 figures; Submitted to ApJ; companion paper Dark Energy Collaboration et al. on consecutive arxiv number 2401.02929

We present cosmological constraints from the sample of Type Ia supernovae (SN Ia) discovered during the full five years of the Dark Energy Survey (DES) Supernova Program. In contrast to most previous cosmological samples, in which SN are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being a SN Ia, we find 1635 DES SNe in the redshift range $0.100.5$ SNe compared to the previous leading compilation of Pantheon+, and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints we combine the DES supernova data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning $0.025

Context. The determination of accurate photometric redshifts (photo-zs) in large imaging galaxy surveys is key for cosmological studies. One of the most common approaches are machine learning techniques. These methods require a spectroscopic or reference sample to train the algorithms. Attention has to be paid to the quality and properties of these samples since they are key factors in the estimation of reliable photo-zs. Aims. The goal of this work is to calculate the photo-zs for the Y3 DES Deep Fields catalogue using the DNF machine learning algorithm. Moreover, we want to develop techniques to assess the incompleteness of the training sample and metrics to study how incompleteness affects the quality of photometric redshifts. Finally, we are interested in comparing the performance obtained with respect to the EAzY template fitting approach on Y3 DES Deep Fields catalogue. Methods. We have emulated -- at brighter magnitude -- the training incompleteness with a spectroscopic sample whose redshifts are known to have a measurable view of the problem. We have used a principal component analysis to graphically assess incompleteness and to relate it with the performance parameters provided by DNF. Finally, we have applied the results about the incompleteness to the photo-z computation on Y3 DES Deep Fields with DNF and estimated its performance. Results. The photo-zs for the galaxies on DES Deep Fields have been computed with the DNF algorithm and added to the Y3 DES Deep Fields catalogue. They are available at https://des.ncsa.illinois.edu/releases/y3a2/Y3deepfields. Some techniques have been developed to evaluate the performance in the absence of "true" redshift and to assess completeness. We have studied... (Partial abstract), Comment: 14 pages, 17 figures

Low Surface Brightness Galaxies (LSBGs) are excellent probes of quenching and other environmental processes near massive galaxies. We study an extensive sample of LSBGs near massive hosts in the local universe that are distributed across a diverse range of environments. The LSBGs with surface-brightness $\mu_{\rm eff,g}> $24.2 mag arcsec$^{-2}$ are drawn from the Dark Energy Survey Year 3 catalog while the hosts with masses $9.0< log(M_{\star}/M_{\odot})< 11.0$ comparable to the Milky Way and the Large Magellanic Cloud are selected from the z0MGS sample. We study the projected radial density profiles of LSBGs as a function of their color and surface brightness around hosts in both the rich Fornax-Eridanus cluster environment and the low-density field. We detect an overdensity with respect to the background density, out to 2.5 times the virial radius for both hosts in the cluster environment and the isolated field galaxies. When the LSBG sample is split by $g-i$ color or surface brightness $\mu_{\rm eff,g}$, we find the LSBGs closer to their hosts are significantly redder and brighter, like their high surface-brightness counterparts. The LSBGs form a clear 'red sequence' in both the cluster and isolated environments that is visible beyond the virial radius of the hosts. This suggests a pre-processing of infalling LSBGs and a quenched backsplash population around both host samples. However, the relative prominence of the 'blue cloud' feature implies that pre-processing is ongoing near the isolated hosts compared to the cluster hosts., Comment: 22 pages, 7 figures

We present a measurement of the Hubble Constant $H_0$ using the gravitational wave event GW190412, an asymmetric binary black hole merger detected by LIGO/Virgo, as a dark standard siren. This event does not have an electromagnetic counterpart, so we use the statistical standard siren method and marginalize over potential host galaxies from the Dark Energy Spectroscopic Instrument (DESI) survey. GW190412 is well-localized to 12 deg$^2$ (90% credible interval), so it is promising for a dark siren analysis. The dark siren value for $H_0=85.4_{-33.9}^{+29.1}$ km/s/Mpc, with a posterior shape that is consistent with redshift overdensities. When combined with the bright standard siren measurement from GW170817 we recover $H_0=77.96_{-5.03}^{+23.0}$ km/s/Mpc, consistent with both early and late-time Universe measurements of $H_0$. This work represents the first standard siren analysis performed with DESI data, and includes the most complete spectroscopic sample used in a dark siren analysis to date., Comment: Submitted to RNAAS

FRB 20220610A is a high-redshift Fast Radio Burst (FRB) that has not been observed to repeat. Here, we present rest-frame UV and optical $\textit{Hubble Space Telescope}$ observations of the field of FRB 20220610A. The imaging reveals seven extended sources, one of which we identify as the most likely host galaxy with a spectroscopic redshift of $z$=1.017. We spectroscopically confirm at least three additional sources to be at the same redshift, and identify the system as a compact galaxy group with possible signs of interaction among group members. We determine the host of FRB 20220610A to be a star-forming galaxy with stellar mass of $\approx10^{9.7}\,M_{\odot}$, mass-weighted age of $\approx2.6$~Gyr, and star formation rate (integrated over the last 100 Myr) of $\approx1.7$~M$_{\odot}$~yr$^{-1}$. These host properties are commensurate with the star-forming field galaxy population at z~1 and trace their properties analogously to the population of low-$z$ FRB hosts. Based on estimates of the total stellar mass of the galaxy group, we calculate a fiducial contribution to the observed Dispersion Measure (DM) from the intragroup medium of $\approx 110-220$ $\rm pc \, cm^{-3}$ (rest-frame). This leaves a significant excess of $500^{+272}_{-109}$ $\rm pc \, cm^{-3}$ (in the observer frame), with additional sources of DM possibly originating from the circumburst environment, host galaxy interstellar medium, and/or foreground structures along the line of sight. Given the low occurrence rates of galaxies in compact groups, the discovery of an FRB in such a group demonstrates a rare and novel environment in which FRBs can occur., Comment: 24 pages, 8 figures, 2 tables, submitted

The current and next observation seasons will detect hundreds of gravitational waves (GWs) from compact binary systems coalescence at cosmological distances. When combined with independent electromagnetic measurements, the source redshift will be known, and we will be able to obtain precise measurements of the Hubble constant $H_0$ via the distance-redshift relation. However, most observed mergers are not expected to have electromagnetic counterparts, which prevents a direct redshift measurement. In this scenario, one of the possibilities is to use the dark sirens method that statistically marginalizes over all the potential host galaxies within the GW location volume to provide a probabilistic redshift to the source. Here we presented $H_{0}$ measurements using two new dark sirens compared to previous analyses using DECam data, GW190924$\_$021846 and GW200202$\_$154313. The photometric redshifts of the possible host galaxies of these two events are acquired from the DECam Local Volume Exploration Survey (DELVE) carried out on the Blanco telescope at Cerro Tololo in Chile. The combination of the $H_0$ posterior from GW190924$\_$021846 and GW200202$\_$154313 together with the bright siren GW170817 leads to $H_{0} = 68.84^{+15.51}_{-7.74}\, \rm{km/s/Mpc}$. Including these two dark sirens improves the 68% confidence interval (CI) by 7% over GW170817 alone. This demonstrates that the inclusion of well-localized dark sirens in such analysis improves the precision with which cosmological measurements can be made. Using a sample containing 10 well-localized dark sirens observed during the third LIGO/Virgo observation run, we determine a measurement of $H_{0} = 76.00^{+17.64}_{-13.45}\, \rm{km /s/Mpc}$., Comment: v2: minor corrections and references added

We use Dark Energy Survey Year 3 (DES Y3) clusters with archival X-ray data from XMM-Newton and Chandra to assess the centering performance of the redMaPPer cluster finder and to measure key richness observable scaling relations. In terms of centering, we find that 10-20% of redMaPPer clusters are miscentered with no significant difference in bins of low versus high richness ($20<\lambda<40$ and $\lambda>40$) or redshift ($0.2

We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibration using DES data. For the DES Year 3 data, we report a systematic uncertainty in weak-lensing mass calibration that increases from 1% at $z=0.25$ to 10% at $z=0.95$, to which we add 2% in quadrature to account for uncertainties in the impact of baryonic effects. We implement an analysis pipeline that joins the cluster abundance likelihood with a multi-observable likelihood for the Sunyaev-Zel'dovich effect, optical richness, and weak-lensing measurements for each individual cluster. We validate that our analysis pipeline can recover unbiased cosmological constraints by analyzing mocks that closely resemble the cluster sample extracted from the SPT-SZ, SPTpol ECS, and SPTpol 500d surveys and the DES Year 3 and HST-39 weak-lensing datasets. This work represents a crucial prerequisite for the subsequent cosmological analysis of the real dataset., Comment: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around $10^5$ clusters with mass and redshift ranges $10^{13.7} < M_{\rm 200m}/M_\odot < 10^{15.5}$ and $0.1 < z < 2$, and the total sky coverage of the maps is $\approx 15,000 \,\,{\rm deg}^2$. We find a clear pressure deficit at $R/R_{\rm 200m}\approx 1.1$ in SZ profiles around both ACT and SPT clusters, estimated at $6\sigma$ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2-5 compared to predictions, and we discuss possible interpretations of this behavior. An oriented stacking of clusters -- where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogs -- shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius., Comment: [v2]: Version accepted to MNRAS

We introduce a software suite developed for galaxy cluster cosmological analysis with the Dark Energy Survey Data. Cosmological analyses based on galaxy cluster number counts and weak-lensing measurements need efficient software infrastructure to explore an increasingly large parameter space, and account for various cosmological and astrophysical effects. Our software package is designed to model the cluster observables in a wide-field optical survey, including galaxy cluster counts, their averaged weak-lensing masses, or the cluster's averaged weak-lensing radial signals. To ensure maximum efficiency, this software package is developed in C++ in the CosmoSIS software framework, making use of the CUBA integration library. We also implement a testing and validation scheme to ensure the quality of the package. We demonstrate the effectiveness of this development by applying the software to the Dark Energy Survey Year 1 galaxy cluster cosmological data sets, and acquired cosmological constraints that are consistent with the fiducial Dark Energy Survey analysis.

Using the full six years of imaging data from the Dark Energy Survey, we study the surface brightness profiles of galaxy cluster central galaxies and intra-cluster light. We apply a ``stacking'' method to over four thousand galaxy clusters identified by the redMaPPer cluster finding algorithm in the redshift range of 0.2 to 0.5. This yields high signal-to-noise radial profile measurements of the central galaxy and intra-cluster light out to 1 Mpc from the cluster center. Using redMaPPer richness as a cluster mass indicator, we find that the intra-cluster light brightness has a strong mass dependence throughout the 0.2 to 0.5 redshift range, and the dependence grows stronger at a larger radius. In terms of redshift evolution, we find some evidence that the central galaxy, as well as the diffuse light within the transition region between the cluster central galaxy and intra-cluster light within 80 kpc from the center, may be growing over time. At larger radii, more than 80 kpc away from the cluster center, we do not find evidence of additional redshift evolution beyond the cluster mass dependence, which is consistent with the findings from the IllustrisTNG hydrodynamic simulation. We speculate that the major driver of intra-cluster light growth, especially at large radii, is associated with cluster mass growth. Finally, we find that the color of the cluster central galaxy and intra-cluster light displays a radial gradient that becomes bluer at a larger radius, which is consistent with a stellar stripping and disruption origin of intra-cluster light as suggested by simulation studies., Comment: Submitted to MNRAS

Fitting model spectral energy distributions (SED) to galaxy photometric data is a widely used method to recover galaxy parameters from galaxy surveys. However, the parameter space used to describe galaxies is wide and interdependent, and distinctions between real and spurious correlations that are found between these parameters can be difficult to discern. In this work, we use the SED fitting code BAGPIPES to investigate degeneracies between galaxy parameters and the effect of the choice of different sets of photometric bands. In particular, we focus on optical to infrared wavelength coverage, and on two parameters describing the galaxies' dust attenuation law: $A_V$ and $\delta$, which characterize dust column density and the slope of a flexible dust attenuation law, respectively. We demonstrate that 1) a degeneracy between the residual (the difference between truth and recovered value) $A_V$ and star formation rate exists, but this is lifted when WISE bands are included; 2) BAGPIPES is able to accurately recover the input $A_V$ and $\delta$ distributions and relations (differences in slope of less than 1.7$\sigma$ for a flat relation, less than 1.2$\sigma$ for an observationally-motivated relation from Salim et al. 2018) and is not introducing spurious correlations between these parameters. Our findings suggest that the information needed to constrain $A_V$ and $\delta$ well enough individually exists in the data, especially when IR is added. This indicates that recent works finding a correlation between $A_V$ and $\delta$ are not being misled by fitting degeneracies from their SED fitting code., Comment: 15 pages

Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at multiple scales as a summary of the galaxy lensing convergence field. Using a suite of N-body lightcone simulations, we show the CDFs' constraining power is modestly better than that of the 2nd and 3rd moments of the field, as they approximately capture the information from all moments of the field in a concise data vector. We then study the practical aspects of applying the CDFs to observational data, using the first three years of the Dark Energy Survey (DES Y3) data as an example, and compute the impact of different systematics on the CDFs. The contributions from the point spread function are 2-3 orders of magnitude below the cosmological signal, while those from reduced shear approximation contribute $\lesssim 1\%$ to the signal. Source clustering effects and baryon imprints contribute $1-10\%$. Enforcing scale cuts to limit systematics-driven biases in parameter constraints degrades these constraints a noticeable amount, and this degradation is similar for the CDFs and the moments. We also detect correlations between the observed convergence field and the shape noise field at $13\sigma$. We find that the non-Gaussian correlations in the noise field must be modeled accurately to use the CDFs, or other statistics sensitive to all moments, as a rigorous cosmology tool., Comment: 21 pages, 12 figures