Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). At $z=0.0866$, ZTF19aawfbtg (SN2019hge) has more than 10 optical spectra at phases from $-41$ to $+103$ days. Most of these spectra match well with that of PTF10hgi, suggesting SN2019hge as the second case of a He-rich SLSN-I. Confirmation comes from its near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features, yielding a lower limit to the He-rich fraction of $\sim9$\%. It is suggested that the excitation of He I atoms requires non-thermal radiation. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. However, how the magnetar energy is transported to the outer layers is unclear. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities, around $-20.06\pm0.16$ mag ($g$-band, rest-frame), in comparison to the full range of $-20$ to $-23$ mag covered by SLSN-I.

We report the discovery and panchromatic followup observations of the young Type Ic supernova, SN2020oi, in M100, a grand design spiral galaxy at a mere distance of $14$ Mpc. We followed up with observations at radio, X-ray and optical wavelengths from only a few days to several months after explosion. The optical behaviour of the supernova is similar to those of other normal Type Ic supernovae. The event was not detected in the X-ray band but our radio observation revealed a bright mJy source ($L_{\nu} \approx 1.2 \times 10^{27} {\rm erg\,s}^{-1} {\rm Hz}^{-1}$). Given, the relatively small number of stripped envelope SNe for which radio emission is detectable, we used this opportunity to perform a detailed analysis of the comprehensive radio dataset we obtained. The radio emitting electrons initially experience a phase of inverse Compton cooling which leads to steepening of the spectral index of the radio emission. Our analysis of the cooling frequency points to a large deviation from equipartition at the level of $\epsilon_e/\epsilon_B \gtrsim 200$, similar to a few other cases of stripped envelope SNe. Our modeling of the radio data suggests that the shockwave driven by the SN ejecta into the circumstellar matter (CSM) is moving at $\sim 3\times 10^{4}\,{\rm km\,s}^{-1}$. Assuming a constant mass-loss from the stellar progenitor, we find that the mass-loss rate is $\dot{M} \approx 1.4\times 10^{-4}\,{M}_{\odot}\,{\rm yr}^{-1}$, for an assumed wind velocity of $1000\,{\rm km\,s}^{-1}$. The temporal evolution of the radio emission suggests a radial CSM density structure steeper than the standard $r^{-2}$.

The true nature of type Ia supernovae (SNe Ia) remains largely unknown, and as survey statistics increase, the question of astrophysical systematic uncertainties rises, notably that of the SN Ia population evolution. In this paper, we study the dependence with redshift of the SN Ia SALT2.4 lightcurve stretch, a purely intrinsic SN property, to probe its potential redshift drift. The SN stretch has been shown to strongly correlate with the SN environment, notably with stellar age tracers. We model the underlying stretch distribution as a function of redshift, using the evolution of fraction of young and old SNe Ia as predicted by Rigault et al. (2018), and assuming constant underlying stretch distribution for each age population made of Gaussian mixtures. We test our prediction against published samples chosen to have negligible magnitude selection effects, so that any observed change is indeed of astrophysical and not observational origin. We clearly demonstrate that the underlying SN Ia stretch distribution is evolving as a function of redshift, and that the young/old drifting model is a much better description of the data than any time-constant model, including the sample-based asymmetric distributions usually used to correct Malmquist bias. The favored underlying stretch model is the bimodal one derived from Rigault et al. (2018): a high-stretch mode shared by both young and old environments, and a low-stretch mode exclusive to old environments. The precise impact of the redshift evolution of the SN Ia population intrinsic properties on cosmology remains to be studied. Yet,the astrophysical drift of the SN stretch distribution does affect current Malmquist bias corrections and thereby distances derived from SN affected by selection effects. We highlight that such a bias will increase with surveys covering increasingly larger redshift ranges, which is particularly important for LSST.

Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently associated one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the association of a radio-emitting tidal disruption event (AT2019dsg) with another high-energy neutrino, identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility (ZTF). The probability of finding any radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. The association suggests that tidal disruption events contribute to the cosmic neutrino flux. Unlike previous work which considered the rare subset of tidal disruption events with relativistic jets, our observations of AT2019dsg suggest an empirical model with a mildly-relativistic outflow.

The Nearby Supernova Factory has made spectrophotometric observations of Type Ia supernovae since $2004$. This work presents an interim version of the data produced, including $210$ supernovae observed between $2004$ and $2013$.

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey.

We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% on the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modelling of these events.

The Zwicky Transient Facility (ZTF) is performing a three-day cadence survey of the visible Northern sky (~3$\pi$). The transient candidates found in this survey are announced via public alerts. As a supplementary product ZTF is also conducting a large spectroscopic campaign: the ZTF Bright Transient Survey (BTS). The goal of the BTS is to spectroscopically classify all extragalactic transients brighter than 18.5 mag at peak brightness and immediately announce those classifications to the public. Extragalactic discoveries from ZTF are predominantly Supernovae (SNe). The BTS is the largest flux-limited SN survey to date. Here we present a catalog of the761 SNe that were classified during the first nine months of the survey (2018 Apr. 1 to 2018 Dec. 31). The BTS SN catalog contains redshifts based on SN template matching and spectroscopic host galaxy redshifts when available. Based on this data we perform an analysis of the redshift completeness of local galaxy catalogs, dubbed as the Redshift Completeness Fraction (RCF; the number of SN host galaxies with known spectroscopic redshift prior to SN discovery divided by the total number of SN hosts). In total, we identify the host galaxies of 512 Type Ia supernovae, 227 of which have known spectroscopic redshifts, yielding an RCF estimate of $44\% \pm1\%$. We find a steady decrease in the RCF with increasing distance in the local universe. For z<0.05, or ~200 Mpc, we find RCF=0.6, which has important ramifications when searching for multimessenger astronomical events. Prospects for dramatically increasing the RCF are limited to new multi-fiber spectroscopic instruments, or wide-field narrowband surveys. We find that existing galaxy redshift catalogs are only $50\%$ complete at $r\approx16.9$ mag. Pushing this limit several magnitudes deeper will pay huge dividends when searching for electromagnetic counterparts to gravitational wave events.

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain their progenitor properties. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point spread function (PSF) photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night ($3g+3r$), all of the 127 SNe Ia are detected in both $g$ and $r$ band more than 10\,d (in the rest frame) prior to the epoch of $g$-band maximum light. The redshifts of these objects range from $z=0.0181$ to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14\,d prior to maximum light (in the rest frame), with a subset of 9 objects being detected more than 17\,d before $g$-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample, including one 02cx-like event ZTF18abclfee (SN\,2018crl), one Ia-CSM SN ZTF18aaykjei (SN\,2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN\,2018eul), ZTF18abdpvnd (SN\,2018dvf), ZTF18aawpcel (SN\,2018cir) and ZTF18abddmrf (SN\,2018dsx).

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

When planning a survey for astronomical transients, many factors such as cadence, filter choice, sky coverage, and depth of observations need to be balanced in order to optimize the scientific gain of the survey. Here we present a software package called $\texttt{simsurvey}$ for simulating the supernova lightcurves that are expected based on a survey strategy, which can then be used to determine the potential for discoveries of each strategy in question. The code is set up in a modular fashion that allows easy modification of small details of the survey and enables the user to adapt it to any survey design and transient template that they wish to use in planning their survey. As an example of its utility, we use $\texttt{simsurvey}$ to simulate the lightcurve of several types of supernovae that the recently started Zwicky Transient Facility (ZTF) is expected to find and compare the results to the discoveries made during its early operations. We conclude that ZTF will find thousands of bright supernovae per year, of which about 10 could potentially be found with two days of explosion. Over the course of three years the survey will obtain lightcurves of about 1800 type Ia supernovae with $z < 0.1$ that can be used as distance indicators in cosmology if they are spectroscopically classified using additional telescopes. In a comparison to detections from the ZTF public survey, we found good agreement with the numbers of detections expected from the simulations.

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model is constructed from SNe Ia spectral properties and spectrophotometric data from The Nearby Supernova Factory collaboration. In a first step, a PCA-like method is used on spectral features measured at maximum light, which allows us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties are used to extract the average extinction curve. Third, an interpolation using Gaussian Processes facilitates using data taken at different epochs during the lifetime of a SN Ia and then projecting the data on a fixed time grid. Finally, the three steps are combined to build the SED model as a function of time and wavelength. This is the SUGAR model. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity, the second is correlated with their calcium lines. The addition of these parameters, as well as the high quality the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like ZTF, LSST or WFIRST.

As a candidate 'super-Chandrasekhar' or 09dc-like Type Ia supernova (SN Ia), SN 2012dn shares many characteristics with other members of this remarkable class of objects but lacks their extraordinary luminosity. Here, we present and discuss the most comprehensive optical data set of this SN to date, comprised of a densely sampled series of early-time spectra obtained within the Nearby Supernova Factory project, plus photometry and spectroscopy obtained at the VLT about 1 yr after the explosion. The light curves, colour curves, spectral time series and ejecta velocities of SN 2012dn are compared with those of other 09dc-like and normal SNe Ia, the overall variety within the class of 09dc-like SNe Ia is discussed, and new criteria for 09dc-likeness are proposed. Particular attention is directed to additional insight that the late-phase data provide. The nebular spectra show forbidden lines of oxygen and calcium, elements that are usually not seen in late-time spectra of SNe Ia, while the ionisation state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterised by an enhanced fading starting ~60 d after maximum and very low luminosities in the nebular phase, which is most readily explained by unusually early formation of clumpy dust in the ejecta. Taken together, these effects suggest a strongly perturbed ejecta density profile, which might lend support to the idea that 09dc-like characteristics arise from a brief episode of interaction with a hydrogen-deficient envelope during the first hours or days after the explosion.

High cadence transient surveys are able to capture supernovae closer to their first light than before. Applying analytical models to such early emission, we can constrain the progenitor stars properties. In this paper, we present observations of SN2018fif (ZTF18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN2018fif was surrounded by relatively small amounts of circumstellar material (CSM) compared to all previous cases. This particularity, coupled with the high cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman (2017). Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model. We find that the progenitor of SN2018fif was a large red supergiant, with a radius of R=1174_{-81}^{+208} solar radii, and an ejected mass of M=5.6_{-1.0}^{+9.1} solar masses. Our model also gives information on the explosion epoch, the progenitor inner structure, the shock velocity and the extinction. The large radius differs from previously modeled objects, and the difference could be either intrinsic or due to the relatively small amount of CSM around SN2018fif, perhaps making it a "cleaner" candidate for applying shock-cooling analytical models.

Both multi-messenger astronomy and new high-throughput wide-field surveys require flexible tools for the selection and analysis of astrophysical transients. We here introduce the Alert Management, Photometry and Evaluation of Lightcurves (AMPEL) system, an analysis framework designed for high-throughput surveys and suited for streamed data. AMPEL combines the functionality of an alert broker with a generic framework capable of hosting user-contributed code, that encourages provenance and keeps track of the varying information states that a transient displays. The latter concept includes information gathered over time and data policies such as access or calibration levels. We describe a novel ongoing real-time multi-messenger analysis using AMPEL to combine IceCube neutrino data with the alert streams of the Zwicky Transient Facility (ZTF). We also reprocess the first four months of ZTF public alerts, and compare the yields of more than 200 different transient selection functions to quantify efficiencies for selecting Type Ia supernovae that were reported to the Transient Name Server (TNS). We highlight three channels suitable for (1) the collection of a complete sample of extragalactic transients, (2) immediate follow-up of nearby transients and (3) follow-up campaigns targeting young, extragalactic transients. We confirm ZTF completeness in that all TNS supernovae positioned on active CCD regions were detected. AMPEL can assist in filtering transients in real time, running alert reaction simulations, the reprocessing of full datasets as well as in the final scientific analysis of transient data. This text introduces how users can design their own channels for inclusion in the AMPEL live instance that parses the ZTF stream and the real-time submission of high quality extragalactic supernova candidates to the TNS.

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

Current time domain facilities are discovering hundreds of new galactic and extra-galactic transients every week. Classifying the ever-increasing number of transients is challenging, yet crucial to further our understanding of their nature, discover new classes, or ensuring sample purity, for instance, for Supernova Ia cosmology. The Zwicky Transient Facility is one example of such a survey. In addition, it has a dedicated very-low resolution spectrograph, the SEDMachine, operating on the Palomar 60-inch telescope. This spectrograph's primary aim is object classification. In practice most, if not all, transients of interest brighter than ~19 mag are typed. This corresponds to approximately 10 to 15 targets a night. In this paper, we present a fully automated pipeline for the SEDMachine. This pipeline has been designed to be fast, robust, stable and extremely flexible. pysedm enables the fully automated spectral extraction of a targeted point source object in less than 5 minutes after the end of the exposure. The spectral color calibration is accurate at the few percent level. In the 19 weeks since pysedm entered production in early August of 2018, we have classified, among other objects, about 400 Type Ia supernovae and 140 Type II supernovae. We conclude that low resolution, fully automated spectrographs such as the `SEDMachine with pysedm' installed on 2-m class telescopes within the southern hemisphere could allow us to automatically and simultaneously type and obtain a redshift for most (if not all) bright transients detected by LSST within z<0.2, notably potentially all Type Ia Supernovae. In comparison to the current SEDM design, this would require higher spectral resolution (R~1000) and slightly improved throughput. With this perspective in mind, pysedm has been designed to easily be adaptable to any IFU-like spectrograph (see https://github.com/MickaelRigault/pysedm).

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

Type Ia supernova cosmology depends on the ability to fit and standardize observations of supernova magnitudes with an empirical model. We present here a series of new models of Type Ia Supernova spectral time series that capture a greater amount of supernova diversity than possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (\textsc{SNEMO}\footnote{https://snfactory.lbl.gov/snemo}). The models are constructed using spectrophotometric time series from $172$ individual supernovae from the Nearby Supernova Factory, comprising more than $2000$ spectra. Using the available observations, Gaussian Processes are used to predict a full spectral time series for each supernova. A matrix is constructed from the spectral time series of all the supernovae, and Expectation Maximization Factor Analysis is used to calculate the principal components of the data. K-fold cross-validation then determines the selection of model parameters and accounts for color variation in the data. Based on this process, the final models are trained on supernovae that have been dereddened using the Fitzpatrick and Massa extinction relation. Three final models are presented here: \textsc{SNEMO2}, a two-component model for comparison with current Type~Ia models; \textsc{SNEMO7}, a seven component model chosen for standardizing supernova magnitudes which results in a total dispersion of $0.100$~mag for a validation set of supernovae, of which $0.087$~mag is unexplained (a total dispersion of $0.113$~mag with unexplained dispersion of $0.097$~mag is found for the total set of training and validation supernovae); and \textsc{SNEMO15}, a comprehensive $15$ component model that maximizes the amount of spectral time series behavior captured.

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SNIa) for cosmology, we have statistically classified a large sample of nearby SNeIa into those located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNeIa in predominantly younger environments are DY=0.163\pm0.029 mag (5.7 sigma) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNeIa in hosts with lower or higher total stellar masses is smaller at DM=0.119\pm0.032 mag (4.5 sigma), for the same set of SNeIa. When fit simultaneously, the environment age offset remains very significant, with DY=0.129\pm0.032 mag (4.0 sigma), while the global stellar mass step is reduced to DM=0.064\pm0.029 mag (2.2 sigma). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Standardization using only the SNeIa in younger environments reduces the total dispersion from 0.142\pm0.008 mag to 0.120\pm0.010 mag. We show that as environment ages evolve with redshift a strong bias on measurement of the dark energy equation of state parameters can develop. Fortunately, data to measure and correct for this effect is likely to be available for many next-generation experiments. [abstract shorten]

We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of 2.8+0.6-0.5 (1.10 +- 0.23 mag)---compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from LambdaCDM simulations---making it the most amplified SN Ia discovered behind a galaxy cluster.

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. To perform such measurements the luminosity and cosmological redshift ($z$) of the SNe Ia have to be determined. The uncertainty on $z$ includes an unknown peculiar velocity, which can be very large for SNe Ia in the virialized cores of massive clusters. We determine which SNe Ia exploded in galaxy clusters. We then study how the correction for peculiar velocities of host galaxies inside the clusters improves the Hubble residuals. Using 145 SNe Ia from the Nearby Supernova Factory we found 11 candidates for membership in clusters. To estimate the redshift of a cluster we applied the bi-weight technique. Then, we use the galaxy cluster redshift instead of the host galaxy redshift to construct the Hubble diagram. For SNe Ia inside galaxy clusters the dispersion around the Hubble diagram when peculiar velocities are taken into account is smaller in comparison with a case without peculiar velocity correction, with a $wRMS=0.130\pm0.038$ mag instead of $wRMS=0.137\pm0.036$ mag. The significance of this improvement is 3.58 $\sigma$. If we remove the very nearby Virgo cluster member SN2006X ($z<0.01$) from the analysis, the significance decreases to 1.34 $\sigma$. The peculiar velocity correction is found to be highest for the SNe Ia hosted by blue spiral galaxies, with high local specific star formation rate and smaller stellar mass, seemingly counter to what might be expected given the heavy concentration of old, massive elliptical galaxies in clusters. As expected, the Hubble residuals of SNe Ia associated with massive galaxy clusters improve when the cluster redshift is taken as the cosmological redshift of the SN. This fact has to be taken into account in future cosmological analyses in order to achieve higher accuracy for cosmological redshift measurements. Here we provide an approach to do so.

Context. Observations of Type Ia supernovae (SNe Ia) can be used to derive accurate cosmological distances through empirical standardization techniques. Despite this success neither the progenitors of SNe Ia nor the explosion process are fully understood. The U-band region has been less well observed for nearby SNe, due to technical challenges, but is the most readily accessible band for high-redshift SNe. Aims. Using spectrophotometry from the Nearby Supernova Factory, we study the origin and extent of U-band spectroscopic variations in SNe Ia and explore consequences for their standardization and the potential for providing new insights into the explosion process. Methods. We divide the U-band spectrum into four wavelength regions {\lambda}(uNi), {\lambda}(uTi), {\lambda}(uSi) and {\lambda}(uCa). Two of these span the Ca H&K {\lambda}{\lambda} 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. (1) The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 $\pm$ 0.011 mag RMS. A traditional SALT2.4 fit on the same sample yields a 0.135 mag RMS. Standardization using uTi also reduces the difference in corrected magnitude between SNe originating from different host galaxy environments. (2) Early U-band spectra can be used to probe the Ni+Co distribution in the ejecta, thus offering a rare window into the source of lightcurve power. (3) The uCa flux further improves standardization, yielding a 0.086 $\pm$ 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach {\chi}$^2$/dof $\sim$ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for Type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova Integral Field Spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. The SCALA calibration unit is now fully deployed at the UH\,88 telescope, and with it color-calibration between 4000 {\AA} and 9000 {\AA} is stable at the percent level over a one-year baseline.

Correction of Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning 3,300 - 9,200 {\AA}, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E(B-V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (RV , RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands, and compare the 5780 {\AA} band with the dust-to-band ratio for the Milky Way. Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the Milky Way. Furthermore we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al. (1989), O'Donnell (1994), and Fitzpatrick (1999), and find the predictions of Fitzpatrick (1999) fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6$\pm$1.1 Mpc. We compare this result with distance measurements in the literature.

We introduce a method for identifying "twin" Type Ia supernovae, and using them to improve distance measurements. This novel approach to Type Ia supernova standardization is made possible by spectrophotometric time series observations from the Nearby Supernova Factory (SNfactory). We begin with a well-measured set of supernovae, find pairs whose spectra match well across the entire optical window, and then test whether this leads to a smaller dispersion in their absolute brightnesses. This analysis is completed in a blinded fashion, ensuring that decisions made in implementing the method do not inadvertently bias the result. We find that pairs of supernovae with more closely matched spectra indeed have reduced brightness dispersion. We are able to standardize this initial set of SNfactory supernovae to 0.083 +/- 0.012 magnitudes, implying a dispersion of 0.072 +/- 0.010 magnitudes in the absence of peculiar velocities. We estimate that with larger numbers of comparison SNe, e.g, using the final SNfactory spectrophotometric dataset as a reference, this method will be capable of standardizing high-redshift supernovae to within 0.06-0.07 magnitudes. These results imply that at least 3/4 of the variance in Hubble residuals in current supernova cosmology analyses is due to previously unaccounted-for astrophysical differences among the supernovae

Context. Our Local Group of galaxies appears to be moving relative to the cosmic microwave background with the source of the peculiar motion still uncertain. While in the past this has been studied mostly using galaxies as distance indicators, the weight of type Ia supernovae (SNe Ia) has increased recently with the continuously improving statistics of available low-redshift supernovae. Aims. We measured the bulk flow in the nearby universe ($0.015 < z < 0.1$) using 117 SNe Ia observed by the Nearby Supernova Factory, as well as the Union2 compilation of SN Ia data already in the literature. Methods. The bulk flow velocity was determined from SN data binned in redshift shells by including a coherent motion (dipole) in a cosmological fit. Additionally, a method of spatially smoothing the Hubble residuals was used to verify the results of the dipole fit. To constrain the location and mass of a potential mass concentration (e.g., the Shapley supercluster) responsible for the peculiar motion, we fit a Hubble law modified by adding an additional mass concentration. Results. The analysis shows a bulk flow that is consistent with the direction of the CMB dipole up to $z \sim 0.06$, thereby doubling the volume over which conventional distance measures are sensitive to a bulk flow. We see no significant turnover behind the center of the Shapley supercluster. A simple attractor model in the proximity of the Shapley supercluster is only marginally consistent with our data, suggesting the need for another, more distant source. In the redshift shell $0.06 < z < 0.1$, we constrain the bulk flow velocity to $< 240~\textrm{km s}^{-1}$ (68% confidence level) for the direction of the CMB dipole, in contradiction to recent claims of the existence of a large-amplitude dark flow.

Previously we used the Nearby Supernova Factory sample to show that SNe~Ia having locally star-forming environments are dimmer than SNe~Ia having locally passive environments.Here we use the \constitution\ sample together with host galaxy data from \GALEX\ to independently confirm that result. The effect is seen using both the SALT2 and MLCS2k2 lightcurve fitting and standardization methods, with brightness differences of $0.094 \pm 0.037\ \mathrm{mag}$ for SALT2 and $0.155 \pm 0.041\ \mathrm{mag}$ for MLCS2k2 with $R_V=2.5$. When combined with our previous measurement the effect is $0.094 \pm 0.025\ \mathrm{mag}$ for SALT2. If the ratio of these local SN~Ia environments changes with redshift or sample selection, this can lead to a bias in cosmological measurements. We explore this issue further, using as an example the direct measurement of $H_0$. \GALEX{} observations show that the SNe~Ia having standardized absolute magnitudes calibrated via the Cepheid period--luminosity relation using {\textit{HST}} originate in predominately star-forming environments, whereas only ~50% of the Hubble-flow comparison sample have locally star-forming environments. As a consequence, the $H_0$ measurement using SNe~Ia is currently overestimated. Correcting for this bias, we find a value of $H_0^{corr}=70.6\pm 2.6\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ when using the LMC distance, Milky Way parallaxes and the NGC~4258 megamaser as the Cepheid zeropoint, and $68.8\pm 3.3\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ when only using NGC~4258. Our correction brings the direct measurement of $H_0$ within $\sim 1\,\sigma$ of recent indirect measurements based on the CMB power spectrum.

We estimate systematic errors due to K-corrections in standard photometric analyses of high redshift Type Ia supernovae. Errors due to K-correction occur when the spectral template model underlying the lightcurve fitter poorly represents the actual supernova spectral energy distribution, meaning that the distance modulus cannot be recovered accurately. In order to quantify this effect, synthetic photometry is performed on artificially redshifted spectrophotometric data from 119 low-redshift supernovae from the Nearby Supernova Factory, and the resulting lightcurves are fit with a conventional lightcurve fitter. We measure the variation in the standardized magnitude that would be fit for a given supernova if located at a range of redshifts and observed with various filter sets corresponding to current and future supernova surveys. We find significant variation in the measurements of the same supernovae placed at different redshifts regardless of filters used, which causes dispersion greater than $\sim0.05$ mag for measurements of photometry using the Sloan-like filters and a bias that corresponds to a $0.03$ shift in $w$ when applied to an outside data set. To test the result of a shift in supernova population or environment at higher redshifts, we repeat our calculations with the addition of a reweighting of the supernovae as a function of redshift and find that this strongly affects the results and would have repercussions for cosmology. We discuss possible methods to reduce the contribution of the K-correction bias and uncertainty.

We develop a new framework for use in exploring Type Ia Supernova (SN Ia) spectra. Combining Principal Component Analysis (PCA) and Partial Least Square analysis (PLS) we are able to establish correlations between the Principal Components (PCs) and spectroscopic/photometric SNe Ia features. The technique was applied to ~120 supernova and ~800 spectra from the Nearby Supernova Factory. The ability of PCA to group together SNe Ia with similar spectral features, already explored in previous studies, is greatly enhanced by two important modifications: (1) the initial data matrix is built using derivatives of spectra over the wavelength, which increases the weight of weak lines and discards extinction, and (2) we extract time evolution information through the use of entire spectral sequences concatenated in each line of the input data matrix. These allow us to define a stable PC parameter space which can be used to characterize synthetic SN Ia spectra by means of real SN features. Using PLS, we demonstrate that the information from important previously known spectral indicators (namely the pseudo-equivalent width (pEW) of Si II 5972 / Si II 6355 and the line velocity of S II 5640 / Si II 6355) at a given epoch, is contained within the PC space and can be determined through a linear combination of the most important PCs. We also show that the PC space encompasses photometric features like B or V magnitudes, B-V color and SALT2 parameters c and x1. The observed colors and magnitudes, that are heavily affected by extinction, cannot be reconstructed using this technique alone. All the above mentioned applications allowed us to construct a metric space for comparing synthetic SN Ia spectra with observations.

Observational cosmology employing optical surveys often require precise flux calibration. In this context we present SNIFS Calibration Apparatus (SCALA), a flux calibration system developed for the SuperNova Integral Field Spectrograph (SNIFS), operating at the University of Hawaii 2.2 m telescope. SCALA consists of a hexagonal array of 18 small parabolic mirrors distributed over the face of, and feeding parallel light to, the telescope entrance pupil. The mirrors are illuminated by integrating spheres and a wavelength-tunable (from UV to IR) light source, generating light beams with opening angles of 1 degree. These nearly parallel beams are flat and flux-calibrated at a subpercent level, enabling us to calibrate our "telescope + SNIFS system" at the required precision.

We present a sample of normal type Ia supernovae from the Nearby Supernova Factory dataset with spectrophotometry at sufficiently late phases to estimate the ejected mass using the bolometric light curve. We measure $^{56}$Ni masses from the peak bolometric luminosity, then compare the luminosity in the $^{56}$Co-decay tail to the expected rate of radioactive energy re- lease from ejecta of a given mass. We infer the ejected mass in a Bayesian context using a semi-analytic model of the ejecta, incorporating constraints from contemporary numerical models as priors on the density structure and distribution of $^{56}$Ni throughout the ejecta. We find a strong correlation between ejected mass and light curve decline rate, and consequently $^{56}$Ni mass, with ejected masses in our data ranging from 0.9-1.4 $M_\odot$. Most fast-declining (SALT2 $x_1 < -1$) normal SNe Ia have significantly sub-Chandrasekhar ejected masses in our fiducial analysis.

Kim et al. (2013) [K13] introduced a new methodology for determining peak-brightness absolute magnitudes of type Ia supernovae from multi-band light curves. We examine the relation between their parameterization of light curves and Hubble residuals, based on photometry synthesized from the Nearby Supernova Factory spectrophotometric time series, with global host-galaxy properties. The K13 Hubble residual step with host mass is $0.013\pm 0.031$ mag for a supernova subsample with data coverage corresponding to the K13 training; at $\ll 1\sigma$, the step is not significant and lower than previous measurements. Relaxing the data coverage requirement the Hubble residual step with host mass is $0.045\pm 0.026$ mag for the larger sample; a calculation using the modes of the distributions, less sensitive to outliers, yields a step of 0.019 mag. The analysis of this article uses K13 inferred luminosities, as distinguished from previous works that use magnitude corrections as a function of SALT2 color and stretch parameters: Steps at $>2\sigma $ significance are found in SALT2 Hubble residuals in samples split by the values of their K13 $x(1)$ and $x(2)$ light-curve parameters. $x(1)$ affects the light-curve width and color around peak (similar to the $\Delta m_{15}$ and stretch parameters), and $x(2)$ affects colors, the near-UV light-curve width, and the light-curve decline 20 to 30 days after peak brightness. The novel light-curve analysis, increased parameter set, and magnitude corrections of K13 may be capturing features of SN~Ia diversity arising from progenitor stellar evolution.

(Abridged) We study the host galaxy regions in close proximity to Type Ia supernovae (SNe Ia) to analyze relations between the properties of SN Ia events and environments most similar to where their progenitors formed. We focus on local H\alpha\ emission as an indicator of young environments. The Nearby Supernova Factory has obtained flux-calibrated spectral timeseries for SNe Ia using integral field spectroscopy, allowing the simultaneous measurement of the SN and its immediate vicinity. For 89 SNe Ia we measure H\alpha\ emission tracing ongoing star formation within a 1 kpc radius around each SN. This constitutes the first direct study of the local environment for a large sample of SNe Ia also having accurate luminosity, color and stretch measurements. We find that SNe Ia with local H\alpha\ emission are redder by 0.036+/-0.017 mag, and that the previously-noted correlation between stretch and host mass is entirely driven by the SNe Ia coming from passive regions. Most importantly, the mean standardized brightness for SNe Ia with local H\alpha\ emission is 0.094+/-0.031 mag fainter than for those without. This offset arises from a bimodal structure in the Hubble residuals, that also explains the previously-known host-mass bias. We combine this bimodality with the cosmic star-formation rate to predict changes with redshift in the mean SN Ia brightness and the host-mass bias. This change is confirmed using high-redshift SNe Ia from the literature. These environmental dependences point to remaining systematic errors in SNe Ia standardization. The observed brightness offset is predicted to cause a significant bias in measurements of the dark energy equation of state. Recognition of these effects offers new opportunities to improve SNe Ia as cosmological probes - e.g. SNe Ia having local H\alpha\ emission are more homogeneous, having a brightness dispersion of 0.105+/-0.012 mag.

We present photometric and spectroscopic observations of galaxies hosting Type Ia supernovae (SNe Ia) observed by the Nearby Supernova Factory (SNfactory). Combining GALEX UV data with optical and near infrared photometry, we employ stellar population synthesis techniques to measure SN Ia host galaxy stellar masses, star-formation rates (SFRs), and reddening due to dust. We reinforce the key role of GALEX UV data in deriving accurate estimates of galaxy SFRs and dust extinction. Optical spectra of SN Ia host galaxies are fitted simultaneously for their stellar continua and emission lines fluxes, from which we derive high precision redshifts, gas-phase metallicities, and Halpha-based SFRs. With these data we show that SN Ia host galaxies present tight agreement with the fiducial galaxy mass-metallicity relation from SDSS for stellar masses log(M_*/M_Sun)>8.5 where the relation is well-defined. The star-formation activity of SN Ia host galaxies is consistent with a sample of comparable SDSS field galaxies, though this comparison is limited by systematic uncertainties in SFR measurements. Our analysis indicates that SN Ia host galaxies are, on average, typical representatives of normal field galaxies.

We examine the relationship between Type Ia Supernova (SN Ia) Hubble residuals and the properties of their host galaxies using a sample of 115 SNe Ia from the Nearby Supernova Factory (SNfactory). We use host galaxy stellar masses and specific star-formation rates fitted from photometry for all hosts, as well as gas-phase metallicities for a subset of 69 star-forming (non-AGN) hosts, to show that the SN Ia Hubble residuals correlate with each of these host properties. With these data we find new evidence for a correlation between SN Ia intrinsic color and host metallicity. When we combine our data with those of other published SN Ia surveys, we find the difference between mean SN Ia brightnesses in low and high mass hosts is 0.077 +- 0.014 mag. When viewed in narrow (0.2 dex) bins of host stellar mass, the data reveal apparent plateaus of Hubble residuals at high and low host masses with a rapid transition over a short mass range (9.8 <= log(M_*/M_Sun) <= 10.4). Although metallicity has been a favored interpretation for the origin of the Hubble residual trend with host mass, we illustrate how dust in star-forming galaxies and mean SN Ia progenitor age both evolve along the galaxy mass sequence, thereby presenting equally viable explanations for some or all of the observed SN Ia host bias.

We present a novel class of models for Type Ia supernova time-evolving spectral energy distributions (SED) and absolute magnitudes: they are each modeled as stochastic functions described by Gaussian processes. The values of the SED and absolute magnitudes are defined through well-defined regression prescriptions, so that data directly inform the models. As a proof of concept, we implement a model for synthetic photometry built from the spectrophotometric time series from the Nearby Supernova Factory. Absolute magnitudes at peak $B$ brightness are calibrated to 0.13 mag in the $g$-band and to as low as 0.09 mag in the $z=0.25$ blueshifted $i$-band, where the dispersion includes contributions from measurement uncertainties and peculiar velocities. The methodology can be applied to spectrophotometric time series of supernovae that span a range of redshifts to simultaneously standardize supernovae together with fitting cosmological parameters.

We present 32 epochs of optical (3300-9700 \AA) spectrophotometric observations of the nearby quintessential "normal" type Ia supernova (SN Ia) SN 2011fe in the galaxy M101, extending from -15 to +97 d with respect to B-band maximum, obtained by the Nearby Supernova Factory collaboration. SN 2011fe is the closest (\mu = 29.04) and brightest (Bmax = 9.94 mag) SN Ia observed since the advent of modern large scale programs for the intensive periodic followup of supernovae. Both synthetic light curve measurements and spectral feature analysis attest to the normality of SN 2011fe. There is very little evidence for reddening in its host galaxy. The homogeneous calibration, intensive time sampling, and high signal-to-noise ratio of the data set make it unique. Thus it is ideal for studying the physics of SN Ia explosions in detail, and for furthering the use of SNe Ia as standardizable candles for cosmology. Several such applications are shown, from the creation of a bolometric light curve and measurement of the 56Ni mass, to the simulation of detection thresholds for unburned carbon, direct comparisons with other SNe Ia, and existing spectral templates.

We present a new atmospheric extinction curve for Mauna Kea spanning 3200--9700 \AA. It is the most comprehensive to date, being based on some 4285 standard star spectra obtained on 478 nights spread over a period of 7 years obtained by the Nearby SuperNova Factory using the SuperNova Integral Field Spectrograph. This mean curve and its dispersion can be used as an aid in calibrating spectroscopic or imaging data from Mauna Kea, and in estimating the calibration uncertainty associated with the use of a mean extinction curve. Our method for decomposing the extinction curve into physical components, and the ability to determine the chromatic portion of the extinction even on cloudy nights, is described and verified over the wide range of conditions sampled by our large dataset. We demonstrate good agreement with atmospheric science data obtain at nearby Mauna Loa Observatory, and with previously published measurements of the extinction above Mauna Kea.

We present optical photometry and spectroscopy of five type Ia supernovae discovered by the Nearby Supernova Factory selected to be spectroscopic analogues of the candidate super-Chandrasekhar-mass events SN 2003fg and SN 2007if. Their spectra are characterized by hot, highly ionized photospheres near maximum light, for which SN 1991T supplies the best phase coverage among available close spectral templates. Like SN 2007if, these supernovae are overluminous (-19.5 < M_V < -20) and the velocity of the Si II 6355 absorption minimum is consistent with being constant in time from phases as early as a week before, and up to two weeks after, $B$-band maximum light. We interpret the velocity plateaus as evidence for a reverse-shock shell in the ejecta formed by interaction at early times with a compact envelope of surrounding material, as might be expected for SNe resulting from the mergers of two white dwarfs. We use the bolometric light curves and line velocity evolution of these SNe to estimate important parameters of the progenitor systems, including nickel-56 mass, total progenitor mass, and masses of shells and surrounding carbon/oxygen envelopes. We find that the reconstructed total progenitor mass distribution of the events (including SN 2007if) is bounded from below by the Chandrasekhar mass, with SN 2007if being the most massive. We discuss the relationship of these events to the emerging class of super-Chandrasekhar-mass SNe Ia, estimate the relative rates, compare the mass distribution to that expected for double-degenerate SN Ia progenitors from population synthesis, and consider implications for future cosmological Hubble diagrams.

The nearby supernova SN 2011fe can be observed in unprecedented detail. Therefore, it is an important test case for Type Ia supernova (SN Ia) models, which may bring us closer to understanding the physical nature of these objects. Here, we explore how available and expected future observations of SN 2011fe can be used to constrain SN Ia explosion scenarios. We base our discussion on three-dimensional simulations of a delayed detonation in a Chandrasekhar-mass white dwarf and of a violent merger of two white dwarfs-realizations of explosion models appropriate for two of the most widely-discussed progenitor channels that may give rise to SNe Ia. Although both models have their shortcomings in reproducing details of the early and near-maximum spectra of SN 2011fe obtained by the Nearby Supernova Factory (SNfactory), the overall match with the observations is reasonable. The level of agreement is slightly better for the merger, in particular around maximum, but a clear preference for one model over the other is still not justified. Observations at late epochs, however, hold promise for discriminating the explosion scenarios in a straightforward way, as a nucleosynthesis effect leads to differences in the 55Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect.

We present convincing evidence of unburned carbon at photospheric velocities in new observations of 5 Type Ia supernovae (SNe Ia) obtained by the Nearby Supernova Factory. These SNe are identified by examining 346 spectra from 124 SNe obtained before +2.5 d relative to maximum. Detections are based on the presence of relatively strong C II 6580 absorption "notches" in multiple spectra of each SN, aided by automated fitting with the SYNAPPS code. Four of the 5 SNe in question are otherwise spectroscopically unremarkable, with ions and ejection velocities typical of SNe Ia, but spectra of the fifth exhibits high-velocity (v > 20,000 km/s) Si II and Ca II features. On the other hand, the light curve properties are preferentially grouped, strongly suggesting a connection between carbon-positivity and broad band light curve/color behavior: Three of the 5 have relatively narrow light curves but also blue colors, and a fourth may be a dust-reddened member of this family. Accounting for signal-to-noise and phase, we estimate that 22 +10/-6% of SNe Ia exhibit spectroscopic C II signatures as late as -5 d with respect to maximum. We place these new objects in the context of previously recognized carbon-positive SNe Ia, and consider reasonable scenarios seeking to explain a physical connection between light curve properties and the presence of photospheric carbon. We also examine the detailed evolution of the detected carbon signatures and the surrounding wavelength regions to shed light on the distribution of carbon in the ejecta. Our ability to reconstruct the C II 6580 feature in detail under the assumption of purely spherical symmetry casts doubt on a "carbon blobs" hypothesis, but does not rule out all asymmetric models. A low volume filling factor for carbon, combined with line-of-sight effects, seems unlikely to explain the scarcity of detected carbon in SNe Ia by itself.