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/

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond., Comment: 32 pages, 18 figures, accepted for publication by ApJ

Improving the use of Type Ia supernovae (SNIa) as standard candles requires a better approach to incorporate the relationship between SNIa and the properties of their host galaxies. Using a spectroscopically-confirmed sample of $\sim$1600 SNIa, we develop the first empirical model of underlying populations for SNIa light-curve properties that includes their dependence on host-galaxy stellar mass. These populations are important inputs to simulations that are used to model selection effects and correct distance biases within the BEAMS with Bias Correction (BBC) framework. Here we improve BBC to also account for SNIa-host correlations, and we validate this technique on simulated data samples. We recover the input relationship between SNIa luminosity and host-galaxy stellar mass (the mass step, $\gamma$) to within 0.004 mags, which is a factor of 5 improvement over the previous method that results in a $\gamma$-bias of ${\sim}0.02$. We adapt BBC for a novel dust-based model of intrinsic brightness variations, which results in a greatly reduced mass step for data ($\gamma = 0.017 \pm 0.008$), and for simulations ($\gamma =0.006 \pm 0.007$). Analysing simulated SNIa, the biases on the dark energy equation-of-state, $w$, vary from $\Delta w = 0.006(5)$ to $0.010(5)$ with our new BBC method; these biases are significantly smaller than the $0.02(5)$ $w$-bias using previous BBC methods that ignore SNIa-host correlations.

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607$\pm$7 days) at 99.95% significance (with a global p-value of ~$10^{-3}$ accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67-002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ~20-yr time baseline spans ~4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302-102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna., Comment: Accepted to MNRAS, 17 pages, 9 figures

Improvement in the precision of measurements of cosmological parameters with Type Ia Supernovae (SNIa) is expected to come from large photometrically identified (photometric) SN samples. Here we re-analyse the SDSS photometric SN sample, with roughly 700 high-quality, likely but unconfirmed SNIa light-curves, to develop new analysis tools aimed at evaluating systematic uncertainties on the dark energy equation-of-state parameter $w$. Since we require a spectroscopically measured host galaxy redshift for each SN, we determine the associated selection efficiency of host galaxies in order to simulate bias corrections. We determine that the mis-association rate of host galaxies is $0.6\%$; ignoring this effect in simulated bias corrections leads to a $w$-bias of $\Delta w = +0.0007$, where $w$ is evaluated from SNIa and priors from measurements of baryon acoustic oscillations and the cosmic microwave background. We assess the uncertainty in our modeling of the host galaxy selection efficiency and find the associated $w$ uncertainty to be $-0.0072$. Finally, we explore new core collapse (CC) models in simulated training samples and find that adjusting the CC luminosity distribution to be in agreement with previous Pan-STARRS analyses yields a better match to the SDSS data. The impact of ignoring this adjustment is $\Delta w = -0.0109$; the impact of replacing the new CC models with those used by Pan-STARRS is $\Delta w = -0.0028$. These systematic uncertainties are subdominant to the statistical constraints from the SDSS sample, but must be considered in future photometric analyses of large SN samples such as those from DES, LSST and WFIRST.

The Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) is an open data challenge to classify simulated astronomical time-series data in preparation for observations from the Large Synoptic Survey Telescope (LSST), which will achieve first light in 2019 and commence its 10-year main survey in 2022. LSST will revolutionize our understanding of the changing sky, discovering and measuring millions of time-varying objects. In this challenge, we pose the question: how well can we classify objects in the sky that vary in brightness from simulated LSST time-series data, with all its challenges of non-representativity? In this note we explain the need for a data challenge to help classify such astronomical sources and describe the PLAsTiCC data set and Kaggle data challenge, noting that while the references are provided for context, they are not needed to participate in the challenge., Comment: Research note to accompany the https://www.kaggle.com/c/PLAsTiCC-2018 challenge

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible., Comment: 312 pages, 90 figures. Browse the current version at https://github.com/LSSTScienceCollaborations/ObservingStrategy, new contributions welcome!

In response to a recently reported observation of evidence for two classes of Type Ia Supernovae (SNe Ia) distinguished by their brightness in the rest-frame near ultraviolet (NUV), we search for the phenomenon in publicly available light-curve data. We use the SNANA supernova analysis package to simulate SN Ia-light curves in the Sloan Digital Sky Survey Supernova Search (SDSS) and the Supernova Legacy Survey (SNLS) with a model of two distinct ultraviolet classes of SNe Ia and a conventional model with a single broad distribution of SN-Ia ultraviolet brightnesses. We compare simulated distributions of rest-frame colors with these two models to those observed in 158 SNe Ia in the SDSS and SNLS data. The SNLS sample of 99 SNe Ia is in clearly better agreement with a model with one class of SN Ia light curves and shows no evidence for distinct NUV sub-classes. The SDSS sample of 59 SNe Ia with poorer color resolution does not distinguish between the two models., Comment: 8 pages, in review at MNRAS

We present a new technique to create a bin-averaged Hubble Diagram (HD) from photometrically identified SN~Ia data. The resulting HD is corrected for selection biases and contamination from core collapse (CC) SNe, and can be used to infer cosmological parameters. This method, called "BBC" (BEAMS with Bias Corrections), includes two fitting stages. The first BBC fitting stage uses a posterior distribution that includes multiple SN likelihoods, a Monte Carlo simulation to bias-correct the fitted SALT-II parameters, and CC probabilities determined from a machine learning technique. The BBC fit determines 1) a bin-averaged HD (average distance vs. redshift), and 2) the nuisance parameters alpha and beta, which multiply the stretch and color (respectively) to standardize the SN brightness. In the second stage, the bin-averaged HD is fit to a cosmological model where priors can be imposed. We perform high precision tests of the BBC method by simulating large (150,000 event) data samples corresponding to the Dark Energy Survey Supernova Program. Our tests include three models of intrinsic scatter, each with two different CC rates. In the BBC fit, the SALT-II nuisance parameters alpha and beta are recovered to within 1% of their true values. In the cosmology fit, we determine the dark energy equation of state parameter w using a fixed value of Omega_matter as a prior: averaging over all six tests based on 6 x 150,000 = 900,000 SNe, there is a small w-bias of 0.006 +- 0.002. Finally, the BBC fitting code is publicly available in the SNANA package.

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey., Comment: 22 pages, 18 figures; Accepted by AJ, revised to incorporate referee comments

Simulations of Type Ia Supernovae (SNIa) surveys are a critical tool for correcting biases in the analysis of SNIa to infer cosmological parameters. Large scale Monte Carlo simulations include a thorough treatment of observation history, measurement noise, intrinsic scatter models and selection effects. In this paper, we improve simulations with a robust technique to evaluate the underlying populations of SNIa color and stretch that correlate with luminosity. In typical analyses, the standardized SNIa brightness is determined from linear `Tripp' relations between the light curve color and luminosity and between stretch and luminosity. However, this solution produces Hubble residual biases because intrinsic scatter and measurement noise result in measured color and stretch values that do not follow the Tripp relation. We find a $10\sigma$ bias (up to 0.3 mag) in Hubble residuals versus color and $5\sigma$ bias (up to 0.2 mag) in Hubble residuals versus stretch in a joint sample of 920 spectroscopically confirmed SNIa from PS1, SNLS, SDSS and several low-z surveys. After we determine the underlying color and stretch distributions, we use simulations to predict and correct the biases in the data. We show that removing these biases has a small impact on the low-z sample, but reduces the intrinsic scatter $\sigma_{\textrm{int}}$ from $0.101$ to $0.083$ in the combined PS1, SNLS and SDSS sample. Past estimates of the underlying populations were too broad, leading to a small bias in the equation-of-state of dark energy $w$ of $\Delta w=0.005$., Comment: Accepted by ApJL