The Multicomponent Fitting to the Star Formation Histories in the TNG Simulation

The star formation history ( SFH ) is a key issue in the evolution of galaxies. In this work, we developed a model based on a Gaussian and gamma function mixture to ﬁ t SFHs with varying numbers of components. Our primary objective was to use this model to reveal the shape of SFHs and the corresponding physical driving factors. Speci ﬁ cally, we applied this model to ﬁ t SFHs from the TNG100-1 simulation. Our study led to the following ﬁ ndings. ( 1 ) Our model ﬁ ts with TNG star formation histories well, especially for high-mass and red galaxies; ( 2 ) a clear relationship exists between the number and shape of ﬁ tted components and the mass and color of galaxies, with notable differences observed between central / isolated and satellite galaxies. ( 3 ) Our model allows us to extract different episodes of star formation within star formation histories with ease and analyze the duration and timing of each star formation episode. Our ﬁ ndings indicate a strong relationship between the timing of each star formation episode and galaxy mass and color.


INTRODUCTION
Star formation is one of the most important processes driving the evolution of galaxies.Well measured star formation history (SFH) can help us understand many aspects in galaxy evolution, such as the underlying physics of boosting and quenching of star formation, and the timing and time scale of these mechanisms (Thomas et al. 2005;Gallazzi et al. 2005;Panter et al. 2007;Choi et al. 2014;Conroy et al. 2014;Pacifici et al. 2016;Carnall et al. 2018;Schreiber et al. 2018).On the other hand, well modeled SFH is important when inferring galaxy properties, such as stellar mass, SFR, metallicity, and dust contents, from the observed spectra or spectral energy distribution (SED) (Conroy 2013;Kauffmann 2014;Janowiecki et al. 2017;Leja et al. 2017;Telles & Melnick 2018;Zhou et al. 2020b).
Parametric models, however, cannot capture all SFHs and may introduce bias if inappropriate functions are utilized (Simha et al. 2014;Carnall et al. 2019).Consequently, researchers have introduced nonparametric models to accommodate the increased complexity and variability of SFH shapes.The simplest nonparametric model employs a piece-wise constant function to fit the  () (Ocvirk et al. 2006;Kelson 2014;Leja et al. 2017;Chauke et al. 2018).The piecewise function method can be enhanced by incorporating adaptive time binning (Tojeiro et al. 2007).Other methods include direct fitting of SFH patterns from theoretical models (Finlator et al. 2007;Pacifici et al. 2012), polynomial expansion of the SFH (Jiménez-López et al. 2022), and Gaussian process regression (Iyer et al. 2019).Although nonparametric methods offer improved flexibility and accuracy, they are computationally demanding and are susceptible to various degeneracies due to the large number of involved free variables (Walcher et al. 2011;Conroy 2013;Lower et al. 2020).
Regardless of the methods employed, the shape of SFHs is always a significant consideration.There hasn't been universal agreement on the best form(s) of SFH models.Furthermore, Carnall et al. (2019) and Leja et al. (2019) have shown that the form of prior significantly impacts the posterior distributions of SFHs.Thus, selecting an appropriate prior, which is based on our comprehension of SFH shapes, is crucial in accurately reconstructing SFHs.
With the emergency of hydrodynamical cosmological simulations (e.g., Illustris (Vogelsberger et al. 2014), Illu-stirsTNG (Nelson et al. 2019;Pillepich et al. 2018a), EA-GLE (Schaye et al. 2015;Crain et al. 2015), FIRE (Wetzel et al. 2016;Hopkins et al. 2018)), it is possible to get numerous realistic and physics motivated SFHs.These simulated SFHs provide references and training sets for studies aimed at recovering SFHs (e.g.Iyer & Gawiser 2017;Iyer et al. 2019).Analyzing simulated SFHs will assist in selecting appropriate priors for observational tasks.Moreover, analyzing the SFHs is necessary to infer the underlying physics of galaxy formation.The simulations encompass numerous processes that play significant roles in galaxy SFHs.These processes include environment-dependent accretion, stochastic variations, minor and major mergers, gas-to-star conversion dictated by physical conditions in the interstellar medium, ejection of gas through galactic winds, and subsequent recycling of ejected material through accretion.Explicit functions can provide insights into the roles and quantities of these physical processes in galaxy formation.For these reasons, the parameterization of SFHs derived from simulations is a topic worthy of investigation (e.g.Simha et al. 2014).
In this study, our approach involves fitting the SFHs from TNG100-1 using a combination of two basis forms: Gaussian distribution and Γ distribution.Our idea aligns closely with the dense basis approach advocated by Iyer & Gawiser (2017), with the distinction that we utilize a reduced number of basis forms.Our goal is to uncover the underlying patterns of SFHs and offer valuable insights for studies focused on recovering SFHs.
The paper is organized as following.In Section 2, we present the dataset employed in this study.In Section 3, we analyze the performance of fitting SFHs.In Section 4, we classify the SFHs based on our fitting results and investigate the relationship between SFH types and galaxy properties.In Section 5, we examine the specific features of SFHs, including peak positions, peak widths, and peak separations.Additionally, we analyze the relationship between SFH characteristics and galaxy properties.Finally, we present our conclusions and engage in additional discussions in Section 6.

DATA AND METHODOLOGY
This section provides a brief introduction to the data used in our study.We explore the star formation history of simulated galaxies using TNG100-1 simulation.

The Simulation and Samples
To conduct our analysis, we select galaxies from the TNG100-1 simulation and reconstruct their SFHs.The TNG100-1 simulation is a large-scale, cosmological simulation which uses the AREPO moving mesh code developed by Springel (2010) and includes the effects of gravity and magnetohydrodynamics.Its cosmological parameters, in agreement with Planck2015 (Collaboration et al. 2016), are Ω Λ = 0.6911, Ω  = 0.3089, Ω  = 0.0486,  8 = 0.8159,   = 0.9667, and ℎ = 0.6774.The simulation covers a box size of 110.7 3 Mpc 3 with 1820 3 dark matter particles and 1820 3 initial hydrodynamic cells.The dark matter particle mass resolution is 7.5 × 10 6  ⊙ and the initial mass resolution of baryons is 1.4 × 10 6  ⊙ .The simulation has 100 output snapshots, covering the redshift range 127 to 0. The TNG100-1 simulation features an updated physical model that encompasses the revised recipes for star formation and evolution, chemical enrichment, cooling, and feedbacks (Weinberger et al. 2017;Pillepich et al. 2018a,b;Nelson et al. 2018a).It also features a modified AGN feedback model to regulate massive galaxies (Weinberger et al. 2017) and a galactic winds model to shape low-mass galaxies (Pillepich et al. 2018a).More details about the simulation and its data release can be found in the introductory paper series of TNG100-1 (Pillepich et al. 2018a;Springel et al. 2018;Nelson et al. 2018a;Naiman et al. 2018;Marinacci et al. 2018) and in the data release itself (Nelson et al. 2019).
A total of 43, 440 galaxies were identified using the SUB-FIND algorithm (Springel et al. 2001).These galaxies satisfy two criteria, with  * ( = 0) > 10 8  ⊙ /ℎ at redshift 0, and a merger tree length constraint such that it is longer than 10 snapshots.We obtained the merger trees of these galaxies from the TNG project supplementary file, which are generated using the SubLink algorithm (Rodriguez-Gomez et al. 2015).
The SFH of a galaxy is determined based on its stellar population.Each stellar particle records its age, measured in terms of lookback time   in units of giga-years (Gyr), as well as its initial stellar mass.We divide their initial stellar masses into different time bins based on their age, using a time bin width of Δ = 0.1 Gyr.The SFH of each galaxy is then calculated as SFR(  ) = ( * , (  − Δ) −  * , (  ))/Δ.This definition of SFH closely resembles those reconstructed through observations, with numerous studies deriving SFR estimates using identical methods (Donnari et al. 2019;Hahn et al. 2019;Matthee & Schaye 2019;Conroy 2013;Johnson et al. 2013).The time bin width of 0.1 Gyr is widely employed in previous works.The chosen of time bin may affect the reconstruction of SFH.We represent the results of tests on different time bin widths in Appendix A. We found that the time bin width of 0.1 Gyr is appropriate for the SFH reconstruction based on the TNG100-1 data.Different from previous studies (Benítez-Llambay et al. 2015;Digby et al. 2019;Joshi et al. 2021), we employ the initial stellar mass rather than the stellar mass at  = 0 to reconstruct the SFH.In the TNG100-1 simulation, the stellar particles return approximately 40% of in their mass to inter stellar medium during their evolution (Pillepich et al. 2018b) .Hence, utilizing the initial stellar mass allows us to account for stellar mass loss and obtain a more accurate estimation of the SFH.Appendix B shows the difference between the SFH derived from the initial stellar mass and the SFH derived from the stellar mass at  = 0.

Fitting Method
We employed a range of combinations of Gaussian and Gamma distributions to establish the best-fit to the SFH curves.The reason we use Gaussian and Gamma distributions as basic functions is because these two patterns have been widely applied in previous works on SFH reconstruction and spectral energy distribution (SED) fitting (Iyer & Gawiser 2017;Zhou et al. 2020a).Moreover, the shapes of these two distributions can respectively fit symmetric and asymmetric distributions, thereby enabling the overall fitting function to have better universality while limiting the number of basis.The general fitting formula was expressed as follows: where   (  ) could be either a Gaussian distribution, or a Gamma distribution, Here,   ,   ,   ,   ,   , and   are fitting parameters.
Also, we utilized the  2 to evaluate the goodness of our fitting through all steps.For a sequence   and its corresponding fitting values    , , the  2 is calculated as: where  represents the mean of the sequence   ,  is the number of data points involved in the fitting, and  is the number of variables, which is set to 1 in this case.The  2 ranges between −∞ and 1, with the results closer to 1 indicating better agreement between data and model.While generally ranging from 0 to 1,  2 may be negative if the fitting is severely poor.The fitting process is given in pseudo-code Algorithm 1.The main idea is to obtain the best fit of each component independently first and then combine them together.Each component is adjusted to match either the local or global curve of the SFH, and both Gaussian and gamma distribution shape of fittings are tested.The best component fit is determined by the maximum  2 value.The overall fit is obtained by summing all the components together.Overall fits with different numbers of maximum components   , ranging from 1 to 6, are created.Then, the best fit is found based on the  2 values of overall fit candidates.We did not test the SFH fittings with   greater than 6 due to the very low fraction of these samples.Our fitting results indicate a significant decrease in the number of SFHs with an increasing number of components.In current set, fraction of SFH fits with 6 components is 4%.The fraction of SFH fits with more than 6 components should make up less than this value.
Figure 1 illustrates the fittings to four galaxies SFHs.We select these four SFHs to cover a range of scenarios with regard to different fitting components and goodness.As it shows, our method is capable of accurately capturing the trends of SFHs, even when the  2 value is low.Furthermore, the bottom right sub-plot reveals that the extremely small  2 value is attributed to the significant fluctuations of the origin SFH.By applying the F-test to determine significance, we find that only 281 SFH fittings fail to reach the 5% significance level (i.e.,  value is greater than 0.05), indicating that only these 281 fittings should be considered unreliable.The  2 values of thess unreliable fittings are from −0.067 to 0.021.However, in order to improve the credibility of following analysis, we further define a group of well-matched samples by selecting SFH fittings that satisfy an arbitrary criterion of  2 ≥ 0.5.There are 23, 229 SFHs remained in this group.In the following sections, we will specify whether we use all samples or well-matched samples for our investigations.
The SFHs of galaxies are classified based on the properties of the fittings.Six categories labeled as "C1" to "C6" are assembled based on the number of fitting components.The SFHs are also divided into three additional categories based on the major shapes of the fitting components.These categories are "G" type SFH with more Gaussian components, "Γ" type SFH with more Gamma components, and "GΓ" type SFH with an equal number of both types of components.Table 1 presents the number of SFHs of each type.
As Table 1 shows, fittings with fewer components are more prevalent than those with multiple components.This is because the number of galaxies increases exponentially as the mass decreases, and smaller mass galaxies dominate in terms of quantity.The star formation process in these smaller mass galaxies is usually shorter, therefore it can generally be described with fewer components.On the other hand, there are also numerical reasons for this, such as our method tending to achieve fitting with the fewest possible components.Table 1 also indicates that SFHs that match our fitting model well contain somewhat more Gaussian components than Gamma components, with the exception of single-component galaxies which exhibit a tendency towards a Gamma-shaped distribution.

GOODNESS OF FIT
Before conducting analysis, it is important to evaluate the ability of our fitting method to accurately describe the SFHs of galaxies.We examined the distribution of the goodness of fit for all samples, which is shown in Figure 2.
For the total trends (grey line in Figure 2), the median of  2 value is 0.53, the peak of  2 distribution is at  2 = 0.70.The median value here is consistent with Iyer & Gawiser (2017).To quantitatively evaluate the performance at reconstructing the SFHs from mock catalogs, e.g., SAMs (Somerville & Davé 2015), hydrodynamic simulation MUFASA (Davé et al. 2017)  Table 1.The numbers of SFHs with different fitting components.The three middle columns show the statistics for all SFHs, while the rightmost columns show the results for SFHs that are well fitted with the Gaussian+Gamma format.Types "C1" to "C6" mean the SFH types classified by component number from 1 to 6. Types "Γ", "GΓ" and "G" mean the SFH types classified by the major component shape, corresponding to SFHs with more Gamma components, with an equal number of two components, or with more Gaussian components.(2017) used  2 , which is in the same definition as this work, to quantify the accuracy of reconstruction.The  2 distribution plot indicates that the ability of reconstruction of our model reaches acceptable accuracy, therefore it can be used to fit the majority of SFHs of galaxies.Nonetheless, the long tail in the  2 < 0 region suggests that there is a small number of SFHs that cannot be fitted with our model.Generally, fittings with multiple components exhibit good accuracy in fitting, where the fitting quality of Gaussian components is slightly superior to that of Gamma components, as shown consistently in the results presented in Table 1.
The shape of the SFH curve of a galaxy is tightly related to its underlying physical processes.In our study, the "goodness of fit" value can indicate to some extent whether the shape of the SFH curve of a galaxy approximates a certain type of function/distribution, specifically a combination of Gaussian and Gamma distributions.Hence, exploring how the  2 value relates to the properties of galaxies at  = 0 will be useful in our investigation of star formation history.  Figure 3 depicts the relationship between the  2 value of SFH fittings and the stellar mass of their corresponding galaxies at  = 0.It is evident from the figure that higher mass galaxies exhibit better fitting quality.This implies that the star formation history of high-mass galaxies conforms more to Gaussian or Gamma distribution shapes.This  2 −  * relation is partially ascribed to the resolution effect present in the simulation.Low-mass galaxies contain a lower number of particles, leading to reduced accuracy and greater noise in reconstructing SFHs.Consequently, the overall precision of the fitting is lowered.Although, we have already set a minimum mass threshold to ensure sufficient galaxy particles to ease the resolution effect.A more comprehensive assessment of this point would necessitate additional simulations with different resolutions, which are beyond the scope of this work.Figure 4 depicts the relationship between the  2 value of SFH fittings and the color of their corresponding galaxies at  = 0.It indicates that the correlation between SFH fitting goodness and color is weak for blue galaxies, where fitting goodness is relatively low.Conversely, at the red end, our model yields better fitting results, as redder colors correspond better with SFH peak functions that resemble Gaussian or Gamma distributions.We attribute this phenomenon to the quenching process that red galaxies underwent during stellar formation history, leading to distinctive SFH with peak- shaped functions that resemble Gaussian or Gamma functions more closely.Figure 5 illustrates the relationship between the  2 values and the age of galaxies.Following Nelson et al. (2018a), the age of a galaxy is defined as the mass weighted average age of all stellar particles within it.It represents the overall distribution of star formation time, as well as a effective length of the SFH Figure 5 exhibits that young and old galaxies have better fitting goodness to our model.This result is straightforward to comprehend.Young galaxies have shorter and simpler shapes, which results from fewer physical mechanisms in the star formation process.Old galaxies corresponds to those massive and red galaxies, which can achieve better fitting results to our model, as have been discussed in previous paragraphs.For both old or young galaxies, their SFHs should have relatively significant peaks at early or late stages, which makes the fittings easier.On the other hand, the mass weighted age reaches half of the SFH length, approximately 6   here, when the SFR distributes evenly accross the whole SFH.In this case, the SFH is hard to be fitted by Gaussian or Gamma distributions, thus results in the extremely low  2 value at  ≃ 6 .In subsequent content, the Figure 9 illustrates that the age distribution of SFHs with more components tends to concentrate around 6  , which also verifies the argument here.
We evaluated the fitting goodness of different galaxy types, as depicted in Figure 6.The galaxies were categorized into two groups: satellite galaxies and central/isolated galaxies.In the simulation, we define the biggest galaxy in a dark matter halo as the central or isolated galaxy, and other galaxies as satellites.The results show that satellite galaxies have a relatively better fitting.We also plot the cumulative density function (CDF) for sub-samples with a stellar mass around 10 9  ⊙ /ℎ (dashed lines) and 10 10  ⊙ /ℎ (dotted lines).The discrepancy between centrals and satellites still exists, negating the possibility of mass effect.A plausible explanation is that satellites are more likely to be quenched in their history, resulting in SFH curves that are similar to Gaussian or Gamma-like peaks.This explanation aligns with the color dependence of  2 .
Generally, only half of the galaxies possess SFH shapes that are consistent with the Gaussian-Gamma combined distribution.High-mass galaxies, red galaxies, and satellite galaxies generally demonstrate a better fitting for these shapes.We hypothesize that these results are mainly due to mechanisms that can shape distinctive peaks in SFHs, such as rapid increases and quenching of star formation.

DEPENDENCY OF SFH TYPES
The SFHs that have been fitted are classified into different types based on the number of components and their major shapes, as explained in Section 2.2.This section aims to investigate how galaxy properties differ depending on the type of SFHs.In order to ensure the validity of our study, we will only investigate well-matched samples of SFHs in this section.This is because it is not convincing to attribute poorly fitted SFHs to specific types as defined by our classification.
Figure 7 shows the stellar mass function of galaxies at  = 0 corresponding to the different types of SFHs.Highmass galaxies tend to have SFHs with more decomposition components, which aligns with our intuition.Concerning the component shapes, high-mass galaxies tend to exhibit a prevalence of Gaussian components, whereas low-mass galaxies tend to be dominated by Gamma components.The first two rows of Figure 11, Figure 12 and Figure 13 show SFH examples of low mass and high galaxies.
Figure 8 shows the distribution of galaxy color at  = 0 corresponding to different types of SFHs.Multiple-component SFHs are predominantly found in blue galaxies, while singlecomponent SFHs are largely associated with red galaxies.The color of galaxies shifts towards the blue area as the number of decomposed components in their history increases.A simple understanding of this is that galaxies with fewer stellar formation epochs are less likely to experience star formation in their later stages, when compared to galaxies with multiple stellar formation epochs.Consequently, the colors of the former tend to be red.Galaxies with more Gamma components tend to have a slightly bluer color than others.This is primarily because the typical pattern of star formation for many blue galaxies is a slow but continuous rise until recent time.
The Gamma distribution provides a good fit for this type of curve.The third and fourth rows of Figure 11, Figure 12 and Figure 13 show SFH examples of red and blue galaxies.
Figure 9 shows the cumulative density distribution of the age of galaxies.The age of a galaxy is defined as the mass weighted age of its all particles.Galaxies with a singlecomponent SFH are concentrated around two area, approximately age of 4 Gyr and 10 Gyr.This corresponds to two types of galaxies.The former represents young galaxies that have just formed and only experienced one period of star for- mation.The latter refers to older quiescent galaxies that have remained in a quiescent state after an initial period of star formation, hence having only one component in their early stage.As the number of components increases, the distribution of galaxy ages gradually converges towards the central position of 6 Gyr.The duration of star formation periods is determined by the underlying physical mechanisms thus will not extend infinitely.Therefore, for SFHs with finite lengths (cannot exceed the age of the universe), there must be more effective star formation times to allow for additional star formation cycles and thus more SFH components.When a galaxy maintains a constant rate of star formation from the beginning of the universe, its average age of stellar particles tends to approach half of the age of the SFH length, approximately 6 Gyr in this  11, Figure 12 and Figure 13 show SFH examples of galaxies with large and small ages.
Figure 10 shows the ratio of the number of central and isolated galaxies to satellite galaxies in galaxies with different SFH types.The horizontal dashed line represents the overall sample ratio.Data points above this line indicate fewer satellite galaxies in this subsample, and vice versa.The proportion of central and isolated galaxies increases with increasing decomposed SFH components, indicating more star formation epochs within central and isolated galaxies.Because most of the satellite galaxies are quenched in late stages, they naturally have fewer star formation epochs.The proportion of satellite galaxies with Gamma SFHs is slightly higher than those with Gaussian SFHs.We believe that it is mainly due to the influence of the cluster stripping, the star formation process in satellite galaxies is quenched more quickly, resulting in the asymmetric gamma distribution of the star formation history curve.
Overall, high-mass, reddish, long-lived galaxies, and central and isolated galaxies tend to have more SFH components, suggesting more star formation periods.Meanwhile, the shapes of the main components of SFHs in these galaxies tend to be slightly Gaussian like.Nevertheless, there is considerable overlap in attribute distributions of galaxies with distinct SFH types, make it a challenging to establish a clear standard for strictly constraining the shape and number of SFH components to apply it in reconstructing the SFH of a given galaxy.

EPISODES OF SFH
Both Gaussian and Gamma distributions have well-defined shapes, hence, the fitting results can be used to characterize the episodes of each star formation history.By exploring the parameters of the fitting function, the peak position and full width at half maximum (FWHM) of each episode can be easily identified and used to further understand the occurrence and time scale of star formation process in galaxies.The position and width of SFH peaks are key parameters in this analysis, as they help us gain insight into the star formation history of galaxies.We continue to constrain our samples with  2 ≥ 0.5 in this section.

Number of episodes
In this section, we recalculated the number of peaks in the SFH to determine the number of star formation episodes present in a given SFH.This is because some fitted components may overlap with others.We combine two fitted components into one peak when the distance between their peaks is smaller than the minimum FWHM of them.The numbers of SFHs with varying numbers of peaks obtained after applying this procedure are presented in the third column of Table 2, with the numbers of SFHs having the same number of components shown in the second column.The table indicates that for each type of multi-component SFH, there are approximately 250 SFHs whose components can be combined, which is quite a small fraction to the total number of SFHs.This suggests that in most cases, the number of fitted components is consistent with the number of star formation episodes in a SFH.
Table 2 indicates that 64.9% of the SFHs have a better fit with two or more peaks.We examined the intrinsic peak numbers of SFHs in the TNG100-1 simulation and confirmed that the reconstructed multi-peaks SFH fraction aligns with the intrinsic fraction.However, these results appear to contradict the findings of Iyer & Gawiser (2017) and Iyer et al. (2019), who reported a fraction of multi-episode SFHs ranging from 10% to 20%.They validated their method using hydro simu-lation and SAMs data, which derived a similar fraction from intrinsic peaks.
The main purpose of Iyer & Gawiser (2017) and Iyer et al. (2019) was to constrain the SFH through SED fitting.It is hard to extract information about all stages in a galaxy's life via SED.The accuracy of fitting the episodes of SFH is highly sensitive to the timing, duration, and strength of the star formation rate, all of which are challenging to reconstruct solely from SED photometry/spectroscopy. Due to the degeneracy of stellar population on spectrum and the noise, the reconstructed SFH from SED is always more smoothedout compared to a direct fit to the original SFH.Another possibility is that the discrepancy arises from the method used to count the intrinsic peaks.However we speculate this possibility.Considering the intrinsic peak number, it is unlikely that the TNG100-1 simulation reproduces SFHs with significantly more peaks than those in the MUFASA simulation and SAMs.We conducted tests using different sets for intrinsic peak counting, and the results are listed in Appendix C. In summary, Iyer & Gawiser (2017) and Iyer et al. (2019) find less peaks, which is reasonable considering the nature of reconstructing SFHs from SED.Take advantages of direct fitting, our approach attempts to distinguish more small peaks compared with the work in Iyer & Gawiser (2017); Iyer et al. (2019).This approach only concerns the performance of direct fitting to SFH.Whether the same performance can be achieved in observations and whether such capability can benefit studies are questions necessitate further exploration and analysis.

Timing and time scale of episodes
Figure 14 shows the distribution of peak positions and widths of SFHs.Each data point in the figure corresponds to one peak.In order to provide readers with a better understanding of the distribution of different intensities of peaks in the multiple-component SFH, we ranked the peaks within the same SFH according to their height and group peaks with same index as sub-samples.Scatters and histograms with different colors in Figure 14 represent the results of peaks with different index.However, we found that the position and width distributions of peaks across different orders are remarkably similar.Therefore, in all subsequent discussions, unless otherwise stated, the subject of the discussion is all peaks.
The distribution of SFH peak position  ,  reveals two conspicuous peaks.The first occurs at approximately  ,  ≈ 0 Gyr, while the second arises at  ,  ≈ 10 Gyr ( ≈ 1.8).The peak occurring at  ≈ 10   is consistent with the general trend of cosmic star formation rate, which found that the star formation rate density peaks at  ≈ 2 (e.g., Behroozi et al. 2013;Vogelsberger et al. 2014;Pillepich et al. 2018b).The peak at  ≈ 0   corresponds to SFH components keep rising.Iyer & Gawiser (2017) likewise documented a significant fraction of recent peaks.
The FWHM of most peaks concentrate in the region of 0   <    < 2  , and reach the maximum of distribution at    ≃ 1.2  .Many works emphasis the importance of SFR regulators in this timescale (e.g., Tacchella et al. 2020;Matthee & Schaye 2019;Katsianis et al. 2020), which correspond to various physical process like gas inflow and outflow, AGN feedbacks, galactic wind, giant molecular cloud life cycles, gas recycling, halo dynamical time-scale, etc. (see Table C1 in Iyer et al. 2020, for a list of time-scale estimates).The peak of FWHM distribution of 1.2   coincides with the time scale of dark matter halo's dynamical timescale (Mo et al. 2010;Lilly et al. 2013).This suggests that the dynamical process of dark matter halos significantly affects the star formation process of galaxies.
The FWHM distribution also includes a part of peaks featuring wider widths that correspond to long-term star formation processes, primarily associated with mergers (Robertson et al. 2006;Jiang et al. 2008;Boylan-Kolchin et al. 2008;Hani et al. 2020), metallicity evolution (Torrey et al. 2018), and galaxy quenching (Sales et al. 2015;Nelson et al. 2018b;Rodríguez Montero et al. 2019;Wright et al. 2019).However, the amount of the wide peaks is 1 to 2 orders of magnitude smaller than that of the narrow peaks.It worth noting that, from the aspect of algorithm, our method is designed to get wider peaks to better represent the whole trend of the SFH of a galaxy.But the FWHM distribution in Figure 14 still prefers narrow peaks.Therefore, it demonstrates the importance of short time-scale processes in shaping the SFH of the galaxy.
The joint distribution of SFH width and position shows three major branches, as indicated by the three thick dashed lines in Figure 14.One of these branches is located vertically near  ,  = 0, which corresponds to galaxies exhibiting a sustained increase in star formation rate until present time.Another branch manifests in the form of a horizontal line with 2   width near    ≃ 1.2  .This branch suggests that a large number of star formation processes occurred within the timescale of 0 to 2  .These processes do not exhibit any clear preference for specific timing of occurring.The third branch occurs like a negative correlation between    and  ,  .For these components, the summation of looking back time of peaks position and    of peak width is close to the length of the SFH they belong to, namely  ,  +    ≈  ℎ .These components are in line with the process of sustained SFR increase from the birth of the galaxy, succeeded by a decline at a specific point in time.In Figure 14, the line    = 13.6 −  ,  is shown as a reference, not meaning the regression for these group of components.
It is important to note that the bimodal distribution of peak positions does not imply the distance between adjacent peaks in a single SFH.Within a single SFH, it is found that adjacent peaks still tend to occur at shorter time intervals.This can be seen in Figure 15.The median peak separation, accompanied by a 1 scatter, is 2.07 +1.76  −0.92  .This value is significantly smaller than the one measured in the study by Iyer et al. (2019).In their study, the peak spacing was determined to be 0.42 +0.15  −0.10   , where   is the age of the universe at the time of observation.Their estimate corresponds to 5.72 +2.04  −1.36   in the context of this work.However, the SFH reconstruction performed by Iyer et al. (2019) used a Gaussian process, which is different from the methodology used in this work.In their study, Iyer et al. (2019) attempted to use the fewest possible number of time points, which resulted in the smoothing of the SFH and the omission of certain shortterm fluctuations.Additionally, as mentioned in the previous section, the SED fitting process naturally leads to fewer peaks being found when reconstructing the SFH.As a consequence, larger intervals between peaks are expected.On the other hand, although similar in methodology, the results from Iyer & Gawiser (2017) are not as convincing for comparison.According to Figure 15 in Iyer & Gawiser (2017), the median peak separation is about 0.07  , which is much smaller than the value found in this work and in Iyer et al. (2019).However, there are only a few tens of SFHs with two episodes in Iyer & Gawiser (2017), which reduces the statistical significance.In summary, comparing peak separation between different studies is challenging due to differences in data and fitting methodologies.These differences warrant further investigation.

Mass dependency
The SFH of a galaxy is strongly associated with the galaxy properties.Consequently, it is reasonable to anticipate the presence of correlations between the extracted SFH features(peak position,    and peaks separation) and the galaxy properties, including stellar mass and color.Figure 16 illustrates the relationship between the stellar mass at  = 0 and the SFH features.For each SFH peak, we record its position,    and the stellar mass of the galaxy it belongs to.Also, we select all neighboring peak pairs from all SFHs and record their separation and the stellar mass of the galaxy they belong to.Then we plot the 2D histograms of these SFH features and the stellar mass of corresponding galaxy.
The peak position of SFH in a galaxy exhibits a U-shaped relationship with the galaxy's mass.The peak of SFH occurs earlier in galaxies with high and low masses, while it occurs later in galaxies with moderate masses.Generally, the formation of high-mass galaxies requires larger dark matter halos with deeper potential wells to attract more gas.The gas accretion into these strong potential wells triggers star formation earlier.As the galaxy evolves, mechanisms such as AGN feedback and metal enrichment gradually strengthen to suppress the process of stellar formation.This suppression effect is more pronounced in higher mass galaxies, preventing their star formation rate from surpassing the early levels in subsequent times.The combination of these mechanisms leads to the earlier appearance of the SFH peaks in higher mass galaxies.On the other hand, for galaxies at the lowest mass end, their limited gas reservoirs result in an early cessation of stellar formation, leading to smaller final masses.Consequently, the peak position of their SFH remains in the early stages.As for satellite galaxies, their stellar formation process stops earlier due to accretion and gas stripping from host galaxies, causing the peak position of SFH in moderate mass galaxies to shift towards earlier times The correlation between the width of SFH peak and galaxy mass is relatively weak.With the exception of galaxies at low and high mass ends, there is a slight decrease in the    as galaxy mass increases.Analyzing the joint distribution of FWHM and mass, we observe a reduction in the number of broader peaks with the mass range of approximately 9 ≲  10 ( * /(ℎ −1  ⊙ )) ≲ 11.This suggests a decreasing influence of the long time-scale star formation driving mechanisms on high-mass galaxies, implying that short time-scale process like feedbacks and gas recycling become more active in high-mass galaxies.Both at the high and low mass ends, the median    experiences an increase, primarily due to the decline in the number of most narrow peaks in the SFHs.The SFH of massive galaxies involves a large number of physical processes, resulting in a more complex SFH shape and more decomposable components.As a result, our method tends to exclude short time-scale components when there are too many pronounced components to be extracted.On the other hand, the evolution process in low-mass galaxies is less complex with fewer fluctuations, leading to fewer instances of fitting with narrow peaks.
Excluding the highest and lowest mass galaxies, the peak separation of SFH decreases slightly with increasing galaxy mass.Based on the context of galaxy SFR rejuvenation (Fang et al. 2012), it can be inferred that larger the galaxy mass results in a shorter time required to return to a state of star formation.This is attributed to high-mass galaxies having stronger gas accretion rates and a higher chance of merging, facilitating the acquisition of necessary fuel for the star formation process.When distinguishing central/isolated galaxies from satellite galaxies, the correlation between peak separation and galaxy mass is weaker for satellite galaxies.We propose two reasons for this.On one hand, when satellite galaxies are influenced by central galaxies and experience gas stripping, the rise in their star formation rate is halted earlier, causing an earlier occurrence of the peak position in the star formation process.This suppressing effect is more pronounced in low-mass galaxies.Therefore, the larger peak separation of SFH at low mass is reduced more, resulting in a flatter relationship curve.On the other hand, the frequency of disturbance and stripping events for satellite galaxies is primarily dependent on the external environment rather than their own properties.This introduces more randomness and The colors correspond to logarithmic galaxy count in the corresponding SFH-properties-stellar-mass bins.The median values of the peak position, FWHM, or peak separation in different stellar mass bins are displayed by the gray lines.The error bars denote the range of the same value from 20% to 80% within the bin.The statistics for all galaxies are presented in the left row, those for central and isolated galaxies are in the middle row, while satellite galaxies are on the right row.
disturbance to their star formation process, thereby weakening the relationship between their SFH characteristics and their own stellar mass.

Color dependency
Figure 17 illustrates the relationship between the galaxy color at  = 0 and the SFH features.For each SFH peak, we record its position,    and the  −  color of the galaxy it belongs to.Also, we select all neighboring peak pairs from all SFHs and record their separation and the  −  color of the galaxy they belong to.Then we plot the 2D histograms of these SFH features(peak position,    and peaks separation) and the color of corresponding galaxies.
According to Figure 14, there exists a clear linear correlation between the color of galaxies and the position of the SFH peak.It indicates that galaxies with a redder color tend to have earlier peaks, which is naturally attributed to the fact that galaxies formed in earlier tend to have older stellar ages, thus displaying a redder color.
The width of the SFH peaks exhibits a tri-modal distribution with the colors of galaxies.At the blue end with  −  < 0.2, the peaks have smaller FWHM.This is because the most blue galaxies are still in the early stage of star formation, with their SFR rapidly increasing and their peak positions mainly located at a redshift of 0. Therefore, the fitted peaks have narrower widths.When the  −  value is within the range of 0.2 to 0.6, the width of the peaks has no correlation with the colors of the galaxies.When the galaxy color exceeds 0.6, a large number of wide peaks appear in satellite galaxies.These are SFH components in satellite galaxies whose star formation processes are regulated by self quenching.This result is consistent with previous work which summarized the The statistics for all galaxies are presented in the left row, those for central and isolated galaxies are in the middle row, while satellite galaxies are on the right row.
quenching time scale in other simulations.Sales et al. (2015) found that the quenching time scale of satellite galaxies is is about 2 to 5 Gyr in Illustris.Wright et al. (2019) found that the quenching time scale is about 2.5 ∼ 3.3 Gyr, and extending out to 10 Gyr, in EAGLE.The peaks separation in the SFH also has a clear relationship with galaxy color.In redder galaxies, the stellar formation process is concentrated in earlier periods, resulting in a more concentrated peak position of star formation.

CONCLUSION
To summarize, this study has undertaken the following objectives and obtained statistical findings concerning the SFHs of galaxies from the TNG100-1 simulation: • We employed a function combined by up to 6 components, described by Gaussian or Gamma distributions, to construct fittings for the SFHs in the TNG100-1 simulation.Our fitting method can well recover the trends of the SFH, though the goodness of fit can depend on galaxy mass, color, history length, and the type of galaxy.
• We analyzed the SFHs that are well-fitted by the Gaussian and Gamma distribution combinations.Our finding indicates that low-mass galaxies, red galaxies, and satellite galaxies tend to have fewer components and to be majored by Gamma distribution, while high-mass galaxies, blue galaxies and central galaxies tend to have more components and to be majored by Gaussian distribution.
• We calculated the distribution of the position, width, and separation of SFH peaks.The peak position distri-bution is concentrated at   = 0   and   = 10  .
The width distribution of peaks is mainly concentrated at    ∼ 0   − 2  .The joint distribution of peak position and width is concentrated in three regions, corresponding to different types of star formation modes.In addition, we found that the separation between adjacent peaks in the TNG simulation is concentrated around 2  , much smaller than previously reconstructed SFHs from observations.
• We analyzed the relationships between three temporal attributes of the SFH, namely peak position, peak width, and peak spacing, and two properties of galaxies, namely stellar mass and galaxy color.Our findings indicate that the peak position and peaks separation are related with stellar mass and galaxy color.The peak width shows limited relationship with stellar mass and galaxy color.This is reasonable because the peak width reflects the time scale of basic physical processes which drives star formation.
Our model can accurately reproduce the trend of galaxies' stellar formation evolution.The method of multi-component decomposition can greatly help us understand the different physical driving factors of star formation in galaxies.Although more than half of the SFHs have low goodness of fit, this is mainly due to the excessive fluctuations in the SFHs, which makes it extremely difficult to fit them for all kinds of functional templates.Of course, the method we used is based on simulation data, the gaps between our method and SFH reconstruction methods suit for observational galaxies can not be neglected.Since it is difficult to give an absolute resolution criterion for the shape types and number of fitted components, our model will require better algorithms and more advanced computing power to determine these factors when applied to observational data, which will be further developed and improved in future work.

APPENDIX A. TIME BIN RESOLUTION IN SFH RECONSTRUCTION
The resolution of the time bin significantly influences the extraction of a galaxy's star formation history and the subsequent fitting work.Figure 18 illustrates the SFHs (represented by the blue line) and the corresponding fitting curves (represented by the orange line) obtained using different time bin resolutions.The figure reveals that a larger time bin results in a smoother SFH curve, albeit at the cost of losing finer details.
However, a smaller time bin is not always the optimal choice.Given the finite number of particles in a galaxy during the simulation, an excessive number of time bins would lead to fewer particles in each bin, thereby causing significant numerical fluctuations.This effect consequently reduces the goodness of fit for the SFHs.As depicted in Figure 18, when the width of the time bin is 10 Myr, the noise in the SFH curve itself is considerably high.Despite this, our fitting method can still identify the main components in the curve, although the goodness of fit is significantly reduced numerically.
Conversely, a larger time bin results in a fitted curve that aligns better with the SFH.This is because using a larger time bin is akin to performing a smoothing process on the curve, which effectively eliminates noise.
In Figure A, we present the statistical goodness of fit of our method for SFHs with different time bins.It can be observed that  2 decreases as the time bin decreases.This relationship is particularly pronounced in low-mass galaxies, primarily due to the numerical noise caused by the insufficient number of particles in these galaxies.
Another point of interest is whether our fitting results would be affected when different time bins are used.According to our inspection results (partially shown in Figure 18), although the number of components and the type of each component (Gaussian or Gamma) may vary, the major peaks can be well restored and maintain stable positions and widths.The main differences lie in the handling of some small components, where different time bins can lead to different results.
We use the fitting results of a time bin of 100 Myr as a benchmark and compare the number of components obtained by fitting the SFHs of other time bins with it.As shown in Figure A, most galaxies' SFH fitting does not change the number of components relative to when the time bin equals 100 Myr ( Δ −  0.1 = 0), where  Δ is the number of fitting components of a single SFH when the time bin is Δ.
For high-mass galaxies, whether the time bin increases or decreases, the change in the number of components is mainly a decrease, i.e., the side where  Δ −  0.1 < 0 is higher than the other side.As mentioned above, this is because the numerical noise brought about by the decrease in time bin can blur some small peaks, while the smoothing effect brought about by the increase in time bin may also erase these peaks.This suggests that 100 Myr is a good choice of time bin for high-mass galaxies in the TNG simulation.
For low-mass galaxies, we can see that the number of components tends to increase when the time bin increases.This suggests that the time bin of 100 Myr may still be too small for these galaxies, and a slightly larger time bin matches their particle number better.
At last, readers should remember that the influence of time bin resolution mainly comes from the mass resolution limit in simulation.In other theoretical galaxy models, i.e., SAMs, or observational data, the influence of mass resolution can be alleviated.

B. STELLAR MASS LOSS
Due to the loss of mass in stellar particles in simulation, the Star Formation Histories (SFHs) constructed based on the mass of particles at z=0 and those based on the initial mass of particles will differ.Figure B illustrates the differences between these two types of SFHs.As the age of each particle remains unchanged, with only the value of its stellar mass varying, the differences between the two SFHs are solely vertical shifts.Consequently, the differences in the fitting results based on the two types of SFHs primarily lie in the amplitude of the components, with no changes in the position and width of each component.

C. FACTORS AFFECTING INTRINSIC SFH PEAK COUNTS
Section 5.1 highlights a discrepancy between the "true fraction" of multiple episodes SFHs in this study and those reported in Iyer & Gawiser (2017); Iyer et al. (2019).To investigate this difference, we emulate the method of identifying intrinsic peaks as described in Iyer & Gawiser (2017) and Iyer et al. (2019), and examine various sets of prominence for peak finding algorithm.As demonstrated in Table 3, a higher prominence in peak finding algorithm results in a smaller fraction of multi-episode SFHs.Based on the value of multi-episode fraction, the SFHs reconstructed by our multi-component method closely aligns with the  The upper panel shows the histograms for high mass galaxies, while the lower panel shows the histograms for low mass galaxies.
Table 3.This table represents the factors that influence the intrinsic peak counts of SFH.The last row provides the reference for the multi-episode fraction of reconstructed SFHs in this work.Note that the multi-component method we proposed in this work does not require prominence parameter to define the peaks.Note that the prominence parameter in this work is the difference of   rather than  , while the latter one is used in Iyer et al. (2019) 2019) is larger than    /4 for the majority of SFHs in the TNG simulation.We apply the prominence from Iyer et al. (2019) to our data and find that the peaks defined by this prominence are more akin to the performance of the prominence of    /2 (refer to Figure 22).

Figure 1 .
Figure 1.Four examples of the fittings to the SFHs.The blue solid line represents the actual SFH data, and the orange dashed line represents the corresponding fitted curve.The dotted lines show the decomposition of components in the fitted curve.The goodness of fit for the overall fitting is displayed in the legend.We illustrates the situations for different level of fitting goodness: good fit containing multiple components (top left), good fit containing one component (top right), median fit (bottom left), and bad fit (bottom right).

Figure 2 .
Figure 2. Histograms of  2 of SFH fittings.The overall distribution is shown as a gray histogram.The histograms of  2 values for SFHs of different types of fittings are shown in different colors.The SFH fittings are classified by the number of components (left plot) or by the shape (Gamma or Gaussian distribution) of dominant component(right plot).

Figure 3 .
Figure3.The two-dimensional distribution of  2 of SFH fittings and stellar mass of corresponding galaxies at  = 0.The gray line displays the median value of  2 for galaxies of different masses, while the error bars represent the range of  2 values from 20% to 80% in each mass bin.

Figure 5 .Figure 6 .
Figure5.The two-dimensional distribution of  2 of SFH fittings and the aga of them.The gray line displays the median value of  2 for galaxies of different SFH lengths, while the error bars represent the range of  2 values from 20% to 80% in each age bin.

PFigure 7 .
Figure 7.The cumulative stellar mass function of galaxies at  = 0. Different colors represent galaxies with different SFH types.The gray line shows the cumulative stellar mass function of all samples.The SFHs types are classified by the number of components (left plot) or by the shapes of dominant components (right plot).In the top sub-panel of each plot, the residuals to the cumulative mass function of all samples for the samples with different SFH types are given.

PFigure 8 .
Figure 8.The cumulative density function of galaxies color ( − ) at  = 0. Different colors represent different SFH types.The gray line shows the color distribution of all samples.The SFHs types are classified by the number of components (left plot) or by the shapes of dominant components (right plot).The top sub-panel of each plot shows the residuals to the color distribution of all samples for samples with different SFH types.

PFigure 9 .
Figure 9.The cumulative density function of average galaxy age with different SFH types.Different colors of histograms represent different SFH types.The gray histogram shows the average galaxy age distribution of all samples.The SFHs types are classified by the number of components (left plot) or by the shapes of dominant components (right plot).The top sub-panel of each plot shows the residuals to the SFH length distribution of all samples for samples with different SFH types.

Figure 11 .Figure 12 .
Figure 11.Examples of C1, C2 and C3 type SFHs with different galaxy mass, color and age.Each column show one kind of SFHs as labeled on the top.Each row represent one galaxy property range as labeled on the left.In each sub figure, the solid blue line represents the origin SFH, and the dashed orange line represents the best fit to it.The colored dashed lines represent decomposed components.

Figure 13 .
Figure 13.Examples of Γ, GΓ and G type SFHs with different galaxy mass, color and age.Each column show one kind of SFHs as labeled on the top.Each row represent one galaxy property range as labeled on the left.In each sub figure, the solid blue line represents the origin SFH, and the dashed orange line represents the best fit to it.The colored dashed lines represent decomposed components.

Figure 14 .
Figure 14.The position and width of peaks of SFHs. ,  is the lookback time of the peaks.   is the full width at half maximum of the peaks.In the middle panel, each point represents data of one peak.The upper panel shows the histogram of the peak positions, and the right panel shows the histogram of the peak widths.For each SFH, we sort the peaks by their height and distinguish the distribution of the peaks in different colors according to their order in the SFH.'Peak 1' is the strongest peak, 'Peak 2' is the second strongest peak, and so on.The three wide dashed lines in the figure correspond to    = 13.6 −  ,  ,  ,  = 0 and    = 1.

Figure 15 .
Figure 15.The histogram of separations between neighboring peaks.

Figure 16 .
Figure16.The dependence of SFH peak position (top row), FWHM (middle row) and neighboring peak separation on stellar mass at  = 0.The colors correspond to logarithmic galaxy count in the corresponding SFH-properties-stellar-mass bins.The median values of the peak position, FWHM, or peak separation in different stellar mass bins are displayed by the gray lines.The error bars denote the range of the same value from 20% to 80% within the bin.The statistics for all galaxies are presented in the left row, those for central and isolated galaxies are in the middle row, while satellite galaxies are on the right row.

Figure 17 .
Figure17.The of SFH peak position (top row), FWHM (middle row) and neighboring peak separation on galaxy color  −  at  = 0.The colors on plots correspond to logarithmic galaxy count in the corresponding SFH-properties-color bins.The median values of the peak position, FWHM, or peak separation in different color bins are displayed by the gray lines.The error bars denote the range of the same value from 20% to 80% within the bin.The statistics for all galaxies are presented in the left row, those for central and isolated galaxies are in the middle row, while satellite galaxies are on the right row.

Figure 18 .
Figure18.The SFHs of three galaxies obtained at different time bin resolutions (blue solid lines), along with the corresponding fitted curve (orange dashed lines) and the decomposition of components (colored dotted lines).Each row contains the results for one time bin.Each column contains the results for one galaxy.

Figure 19 .
Figure19.The histograms of fitting goodness  2 .The histograms of different colors represent the results of the SFHs obtained at different time bin resolutions.The upper panel shows the histograms for high mass galaxies, while the lower panel shows the histograms for low mass galaxies.

Figure 20 .Figure 21 .
Figure20.The histogram of the difference in the number of components in the SFH obtained at different time bin resolutions.  represents the number of components in an SFH.The x-axis shows the difference between the number of components at four other time bin resolutions and the number of components at a resolution of 100 Myr.The y-axis represents the corresponding number of galaxies.The upper panel shows the results for high-mass galaxies, while the lower panel shows the results for low-mass galaxies.

Gawiser Algorithm 1: Fitting algorithm for one SFH Input: 𝑥
and stochastic SFHs (Kelson 2014), Iyer & =   ,  =  (  ) the priors for Gamma fitting are   =   ,   = 1,   = 1 * / / * the ranges of parameters are   ∈ [0, 14],   ∈ [0.1, +∞],   ∈ [0, 100] * / 10 calculate  2 of  , ,  , , , and  , ;( , ,  , ,  , ,  , ) ; ) to get tuned ŷ′′   ; 16 calculate  2 of ŷ′′   ; 17 select the best fitting ŷ = arg max 2 , ..., ŷ′′ 6 ) The ratio of central and isolated galaxies over satellite galaxies for different types of SFHs.The horizon line indicate the ratio for the whole sample.case.So, more components of SFHs indicates more effective star formation time across a galaxy's life.The galaxy age does not show a strong correlation with the shape (Gaussian or Gamma) of the SFH components.It implies that different physical drivers of star formation have relatively weak correlations with the galaxy age, indicating that they do not have a particular preference for the timing of occurrence.The last two rows of Figure sat N cen, total /N sat, total Figure 10.

Table 2 .
Number of SFHs with different numbers of fitting components (left) and actual peaks (right), limited to the sample with fitting goodness of  2 > 0.5.
. We use the savgol_filter function from scipy.signal module.bThisprominence is the same as that in the line ofIyer et al. 2019.c  is the maximum value of each smoothed SFH. a Iyer et al. (2019)17) of the effects of prominence.In each figure, blue solid line represents a origin SFH from the TNG simulation.The orange dashed line represents the smoothed curve using the scipy.signal.savgol_filterfunctionwith a window of 20 points and polyorder 3. The peaks found with a prominence equivalent to that inIyer et al. (2019)are indicated by green circles.Purple squares indicate the peaks found with a prominence of    /2, while red triangles indicate the peaks found with a prominence of    /4.To ensure the clarity of the markers, they have been slightly shifted upwards.In the left figure, three different sets of markers can identify the same peaks in the smoothed SFH.However, in the right figure, only the smallest prominence is able to distinguish the smaller peak.intrinsicSFHsdefined with with a prominence of    −   , >    /4, where    is the max SFR value of one smoothed SFH.Although the smoothing method and prominence in our test can not be aligned withIyer & Gawiser (2017)andIyer et al. (2019)exactly, we can confirm that the prominence accepted inIyer et al. (