Marcelle Soares Santos

We present constraints on ultralight bosons using binary black hole mergers observed in the second and third parts of the fourth LIGO-Virgo-KAGRA observing run. Directed searches are conducted for long-transient gravitational waves from ultralight vector boson clouds around merger remnants, using a hidden-Markov-model (HMM) tracking scheme. We target the remnant black holes formed in the binary coalescences that produced GW250114 and GW250207. We find no evidence for such signals from either target. Estimating our search sensitivity at a threshold corresponding to a 1% false alarm probability, we thus disfavor vector boson masses in the range of [ 2.80 , 3.95 ] × 10 − 13 eV with greater than 90% confidence. In addition, we derive constraints on ultralight scalar and vector bosons from the inferred high spins of the constituent black holes in three binaries, using events GW240515, GW241113, and GW241225_08. The excluded mass ranges in this approach depend on the assumed black-hole ages. At 10 5 years , corresponding to typical dynamically formed binaries, we exclude scalar and vector bosons in the ranges [ 1.39 , 6.94 ] × 10 − 13 eV and [ 0.32 , 14.4 ] × 10 − 13 eV at 90% confidence, respectively.

The recent candidate gravitational-wave (GW) alert from a compact object merger involving at least one subsolar mass (SSM) object has prompted questions about their origins. S251112cm is reported by LIGO/Virgo with a false-alarm rate of 1 per 6.2 yr, nearby luminosity distance of 93 ± 27 Mpc, and probability of containing an SSM object of 100%. Such a system, if astrophysical, likely did not involve the supersolar neutron stars or black holes invoked to explain kilonovae. One must then also invoke hitherto unobserved and speculative models to produce SSM mergers that may have electromagnetic (EM) counterparts. We introduce a framework that vets and scores candidate counterparts to SSM GW events to inform follow-up in search of any among the zoo of potential EM transients: kilonovae, kilonovae-within-supernovae, superkilonovae, or active galactic nucleus (AGN) flares from binary black hole mergers. We use a suite of telescopes to perform tiling, galaxy-targeted observations, and photometric/spectroscopic follow-up of promising candidates. In near-real time, we ingest candidates reported by the community, including some of the first observations reported by the Vera C. Rubin Observatory. We vet and score a total of 456 candidates, including 67 from Rubin, but find no likely counterpart. We nonetheless highlight candidates that demonstrate the ability of our framework to distinguish between different transient types and describe strategies to maximize the chances of detecting a counterpart to the next SSM event. Our framework will be implemented in the forthcoming Multimessenger Tool for Rapid Object Vetting and Examination.

We present a new measurement of the Hubble constant ( H ₀ ) resulting from the first joint analysis of standard sirens with weak gravitational lensing and galaxy clustering observables comprising three two-point correlation functions (3 × 2 pt). For the 3 × 2 pt component of the analysis, we used data from the Dark Energy Survey (DES) Year 3 release. For the standard siren component, we used data from the Gravitational-Wave Transient Catalog 4.0 released by the LIGO-Virgo-KAGRA (LVK) Collaboration. For GW170817, the only standard siren for which extensive electromagnetic follow-up observations exist, we also used measurements of the host galaxy redshift and inclination angle estimates derived from observations of a superluminal jet from its remnant. Assuming a flat Λ cold dark matter model, our joint analysis yields H ₀ = 67.94 ^+4.40 _−4.34 km s ⁻¹ Mpc ⁻¹ , a 6.4% measurement, while improving the DES constraint on the total abundance of matter Ω _m by 22%. When the jet information is removed, the H ₀ precision decreases to 9.9%. The measurement of H ₀ remains a central problem in cosmology, with a multitude of approaches being vigorously pursued with the aim to reconcile significantly discrepant measurements at the percent level. In light of the impending new data releases from DES and LVK, and anticipating a much higher constraining power from 3 × 2 pt observables using newly commissioned survey instruments, we demonstrate that incorporating standard sirens in the cosmology framework of large cosmic surveys is a viable route towards that goal.

We describe the photometric dataset assembled from the full 6 yr of observations by the Dark Energy Survey (DES) in support of static-sky cosmology analyses. DES Y6 Gold is a curated dataset derived from DES Data Release 2 (DR2) that incorporates improved measurement, photometric calibration, object classification and value-added information. Y6 Gold comprises nearly 5000 deg ² of grizY imaging in the south Galactic cap and includes 669 million objects with a depth of i _AB ∼ 23.4 mag at a signal-to-noise ratio ∼ 10 for extended objects and a top-of-the-atmosphere photometric uniformity i _AB < 22.5. Additionally, it includes per-object quality information, and accompanying maps of the footprint coverage, masked regions, imaging depth, survey conditions, and astrophysical foregrounds that are used for cosmology analyses. After quality selections, benchmark samples contain 448 million galaxies and 120 million stars. This publication is complemented by data access and documentation.

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 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 (⁠|$\rm \mu _{\star }$|⁠) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation at low redshift for a LSST type survey. These results indicate the potential of Copacabana and |$\rm \mu _{\star }$| to be used in cosmological analyses of optically selected clusters in the future.

We present cosmological constraints from the sample of Type Ia supernovae (SNe Ia) discovered and measured during the full 5 yr of the Dark Energy Survey (DES) SN program. In contrast to most previous cosmological samples, in which SNe 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 an SN Ia, we find 1635 DES SNe in the redshift range 0.10 < z < 1.13 that pass quality selection criteria sufficient to constrain cosmological parameters. This quintuples the number of high-quality z > 0.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 SN data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning 0.025 < z < 0.10. Using SN data alone and including systematic uncertainties, we find Ω_M = 0.352 ± 0.017 in flat ΛCDM. SN data alone now require acceleration (q ₀ < 0 in ΛCDM) with over 5σ confidence. We find ( Ω M , w ) = ( 0.264 − 0.096 + 0.074 , − 0.80 − 0.16 + 0.14 ) in flat wCDM. For flat w ₀ w _a CDM, we find ( Ω M , w 0 , w a ) = ( 0.495 − 0.043 + 0.033 , − 0.36 − 0.30 + 0.36 , − 8.8 − 4.5 + 3.7 ) , consistent with a constant equation of state to within ∼2σ. Including Planck cosmic microwave background, Sloan Digital Sky Survey baryon acoustic oscillation, and DES 3 × 2pt data gives (Ω_M, w) = (0.321 ± 0.007, −0.941 ± 0.026). In all cases, dark energy is consistent with a cosmological constant to within ∼2σ. Systematic errors on cosmological parameters are subdominant compared to statistical errors; these results thus pave the way for future photometrically classified SN analyses.

Redshift measurements, primarily obtained from host galaxies, are essential for inferring cosmological parameters from type Ia supernovae (SNe Ia). Matching SNe to host galaxies using images is nontrivial, resulting in a subset of SNe with mismatched hosts and thus incorrect redshifts. We evaluate the host galaxy mismatch rate and resulting biases on cosmological parameters from simulations modeled after the Dark Energy Survey 5 Yr (DES-SN5YR) photometric sample. For both DES-SN5YR data and simulations, we employ the directional light radius method for host galaxy matching. In our SN Ia simulations, we find that 1.7% of SNe are matched to the wrong host galaxy, with redshift differences between the true and matched hosts of up to 0.6. Using our analysis pipeline, we determine the shift in the dark energy equation of state parameter (Δw) due to including SNe with incorrect host galaxy matches. For SN Ia–only simulations, we find Δw = 0.0013 ± 0.0026 with constraints from the cosmic microwave background. Including core-collapse SNe and peculiar SNe Ia in the simulation, we find that Δw ranges from 0.0009 to 0.0032, depending on the photometric classifier used. This bias is an order of magnitude smaller than the expected total uncertainty on w from the DES-SN5YR sample of ∼0.03. We conclude that the bias on w from host galaxy mismatch is much smaller than the uncertainties expected from the DES-SN5YR sample, but we encourage further studies to reduce this bias through better host-matching algorithms or selection cuts.

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