-
pop-cosmos: Scaleable inference of galaxy properties and redshifts with a data-driven population model
Authors:
Stephen Thorp,
Justin Alsing,
Hiranya V. Peiris,
Sinan Deger,
Daniel J. Mortlock,
Boris Leistedt,
Joel Leja,
Arthur Loureiro
Abstract:
We present an efficient Bayesian method for estimating individual photometric redshifts and galaxy properties under a pre-trained population model (pop-cosmos) that was calibrated using purely photometric data. This model specifies a prior distribution over 16 stellar population synthesis (SPS) parameters using a score-based diffusion model, and includes a data model with detailed treatment of neb…
▽ More
We present an efficient Bayesian method for estimating individual photometric redshifts and galaxy properties under a pre-trained population model (pop-cosmos) that was calibrated using purely photometric data. This model specifies a prior distribution over 16 stellar population synthesis (SPS) parameters using a score-based diffusion model, and includes a data model with detailed treatment of nebular emission. We use a GPU-accelerated affine invariant ensemble sampler to achieve fast posterior sampling under this model for 292,300 individual galaxies in the COSMOS2020 catalog, leveraging a neural network emulator (Speculator) to speed up the SPS calculations. We apply both the pop-cosmos population model and a baseline prior inspired by Prospector-$α$, and compare these results to published COSMOS2020 redshift estimates from the widely-used EAZY and LePhare codes. For the $\sim 12,000$ galaxies with spectroscopic redshifts, we find that pop-cosmos yields redshift estimates that have minimal bias ($\sim10^{-4}$), high accuracy ($σ_\text{MAD}=7\times10^{-3}$), and a low outlier rate ($1.6\%$). We show that the pop-cosmos population model generalizes well to galaxies fainter than its $r<25$ mag training set. The sample we have analyzed is $\gtrsim3\times$ larger than has previously been possible via posterior sampling with a full SPS model, with average throughput of 15 GPU-sec per galaxy under the pop-cosmos prior, and 0.6 GPU-sec per galaxy under the Prospector prior. This paves the way for principled modeling of the huge catalogs expected from upcoming Stage IV galaxy surveys.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Variational Inference for Acceleration of SN Ia Photometric Distance Estimation with BayeSN
Authors:
Ana Sofía M. Uzsoy,
Stephen Thorp,
Matthew Grayling,
Kaisey S. Mandel
Abstract:
Type Ia supernovae (SNe Ia) are standarizable candles whose observed light curves can be used to infer their distances, which can in turn be used in cosmological analyses. As the quantity of observed SNe Ia grows with current and upcoming surveys, increasingly scalable analyses are necessary to take full advantage of these new datasets for precise estimation of cosmological parameters. Bayesian in…
▽ More
Type Ia supernovae (SNe Ia) are standarizable candles whose observed light curves can be used to infer their distances, which can in turn be used in cosmological analyses. As the quantity of observed SNe Ia grows with current and upcoming surveys, increasingly scalable analyses are necessary to take full advantage of these new datasets for precise estimation of cosmological parameters. Bayesian inference methods enable fitting SN Ia light curves with robust uncertainty quantification, but traditional posterior sampling using Markov Chain Monte Carlo (MCMC) is computationally expensive. We present an implementation of variational inference (VI) to accelerate the fitting of SN Ia light curves using the BayeSN hierarchical Bayesian model for time-varying SN Ia spectral energy distributions (SEDs). We demonstrate and evaluate its performance on both simulated light curves and data from the Foundation Supernova Survey with two different forms of surrogate posterior -- a multivariate normal and a custom multivariate zero-lower-truncated normal distribution -- and compare them with the Laplace Approximation and full MCMC analysis. To validate of our variational approximation, we calculate the pareto-smoothed importance sampling (PSIS) diagnostic, and perform variational simulation-based calibration (VSBC). The VI approximation achieves similar results to MCMC but with an order-of-magnitude speedup for the inference of the photometric distance moduli. Overall, we show that VI is a promising method for scalable parameter inference that enables analysis of larger datasets for precision cosmology.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
JWST Photometric Time-Delay and Magnification Measurements for the Triply-Imaged Type Ia "Supernova H0pe" at z = 1.78
Authors:
J. D. R. Pierel,
B. L. Frye,
M. Pascale,
G. B. Caminha,
W. Chen,
S. Dhawan,
D. Gilman,
M. Grayling,
S. Huber,
P. Kelly,
S. Thorp,
N. Arendse,
S. Birrer,
M. Bronikowski,
R. Canameras,
D. Coe,
S. H. Cohen,
C. J. Conselice,
S. P. Driver,
J. C. J. Dsilva,
M. Engesser,
N. Foo,
C. Gall,
N. Garuda,
C. Grillo
, et al. (38 additional authors not shown)
Abstract:
Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is locate…
▽ More
Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is located at $z=1.783$, and is the first SN Ia with sufficient light curve sampling and long enough time delays for an $H_0$ inference. Here we present photometric time-delay measurements and SN properties of SN H0pe. Using JWST/NIRCam photometry we measure time delays of $Δt_{ab}=-116.6^{+10.8}_{-9.3}$ and $Δt_{cb}=-48.6^{+3.6}_{-4.0}$ observer-frame days relative to the last image to arrive (image 2b; all uncertainties are $1σ$), which corresponds to a $\sim5.6\%$ uncertainty contribution for $H_0$ assuming $70 \rm{km s^{-1} Mpc^{-1}}$. We also constrain the absolute magnification of each image to $μ_{a}=4.3^{+1.6}_{-1.8}$, $μ_{b}=7.6^{+3.6}_{-2.6}$, $μ_{c}=6.4^{+1.6}_{-1.5}$ by comparing the observed peak near-IR magnitude of SN H0pe to the non-lensed population of SNe Ia.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Using Rest-Frame Optical and NIR Data from the RAISIN Survey to Explore the Redshift Evolution of Dust Laws in SN Ia Host Galaxies
Authors:
Stephen Thorp,
Kaisey S. Mandel,
David O. Jones,
Robert P. Kirshner,
Peter M. Challis
Abstract:
We use rest-frame optical and near-infrared (NIR) observations of 42 Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project at low-$z$ and 37 from the RAISIN Survey at high-$z$ to investigate correlations between SN Ia host galaxy dust, host mass, and redshift. This is the first time the SN Ia host galaxy dust extinction law at high-$z$ has been estimated using combined optical and rest-f…
▽ More
We use rest-frame optical and near-infrared (NIR) observations of 42 Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project at low-$z$ and 37 from the RAISIN Survey at high-$z$ to investigate correlations between SN Ia host galaxy dust, host mass, and redshift. This is the first time the SN Ia host galaxy dust extinction law at high-$z$ has been estimated using combined optical and rest-frame NIR data ($YJ$-band). We use the BayeSN hierarchical model to leverage the data's wide rest-frame wavelength range (extending to $\sim$1.0-1.2 microns for the RAISIN sample at $0.2\lesssim z\lesssim0.6$). By contrasting the RAISIN and CSP data, we constrain the population distributions of the host dust $R_V$ parameter for both redshift ranges. We place a limit on the difference in population mean $R_V$ between RAISIN and CSP of $-1.16<Δμ(R_V)<1.38$ with 95% posterior probability. For RAISIN we estimate $μ(R_V)=2.58\pm0.57$, and constrain the population standard deviation to $σ(R_V)<0.90~[2.42]$ at the 68 [95]% level. Given that we are only able to constrain the size of the low- to high-$z$ shift in $μ(R_V)$ to $\lesssim1.4$ - which could still propagate to a substantial bias in the equation of state parameter $w$ - these and other recent results motivate continued effort to obtain rest-frame NIR data at low and high redshifts (e.g. using the Roman Space Telescope).
△ Less
Submitted 23 April, 2024; v1 submitted 28 February, 2024;
originally announced February 2024.
-
pop-cosmos: A comprehensive picture of the galaxy population from COSMOS data
Authors:
Justin Alsing,
Stephen Thorp,
Sinan Deger,
Hiranya Peiris,
Boris Leistedt,
Daniel Mortlock,
Joel Leja
Abstract:
We present pop-cosmos: a comprehensive model characterizing the galaxy population, calibrated to $140,938$ galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in $26$ bands from the ultra-violet to the infra-red. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (par…
▽ More
We present pop-cosmos: a comprehensive model characterizing the galaxy population, calibrated to $140,938$ galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in $26$ bands from the ultra-violet to the infra-red. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (parameterized by a flexible score-based diffusion model); a state-of-the-art stellar population synthesis (SPS) model connecting galaxies' instrinsic properties to their photometry; and a data-model for the observation, calibration and selection processes. By minimizing the optimal transport distance between synthetic and real data we are able to jointly fit the population- and data-models, leading to robustly calibrated population-level inferences that account for parameter degeneracies, photometric noise and calibration, and selection effects. We present a number of key predictions from our model of interest for cosmology and galaxy evolution, including the mass function and redshift distribution; the mass-metallicity-redshift and fundamental metallicity relations; the star-forming sequence; the relation between dust attenuation and stellar mass, star formation rate and attenuation-law index; and the relation between gas-ionization and star formation. Our model encodes a comprehensive picture of galaxy evolution that faithfully predicts galaxy colors across a broad redshift ($z<4$) and wavelength range.
△ Less
Submitted 1 February, 2024;
originally announced February 2024.
-
Data-Space Validation of High-Dimensional Models by Comparing Sample Quantiles
Authors:
Stephen Thorp,
Hiranya V. Peiris,
Daniel J. Mortlock,
Justin Alsing,
Boris Leistedt,
Sinan Deger
Abstract:
We present a simple method for assessing the predictive performance of high-dimensional models directly in data space when only samples are available. Our approach is to compare the quantiles of observables predicted by a model to those of the observables themselves. In cases where the dimensionality of the observables is large (e.g. multiband galaxy photometry), we advocate that the comparison is…
▽ More
We present a simple method for assessing the predictive performance of high-dimensional models directly in data space when only samples are available. Our approach is to compare the quantiles of observables predicted by a model to those of the observables themselves. In cases where the dimensionality of the observables is large (e.g. multiband galaxy photometry), we advocate that the comparison is made after projection onto a set of principal axes to reduce the dimensionality. We demonstrate our method on a series of two-dimensional examples. We then apply it to results from a state-of-the-art generative model for galaxy photometry (pop-cosmos; arXiv:2402.00935) that generates predictions of colors and magnitudes by forward simulating from a 16-dimensional distribution of physical parameters represented by a score-based diffusion model. We validate the predictive performance of this model directly in a space of nine broadband colors. Although motivated by this specific example, the techniques we present will be broadly useful for evaluating the performance of flexible, non-parametric population models of this kind, and can be readily applied to any setting where two sets of samples are to be compared.
△ Less
Submitted 16 February, 2024; v1 submitted 1 February, 2024;
originally announced February 2024.
-
Scalable hierarchical BayeSN inference: Investigating dependence of SN Ia host galaxy dust properties on stellar mass and redshift
Authors:
Matthew Grayling,
Stephen Thorp,
Kaisey S. Mandel,
Suhail Dhawan,
Ana Sofia M. Uzsoy,
Benjamin M. Boyd,
Erin E. Hayes,
Sam M. Ward
Abstract:
We apply the hierarchical probabilistic SED model BayeSN to analyse a sample of 475 SNe Ia (0.015 < z < 0.4) from Foundation, DES3YR and PS1MD to investigate the properties of dust in their host galaxies. We jointly infer the dust law $R_V$ population distributions at the SED level in high- and low-mass galaxies simultaneously with dust-independent, intrinsic differences. We find an intrinsic mass…
▽ More
We apply the hierarchical probabilistic SED model BayeSN to analyse a sample of 475 SNe Ia (0.015 < z < 0.4) from Foundation, DES3YR and PS1MD to investigate the properties of dust in their host galaxies. We jointly infer the dust law $R_V$ population distributions at the SED level in high- and low-mass galaxies simultaneously with dust-independent, intrinsic differences. We find an intrinsic mass step of $-0.049\pm0.016$ mag, at a significance of 3.1$σ$, when allowing for a constant intrinsic, achromatic magnitude offset. We additionally apply a model allowing for time- and wavelength-dependent intrinsic differences between SNe Ia in different mass bins, finding $\sim$2$σ$ differences in magnitude and colour around peak and 4.5$σ$ differences at later times. These intrinsic differences are inferred simultaneously with a difference in population mean $R_V$ of $\sim$2$σ$ significance, demonstrating that both intrinsic and extrinsic differences may play a role in causing the host galaxy mass step. We also consider a model which allows the mean of the $R_V$ distribution to linearly evolve with redshift but find no evidence for any evolution - we infer the gradient of this relation $η_R = -0.38\pm0.70$. In addition, we discuss in brief a new, GPU-accelerated Python implementation of BayeSN suitable for application to large surveys which is publicly available and can be used for future cosmological analyses; this code can be found here: https://github.com/bayesn/bayesn.
△ Less
Submitted 29 April, 2024; v1 submitted 16 January, 2024;
originally announced January 2024.
-
GausSN: Bayesian Time-Delay Estimation for Strongly Lensed Supernovae
Authors:
Erin E. Hayes,
Stephen Thorp,
Kaisey S. Mandel,
Nikki Arendse,
Matthew Grayling,
Suhail Dhawan
Abstract:
We present GausSN, a Bayesian semi-parametric Gaussian Process (GP) model for time-delay estimation with resolved systems of gravitationally lensed supernovae (glSNe). GausSN models the underlying light curve non-parametrically using a GP. Without assuming a template light curve for each SN type, GausSN fits for the time delays of all images using data in any number of wavelength filters simultane…
▽ More
We present GausSN, a Bayesian semi-parametric Gaussian Process (GP) model for time-delay estimation with resolved systems of gravitationally lensed supernovae (glSNe). GausSN models the underlying light curve non-parametrically using a GP. Without assuming a template light curve for each SN type, GausSN fits for the time delays of all images using data in any number of wavelength filters simultaneously. We also introduce a novel time-varying magnification model to capture the effects of microlensing alongside time-delay estimation. In this analysis, we model the time-varying relative magnification as a sigmoid function, as well as a constant for comparison to existing time-delay estimation approaches. We demonstrate that GausSN provides robust time-delay estimates for simulations of glSNe from the Nancy Grace Roman Space Telescope and the Vera C. Rubin Observatory's Legacy Survey of Space and Time (Rubin-LSST). We find that up to 43.6% of time-delay estimates from Roman and 52.9% from Rubin-LSST have fractional errors of less than 5%. We then apply GausSN to SN Refsdal and find the time delay for the fifth image is consistent with the original analysis, regardless of microlensing treatment. Therefore, GausSN maintains the level of precision and accuracy achieved by existing time-delay extraction methods with fewer assumptions about the underlying shape of the light curve than template-based approaches, while incorporating microlensing into the statistical error budget rather than requiring post-processing to account for its systematic uncertainty. GausSN is scalable for time-delay cosmography analyses given current projections of glSNe discovery rates from Rubin-LSST and Roman.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
Bird-Snack: Bayesian Inference of dust law $R_V$ Distributions using SN Ia Apparent Colours at peaK
Authors:
Sam M. Ward,
Suhail Dhawan,
Kaisey S. Mandel,
Matthew Grayling,
Stephen Thorp
Abstract:
To reduce systematic uncertainties in Type Ia supernova (SN Ia) cosmology, the host galaxy dust law shape parameter, $R_V$, must be accurately constrained. We thus develop a computationally-inexpensive pipeline, Bird-Snack, to rapidly infer dust population distributions from optical-near infrared SN colours at peak brightness, and determine which analysis choices significantly impact the populatio…
▽ More
To reduce systematic uncertainties in Type Ia supernova (SN Ia) cosmology, the host galaxy dust law shape parameter, $R_V$, must be accurately constrained. We thus develop a computationally-inexpensive pipeline, Bird-Snack, to rapidly infer dust population distributions from optical-near infrared SN colours at peak brightness, and determine which analysis choices significantly impact the population mean $R_V$ inference, $μ_{R_V}$. Our pipeline uses a 2D Gaussian process to measure peak $BVriJH$ apparent magnitudes from SN light curves, and a hierarchical Bayesian model to simultaneously constrain population distributions of intrinsic and dust components. Fitting a low-to-moderate-reddening sample of 65 low-redshift SNe yields $μ_{R_V}=2.61^{+0.38}_{-0.35}$, with $68\%(95\%)$ posterior upper bounds on the population dispersion, $σ_{R_V}<0.92(1.96)$. This result is robust to various analysis choices, including: the model for intrinsic colour variations, fitting the shape hyperparameter of a gamma dust extinction distribution, and cutting the sample based on the availability of data near peak. However, these choices may be important if statistical uncertainties are reduced. With larger near-future optical and near-infrared SN samples, Bird-Snack can be used to better constrain dust distributions, and investigate potential correlations with host galaxy properties. Bird-Snack is publicly available; the modular infrastructure facilitates rapid exploration of custom analysis choices, and quick fits to simulated datasets, for better interpretation of real-data inferences.
△ Less
Submitted 11 October, 2023;
originally announced October 2023.
-
Keck Infrared Transient Survey I: Survey Description and Data Release 1
Authors:
S. Tinyanont,
R. J. Foley,
K. Taggart,
K. W. Davis,
N. LeBaron,
J. E. Andrews,
M. J. Bustamante-Rosell,
Y. Camacho-Neves,
R. Chornock,
D. A. Coulter,
L. Galbany,
S. W. Jha,
C. D. Kilpatrick,
L. A. Kwok,
C. Larison,
J. R. Pierel,
M. R. Siebert,
G. Aldering,
K. Auchettl,
J. S. Bloom,
S. Dhawan,
A. V. Filippenko,
K. D. French,
A. Gagliano,
M. Grayling
, et al. (13 additional authors not shown)
Abstract:
We present the Keck Infrared Transient Survey (KITS), a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman…
▽ More
We present the Keck Infrared Transient Survey (KITS), a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roma} SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. The first data release includes data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients $<$17 mag in a red optical band (usually ZTF r or ATLAS o bands); a volume-limited sample including all transients within redshift $z < 0.01$ ($D \approx 50$ Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete to $z = 0.005$. All completeness numbers will rise with the inclusion of data from other telescopes in future data releases. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events (TDEs), luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction using Pypeit, which requires minimal human interaction to ensure reproducibility.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
The Magnificent Five Images of Supernova Refsdal: Time Delay and Magnification Measurements
Authors:
Patrick L. Kelly,
Steven Rodney,
Tommaso Treu,
Simon Birrer,
Vivien Bonvin,
Luc Dessart,
Ryan J. Foley,
Alexei V. Filippenko,
Daniel Gilman,
Saurabh Jha,
Jens Hjorth,
Kaisey Mandel,
Martin Millon,
Justin Pierel,
Stephen Thorp,
Adi Zitrin,
Tom Broadhurst,
Wenlei Chen,
Jose M. Diego,
Alan Dressler,
Or Graur,
Mathilde Jauzac,
Matthew A. Malkan,
Curtis McCully,
Masamune Oguri
, et al. (6 additional authors not shown)
Abstract:
In late 2014, four images of Supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8 arcseconds away in the field. The observed reappearance in late 2015 makes it possible to carry out R…
▽ More
In late 2014, four images of Supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8 arcseconds away in the field. The observed reappearance in late 2015 makes it possible to carry out Refsdal's (1964) original proposal to use a multiply imaged SN to measure the Hubble constant H0, since the time delay between appearances should vary inversely with H0. Moreover, the position, brightness, and timing of the reappearance enable a novel test of the blind predictions of galaxy-cluster models, which are typically constrained only by the positions of multiply imaged galaxies. We have developed a new photometry pipeline that uses DOLPHOT to measure the fluxes of the five images of SN Refsdal from difference images. We apply four separate techniques to perform a blind measurement of the relative time delays and magnification ratios (mu_i/mu_1) between the last image SX and the earlier images S1-S4. We measure the relative time delay of SX-S1 to be 376.0+5.6-5.5 days and the relative magnification to be 0.30+0.05-0.03. This corresponds to a 1.5% precision on the time delay and 17% precision for the magnification ratios, and includes uncertainties due to millilensing and microlensing. In an accompanying paper, we place initial and blind constraints on the value of the Hubble constant.
△ Less
Submitted 10 May, 2023;
originally announced May 2023.
-
Constraints on the Hubble constant from Supernova Refsdal's reappearance
Authors:
Patrick L. Kelly,
Steven Rodney,
Tommaso Treu,
Masamune Oguri,
Wenlei Chen,
Adi Zitrin,
Simon Birrer,
Vivien Bonvin,
Luc Dessart,
Jose M. Diego,
Alexei V. Filippenko,
Ryan J. Foley,
Daniel Gilman,
Jens Hjorth,
Mathilde Jauzac,
Kaisey Mandel,
Martin Millon,
Justin Pierel,
Keren Sharon,
Stephen Thorp,
Liliya Williams,
Tom Broadhurst,
Alan Dressler,
Or Graur,
Saurabh Jha
, et al. (5 additional authors not shown)
Abstract:
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 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 measuremen…
▽ More
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 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 (H0). Using eight cluster lens models, we infer H0 = 64.8 +4.4-4.3 km / s / Mpc, where Mpc is the megaparsec. Using the two models most consistent with the observations, we find H0 = 66.6 +4.1-3.3 km / s / Mpc. The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster.
△ Less
Submitted 2 September, 2023; v1 submitted 10 May, 2023;
originally announced May 2023.
-
Graphene Oxide Photoreduction Recovers Graphene Hot Electron Cooling Dynamics
Authors:
Alden N. Bradley,
Spencer G. Thorp,
Gina Mayonado,
Edward Elliott,
Matt W. Graham
Abstract:
Reduced graphene oxide (rGO) is a bulk-processable quasi-amorphous 2D material with broad spectral coverage and fast electronic response. rGO sheets are suspended in a polymer matrix and sequentially photoreduced while measuring the evolving optical spectra and ultrafast electron relaxation dynamics. Photoreduced rGO yields optical absorption spectra that fit with the same Fano lineshape parameter…
▽ More
Reduced graphene oxide (rGO) is a bulk-processable quasi-amorphous 2D material with broad spectral coverage and fast electronic response. rGO sheets are suspended in a polymer matrix and sequentially photoreduced while measuring the evolving optical spectra and ultrafast electron relaxation dynamics. Photoreduced rGO yields optical absorption spectra that fit with the same Fano lineshape parameters as monolayer graphene. With increasing photoreduction time, rGO transient absorption kinetics accelerate monotonically, reaching an optimal point that matches the hot electron cooling in graphene. All stages of rGO ultrafast kinetics are simulated with a hot-electron cooling model mediated by disorder-assisted supercollisions. While the rGO room temperature 0.31 ps$^{-1}$ electronic cooling rate matches monolayer graphene, subsequent photoreduction can rapidly increase the rate by ~10-12$\times$. Such accelerated supercollision rates imply a reduced mean-free scattering length caused by photoionized point-defects on the rGO sp$^2$ sub-lattice. For visible range excitations of rGO, photoreduction shows three increasing spectral peaks that match graphene quantum dot (GQD) transitions, while a broad peak from oxygenated defect edge states shrinks. These three confined GQD states donate their hot carriers to the graphene sub-lattice with a 0.17 ps rise-time that accelerates with photoreduction. Collectively, many desirable photophysical properties of 2D graphene are replicated through selectively reducing rGO scaffolded within a 3D bulk polymeric network.
△ Less
Submitted 30 January, 2023;
originally announced January 2023.
-
The DEHVILS Survey Overview and Initial Data Release: High-Quality Near-Infrared Type Ia Supernova Light Curves at Low Redshift
Authors:
Erik R. Peterson,
David O. Jones,
Daniel Scolnic,
Bruno O. Sánchez,
Aaron Do,
Adam G. Riess,
Sam M. Ward,
Arianna Dwomoh,
Thomas de Jaeger,
Saurabh W. Jha,
Kaisey S. Mandel,
Justin D. R. Pierel,
Brodie Popovic,
Benjamin M. Rose,
David Rubin,
Benjamin J. Shappee,
Stephen Thorp,
John L. Tonry,
R. Brent Tully,
Maria Vincenzi
Abstract:
While the sample of optical Type Ia Supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good "standard candles," is still $\lesssim$ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from $0.002 < z < 0.09$ as a part of the Dark En…
▽ More
While the sample of optical Type Ia Supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good "standard candles," is still $\lesssim$ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from $0.002 < z < 0.09$ as a part of the Dark Energy, H$_0$, and peculiar Velocities using Infrared Light from Supernovae (DEHVILS) survey. Observations are taken using UKIRT's WFCAM, where the median depth of the images is 20.7, 20.1, and 19.3 mag (Vega) for $Y$, $J$, and $H$-bands, respectively. The median number of epochs per SN Ia is 18 for all three bands ($YJH$) combined and 6 for each band individually. We fit 47 SN Ia LCs that pass strict quality cuts using three LC models, SALT3, SNooPy, and BayeSN and find scatter on the Hubble diagram to be comparable to or better than scatter from optical-only fits in the literature. Fitting NIR-only LCs, we obtain standard deviations ranging from 0.128-0.135 mag. Additionally, we present a refined calibration method for transforming 2MASS magnitudes to WFCAM magnitudes using HST CALSPEC stars that results in a 0.03 mag shift in the WFCAM $Y$-band magnitudes.
△ Less
Submitted 10 April, 2023; v1 submitted 27 January, 2023;
originally announced January 2023.
-
A BayeSN Distance Ladder: $H_0$ from a consistent modelling of Type Ia supernovae from the optical to the near infrared
Authors:
Suhail Dhawan,
Stephen Thorp,
Kaisey S. Mandel,
Sam M. Ward,
Gautham Narayan,
Saurabh W. Jha,
Thaisen Chant
Abstract:
The local distance ladder estimate of the Hubble constant ($H_0$) is important in cosmology, given the recent tension with the early universe inference. We estimate $H_0$ from the Type Ia supernova (SN~Ia) distance ladder, inferring SN~Ia distances with the hierarchical Bayesian SED model, BayeSN. This method has a notable advantage of being able to continuously model the optical and near-infrared…
▽ More
The local distance ladder estimate of the Hubble constant ($H_0$) is important in cosmology, given the recent tension with the early universe inference. We estimate $H_0$ from the Type Ia supernova (SN~Ia) distance ladder, inferring SN~Ia distances with the hierarchical Bayesian SED model, BayeSN. This method has a notable advantage of being able to continuously model the optical and near-infrared (NIR) SN~Ia light curves simultaneously. We use two independent distance indicators, Cepheids or the tip of the red giant branch (TRGB), to calibrate a Hubble-flow sample of 67 SNe~Ia with optical and NIR data. We estimate $H_0 = 74.82 \pm 0.97$ (stat) $\pm\, 0.84$ (sys) km\,s$^{-1}$\,Mpc$^{-1}$ when using the calibration with Cepheid distances to 37 host galaxies of 41 SNe~Ia, and $70.92 \pm 1.14$ (stat) $\pm\,1.49$ (sys) km\,s$^{-1}$\,Mpc$^{-1}$ when using the calibration with TRGB distances to 15 host galaxies of 18 SNe~Ia. For both methods, we find a low intrinsic scatter $σ_{\rm int} \lesssim 0.1$ mag. We test various selection criteria and do not find significant shifts in the estimate of $H_0$. Simultaneous modelling of the optical and NIR yields up to $\sim$15\% reduction in $H_0$ uncertainty compared to the equivalent optical-only cases. With improvements expected in other rungs of the distance ladder, leveraging joint optical-NIR SN~Ia data can be critical to reducing the $H_0$ error budget.
△ Less
Submitted 2 August, 2023; v1 submitted 14 November, 2022;
originally announced November 2022.
-
The Young Supernova Experiment Data Release 1 (YSE DR1): Light Curves and Photometric Classification of 1975 Supernovae
Authors:
P. D. Aleo,
K. Malanchev,
S. Sharief,
D. O. Jones,
G. Narayan,
R. J. Foley,
V. A. Villar,
C. R. Angus,
V. F. Baldassare,
M. J. Bustamante-Rosell,
D. Chatterjee,
C. Cold,
D. A. Coulter,
K. W. Davis,
S. Dhawan,
M. R. Drout,
A. Engel,
K. D. French,
A. Gagliano,
C. Gall,
J. Hjorth,
M. E. Huber,
W. V. Jacobson-Galán,
C. D. Kilpatrick,
D. Langeroodi
, et al. (58 additional authors not shown)
Abstract:
We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from…
▽ More
We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from young and fast-rising supernovae (SNe) to transients that persist for over a year, with a redshift distribution reaching z~0.5. We present relative SN rates from YSE's magnitude- and volume-limited surveys, which are consistent with previously published values within estimated uncertainties for untargeted surveys. We combine YSE and ZTF data, and create multi-survey SN simulations to train the ParSNIP and SuperRAENN photometric classification algorithms; when validating our ParSNIP classifier on 472 spectroscopically classified YSE DR1 SNe, we achieve 82% accuracy across three SN classes (SNe Ia, II, Ib/Ic) and 90% accuracy across two SN classes (SNe Ia, core-collapse SNe). Our classifier performs particularly well on SNe Ia, with high (>90%) individual completeness and purity, which will help build an anchor photometric SNe Ia sample for cosmology. We then use our photometric classifier to characterize our photometric sample of 1483 SNe, labeling 1048 (~71%) SNe Ia, 339 (~23%) SNe II, and 96 (~6%) SNe Ib/Ic. YSE DR1 provides a training ground for building discovery, anomaly detection, and classification algorithms, performing cosmological analyses, understanding the nature of red and rare transients, exploring tidal disruption events and nuclear variability, and preparing for the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time.
△ Less
Submitted 21 February, 2023; v1 submitted 14 November, 2022;
originally announced November 2022.
-
Relative intrinsic scatter in hierarchical Type Ia supernova siblings analyses: Application to SNe 2021hpr, 1997bq & 2008fv in NGC 3147
Authors:
Sam M. Ward,
Stephen Thorp,
Kaisey S. Mandel,
Suhail Dhawan,
David O. Jones,
Kirsty Taggart,
Ryan J. Foley,
Gautham Narayan,
Kenneth C. Chambers,
David A. Coulter,
Kyle W. Davis,
Thomas de Boer,
Kaylee de Soto,
Nicholas Earl,
Alex Gagliano,
Hua Gao,
Jens Hjorth,
Mark E. Huber,
Luca Izzo,
Danial Langeroodi,
Eugene A. Magnier,
Peter McGill,
Armin Rest,
César Rojas-Bravo,
Radosław Wojtak
Abstract:
We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian fr…
▽ More
We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian framework for analyzing SN Ia siblings. At its core is the cosmology-independent relative intrinsic scatter parameter, $σ_{Rel}$: the dispersion of siblings distance estimates relative to one another within a galaxy. It quantifies the contribution towards the total intrinsic scatter, $σ_0$, from within-galaxy variations about the siblings' common properties. It also affects the combined-distance uncertainty. We present analytic formulae for computing a $σ_{Rel}$-posterior from individual siblings distances (estimated using any SN-model). Applying a newly trained BayeSN model, we fit the light curves of each sibling in NGC 3147 individually, to yield consistent distance estimates. However, the wide $σ_{Rel}$-posterior means $σ_{Rel}\approxσ_0$ is not ruled out. We thus combine the distances by marginalizing over $σ_{Rel}$ with an informative prior: $σ_{Rel}\sim U(0,σ_0)$. Simultaneously fitting the trio's light curves improves constraints on distance, and each sibling's individual dust parameters, compared to individual fits. Higher correlation also tightens dust parameter constraints. Therefore, $σ_{Rel}$-marginalization yields robust estimates of siblings distances for cosmology, and dust parameters for siblings-host correlation studies. Incorporating NGC 3147's Cepheid-distance yields $H_0=78.4\pm 6.5\,$km/s/Mpc. Our work motivates analyses of homogeneous siblings samples, to constrain $σ_{Rel}$, and its SN-model dependence.
△ Less
Submitted 1 September, 2023; v1 submitted 21 September, 2022;
originally announced September 2022.
-
Constraining the SN Ia Host Galaxy Dust Law Distribution and Mass Step: Hierarchical BayeSN Analysis of Optical and Near-Infrared Light Curves
Authors:
Stephen Thorp,
Kaisey S. Mandel
Abstract:
We use the BayeSN hierarchical probabilistic SED model to analyse the optical-NIR ($BVriYJH$) light curves of 86 Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project to investigate the SN Ia host galaxy dust law distribution and correlations between SN Ia Hubble residuals and host mass. Our Bayesian analysis simultaneously constrains the mass step and dust $R_V$ population distribution…
▽ More
We use the BayeSN hierarchical probabilistic SED model to analyse the optical-NIR ($BVriYJH$) light curves of 86 Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project to investigate the SN Ia host galaxy dust law distribution and correlations between SN Ia Hubble residuals and host mass. Our Bayesian analysis simultaneously constrains the mass step and dust $R_V$ population distribution by leveraging optical-NIR colour information. We demonstrate how a simplistic analysis where individual $R_V$ values are first estimated for each SN separately, and then the sample variance of these point estimates is computed, overestimates the $R_V$ population variance $σ_R^2$. This bias is exacerbated when neglecting residual intrinsic colour variation beyond that due to light curve shape. Instead, Bayesian shrinkage estimates of $σ_R$ are more accurate, with fully hierarchical analysis of the light curves being ideal. For the 75 SNe with low-to-moderate reddening (peak apparent $B-V\leq0.3$), we estimate an $R_V$ distribution with population mean $μ_R=2.59\pm0.14$, and standard deviation $σ_R=0.62\pm0.16$. Splitting this subsample at the median host galaxy mass ($10^{10.57}~\mathrm{M}_\odot$) yields consistent estimated $R_V$ distributions between low- and high-mass galaxies, with $μ_R=2.79\pm0.18$, $σ_R=0.42\pm0.24$, and $μ_R=2.35\pm0.27$, $σ_R=0.74\pm0.36$, respectively. When estimating distances from the full optical-NIR light curves while marginalising over various forms of the dust $R_V$ distribution, a mass step of $\gtrsim0.06$ mag persists in the Hubble residuals at the median host mass.
△ Less
Submitted 21 September, 2022;
originally announced September 2022.
-
A graph-based spectral classification of Type II supernovae
Authors:
Rafael S. de Souza,
Stephen Thorp,
Lluís Galbany,
Emille E. O. Ishida,
Santiago González-Gaitán,
Morgan A. Schmitz,
Alberto Krone-Martins,
Christina Peters
Abstract:
Given the ever-increasing number of time-domain astronomical surveys, employing robust, interpretative, and automated data-driven classification schemes is pivotal. Based on graph theory, we present new data-driven classification heuristics for spectral data. A spectral classification scheme of Type II supernovae (SNe II) is proposed based on the phase relative to the maximum light in the $V$ band…
▽ More
Given the ever-increasing number of time-domain astronomical surveys, employing robust, interpretative, and automated data-driven classification schemes is pivotal. Based on graph theory, we present new data-driven classification heuristics for spectral data. A spectral classification scheme of Type II supernovae (SNe II) is proposed based on the phase relative to the maximum light in the $V$ band and the end of the plateau phase. We utilize a compiled optical data set that comprises 145 SNe and 1595 optical spectra in 4000-9000 $\overset{\circ}{\mathrm {A}}$. Our classification method naturally identifies outliers and arranges the different SNe in terms of their major spectral features. We compare our approach to the off-the-shelf umap manifold learning and show that both strategies are consistent with a continuous variation of spectral types rather than discrete families. The automated classification naturally reflects the fast evolution of Type II SNe around the maximum light while showcasing their homogeneity close to the end of the plateau phase. The scheme we develop could be more widely applicable to unsupervised time series classification or characterisation of other functional data.
△ Less
Submitted 1 June, 2023; v1 submitted 28 June, 2022;
originally announced June 2022.
-
Cosmological Results from the RAISIN Survey: Using Type Ia Supernovae in the Near Infrared as a Novel Path to Measure the Dark Energy Equation of State
Authors:
D. O. Jones,
K. S. Mandel,
R. P. Kirshner,
S. Thorp,
P. M. Challis,
A. Avelino,
D. Brout,
C. Burns,
R. J. Foley,
Y. -C. Pan,
D. M. Scolnic,
M. R. Siebert,
R. Chornock,
W. L. Freedman,
A. Friedman,
J. Frieman,
L. Galbany,
E. Hsiao,
L. Kelsey,
G. H. Marion,
R. C. Nichol,
P. E. Nugent,
M. M. Phillips,
A. Rest,
A. G. Riess
, et al. (4 additional authors not shown)
Abstract:
Type Ia supernovae (SNe Ia) are more precise standardizable candles when measured in the near-infrared (NIR) than in the optical. With this motivation, from 2012-2017 we embarked on the RAISIN program with the Hubble Space Telescope (HST) to obtain rest-frame NIR light curves for a cosmologically distant sample of 37 SN Ia ($0.2 \lesssim z \lesssim 0.6$) discovered by Pan-STARRS and the Dark Energ…
▽ More
Type Ia supernovae (SNe Ia) are more precise standardizable candles when measured in the near-infrared (NIR) than in the optical. With this motivation, from 2012-2017 we embarked on the RAISIN program with the Hubble Space Telescope (HST) to obtain rest-frame NIR light curves for a cosmologically distant sample of 37 SN Ia ($0.2 \lesssim z \lesssim 0.6$) discovered by Pan-STARRS and the Dark Energy Survey. By comparing higher-$z$ HST data with 42 SN Ia at $z<0.1$ observed in the NIR by the Carnegie Supernova Project, we construct a Hubble diagram from NIR observations (with only time of maximum light and some selection cuts from optical data) to pursue a unique avenue to constrain the dark energy equation of state parameter, $w$. We analyze the dependence of the full set of Hubble residuals on the SN Ia host galaxy mass and find Hubble residual steps of size $\sim$0.06-0.1~mag with 1.5- to 2.5-$σ$ significance depending on the method and step location. Combining our NIR sample with CMB constraints, we find $1+w=-0.17\pm0.12$ (stat$+$syst). The largest systematic errors are the redshift-dependent SN selection biases and the properties of the NIR mass step. We also use these data to measure $H_0=75.9\pm 2.2$ km s$^{-1}$ Mpc$^{-1}$ from stars with geometric distance calibration in the hosts of 8 SNe Ia observed in the NIR versus $H_0=71.2\pm3.8$ km s$^{-1}$ Mpc$^{-1}$ using an inverse distance ladder approach tied to Planck. Using optical data we find $1+w=-0.10\pm0.09$ and with optical and NIR data combined, we find $1+w=-0.06\pm0.07$; these shifts of up to 0.11 in $w$ could point to inconsistency in optical versus NIR SN models. There will be many opportunities to improve this NIR measurement and better understand systematic uncertainties through larger low-$z$ samples, new light-curve models, calibration improvements, and by building high-$z$ samples from the Roman Space Telescope.
△ Less
Submitted 20 July, 2022; v1 submitted 19 January, 2022;
originally announced January 2022.
-
A Reference Survey for Supernova Cosmology with the Nancy Grace Roman Space Telescope
Authors:
B. M. Rose,
C. Baltay,
R. Hounsell,
P. Macias,
D. Rubin,
D. Scolnic,
G. Aldering,
R. Bohlin,
M. Dai,
S. E. Deustua,
R. J. Foley,
A. Fruchter,
L. Galbany,
S. W. Jha,
D. O. Jones,
B. A. Joshi,
P. L. Kelly,
R. Kessler,
R. P. Kirshner,
K. S. Mandel,
S. Perlmutter,
J. Pierel,
H. Qu,
D. Rabinowitz,
A. Rest
, et al. (11 additional authors not shown)
Abstract:
This note presents an initial survey design for the Nancy Grace Roman High-latitude Time Domain Survey. This is not meant to be a final or exhaustive list of all the survey strategy choices, but instead presents a viable path towards achieving the desired precision and accuracy of dark energy measurements using Type Ia supernovae (SNe Ia). We describe a survey strategy that use six filters (RZYJH…
▽ More
This note presents an initial survey design for the Nancy Grace Roman High-latitude Time Domain Survey. This is not meant to be a final or exhaustive list of all the survey strategy choices, but instead presents a viable path towards achieving the desired precision and accuracy of dark energy measurements using Type Ia supernovae (SNe Ia). We describe a survey strategy that use six filters (RZYJH and F) and the prism on the Roman Wide Field Instrument. This survey has two tiers, one "wide" which targets SNe Ia at redshifts up to 1 and one "deep" targeting redshifts up to 1.7; for each, four filters are used (with Y and J used in both tiers). We propose one field each in the north and south continuous viewing zones, and expect to obtain high-quality distances of $\sim$12,000 SNe Ia with $\sim$5,000 at z > 1. We propose a wide-tier area of $\sim$19 deg$^2$ and a deep tier of $\sim$5 deg$^2$. Exposure times range from 100 s to 900 s for imaging and 900 s to 3600 s for the prism. These exposure times would reach $\sim$25.5 mag and $\sim$26.5 mag for the wide and deep tiers respectively, with deep co-add stacks reaching $\sim$28 mag and $\sim$29 mag. The total survey spans two years, with a total allocation time of six months, and a cadence of $\sim$5 days.
△ Less
Submitted 4 November, 2021;
originally announced November 2021.
-
Testing the Consistency of Dust Laws in SN Ia Host Galaxies: A BayeSN Examination of Foundation DR1
Authors:
Stephen Thorp,
Kaisey S. Mandel,
David O. Jones,
Sam M. Ward,
Gautham Narayan
Abstract:
We apply BayeSN, our new hierarchical Bayesian model for the SEDs of Type Ia supernovae (SNe Ia), to analyse the $griz$ light curves of 157 nearby SNe Ia ($0.015<z<0.08$) from the public Foundation DR1 dataset. We train a new version of BayeSN, continuous from 0.35--0.95 $μ$m, which we use to model the properties of SNe Ia in the rest-frame $z$-band, study the properties of dust in their host gala…
▽ More
We apply BayeSN, our new hierarchical Bayesian model for the SEDs of Type Ia supernovae (SNe Ia), to analyse the $griz$ light curves of 157 nearby SNe Ia ($0.015<z<0.08$) from the public Foundation DR1 dataset. We train a new version of BayeSN, continuous from 0.35--0.95 $μ$m, which we use to model the properties of SNe Ia in the rest-frame $z$-band, study the properties of dust in their host galaxies, and construct a Hubble diagram of SN Ia distances determined from full $griz$ light curves. Our $griz$ Hubble diagram has a low total RMS of 0.13 mag using BayeSN, compared to 0.16 mag using SALT2. Additionally, we test the consistency of the dust law $R_V$ between low- and high-mass host galaxies by using our model to fit the full time- and wavelength-dependent SEDs of SNe Ia up to moderate reddening (peak apparent $B-V \lesssim 0.3$). Splitting the population at the median host mass, we find $R_V=2.84\pm0.31$ in low-mass hosts, and $R_V=2.58\pm0.23$ in high-mass hosts, both consistent with the global value of $R_V=2.61\pm0.21$ that we estimate for the full sample. For all choices of mass split we consider, $R_V$ is consistent across the step within $\lesssim1.2σ$. Modelling population distributions of dust laws in low- and high-mass hosts, we find that both subsamples are highly consistent with the full sample's population mean $μ(R_V) = 2.70\pm0.25$ with a 95% upper bound on the population $σ(R_V) < 0.61$. The $R_V$ population means are consistent within $\lesssim1.2σ$. We find that simultaneous fitting of host-mass-dependent dust properties within our hierarchical model does not account for the conventional mass step.
△ Less
Submitted 29 October, 2021; v1 submitted 10 February, 2021;
originally announced February 2021.
-
A Hierarchical Bayesian SED Model for Type Ia Supernovae in the Optical to Near-Infrared
Authors:
Kaisey S. Mandel,
Stephen Thorp,
Gautham Narayan,
Andrew S. Friedman,
Arturo Avelino
Abstract:
While conventional Type Ia supernova (SN Ia) cosmology analyses rely primarily on rest-frame optical light curves to determine distances, SNe Ia are excellent standard candles in near-infrared (NIR) light, which is significantly less sensitive to dust extinction. A SN Ia spectral energy distribution (SED) model capable of fitting rest-frame NIR observations is necessary to fully leverage current a…
▽ More
While conventional Type Ia supernova (SN Ia) cosmology analyses rely primarily on rest-frame optical light curves to determine distances, SNe Ia are excellent standard candles in near-infrared (NIR) light, which is significantly less sensitive to dust extinction. A SN Ia spectral energy distribution (SED) model capable of fitting rest-frame NIR observations is necessary to fully leverage current and future SN Ia datasets from ground- and space-based telescopes including HST, LSST, JWST, and RST. We construct a hierarchical Bayesian model for SN Ia SEDs, continuous over time and wavelength, from the optical to NIR ($B$ through $H$, or $0.35 -1.8\, μ$m). We model the SED as a combination of physically-distinct host galaxy dust and intrinsic spectral components. The distribution of intrinsic SEDs over time and wavelength is modelled with probabilistic functional principal components and the covariance of residual functions. We train the model on a nearby sample of 79 SNe Ia with joint optical and NIR light curves by sampling the global posterior distribution over dust and intrinsic latent variables, SED components, and population hyperparameters. The photometric distances of SNe Ia with NIR data near maximum light obtain a total RMS error of 0.10 mag with our BayeSN model, compared to 0.14 mag with SNooPy and SALT2 for the same sample. Jointly fitting the optical and NIR data of the full sample for a global host dust law, we find $R_V = 2.9 \pm 0.2$, consistent with the Milky Way average.
△ Less
Submitted 24 August, 2020; v1 submitted 17 August, 2020;
originally announced August 2020.
-
Tidal disruption events from massive black hole binaries: predictions for ongoing and future surveys
Authors:
Stephen Thorp,
Eli Chadwick,
Alberto Sesana
Abstract:
We compute the expected cosmic rates of tidal disruption events induced by individual massive black holes (MBHs) and by MBH binaries (MBHBs) - with a specific focus on the latter class - to explore the potential of TDEs to probe the cosmic population of sub-pc MBHBs. Rates are computed by combining MBH and MBHB population models derived from large cosmological simulations with estimates of the ind…
▽ More
We compute the expected cosmic rates of tidal disruption events induced by individual massive black holes (MBHs) and by MBH binaries (MBHBs) - with a specific focus on the latter class - to explore the potential of TDEs to probe the cosmic population of sub-pc MBHBs. Rates are computed by combining MBH and MBHB population models derived from large cosmological simulations with estimates of the induced TDE rates for each class of objects. We construct empirical TDE spectra that fit a large number of observations in the optical, UV and X-ray and consider their observability by current and future survey instruments. Consistent with results in the literature, and depending on the detailed assumption of the model, we find that LSST and Gaia in optical and eROSITA in X-ray will observe a total of 3000-6000, 80-180 and 600-900 TDEs per year, respectively. Depending on the survey, one to several percent of these are prompted by MBHBs. In particular both LSST and eROSITA are expected to see 150-450 MBHB induced TDEs in their respective mission lifetimes, including 5-100 repeated flares. The latter provide an observational sample of binary candidates with relatively low contamination and have the potential of unveiling the sub-pc population of MBHBs in the mass range $10^5M_\odot<M<10^7M_\odot$, thus informing future low frequency gravitational wave observatories.
△ Less
Submitted 15 July, 2019; v1 submitted 5 November, 2018;
originally announced November 2018.
-
Transient Stability Assessment of Cascade Tripping of Renewable Sources Using SOS
Authors:
Chetan Mishra,
James S. Thorp,
Anamitra Pal,
Virgilio A. Centeno
Abstract:
There has been significant increase in penetration of renewable generation (RG) sources all over the world. Localized concentration of many such generators could initiate a cascade tripping sequence that might threaten the stability of the entire system. Understanding the impact of cascade tripping process would help the system planner identify trip sequences that must be blocked in order to incre…
▽ More
There has been significant increase in penetration of renewable generation (RG) sources all over the world. Localized concentration of many such generators could initiate a cascade tripping sequence that might threaten the stability of the entire system. Understanding the impact of cascade tripping process would help the system planner identify trip sequences that must be blocked in order to increase stability. In this work, we attempt to understand the consequences of cascade tripping mechanism through a Lyapunov approach. A conservative definition for the stability region (SR) along with its estimation for a given cascading sequence using sum of squares (SOS) programming is proposed. Finally, a simple probabilistic definition of the SR is used to visualize the risk of instability and understand the impact of blocking trip sequences. A 3-machine system with significant RG penetration is used to demonstrate the idea.
△ Less
Submitted 12 February, 2018;
originally announced February 2018.
-
Estimating Relevant Portion of Stability Region using Lyapunov Approach and Sum of Squares
Authors:
Chetan Mishra,
James S. Thorp,
Virgilio A. Centeno,
Anamitra Pal
Abstract:
Traditional Lyapunov based transient stability assessment approaches focus on identifying the stability region (SR) of the equilibrium point under study. When trying to estimate this region using Lyapunov functions, the shape of the final estimate is often limited by the degree of the function chosen, a limitation that results in conservativeness in the estimate of the SR. More conservative the es…
▽ More
Traditional Lyapunov based transient stability assessment approaches focus on identifying the stability region (SR) of the equilibrium point under study. When trying to estimate this region using Lyapunov functions, the shape of the final estimate is often limited by the degree of the function chosen, a limitation that results in conservativeness in the estimate of the SR. More conservative the estimate is in a particular region of state space, smaller is the estimate of the critical clearing time for disturbances that drive the system towards that region. In order to reduce this conservativeness, we propose a methodology that uses the disturbance trajectory data to skew the shape of the final Lyapunov based SR estimate. We exploit the advances made in the theory of sum of squares decomposition to algorithmically estimate this region. The effectiveness of this technique is demonstrated on a power systems classical model.
△ Less
Submitted 20 October, 2017;
originally announced October 2017.
-
Stability Region Estimation Under Low Voltage Ride Through Constraints using Sum of Squares
Authors:
Chetan Mishra,
James S. Thorp,
Virgilio A. Centeno,
Anamitra Pal
Abstract:
The increasing penetration of inverter based renewable generation (RG) in the form of solar photo-voltaic (PV) or wind has introduced numerous operational challenges and uncertainties. According to the standards, these generators are made to trip offline if their operating requirements are not met. In an RG-rich system, this might alter the system dynamics and/or cause shifting of the equilibrium…
▽ More
The increasing penetration of inverter based renewable generation (RG) in the form of solar photo-voltaic (PV) or wind has introduced numerous operational challenges and uncertainties. According to the standards, these generators are made to trip offline if their operating requirements are not met. In an RG-rich system, this might alter the system dynamics and/or cause shifting of the equilibrium points to the extent that a cascaded tripping scenario is manifested. The present work attempts at avoiding such scenarios by estimating the constrained stability region (CSR) inside which the system must operate using maximal level set of a Lyapunov function estimated through sum of squares (SOS) technique. A time-independent conservative approximation of the LVRT constraint is initially derived for a classical model of the power system. The proposed approach is eventually validated by evaluating the stability of a 3 machine test system with trip-able RG.
△ Less
Submitted 28 June, 2017;
originally announced June 2017.
