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

In 2016, the first strongly lensed Type Ia supernova, iPTF16geu at redshift $z=0.409$ with four resolved images arranged symmetrically around the lens galaxy at $z=0.2163$, was discovered. Here, refined observations of iPTF16geu, including the time delay between images, are used to decrease uncertainties in the lens model, including the the slope of the projected surface density of the lens galaxy, $\Sigma\propto r^{1-\eta}$, and to constrain the universal expansion rate $H_0$. Imaging with HST provides an upper limit on the slope $\eta$, in slight tension with the steeper density profiles indicated by imaging with Keck after iPTF16geu had faded, potentially due to dust extinction not corrected for in host galaxy imaging. Since smaller $\eta$ implies larger magnifications, we take advantage of the standard candle nature of Type Ia supernovae constraining the image magnifications, to obtain an independent constraint of the slope. We find that a smooth lens density fails to explain the iPTF16geu fluxes, regardless of the slope, and additional sub-structure lensing is needed. The total probability for the smooth halo model combined with star microlensing to explain the iPTF16geu image fluxes is maximized at $12\,\%$ for $\eta\sim 1.8$, in excellent agreement with Keck high spatial resolution data, and flatter than an isothermal halo. It also agrees perfectly with independent constraints on the slope from lens velocity dispersion measurements. Combining with the observed time delays between the images, we infer a lower bound on the Hubble constant, $H_0 \gtrsim 40\,{\rm km \ s^{-1} Mpc^{-1}}$ at $68.3\,\%$ confidence level.

Colors of Type Ia supernovae in the first few days after explosion provide a potential discriminant between different models. In this paper, we present $g-r$ colors of 65 Type Ia supernovae discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that $g-r$ colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties ($\sigma_\mathrm{noise}\sim\sigma_\mathrm{int}\sim$ 0.18 mag). Colors are nearly constant starting from 6 days after first light ($g-r\sim-0.15$ mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of $\sim-0.25$ mag day$^{-1}$) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of $^{56}$Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers ($M_\mathrm{He}=0.01\,M_\odot$) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by "double-detonation" models with larger helium masses. We finally identify a significant correlation between the early-time $g-r$ slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value=0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on $B-V$ colors.

We report the results from spectroscopic observations of the multiple images of the strongly lensed Type Ia supernova (SN Ia), iPTF16geu, obtained with ground based telescopes and the Hubble Space Telescope (HST). From a single epoch of slitless spectroscopy with HST, we can resolve spectra of individual lensed supernova images for the first time. This allows us to perform an independent measurement of the time-delay between the two brightest images, $\Delta t = 1.4 \pm 5.0$ days, which is consistent with the time-delay measured from the light-curves. We also present measurements of narrow emission and absorption lines characterizing the interstellar medium in the host galaxy at z=0.4087, as well as in the foreground lensing galaxy at z=0.2163. We detect strong Na ID absorption in the host galaxy, indicating that iPTF16geu belongs to a subclass of SNe Ia displaying "anomalously" large Na ID column densities in comparison to the amount of dust extinction derived from their light curves. For the deflecting galaxy, we refine the measurement of the velocity dispersion, $\sigma = 129 \pm 4$ km/s, which significantly constrains the lens model. Since the time-delay between the SN images is negligible, we can use unresolved ground based spectroscopy, boosted by a factor ~70 from lensing magnification, to study the properties of a high-z SN Ia with unprecedented signal-to-noise ratio. The spectral properties of the supernova, such as pseudo-Equivalent widths of several absorption features and velocities of the Si II-line indicate that iPTF16geu, besides being lensed, is a normal SN Ia, indistinguishable from well-studied ones in the local universe, providing support for the use of SNe Ia in precision cosmology. We do not detect any significant deviations of the SN spectral energy distribution from microlensing of the SN photosphere by stars and compact objects in the lensing galaxy.

The observed tension ($\sim 9\%$ difference) between the local distance ladder measurement of the Hubble constant, $H_0$, and its value inferred from the cosmic microwave background (CMB) could hint at new, exotic, cosmological physics. We test the impact of the assumption about the expansion history of the universe ($0.01<z<2.3$) on the local distance ladder estimate of $H_0$. In the fiducial analysis, the Hubble flow Type Ia supernova (SN~Ia) sample is truncated to $z < 0.15$ and the deceleration parameter ($q_0$) fixed to -0.55. We create realistic simulations of the calibrator and Pantheon samples and account for a full systematics covariance between these two sets. We fit several physically motivated dark energy models and derive combined constraints from calibrator and Pantheon SNe~Ia and simultaneously infer $H_0$ and dark energy properties. We find that the assumption on the dark energy model does not significantly change the local distance ladder value of $H_0$, with a maximum difference ($\Delta H_0$) between the inferred value for different models of 0.47 km$^{-1}$ s$^{-1}$ Mpc $^{-1}$, i.e. a 0.6$\%$ shift in $H_0$, significantly smaller than the observed tension. Additional freedom in the dark energy models does not increase the error in the inferred value of $H_0$. Including systematics covariance between the calibrators, low redshift SNe, and high redshift SNe can induce small shifts in the inferred value for $H_0$. The SN~Ia systematics in this study contribute $\lesssim 0.8 \%$ to the total uncertainty on $H_0$.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

There is a strong degeneracy between the luminosity distance ($D_L$) and the observer viewing angle ($\theta_{\rm obs}$; hereafter viewing angle) of the gravitational wave (GW) source with an electromagnetic counterpart, GW170817. Here, for the first time, we present independent constraints on $\theta_{\rm obs} = 32.5 ^{\circ +11.7}_{-9.7}$ from broad-band photometry of the kilonova (kN) AT2017gfo associated with GW170817. These constraints are consistent with independent results presented in the literature using the associated gamma ray burst GRB170817A. Combining the constraints on $\theta_{\rm obs}$ with the GW data, we find an improvement of 24$\%$ on $H_0$. The observer angle constraints are insensitive to other model parameters, e.g. the ejecta mass, half-opening angle of the lanthanide-rich region and the temperature. A broad wavelength coverage extending to the near infrared is helpful to robustly constrain $\theta_{\rm obs}$. While the improvement on $H_0$ presented here is smaller than the one from high angular resolution imaging of the radio counterpart of GW170817, kN observations are significantly more feasible at the typical distances of such events from current and future LIGO-Virgo Collaboration observing runs ($D_L \sim 100$ Mpc). Our results are insensitive to the assumption on the peculiar velocity of the kN host galaxy.

We report lensing magnifications, extinction, and time-delay estimates for the first resolved, multiply-imaged Type Ia supernova iPTF16geu, at $z = 0.409$, using $Hubble\,Space\,Telescope$ ($HST$) observations in combination with supporting ground-based data. Multi-band photometry of the resolved images provides unique information about the differential dimming due to dust in the lensing galaxy. Using $HST$ and Keck AO reference images taken after the SN faded, we obtain a total lensing magnification for iPTF16geu of $\mu = 67.8^{+2.6}_{-2.9}$, accounting for extinction in the host and lensing galaxy. As expected from the symmetry of the system, we measure very short time-delays for the three fainter images with respect to the brightest one: -0.23 $\pm$ 0.99, -1.43 $\pm$ 0.74 and 1.36 $\pm$ 1.07 days. Interestingly, we find large differences between the magnifications of the four supernova images, even after accounting for uncertainties in the extinction corrections: $\Delta m_1 = -3.88^{+0.07}_{-0.06}$, $\Delta m_2 = -2.99^{+0.09}_{-0.08}$, $\Delta m_3 = -2.19^{+0.14}_{-0.15}$ and $\Delta m_4 = -2.40^{+0.14}_{-0.12}$ mag, discrepant with model predictions suggesting similar image brightnesses. A possible explanation for the large differences is gravitational lensing by substructures, micro- or millilensing, in addition to the large scale lens causing the image separations. We find that the inferred magnification is insensitive to the assumptions about the dust properties in the host and lens galaxy.

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

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

A non-local-thermodynamic-equilibrium (NLTE) level population model of the first and second ionisation stages of iron, nickel and cobalt is used to fit a sample of XShooter optical + near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia). From the ratio of the NIR lines to the optical lines limits can be placed on the temperature and density of the emission region. We find a similar evolution of these parameters across our sample. Using the evolution of the Fe II 12$\,$570$\,\mathring{A}\,$to 7$\,$155$\,\mathring{A}\,$line as a prior in fits of spectra covering only the optical wavelengths we show that the 7200$\,\mathring{A}\,$feature is fully explained by [Fe II] and [Ni II] alone. This approach allows us to determine the abundance of Ni II$\,$/$\,$Fe II for a large sample of 130 optical spectra of 58 SNe Ia with uncertainties small enough to distinguish between Chandrasekhar mass (M$_{\text{Ch}}$) and sub-Chandrasekhar mass (sub-M$_{\text{Ch}}$) explosion models. We conclude that the majority (85$\%$) of normal SNe Ia have a Ni/Fe abundance that is in agreement with predictions of sub-M$_{\text{Ch}}$ explosion simulations of $\sim Z_\odot$ progenitors. Only a small fraction (11$\%$) of objects in the sample have a Ni/Fe abundance in agreement with M$_{\text{Ch}}$ explosion models.

An increase in the number of studied Type Ia supernovae (SNe~Ia) has demonstrated that this class of explosions has a greater diversity in its observables than was previously assumed. The reasons (e.g. the explosion mechanism, progenitor system) for such a diversity remain unknown. Here, we analyse a sample of $r$-band light curves of SNe~Ia, focusing on their behaviour $\sim$ 2-4 weeks after maximum light, i.e. the second maximum. We characterise the second maximum by its timing ($t_{r_2}$) and the integrated flux ($\overline{\mathcal{F}}_{r_2}$). We find that $t_{r_2}$ correlates with the "colour-stretch" parameter s$_{BV}$, which can be used as a proxy for $^{56}$Ni mass, and $\overline{\mathcal{F}}_{r_2}$, correlates with the transparency timescale, t$_0$. Using $\overline{\mathcal{F}}_{r_2}$, for a sample of 199 SNe from the Palomar Transient Factory and intermediate Palomar Transient Factory, we evaluate a distribution on t$_0$ for a sample of SNe~Ia found in an "untargeted" survey. Comparing this distribution to the predictions of t$_0$ ranges from models we find that the largest overlap in t$_0$ values between models and observations is for the sub-Chandrasekhar double detonation models. We also compare our relations between t$_0$ and $\overline{\mathcal{F}}_{r_2}$, with that from 1-D explosion models of \citet{GK18} and confirm that $\overline{\mathcal{F}}_{r_2}$, can be used as a diagnostic of the total ejecta mass.

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

We present the best 265 sampled R-band light curves of spectroscopically identified Type Ia supernovae (SNe) from the Palomar Transient Factory (PTF; 2009-2012) survey and the intermediate Palomar Transient Factory (iPTF; 2013-2017). A model-independent light curve template is built from our data-set with the purpose to investigate average properties and diversity in our sample. We searched for multiple populations in the light curve properties using machine learning tools. We also utilised the long history of our light curves, up to 4000 days, to exclude any significant pre- or post- supernova flares. From the shapes of light curves we found the average rise time in the R band to be $16.8^{+0.5}_{-0.6}$ days. Although PTF/iPTF were single-band surveys, by modelling the residuals of the SNe in the Hubble-Lema\^{i}tre diagram, we estimate the average colour excess of our sample to be $<$E$($B$-$V$)> \approx 0.05(2)$ mag and thus the mean corrected peak brightness to be $M_R = -19.02\pm0.02$ $+5 \log( {\rm H}_0 [{\rm km} \cdot{\rm s}^{-1} {\rm Mpc}^{-1}]/70)$ mag with only weakly dependent on light curve shape. The intrinsic scatter is found to be $\sigma_R = 0.186 \pm 0.033$ mag for the redshift range $0.05<z<0.1$, without colour corrections of individual SNe. Our analysis shows that Malmquist bias becomes very significant at z=0.13. A similar limitation is expected for the ongoing Zwicky Transient Facility (ZTF) survey using the same telescope, but new camera expressly designed for ZTF.

We obtained optical and near-infrared spectra of Type$\,$Ia supernovae (SNe$\,$Ia) at epochs ranging from 224 to 496 days after the explosion. The spectra show emission lines from forbidden transitions of singly ionised iron and cobalt atoms. We used non-local thermodynamic equilibrium (NLTE) modelling of the first and second ionisation stages of iron, nickel, and cobalt to fit the spectra using a sampling algorithm allowing us to probe a broad parameter space. We derive velocity shifts, line widths, and abundance ratios for iron and cobalt. The measured line widths and velocity shifts of the singly ionised ions suggest a shared emitting region. Our data are fully compatible with radioactive $^{56}$Ni decay as the origin for cobalt and iron. We compare the measured abundance ratios of iron and cobalt to theoretical predictions of various SN$\,$Ia explosion models. These models include, in addition to $^{56}$Ni, different amounts of $^{57}$Ni and stable $^{54,56}$Fe. We can exclude models that produced only $^{54,56}$Fe or only $^{57}$Ni in addition to $^{56}$Ni. If we consider a model that has $^{56}$Ni, $^{57}$Ni, and $^{54,56}$Fe then our data imply that these ratios are $^{54,56}$Fe / $^{56}$Ni $=0.272\pm0.086$ and $^{57}$Ni / $^{56}$Ni $=0.032\pm0.011$.

We explore cosmological constraints on the sum of the three active neutrino masses $M_{\nu}$ in the context of dynamical dark energy (DDE) models with equation of state (EoS) parametrized as a function of redshift $z$ by $w(z)=w_0+w_a\,z/(1+z)$, and satisfying $w(z)\geq-1$ for all $z$. We perform a Bayesian analysis and show that, within these models, the bounds on $M_{\nu}$ \textit{do not degrade} with respect to those obtained in the $\Lambda$CDM case; in fact the bounds are slightly tighter, despite the enlarged parameter space. We explain our results based on the observation that, for fixed choices of $w_0\,,w_a$ such that $w(z)\geq-1$ (but not $w=-1$ for all $z$), the upper limit on $M_{\nu}$ is tighter than the $\Lambda$CDM limit because of the well-known degeneracy between $w$ and $M_{\nu}$. The Bayesian analysis we have carried out then integrates over the possible values of $w_0$-$w_a$ such that $w(z)\geq-1$, all of which correspond to tighter limits on $M_{\nu}$ than the $\Lambda$CDM limit. We find a 95\% confidence level (C.L.) upper bound of $M_{\nu}<0.13\,\mathrm{eV}$. This bound can be compared with $M_{\nu}<0.16\,\mathrm{eV}$ at 95\%~C.L., obtained within the $\Lambda$CDM model, and $M_{\nu}<0.41\,\mathrm{eV}$ at 95\%~C.L., obtained in a DDE model with arbitrary EoS (which allows values of $w < -1$). Contrary to the results derived for DDE models with arbitrary EoS, we find that a dark energy component with $w(z)\geq-1$ is unable to alleviate the tension between high-redshift observables and direct measurements of the Hubble constant $H_0$. Finally, in light of the results of this analysis, we also discuss the implications for DDE models of a possible determination of the neutrino mass hierarchy by laboratory searches. (abstract abridged)

Key information about the progenitor system and the explosion mechanism of Type Ia supernovae (SNe~Ia) can be obtained from early observations, within a few days from explosion. iPTF16abc was discovered as a young SN~Ia with excellent early time data. Here, we present photometry and spectroscopy of the SN in the nebular phase. A comparison of the early time data with a sample of SNe~Ia shows distinct features, differing from normal SNe~Ia at early phases but similar to normal SNe~Ia at a few weeks after maximum light (i.e. the transitional phase) and well into the nebular phase. The transparency timescales ($t_0$) for this sample of SNe~Ia range between $\sim$ 25 and 41 days indicating a diversity in the ejecta masses. $t_0$ also weakly correlates with the peak bolometric luminosity, consistent with the interpretation that SNe with higher ejecta masses would produce more $^{56}$Ni. Comparing the $t_0$ and the maximum luminosity, L$_{max}$\, distribution of a sample of SNe~Ia to predictions from a wide range of explosion models we find an indication that the sub-Chandrasekhar mass models span the range of observed values. However, the bright end of the distribution can be better explained by Chandrasekhar mass delayed detonation models, hinting at multiple progenitor channels to explain the observed bolometric properties of SNe~Ia. iPTF16abc appears to be consistent with the predictions from the M$_{ch}$ models.

Constraints on models of the late time acceleration of the universe assume the cosmological principle of homogeneity and isotropy on large scales. However, small scale inhomogeneities can alter observational and dynamical relations, affecting the inferred cosmological parameters. For precision constraints on the properties of dark energy, it is important to assess the potential systematic effects arising from these inhomogeneities. In this study, we use the Type Ia supernova magnitude-redshift relation to constrain the inhomogeneities as described by the Dyer-Roeder distance relation and the effect they have on the dark energy equation of state ($w$), together with priors derived from the most recent results of the measurements of the power spectrum of the Cosmic Microwave Background and Baryon Acoustic Oscillations. We find that the parameter describing the inhomogeneities ($\eta$) is correlated with $w$. The best fit values $w = -0.933 \pm 0.065$ and $\eta = 0.61 \pm 0.37$ are consistent with homogeneity at $< 2 \sigma$ level. Assuming homogeneity ($\eta =1$), we find $w = -0.961 \pm 0.055$, indicating only a small change in $w$. For a time-dependent dark energy equation of state, $w_0 = -0.951 \pm 0.112$ and $w_a = 0.059 \pm 0.418$, to be compared with $w_0 = -0.983 \pm 0.127$ and $w_a = 0.07 \pm 0.432$ in the homogeneous case, which is also a very small change. Current data are not sufficient to constrain the fraction of dark matter (DM) in compact objects, $f_p$ at the 95$\%$ C.L., however at 68$\%$ C.L. $f_p < 0.73$. Future supernova surveys will improve the constraints on $\eta$, and $f_p$, by a factor of $\sim$ 10.

The colour evolution of reddened Type Ia supernovae can place strong constraints on the location of dust and help address the question of whether the observed extinction stems from the interstellar medium or from circumstellar material surrounding the progenitor. Here we analyse BV photometry of 48 reddened Type Ia supernovae from the literature and estimate the dust location from their $B-V$ colour evolution. We find a time-variable colour excess $E(B-V)$ for 15 supernovae in our sample and constrain dust to distances between 0.013 and 45 pc ($4\times10^{16}-10^{20}$ cm). For the remaining supernovae, we obtain a constant $E(B-V)$ evolution and place lower limits on the dust distance from the explosion. In all the 48 supernovae, the inferred dust location is compatible with an interstellar origin for the extinction. This is corroborated by the observation that supernovae with relatively nearby dust ($\lesssim$ 1 pc) are located close to the center of their host galaxy, in high-density dusty regions where interactions between the supernova radiation and interstellar clouds close by are likely to occur. For supernovae showing time-variable $E(B-V)$, we identify a potential preference for low $R_\mathrm{V}$ values, unusually strong sodium absorption and blue-shifted and time-variable absorption features. Within the interstellar framework, this brings evidence to a proposed scenario where cloud-cloud collisions induced by the supernova radiation pressure can shift the grain size distribution to smaller values and enhance the abundance of sodium in the gaseous phase.

Early observations of Type Ia supernovae (SNe Ia) provide a unique probe of their progenitor systems and explosion physics. Here we report the intermediate Palomar Transient Factory (iPTF) discovery of an extraordinarily young SN Ia, iPTF 16abc. By fitting a power law to our early light curve, we infer that first light for the SN, that is when the SN could have first been detected by our survey, occurred only $0.15\pm_{0.07}^{0.15}$ days before our first detection. In the $\sim$24 hr after discovery, iPTF 16abc rose by $\sim$2 mag, featuring a near-linear rise in flux for $\gtrsim$3 days. Early spectra show strong C II absorption, which disappears after $\sim$7 days. Unlike the extensivelyobserved SN Ia SN 2011fe, the $(B-V)_0$ colors of iPTF 16abc are blue and nearly constant in the days after explosion. We show that our early observations of iPTF 16abc cannot be explained by either SN shock breakout and the associated, subsequent cooling or the SN ejecta colliding with a stellar companion. Instead, we argue that the early characteristics of iPTF 16abc, including (i) the rapid, near-linear rise, (ii) the nonevolving blue colors, and (iii) the strong C II absorption, are the result of either ejecta interaction with nearby, unbound material or vigorous mixing of radioactive $^{56}$Ni in the SN ejecta, or a combination of the two. In the next few years, dozens of very young \textit{normal} SNe Ia will be discovered, and observations similar to those presented here will constrain the white dwarf explosion mechanism.

We present near infrared (NIR) spectroscopy of the nearby supernova 2014J obtained $\sim$450 d after explosion. We detect the [Ni II] 1.939 $\mu$m line in the spectra indicating the presence of stable $^{58}$Ni in the ejecta. The stable nickel is not centrally concentrated but rather distributed as the iron. The spectra are dominated by forbidden [Fe II] and [Co II] lines. We use lines, in the NIR spectra, arising from the same upper energy levels to place constraints on the extinction from host galaxy dust. We find that that our data are in agreement with the high $A_V$ and low $R_V$ found in earlier studies from data near maximum light. Using a $^{56}$Ni mass prior from near maximum light $\gamma$-ray observations, we find $\sim$0.05 M$_\odot$ of stable nickel to be present in the ejecta. We find that the iron group features are redshifted from the host galaxy rest frame by $\sim$600 km s$^{-1}$.

Observations of high-redshift Type Ia supernovae (SNe~Ia) are used to study the cosmic transparency at optical wavelengths. Assuming a flat $\Lambda$CDM cosmological model based on BAO and CMB results, redshift dependent deviations of SN~Ia distances are used to constrain mechanisms that would dim light. The analysis is based on the most recent Pantheon SN compilation, for which there is a $0.03\pm0.01 {\textrm \,(\rm stat)}$ mag discrepancy in the distant supernova distance moduli relative to the $\Lambda$CDM model anchored by supernovae at $z<0.05$. While there are known systematic uncertainties that combined could explain the observed offset, here we entertain the possibility that the discrepancy may instead be explained by scattering of supernova light in the intergalactic medium (IGM). We focus on two effects: Compton scattering by free electrons and extinction by dust in the IGM. We find that if the discrepancy is due entirely to dimming by dust, the measurements can be modeled with a cosmic dust density $\Omega_{\rm IGM}^{\rm dust} = 8 \cdot 10^{-5} (1+z)^{-1}$, corresponding to an average attenuation of $2\cdot 10^{-5}$ mag Mpc$^{-1}$ in V-band. Forthcoming SN~Ia studies may provide a definitive measurement of the IGM dust properties, while still providing an unbiased estimate of cosmological parameters by introducing additional parameters in the global fits to the observations.

The most precise local measurements of $H_0$ rely on observations of Type Ia supernovae (SNe Ia) coupled with Cepheid distances to SN Ia host galaxies. Recent results have shown tension comparing $H_0$ to the value inferred from CMB observations assuming $\Lambda$CDM, making it important to check for potential systematic uncertainties in either approach. To date, precise local $H_0$ measurements have used SN Ia distances based on optical photometry, with corrections for light curve shape and colour. Here, we analyse SNe Ia as standard candles in the near-infrared (NIR), where intrinsic variations in the supernovae and extinction by dust are both reduced relative to the optical. From a combined fit to 9 nearby calibrator SNe with host Cepheid distances from Riess et al. (2016) and 27 SNe in the Hubble flow, we estimate the absolute peak $J$ magnitude $M_J = -18.524\;\pm\;0.041$ mag and $H_0 = 72.8\;\pm\;1.6$ (statistical) $\pm$ 2.7 (systematic) km s$^{-1}$ Mpc$^{-1}$. The 2.2 $\%$ statistical uncertainty demonstrates that the NIR provides a compelling avenue to measuring SN Ia distances, and for our sample the intrinsic (unmodeled) peak $J$ magnitude scatter is just $\sim$0.10 mag, even without light curve shape or colour corrections. Our results do not vary significantly with different sample selection criteria, though photometric calibration in the NIR may be a dominant systematic uncertainty. Our findings suggest that tension in the competing $H_0$ distance ladders is likely not a result of supernova systematics that could be expected to vary between optical and NIR wavelengths, like dust extinction. We anticipate further improvements in $H_0$ with a larger calibrator sample of SNe Ia with Cepheid distances, more Hubble flow SNe Ia with NIR light curves, and better use of the full NIR photometric data set beyond simply the peak $J$-band magnitude.

The $\Lambda$ Cold Dark Matter model ($\Lambda$CDM) represents the current standard model in cosmology. Within this, there is a tension between the value of the Hubble constant, $H_0$, inferred from local distance indicators and the angular scale of fluctuations in the Cosmic Microwave Background (CMB). We investigate whether the tension is significant enough to warrant new physics in the form of modifying or adding energy components to the standard cosmological model. We find that late time dark energy explanations are slightly disfavoured whereas a pre-CMB decoupling extra dark energy component has a marginally positive Bayesian evidence. A constant equation of state of the additional early energy density is constrained to 0.086$^{+0.04}_{-0.03}$. Although this value deviates significantly from 1/3, valid for dark radiation, the latter is not disfavoured based on the Bayesian evidence. If the tension persists, future estimates of $H_0$ at the 1$\%$ level will be able to decisively determine which of the proposed explanations is favoured.

Recent re-calibration of the Type Ia supernova (SNe~Ia) magnitude-redshift relation combined with cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) data have provided excellent constraints on the standard cosmological model. Here, we examine particular classes of alternative cosmologies, motivated by various physical mechanisms, e.g. scalar fields, modified gravity and phase transitions to test their consistency with observations of SNe~Ia and the ratio of the angular diameter distances from the CMB and BAO. Using a model selection criterion for a relative comparison of the models (the Bayes Factor), we find moderate to strong evidence that the data prefer flat $\Lambda$ CDM over models invoking a thawing behaviour of the quintessence scalar field. However, some exotic models like the growing neutrino mass cosmology and vacuum metamorphosis still present acceptable evidence values. The bimetric gravity model with only the linear interaction term can be ruled out by the combination of SNe~Ia and CMB/BAO datasets whereas the model with linear and quadratic interaction terms has a comparable evidence value to standard $\Lambda$ CDM. Thawing models are found to have significantly poorer evidence compared to flat $\Lambda$ CDM cosmology under the assumption that the CMB compressed likelihood provides an adequate description for these non-standard cosmologies. We also present estimates for constraints from future data and find that geometric probes from oncoming surveys can put severe limits on non-standard cosmological models.

Fast-declining Type Ia supernovae (SN Ia) separate into two categories based on their bolometric and near-infrared (NIR) properties. The peak bolometric luminosity ($\mathrm{L_{max}}$), the phase of the first maximum relative to the optical, the NIR peak luminosity and the occurrence of a second maximum in the NIR distinguish a group of very faint SN Ia. Fast-declining supernovae show a large range of peak bolometric luminosities ($\mathrm{L_{max}}$ differing by up to a factor of $\sim$ 8). All fast-declining SN Ia with $\mathrm{L_{max}} < 0.3 \cdot$ 10$^{43}\mathrm{erg s}^{-1}$ are spectroscopically classified as 91bg-like and show only a single NIR peak. SNe with $\mathrm{L_{max}} > 0.5 \cdot$ 10$^{43}\mathrm{erg s}^{-1}$ appear to smoothly connect to normal SN Ia. The total ejecta mass (M$_{ej}$) values for SNe with enough late time data are $\lesssim$1 $M_{\odot}$, indicating a sub-Chandrasekhar mass progenitor for these SNe.

Accurate standardisation of Type Ia supernovae (SNIa) is instrumental to the usage of SNIa as distance indicators. We analyse a homogeneous sample of 22 low-z SNIa, observed by the Carnegie Supernova Project (CSP) in the optical and near infra-red (NIR). We study the time of the second peak in the NIR band due to re-brightening, t2, as an alternative standardisation parameter of SNIa peak brightness. We use BAHAMAS, a Bayesian hierarchical model for SNIa cosmology, to determine the residual scatter in the Hubble diagram. We find that in the absence of a colour correction, t2 is a better standardisation parameter compared to stretch: t2 has a 1 sigma posterior interval for the Hubble residual scatter of [0.250, 0.257] , compared to [0.280, 0.287] when stretch (x1) alone is used. We demonstrate that when employed together with a colour correction, t2 and stretch lead to similar residual scatter. Using colour, stretch and t2 jointly as standardisation parameters does not result in any further reduction in scatter, suggesting that t2 carries redundant information with respect to stretch and colour. With a much larger SNIa NIR sample at higher redshift in the future, t2 could be a useful quantity to perform robustness checks of the standardisation procedure.

The increase in the number of Type Ia supernovae (SNe\,Ia) has demonstrated that the population shows larger diversity than has been assumed in the past. The reasons (e.g. parent population, explosion mechanism) for this diversity remain largely unknown. We have investigated a sample of SNe\,Ia near-infrared light curves and have correlated the phase of the second maximum with the bolometric peak luminosity. The peak bolometric luminosity is related to the time of the second maximum (relative to the {\it B} light curve maximum) as follows : $L_{max}(10^{43} erg s^{-1}) = (0.039 \pm 0.004) \times t_2(J)(days) + (0.013 \pm 0.106)$. $^{56}$Ni masses can be derived from the peak luminosity based on Arnett's rule, which states that the luminosity at maximum is equal to instantaneous energy generated by the nickel decay. We check this assumption against recent radiative-transfer calculations of Chandrasekhar-mass delayed detonation models and find this assumption is valid to within 10\% in recent radiative-transfer calculations of Chandrasekhar-mass delayed detonation models. The $L_{max}$ vs. $t_2$ relation is applied to a sample of 40 additional SNe\,Ia with significant reddening ($E(B-V) >$ 0.1 mag) and a reddening-free bolometric luminosity function of SNe~Ia is established. The method is tested with the $^{56}$Ni mass measurement from the direct observation of $\gamma-$rays in the heavily absorbed SN 2014J and found to be fully consistent. Super-Chandrasekhar-mass explosions, in particular SN\,2007if, do not follow the relations between peak luminosity and second IR maximum. This may point to an additional energy source contributing at maximum light. The luminosity function of SNe\,Ia is constructed and is shown to be asymmetric with a tail of low-luminosity objects and a rather sharp high-luminosity cutoff, although it might be influenced by selection effects.

Type Ia supernovae have been proposed to be much better distance indicators at near-infrared compared to optical wavelengths -- the effect of dust extinction is expected to be lower and it has been shown that SNe Ia behave more like `standard candles' at NIR wavelengths. To better understand the physical processes behind this increased uniformity, we have studied the $Y$, $J$ and $H$-filter light curves of 91 SNe Ia from the literature. We show that the phases and luminosities of the first maximum in the NIR light curves are extremely uniform for our sample. The phase of the second maximum, the late-phase NIR luminosity and the optical light curve shape are found to be strongly correlated, in particular more luminous SNe Ia reach the second maximum in the NIR filters at a later phase compared to fainter objects. We also find a strong correlation between the phase of the second maximum and the epoch at which the SN enters the Lira law phase in its optical colour curve (epochs $\sim$ 15 to 30 days after $B$ band maximum). The decline rate after the second maximum is very uniform in all NIR filters. We suggest that these observational parameters are linked to the nickel and iron mass in the explosion, providing evidence that the amount of nickel synthesised in the explosion is the dominating factor shaping the optical and NIR appearance of SNe Ia.

A key tracer of the elusive progenitor systems of Type Ia supernovae (SNe Ia) is the detection of narrow blueshifted time-varying Na I D absorption lines, interpreted as evidence of circumstellar material (CSM) surrounding the progenitor system. The origin of this material is controversial, but the simplest explanation is that it results from previous mass loss in a system containing a white dwarf and a non-degenerate companion star. We present new single-epoch intermediate-resolution spectra of 17 low-redshift SNe Ia taken with XShooter on the ESO Very Large Telescope. Combining this sample with events from the literature, we confirm an excess (~20 per cent) of SNe Ia displaying blueshifted narrow Na I D absorption features compared to non-blueshifted Na I D features. The host galaxies of SNe Ia displaying blueshifted absorption profiles are skewed towards later-type galaxies, compared to SNe Ia that show no Na I D absorption, and SNe Ia displaying blueshifted narrow Na I D absorption features have broader light curves. The strength of the Na I D absorption is stronger in SNe Ia displaying blueshifted Na I D absorption features than those without blueshifted features, and the strength of the blueshifted Na I D is correlated with the B-V colour of the SN at maximum light. This strongly suggests the absorbing material is local to the SN. In the context of the progenitor systems of SNe Ia, we discuss the significance of these findings and other recent observational evidence on the nature of SN Ia progenitors. We present a summary that suggests there are at least two distinct populations of normal, cosmologically useful SNe Ia.