R. B. Koshi, M. Doi, S. Sako, K. Maeda, M. Tanaka, N. Takanashi

Attempts to reveal the spectroscopic diversity of Type Ia supernovae (SNe Ia) have led to subclassification schemes such as the Branch system, which classifies SNe Ia into four categories: core normal (CN), broad line (BL), cool (CL), and shallow silicon (SS). The physical origin of these spectroscopic differences, including progenitor channels, explosion mechanisms, or other parameters, however, remains unclear. Moreover, previous work has concentrated primarily on properties near peak luminosity, yielding limited insight into their behavior at later epochs. In this study, we compile $UBVRI$ photometry for 109 SNe Ia and construct the first set of average light curves for each Branch subgroup, spanning from pre-maximum through the late tail. We find pronounced diversity in the $I$-band, especially in the timing of the secondary maximum across subgroups and in the late-time decline of CL events. After correcting for light curve stretch, which reflects the combined influence of ${}^{56}$Ni and ejecta masses, we show that the secondary maximum is powered by Fe II recombination, and its timing is particularly sensitive to the amount of stable iron-group elements (IGEs) synthesized in the explosion. This implies an anti-correlation between the mass of stable IGEs and ${}^{56}$Ni: CL (and possibly BL) events have a larger mass ratio of IGEs/$^{56}$Ni resulting in earlier secondary maxima, while SS events have a smaller ratio and thus later secondary maxima. This trend is naturally explained within the near-$M_{\rm Ch}$ delayed-detonation scenario, whereas it is inconsistent with the positive correlation predicted by the sub-$M_{\rm Ch}$ double-detonation scenario. Finally, we show that stretch-corrected late-time slopes provide a practical diagnostic for CL events, likely linked to an emission feature around $7,200$ Angstroms.