Draft version January 8, 2020 Typeset using LATEX twocolumn style in AASTeX63 ZTF Early Observations of Type Ia Supernovae II: First Light, the Initial Rise, and Time to Reach Maximum Brightness A. A. Miller,1, 2 Y. Yao,3 M. Bulla,4 C. Pankow,1 E. C. Bellm,5 S. B. Cenko,6, 7 R. Dekany,8 C. Fremling,3 M. J. Graham,9 T. Kupfer,10 R. R. Laher,11 A. A. Mahabal,9, 12 F. J. Masci,11 P. E. Nugent,13, 14 R. Riddle,8 B. Rusholme,11 R. M. Smith,8 D. L. Shupe,11 J. van Roestel,9 and S. R. Kulkarni3 1Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA) and Department of Physics and Astronomy, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA 2The Adler Planetarium, Chicago, IL 60605, USA 3Cahill Center for Astrophysics, California Institute of Technology, 1200 E. California Boulevard, Pasadena, CA 91125, USA 4Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden 5DIRAC Institute, Department of Astronomy, University of Washington, 3910 15th Avenue NE, Seattle, WA 98195, USA 6Astrophysics Science Division, NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD 20771, USA 7Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA 8Caltech Optical Observatories, California Institute of Technology, Pasadena, CA 91125, USA 9Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 10Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 11IPAC, California Institute of Technology, 1200 E. California Blvd, Pasadena, CA 91125, USA 12Center for Data Driven Discovery, California Institute of Technology, Pasadena, CA 91125, USA 13Computational Cosmology Center, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 14Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA (Received January 8, 2020) Submitted to ApJ Abstract While it is clear that Type Ia supernovae (SNe) are the result of thermonuclear explosions in C/O white dwarfs (WDs), a great deal remains uncertain about the binary companion that facilitates the explosive disruption of the WD. Here, we present a comprehensive analysis of a unique, and large, data set of 127 SNe Ia with exquisite coverage by the Zwicky Transient Facility (ZTF). High-cadence (6 observations per night) ZTF observations allow us to measure the SN rise time and examine its initial evolution. We develop a Bayesian framework to model the early rise as a power-law in time, which enables the inclusion of priors in our model. For a volume-limited subset of normal SNe Ia, we find the mean power-law index is consistent with 2 in the rZTF-band (αr = 2:01 ± 0:02), as expected in the expanding fireball model. There are, however, individual SNe that are clearly inconsistent with αr = 2. We estimate a mean rise time of 18.5 d (with a range extending from ∼15{22 d), though this is subject to the adopted prior. We identify an important, previously unknown, bias whereby the rise arXiv:2001.00598v2 [astro-ph.HE] 7 Jan 2020 times for higher redshift SNe within a flux-limited survey are systematically underestimated. This effect can be partially alleviated if the power-law index is fixed to α = 2, in which case we estimate a mean rise time of 21.0 d (with a range from ∼18{23 d). The sample includes a handful or rare and peculiar SNe Ia. Finally,DRAFT we conclude with a discussion of lessons learned from the ZTF sample that can eventually be applied to Large Synoptic Survey Telescope observations. [email protected] 2 Miller et al. Keywords: supernovae: general | methods: observational | surveys | catalogs | supernovae: individual (ZTF18abclfee/SN 2018cxk) 1. Introduction bink 1984). While the debate long focused on which of The fact that supernovae (SNe) of Type Ia can be these two scenarios is correct, there is now strong em- empirically calibrated as standardizable candles makes pirical evidence in support of both channels. PTF 11kx, them arguably the most important tool in all of physics an extreme example of a SN Ia, showed clear evidence for the past ∼two decades. By unlocking our ability of multiple shells of circumstellar material (Dilday et al. to accurately measure distances at high redshift, SNe 2012), which is precisely what one would expect in a Ia have completely revolutionized our understanding of WD+red giant system that had undergone multiple no- the Universe (Riess et al. 1998; Perlmutter et al. 1999). vae prior to the final, fatal thermonuclear explosion. On While it is all but certain that SNe Ia are the result of the other hand, hypervelocity WDs discovered by Gaia thermonuclear explosions in carbon-oxygen (C/O) white are likely the surviving companions of DD explosions dwarfs (WDs) in binary star systems (see Maoz et al. (Shen et al. 2018). 2014; Livio & Mazzali 2018 for recent reviews), there There is also emerging evidence that WDs can explode remains a great deal about SNe Ia progenitors and the prior to reaching MCh. Detailed modeling of SNe Ia precise explosion mechanism that we do not know. This light curves (Scalzo et al. 2014) and a blue-to-red-to-blue leads to the tantalizing hope that an improved under- color evolution observed in a few SNe (Jiang et al. 2017; standing of the binary companions or explosion may im- Noebauer et al. 2017; Polin et al. 2019; De et al. 2019), prove our ability to calibrate these explosions as stan- point to sub-MCh mass explosions. Such explosions are dardizable candles. possible if a C/O WD accretes and retains a thick He One clear avenue for better understanding the progen- shell. A detonation in this shell can trigger an explosion itors of SNe Ia is to obtain observations in the hours to in the C/O core of the WD (e.g., Nomoto 1982a,b). days after explosion (Maoz et al. 2014). Such detections Now the most pressing questions are the following: provide an opportunity to probe the progenitor environ- which scenario, SD or DD, and which mass explosion, ment and binary companion, which is simply not pos- MCh or sub-MCh, is dominant? sible once the SN evolves well into the expansion phase Here too, early observations should prove useful. In (they are standardizable precisely because they are all the SD scenario the SN ejecta will collide with the non- nearly identical at maximum light). Indeed, in the land- degenerate companion creating a shock that gives rise mark discovery of SN 2011fe by the Palomar Transient to a bright ultraviolet/optical flash in the days after Factory (PTF; Law et al. 2009; Rau et al. 2009), Nugent explosion (Kasen 2010). To date, the search for such et al.(2011) were able to constrain the time of explo- a signature in large samples typically results in a non- sion to ±20 min (though see Piro & Nakar 2013 for an detection (e.g., Hayden et al. 2010; Ganeshalingam et al. 1 explanation of a potential early \dark phase"). Bloom 2011; Bianco et al. 2011). For the DD scenario, some et al.(2012) would later combine the PTF observations sub-MCh DD explosions exhibit a highly unusual color with an early non-detection while comparing the limits evolution, though this can only be observed a few days to shock-breakout models to constrain the size of the after explosion (Noebauer et al. 2017; Polin et al. 2019). Measurements of the rise time, i.e., the time it takes progenitor to be . 0:01R , providing the most direct evidence to date that SNe Ia come from WDs. These the SN to evolve from first light, the moment when op- critical findings have demonstrated the importance of tical photons diffuse out of the photosphere, to max- high-cadence time-domain surveys for discovering SNe imum brightness, can also play a role in constraining shortly after explosion. the progenitor systems of SNe Ia. Initial work to esti- Pinning down the binary companion to the exploding WD remains a challenge. There are likely two dominant 1 There are claims of companion interaction based on short-lived pathways towards explosion. In the first, the WD ac- optical \bumps" in the early light curves for individual SNe (e.g., cretes H from a non-degenerate companion and eventu- Cao et al. 2015; Marion et al. 2016; Hosseinzadeh et al. 2017; Dimitriadis et al. 2019). Alternative models (e.g., Dessart et al. ally explodes as it approaches the Chandrasekhar mass 2014; Piro & Morozova 2016) utilizing different physical scenarios (MCh; known as the single degenerate, or SD, scenario; can produce similar bumps, leading many (e.g., Kromer et al. Whelan & Iben 1973), while in the second the explo- 2016; Shappee et al. 2016a; Noebauer et al. 2017; Miller et al. 2018; Shappee et al. 2019) to appeal to alternative explanations sion follows the interaction or merger of two WD stars than ejecta-companion interaction. (known as the double degenerate, or DD, scenario; Web- The Rise of ZTF SNe Ia 3 mate the rise times of SNe Ia clearly demonstrated that servations discovered by the Zwicky Transient Facility early efforts to model WD explosions significantly un- (ZTF; Bellm et al. 2019a; Graham et al. 2019) in 2018, derestimated the opacities in the SN ejecta (e.g., Riess we examine the rise time of SNe Ia and whether or not et al. 1999). Furthermore, while the famous luminosity- their early emission can be characterized as a simple decline relationship for SNe Ia makes them standardiz- power law. Paper I (Yao et al. 2019) describes the sam- able (Phillips 1993), recent evidence suggests that the ple, while Paper III (M.