Search

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., Comment: 27 pages, 24 figures, accepted by The Astrophysical Journal

Future constraints of cosmological parameters from Type Ia supernovae (SNe Ia) will depend on the use of photometric samples, those samples without spectroscopic measurements of the SNe Ia. There is a growing number of analyses that show that photometric samples can be utilised for precision cosmological studies with minimal systematic uncertainties. To investigate this claim, we perform the first analysis that combines two separate photometric samples, SDSS and Pan-STARRS, without including a low-redshift anchor. We evaluate the consistency of the cosmological parameters from these two samples and find they are consistent with each other to under $1\sigma$. From the combined sample, named Amalgame, we measure $\Omega_M = 0.328 \pm 0.024$ with SN alone in a flat $\Lambda$CDM model, and $\Omega_M = 0.330 \pm 0.018$ and $w = -1.016^{+0.055}_{-0.058}$ when combining with a Planck data prior and a flat $w$CDM model. These results are consistent with constraints from the Pantheon+ analysis of only spectroscopically confirmed SNe Ia, and show that there are no significant impediments to analyses of purely photometric samples of SNe Ia., Comment: Submitting to MNRAS; comments welcome

Bayesian Estimation Applied to Multiple Species (BEAMS) is implemented in the BEAMS with Bias Corrections (BBC) framework to produce a redshift-binned Hubble diagram (HD) for Type Ia supernovae (SNe Ia). BBC corrects for selection effects and non-SNIa contamination, and systematic uncertainties are described by a covariance matrix with dimension matching the number of BBC redshift bins. For spectroscopically confirmed SN Ia samples, a recent "Binning is Sinning" article (BHS21, arxiv:2012.05900) showed that an unbinned HD and covariance matrix reduces the systematic uncertainty by a factor of ~1.5 compared to the binned approach. Here we extend their analysis to obtain an unbinned HD for a photometrically identified sample processed with BBC. To test this new method, we simulate and analyze 50 samples corresponding to the Dark Energy Survey (DES) witha low-redshift anchor; the simulation includes SNe Ia, and contaminants from core-collapse SNe and peculiar SNe Ia. The analysis includes systematic uncertainties for calibration, and measures the dark energy equation of state parameter (w). Compared to a redshift-binned HD, the unbinned HD with nearly 2000 events results in a smaller systematic uncertainty, in qualitative agreement with BHS21, and averaging results among the 50 samples we find no evidence for a w-bias. To reduce computation time for fitting an unbinned HD with large samples, we propose an HD-rebinning method that defines the HD in bins of redshift, color, and stretch; the rebinned HD results in similar uncertainty as the unbinned case, and shows no evidence for a w-bias.

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will discover an unprecedented number of supernovae (SNe), making spectroscopic classification for all the events infeasible. LSST will thus rely on photometric classification, whose accuracy depends on the not-yet-finalized LSST observing strategy. In this work, we analyze the impact of cadence choices on classification performance using simulated multi-band light curves. First, we simulate SNe with an LSST baseline cadence, a non-rolling cadence, and a presto-color cadence which observes each sky location three times per night instead of twice. Each simulated dataset includes a spectroscopically-confirmed training set, which we augment to be representative of the test set as part of the classification pipeline. Then, we use the photometric transient classification library snmachine to build classifiers. We find that the active region of the rolling cadence used in the baseline observing strategy yields a 25% improvement in classification performance relative to the background region. This improvement in performance in the actively-rolling region is also associated with an increase of up to a factor of 2.7 in the number of cosmologically-useful Type Ia supernovae relative to the background region. However, adding a third visit per night as implemented in presto-color degrades classification performance due to more irregularly sampled light curves. Overall, our results establish desiderata on the observing cadence related to classification of full SNe light curves, which in turn impacts photometric SNe cosmology with LSST., Comment: 22 pages, 14 figures. Changed to match version accepted by the Astrophysical Journal Supplement Series (accepted 06/02/2023)

We perform a rigorous cosmology analysis on simulated type Ia supernovae (SN~Ia) and evaluate the improvement from including photometric host-galaxy redshifts compared to using only the "zspec" subset with spectroscopic redshifts from the host or SN. We use the Deep Drilling Fields (~50 deg^2) from the Photometric LSST Astronomical Time-Series Classification Challenge (PLaSTiCC), in combination with a low-z sample based on Data Challenge2 (DC2). The analysis includes light curve fitting to standardize the SN brightness, a high-statistics simulation to obtain a bias-corrected Hubble diagram, a statistical+systematics covariance matrix including calibration and photo-z uncertainties, and cosmology fitting with a prior from the cosmic microwave background. Compared to using the zspec subset, including events with SN+host photo-z results in i) more precise distances for z>0.5, ii) a Hubble diagram that extends 0.3 further in redshift, and iii) a 50 % increase in the Dark Energy Task Force figure of merit (FoM) based on the w0-wa CDM model. Analyzing 25 simulated data samples, the average bias on w0 and wa is consistent with zero. The host photo-z systematic of 0.01 reduces FoM by only 2 % because i) most z<0.5 events are in the zspec subset, ii) the combined SN+host photo-z has X 2 smaller bias, and iii) the anti-correlation between fitted redshift and color self corrects distance errors. To prepare for analysing real data, the next SNIa-cosmology analysis with photo-z's should include non SN-Ia contamination and host galaxy mis-associations., Comment: 13 pages, 9 figures, 7 tables. Submitted to ApJ

As we observe a rapidly growing number of astrophysical transients, we learn more about the diverse host galaxy environments in which they occur. Host galaxy information can be used to purify samples of cosmological Type Ia supernovae, uncover the progenitor systems of individual classes, and facilitate low-latency follow-up of rare and peculiar explosions. In this work, we develop a novel data-driven methodology to simulate the time-domain sky that includes detailed modeling of the probability density function for multiple transient classes conditioned on host galaxy magnitudes, colours, star formation rates, and masses. We have designed these simulations to optimize photometric classification and analysis in upcoming large synoptic surveys. We integrate host galaxy information into the SNANA simulation framework to construct the Simulated Catalogue of Optical Transients and Correlated Hosts (SCOTCH), a publicly-available catalogue of 5 million idealized transient light curves in LSST passbands and their host galaxy properties over the redshift range $0

We present results from simulating slitless spectroscopic observations with the Nancy Grace Roman Space Telescope's (Roman) Wide-Field Instrument (WFI) P127 prism spanning 0.75 $\mu m$ to 1.8 $\mu m$. We quantify the efficiency of recovered Type Ia supernovae (SNe Ia) redshifts, as a function of P127 prism exposure time, to guide planning for future observing programs with the Roman prism. Generating the two-dimensional dispersed images and extracting one-dimensional spectra is done with the slitless spectroscopy package pyLINEAR along with custom-written software. From the analysis of 1698 simulated SN Ia P127 prism spectra, we show the efficiency of recovering SN redshifts to $z\lesssim3.0$, highlighting the exceptional sensitivity of the Roman P127 prism. Redshift recovery is assessed by setting a requirement of $\sigma_z = (\left|z - z_\mathrm{true} \right|)/(1+z) \leq 0.01$. We find that 3 hr exposures are sufficient for meeting this requirement, for $\gtrsim 50\%$ of the sample of mock SNe Ia at $z\approx2$ and within $\pm5$ days of rest-frame maximum light in the optical. We also show that a 1 hr integration of Roman can achieve the same precision in completeness to a depth of $24.4 \pm 0.06$ AB mag (or $z\lesssim 1$). Implications for cosmological studies with Roman P127 prism spectra of SNe Ia are also discussed., Comment: ApJ accepted. Updated to match accepted version

Recent studies have shown that the observed colour distributions of Type Ia SNe (SNIa) are well-described by a combination of distributions from dust and intrinsic colour. Here we present a new forward-modeling fitting method (Dust2Dust) to measure the parent dust and colour distributions, including their dependence on host-galaxy mass. At each fit step, the SNIa selection efficiency is determined from a large simulated sample that is re-weighted to reflect the proposed distributions. We use five separate metrics to constrain the Dust2Dust parameters: distribution of fitted light-curve colour $c$, cosmological residual trends with $c$, cosmological residual scatter with $c$, fitted colour-luminosity relationship $\beta_{\rm SALT2}$, and intrinsic scatter $\sigma_{\rm int}$. Using the Pantheon+ data sample, we present results for a Dust2Dust fit that includes 4 parameters describing intrinsic colour variations and 8 parameters describing dust. Furthermore, we propagate the Dust2Dust parameter uncertainties and covariance to the dark energy equation-of-state $w$ and Hubble constant H$_0$: we find $\sigma_w = 0.005$ and $\sigma_{\textrm{H}_0} = 0.145~$km/s/Mpc. The Dust2Dust code is publically available.

We examine the redshifts of a comprehensive set of published Type Ia supernovae, and provide a combined, improved catalogue with updated redshifts. We improve on the original catalogues by using the most up-to-date heliocentric redshift data available; ensuring all redshifts have uncertainty estimates; using the exact formulae to convert heliocentric redshifts into the Cosmic Microwave Background (CMB) frame; and utilising an improved peculiar velocity model that calculates local motions in redshift-space and more realistically accounts for the external bulk flow at high-redshifts. In total we reviewed 2821 supernova redshifts; 534 are comprised of repeat-observations of the same supernovae and 1764 pass the cosmology sample quality cuts. We found 5 cases of missing or incorrect heliocentric corrections, 44 incorrect or missing supernova coordinates, 230 missing heliocentric or CMB frame redshifts, and 1200 missing redshift uncertainties. Of the 2287 unique Type Ia supernovae in our sample (1594 of which satisfy cosmology-sample cuts) we updated 990 heliocentric redshifts. The absolute corrections range between $10^{-8} \leq \Delta z \leq 0.038$, and RMS$(\Delta z) \sim 3\times 10^{-3}$. The sign of the correction was essentially random, so the mean and median corrections are small: $4\times 10^{-4}$ and $4\times 10^{-6}$ respectively. We examine the impact of these improvements for $H_0$ and the dark energy equation of state $w$ and find that the cosmological results change by $\Delta H_0 = -0.11$ km s$^{-1}$ Mpc$^{-1}$ and $\Delta w = -0.001$, both significantly smaller than previously reported uncertainties for $H_0$ of 1.4 km s$^{-1}$ Mpc$^{-1}$ and $w$ of 0.04 respectively., Comment: 24 pages, 13 figures, 7 tables. Updated to match version published in PASA. Data and papers available at https://pantheonplussh0es.github.io/