Steven A. Rodney

The James Webb Space Telescope (JWST) North Ecliptic Pole (NEP) Time-domain Field (TDF) is a >14′ diameter field optimized for multiwavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over three cycles (the “TREASUREHUNT” program), deep images were obtained with the Wide Field Camera on the Advanced Camera for Surveys in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ∼88 arcmin², with time intervals between individual epochs that range between 1 day and 4+ yr. To a depth of m _AB ≃ 29.5 mag (F606W), we present the discovery of 12 transients and 190 variable candidates. For the variable candidates, we demonstrate that Gaussian statistics are applicable and estimate that ∼80 are false positives. The majority of the transients will be supernovae, although at least two are likely quasars. Most variable candidates are active galactic nuclei (AGNs), where we find 0.42% of the general z ≲ 6 field galaxy population to vary at the ∼3σ level. Based on a 5 yr time frame, this translates into a random supernova areal density of up to ∼0.07 transients arcmin⁻² (∼245 deg⁻²) per epoch and a variable AGN areal density of ∼1.25 variables arcmin⁻² (∼4500 deg⁻²) to these depths.

The gravitationally lensed supernova Refsdal appeared in multiple images produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted that an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measurement of the expansion rate of the Universe, quantified by the Hubble constant ( H ₀ ). Using eight cluster lens models, we infer H 0 = 64.8 − 4.3 + 4.4  kilometers per second per megaparsec . Using the two models most consistent with the observations, we find H 0 = 66.6 − 3.3 + 4.1  kilometers per second per megaparsec . The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster.

A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models—including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is 2.5× larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range 2000–11,000 Å (1800 Å redder) and reduced uncertainties compared to SALT2, enabling accurate use of low-z I and iz photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low-z Foundation and CfA3 samples by 15% and 10%, respectively. To check for potential systematic uncertainties, we compare distances of low (0.01 < z < 0.2) and high (0.4 < z < 0.6) redshift SNe in the training compilation, finding an insignificant 3 ± 14 mmag shift between SALT2.4 and SALT3.K21. While the SALT3.K21 model was trained on optical data, our method can be used to build a model for rest-frame NIR samples from the Roman Space Telescope. Our open-source training code, public training data, model, and documentation are available at https://saltshaker.readthedocs.io/en/latest/, and the model is integrated into the sncosmo and SNANA software packages.

In the next decade, transient searches from the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope will increase the sample of known Type Ia supernovae (SNe Ia) from ∼10³ to 10⁵. With this reduction of statistical uncertainties on cosmological measurements, new methods are needed to reduce systematic uncertainties. Characterizing the underlying spectroscopic evolution of SN Ia remains a major systematic uncertainty in current cosmological analyses, motivating a new simulation tool for the next era of SN Ia cosmology: Build Your Own Spectral Energy Distribution (byosed). byosed is used within the SNANA framework to simulate light curves by applying spectral variations to model SEDs, enabling flexible testing of possible systematic shifts in SN Ia distance measurements. We test the framework by comparing a nominal Roman SN Ia survey simulation using a baseline SED model to simulations using SEDs perturbed with byosed, and investigating the impact of ignoring specific SED features in the analysis. These features include semiempirical models of two possible, predicted relationships: between SN ejecta velocity and light-curve observables, and a redshift-dependent relationship between SN Hubble residuals and host-galaxy mass. We analyze each byosed simulation using the SALT2 and BEAMS with Bias Corrections framework, and estimate changes in the measured value of the dark-energy equation-of-state parameter, w. We find a difference of Δw = −0.023 for SN velocity and Δw = 0.021 for redshift-evolving host mass when compared to simulations without these features. By using byosed for SN Ia cosmology simulations, future analyses (e.g., the Rubin and Roman SN Ia samples) will have greater flexibility to constrain or reduce such SN Ia modeling uncertainties.

One of the primary mission objectives of the Roman Space Telescope is to investigate the nature of dark energy with a variety of methods. Observations of Type I supernovae (SNe Ia) will be one of the principal anchors of the Roman cosmology program through traditional luminosity distance measurements. This SNe Ia cosmology program can provide another valuable cosmological probe, without altering the strategy of the mission: time delay cosmography with gravitationally lensed supernova (SN). In this work, we forecast lensed SN cosmology constraints with the Roman Space Telescope, while providing useful tools for future work. Using the anticipated characteristics of the Roman SNe Ia survey, we have constructed mock catalogs of expected resolved lensing systems, as well as strongly lensed Type Ia and core-collapse (CC) SN light curves, including microlensing effects. We predict Roman will find ∼11 lensed SNe Ia and ∼20 CCSNe, depending on the survey strategy. Next, we estimate the time delay precision obtainable with Roman (Ia: ∼2 days, CC: ∼3 days), and use a Fisher matrix analysis to derive projected constraints on H ₀,Ω _m , and the dark energy equation of state, w, for each SNe Ia survey strategy. A strategy optimized for the discovery of high-redshift SNe Ia is preferred when considering the constraints possible from both SNe Ia and lensed SN cosmology, also delivering ∼1.5 times more lensed SNe than other proposed survey strategies.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope (HST) Treasury program taking data from 2018 to 2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in Wide Field Camera 3/IR F105W, F125W, and F160W and Advanced Camera for Surveys/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multiwavelength data sets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies that would otherwise lie below HST detection limits and model foreground clusters to study the properties of dark matter and galaxy assembly. The expanded area will provide the first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient follow-up. BUFFALO data including mosaics, value-added catalogs, and cluster mass distribution models will be released via MAST on a regular basis as the observations and analysis are completed for the six individual clusters.

We present analytical reconstructions of SN Ia delay time distributions (DTDs) by way of two independent methods: by a Markov Chain Monte Carlo best-fit technique comparing the volumetric SN Ia rate history to today’s compendium cosmic star formation history, and second through a maximum likelihood analysis of the star formation rate histories of individual galaxies in the GOODS/CANDELS field, in comparison to their resultant SN Ia yields. We adopt a flexible skew-normal DTD model, which could match a wide range of physically motivated DTD forms. We find a family of solutions that are essentially exponential DTDs, similar in shape to the β ≈ −1 power-law DTDs, but with more delayed events (>1 Gyr in age) than prompt events (<1 Gyr). Comparing these solutions to delay time measures separately derived from field galaxies and galaxy clusters, we find the skew-normal solutions can accommodate both without requiring a different DTD form in different environments. These model fits are generally inconsistent with results from single-degenerate binary population synthesis models, and are seemingly supportive of double-degenerate progenitors for most SN Ia events.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ≈200 arcmin². These clusters were selected to be excellent lenses, and we find similar high-redshift sample sizes and magnitude distributions as the Cluster Lensing And Supernova survey with Hubble (CLASH). We discover 257, 57, and eight candidate galaxies at z ∼ 6, 7, and 8 respectively, (322 in total). The observed (lensed) magnitudes of the z ∼ 6 candidates are as bright as AB mag ∼23, making them among the brightest known at these redshifts, comparable with discoveries from much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.