Outside the Wall: Hydrodynamics of Type I Supernovae Interacting with a Partially Swept-Up Circumstellar Medium

Explaining the observed diversity of supernovae (SNe) and the physics of explosion requires knowledge of their progenitor stars, which can be obtained by constraining the circumstellar medium (CSM). Models of the SN ejecta colliding with CSM are necessary to infer the structure of the CSM and tie it back to a progenitor model. Recent SNe I revealed CSM concentrated at a distance $r\sim10^16$ cm, for which models of SN interaction are extremely limited. In this paper, we assume the concentrated region is a "wall" representing swept-up material, and unswept material lies outside the wall. We simulate one-dimensional hydrodynamics of SNe Ia & Ib impacting 300 unique CSM configurations using RT1D, which captures the Rayleigh-Taylor instability. We find that the density ratio between the wall and ejecta -- denoted $A_0$ or "wall height" -- is key, and higher walls deviate more from self-similar evolution. Functional fits accounting for $A_0$ are presented for the forward shock radius evolution. We show that higher walls have more degeneracy between CSM properties in the deceleration parameter, slower shocks, deeper-probing reverse shocks, slower shocked ejecta, less ejecta mass than CSM in the shock, and more mixing of ejecta into the CSM at early times. We analyze observations of SN 2014C (Type Ib) and suggest that it had a moderately high wall ($10 < A_0 < 200$) and wind-like outer CSM. We also postulate an alternate interpretation for the radio data of SN 2014C, that the radio rise occurs in the wind rather than the wall. Finally, we find that hydrodynamic measurements at very late times cannot distinguish the presence of a wall, except perhaps as an anomalously wide shock region.Comment: 17 pages, 13 figures, accepted to Ap

Skill set profile clustering based on student capability vectors computed from online tutoring data

In educational research, a fundamental goal is identifying which skills students have mastered, which skills they have not, and which skills they are in the process of mastering. As the number of examinees, items, and skills increases, the estimation of even simple cognitive diagnosis models becomes difficult. To address this, we introduce a capability matrix showing for each skill the proportion correct on all items tried by each student involving that skill. We apply variations of common clustering methods to this matrix and discuss conditioning on sparse subspaces. We demonstrate the feasibility and scalability of our method on several simulated datasets and illustrate the difficulties inherent in real data using a subset of online mathematics tutor data. We also comment on the interpretability and application of the results for teachers

Skill set profile clustering: the empty K-means algorithm with automatic specification of starting cluster centers

While students’ skill set profiles can be estimated with formal cognitive diagnosis models [8], their computational complexity makes simpler proxy skill estimates attractive [1, 4, 6]. These estimates can be clustered to generate groups of similar students. Often hierarchical agglomerative clustering or k-means clustering is utilized, requiring, for K skills, the specification of 2^K clusters. The number of skill set profiles/clusters can quickly become computationally intractable. Moreover, not all profiles may be present in the population. We present a flexible version of k-means that allows for empty clusters. We also specify a method to determine efficient starting centers based on the Q-matrix. Combining the two substantially improves the clustering results and allows for analysis of data sets previously thought impossible

Intermediate Palomar Transient Factory: Realtime Image Subtraction Pipeline

A fast-turnaround pipeline for realtime data reduction plays an essential role in discovering and permitting follow-up observations to young supernovae and fast-evolving transients in modern time-domain surveys. In this paper, we present the realtime image subtraction pipeline in the intermediate Palomar Transient Factory. By using high-performance computing, efficient database, and machine learning algorithms, this pipeline manages to reliably deliver transient candidates within ten minutes of images being taken. Our experience in using high performance computing resources to process big data in astronomy serves as a trailblazer to dealing with data from large-scale time-domain facilities in near future.Comment: 18 pages, 6 figures, accepted for publication in PAS

Rates and Properties of Strongly Gravitationally Lensed Supernovae and their Host Galaxies in Time-Domain Imaging Surveys

Supernovae that are strongly gravitationally lensed (gLSNe) by galaxies are powerful probes of astrophysics and cosmology that will be discovered systematically by next-generation wide-field, high-cadence imaging surveys such as the Zwicky Transient Facility (ZTF) and the Large Synoptic Survey Telescope (LSST). Here we use pixel-level simulations that include dust, observing strategy, and multiple supernova subtypes to forecast the rates and properties of gLSNe that ZTF and LSST will find. Applying the resolution-insensitive discovery strategy of Goldstein et al. (2018), we forecast that ZTF (LSST) can discover 0.02 (0.79) 91bg-like, 0.17 (5.92) 91T-like, 1.22 (47.84) Type Ia, 2.76 (88.51) Type IIP, 0.31 (12.78) Type IIL, and 0.36 (15.43) Type Ib/c gLSNe per year. We also forecast that the surveys can discover at least 3.75 (209.32) Type IIn gLSNe per year, for a total of at least 8.60 (380.60) gLSNe per year under fiducial observing strategies. ZTF gLSNe have a median $z_s=0.9$, $z_l=0.35$, $\mu_\mathrm{tot}=30$, $\Delta t_\mathrm{max}= 10$ days, $\min(\theta)= 0.25^{\prime\prime}$, and $N_\mathrm{img} = 4$. LSST gLSNe are less compact and less magnified, with a median $z_s=1.0$, $z_l=0.4$, $\mu_\mathrm{tot}\approx6$, $\Delta t_\mathrm{max} = 25$ days, $\min(\theta)=0.6^{\prime\prime}$, and $N_\mathrm{img} = 2$. As the properties of lensed host galaxy arcs provide critical information for lens mass modeling, we develop a model of the supernova--host galaxy connection and use it to simulate realistic images of the supernova--host--lens systems. We find that the vast majority of gLSN host galaxies will be multiply imaged, enabling detailed constraints on lens models with sufficiently deep high-resolution imaging taken after the supernova has faded. We release the results of our simulations to the public as catalogs at this URL: http://portal.nersc.gov/project/astro250/glsne/.Comment: 57 pages, 66 equations, 36 figures, 4 tables, Submitted to ApJS, comments welcome, v2 replaced some figures with rasterized versions to reduce load on PDF viewer

The volumetric rate of calcium-rich transients in the local universe

We present a measurement of the volumetric rate of `calcium-rich' optical transients in the local universe, using a sample of three events from the Palomar Transient Factory (PTF). This measurement builds on a detailed study of the PTF transient detection efficiencies, and uses a Monte Carlo simulation of the PTF survey. We measure the volumetric rate of calcium-rich transients to be higher than previous estimates: $1.21^{+1.13}_{-0.39}\times10^{-5}$ events yr$^{-1}$ Mpc$^{-3}$. This is equivalent to 33-94% of the local volumetric type Ia supernova rate. This calcium-rich transient rate is sufficient to reproduce the observed calcium abundances in galaxy clusters, assuming an asymptotic calcium yield per calcium-rich event of ~0.05$\mathrm{M}_{\odot}$. We also study the PTF detection efficiency of these transients as a function of position within their candidate host galaxies. We confirm as a real physical effect previous results that suggest calcium-rich transients prefer large physical offsets from their host galaxies.Comment: Accepted for publication in ApJ. 9 pages, 5 figure

Against the Wind: Radio Light Curves of Type Ia Supernovae Interacting with Low-Density Circumstellar Shells

For decades, a wide variety of observations spanning the radio through optical and on to the x-ray have attempted to uncover signs of type Ia supernovae (SNe Ia) interacting with a circumstellar medium (CSM). The goal of these studies is to constrain the nature of the hypothesized SN Ia mass-donor companion. A continuous CSM is typically assumed when interpreting observations of interaction. However, while such models have been successfully applied to core-collapse SNe, the assumption of continuity may not be accurate for SNe Ia, as shells of CSM could be formed by pre-supernova eruptions (novae). In this work, we model the interaction of SNe with a spherical, low density, finite-extent CSM and create a suite of synthetic radio synchrotron light curves. We find that CSM shells produce sharply peaked light curves, and identify a fiducial set of models that all obey a common evolution and can be used to generate radio light curves for interaction with an arbitrary shell. The relations obeyed by the fiducial models can be used to deduce CSM properties from radio observations; we demonstrate this by applying them to the non-detections of SN 2011fe and SN 2014J. Finally, we explore a multiple shell CSM configuration and describe its more complicated dynamics and resultant radio light curves.Comment: 15 pages, 11 figures, ApJ accepte

PTF11kx: A Type Ia Supernova with Hydrogen Emission Persisting After 3.5 Years

The optical transient PTF11kx exhibited both the characteristic spectral features of Type Ia supernovae (SNe Ia) and the signature of ejecta interacting with circumstellar material (CSM) containing hydrogen, indicating the presence of a nondegenerate companion. We present an optical spectrum at $1342$ days after peak from Keck Observatory, in which the broad component of H$\alpha$ emission persists with a similar profile as in early-time observations. We also present $Spitzer$ IRAC detections obtained $1237$ and $1818$ days after peak, and an upper limit from $HST$ ultraviolet imaging at $2133$ days. We interpret our late-time observations in context with published results - and reinterpret the early-time observations - in order to constrain the CSM's physical parameters and compare to theoretical predictions for recurrent nova systems. We find that the CSM's radial extent may be several times the distance between the star and the CSM's inner edge, and that the CSM column density may be two orders of magnitude lower than previous estimates. We show that the H$\alpha$ luminosity decline is similar to other SNe with CSM interaction, and demonstrate how our infrared photometry is evidence for newly formed, collisionally heated dust. We create a model for PTF11kx's late-time CSM interaction and find that X-ray reprocessing by photoionization and recombination cannot reproduce the observed H$\alpha$ luminosity, suggesting that the X-rays are thermalized and that H$\alpha$ radiates from collisional excitation. Finally, we discuss the implications of our results regarding the progenitor scenario and the geometric properties of the CSM for the PTF11kx system.Comment: 15 pages, 8 figures, 3 tables; submitted to Ap