## CONSTRAINTS ON THE PHYSICAL PROPERTIES OF TYPE IA SUPERNOVAE FROM PHOTOMETRY

R. K¨nyves-T´th, 1 o o J. Vink´, 1,2 o A. Ordasi, 1 K. S´rneczky, 1 a A. B´di, 1,3 o B. Cseh, 1 G. Cs¨rnyei, 1 o Z. Dencs, 1 O. Hanyecz, 1 B. Ign´cz, 1 a Cs. Kalup, 1 L. Kriskovics, 1 A. P´l, 1 a B. Seli, 1 ´ . A S´ odor, 1 R. Szak´ts, 1 a P. Sz´kely, 4 e E. Varga-Vereb´ elyi, 1 K. Vida, 1 and G. Zsidi 1

- 1 Konkoly Observatory, CSFK, Konkoly-Thege M. ut 15-17, Budapest, Hungary
- 2 Department of Optics and Quantum Electronics, University of Szeged, Hungary
- 3 CSFK Lend¨ ulet Near-Field Cosmology Research Group
- 4 Department of Experimental Physics,University of Szeged, D´m t´r 9., Szeged, Hungary o e

## ABSTRACT

We present a photometric study of 17 Type Ia supernovae (SNe) based on multi-color (Johnson-Cousins-Bessell BVRI ) data taken at Piszk´stet˝ mountain station of Konkoly Observatory, Hungary between 2016 and 2018. e o We analyze the light curves (LCs) using the publicly available LC-fitter SNooPy2 to derive distance and reddening information. The bolometric LCs are fit with a radiation-diffusion Arnett model to get constraints on the physical parameters of the ejecta: the optical opacity, the ejected mass and the initial nickel mass in particular. We also study the pre-maximum, de-reddened ( B -V ) 0 color evolution by comparing our data with standard delayed detonation and pulsational delayed detonation models, and show that the 56 Ni masses of the models that fit the ( B -V ) 0 colors are consistent with those derived from the bolometric LC fitting. We find similar correlations between the ejecta parameters (e.g. ejecta mass, or 56 Ni mass vs decline rate) as published recently by Scalzo et al. (2019).

Keywords: supenovae: general - supernovae: individual (Gaia16alq, SN 2016asf, SN 2016bln, SN 2016coj, SN 2016eoa, SN 2016ffh, SN 2016gcl, SN 2016gou, SN 2016ixb, SN 2017cts, SN 2017drh, SN 2017erp, SN 2017fgc, SN 2017fms, SN 2017hjy, SN 2017igf, SN 2018oh)

## 1. INTRODUCTION

Type Ia supernovae (SNe) are especially important objects for measuring extragalactic distances as their peak absolute magnitudes can be inferred via fitting their observed, multi-color LCs. Normal SN Ia events obey the empirical Phillips-relation (Pskovskii 1977; Phillips 1993), which states that the LCs of intrinsically fainter objects decline faster than those of brighter ones. The decline rate is often parametrized by ∆ m 15 , i.e. the magnitude difference between the peak and the one measured at 15 days after maximum in a given (often the B ) band. For example, the earlier version of the SNooPy code (Burns et al. 2011) applied ∆ m 15 as a fitting parameter for the decline rate. In the new version, SNooPy2 (Burns et al. 2014, 2018), a new parameter ( s BV ) that measures the time difference between the maxima of the B -band light curve and the B -V color curve, was introduced. Other parametrizations also exist: for example the SALT2 code (Guy et al. 2007, 2010; Betoule et al. 2014) applies the x 1 (stretch) parameter, while MLCS2k2 (Riess et al. 1998; Jha et al. 1999, 2007) uses ∆ that corresponds to the magnitude difference between the actual V -band maximum brightness and that of a fiducial Type Ia SN. All of them are based on the same Phillips-relation, thus, ∆ m 15 , s BV , x 1 or ∆ are related to each other.

Studying SNe Ia opens a door for constraining the Hubble-parameter H 0 (Riess et al. 2012, 2016; Dhawan et al. 2018a) by getting accurate distances to their host galaxies. Such absolute distances are the quintessential cornerstones of the cosmic distance ladder. Via constraining H 0 , SNe Ia play a major role in investigating the expansion of the Universe (Riess et al. 1998; Perlmutter et al. 1999; Astier et al. 2006; Riess et al. 2007; Wood-Vasey et al. 2007; Kessler et al. 2009; Guy et al. 2010; Conley et al. 2011; Betoule et al. 2014; Rest et al. 2014; Scolnic et al. 2014; Bengaly et al. 2015; Jones et al. 2015; Li et al. 2016; Zhang et al. 2017) and testing the most recent cosmological models (e.g. Benitez-Herrera et al. 2013; Betoule et al. 2014).

Even though they are extensively used to estimate distances, the improvement of the precision as well as the accuracy of the method is still a subject of recent studies (see e.g. Vink´ et al. 2018). o In order to achieve the desired 1% accuracy, it is important to understand the physical properties of the progenitor system and the explosion mechanism better.

The actual progenitor that explodes as a SN Ia, as well as the explosion mechanism, is still an issue. There are two main proposed progenitor scenarios: singledegenerate (SD) (Whelan &amp; Iben 1973) and double- degenerate (DD) (Iben &amp; Tutukov 1984). The SD scenario presumes that a carbon-oxygen white dwarf (C/O WD) has a non-degenerate companion star, e.g. a red giant, which, after overflowing its Roche-lobe, transfers mass to the WD. When the WD approaches the Chandrasekhar mass ( M Ch ), spontaneous fusion of C/O to 56 Ni develops that quickly engulfs the whole WD, leading to a thermonuclear explosion.

The details on the onset and the progress of the C/O fusion is still debated, and many possible mechanisms have been proposed in the literature. The most successful one is the delayed detonation explosion (DDE) model, in which the burning starts as a deflagration, but later it turns into a detonation wave (Nomoto et al. 1984; Khokhlov 1991; Dessart et al. 2014; Maoz et al. 2014). A variant of that is the pulsational delayed detonation explosion (PDDE): during the initial deflagration phase the expansion of the WD expels some material from its outmost layers, which pulsates, expands and avoids burning. After that, the bound material falls back to the WD that leads to a subsequent detonation (Dessart et al. 2014).

There is a theoretical possibility for a sub-Chandrasekhar double-detonation scenario, where the WD accretes a thin layer of helium onto its surface, which is compressed by its own mass that leads to He-detonation This triggers the thermonuclear explosion of the underlying C/O WD (Woosley &amp; Weaver 1994; Fink et al. 2010; Kromer et al. 2010; Sim et al. 2010, 2012).

In the DD scenario two WDs merge or collide that results in a subsequent explosion (Maoz et al. 2014; van Rossum et al. 2016).

It may be possible to distinguish between these scenarios e.g. by constraining the mass of the progenitor. While the traditional SD and DD scenarios require nearly M Ch WD masses, explosion of subM Ch WDs may be possible via the 'violent merger' or the double detonation mechanisms (see e.g. Maoz et al. 2014, for a detailed review). Thus, the ejecta mass is an extremely important physical quantity, which can be inferred by fitting LC models to the observations.

The idea that the bolometric LC of SNe Ia can be used to infer the ejecta mass via a semi-analytical model, was introduced by Arnett (1982) and developed further by Jeffery (1999). Arnett (1982) showed that the ejecta mass correlates with the rise time to maximum light, provided the expansion velocity and the mean optical opacity of the ejecta are known. Later, Jeffery (1999) suggested the usage of the rate of the deviation of the observed LC from the rate of the Co-decay during the early nebular phase (i.e. the transparency timescale, t γ ), which measures the leakage of γ -photons from the

diluting ejecta. The advantage of using t γ for constraining the ejecta mass is that t γ is proportional to the gamma-ray opacity, which is much better known than the mean optical opacity. This technique was applied to real data by Stritzinger et al. (2006), Scalzo et al. (2014) and more recently by Scalzo et al. (2019). The conclusion of all these studies was that most SNe Ia seem to have subM Ch ejecta, which was also confirmed recently by theoretical models (Dhawan et al. 2018b; Goldstein, &amp; Kasen 2018; Wilk et al. 2018; Papadogiannakis et al. 2019). Other studies lead to similar conclusions: for example, Dhawan et al. (2016, 2017) studied the correlation of the bolometric peak luminosity with the phase of the second maximum in NIR, and found that subM Ch explosions can be separated into two categories in respect of this property, however superM Ch explosions do not follow such a relation.

The usage of the transparency timescale as a proxy for the ejecta mass has some caveats, though. t γ also depends on the characteristic velocity ( v e ) of the expanding ejecta, which is not easy to constrain as it is related to the velocity of a layer deep inside the ejecta. Stritzinger et al. (2006), for example, assumed that v e is uniform for all SNe Ia (they adopted v e ∼ 3000 km s -1 ), which may not be true in reality, because it is known that a diversity in expansion velocities exists for most SNe Ia (Wang et al. 2013). Another issue is the distribution of the radioactive 56 Ni, encoded by the q parameter by Jeffery (1999), which is usually taken from models: q ∼ 1 / 3, as given by the W7 model, is often assumed. These assumptions, although may not be too far from reality, might introduce some sort of systematic uncertainties in the inferred ejecta masses, which may be worth for further studies.

Additionally, Levanon, &amp; Soker (2019) investigated the influence of asymmetric 56 Ni distribution to t γ , and found that the full range of t γ -values cannot be explained by variations in this asymmetry. Their suggested that both M Ch and subM Ch explosions are needed to account for the diversity of the observed physical parameters of SN Ia explosions.

The main motivation of the present paper is to give constraints on the ejecta mass and some other physical parameters for a sample of 17 recent SNe Ia (Figure 1 and 2) observed from Piszk´ estet˝ station of Konkoly o Observatory, Hungary. We generalize the prescription of inferring the ejecta masses by combining the LC rise time and the transparency timescale within the framework of the constant-density Arnett model.

In the following we present the description of the photometric sample (Section 2), then we show the results from multi-color LC modeling (Section 3). We construct and fit the bolometric LCs in order to derive the ejecta mass, and other parameters such as the diffusion- and gamma-leakage timescales, the optical opacity and the expansion velocity.

In Section 4, we first discuss the early (de-reddened) ( B -V ) 0 color evolution, which might also provide some constraints on the explosion mechanism and the progenitor system (e.g. Hosseinzadeh et al. 2017; Miller et al. 2018; Stritzinger et al. 2018). There is a growing number of evidence for the appearance of blue excess light during the earliest phase of some SNe Ia (Marion et al. 2016; Hosseinzadeh et al. 2017; Dimitriadis et al. 2019; Li et al. 2019; Shappee et al. 2019; Stritzinger et al. 2018). At present the cause of this excess emission is debated, and a number of possible explanations were proposed recently (Hosseinzadeh et al. 2017; Miller et al. 2018; Stritzinger et al. 2018). These include i ) SN shock cooling, ii ) interaction with a non-degenerate companion, iii ) presence of high velocity 56 Ni in the outer layers of the ejecta, iv ) interaction with the circumstellar matter (CSM) and v ) differences in the composition or variable opacity.

Furthermore, we compare our measured ( B -V ) 0 colors with the predictions of various explosion models (Dessart et al. 2014), and examine the possible correlations between the derived physical parameters following Scalzo et al. (2014), Scalzo et al. (2019) and Khatami &amp; Kasen (2018).

Finally Section 5 summarizes the results of this paper.

## 2. OBSERVATIONS

We obtained multi-band photometry for 17 bright Type Ia SNe from the Piszk´ estet˝ station of o Konkoly Observatory, Hungary between 2016 and 2018. The selection criteria for the sample were as follows: i ) accessibility from the site (i.e. declination above -15 degree), ii ) sufficiently early (pre-maximum) discovery, iii ) ability for follow-up beyond t ∼ +40 days after maximum, and iv ) low redshift ( z /lessorsimilar 0 05). .

All data were taken with the 0.6/0.9 m Schmidttelescope equipped with a 4096 × 4096 FLI CCD and Bessell BVRI filters, thus, providing a homogeneous, high signal-to-noise data sample of nearby SNe Ia.

Data reduction and photometry was done the same way as described in Vink´ et al. (2018). o The raw data were reduced using IRAF 1 (Image Reduction and Analysis Facility) by completing bias, dark and flatfield

1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. http://iraf.noao.edu

SN 2016coj

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-5a9ec25c-f356-4e0b-9655-e4095274d872.jpg)

SN 2016bln

SN 2016ffh

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-562a03ef-82f6-483d-8961-21bdf09d8d66.jpg)

SN 2016asf

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-51289b13-dc13-499b-827f-c70172fa2caf.jpg)

SN 2016eoa

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-1fb330fd-5aad-4b46-8b41-c5c28adc4119.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-8a22890a-1f6d-4439-aeb8-3edb8397be89.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-5ced3c18-cb73-44ee-b3ca-30793467a1e6.jpg)

SN 2016gou

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-632d631f-56cd-4249-8483-762bfb821567.jpg)

SN 2016ixb

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-c35d0009-3fed-472e-b7a4-77ee64624cee.jpg)

Figure 1. Images of the program SNe observed in 2016. The size of each subframe is 1 7 . × 1 7 arcmin . 2 . The supernova is the central object, while North is up and East is to the left.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-5f09724d-ee19-4b77-a52e-9dfb91740c1d.jpg)

corrections. Geometric registration of the sky frames was made in two steps. First, we used SExtractor (Bertin &amp; Arnouts 1996) for identifying point sources on each frame. Second, the imwcs routine from the wcstools 2 package was applied to assign R.A. and Dec. coordinates to pixels on the CCD frames.

Photometry on each SN and several other local comparison stars was made via PSF-fitting with DAOPHOT in IRAF. Note that due to the strong, variable background from the host galaxy, image subtraction was unavoidable in the case of SN 2016coj, 2016gcl, 2016ixb, 2017drh and 2017hjy. For subtraction we used template frames taken with the same telescope and instrumental setup more than 1 year after the discovery of the SN. In these cases the photometry of the comparison stars was

2 http://tdc-www.harvard.edu/wcstools/

computed on the unsubtracted frames, while for the SN it was done on the host-subtracted frames. Particular attention was paid to keep the flux zero point of the subtracted frame the same as that of the unsubtracted one, thus, getting consistent photometry from both frames. Simple PSF fitting gave acceptable results in the case of the other SNe that suffered less severe contamination from their hosts.

The magnitudes of the local comparison stars were determined from their PS1-photometry 3 after transforming the PS1 g P , r P , i P magnitudes to the JohnsonCousins BVRI system. The zero points of the standard transformation were tied to these magnitudes.

Uncertainties of the magnitudes in each filter are estimated by combining the photometric errors reported

3 https://archive.stsci.edu/panstarrs/

SN 2017cts

SN 2017fgc

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-3dc0dca6-3026-4cd0-8213-d04d4eb50a0d.jpg)

SN 2017erp

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-66d21835-dd1c-438d-a8af-917fd4296e7b.jpg)

SN 2017hjy

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-35b99dc2-6fa3-4bd9-8dee-41d0f3160171.jpg)

SN 2017drh

SN 2017fms

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-e4b8ffa3-2f9c-4d2f-a79a-ad957fee1a1f.jpg)

SN 2017igf

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-16535a97-f7ad-454a-ad8f-9864e82ed025.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-fec9ce21-aa6d-42b1-b33f-62f93f71de6b.jpg)

SN 20180h

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-2da373c2-689d-437c-9ec3-1b028312c18e.jpg)

Figure 2. The same as Figure 1 but for the SNe observed in 2017 and 2018.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-625a0e34-91f5-4a2f-853c-1049a8c81ac4.jpg)

by DAOPHOT with the uncertainties of the standard transformation represented by the rms of the residuals between the measured and catalogued magnitudes of the local comparison stars. The latter is found as the leading source of error in most cases, resulting in 0.071, 0.045, 0.048 and 0.053 mag for the median uncertainty in B -, V -, R - and I -bands, respectively.

The basic data of the observed SNe are collected in Table 1. Plots of the V -band sub-frames centered on the SNe are shown in Fig 1 and 2. After acceptance, all photometric data will be made available via the Open Supernova Catalog 4 .

To fit and analyze the observed LCs, we used the SNooPy2 5 (Burns et al. 2011, 2014) LC fitter, which is based on the Phillips-relation (Phillips 1993). It can be categorized as a 'distance calculator' (Conley et al. 2008), since it provides the absolute distance as a fitting parameter. This was an important criterion for choosing SNooPy2 , because precise distance estimates were necessary for the goals of this paper.

## 3. ANALYSIS

3.1. Distance and reddening estimates via multi-color light curve modeling

4 https://sne.space

We applied both the EBV-model that fits the template LCs by Prieto et al. (2006) to the data in BVRI filters, and the EBV2-model that is based on the uBVgriY JH light curves by Burns et al. (2011). Both of these models relate the distance modulus of a normal SN Ia to its decline rate and color as

$$m _ { X } ( \varphi _ { o } ) \ = \ T _ { Y } ( \varphi, p ) + M _ { Y } ( p ) + \mu _ { 0 } + K _ { X Y } + \\ + R _ { X } \cdot E ( B - V ) _ { M W } + R _ { Y } \cdot E ( B - V ) _ { h o s t }, \ ( 1 )$$

5 http://csp.obs.carnegiescience.edu/data/snpy/snpy

Table 1. Basic data of the observed SNe

| SN name    | Type    | R.A.        | Dec.         | Host galaxy                |        z | Discovery   | Ref.   |
|------------|---------|-------------|--------------|----------------------------|----------|-------------|--------|
| Gaia16alq  | Ia-norm | 18:12:29.36 | +31:16:47.32 | PSO J181229.441+311647.834 | 0.023    | 2016-04-21  | a      |
| SN 2016asf | Ia-norm | 06:50:36.73 | +31:06:45.36 | KUG 0647+311               | 0.021    | 2016-03-06  | b      |
| SN 2016bln | Ia-91T  | 13:34:45.49 | +13:51:14.30 | NGC 5221                   | 0.0235   | 2016-04-04  | c      |
| SN 2016coj | Ia-norm | 12:08:06.80 | +65:10:38.24 | NGC 4125                   | 0.005    | 2016-05-28  | d      |
| SN 2016eoa | Ia-91bg | 00:21:23.10 | +22:26:08.30 | NGC 0083                   | 0.021    | 2016-08-02  | e      |
| SN 2016ffh | Ia-norm | 15:11:49.48 | +46:15:03.22 | MCG +08-28-006             | 0.018    | 2016-08-17  | f      |
| SN 2016gcl | Ia-91T  | 23:37:56.62 | +27:16:37.73 | AGC 331536                 | 0.028    | 2016-09-08  | g      |
| SN 2016gou | Ia-norm | 18:08:06.50 | +25:24:31.32 | PSO J180806.461+252431.916 | 0.016    | 2016-09-22  | h      |
| SN 2016ixb | Ia-91bg | 04:54:00.04 | +01:57:46.62 | NPM1G +01.0158             | 0.028343 | 2016-12-17  | i      |
| SN 2017cts | Ia-norm | 17:03:11.76 | +61:27:26.06 | CGCG 299-048               | 0.02     | 2017-04-02  | j      |
| SN 2017drh | Ia-norm | 17:32:26.05 | +07:03:47.52 | NGC 6384                   | 0.005554 | 2017-05-03  | k      |
| SN 2017erp | Ia-norm | 15:09:14.81 | -11:20:03.20 | NGC 5861                   | 0.0062   | 2017-06-13  | l      |
| SN 2017fgc | Ia-norm | 01:20:14.44 | +03:24:09.96 | NGC 0474                   | 0.008    | 2017-07-11  | m      |
| SN 2017fms | Ia-91bg | 21:20:14.60 | -04:52:51.30 | IC 1371                    | 0.031    | 2017-07-17  | n      |
| SN 2017hjy | Ia-norm | 02:36:02.56 | +43:28:19.51 | PSO J023602.146+432817.771 | 0.007    | 2017-10-14  | o      |
| SN 2017igf | Ia-91bg | 11:42:49.85 | +77:22:12.94 | NGC 3901                   | 0.006    | 2017-11-18  | p      |
| SN 2018oh  | Ia-norm | 09:06:39.54 | +19:20:17.77 | UGC 04780                  | 0.012    | 2018-02-04  | q      |

Note -a: Piascik, &amp; Steele (2016); b: Cruz et al. (2016); c: Miller et al. (2016); d: Zheng et al. (2017); e: Gagliano et al. (2016); f: Tonry et al. (2016); g: Brown (2016); h: Tonry et al. (2016); i: Stanek (2016); j: Brimacombe et al. (2017); k: Valenti et al. (2017); l: Itagaki (2017); m: Sand et al. (2017); n: Gagliano et al. (2017); o: Tonry et al. (2017); p: Stanek (2017); q: Stanek (2018)

where X and Y represent the filter of observed data and the template LC, ϕ is the rest-frame phase of the SN (defined as rest-frame days since the moment of B -band maximum light), ϕ o = ϕ · (1 + z ) is the phase in the observer's frame, m X ( ϕ o ) is the observed magnitude in filter X, T Y ( ϕ, p ) is the template LC as a function of phase, p is a generalized decline rate parameter ( p = ∆ m 15 or p = s BV in the EBV and EBV2 models, respectively), M Y ( p ) is the peak absolute magnitude of the SN having p decline-rate, µ 0 is the reddening-free distance modulus in magnitudes, E B ( -V ) MW is the color excess due to interstellar extinction in the Milky Way (MW), ( E B ( -V ) host in the host galaxy), R X,Y are the reddening slopes in filter X or Yand K XY = K XY ( ϕ, z, E ( B -V ) MW , E ( B -V ) host ) is the cross-band K-correction that matches the observed broad-band magnitudes of a redshifted SN taken with filter X to a template SN LC taken in filter Y. Since our SNe have very low redshifts (see Table 1), K-corrections are often negligible compared to the observational uncertainties.

The motivation behind using the EBV2 model, despite the difference between the photometric system of its LC templates and our data, is that the EBV2 model allows the fitting of the s BV color-stretch parameter, which is more suitable for describing fast-decliner (91bg-like) SNe Ia than the canonical ∆ m 15 (Burns et al. 2014) 6 . Also, the EBV2 model is based on a more recent, better photometric calibration than the EBV model.

First, we fit the observed BVRI LCs, adopting R V = 3 1 . (utilizing the calibration=1 mode) for the reddening slope both in the MW and in the host , with χ 2 -minimization using the built-in MCMC routine in SNooPy2 . SNooPy2 has a capability of K-correcting the data while fitting them with templates from a different photometric system, which was also applied in this case. The inferred parameters were the following ones:

- · E B ( -V ) host : interstellar extinction of the host galaxy (in magnitude);
- · T max : moment of the maximum light in the Bband (in MJD);
- · µ 0 : extinction-free distance modulus (in magnitude);

6 We thank the referee for this suggestion.

Table 2. Best-fit parameter from SNooPy2 assuming R V = 3 1 .

| Name       |   E ( B - V ) MW mag | E ( B - V ) host mag   | T max MJD       |     µ | 0 mag   | s BV          |   ∆ m mag | 15      |   χ 2 |
|------------|----------------------|------------------------|-----------------|-------|---------|---------------|-----------|---------|-------|
| SN 2011fe  |               0.0075 | 0.048 (0.060)          | 55815.31 (0.06) | 29.08 | (0.082) | 0.940 (0.031) |     1.175 | (0.008) | 1.058 |
| Gaia16alq  |               0.0576 | 0.089 (0.061)          | 57508.11 (0.14) | 35.05 | (0.084) | 1.132 (0.033) |     0.939 | (0.024) | 1.358 |
| SN 2016asf |               0.1149 | 0.076 (0.097)          | 57464.66 (0.11) | 34.69 | (0.083) | 1.001 (0.031) |     1.102 | (0.023) | 2.722 |
| SN 2016bln |               0.0249 | 0.213 (0.061)          | 57499.40 (0.34) | 34.78 | (0.116) | 1.058 (0.032) |     1.005 | (0.024) | 3.839 |
| SN 2016coj |               0.0163 | 0.000 (0.061)          | 57548.30 (0.34) | 32.08 | (0.083) | 0.788 (0.030) |     1.438 | (0.020) | 1.954 |
| SN 2016eoa |               0.0633 | 0.242 (0.062)          | 57615.66 (0.35) | 34.24 | (0.101) | 0.775 (0.032) |     1.486 | (0.036) | 1.867 |
| SN 2016ffh |               0.0239 | 0.198 (0.061)          | 57630.48 (0.36) | 34.61 | (0.094) | 0.926 (0.032) |     1.082 | (0.018) | 2.392 |
| SN 2016gcl |               0.063  | 0.056 (0.061)          | 57649.48 (0.36) | 35.45 | (0.090) | 1.155 (0.043) |     0.901 | (0.031) | 2.792 |
| SN 2016gou |               0.1095 | 0.258 (0.060)          | 57666.40 (0.34) | 34.25 | (0.104) | 0.984 (0.031) |     0.925 | (0.032) | 6.56  |
| SN 2016ixb |               0.052  | 0.077 (0.061)          | 57745.11 (0.37) | 35.62 | (0.087) | 0.758 (0.032) |     1.657 | (0.039) | 2.072 |
| SN 2017cts |               0.0265 | 0.158 (0.061)          | 57856.83 (0.35) | 34.56 | (0.088) | 0.931 (0.030) |     1.208 | (0.020) | 1.394 |
| SN 2017drh |               0.109  | 1.396 (0.062)          | 57890.94 (0.36) | 32.29 | (0.326) | 0.838 (0.031) |     1.352 | (0.065) | 2.075 |
| SN 2017erp |               0.0928 | 0.210 (0.061)          | 57934.40 (0.35) | 32.34 | (0.097) | 1.174 (0.030) |     1.129 | (0.011) | 1.627 |
| SN 2017fgc |               0.0294 | 0.162 (0.061)          | 57959.77 (0.37) | 32.61 | (0.092) | 1.137 (0.037) |     1.086 | (0.027) | 1.761 |
| SN 2017fms |               0.0568 | 0.022 (0.061)          | 57960.09 (0.35) | 35.52 | (0.084) | 0.746 (0.033) |     1.425 | (0.026) | 1.875 |
| SN 2017hjy |               0.0768 | 0.211 (0.061)          | 58056.02 (0.35) | 34.05 | (0.096) | 0.949 (0.031) |     1.137 | (0.011) | 1.519 |
| SN 2017igf |               0.0456 | 0.158 (0.061)          | 58084.61 (0.36) | 32.8  | (0.096) | 0.608 (0.032) |     1.757 | (0.006) | 3.461 |
| SN 2018oh  |               0.0382 | 0.071 (0.061)          | 58162.96 (0.37) | 33.4  | (0.083) | 1.089 (0.034) |     0.989 | (0.013) | 1.231 |

- · s BV : color-stretch parameter in the EBV2 model;
- · ∆ m 15 : decline rate parameter in the EBV model (in magnitude).

We also attempted to include R V as a fitting parameter, but the results were close the original R V = 3 1 value, . indicating that our photometry is not suitable for constraining this parameter. This is not surprising, given the lack of UV- or NIR-data in our sample.

Note that since the EBV2 model is based on a different photometric system than our Johnson-Cousins BVRI data, the inferred best-fit values for E B ( -V ) host and µ 0 are expected to be slightly different from those obtained from the EBV model that is based on BVRI data. After comparing the best-fit parameters taken from both models, it is found that the E B ( -V ) host values are consistent with each other within 1 σ , while there is a systematic shift between the distance moduli of µ 0 ( EBV ) -µ 0 ( EBV 2) ∼ 0 1 mag. . For consistency, we decided to adopt the E B ( -V ) host and µ 0 parameters from the EBV2 model as final, but added ∼ 0 1 mag . systematic uncertainty to the distance modulus (corresponding to ∼ 5 percent relative error in the distances).

The final best-fit parameters are shown in Table 2. The reported errors include the systematic uncertainties of the template vectors as given by SNooPy2 . Plots of the observed LCs and their best-fit SNooPy2 templates can be found in the Appendix. The overall fitting quality is good; most of the reduced χ 2 values (column 8 in Table 2) are in between 1 and 2, and the highest χ 2 ∼ 6 6 . is that of SN 2016gou. As Figure 12 of the Appendix demonstrates, the best model fits the near-maximum data of SN 2016gou adequately, and the relatively high χ 2 value is likely caused by the scattering in the last few data points after +40 days.

As seen in Table 2, the host reddening of SN 2017drh turned out to be extremely high ( E B ( -V ) host ∼ 1 4 . mag) compared to the rest of the sample. Since this may add an extra uncertainty in the derived distance and other parameters, SN 2017drh was excluded from the sample and was not analyzed further. From the remaining 17 SNe, 8 of them have low E B ( -V ) host values that do not exceed their 2 σ uncertainties ( ∼ 0 12 . mag for most of them), but their 3 σ upper limits still indicate low reddening for them. This suggests that E B ( -V ) host cannot be constrained for these SNe from our data due to the photometric uncertainties. On the other hand, none of the other SNe in our sample show E B ( -V ) host /greaterorsimilar 0 3 mag. . Thus, even if we adopt the upper limits for the 8 low-reddened objects, our sample is only moderately reddened, which significantly re-

Figure 3. Differences between the best-fit SNooPy2 parameters from the R V = 1 5 model and those from the . R V = 3 1 . model. Only the distance modulus (bottom left panel) shows a slight trend against reddening, but is statistically insignificant (see text).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-f9583cf4-98ad-4822-b5b9-844b9a5bb04f.jpg)

duces the extinction-related issues while converting the observed data to physical quantities.

Nevertheless, we also tested the effect of the adopted reddening law on the inferred LC parameters. The motivation for this kind of test was the fact that the validity of the R V = 3 1 . reddening law have been questioned for larger samples of SNe Ia in several studies (e.g. Reindl et al. 2005; Folatelli et al. 2010; Amanullah et al. 2015), but see Scolnic et al. (2014) for a different conclusion. In order to test the bias caused by the assumption of the R V = 3 1 reddening . law on the inferred SNooPy2 parameters, especially on the distance, all fits were re-computed with the built-in calibration = 3 mode (corresponding to R V = 1 5) in . SNooPy2. Since our sample contains only moderately reddened SNe, no major differences in the fit parameters are expected. Indeed, it was found that all the fit parameters from the R V = 1 5 model are consistent . with their counterparts from the R V = 3 1 model within . 1 σ . This is illustrated in Figure 3 where the differences between the best-fit parameters for the R V = 1 5 model . and those of the R V = 3 1 model are plotted against . E B ( -V ) host . Disregarding a single deviant point, the inferred E B ( -V ) host , T max and s BV parameters do not show any trend with the E B ( -V ) host parameter, suggesting negligible dependence of these fit parameters on the assumed reddening law. However, the ∆ µ 0 distance modulus difference ( µ 0 (1 5) . -µ 0 (3 1)) exhibits a . slight trend with reddening as ∆ µ 0 = -0 186( . ± 0 062) . · E B ( -V ) host + 0 008( . ± 0 010), but none of the ∆ . µ 0 residuals exceed their uncertainties more than 1 σ . Thus, we adopt the parameters from the R V = 3 1 model for .

further studies, but we also investigate the effect of the lower reddening slope on the bolometric LC and the final inferred physical parameters (see Sections 3.2, 3.3 and 4).

## 3.2. Construction of the bolometric light curve

While constructing the bolometric LCs we followed the same procedure as applied recently by Li et al. (2019) for SN 2018oh. Briefly, after correcting for extinction within the Milky Way and the host galaxy (Table 1 and 2), the observed BVRI magnitudes were converted to physical fluxes via the calibration of Bessell et al. (1998). Uncertainties of the fluxes were estimated in the same way using the magnitude uncertainties obtained via PSF fitting (Sect. 2). Note that because of the relatively low reddening of the sample, the conversion of broadband magnitudes to quasimonochromatic fluxes is not expected to be significantly affected by the alteration of the spectral energy distribution (SED) due to interstellar extinction.

After that, the fluxes were integrated against wavelength via the trapezoidal rule to get a pseudobolometric ( L BVRI ) LC.

Since the UV- and IR regimes are known to have significant contribution to the true bolometric flux of SNe Ia, they cannot be neglected. Because neither the UV, nor the IR bands were covered by our observations, they were estimated by extrapolations in the following way.

In the UV regime the flux was assumed to decrease linearly between 2000 ˚ and A λ B , and the UV-contribution was estimated from the extinction-corrected B -band flux f B as f UV bol = 0 5 . f B ( λ B -2000) (Li et al. 2019). This estimate is tied directly to the B -band flux, and its validity is tested below.

In the IR, to get an estimate for the contribution of the unobserved flux, a Rayleigh-Jeans tail was fit to the corrected I -band fluxes ( f I ) and integrated between λ I and infinity. This resulted in the formula of f IR bol = 1 3 . f λ / I I 3 that we adopted as the estimate of the contribution of the missing IR-bands. The factor 1.3 was found by matching the IR contribution estimated this way to the one that was obtained via direct integration of real observed near-IR JHK photometry taken from literature for a few well-observed objects (see below).

Finally, the bolometric fluxes are corrected for distances inferred from the distance moduli from the SNooPy2 fits (Section 3.1) as D (Mpc) = 10 0 2 ( . · µ 0 -25)

.

This procedure was validated by comparing the estimated UV and IR-contributions to those calculated from existing data for three well-observed normal Type Ia SNe: 2011fe (Brown et al. 2014; Vink´ o et al. 2012;

Matheson et al. 2012), 2017erp (Brown et al. 2018) and 2018oh (Li et al. 2019). All these data were downloaded directly from the Open Supernova Catalog 7 (Guillochon et al. 2017).

Fig. 4 shows the comparison of the BVRI -extrapolated pseudo-bolometric fluxes (colored symbols) with the ones obtained by direct integrations from the UV to the NIR. In the latter case, the integral of a Rayleigh-Jeans tail was also added to the final bolometric flux, but the tail was fit to the K -band flux ( f K ) instead of f I and the factor 1.3 was not applied.

The top-left panel exhibits the comparison of the flux ratio f BUV to f bol as a function of phase, where f BUV is the integrated flux between the B -band and 2000 ˚ A and f bol is the total bolometric flux. The filled circles are based on extrapolated fluxes (as described above), while the lines are from direct integration using the UV data from the Neil Gehrels Swift Observatory . It is seen that there is reasonable agreement between the extrapolated and the directly integrated blue-UV fluxes: the differences between the f BUV /f bol flux ratios from the extrapolated and the directly integrated data (marked by circles and lines in Fig.4, respectively) do not exceed ± 5 percent for most of the data, which is less than the estimated uncertainty of the total bolometric flux (see below).

The bottom-left panel shows the same flux ratio computed from the interpolated data but for all SNe in our sample. This illustrate that there is no major difference in the flux ratios between the slowand fast-decliner SNe Ia around maximum light, despite a ∼ 5 percent relative flux uncertainty (estimated from the scattering of the data, see above) that we include in the final flux uncertainty estimate.

In the top-right panel the full bolometric LC from direct integration (solid line) and the one based on extrapolation (circles) is shown for the extremely well-observed SN 2011fe. The insets illustrate the same but only for the UV (left) and NIR (right) regimes. Again, there seems to be good agreement between the directly integrated and the extrapolated bolometric fluxes. In the bottom-right panel the calculated luminosity evolution is plotted for all SNe in our sample.

It is concluded that the pseudo-bolometric LCs obtained from extrapolations described above are reliable representations of the true bolometric data, and the systematic errors due to the missing bands do not exceed ∼ 5 percent. Together with the errors due to uncertainties in the distances (see above), the final relative un-

7 https://sne.space

Figure 4. Top left panel: the ratio of BUV flux to the total bolometric flux as a function of time since explosion for SN 2018oh (magenta), SN 2011fe (blue) and SN 2017erp (red). Data obtained by direct integration of Swift fluxes are plotted with lines, while the symbols correspond to data estimated from BVRI observations (see text). Top right panel: Comparison of the pseudobolometric LCs of SN 2011fe derived from the extrapolated BVRI SED (symbols) and by direct integration of the observed UV + optical + NIR data (black line). The two insets show the UV and the IR contributions separately. Bottom left panel: the same as the top left panel but for the whole observed sample. Colors code the different SN subtypes as indicated in the legend. Bottom right panel: the derived pseudo-bolometric light curves for the sample SNe.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-7f3f01dd-3b42-480d-bad9-7094d914c893.jpg)

certainty of the bolometric fluxes are estimated as ∼ 10 percent.

## 3.3. Fitting the bolometric light curve

In order to cross-check the sensitivity of the bolometric LCs to the assumed reddening law, all LCs were recomputed using R V = 1 5 for estimating the extinction . within the host galaxy (cf. Section 3). As expected, this resulted in bolometric LCs having reduced amplitudes with respect to the ones calculated by assuming R V = 3 1. . The amplitude ratio, L (1 5) . /L (3 1), . as a function of E B ( -V ) host turned out to be L (1 5) . /L (3 1) = . 1 -1 349( . ± 0 046) . · E B ( -V ) host , which gives L (1 5) . ≈ 0 65 . · L (3 1) for SN 2016gou that had the highest . E B ( -V ) host ∼ 0 26 mag in our sample (see Table 2). . This directly influences the derived 56 Ni masses, predicting lower values for those SNe that suffered higher reddening within their host (cf. Section 3.3).

We estimated the physical parameters of the SN ejecta via applying the radiation-diffusion model of Arnett (1982) by fitting the bolometric LCs with the Minim code (Chatzopoulos et al. 2013). This is a Monte-Carlo code utilizing the Price-algorithm that intends to find the position of the absolute minimum on the χ 2 -hypersurface within the allowed volume of the parameter space. Parameter uncertainties are estimated from the final distribution of N = 200 test points that probe the parameter space around the χ 2 minimum where ∆ χ 2 ≤ 1. See Chatzopoulos et al. (2013) for more details.

The fitted parameters were the following: the epoch of first light with respect to the moment of the B -band maximum ( t 0 ), the light curve time scale ( t lc ), the gamma-ray leaking time scale ( t γ ), and the initial nickel mass ( M Ni ). These parameters can be found in Table

Table 3. Best-fit and inferred parameters from the bolometric LC fitting using R V = 3 1 .

.

| Name      | t 0           | t lc          | t γ            | M Ni            | κ - (cm   | κ + - 1   | M ej            | v exp - 1    | E kin (10 51 er   |
|-----------|---------------|---------------|----------------|-----------------|-----------|-----------|-----------------|--------------|-------------------|
|           | (day)         | (day)         | (day)          | (M /circledot ) | 2 g )     | 2 g )     | (M /circledot ) | (km s )      |                   |
| 2011fe    | 16.59 (0.061) | 14.87 (0.321) | 37.60 (0.670)  | 0.567 (0.042)   | 0.166     | 0.232     | 1.002 (0.070)   | 11660 (542)  | 0.817 (0.4        |
| Gaia16alq | 19.92 (0.418) | 13.78 (0.899) | 46.669 (0.717) | 0.744 (0.055)   | 0.115     | 0.129     | 1.274 (0.330)   | 10594 (1421) | 0.858 (0.2        |
| 2016asf   | 15.08 (2.579) | 11.35 (1.140) | 39.192 (1.329) | 0.597 (0.149)   | 0.093     | 0.124     | 1.051 (0.430)   | 11459 (2429) | 0.828 (0.4        |
| 2016bln   | 17.42 (0.342) | 14.00 (1.219) | 44.508 (1.125) | 0.789 (0.097)   | 0.124     | 0.146     | 1.211 (0.430)   | 10833 (1964) | 0.852 (0.3        |
| 2016coj   | 14.17 (0.261) | 10.57 (0.367) | 32.967 (0.863) | 0.401 (0.053)   | 0.095     | 0.152     | 0.856 (0.130)   | 12298 (1070) | 0.777 (0.5        |
| 2016eoa   | 14.07 (3.053) | 10.55 (1.092) | 39.038 (0.935) | 0.482 (0.103)   | 0.080     | 0.108     | 1.047 (0.440)   | 11483 (2439) | 0.828 (0.4        |
| 2016ffh   | 14.03 (1.16)  | 9.736 (0.521) | 40.521 (0.926) | 0.573 (0.078)   | 0.067     | 0.088     | 1.035 (0.230)   | 11298 (1316) | 0.836 (0.3        |
| 2016gcl   | 17.79 (2.184) | 15.85 (1.177) | 43.623 (1.182) | 0.689 (0.164)   | 0.162     | 0.195     | 1.185 (0.360)   | 10931 (1728) | 0.849 (0.3        |
| 2016gou   | 15.03 (0.525) | 11.05 (0.798) | 45.793 (1.412) | 0.678 (0.063)   | 0.075     | 0.086     | 1.250 (0.30)    | 10696 (1681) | 0.856 (0.3        |
| 2016ixb   | 15.64 (2.088) | 13.12 (2.364) | 30.520 (1.345) | 0.483 (0.064)   | 0.159     | 0.274     | 0.777 (0.30)    | 12653 (3050) | 0.747 (0.5        |
| 2017cts   | 14.23 (1.101) | 9.97 (1.259)  | 42.485 (0.959) | 0.539 (0.063)   | 0.066     | 0.081     | 1.151 (0.58)    | 11062 (2838) | 0.845 (0.5        |
| 2017erp   | 17.96 (0.092) | 17.40 (0.518) | 37.603 (1.084) | 0.975 (0.083)   | 0.227     | 0.317     | 1.002 (0.13)    | 11661 (966)  | 0.817 (0.4        |
| 2017fgc   | 16.21 (0.239) | 12.50 (0.386) | 45.398 (0.941) | 0.692 (0.047)   | 0.097     | 0.112     | 1.237 (0.16)    | 10734 (799)  | 0.855 (0.2        |
| 2017fms   | 14.04 (0.504) | 10.58 (0.575) | 34.612 (0.731) | 0.360 (0.029)   | 0.091     | 0.138     | 0.909 (0.20)    | 12066 (1407) | 0.794 (0.5        |
| 2017hjy   | 16.29 (0.492) | 12.83 (0.823) | 39.484 (0.840) | 0.688 (0.057)   | 0.117     | 0.156     | 1.060 (0.270)   | 11425 (1545) | 0.830 (0.4        |
| 2017igf   | 19.58 (0.933) | 15.30 (1.353) | 34.554 (1.193) | 0.420 (0.051)   | 0.191     | 0.291     | 0.906 (0.320)   | 12070 (2291) | 0.792 (0.5        |
| 2018oh    | 14.86 (0.864) | 11.17 (1.098) | 44.654 (0.928) | 0.598 (0.059)   | 0.078     | 0.092     | 1.217 (0.480)   | 10824 (2175) | 0.856 (0.3        |

3. The final uncertainty of the nickel mass also contains the error of the distance (as given by SNooPy2 ) which is added to the fitting uncertainty reported by Minim in quadrature.

Plots of the the best-fit bolometric LCs corresponding to the smallest χ 2 are shown in the Appendix.

The crucial parameter in the semi-analytic LC codes is the effective optical opacity ( κ ) that is assumed to be constant both in space and time. To estimate the effective optical opacity for our sample, we applied the same technique as done by Li et al. (2019) for SN 2018oh recently. This technique is based on the combination of the light curve time scale, t lc and that of the gamma-ray leakage, t γ . These parameters can be expressed with the physical parameters of the ejecta as

$$t _ { l c } ^ { 2 } \, = \frac { 2 \kappa M _ { e j } } { \beta c v _ { e x p } } \ \text{ and } \ t _ { \gamma } ^ { 2 } \, = \frac { 3 \kappa _ { \gamma } M _ { e j } } { 4 \pi v _ { e x p } ^ { 2 } } \quad ( 2 ) \quad \ \text{$\quad$\cdot$}$$

(Arnett 1982; Clocchiatti &amp; Wheeler 1997; Valenti et al. 2008; Chatzopoulos et al. 2012; Li et al. 2019), where M ej is the ejecta mass, β = 13 8 is a fixed LC parameter . related to the density distribution (Arnett 1982), v exp is the expansion velocity and κ γ = 0 03 cm g . 2 -1 is the opacity for γ -rays (Wheeler et al. 2015).

Since t lc and t γ are measured quantities, the two formulae in Equation 2 contains three unknowns: M ej , v exp , and κ . Following Li et al. (2019), we apply two additional constraints for M ej and v exp to get upper and lower limits for κ . Assuming that the ejecta mass cannot exceed the Chandrasekhar limit ( M ej ≤ M Ch ) we get a lower limit for the optical opacity, κ -, while assuming a lower limit for v exp as v exp ≥ 10 000 km s , -1 we get an upper limit, κ + (see Equation 2.).

In the first case the lower limit for κ can be calculated as

$$\kappa ^ { - } \ = \sqrt { \frac { 3 \kappa _ { \gamma } t _ { l c } ^ { 4 } \beta ^ { 2 } c ^ { 2 } } { 1 6 \pi t _ { \gamma } ^ { 2 } M _ { C h } } }. \quad \quad ( 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Second, the lower limit for the expansion velocity, v exp = 10 000 km s , -1 , implies

$$\kappa ^ { + } \ = \ \frac { 3 \kappa _ { \gamma } t _ { l c } ^ { 2 } \beta c } { 8 \pi v _ { e x p } t _ { \gamma } ^ { 2 } }.$$

Finally we estimate κ as the average of κ -and κ + , and derive M ej and v exp via the following expressions

$$M _ { e j } \ = \ \frac { 3 \kappa _ { \gamma } t _ { l c } ^ { 4 } \beta ^ { 2 } c ^ { 2 } } { 1 6 \pi t _ { \gamma } ^ { 2 } \kappa ^ { 2 } } \ \text{ and } \ v _ { e x p } \ = \ \frac { 3 \kappa _ { \gamma } t _ { l c } ^ { 2 } \beta c } { 8 \pi \kappa t _ { \gamma } ^ { 2 } }. \ \ ( 5 ) \quad \text{$\text{$As$}$} \\ \quad \text{$the$}$$

Having M ej and v exp evaluated, we express the kinetic energy of ejecta as E kin = 0 3 . · M ej · v 2 exp (Arnett 1982; Chatzopoulos et al. 2012).

The results of these calculations are collected in Table 3, where the best-fit LC parameters for the R V = 3 1 . reddening model are shown. Uncertainties (given in parentheses, as previously) are calculated via error propagation taking into account the uncertainties of the fitted timescales and the mean optical opacity, the latter approximated as σ κ ≈ 0 5( . κ + -κ -).

The same fitting process was repeated on the set of bolometric LCs computed with the R V = 1 5 reddening . law (Section 3.2). Since only the amplitude of the LCs are affected by the choice of the reddening law, and the timescales are not, it is expected that only the 56 Ni values will be different from those listed in Table 3. Indeed, the resulting best-fit parameters remained almost the same, i.e. within their 1 σ uncertainties, except for the 56 Ni masses that were found to be systematically lower. The nickel mass ratio between the two different reddening models had the same dependence on E B ( -V ) host as the luminosity ratio shown in Section 3.2, which is expected since the amplitude of the LC is directly related to the calculated nickel mass.

Note that the lack of spectroscopic (i.e. velocity) data makes the kinetic energies poorly constrained. This can be seen from the large relative errors of the E kin values given in Table 3 that may exceed ∼ 50 percent in some cases. Therefore, we do not use the kinetic energies derived from the LC fitting for constraining the physics in these Type Ia SNe.

Instead, in this paper we focus on constraining the ejecta masses. In order to illustrate that the ejecta mass can be inferred from the combination of t lc and t γ via Eq. 5, we constructed model LCs using the same Arnett model as above, assuming various ejecta parameters: progenitor radius ( R 0 ), ejecta mass ( M ej ), initial nickel mass ( M Ni ), expansion velocity ( v exp ) and optical opacity ( κ ). The left panel of Fig. 5 shows some of them corresponding to κ = 0 1 cm g . 2 -1 , v exp = 11 000 , km s -1 , and M ej in between 0.5 and 2.1 M /circledot . It can be seen that larger M ej implies both longer rise time and slower decline rate consistently with Eq. 2. The right panel indicates the correlation between the ejecta mass and the two timescale parameters, t lc and t γ . Here t γ is plotted with squares while t lc with triangles, and different colors code different physical parameters as follows: blue means v exp = 11 000 km s , -1 , κ = 0 1 cm g . 2 -1 , green denotes v exp = 15 000 km s , -1 , κ = 0 1 cm g . 2 -1 , and orange is v exp = 11 000 km s , -1 , κ = 0 2 cm g . 2 -1 . As it is expected, the shorter the time scales, the lower the model ejecta mass, suggesting that the combination of t lc and t γ outlined above may indeed provide realistic estimates for M ej within the framework of the Arnett model.

Figure 5. Left panel: bolometric LCs for models having M ej in between 0.5 and 2.1 M /circledot . All other input parameters for the models were fixed as R 0 = 0 01 . R /circledot , v exp = 11 000 km s -1 , κ = 0.1 cm g 2 -1 , and M Ni = 0 6 M . /circledot . Right panel: t lc (triangles) and t γ (squares) as a function of M ej . The values corresponding to v exp = 11000 km s -1 and κ = 0.1 cm g 2 -1 are plotted with blue, v exp = 15000 km s -1 and κ = 0 1 cm g . 2 -1 with green, and v exp = 11000 km/s, κ = 0 2 cm g . 2 -1 with orange symbols.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-169c581f-862d-414b-b38f-e2ea59c9a1ef.jpg)

We use two SNe as reference objects in order to test the consistency of our LC modeling described above with those presented in other studies. SN 2011fe is chosen as the first test object due to the availability of precise, high-cadence observations spanning from the near-UV to the near-IR regimes (see Section 3.2). Scalzo et al. (2014) modeled the bolometric LC of SN 2011fe and obtained M ej = 1 19 . ± 0 12 M . /circledot and M Ni = 0 42 . ± 0 08 . M /circledot for the ejected mass and the nickel mass, respectively. Using R V = 3 1, . our best-fit results are M ej = 1 00 . ± 0 070 M . /circledot and M Ni = 0 567 . ± 0 042 M . /circledot (see Table 3). It is seen that these two estimates are only marginally consistent: the difference between the two ejecta masses slightly exceeds 1 σ , while the 56 Ni -masses differ by ∼ 2 σ . In addition to the sensitivity of the 56 Ni -mass to the uncertainties in the distance, this highlights the possible systematic differences between the two modeling schemes applied by Scalzo et al. (2014) and in this paper: Scalzo et al. (2014) used only the late-time bolometric LC to constrain the ejecta and nickel masses via t γ based on the method of Jeffery (1999), while we fit the full Arnett model to the entire LC. Note that the 56 Ni mass of SN 2011fe has also been determined in several other papers, including Pereira et al. (2013) (0 53 . ± 0 11 M . /circledot ), Mazzali et al. (2015) (0 47 . ± 0 07 M . /circledot ) and Zhang et al. (2016) (0 57 . M ). The scattering of these various estimates suggest /circledot a value of M Ni ∼ 0 5 . ± 0 1 M . /circledot for SN 2011fe, which makes both our result and that of Scalzo et al. (2014) consistent.

On the other hand, good agreement is found between the parameters of the other test object, SN 2018oh. Recently Li et al. (2019) derived M ej = 1 27 . ± 0 15 M . /circledot and M Ni = 0 55 . ± 0 04 M . /circledot , which are very similar to our best-fit values (Table 3), M ej = 1 22 . ± 0 48 M . /circledot and M Ni = 0 60 . ± 0 06 M . /circledot . Even though Li et al. (2019) applied the same method for the LC fitting as we use in this paper, their bolometric LC was assembled from a much denser, more extended dataset, including observed near-UV and near-IR photometry. The good agreement between our best-fit parameters and theirs suggests that our parameters are not unrealistic, and probably do not suffer from severe systematic errors. Note, however, that our estimate for the ejecta mass of SN 2018oh has higher uncertainty ( ∼ 40%) than that of Li et al. (2019), thus, our data are less constraining regarding the ejecta mass. This is true for most SNe listed in Table 3.

## 4. DISCUSSION

## 4.1. Early color evolution

In Fig. 6 we plot the early ( B -V ) 0 colors, corrected for both Milky Way and host galaxy reddening (Table 2), of our sample together with the data for several other well-observed SNe Ia: SN 2017cbv (Hosseinzadeh et al. 2017), SN 2011fe (Vink´ et al. o 2012), SN 20112cg (Vink´ et al. 2018), o SN 2017fr (Contreras et al. 2018), SN 2009ig (Marion et al. 2016; Foley et al. 2012), iPTF16abc (Miller et al. 2018), SN 2012ht (Vink´ et al. o 2018) and SN 2013dy (Vink´ o et al. 2018)). In the following we investigate the pre-maximum color evolution of our observed sample.

Early-phase ( B -V ) 0 observations of SNe Ia suggest that they can be divided into two categories: earlyred and early-blue type (e.g. Stritzinger et al. 2018). The cause of this dichotomy is still debated (see Section 1). For example, Miller et al. (2018) proposed phys-

Figure 6. The evolution of the reddening-corrected ( B -V ) 0 colors for the sample SNe (colored symbols). Left panel: the pre-maximum color evolution compared to other well-observed SNe collected from literature (colored lines). Right panel: B -V 0 colors of the sample plotted together with Hsiao template (green curve) up to 90 days after maximum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-3641aa64-29c4-46e6-95ac-d6416f9d1cab.jpg)

ical models for the progenitor system and the explosion of iPTF16abc. This SN Ia showed blue, nearly constant ( B -V ) 0 color starting from t ∼ -10 days, which was thought to be caused by strong 56 Ni-mixing in the ejecta.

SNe Ia experience a dark phase after shock breakout (SB), before the heating from radioactive decay diffuses through the photosphere. The duration of this dark phase depends on how much 56 Ni is mixed into the outer layers of the ejecta. If 56 Ni is confined to the innermost layers, the dark phase lasts for a few days, so that weak mixing leads to redder colors and moderate luminosity rise. On the contrary, strong mixing results in higher luminosities and bluer optical colors. In the latter case the dark phase does not exist or very short, because the γ -photons originating from the Ni-decay rapidly diffuse out from the ejecta. Shappee et al. (2019) found that strong mixing accounts for the early excess light in the LC of SN 2018oh observed by Kepler during the K2-C16 campaign, although its color could not be constrained by the unfiltered K2 observations.

for our sample SNe Ia. Fig.6 illustrates that our data are consistent with the colors of other SNe Ia collected from recent literature (plotted as continuous lines in the left panel of Fig. 6). Also shown is the prediction based on the Hsiao-template (plotted with a green line in the right panel) that represents an empirical description of the time-evolving spectral energy distribution of a fiducial SN Ia (Hsiao et al. 2007). In this case the colors are derived by synthetic photometry using Bessell B and V filter functions (Bessell 1990) on the spectrum templates given by Hsiao et al. (2007).

Another possibility for the early blue flux is the collision of the ejecta into a nearby companion star or some kind of a circumstellar envelope. In this case a strong, quickly declining ultraviolet pulse is thought to be the root cause of the excess blue emission during the earliest phases. However, the observability of this emission requires a favorable geometric configuration, i.e. the companion being in front of the SN toward the observer, so it is expected to occur in less than 10% of the actually observed SNe Ia (e.g. Kasen 2010).

With the use of our new photometric data we attempt to investigate the evolution of the early ( B -V ) 0 color

In order to classify SNe Ia into the early-red and earlyblue groups, photometric data taken between -20 and -10 days before maximum is necessary. Between t = -10 days and t max the ( B -V ) 0 colors are so similar for most SNe Ia that it is almost impossible to make such a distinction. Unfortunately, most of the SNe Ia in our sample do not have photometry taken early enough for this purpose. During our campaign this very early phase ( -20 &lt; t &lt; -10 days) have been observed only in two cases: SN 2016bln and 2017erp (see the left panel in Figure 6).

SN 2016bln is a 1991T-like or slow-decliner Ia belonging to the Type SS (Shallow Silicon) subclass on the Branch-diagram (Cenko et al. 2016). Such SNe Ia are recently found to be associated with the early-blue group by Stritzinger et al. (2018): they pointed out that such SNe tend to be located in between the Core Normal and the Shallow Silicon types on the Branch diagram, similar to 1991T-like events. Their finding suggests that SN 2016bln may also belong to the early-blue group. As seen on the left panel of Figure 6, the earliest data on SN 2016bln are indeed close to those of iPTF16abc

Table 4. Parameters of the DDE and PDDE models by Dessart et al. (2014)

| Model   |   E kin (10 51 erg) |   M Ni (M /circledot ) | Model   | E kin (10 51 erg)   | M N i (M /circledot )   |
|---------|---------------------|------------------------|---------|---------------------|-------------------------|
| DDC0    |               1.573 |                  0.869 | PDDEL1  | 1.398               | 0.758                   |
| DDC6    |               1.53  |                  0.722 | PDDEL3  | 1.353               | 0.685                   |
| DDC10   |               1.52  |                  0.623 | PDDEL7  | 1.336               | 0.604                   |
| DDC15   |               1.465 |                  0.511 | PDDEL4  | 1.344               | 0.529                   |
| DDC17   |               1.459 |                  0.412 | PDDEL9  | 1.342               | 0.408                   |
| DDC20   |               1.442 |                  0.3   | PDDEL11 | 1.236               | 0.312                   |
| DDC22   |               1.345 |                  0.211 | PDDEL12 | 1.262               | 0.268                   |
| DDC25   |               1.185 |                  0.119 |         |                     |                         |

and SN 2017cbv, two well-known members of the earlyblue group, thus, they are consistent with the expected early-color evolution. However, our data are too sparse to draw a more definite conclusion.

The other object in our sample that was observed sufficiently early, SN 2017erp, shows an early ( B -V ) 0 color that is similar to that of SN 2011fe, thus, it belongs to the early-red group. It is interesting that Brown et al. (2018) found that SN 2017erp was a near-UV (NUV)red object, while SN 2011fe was a NUV-blue event. This difference between the NUV-colors seems to be independent from the early-phase optical colors, because 2011fe and 2017erp both had similarly red early ( B -V ) 0 color.

This result may suggest that the observed spread in the early NUV- and optical ( B -V ) 0 colors of SNe Ia has different physical reasons. Brown et al. (2018) concluded that the diversity in the NUV colors is likely due to the metallicity of the progenitor that affects the NUV-continuum and the strength of the Ca H&amp;K features. The fact that SN 2017erp and SN 2011fe have similar early ( B -V ) 0 but different NUV colors may suggest that the early ( B -V ) 0 diversity might not be directly related to the progenitor metallicity.

## 4.2. Comparison with explosion models

In this subsection we compare parameters derived from the bolometric LC fitting, in particular the nickel mass, to those taken from several explosion models. First, we consider the DDE and PDDE models computed by Dessart et al. (2014), as listed in Table 4 (we keep the original naming scheme of these models, such as DDC nn models use the DDE scenario, while the PDDEL nn models assume the PDDE mechanism). It is seen that the Ni-masses of these models span the same range than the ones inferred from the bolometric LC fit-

Table 5. 56 Ni masses of the best-fit explosion models. The second column contains the 56 Ni masses inferred from the bolometric LCs assuming R V = 3 1 reddening law. . Uncertainties are given in parentheses.

| Name      | M bol Ni         | M DDE Ni         | M PDDE Ni        |
|-----------|------------------|------------------|------------------|
|           | ( M /circledot ) | ( M /circledot ) | ( M /circledot ) |
| Gaia16alq | 0.744 (0.055)    | 0.623            | 0.604            |
| 2016asf   | 0.597 (0.149)    | 0.623            | 0.604            |
| 2016bln   | 0.789 (0.097)    | 0.869            | 0.758            |
| 2016coj   | 0.401 (0.053)    | 0.412            | 0.604            |
| 2016eoa   | 0.482 (0.103)    | 0.412            | 0.604            |
| 2016ffh   | 0.573 (0.078)    | 0.511            | 0.529            |
| 2016gcl   | 0.689 (0.164)    | 0.623            | 0.758            |
| 2016gou   | 0.678 (0.063)    | 0.623            | 0.758            |
| 2016ixb   | 0.483 (0.064)    | 0.511            | 0.685            |
| 2017cts   | 0.539 (0.063)    | 0.511            | 0.685            |
| 2017erp   | 0.975 (0.083)    | 0.722            | 0.685            |
| 2017fgc   | 0.692 (0.047)    | 0.623            | 0.604            |
| 2017fms   | 0.360 (0.029)    | 0.511            | 0.685            |
| 2017hjy   | 0.688 (0.057)    | 0.623            | 0.685            |
| 2017igf   | 0.420 (0.051)    | 0.300            | 0.408            |
| 2018oh    | 0.598 (0.059)    | 0.623            | 0.685            |

ting (see Table 3), but the kinetic energies of the models are higher by about a factor of ∼ 2.

PDDE models exhibit strong C II lines shortly after explosion, which are formed in the outer, unburned material. Furthermore, the collision with the previously ejected unbound material surrounding the WD results in the heat-up of the outer layers of the ejecta. Thus, in the PDDE scenario the early color of the SN is bluer, and the luminosity rises faster than in the conventional DDE models.

On the contrary, standard DDE models typically leave no unburned material. Instead, at 1-2 days after explosion they show red optical colors (( B -V ) 0 ∼ 1 mag ), which gets bluer continuously as the SN evolves toward maximum light. After maximum both the DDE and the PDDE scenarios show nearly the same ( B -V ) 0 color.

We compare the observed, reddening-corrected ( B -V ) 0 colors of our sample with the predictions from these explosion models. The colors from the models were derived via synthetic photometry applying the standard Bessell B and V filters (Bessell 1990), as above. Note that the in the redshift range of the observed SNe ( z ≤ 0 031) the . K-corrections for the ( B -V ) 0 color indices does not exceed 0.06 mag, which is comparable

Figure 7. Comparison of the Ni-masses derived from the bolometric LCs (plotted on the horizontal axis) with those from DDE (left panel) and PDDE (right panel) models. Data inferred by assuming R V = 3 1 reddening law are plotted with filled red . symbols, while those from the R V = 1 5 reddening law are shown by open blue ones. . Solid lines represents the expected 1:1 relation on both panels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-86798042-29ee-4ba5-94e3-5cd3457a4c5d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-42c5714f-ca76-40cb-9a5e-0a777816dd11.jpg)

Table 6. The same as Table 5 but assuming R V = 1 5 . reddening law.

| Name      | M bol Ni         | M DDE Ni         | M PDDE Ni        |
|-----------|------------------|------------------|------------------|
|           | ( M /circledot ) | ( M /circledot ) | ( M /circledot ) |
| Gaia16alq | 0.651 (0.055)    | 0.623            | 0.604            |
| 2016asf   | 0.492 (0.149)    | 0.511            | 0.604            |
| 2016bln   | 0.560 (0.097)    | 0.869            | 0.758            |
| 2016coj   | 0.397 (0.053)    | 0.412            | 0.604            |
| 2016eoa   | 0.333 (0.103)    | 0.412            | 0.604            |
| 2016ffh   | 0.410 (0.078)    | 0.511            | 0.604            |
| 2016gcl   | 0.641 (0.164)    | 0.623            | 0.758            |
| 2016gou   | 0.450 (0.063)    | 0.623            | 0.685            |
| 2016ixb   | 0.417 (0.064)    | 0.511            | 0.685            |
| 2017cts   | 0.411 (0.063)    | 0.511            | 0.685            |
| 2017erp   | 0.686 (0.083)    | 0.722            | 0.685            |
| 2017fgc   | 0.538 (0.047)    | 0.623            | 0.604            |
| 2017fms   | 0.338 (0.029)    | 0.511            | 0.685            |
| 2017hjy   | 0.489 (0.057)    | 0.623            | 0.685            |
| 2017igf   | 0.327 (0.051)    | 0.300            | 0.408            |
| 2018oh    | 0.489 (0.059)    | 0.623            | 0.685            |

to the uncertainty of our color measurements in these bands. Thus, the K-corrections were neglected when comparing the observed ( B -V ) 0 colors to the synthetic colors inferred from the the explosion models. A plot comparing the observed and synthetic color curves, for both DDE and PDDE models, can be found in the Appendix.

After computing the synthetic B -V curves as a function of phase, the one that has the lowest χ 2 with respect to the observed ( B -V ) 0 colors was chosen as the most probable explosion model that describes the observed SN. These best-fit models and their nickel masses are shown in Table 5 and 6 for the R V = 3 1 and . R V = 1 5 . reddening laws, respectively. Figure 7 compares the Nimasses of the best-fit models to those derived directly from bolometric LC fitting (Table 3).

Although the grid resolution of the models is inferior, it is seen that the nickel masses from the DDE models nicely correlate with the ones from the bolometric LC fitting; the only outlier, SN 2017erp, may have an overestimated M Ni ∼ 1 M /circledot due to a reddening issue (Brown et al. 2018), even though its early ( B -V ) 0 color was similar to that of SN 2011fe, cf. Section 4.1). The agreement is worse for the PDDE models: most of the best-fit models have practically the same Ni-mass, M Ni ∼ 0.6 - 0.7 M /circledot .

In Figure 8 a similar comparison of the synthetic ( B -V ) colors from different explosion models with observations is shown, but only for the pre-maximum phases. Beside the best-fit DDE and PDDE models by Dessart et al. (2014), synthetic colors from two other theoretical models are also plotted: the N100 explosion model (DDE in a Chandrasekhar-mass WD) by Seitenzahl et al. (2013) and the Violent Merger (VM) model by Pakmor et al. (2012). The color evolution from the empirical Hsiao template (thick grey curve; Hsiao et al. 2007) is also shown for comparison.

It is seen that three of these model families (N100, PDDEL and VM) show similar pre-maximum ( B -V ) colors to the observed data at the earliest ( t &lt; -14

Figure 8. Comparison of the early (dereddened) ( B -V ) 0 colors of SN 2017erp (red circles) and 2016bln (blue triangles) with the synthetic colors predicted by several Ia explosion models (DDE, PDDE and Violent Merger mechanisms, see text). The color evolution from the empirical Hsiao template is also shown with a thick grey line.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-733d9d7a-2169-4725-b0d9-07f566b1a02d.jpg)

days) phases, although the PDDEL models are somewhat redder than the observations after -10 days. The DDC0 model fits the data of SN 2016bln very well. The DDC6 model by Dessart et al. (2014) predicts too red ( B -V ) color at the earliest phases, while those from the Hsiao template look being too blue. Although the models shown here do not seem to constrain the observed early-blue and early-red events, at least the two SNe in our sample (2016bln and 2017erp) that have been sampled at t &lt; -14d have similar ( B -V ) 0 colors to those of the models considered here.

It is concluded that the observed, de-reddened ( B -V ) 0 color evolution of SNe Ia seems to be more-orless reproduced by current models of DDE and/or VM mechanisms. The Ni-masses of the DDE models by Dessart et al. (2014) that match the observed ( B -V ) 0 colors are consistent with the Ni-masses inferred from the bolometric LC fitting assuming R V = 3 1 for the . reddening law, but they are systematically higher than those inferred by applying the R V = 1 5 reddening law. . This is not true for the PDDE models, since they predict too high nickel masses for SNe that have M Ni &lt; 0 6 . M /circledot from their bolometric LC fitting, regardless of the adopted reddening law.

## 4.3. Ejecta parameters

In this section we examine the relations between the inferred ejecta parameters (Table 3) following Scalzo et al. (2014) and Scalzo et al. (2019). First, we use the parameters derived from the R V = 3 1 reddening .

law, then we discuss how the assumption of a different reddening law ( R V = 1 5) affects the conclusions. .

## 4.3.1. Comparison with Scalzo et al. (2014, 2019)

The top left panel of Figure 9 shows the dependence between the SNooPy2 color-stretch parameter, s BV and the ejecta mass. Different colors represent different SN Ia subtypes: normal SNe are plotted with green symbols, while 91T-like events are blue and 91bg-like objects are red.

The dashed line indicates the correlation found by Scalzo et al. (2019) between M ej and s BV : M ej = (1 253 . ± 0 021) + (1 036 . . ± 0 095) . · ( s BV -1). It is seen that their finding is consistent with our data, at least for those that have s BV /lessorsimilar 1.

Fitting a straight line to the observed data resulted in the following empirical relation

$$\underset { n - } { \overset { \iota \sigma } { m } }, \quad M _ { e j } \, = \, ( 1. 1 0 2 \pm 0. 0 8 7 ) + ( 0. 6 6 1 \pm 0. 4 5 5 ) \cdot ( s _ { B V } - 1 ), \, ( 6 )$$

which is also plotted in Figure 9 as a continuous line. It is seen that even though the new fitting parameters predict a less steep correlation between s BV and M ej , they are consistent in ± 1 σ with the ones given by Scalzo et al. (2019). As shown by Figure 9, both the continuous and the dashed lines run within the errorbars of the data in the observed range of the s BV parameter.

The statistical significance of the correlation is investigated via deriving the Pearson correlation coefficient from all data shown in the top left panel of Figure 9. Its value is r = 0 751, suggesting at first glance that . s BV and M ej might be correlated. However, this estimate does not take into account the relatively high errorbars of our data, which likely make the real correlation statistically less significant. In order to test the effect of errors on the correlation coefficient, we computed 5,000 random realizations of the observed data by adding Gaussian random noise (having FWHM values equal to the errors) to each data point, then the Pearson correlation coefficient was computed for each random sample. The final value of the correlation coefficient and its uncertainty was estimated as the average and the standard deviation of the r-values from all of these randomized samples. This resulted in r ′ = 0 301 . ± 0 200, which is . significantly less than the value above that does not take into account the uncertainties. Since the uncertainty of the correlation coefficient is similar to the value of the r ′ parameter itself, it is concluded that the correlation between s BV and M ej is not significant in the present sample due to the relatively high uncertainties of the inferred ejecta masses.

The top right panel of Figure 9 plots M Ni versus s BV . The dashed and continuous lines show the same

Figure 9. Correlation between s BV and M ej (top left panel), s BV and M Ni (top right panel), s BV and t γ (bottom left panel), and M Ni and M ej (bottom right panel). Filled dots correspond to parameters derived from the R V = 3 1 reddening law, while . open triangles denote the results from the R V = 1 5 reddening law. . Normal SNe are plotted with green symbols, while blue and red symbols denote 91T-like and 91bg-like SNe, respectively. Dashed lines in the top panels represent the relations found by Scalzo et al. (2019), while continuous ( R V = 3 1) and dotted ( . R V = 1 5) lines indicate fits to the present data. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-fc919fac-9aab-462e-999f-6cc4660f5a8b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-8503cb7b-77b9-4ec3-bb39-affda6fa34ed.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-f6486e8f-8400-46eb-a354-8d94be413fac.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-6fec36b6-69a4-4382-8700-ef001ae8e1dd.jpg)

correlations as found by Scalzo et al. (2019) and this paper, respectively. The former was given as M Ni = (0 718 . ± 0 027) + (0 903 . . ± 0 108) . · ( s BV -1), while the latter can be expressed as

$$M _ { N i } = ( 0. 6 4 3 \pm 0. 0 2 3 ) + ( 0. 7 6 8 \pm 0. 1 2 2 ) \cdot ( s _ { B V } - 1 ). \ ( 7 ) \quad \text{-k} i$$

The Pearson correlation coefficient, after taking into account the uncertainties as above, is r ′ = 0 745 . ± 0 083. . Its conventional value (from the data alone without the errorbars) would be r = 0 852. . The correlation coefficient suggests that M Ni indeed correlates with s BV , even though the parameters given in Equation 7 are only marginally consistent with those found by Scalzo et al. (2019).

Equation 7 suggests that 91T-like objects with slower decline rate (i.e. higher s BV ) tend to have larger Nimasses, while 91bg-like SNe having lower s BV show smaller M Ni .

The bottom left panel in Figure 9 illustrates the dependence between s BV and t γ that is similar to the one found by Scalzo et al. (2014) between their 'trans- parency time scale' and the color-stretch parameter. The line represents the fit to the data as

$$t _ { \gamma } \ = \ ( 4 1. 0 4 \pm 0. 8 5 ) \ + \ ( 2 1. 7 9 \pm 4. 9 8 ) \cdot ( s _ { B V } - 1 ), \ ( 8 )$$

while the Pearson correlation coefficient is r ′ = 0 722 . ± 0 048 ( , r = 0 749 without the uncertainties). .

The results above suggests that the nickel masses also depend on the ejecta masses. The correlation between the derived M Ni and M ej masses is shown in the bottom right panel in Figure 9. Except one outlier (SN 2017erp that likely has an overestimated M Ni due to its overestimated reddening), these two quantities seem to be connected, as the SNe with higher M Ni tend to have larger M ej as well. However, the Pearson correlation coefficient, after removing the outlying SN 2017erp, but taking into account the uncertainties as above, is only r ′ = 0 280 . ± 0 224, which suggests that the correlation . between the data shown is not statistically significant, even though the correlation coefficient without the errorbars would be r = 0 822. .

Nevertheless, after fitting a straight line to the remaining 16 data points, the result can be expressed as

$$M _ { e j } \ = \ ( 0. 7 2 8 \pm 0. 5 8 6 ) \cdot M _ { N i } \ + \ ( 0. 6 5 1 \pm 0. 3 4 3 ). \ ( 9 )$$

Since the uncertainty of the slope parameter is comparable to its best-fit value itself ( ± 0 586 vs. . 0 728), this . result is also consistent with the low value of the r ′ correlation coefficient determined above. We conclude that although the correlation between M Ni and M ej cannot be ruled out, it cannot be constrained from the present dataset due to mainly the relatively high errorbars of these two inferred parameters.

From Figure 9 it is seen that the parameters derived from the R V = 1 5 reddening law (plotted with open . triangles) do not differ significantly from the ones inferred from the R V = 3 1 model. . The only exception is the 56 Ni-mass (top right panel in Figure 9) that is directly related to the systematic shift of the distance due to the different reddening law. In this case the relation between s BV and M Ni (again, taking the errorbars into account) is found as

$$M _ { N i } \ = \ ( 0. 5 0 6 \pm 0. 0 2 3 ) \ + \ ( 0. 5 5 4 \pm 0. 1 2 5 ) \cdot ( s _ { B V } - 1 ), \quad \begin{smallmatrix} \ u a \\ \text{on} \ \ \\ ( 1 0 ) \ \ \\ \text{str} \ \end{smallmatrix} \.$$

which predicts systematically lower nickel masses for the same s BV compared to either Equation 7 or the relation given by Scalzo et al. (2019).

Similarly, the relation between M ej and M Ni changes to

$$M _ { e j } \ = \ ( 0. 4 7 2 \pm 0. 8 6 8 ) \cdot M _ { N i } \ + \ ( 0. 8 6 6 \pm 0. 4 4 5 ). \ ( 1 1 ) \quad \text{me}$$

As above, the high uncertainty of the slope parameter indicates that the correlation between M Ni and M ej , if exists, is not statistically significant from the present sample, because of the too large uncertainties of the inferred mass parameters. Nevertheless, the trend that higher ejecta masses seem to be coupled with higher nickel masses remains visible in the bottom right panel of Figure 9. A plot published by Scalzo et al. (2014) (see their Fig. 10) also shows a similar increase in M ej for increasing M Ni , even though Scalzo et al. (2014) did not present any analytical fit to their data.

Scalzo et al. (2019) argued that within the framework of the Arnett model the ratio of t lc and t γ ( τ m /t 0 in their nomenclature) should be nearly constant, at least for SNe Ia having M ej &lt; M Ch , i.e. t γ ∼ t lc . Figure 10 shows the dependence between our best-fit t γ and t lc values (Table 3) together with the simple linear relation suggested by Scalzo et al. (2019) (plotted as a dashed line). It is seen that the parameters inferred directly from the bolometric LC fitting do not follow the linear trend proposed by Scalzo et al. (2019). Instead, t γ

Figure 10. The best-fit t lc and t γ values plotted together with the linear relation suggested by Scalzo et al. (2019). Colors and symbol shapes

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-8fd71c59-b919-4962-9160-48bb49a10b3e.jpg)

have the same meaning as in Fig.9.

seems to be nearly independent of t lc . In fact, this finding agrees with the assumption by Scalzo et al. (2019) that the ejecta mass can be reliably estimated from t γ only. It suggests that t γ is a better parameter for constraining the ejecta parameters. However, since t γ is more difficult to measure than t lc , as it requires precise photometry extending up to several months after maximum light, the combination of the two timescales, as shown in this paper, might be useful to reduce the systematic errors that might occur during the inference of the ejecta parameters solely from a single timescale.

## 4.3.2. Comparison with Khatami &amp; Kasen (2018)

Khatami &amp; Kasen (2018) introduced a new analytic relation between peak time and luminosity. They also inferred a new equation connecting the peak time ( t peak ) of the bolometric LC to the diffusion timescale ( t d ) of the Arnett model defined similarly to t lc (see Equation 2). For a centrally located 56 Ni distribution,

$$\frac { t _ { p e a k } } { t _ { d } } \ = \ 0. 1 1 \cdot \ln ( 1 + \frac { 9 t _ { s } } { t _ { d } } ) + 0. 3 6, \quad ( 1 2 )$$

where t d = ( κ · M / v ej ( · c )) 1 / 2 and t s = t Ni = 8 8 day is . the nickel decay time scale.

We can test whether Equation 12 were applicable in the case of our SNe by comparing the best-fit t 0 parameter (i.e. the time between the moment of first light and the epoch of B-band maximum) from the Arnett model (Table 3) to the t peak values inferred from Equation 12. The left panel of Figure 11 shows t 0 as a function of the corresponding t peak values. The solid line indicates the 1:1 relation, which suggests that the t peak values given

Figure 11. Left panel: The LC rise time from the Arnett model ( t 0 in Table 3) versus t peak inferred from Equation 12. Right panel: M Ni inferred from Eq. 13. vs. M Ni from the Arnett model. The black line represents the 1:1 relation in both panels. Colors have the same meaning as in Figure 9.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-0f289fe5-09aa-49f7-b0fe-cdac8cb4a935.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-59eb61a3-d3c2-47e4-8ee8-3f2daa68892d.jpg)

by Equation 12 are consistent with the best-fit t 0 parameters.

Khatami &amp; Kasen (2018) also derived the peak luminosity ( L peak ) as

$$L _ { p e a k } \, = \, \frac { 2 \, \epsilon _ { N i } \cdot M _ { N i } t _ { s } ^ { 2 } } { \beta _ { K } ^ { 2 } t _ { p e a k } ^ { 2 } } [ 1 - ( 1 + \beta _ { K } t _ { p e a k } / t _ { s } ) e ^ { - \beta _ { K } t _ { p e a k } / t _ { s } } ], \quad \, \tan \cdot \\ ( 1 3 ) \quad \, \text{ via } \\ \dots \, \dots \, - \, 9 \, 0 \, \ 1 0 1 0 \, \max \, \leftarrow - 1 \, \leftarrow \, \leftarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \leftarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightoration \, \leftarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightoration \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightoration \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightoration \, \leftarrow \, \leftarrow \, \leftarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightrop \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightoration \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightoration \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightrop \, \rightarrow \, \rightarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightrop \, \rightarrow \, \rightarrow \, \leftarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightarrow \, \rightrop \, \right$$

where /epsilon1 Ni = 3 9 . · 10 10 erg g -1 s -1 is the heating rate of Ni-decay, and β K is a dimensionless parameter, whose value depends on the heating mechanism powering the LC. introduced by Khatami &amp; Kasen (2018) (not related to the β ∼ 13 8 . density distribution parameter in the Arnett model). Khatami &amp; Kasen (2018) showed that β K ∼ 1 for a centrally located heating source (i.e. 56 Ni), while mixing the radioactive Ni toward the outer parts of the ejecta tends to increase the value of β K .

telescope at Piszk´stet˝ station of Konkoly Observatory, e o Hungary. The reduced BVRI LCs were analyzed using the SNooPy2 public LC-fitter code. The reddening of the host galaxy ( E B ( -V ) host ), the moment of the B-band maximum light ( T max ), the extinction-free distance modulus ( µ 0 ), and color-stretch parameter and the LC decline rate ( s BV , ∆ m 15 ) were were obtained via LC-fitting.

In the right panel of Figure 11 we plot the nickel masses calculated via Equation 13 using the observed L peak values from the assembled bolometric light curves and choosing β K = 1, versus the best-fit M Ni parameters from Table 3 assuming R V = 3 1. . It is seen that these two parameters are nicely consistent. Overall, Figure 11 illustrates that the LC timescales and Nimasses inferred by the bolometric LC fits in this paper are in very good agreement with those resulting from the application of the new theoretical relations given by Khatami &amp; Kasen (2018). This supports the validity of using best-fit parameters from our Arnett models for representing the real ejecta parameters in the studied SNe Ia.

## 5. SUMMARY

We presented a photometric investigation of 17 Type Ia supernovae observed with the 0.6/0.9 m Schmidt-

After correcting for the extinction in the Milky Way and the host galaxy, the fluxes of the missing UV- and IR-bands were estimated by extrapolations. The bolometric LCs were constructed applying the trapezodial integration rule, and our methods were tested and validated with NIR and UV data of three well-observed normal Ia SNe: SN 2011fe, SN 2017erp and SN 2018oh. The integrated optical fluxes supplemented by extrapolations into the unobserved UV and IR-bands are found to be reliable representations of the true bolometric data, thus the systematic errors caused by the missing bands should be negligible. Finally, bolometric luminosities were calculated from the integrated bolometric fluxes and the distances derived via SNooPy2 .

We applied the Minim code (Chatzopoulos et al. 2012; Li et al. 2019) to fit the bolometric LCs with the radiation diffusion model of Arnett (1982). The optimized parameters of this model were the moment of first light ( t 0 ), the LC time scale ( t lc ), the γ -ray leakage time scale ( t γ ), and the initial nickel-mass ( M Ni ). It was found that the fitting results are not particularly sensitive to the assumed reddening law within the host galaxy, except for the inferred Ni-masses. Adopting R V = 1 5 instead . of R V = 3 1 in . the hosts resulted in systematically lower Ni-masses depending on the E B ( -V ) host parameter: M Ni (1 5) . /M Ni (3 1) . = 1 -(1 349 . ± 0 046) . · E B ( -V ) host .

One of the critical parameters of the Arnett model is the value of the optical opacity ( κ ), which is approximated as a constant in both space and time. Upper and lower limits for the optical opacity were estimated by using the same method as in Li et al. (2019). This method combines the t lc and t γ parameters, and may give a reliable evaluation of the ejecta mass ( M ej ) and the expansion velocity ( v exp ). As above, comparing the inferred parameters of SN 2011fe and SN 2018oh to those published by Scalzo et al. (2014) and Li et al. (2019), reasonable agreement has been found.

The evolution of the reddening-corrected pre-maximum ( B -V ) 0 color was also studied for the sample SNe in order to decide whether they belong to the early red or early blue group defined recently by Stritzinger et al. (2018). Even though our data are consistent with the expected color evolution of SNe Ia, only two SNe in our sample, SN 2016bln and SN 2017erp, were observed at sufficiently early epochs ( -10 --20 days before maximum) for this purpose. The early ( B -V ) 0 color of SN 2016bln looks like similar to those of the early blue group, which is consistent with its 91T/SS spectral type (Stritzinger et al. 2018). SN 2017erp, a NUV-red object (Brown et al. 2018), however, seems to belong to the early red group, although our data are sparse and have higher uncertainty that prevents drawing a more definite conclusion.

ors and the Ni-masses inferred from the bolometric LCs assuming R V = 3 1 reddening law. . The agreement is worse for the PDDE models, since those models that have synthetic colors most similar to the observed ones have nearly the same Ni-mass, M Ni ∼ 0 65 . ± 0 1 M . /circledot . Also, using the R V = 1 5 reddening law made the agree-. ment between the observed and model Ni-masses worse, since in this case the observed Ni-masses were found to be systematically lower.

Finally, we examined the possible correlations between various physical parameters of the ejecta, such as s BV vs M ej , M Ni , and t γ , as well as M ej vs M Ni . Similar correlations were found as published recently by Scalzo et al. (2019). Our results also turned out to be consistent with the predictions from the new formalism proposed by Khatami &amp; Kasen (2018), again, while assuming R V = 3 1 in the host galaxies. .

The early-phase ( B -V ) 0 colors were also compared to the synthetic colors from DDE, PDDE (Dessart et al. 2014) and other explosion models (e.g. N100 from Noebauer et al. (2017)) in order to test whether the nickel masses in these models were consistent with those derived from the bolometric LC modeling. We found good agreement between the Ni-masses of DDE models whose color curves match the observed ( B -V ) 0 col-

We express our thanks to an anonymous referee whose criticism and helpful suggestions led to a significant improvement of this paper. This work is part of the project 'Transient Astrophysical Objects' GINOP 2.3.2-152016-00033 of the National Research, Development and Innovation Office (NKFIH), Hungary, funded by the European Union. This study has also been partly supported by the Lend¨ ulet LP2018-7/2019 grant of the Hungarian Academy of Sciences, and the K-130405 and K-131508 grants of the Hungarian National Research, Development and Innovation Office. KV is supported by the Bolyai J´nos Research Scholarship of the Huna garian Academy of Sciences. RKT is also supported by the UNKP-19-02 New National Excellence Program of ´ the Ministry for Innovation and Technology.

Software: IRAF (Tody 1986, 1993), SExtractor (Bertin &amp; Arnouts 1996), wcstools (Mink 1997), SNooPy2 (Burns et al. 2011, 2014), Minim (Chatzopoulos et al. 2013)

## REFERENCES

```
                                                                                                                                                                                                        REFERENCIF                                                                                                                                                                                                       
                                                                                                                                                                                                       

                                                                                                                                                                                                        <?                                                                                                                                                                                                        MNRAS, 453, 3300                                                                                                                                                                                                        1                                                                                                                                                                                                        2                                                                                                                                                                                                        3                                                                                                                                                                                                        4                                                                                                                                                                                                        0                                                                                                                                                                                                        6                                                                                                                                                                                                        5                                                                                                                                                                                                        7                                                                                                                                                                                                        8                                                                                                                                                                                                        9                                                                                                                                                                                                        10                                                                                                                                                                                                        11                                                                                                                                                                                                        19
```

Betoule, M., Kessler, R., Guy, J., et al. 2014, A&amp;A, 568,

Brimacombe, J., Stone, G., Post, R., et al. 2017, Transient Name Server Discovery Report 2017-391, 1.

Brown, J. 2016, Transient Name Server Discovery Report

Brown, P. J., Breeveld, A. A., Holland, S., Kuin, P., &amp; Pritchard, T. 2014, Ap&amp;SS, 354, 89

Brown, P. J., Hosseinzadeh, G., Jha, S. W., et al. 2018,

Burns, C. R., Stritzinger, M., Phillips, M. M., et al. 2011,

```
<_Erlang_>
```

| Burns, C. R., Stritzinger, M., Phillips, M. M., et al. 2014, ApJ, 789, 32                              |
|--------------------------------------------------------------------------------------------------------|
| Burns, C. R., Parent, E., Phillips, M. M., et al. 2018, ApJ, 869, 56                                   |
| Cenko, S. B., Cao, Y., Kasliwal, M., et al. 2016, The Astronomer's Telegram, 8909,                     |
| Chatzopoulos, E., Wheeler, J. C., & Vinko, J. 2012, ApJ, 746, 121                                      |
| Chatzopoulos, E., Wheeler, J. C., Vinko, J., Horvath, Z. L., & Nagy, A. 2013, ApJ, 773, 76             |
| Clocchiatti, A., & Wheeler, J. C. 1997, ApJ, 491, 375                                                  |
| Conley, A., Sullivan, M., Hsiao, E. Y., et al. 2008, ApJ, 681, 482                                     |
| Conley, A., Guy, J., Sullivan, M., et al. 2011, ApJS, 192, 1                                           |
| Contreras, C., Phillips, M. M., Burns, C. R., et al. 2018, ApJ, 859, 24                                |
| Cruz, I., Holoien, T. W.-S., Stanek, K. Z., et al. 2016, The Astronomer's Telegram, 8784, 1.           |
| Dessart, L., Blondin, S., Hillier, D. J., & Khokhlov, A. 2014, MNRAS, 441, 532                         |
| Dhawan, S., Leibundgut, B., Spyromilio, J., et al. 2016, A&A, 588, A84                                 |
| Dhawan, S., Leibundgut, B., Spyromilio, J., et al. 2017, A&A, 602, A118                                |
| Dimitriadis, G., Foley, R. J., Rest, A., et al. 2019, ApJL, 870, L1                                    |
| Dhawan, S., Jha, S. W., & Leibundgut, B. 2018, A&A, 609, A72                                           |
| Dhawan, S., Bulla, M., Goobar, A., et al. 2018, MNRAS, 480, 1445                                       |
| Fink, M., R¨pke, o F. K., Hillebrandt, W., et al. 2010, A&A, 514, A53                                  |
| Folatelli, G., Phillips, M. M., Burns, C. R., et al. 2010, AJ, 139, 120                                |
| Foley, R. J., Challis, P. J., Filippenko, A. V., et al. 2012, ApJ, 744, 38                             |
| Gagliano, R., Post, R., Weinberg, E., et al. 2016, Transient Name Server Discovery Report 2016-516, 1. |
| Gagliano, R., Post, R., Weinberg, E., et al. 2017, Transient Name Server Discovery Report 2017-774, 1. |
| Goldstein, D. A., & Kasen, D. 2018, ApJL, 852, L33                                                     |
| Guillochon, J., Parrent, J., Kelley, L. Z., & Margutti, R. 2017, ApJ, 835, 64                          |
| Guy, J., Astier, P., Baumont, S., et al. 2007, A&A, 466, 11                                            |
| Hosseinzadeh, G., Sand, D. J., Valenti, S., et al. 2017, ApJL, 845, L11                                |
| Hsiao, E. Y., Conley, A., Howell, D. A., et al. 2007, ApJ,                                             |
| 663, 1187                                                                                              |
| Iben, I., Jr., & Tutukov, A. V. 1984, ApJS, 54, 335                                                    |

| Itagaki, K. 2017, Transient Name Server Discovery Report 2017-647, 1.                                      |
|------------------------------------------------------------------------------------------------------------|
| Jeffery, D. J. 1999, arXiv e-prints, astro-ph/9907015                                                      |
| Jha, S., Garnavich, P. M., Kirshner, R. P., et al. 1999, ApJS, 125, 73                                     |
| Jones, D. O., Riess, A. G., & Scolnic, D. M. 2015, ApJ, 812, 31                                            |
| Kasen, D. 2010, ApJ, 708, 1025                                                                             |
| Kessler, R., Becker, A. C., Cinabro, D., et al. 2009, ApJS, 185, 32                                        |
| Khatami, D. K., & Kasen, D. N. 2018, arXiv:1812.06522                                                      |
| Kromer, M., Sim, S. A., Fink, M., et al. 2010, ApJ, 719, 1067                                              |
| Levanon, N., & Soker, N. 2019, MNRAS, 486, 5528 Li, Z., Gonzalez, J. E., Yu, H., Zhu, Z.-H., & Alcaniz, J. |
| S. 2016, PhRvD, 93, 043014 Li, W., Wang, X., Vink´, o J., et al. 2019, ApJ, 870, 12                        |
| Maoz, D., Mannucci, F., & Nelemans, G. 2014, ARA&A, 52, 107                                                |
| Marion, G. H., Brown, P. J., Vink´, o J., et al. 2016, ApJ, 820, 92                                        |
| Matheson, T., Joyce, R. R., Allen, L. E., et al. 2012, ApJ, 754, 19                                        |
| Mazzali, P. A., Sullivan, M., Filippenko, A. V., et al. 2015, MNRAS, 450, 2631                             |
| Miller, A. A., Cao, Y., Piro, A. L., et al. 2018, ApJ, 852, 100                                            |
| Miller, A. A., Laher, R., Masci, F., et al. 2016, The Astronomer's Telegram, 8907, 1.                      |
| Systems VI, 249 Noebauer, U. M., Kromer, M., Taubenberger, S., et al.                                      |
| 2017, MNRAS, 472, 2787                                                                                     |
| Nomoto, K., Thielemann, F.-K., & Yokoi, K. 1984, ApJ, 286, 644                                             |
| Pakmor, R., Kromer, M., Taubenberger, S., et al. 2012, ApJL, 747, L10                                      |
| Papadogiannakis, S., Dhawan, S., Morosin, R., et al. 2019, MNRAS, 485, 2343                                |
| Pereira, R., Thomas, R. C., Aldering, G., et al. 2013, A&A, 554, A27                                       |
| Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565                                    |
| Piascik, A. S., & Steele, I. A. 2016, The Astronomer's Telegram, 8991, 1.                                  |
| Prieto, J. L., Rest, A., & Suntzeff, N. B. 2006, ApJ,                                                      |
| 501                                                                                                        |
| 647,                                                                                                       |
| Phillips, M. M. 1993, ApJL, 413, L105                                                                      |

| Pskovskii, I. P. 1977, AZh, 54, 1188                                                       |
|--------------------------------------------------------------------------------------------|
| Reindl, B., Tammann, G. A., Sandage, A., et al. 2005, ApJ, 624, 532                        |
| Rest, A., Scolnic, D., Foley, R. J., et al. 2014, ApJ, 795, 44                             |
| Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009                   |
| Riess, A. G., Macri, L., Casertano, S., et al. 2012, ApJ, 752, 76                          |
| Riess, A. G., Macri, L. M., Hoffmann, S. L., et al. 2016, ApJ, 826, 56                     |
| Riess, A. G., Strolger, L.-G., Casertano, S., et al. 2007, ApJ, 659, 98                    |
| Sand, D. J., Valenti, S., Tartaglia, L., et al. 2017, The Astronomer's Telegram, 10569, 1. |
| Scalzo, R., Aldering, G., Antilogus, P., et al. 2014, MNRAS, 440, 1498                     |
| Scalzo, R. A., Parent, E., Burns, C., et al. 2019, MNRAS, 483, 628                         |
| Scolnic, D., Rest, A., Riess, A., et al. 2014, ApJ, 795, 45                                |
| Seitenzahl, I. R., Ciaraldi-Schoolmann, F., R¨pke, o F. K., et al. 2013, MNRAS, 429, 1156  |
| Shappee, B. J., Holoien, T. W.-S., Drout, M. R., et al. 2019, ApJ, 870, 13                 |
| Sim, S. A., R¨pke, o F. K., Hillebrandt, W., et al. 2010, ApJL, 714, L52                   |
| Sim, S. A., Fink, M., Kromer, M., et al. 2012, MNRAS, 420, 3003                            |
| Stanek, K. Z. 2016, Transient Name Server Discovery Report 2016-1049, 1                    |
| Stanek, K. Z. 2017, Transient Name Server Discovery Report 2017-1273, 1.                   |
| Stanek, K. Z. 2018, Transient Name Server Discovery Report 2018-150, 1.                    |
| Stritzinger, M., Leibundgut, B., Walch, S., et al. 2006, A&A, 450, 241                     |
| Stritzinger, M. D., Shappee, B. J., Piro, A. L., et al. 2018, arXiv:1807.07576             |

| Tody, D. 1986, Proc. SPIE, 733                                                                                |
|---------------------------------------------------------------------------------------------------------------|
| Tody, D. 1993, Astronomical Data Analysis Software and Systems II, 173                                        |
| Tonry, J., Denneau, L., Stalder, B., et al. 2016, Transient Name Server Discovery Report 2016-583, 1.         |
| Tonry, J., Denneau, L., Stalder, B., et al. 2016, Transient Name Server Discovery Report 2016-718, 1.         |
| Tonry, J., Stalder, B., Denneau, L., et al. 2017, Transient Name Server Discovery Report 2017-1123, 1.        |
| Valenti, S., Benetti, S., Cappellaro, E., et al. 2008, MNRAS, 383, 1485                                       |
| Valenti, S., Sand, D. J., & Tartaglia, L. 2017, Transient Name Server Discovery Report 2017-513, 1.           |
| van Rossum, D. R., Kashyap, R., Fisher, R., et al. 2016, ApJ, 827, 128                                        |
| Vink´, o J., S´rneczky, a K., Tak´ts, a K., et al. 2012, A&A, 546, A12                                        |
| Vink´, o J., Ordasi, A., Szalai, T., et al. 2018, PASP, 130, 064101                                           |
| Wang, X., Wang, L., Filippenko, A. V., et al. 2013, Science, 340, 170                                         |
| Wheeler, J. C., Johnson, V., & Clocchiatti, A. 2015, MNRAS, 450, 1295                                         |
| Whelan, J., & Iben, I., Jr. 1973, ApJ, 186, 1007                                                              |
| Wilk, K. D., Hillier, D. J., & Dessart, L. 2018, MNRAS, 474, 3187                                             |
| Wood-Vasey, W. M., Miknaitis, G., Stubbs, C. W., et al. 2007, ApJ, 666, 694                                   |
| Woosley, S. E., & Weaver, T. A. 1994, ApJ, 423, 371 Zhang, K., Wang, X., Zhang, J., et al. 2016, ApJ, 820, 67 |
| Zhang, B. R., Childress, M. J., Davis, T. M., et al. 2017, MNRAS, 471, 2254                                   |
| Zheng, W., Filippenko, A. V., Mauerhan, J., et al. 2017, ApJ, 841, 64.                                        |

## 6. APPENDIX

In the following we present plots of the best-fit SNooPy2 templates to the multi-color data of our SN sample, the best-fit Arnett models to the assembled bolometric LCs ( Appendix) and the comparison between the dereddened ( B -V ) 0 color curves and the synthetic colors computed from the DDE and PDDE models by Dessart et al. (2014).

Figure 12. BVRI LC-fitting with SNooPy2 . The curves corresponding to different filters are shifted vertically for better visibility.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-ad280044-615a-4634-95f5-b6820d176cc8.jpg)

Figure 13. The best-fit bolometric LC-s computed with Minim .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-5b53c15f-4fed-44d9-9446-f96988429f4f.jpg)

Figure 14. Comparison of the observed, de-reddened ( B -V ) 0 colors (filled symbols) with synthetic colors from DDE and PDDE models by Dessart et al. (2014) (colored curves). The black curve corresponds to the synthetic colors inferred from the empirical Hsiao-template (Hsiao et al. 2007).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-cbe23fdb-1b4c-43c3-802c-8b6b5575c883.jpg)

## K¨ onyves-T´th et al. o

Figure 15. The same as Fig.14 but for additional SNe.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00582v2-figure-d2e50ff1-5f61-4ace-b2cc-c5424821dcc7.jpg)
R. Konyves-Toth J. Vinko A. Ordasi K. Sarneczky A. Bodi B. Cseh G. Csornyei Z. Dencs O. Hanyecz B. Ignacz Cs. Kalup L. Kriskovics A. Pal B. Seli A. Sodor R. Szakats P. Szekely E. Varga-Verebelyi K. Vida G. Zsidi
Constraints on the physical properties of Type Ia supernovae from photometry
We present a photometric study of 17 Type Ia supernovae (SNe) based on
multi-color (Bessell BVRI) data taken at Piszkesteto mountain station of
Konkoly Observatory, Hungary between 2016 and 2018. We analyze the light curves
(LCs) using the publicly available LC-fitter SNooPy2 to derive distance and
reddening information. The bolometric LCs are fit with a radiation-diffusion
Arnett-model to get constraints on the physical parameters of the ejecta: the
optical opacity, the ejected mass and the expansion velocity in particular. We
also study the pre-maximum (B-V) color evolution by comparing our data with
standard delayed detonation and pulsational delayed detonation models, and show
that the Ni56 masses of the models that fit the (B-V) colors are consistent
with those derived from the bolometric LC fitting. We find similar correlations
between the ejecta parameters (e.g. ejecta mass, or Ni56 mass vs decline rate)
as published recently by Scalzo et al. (2019).