Researcher profile

David Rubin

David Rubin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2026arXiv

Three-dimensional scene reconstruction using Roman slitless spectra

The Nancy Grace Roman Space Telescope will carry out a wide-field imaging and slitless spectroscopic survey of Type Ia Supernovae to improve our understanding of dark energy. Crucial to this endeavor is obtaining supernova spectra uncontaminated by light from their host galaxies. However, obtaining such spectra is made more difficult by the inherent problem in wide-field slitless spectroscopic surveys: the blending of spectra of close objects. The spectrum of a supernova will blend with the host galaxy, even from regions distant from the supernova on the sky. If not properly removed, this contamination will introduce systematic bias when the supernova spectra are later used to determine intrinsic supernova parameters and to infer the parameters of dark energy. To address this problem we developed an algorithm that makes use of the spectroscopic observations of the host galaxy at all available observatory roll angles to reconstruct a three-dimensional (3d; 2d spatial, 1d spectral) representation of the underlying host galaxy that accurately matches the 2d slitless spectrum of the host galaxy when projected to an arbitrary rotation angle. We call this ``scene reconstruction''. The projection of the reconstructed scene can be subtracted from an observation of a supernova to remove the contamination from the underlying host. Using simulated Roman data, we show that our method has extremely small systematic errors and significantly less random noise than if we subtracted a single perfectly aligned spectrum of the host obtained before or after the supernova was visible.

preprint2022arXiv

Evaluating and Optimizing a Slitless Prism for Nancy Grace Roman Space Telescope SN Cosmology

This work presents a set of studies addressing the use of the low-dispersion slitless prism on Roman for SN spectroscopy as part of the Roman High Latitude Time Domain Survey (HLTDS). We find SN spectral energy distributions including prism data carry more information than imaging alone at fixed total observing time, improving redshift measurements and sub-typing of SNe. The Roman field of view will typically include ~ 10 SNe Ia at observable redshifts at a range of phases (the multiplexing of host galaxies is much greater as they are always present), building up SN spectral time series without targeted observations. We show that fitting these time series extracts more information than stacking the data over all the phases, resulting in a large improvement in precision for SN Ia subclassification measurements. A prism on Roman thus significantly enhances scientific opportunities for the mission, and is particularly important for the Roman SN cosmology program to provide the systematics-controlled measurement that is a focus of the Roman dark energy mission. Optimizing the prism parameters, we conclude that the blue cutoff should be set as blue as the prism image quality allows (~ 7500A), the red cutoff should be set to ~ 18000A to minimize thermal background, and the two-pixel dispersion should be >~ 70.

preprint2022arXiv

Uniform Recalibration of Common Spectrophotometry Standard Stars onto the CALSPEC System using the SuperNova Integral Field Spectrograph

We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the HST-based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300 A to 9400 A as calibration for the Nearby Supernova Factory cosmology experiment. In total, this analysis used 4289 standard-star spectra taken on photometric nights. As a modern cosmology analysis, all pre-submission methodological decisions were made with the flux scale and external comparison results blinded. The large number of spectra per star allows us to treat the wavelength-by-wavelength calibration for all nights simultaneously with a Bayesian hierarchical model, thereby enabling a consistent treatment of the Type Ia supernova cosmology analysis and the calibration on which it critically relies. We determine the typical per-observation repeatability (median 14 mmag for exposures >~ 5 s), the Maunakea atmospheric transmission distribution (median dispersion of 7 mmag with uncertainty 1 mmag), and the scatter internal to our CALSPEC reference stars (median of 8 mmag). We also check our standards against literature filter photometry, finding generally good agreement over the full 12-magnitude range. Overall, the mean of our system is calibrated to the mean of CALSPEC at the level of ~ 3 mmag. With our large number of observations, careful crosschecks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys.

preprint2020arXiv

Accelerated Bayesian Optimization throughWeight-Prior Tuning

Bayesian optimization (BO) is a widely-used method for optimizing expensive (to evaluate) problems. At the core of most BO methods is the modeling of the objective function using a Gaussian Process (GP) whose covariance is selected from a set of standard covariance functions. From a weight-space view, this models the objective as a linear function in a feature space implied by the given covariance K, with an arbitrary Gaussian weight prior ${\bf w} \sim \mathcal{N} ({\bf 0}, {\bf I})$. In many practical applications there is data available that has a similar (covariance) structure to the objective, but which, having different form, cannot be used directly in standard transfer learning. In this paper we show how such auxiliary data may be used to construct a GP covariance corresponding to a more appropriate weight prior for the objective function. Building on this, we show that we may accelerate BO by modeling the objective function using this (learned) weight prior, which we demonstrate on both test functions and a practical application to short-polymer fibre manufacture.

preprint2020arXiv

Does Gravity Fall Down? Evidence for Gravitational Wave Deflection Along the Line of Sight to GW 170817

We present a novel test of general relativity (GR): measuring the geometric component of the time delay due to gravitational lensing. GR predicts that photons and gravitational waves follow the same geodesic paths and thus experience the same geometric time delay. We show that for typical systems, the time delays are tens of seconds, and thus can dominate over astrophysical delays in the timing of photon emission. For the case of GW 170817, we use a multi-plane lensing code to evaluate the time delay due to four massive halos along the line of sight. From literature mass and distance measurements of these halos, we establish at high confidence (significantly greater than 5 sigma) that the gravitational waves of GW 170817 underwent gravitational deflection to arrive within 1.7 seconds of the photons.

preprint2020arXiv

Evaluating the Calibration of SN Ia Anchor Datasets with a Bayesian Hierarchical Model

Inter-survey calibration remains an important systematic uncertainty in cosmological studies using type Ia supernova (SNe Ia). Ideally, each survey would measure its system throughputs, for instance with bandpass measurements combined with observations of well-characterized spectrophotometric standard stars; however, many important nearby-SN surveys have not done this. We recalibrate these surveys by tying their tertiary survey stars to Pan-STARRS1 g, r, and i, and SDSS/CSP u. This improves upon previous recalibration efforts by taking the spatially variable zeropoints of each telescope/camera into account, and applying improved color transformations in the surveys' natural instrumental photometric systems. Our analysis uses a global hierarchical model of the data which produces a covariance matrix of magnitude offsets and bandpass shifts, quantifying and reducing the systematic uncertainties in the calibration. We call our method CROSS-CALIBration with a Uniform Reanalysis (X-CALIBUR). This approach gains not only from a sophisticated analysis, but also from simply tying our calibration to more color calibrators, rather than just the one color calibrator (BD+17 4708) as many previous efforts have done. The results presented here have the potential to help understand and improve calibration uncertainties upcoming SN Ia cosmological analyses.

preprint2020arXiv

Is the expansion of the universe accelerating? All signs still point to yes a local dipole anisotropy cannot explain dark energy

Type Ia supernovae (SNe Ia) provided the first strong evidence that the expansion of the universe is accelerating. With SN samples now more than ten times larger than those used for the original discovery and joined by other cosmological probes, this discovery is on even firmer ground. Two recent, related studies (Nielsen et al. 2016 and Colin et al. 2019, hereafter N16 and C19, respectively) have claimed to undermine the statistical significance of the SN Ia constraints. Rubin & Hayden (2016) (hereafter RH16) showed N16 made an incorrect assumption about the distributions of SN Ia light-curve parameters, while C19 also fails to remove the impact of the motion of the solar system from the SN redshifts, interpreting the resulting errors as evidence of a dipole in the deceleration parameter. Building on RH16, we outline the errors C19 makes in their treatment of the data and inference on cosmological parameters. Reproducing the C19 analysis with our proposed fixes, we find that the dipole parameters have little effect on the inferred cosmological parameters. We thus affirm the conclusion of RH16: the evidence for acceleration is secure.

preprint2020arXiv

SN2019dge: a Helium-rich Ultra-Stripped Envelope Supernova

We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture.

preprint2019arXiv

Constraining the Dimensionality of SN Ia Spectral Variation with Twins

SNe Ia continue to play a key role in cosmological measurements. Their interpretation over a range in redshift requires a rest-frame spectral energy distribution model. For practicality, these models are parameterized with a limited number of parameters and are trained using linear or nonlinear dimensionality reduction. This work focuses on the related problem of estimating the number of parameters underlying SN Ia spectral variation (the dimensionality). I present a technique for using the properties of high-dimensional space and the counting statistics of &#34;twin&#34; SNe Ia to estimate this dimensionality. Applying this method to the supernova pairings from Fakhouri et al. (2015) shows that a modest number of parameters (three to five, not including extinction) explain those data well. The analysis also finds that the intrinsic parameters are approximately Gaussian-distributed. The limited number of parameters hints that improved SED models are possible that may enable substantial reductions in SN cosmological uncertainties with current and near-term datasets.

preprint2019arXiv

Precise Mass Determination of SPT-CL J2106-5844, the Most Massive Cluster at z>1

We present a detailed high-resolution weak-lensing (WL) study of SPT-CL J2106-5844 at z=1.132, claimed to be the most massive system discovered at z > 1 in the South Pole Telescope Sunyaev-Zel&#39;dovich (SPT-SZ) survey. Based on the deep imaging data from the Advanced Camera for Surveys and Wide Field Camera 3 on-board the Hubble Space Telescope, we find that the cluster mass distribution is asymmetric, composed of a main clump and a subclump ~640 kpc west thereof. The central clump is further resolved into two smaller northwestern and southeastern substructures separated by ~150 kpc. We show that this rather complex mass distribution is more consistent with the cluster galaxy distribution than a unimodal distribution as previously presented. The northwestern substructure coincides with the BCG and X-ray peak while the southeastern one agrees with the location of the number density peak. These morphological features and the comparison with the X-ray emission suggest that the cluster might be a merging system. We estimate the virial mass of the cluster to be $M_{200c} = (10.4^{+3.3}_{-3.0}\pm1.0)~\times~10^{14}~M_{\odot}$, where the second error bar is the systematic uncertainty. Our result confirms that the cluster SPT-CL J2106-5844 is indeed the most massive cluster at z>1 known to date. We demonstrate the robustness of this mass estimate by performing a number of tests with different assumptions on the centroids, mass-concentration relations, and sample variance.