Open Access
Issue
A&A
Volume 691, November 2024
Article Number A56
Number of page(s) 15
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202349139
Published online 29 October 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

It has been shown that the global star formation rate (SFR) has steadily decreased over cosmic time from z ∼ 2−3 (Madau & Dickinson 2014), so that galaxies today generally form fewer stars than in the past. However, this decline did not occur at the same rate for all galaxies. On average, the SFR in massive galaxies reached their highest value at a higher redshift and declined faster than in less massive galaxies, many of which are still forming stars (e.g. Speagle et al. 2014; Rodríguez-Puebla et al. 2017; Sánchez et al. 2019). These differences lead to a number of correlations between the stellar mass (M*) and the mean age (e.g. Gallazzi et al. 2005; Blanton & Moustakas 2009; Sánchez et al. 2018) and metallicity (e.g. Trager et al. 2000; Savaglio et al. 2005; Erb et al. 2006; Sánchez-Blázquez et al. 2006a,b; Maiolino et al. 2008; Camps-Fariña et al. 2021, 2022). Moreover, the M* function of galaxies (when observed at different redshifts) shows much smaller changes in the high M* range than at the low end of the scale (e.g. Fontana et al. 2004; Bundy et al. 2005; Pozzetti et al. 2010). These results are known as the downsizing phenomenon (Cowie et al. 1996), which generally suggests that more massive galaxies assemble their mass in shorter timescales. Despite all of the advances made in understanding these processes, we still do not know the details of the transition between star-forming and retired galaxies, what triggers that transformation, when it happens, and what regulates it.

To explore this issue, two different methods could be adopted: (i) study the SFR and M* of different cosmological samples at a different redshift and compare them or (ii) derive the star formation histories (SFHs) of individual galaxies based on their current observed properties – for example, spectral energy distribution (SED) fitting – and infer their properties at different cosmological times. The method we used to calculate the SFH for the considered sample of analysed galaxies is full spectral fitting: a technique that allows stellar populations to be dated within a galaxy (e.g. see Sánchez 2020; Sánchez et al. 2021a). Early studies used the central spectra of galaxies from the Sloan Digital Sky Survey (SDSS) (e.g. York et al. 2000), but in recent years the advent of integral field unit (IFU) surveys such as CALIFA (Sánchez et al. 2012), SAMI (Bryant et al. 2015), and MaNGA (Bundy et al. 2015) has allowed for more robust determinations using the integrated spectrum of the galaxies (e.g. Sánchez-Blázquez et al. 2011, 2014; Cid Fernandes et al. 2013; Ibarra-Medel et al. 2016; Goddard et al. 2017). Sánchez et al. (2019) studied the star formation main sequence (SFMS, Brinchmann et al. 2004; Noeske et al. 2007; Elbaz et al. 2011) over cosmic time and reproduced the results of Madau & Dickinson (2014) using stellar population synthesis techniques.

The way that galaxies transition from being in the main sequence of star formation until they become retired galaxies is still not well understood. There are many different mechanisms that have been proposed and appear to correlate well with the quenching of the galaxies, such as feedback from active galactic nuclei (AGN) or galactic outflows (e.g. Dekel & Birnboim 2006; Sturm et al. 2011; Gaspari et al. 2012; Cicone et al. 2014; Sánchez et al. 2018; Dekel et al. 2019; Trussler et al. 2020), morphological quenching (e.g. Martig et al. 2009; Schawinski et al. 2014; Smethurst et al. 2015), or environment (e.g. Wetzel et al. 2012; Boselli et al. 2014), as the main actors. It is thus very likely that quenching effects act on global scales (e.g. Belfiore et al. 2017, 2018; Trussler et al. 2020) and involve a loss of gas supply (e.g. Boselli et al. 2006; Dekel & Birnboim 2006; Schaye et al. 2010; Dekel et al. 2019; Trussler et al. 2020).

As indicated before, global trends in SFHs can also be studied using cosmic surveys in which galaxies are observed at various redshifts. While these studies lose the advantage of measuring the evolution of the SFR on the same sample of galaxies, they also avoid the uncertainties and degeneracies present in spectral fitting techniques (Worthey 1999; Kuntschner et al. 2001; Sánchez-Blázquez et al. 2011; Walcher et al. 2011; Conroy 2013; Cid Fernandes et al. 2014). Deep observations of cosmological fields on the sky, in particular the Hubble Deep Field and the Great Observatories Origins Deep Survey fields (GOODS; Giavalisco et al. 2004; Dickinson et al. 2003) with the Hubble Space Telescope (HST), have enabled a large number of surveys to measure the SFR and M* of galaxies over a wide range of redshifts (e.g. Salim et al. 2007; Noeske et al. 2007; Daddi et al. 2007; Pérez-González et al. 2008a; Stark et al. 2009; Karim et al. 2011; Grogin et al. 2011; Schreiber et al. 2015; Ilbert et al. 2015; Salmon et al. 2015; Tomczak et al. 2016; Santini et al. 2017), which are currently being dramatically expanded with the introduction of the James Webb Space Telescope (Gardner et al. 2006).

These surveys generally agree well with the general picture of the evolution of the global SFR and SFMS described above (Speagle et al. 2014), but they have significant uncertainties in the absolute values of the SFR and M*. As cosmic surveys employ several methods and they measure galaxies, by design, at vastly different distances, it is also difficult to distinguish features of the SFR- or sSFR–M* relations beyond the slope and zero-point and how these evolve through time. Morphology is also a complicated parameter to include, both due to observational constraints in spatial resolution given the distances involved, and the fact that it is not constant through time (Conselice 2014). This is precisely where the full spectral fitting results shown in this paper can contribute to the understanding of how galaxies grow their mass, as the fact that our sample is composed of the same galaxies at all cosmic times allows us to probe qualitative trends intrinsic to the galaxies regarding how the sSFR–M* relation has evolved in cosmic time.

Specifically, we investigate how star formation is regulated in galaxies by examining which parameters determine the shape and evolution of the relation between specific star formation and stellar mass (sSFR–M*) with cosmic time. Using the MaNGA galaxy sample and applying stellar population synthesis techniques, we could track the evolution of the same galaxies over cosmic time and partition them into current morphology bins that are consistent over cosmic time because they contain the same galaxies that end up in a given morphology bin at all times. We also checked whether galaxies that are currently forming stars show a different evolution than galaxies that are retired. The results are compared with several cosmic surveys adapted to the format of our results.

The text is organised as follows: In Sect. 2 we describe the two types of data used, IFU observations, and a spectrophotometric survey; in Sect. 3 we describe the sample employed and how it was selected; in Sect. 4 we describe the reduction and analysis procedures used to obtain the relevant quantities; and in Sect. 5 we show the resulting relations between the sSFR(t), M*(t), current SFR, and current morphology. In Sect. 6 we discuss the implications of the results and their validity before presenting the conclusions in Sect. 7.

2. Data

Mapping Nearby Galaxies at APO (MaNGA) (Bundy et al. 2015) is a project in the fourth iteration of SDSS (York et al. 2000) to produce maps of the ionised gas and stellar properties for a large number of galaxies in the Local Universe (⟨z⟩∼0.03) as a statistically significant sample to study how galaxies have evolved over their lifetime. Observations were made using a set of fibre bundles fed into the Baryon Oscillation Spectroscopic Survey (BOSS) spectrographs (Smee et al. 2013) to map the spectral emission of the galaxies at each spatial position, aiming to cover most of the area of each galaxy. The coverage primary sample (66%) galaxies includes at least 1.5 effective radii (Re) and for the secondary sample (33%) it includes 2.5 Re. The spectrographs are fed by optical fibres which are bundled to produce integral field units (IFUs) that provide spatially resolved spectral information. The IFUs consist of a primary bundle that can range from 19 (with a diameter of 12″) to 127 fibres (with a diameter of 32″), depending on the object, and 12 sets of seven fibre bundles for flux calibration and 92 fibres for sky subtraction (Drory et al. 2015; Yan et al. 2016). The fibre bundles are inserted into precision perforated plates prepared for each night’s observing targets and then mounted on the Apache Point Observatory 2.5-m Sloan Telescope at Apache Point Observatory (Gunn et al. 2006).

The raw data from the telescope are processed with the MaNGA data reduction pipeline (DRP, Law et al. 2016) into the project’s core products (the IFU data cubes that we use here). Each data cube has a spectral range between 3600 Å and 10 300 Å at a resolution of R ∼ 2000, with an effective spatial resolution after reduction of about 2.5″/FWHM.

3. Sample

The MaNGA sample consists of ∼10 000 galaxies selected in the Local Universe with the primary criteria being that (i) they are a representative sample of the local galaxy population and (ii) they are well covered by the field of view (FOV) of the instrument. The latter is necessary to obtain 2D maps on several properties such as M* density, mean stellar age, SFR, metallicity, etc.

We apply some selection criteria to the parent sample to obtain a sub-sample suitable for our study, removing galaxies with an inclination above 70 Å and those in which an AGN has been detected (from Lacerda et al. 2020). Measurements in highly inclined galaxies are heavily affected by dust lanes and it is generally harder to obtain quality SFHs, while wide emission lines from AGN can also affect the analysis whilst being a relatively small number of galaxies. The refined sample consists of 9087 galaxies and is the same as that used in Camps-Fariña et al. (2022).

We use the Vázquez-Mata et al. (2022) morphological classification (v2.0.1), which is available as a Value Added Catalog1. The classification procedure uses r-band SDSS and Dark Energy Spectroscopic Instrument (DESI, Dey et al. 2019) images which are processed with a sharpening filter and have surface brightness models subtracted to detect sub-structures such as bars.

The main morphology bins used in this article can be reproduced by selecting the following morphology labels from the Value-Added Catalog:

  • E-S0: E, E(dSph), S0, S0(dwarf), S0a, SAB0, SAB0a, SB0, SB0a.

  • Sa-Sb: SABa, SABab, SBa, SBab, Sa, Sab, SABb, SABbc, SBb, SBbc, Sb, Sbc.

  • Sc-Irr: SABc, SABcd, SBc, SBcd, Sc, Scd, Irr, IrrAB, IrrB, SABd, SABdm, SBd, SBdm, Sd, Sdm, dIrr, dSph, dwarf.

In other words, we segregate galaxies into: (i) early-type ones (bulge dominated), (ii) early spirals (galaxies with both bulges and discs), and (iii) late spirals (galaxies without a significant bulge, mostly disc dominated).

4. Analysis

The composite optical-near-infrared stellar spectrum of a given region in a galaxy varies as a function of the SFH, the metallicity of the stars, and the present-day dust extinction. The process of fitting a library of stellar population spectra to observations to determine the composition of the underlying stellar population allows us to study the past of a galaxy using features left behind in the current stellar populations.

4.1. pyPipe3D

pyPipe3D (Lacerda et al. 2022) is an analysis pipeline used to obtain the ionised gas and stellar population parameters of IFU observations such as emission line fluxes, velocities and widths and stellar population ages and metallicities, among others. It also provides us with the relevant 2D maps of these parameters. The stellar spectra are fitted using templates from the library of simple stellar populations (SSPs) specified by the user, returning the light fractions of the populations for each given spectrum. It is an updated version of the Pipe3D analysis pipeline (Sánchez 2006; Sánchez et al. 2016a,b) rewritten in Python, with added modularity and ease of use, and a significant increase in computational speed. It separates the emission and absorption components and fits them separately to measure and generate spatially resolved maps of the relevant physical quantities of the galaxy. Internal dust attenuation is determined using a two-step procedure as extensively described in Lacerda et al. (2022) and Sánchez et al. (2016a). The procedure performs an estimation of the dust attenuation by doing a brute force exploration in a fixed range of values, fitting the stellar population using a limited subset of the stellar template that limits the possible degeneracies between properties. A Cardelli et al. (1989) extinction law is applied.

4.2. Obtaining the star formation relations

The M* of a galaxy can be determined from the light fractions provided by pyPipe3D as follows:

M = i N age j N Z w i , j Υ i , j L V , tot $$ \begin{aligned} M_* = \sum _{i}^{N_{\rm age}}\sum _{j}^{N_Z} \, { w}_{i,j} \, \Upsilon _{i,j} \, L_{V,\mathrm{tot}} \end{aligned} $$(1)

where wi, j and Υi, j are the light fraction and the mass-to-light ratio of the SSP template with a given age and metallicity, respectively. We can obtain the mass assembly history (MAH), which is the mass of all the stars present at a given age of the Universe (tU = 13.7 − t) as:

M ( t ) = i > t N age j N Z ( 1 + F loss ( i t ) ) w i , j Υ i , j L V , tot $$ \begin{aligned} M_*(t) = \sum _{i > t}^{N_{\rm age}}\sum _{j}^{N_Z} \, \left(1+F_{\rm loss}(i - t)\right) \, { w}_{i,j} \, \Upsilon _{i,j} \, L_{V,\mathrm{tot}} \end{aligned} $$(2)

where Floss(i − t) is the fraction of M* returned to the ISM by a SSP of a given age and metallicity in the time interval between t and now (t = 0). The sSFR is defined as the SFR divided by the M* so we can trivially obtain the sSFR history with the SFH and the M*(t).

Due to the nature of the sample selection performed for MaNGA, which prioritised a good coverage of the FOV, there is a significant redshift range for these objects (z ∼ 0.0002−0.15) which correlates with M*. This results in more massive galaxies being observed at slightly earlier epochs compared to less massive ones.

We can consider the SFR for each look-back time (LBT) as equivalent to a measurement of a single object at a given redshift in a cosmic survey, similar to the methodology used in Sánchez et al. (2019). We thus select several discrete LBT values and measure the SFR and M* of each galaxy at each time point. We do this by interpolating the SFH and M*(t) curves, since the SFH of each galaxy is shifted by the value of its redshift and therefore, the ages of the templates correspond to slightly different LBT values. The lowest LBT value for which we measure the sSFR–M* relation is 1 Gyr. This is due to the aforementioned redshift range in the sample, which makes it so that for very low LBT values many galaxies are not observed as the light travel time is longer than 1 Gyr for the distance they are located at. In other words, for some galaxies the light we currently measure was emitted longer than 1 Gyr ago due to their distance. At 1 Gyr 87% of the sample is included, which we considered a good compromise for picking the lowest LBT value.

5. Results

In Fig. 1 we show the sSFR–M* relation for the galaxies in our sample. The SFR is measured by averaging it over the last 30 Myr which is a typical value for the lifetime of an HII region and this is also the timescale at which the SFR measured using Hα and full spectral fitting coincide (Asari et al. 2007; González Delgado et al. 2016). Similarly to the SFMS, there is a bi-modal distribution with a fairly flat upper distribution of star-forming galaxies and a lower distribution of retired galaxies, which are connected by the Green Valley.

thumbnail Fig. 1.

Currently observed sSFR of MaNGA galaxies compared to their M* (bottom left). The black contours enclose 95%, 65% and 35% of the sample respectively. The green line is the median sSFR within 0.25 dex wide M* bins. In the top and right panels, the distributions in M* and sSFR, respectively, are shown as histograms, with the three morphological bins shown in lines of the corresponding colour.

In Fig. 2 we show the sSFR–M* distribution measured at different LBT by evaluating the SFH and MAH at these LBT values. At 1 Gyr, the relation is very similar to that shown in Fig. 1, with a general negative trend between the two parameters which steepens at about log M*/M = 10.5. The retired galaxies no longer have a clearly defined separate distribution, likely the effect of measuring at a specific LBT rather than the last 30 Myr, meaning that for a portion of these galaxies the SFH is measured at the time when they are still becoming quenched, widening the distribution.

thumbnail Fig. 2.

Evolution of the sSFR over cosmic time for the galaxies in the MaNGA sample. In each panel the M* and sSFR of the galaxies measured at different LBT are shown. The colour of the data points corresponds to the currently observed morphology of the galaxies. The black contours enclose 95%, 65% and 35% of the sample, respectively. The cyan and violet dots and lines show the median sSFR within 0.5 dex wide M* bins, selecting galaxies above and below 6 Å in EWHα, respectively. The error bars for the dots show the 25th and 75th percentile of the distribution within each M* bin.

We show averaged sSFR values of the distribution within M* and current star formation bins as large dots on Fig. 2. The M* bins are selected at each LBT but the star-forming (SFG) and retired (RG) galaxy bins are selected based on whether the current equivalent width in Hα (EWHα) of the galaxies is higher or lower, respectively, than 6 Å (Cano-Díaz et al. 2016). Each bin that is represented has at least ten objects in it. SFG and RG are still segregated at 1 Gyr and they also have different sSFR–M* slopes such that at higher M* RG have much lower sSFR compared to SFG. The bend at log M*/M = 10.5 is mainly observed in SFG rather than RG. The sSFR–M* relation for both star formation bins appears to converge at low masses such that they have similar sSFR at 1 Gyr in LBT, despite having different EWHα values. The most likely explanation for this is that low M* galaxies have their star formation in short, intermittent bursts (e.g. Lee et al. 2009; Domínguez et al. 2015; Atek et al. 2022). As such, the emission line based separation using EWHα would not correlate with a significant separation at 1 Gyr. Galaxies which are currently star-forming are not as statistically likely to be so once we look at the recent past if they tend to have short bursts of star formation.

The relation between the sSFR and M* becomes flatter at higher redshifts, indicating that the SFR used to be roughly proportional to the M* at earlier times, which creates a flat relation between M* and the sSFR. A flat correlation means that galaxies double their M* on similar timescales and the currently observed negative correlation means that more massive galaxies take longer to double their M*. As such, we propose that the phenomenon of downsizing is not best described as more massive galaxies having a quicker evolution than less massive galaxies. Instead, the self-regulation of the SFR means that galaxies grow at a similar pace (similar M* doubling times) intrinsically but once galaxy-wide quenching mechanisms become efficient their growth is stopped, a phenomenon that happens earlier in more massive galaxies. Therefore, it is not that the evolution of the latter is intrinsically quicker (defined by M*-doubling times) but that their steady-state of star formation was shut down at an earlier time than less massive galaxies. As such, downsizing is defined by the age at which quenching mechanisms become efficient.

The difference between star-forming and retired galaxies (in the Local Universe) quickly shrinks with increasing LBT, showing once more that these galaxies used to form stars at a similar pace before some of them shut off further star formation. The timing for currently retired galaxies to become quenched depends on M*, with less massive galaxies reaching similar sSFR values at 3 Gyr in LBT (z ∼ 0.2) but more massive ones only do so at 5 Gyr (z ∼ 0.4) and earlier. Interestingly, there is a small but consistent gap between RG and SFG at the highest LBT of about 0.11−0.23 dex, showing that currently retired galaxies had slightly higher sSFR at the beginning. This implies that galaxies that are retired in the present day are statistically more likely to have been highly efficient at forming stars in the early universe. As such, there is a degree of intrinsic downsizing such that some galaxies did have a quicker evolution, but effect does not appear to be strongly tied to M*, unlike the global quenching described above.

At 3−7 Gyr in LBT (z ∼ 0.2−0.8) the lowest M* galaxies show a deviation from the general trend such that the relation between sSFR and M* significantly steepens, especially at 5 Gyr (z ∼ 0.4). The number of galaxies involved is small relative to the full sample except at 5 Gyr, where the change in slope occurs at around log M*/M ∼ 9 − 9.5. Were this a real feature, it would imply that galaxies with log M*/M ≲ 8.5 − 9 have the opposite trend to that observed for galaxies above this M*, that is, that for earlier LBT the relation between sSFR–M* becomes steeper instead of shallower.

5.1. Morphology

In this section we explore if we can distinguish galaxies with different morphologies according to their SFH and how far into the past this segregation is maintained. The morphology of galaxies in the local universe is closely related to their SFH (e.g. Camps-Fariña et al. 2021, 2022). Thus, to some extent, we should expect to find additional correlations if we divide the sample into morphology bins before computing the evolution of the sSFR–M* relationship. At the very least, we should find a more pronounced decline in sSFR over cosmic time for early-type galaxies, whose retired fraction is much higher than for late-type galaxies, which tend to have ongoing star formation.

A crucial difference for the comparison between morphology bins is that, unlike M* bins, we are not able to distinguish the morphology of galaxies at different cosmic times. For this reason, these morphology bins must be interpreted as representing galaxies that end up in a particular morphological bin. In addition, the current morphology classification is unlikely to be representative of the galaxy population at all cosmic times (e.g. Conselice 2014).

In Fig. 3, we show the same data points as in Fig. 2 but we overlay the contours of galaxies divided in three morphological bins: early type (E and S0), intermediate type (Sa and Sb), and late type (Sc, Sd, and Irr). At 1 Gyr in LBT the morphological types have clear differences in both the extent of the M* range that they cover, their average sSFR and the slope of the correlation between M* and sSFR. Early-type galaxies are widely distributed in M* and have lower sSFR in general as well as having a steep negative correlation such that more massive galaxies have much lower sSFR than less massive ones. On the other hand, the M* distribution of late-type galaxies is narrower and centred at lower values, and the correlation becomes flatter. Late-type galaxies show a fairly flat distribution while intermediate ones have a clear bend in the relation such that the slope increases above log M*/M ∼ 10−10.5. The difference in slope between morphological bins results in a convergence of sorts at lower M* such that at log M*/M ∼ 9−10 the differences in sSFR due to morphology are relatively small.

thumbnail Fig. 3.

Evolution of the sSFR over cosmic time for the galaxies in the MaNGA sample. In each panel the M* and sSFR of the galaxies measured at different LBT are shown. There are three sets of contours whose colours correspond to the current morphology of the galaxies and whose levels enclose 95% and 50% of each sub-sample.

The differences between morphological types at 1 Gyr are similar to those between SFG and RG seen in Fig. 2 as would be expected due to the correlation between star formation status and morphology. The relation between sSFR and M* becomes shallower at earlier LBTs for all morphology bins, most clearly for E+S0 and Sa+Sb but even Sc-Irr (whose distribution is fairly flat at all times) can be seen to have a shallower relation at 11 Gyr compared to 1 Gyr. The main change in the latter morphological bin is that its distribution greatly widens between 3−7 Gyr, perhaps indicating a period of increased burstiness. Early-type galaxies maintain a shallow relation between sSFR and M* for LBT ∼ 11 to 7 Gyr, but at 7 Gyr the centre of the distribution has already started to drop and by 5 Gyr it is clearly located at lower sSFR compared to the other morphological bins. This indicates that their star formation slows down at around 7 Gyr compared to galaxies in the other morphological bins before the current negative correlation between sSFR and M* was established.

At high LBT we see that the galaxies in each of the current morphology bins have similar distributions, but they are distinctly segregated in M*, even more so than at 1 Gyr and especially at the lower M* end of the distributions. Such a clear difference between the galaxies at such an early time despite the bins corresponding to only current properties of the galaxies shows that current morphology is determined by the entire history of a galaxy’s evolution starting at the very beginning of its lifetime. It also suggests that currently retired galaxies are not, on average, simply star-forming galaxies with a recent and rapid quenching but that they tend to have different SFHs even at early times.

There is also a small (0.1 dex) separation between the average sSFR between early- and late-type galaxies, with early-type galaxies having slightly higher sSFR. We expect this segregation to be related to the one found at this same LBT between RG and SFG in Fig. 2; however, the segregation is significantly smaller between the morphology bins despite us measuring the distance between the two extremes of three bins instead of the difference between two bins. As such, we can conclude that the segregation in sSFR at early times is more likely to be related to whether galaxies continue to have star formation in the present than their morphology.

Of course, the opposite interpretation is also valid: More massive galaxies which have high star formation efficiency (relative to M*) at early times are less likely to end up in later type morphologies compared to galaxies of similar M* but with lower early sSFR. Either way, a clear link between M* at early times and present-day morphology is established. This result points towards a possible link between large scale mechanisms that make a galaxy grow quicker than average in the early universe and the eventual morphology of the galaxies. One possible explanation would be that galaxies which are formed in deeper dark matter potential wells accrete more gas at early times producing an early growth, but it also makes them more likely to attract other galaxies to interact with, eventually favouring an early-type morphology as a result of the interactions. The existence of long-term differentiation of SFHs between galaxies has been proposed by Caplar & Tacchella (2019), Matthee & Schaye (2019) from results in simulations.

The fact that the distributions are segregated more at low LBT than at high ones suggests a ‘convergence’ of sorts in the MAHs of galaxies of different morphology. In other words, galaxies with different morphology which have different M* at very early times end up, on average, with M* values which are more similar in later times. This implies that late-type galaxies have formed a higher percentage of their M* in the 1−11 Gyr interval in LBT compared to early-type ones, which could be an effect of gas accretion. Cosmic gas is expected to (i) be preferentially accreted at the outskirts of galaxies (Sánchez Almeida et al. 2014; Sánchez Almeida 2017) and (ii) to form discs, making so that galaxies with lower initial M* in our sample are statistically more likely to develop discs and therefore they are more likely to be late type currently. Camps-Fariña et al. (2023) measures the gas accretion history of galaxies and finds a correlation such that late-type galaxies have higher gas accretion rates in their recent history, in line with this interpretation.

5.2. Statistical significance of the populations

The previous figures show how the different populations of galaxies in terms of M*, morphology and present-day star formation have different evolutionary paths well into their past. The differences appear to be fairly consistent between ages for all the galaxy populations, but it is still important to use proper statistical metrics to determine how significant their differences are and how far into the past do these populations remain distinct.

For this purpose, we have performed a Kolmogorov-Smirnov (KS) two-sample test for each relevant pair of distributions to see how likely it is that the two populations have the same origin. A p-value lower than 0.05 tells us that it is extremely unlikely that the two samples are the same. The relevant pairs are: one for the difference between current SFG and RG, and three pairs for all the combinations between early-, intermediate-, and late-type galaxies. The distribution in log M*/M within each population is different and therefore comparing the entire population could lead us to overestimate how different their distributions in sSFR are. As an example, consider a single sample of galaxies which have a correlation between log M*/M and sSFR and no secondary correlations. If we divide the sample into a high-M* one and a low-M* one their average sSFR will be different due to the correlation with M* and we would obtain a false result of two different samples.

To avoid this, we have performed the KS test within the 0.5 dex wide M* bins employed in the averaged SFG and RG populations in Fig. 2. This removes any dependence on log M*/M as now we measure the p-value for the distributions of sSFR for each population within a given M* bin. The resulting p-value vs log M*/M inform us of which M* ranges show significant differences in the distribution of sSFR values of the populations considered. The p-values for each M* and LBT bin in Fig. 4 prove that, for the M*-range where we have good statistical sampling, log M*/M ∼ 9 − 11 the samples considered are clearly distinct, with values much lower than the 0.05 threshold. In fact, even at 11 Gyr in LBT the median p-value between log M*/M = 9 − 11 for all pairs of populations is 1.6 ⋅ 10−15 indicating that, even at the highest age we consider, these populations remain very significantly distinct.

thumbnail Fig. 4.

Two-sample Kolmogorov-Smirnov test for each pair of galaxy populations considered in the article. In each panel we show the p-value for each pair of populations performed individually at each of the 0.5 dex wide M* bins. The dotted line shows the limit of p-value = 0.05, below which the samples are considered to be statistically significantly different. The colours identify the pairs of populations used.

There is also a possible interesting result showing that at log M*/M = 10.5 there are two LBT values with an unusually high similarity in one of the morphology pairs. This is the M* where the Green Valley in the SFMS appears and is considered a watershed mass for galaxies to start to quench (e.g. Dekel & Birnboim 2006; Dekel et al. 2019). The fact that we observe first an isolated very high similarity between the early and intermediate and then another one between intermediate and late (as seen in the present) could indicate the moment in which the currently early and the currently intermediate samples started to quench. Comparing this to Fig. 3 we do see that at these times and M* the relevant pair of distributions appear to have similar values as the sSFR starts to drop for earlier-type morphology. This is not seen in the SFG-RG distributions and the morphology bins become distinct again in the next LBT, pointing to a transitional period with a ≲2 Gyr timescale.

5.3. Comparison with cosmological surveys

In this section we compare our results with those of other studies which measure the evolution of either the SFR–M* or the sSFR–M* relations by measuring these parameters at different redshift values based of different galaxy surveys. Given the drastically different methodology employed by these studies and our own, it is very useful to compare them to see if our results agree with theirs.

In Fig. 5 we show the evolution of the sSFR–M* relation using eight different cosmic surveys:

thumbnail Fig. 5.

sSFR–M* relation for eight cosmic surveys. The data have been adapted into the M*–sSFR plane in some cases. The colour of each line corresponds to the average redshift of its bin and the redshift colour scale is the same for all panels. The dotted lines correspond to M* intervals where sample completeness is not guaranteed as reported in each work. The lower right corner of each field indicates the median scatter of the data.

  1. Pérez-González et al. (2008a) (PG08) selects galaxies from the Hubble Deep Field-North (HDF-N; Williams et al. 1996), Chandra Deep Field-South (CDF-S; Szokoly et al. 2004) and the Lockman Hole Field (LHF; Lockman et al. 1986) using Spitzer-IRAC 3.6 μm and 4.5 μm bands. These observations are complemented by I-band ground-based observations to add galaxies which are faint in NIR. The redshift, M* and SFR were calculated via spectral energy distribution (SED) fitting of photometric measurements in the ultraviolet (UV), optical and near- and mid-infrared (NIR/MIR). In their Fig. 8, they show the redshift evolution of the sSFR for five M* bins, which we transformed into the sSFR–M* relation measured at six redshift values.

  2. Karim et al. (2011) (K11) selects galaxies in the Cosmic Evolution Survey (COSMOS; Scoville et al. 2007) fields from their Spitzer-IRAC 3.6 μm emission. The redshift is taken from the catalogue in Ilbert et al. (2009), determined using far-UV (FUV) to MIR photometry, the M* was estimated with a parametric SED fit to the photometry data, and the SFR was measured from the 1.4 GHz luminosity of the sources, measured with the Very Large Array (VLA). Their Fig. 5 shows the evolution of the sSFR–M* relation, so no conversion was required, and we selected only half of the redshift bins to avoid clutter in the figure due to having too many lines.

  3. Schreiber et al. (2015) (Sc15) selects galaxies from GOODS, Ultra Deep Survey (UDS; Galametz et al. 2013) and COSMOS fields, using H-band images except for GOODS-N which uses Ks band images. The redshift is calculated using photometry from NUV-to-8 μm bands, the M* is measured using parametric SED fitting and the SFR is calculated as the sum of the SFR obtained from UV and IR luminosity values. In their Fig. 10, they show the SFR–M* relation for different redshift values, which we converted to sSFR values by dividing them by the M*.

  4. Salmon et al. (2015) (Sa15) selects galaxies from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS; Grogin et al. 2011) GOODS-S field in H-band images. The redshift is calculated with a SED fit to FUV-MIR photometric measurements from the HST, Spitzer and ground-based observations, while the M*, dust attenuation, stellar ages and SFH shape are calculated with a Bayesian approach to SED fitting to the same data. The attenuation is used to correct the UV emission to calculate the SFR with the stellar ages and SFH used to further refine the values. We adapt the M* and SFR values shown in their Fig. 11.

  5. Tomczak et al. (2016) (T16) uses the sample from the FourStar Galaxy Evolution Survey (ZFOURGE) IR survey (Straatman et al. 2016) which selects galaxies on the CDF-S, COSMOS, and UDS fields in the Ks band. The redshift and M* of the galaxies are obtained with separate SED fits to the five 0.3−8 μm photometry bands using different codes. The SFR is measured using archived Spitzer and Herschel measurements and calculated combining rest-frame UV and IR emission. In their Fig. 3 they show the SFR–M* relation.

  6. Santini et al. (2017) (S17) uses data from the HST Frontier Fields (Lotz et al. 2017) of four galaxy clusters, selecting the galaxies in the H-band image as well as adding faint detections performed on a weighted average of the Y, J, JH, and H bands. The redshift is measured spectroscopically for a small portion of the sample and via SED fitting of optical-NIR photometry. The M* was derived via parametric SED fitting and the SFR by correcting the UV emission for dust attenuation. In their Fig. 2 they show the SFR–M* relation.

  7. Barro et al. (2019) (B19) selects galaxies from the GOODS-N field in HST H-band images. The redshift is calculated via SED fitting using Survey of High-z Absorption Red and Dead Sources (SHARDS; Pérez-González et al. 2013) photometry and HST grism spectroscopy, which is roughly equivalent to low-resolution spectroscopy of R ∼ 50 − 210. The M* is calculated via separate SED fits using UV-FIR photometry and the SFR is primarily based on rest-frame UV and IR emission which is corrected for dust-attenuation via information from several estimators depending on redshift. They give the values of SFR–M* relation calculated with running medians in their Table 15.

  8. Mérida et al. (2023) (M23) combines the CANDELS data-set obtained in the five CANDELS fields, namely, the source catalogues from GOODS-N (B19), GOODS-S (Guo et al. 2013), COSMOS (Nayyeri et al. 2017), UDS, Extended Groth Strip (EGS; Stefanon et al. 2017), with a sample of faint sources that they call SHARDS/CANDELS faint, performed in GOODS-N. They obtain the latter by combining observations from HST and SHARDS, looking for fainter sources via image stacking in the optical. Redshift, M* and SFR values are computed similarly to B19. The sSFR–M* relation we show in Fig. 5 uses the M*, SFR values from their Table 6, which considers only the fraction of galaxies within their mass-complete limits.

The main results in this article that we can compare with these studies are (i) the progressive steepening of the sSFR–M* relation over time and (ii) the faster steepening above log M*/M ∼ 10.5. The results related to morphology and present-day star formation cannot be easily reproduced using cosmic surveys. As for (i), PG08, K11, Sc15 and T16 show this quite clearly, although only PG08, K11 and Sc15 clearly show (ii), the high M* inflexion point at log M*/M ∼ 10.5. PG08 shows the clearest change in slope at 10.5 in log M*/M, while K11 shows more of a peak in sSFR at this M* for the earliest redshift, moving to lower M* over time. Sc15 shows no clear transition point except for the two lowest redshift bins, where the relation becomes noticeably steeper at log M*/M ∼ 10.5.

Regarding morphology as a key parameter for M* growth, a companion publication to PG08 (Pérez-González et al. 2008b) uses infrared data to measure the mass assembly rate of galaxies up to z ∼ 2 (∼10 Gyr), and separates the sample into spheroidal and disc galaxies as proxies for early and late-type galaxies. They found a general decrease in sSFR over time with a difference between the two morphological groups, so that at the highest redshift values their sSFR is practically the same (within the uncertainty range), but they differ in more recent times, with the spheroidal group decreasing faster. There is a correlation between M* and galaxy morphology, so this trend cannot be attributed to morphology alone. However, the general conclusion is that more massive earlier-type galaxies show a faster and stronger decline in sSFR compared to less massive later-type galaxies. This fits well with our results, although our sample and method allow us to break the correlation between M* and morphology.

6. Discussion

The results shown here provide a general picture on how galaxies have grown their M* over cosmic time and how this growth is regulated. Using the sSFR rather than the SFR makes it more easy to directly compare the SFH between galaxies as we are controlling for a general indicator of how ‘big’ a galaxy is: in statistical terms, we expect galaxies with higher M* to have higher gas and halo masses, and also to be larger in physical size and have deeper potential wells.

It is argued that the SFMS is not a physical relation, but rather the result of other, more fundamental relations (Lin et al. 2019; Sánchez et al. 2021b) such as the Schmidt-Kennicutt law (Schmidt 1959; Kennicutt 1989) between the surface densities of molecular gas and SFR, and the molecular gas main sequence (Wong et al. 2013; Barrera-Ballesteros et al. 2020) between the surface densities of M* and molecular gas. Thus, to interpret the results, it is important to keep in mind that the correlations we detect do not necessarily imply a direct causal link between the parameters. Instead, they could arise from a combination of other, unobserved, correlations between more fundamental parameters in determining how galaxies grow. For example, many of the correlations with M* are likely correlations with halo mass, but we still detect them due to the intrinsic correlation between halo and M*. Alternatively, it is known that environment is a strong driver of galaxy evolution (e.g. Wetzel et al. 2012; Boselli et al. 2014) and strongly correlates with the morphology of the galaxies (e.g. Martig et al. 2009; Schawinski et al. 2014; Smethurst et al. 2015). It is quite feasible that environmental effects are at least part of the reason that, on average, earlier-type galaxies have a different evolution, as galaxies in more dense environments are both more likely to be early type and massive, and also have quenched their star formation earlier in their lifetime.

The accuracy of the SFH determination correlates with the age of the populations. The older a population is, the smaller the changes in its integrated spectrum due to the progressive loss of stars in a given interval. The most massive stars produce most of the luminosity, although they account for a much smaller fraction of the total M* and they also die much faster than less massive stars. While the representative spectrum of a given population changes greatly between 0 and 500 Myr in age, the spectra at 10 Gyr and 10.5 Gyr in age, for example, are virtually indistinguishable. The inescapable conclusion is that stellar population templates require wider sampling for older ages (see Conroy 2013; Walcher et al. 2011; Ibarra-Medel et al. 2016). As a result, SFHs derived from stellar population synthesis become less reliable for the oldest ages and are generally unable to resolve the initial stages of galaxy formation. At the old end of the age range, the sampling of the SSP we employ has a gap of 1.7 Gyr between ages for the three oldest ones, which becomes a much higher ‘error’ in terms of redshift. The difference between 11 Gyr and 12.7 Gyr in LBT is 2.5 to 5 in redshift, highlighting one of the most striking differences between cosmic surveys and full spectral fitting regarding the measurement of when phenomena happen. Cosmic surveys lose sample completeness at higher redshift but in exchange gain ever increasing time resolution, while full spectral fitting always maintains complete sample integrity, but age measurements at high z are highly unreliable. This is compounded with the fact that older populations have much smaller light fractions, further adding to the uncertainty in age determinations. Despite this, full spectral fitting works quite well at finding differences in age between spectra assuming that the signal to noise, spectral resolution and spectral range are adequate. For this reason, while the absolute values of z at which galaxies quench are not robust, the observed differences between bins should be reliable.

For this reason, we consider the fact that several of the surveys show the same basic behaviour of the sSFR–M* relation becoming flatter at higher LBT to be an agreement with our results, despite the fact that the quenching phase of the high-mass galaxies in the cosmic surveys occurs at significantly higher z than in our results. It is also worth pointing out that the surveys show significant differences between them in this regard, In PG08 the quenching occurs between z ∼ 2.5 and 1.75 whereas in K11 it occurs between z ∼ 1.8 and 1.1. Sc15 and T16 shows a steady steepening with no clear drop-off time. Our results suggest z ∼ 0.8, lower than either of the first two but not far from the lower value of 1.1 in K11. The inflection point at log M*/M = 10.5 that we see coincides very well with that of PG08 and K11, but it is not clearly observed in any other survey. M23 is exempt as it does not reach beyond 10.5. We also observe a general offset in the sSFR values between our results and that of cosmic surveys, likely related to some systematic offset in our measurements and at least partly explained by the fact that we do not have full coverage of the galaxies within the FoV of MaNGA (see Appendix B).

Despite this, even among the cosmic surveys there are significant differences in the sSFR. At z ∼ 0.5 and 1.2, for log M*/M = 9.5, it varies between surveys up to 0.5 dex, while at higher z ∼ 4 it varies 1 dex or more between surveys. These differences should not be construed as these surveys being wrong, high-z measurements are difficult to constrain well and depend on the tracers of SFR employed. This is where our results can be useful despite the unreliability of the age determinations at high z and loss of resolution past z ∼ 1 − 2, as we can measure slight differences between sub-samples of galaxies consistently.

Compared to cosmological simulations, Behroozi et al. (2019) in their Fig. 18 shows the evolution of the sSFR over cosmic time for different halo mass bins, where initially all mass bins have a similar sSFR value and evolve in parallel, but the more massive galaxies gradually drop in sSFR. This is consistent with our results regarding the fact that galaxies used to grow apace but more massive ones dropped off the SFMS. Davé et al. (2016) similarly shows in their Fig. 6 the evolution of the sSFR–M* relation in MUFASA simulations which goes from a negative relation with a bend at log M*/M ∼ 10.5 to a fairly flat (slightly positive correlation even) at high redshift. In the same figure they find the same behaviour comparing to EAGLE (Crain et al. 2015; Schaye et al. 2015) and ILLUSTRIS (Genel et al. 2014; Vogelsberger et al. 2014) simulations. Like the comparison to the cosmic surveys, the redshift values at which the quenching occurs does not match our measurements well, which coincides much better with the results from cosmic surveys. This once again points towards the cosmic clock in our results being unreliable compared to cosmic surveys.

6.1. The correlation at low M*

In the low end of the log M* range we find a steeper correlation between M* and sSFR below log M*/M ∼ 9 for the intermediate LBT values 3−7 Gyr. This is an interesting feature if real, as it suggests that these very low M* galaxies have a difference in the pace at which they grow their M*. It could be the result of either an increase in star formation efficiency for lower M* galaxies or an increase in the ratio between molecular gas mass and stellar gas mass for lower M* galaxies, that is, a break in the molecular gas main sequence (MGMS). For the latter case, it could once again arise from either a generally larger amount of gas or a higher efficiency in converting HI to H2. Lower M* galaxies are expected to be heavily affected by stellar feedback expelling and heating the gas present in them so the explanation that these galaxies have higher molecular gas mass ratios fits better with the current paradigm on galaxy evolution models.

In Fig. 5 we showed the results from surveys and in several of them (Sa15, Sa17 and B19) this steepening of the relation at early times can be observed, especially for B19. It generally occurs at higher redshift values compared to our results (z ∼ 0.2 − 0.8), B19 shows it most clearly at z ∼ 2 − 4, though at z = 0.25 there is a slight steepening between log M*/M ∼ 8 − 8.75 with an increase of about 0.27 dex. This feature, however, is prominent at the M* ranges where the samples are not complete and it is typically explained away (and rightly so) as a selection effect arising from the fact that more star-forming galaxies tend to be brighter and as such in the fainter end of the relation they are preferentially detected. This selection effect grows for lower M*, fainter, and higher redshift galaxies and the negative correlation which becomes more prominent at early times arises as such.

For M23 the sSFR–M* values we show correspond to the M* range where completeness was guaranteed, but below the completeness limit many low to very low M* galaxies are also reported. In their Figure 13, the distribution in M*–SFR for the lowest redshift appears bimodal, with one population following the SFMS and another parallel but above it. The double peak in the distribution quickly fades at higher redshift values, leaving only the distribution above the SFMS. Including galaxies from their parent sample that are less massive than those included between their M* limits would give similar results of a flattening of the SFMS or steepening of the sSFR–M* relation. It bears mention, however, that this bimodality is only reported above z = 1 and therefore it again does not exactly coincide with our results regarding when it is observed. Given that z = 1 is the lowest z bin it is possible that it persists at later times, but it is not guaranteed.

Our results, by contrast, are free of the selection effects present in cosmic surveys as all galaxies are selected in the Local Universe. The MaNGA sample is reported to be mass-complete above 5 × 108 M (Wake et al. 2017). Even if the mass-completeness limit was underestimated and there were SFR selection effects in the MaNGA sample, our result at low M* cannot be reproduced. This is because the selection effect only applies to the most recent sSFR value for each galaxy rather than at specific LBT values. Despite this, the spectral fitting procedure we employ is affected by degeneracies and method inaccuracies as described in the previous sub-section. It is therefore possible, especially as these galaxies have lower observed flux due to being in the lower M* end, for this feature to be method-related, meaning we have two scenarios:

  • We observe a coincidental agreement between the effects of sample selection in cosmic surveys and method degeneracies from our spectral fitting.

  • There is an unreported correlation for low-M* galaxies which has so far been attributed entirely to selection effects intrinsic to cosmic surveys.

Further explorations are required to confirm whether this feature is real and understand the underlying reason for what would effectively be a flattening of the SFMS at very low M*. Observations with JWST and in the near future with instruments such as MOSAIC in the Extremely Large Telescope will provide constraints for the evolution of low M* galaxies in cosmic time.

6.2. Compatibility with theoretical models

6.2.1. The golden mass

The golden mass is an explanation proposed by Dekel & Birnboim (2006), Dekel et al. (2019) for the existence of a bi-modality in several parameters of galaxies mainly related to distribution functions and star formation, such as the one we observe in the SFMS and sSFR–M* relation currently (see Fig. 1). The basic idea is that star formation in galaxies is supported by the accretion of gas from the surroundings, but above ∼1012 M in halo mass (∼1010.5 M in M*) the halo becomes virialised and the infalling gas is shock-heated as it approaches the galaxy, dramatically slowing down the accretion and conversely the SFR in galaxies.

This is supported by results from simulations (Birnboim & Dekel 2003; Kereš et al. 2005; van de Voort & Schaye 2012) and the bi-modality is explained as a transition in the dominant mechanism that regulates star formation inside the galaxies. Below log M*/M ∼ 10.5, stellar winds mainly from supernovae (SNe) suppress star formation in an essentially self-regulating mechanism which becomes less efficient for higher masses as the growing gravitational potential makes it harder for the winds to expel gas from the galaxy (e.g. Silk 1997; Hopkins et al. 2011, 2014). Conversely, AGN feedback starts to become efficient at about the same M* (∼1010.5 M) (Dekel & Birnboim 2006) as, statistically, the central black hole mass grows with that of the host galaxy. The outflow efficiency bimodality has been clearly detected and attributed to these two mechanisms (SNe and AGN) in simulations (Mitchell et al. 2020).

In our results we observe the bimodality clearly in the current properties, but we also see the signature of the golden mass in the fact that at precisely this mass value we see a bend in the relation for most LBT values. This bend appears even for currently star-forming galaxies and is observed even at the highest LBT we consider, albeit to a lesser degree than at the latest LBT such as 1 and 3 Gyr (see Fig. 2).

The difference in the extent to which the bi-modality is observed in our results could be due to either a loss of resolution in our results for earlier times, which is expected, or a sign that the virialisation of the halo is a process which takes time to have an effect on the galaxies. Overall, our results largely support a scenario where galaxies are initially self-regulated such that their sSFR is similar, regardless of their M* until AGN feedback quenches galaxies above the golden mass. The previously discussed feature of a negative correlation for low M* galaxies is, however, at odds with the interpretation that below ∼1010.5 M in M*, galaxies are dominated by stellar feedback. This feedback is expected to be especially efficient for very low M* galaxies with shallow potential wells (e.g. Veilleux et al. 2005; Bland-Hawthorn et al. 2007; Mitchell et al. 2020) and therefore we would expect lower, not higher, sSFR for them. As mentioned above; to obtain a complete picture of this feature, we need to obtain measurements of the actual gas present in the galaxies as well as confirm the trend with deeper observations.

6.2.2. Morphological quenching

To collapse and form stars, the gas must cool and dissipate kinematic energy (especially turbulence). Consequently, any mechanism that adds turbulence to the ISM can suppress the ability of clouds to form stars and thus reduce the efficiency of star formation. AGN and stellar feedback can directly introduce turbulence into the ISM via injection of kinetic energy, but galaxy dynamics can also create it through shear.

The most important of these mechanisms is the effect of the bulge on the disc, which can transfer a large amount of kinetic energy as turbulence and stabilise the disc to suppress star formation as it grows: the so-called morphological quenching (Martig et al. 2009). Bars can be considered as a kind of extension of the bulge into the disc and are expected to help amplify turbulence in the ISM around them. Apart from this, simulations show that the bars are able to funnel gas out of the disc into the centre of the galaxy, promoting star formation there, but ultimately also contributing to extinction by accelerating the time scale for the depletion of the gas (Binney et al. 1991; Athanassoula 1992; Boone et al. 2007; Sormani et al. 2015).

Our results are very much in line with the idea that quenching and the way it is implemented in galaxies correlates with morphological type. We find differences in the number of galaxies that form stars, the extent to which the sSFR has decreased over cosmic time, and the timing of quenching.

6.2.3. Galaxy interactions

Mergers and interactions can have dramatic effects on the evolution of galaxies and will therefore be a key factor for it, especially in hierarchical galaxy formation models where today’s galaxies are the result of many merger processes (e.g. White & Rees 1978; White & Frenk 1991; Cole et al. 2000).

Unfortunately, our results are largely blind to the effects of mergers because we measure the history of galaxies in terms of their current stellar composition. Any object in our sample that resulted from a previous merger of two galaxies has an SFH that is the sum of the two galaxies at any point in time. As a result, any such galaxy will appear in our figures at both a higher M* and a higher SFR, keeping it within the global relation but extending it to higher M*. In general, we do not expect this to produce significant differences in our results, with one possible exception. Major mergers between galaxies tend to produce a galaxy that is more early type than its progenitors (e.g. Mihos & Hernquist 1996; Lotz et al. 2008). Additionally, given that a merged galaxy will be measured at higher M* in the past than it should be (as two smaller galaxies), mergers are a possible explanation for the M* segregation with morphology observed at very high LBT. The extent to which the segregation is due to mergers cannot be readily deduced from our data. Sánchez et al. (2019) measures the cosmic density of SFR at different redshifts using similar methodology to ours and finds that mergers do not significantly affect the results.

7. Summary and conclusions

In this paper, we have studied how the current correlation between M* and sSFR came to be by measuring both parameters at a different LBT using full spectral fitting. We have also studied whether and how the current star formation and morphology of the galaxies is tied to their past sSFR history. Lastly, we compared our results with several cosmic surveys to see if we can find an agreement with our findings. Our main conclusions are:

  • Galaxies in the past used to have similar mass-doubling times regardless of M*. Earlier-type and more massive galaxies then start to drop from the sequence around 7 Gyr ago in our results, producing the currently observed downsizing.

  • Galaxies that end up in different morphology bins in the present are segregated in the sSFR–M* relation even at very early times. They are mainly separated in terms of their M*, but earlier-type galaxies also have a slightly higher sSFR high z; we measured a small (∼0.1 dex) yet significant (median p-value for all M* < 10−5) offset between early- and late-type galaxies.

  • The existence of a high-M* inflection point is similar to predictions from simulations, though the z at which it occurs is much later in our results, likely due to methodology limitations.

  • The differences between morphology bins agree well with predictions for morphological quenching.

  • Comparing to cosmic surveys, we find agreement in several of them for some of the features we describe, though again the z values differ significantly.

These results show how SFHs derived through full spectral fitting can be used as a complementary study to cosmic surveys to probe cosmic star formation if its caveats are considered. The specific advantages of the method allow for a consistent measurement of the differences in the growth rate of distinct bins of galaxies. This allowed us to see how the morphology of the galaxies is determined by their past history starting from > 10 Gyr ago.


Acknowledgments

We thank the referee for useful comments which improved the robustness of the results shown here. ACF and PSB acknowledge financial support by the Spanish Ministry of Science and Innovation through the research grant PID2019-107427-GB-31 and PID2021-123417OB-I00. RMM acknowledges support from Spanish Ministerio de Ciencia e Innovación MCIN/AEI/10.13039/501100011033 through grant PGC2018-093499-B-I00 and MDM-2017-0737 Unidad de Excelencia “Maria de Maeztu”-Centro de Astrobiología (CAB), INTA-CSIC, as well as by “ERDF A way of making Europe”. RMM acknowledges the support from the Instituto Nacional de Técnica Aeroespacial SHARDSJWST project through the PRE-SHARDSJWST/2020 PhD fellowship, and by “ESF Investing in your future”. SFS acknowledges funding from PAPIIT-DGAPA-AG100622 (UNAM) project. SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, Center for Astrophysics | Harvard & Smithsonian (CfA), the Chilean Participation Group, the French Participation Group, Instituto de Astrofísica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe (IPMU)/University of Tokyo, the Korean Participation Group, Lawrence Berkeley National Laboratory, Leibniz Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astrophysik (MPA Garching), Max-Planck-Institut für Extraterrestrische Physik (MPE), National Astronomical Observatories of China, New Mexico State University, New York University, University of Notre Dame, Observatório Nacional/MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Autónoma de México, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University.

References

  1. Asari, N. V., Cid Fernandes, R., Stasińska, G., et al. 2007, MNRAS, 381, 263 [Google Scholar]
  2. Atek, H., Furtak, L. J., Oesch, P., et al. 2022, MNRAS, 511, 4464 [CrossRef] [Google Scholar]
  3. Athanassoula, E. 1992, MNRAS, 259, 345 [Google Scholar]
  4. Barrera-Ballesteros, J. K., Utomo, D., Bolatto, A. D., et al. 2020, MNRAS, 492, 2651 [NASA ADS] [CrossRef] [Google Scholar]
  5. Barro, G., Pérez-González, P. G., Cava, A., et al. 2019, ApJS, 243, 22 [NASA ADS] [CrossRef] [Google Scholar]
  6. Behroozi, P., Wechsler, R. H., Hearin, A. P., & Conroy, C. 2019, MNRAS, 488, 3143 [NASA ADS] [CrossRef] [Google Scholar]
  7. Belfiore, F., Maiolino, R., Maraston, C., et al. 2017, MNRAS, 466, 2570 [Google Scholar]
  8. Belfiore, F., Maiolino, R., Bundy, K., et al. 2018, MNRAS, 477, 3014 [Google Scholar]
  9. Binney, J., Gerhard, O. E., Stark, A. A., Bally, J., & Uchida, K. I. 1991, MNRAS, 252, 210 [Google Scholar]
  10. Birnboim, Y., & Dekel, A. 2003, MNRAS, 345, 349 [Google Scholar]
  11. Bland-Hawthorn, J., Veilleux, S., & Cecil, G. 2007, Ap&SS, 311, 87 [NASA ADS] [CrossRef] [Google Scholar]
  12. Blanton, M. R., & Moustakas, J. 2009, ARA&A, 47, 159 [Google Scholar]
  13. Boone, F., Baker, A. J., Schinnerer, E., et al. 2007, A&A, 471, 113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Boselli, A., Boissier, S., Cortese, L., et al. 2006, ApJ, 651, 811 [Google Scholar]
  15. Boselli, A., Voyer, E., Boissier, S., et al. 2014, A&A, 570, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151 [Google Scholar]
  17. Bryant, J. J., Owers, M. S., Robotham, A. S. G., et al. 2015, MNRAS, 447, 2857 [Google Scholar]
  18. Bundy, K., Ellis, R. S., & Conselice, C. J. 2005, ApJ, 625, 621 [NASA ADS] [CrossRef] [Google Scholar]
  19. Bundy, K., Bershady, M. A., Law, D. R., et al. 2015, ApJ, 798, 7 [Google Scholar]
  20. Camps-Fariña, A., Sanchez, S. F., Lacerda, E. A. D., et al. 2021, MNRAS, 504, 3478 [CrossRef] [Google Scholar]
  21. Camps-Fariña, A., Sánchez, S. F., Mejía-Narváez, A., et al. 2022, ApJ, 933, 44 [CrossRef] [Google Scholar]
  22. Camps-Fariña, A., Sánchez-Blázquez, P., Roca-Fàbrega, S., & Sánchez, S. F. 2023, A&A, 678, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Cano-Díaz, M., Sánchez, S. F., Zibetti, S., et al. 2016, ApJ, 821, L26 [Google Scholar]
  24. Caplar, N., & Tacchella, S. 2019, MNRAS, 487, 3845 [NASA ADS] [CrossRef] [Google Scholar]
  25. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
  26. Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&A, 562, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Cid Fernandes, R., Pérez, E., García Benito, R., et al. 2013, A&A, 557, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Cid Fernandes, R., González Delgado, R. M., García Benito, R., et al. 2014, A&A, 561, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, MNRAS, 319, 168 [Google Scholar]
  30. Conroy, C. 2013, ARA&A, 51, 393 [NASA ADS] [CrossRef] [Google Scholar]
  31. Conselice, C. J. 2014, ARA&A, 52, 291 [CrossRef] [Google Scholar]
  32. Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 [Google Scholar]
  33. Crain, R. A., Schaye, J., Bower, R. G., et al. 2015, MNRAS, 450, 1937 [NASA ADS] [CrossRef] [Google Scholar]
  34. Daddi, E., Dickinson, M., Morrison, G., et al. 2007, ApJ, 670, 156 [NASA ADS] [CrossRef] [Google Scholar]
  35. Davé, R., Thompson, R., & Hopkins, P. F. 2016, MNRAS, 462, 3265 [Google Scholar]
  36. Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 [NASA ADS] [CrossRef] [Google Scholar]
  37. Dekel, A., Sarkar, K. C., Jiang, F., et al. 2019, MNRAS, 488, 4753 [CrossRef] [Google Scholar]
  38. Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168 [Google Scholar]
  39. Dickinson, M., Bergeron, J., Casertano, S., et al. 2003, Great Observatories Origins Deep Survey (GOODS) Validation Observations, 196 [Google Scholar]
  40. Domínguez, A., Siana, B., Brooks, A. M., et al. 2015, MNRAS, 451, 839 [CrossRef] [Google Scholar]
  41. Drory, N., MacDonald, N., Bershady, M. A., et al. 2015, AJ, 149, 77 [CrossRef] [Google Scholar]
  42. Elbaz, D., Dickinson, M., Hwang, H. S., et al. 2011, A&A, 533, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Erb, D. K., Shapley, A. E., Pettini, M., et al. 2006, ApJ, 644, 813 [NASA ADS] [CrossRef] [Google Scholar]
  44. Fontana, A., Pozzetti, L., Donnarumma, I., et al. 2004, A&A, 424, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Galametz, A., Grazian, A., Fontana, A., et al. 2013, ApJS, 206, 10 [Google Scholar]
  46. Gallazzi, A., Charlot, S., Brinchmann, J., White, S. D. M., & Tremonti, C. A. 2005, MNRAS, 362, 41 [Google Scholar]
  47. Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Sci. Rev., 123, 485 [Google Scholar]
  48. Gaspari, M., Brighenti, F., & Temi, P. 2012, MNRAS, 424, 190 [NASA ADS] [CrossRef] [Google Scholar]
  49. Genel, S., Vogelsberger, M., Springel, V., et al. 2014, MNRAS, 445, 175 [Google Scholar]
  50. Giavalisco, M., Ferguson, H. C., Koekemoer, A. M., et al. 2004, ApJ, 600, L93 [NASA ADS] [CrossRef] [Google Scholar]
  51. Goddard, D., Thomas, D., Maraston, C., et al. 2017, MNRAS, 466, 4731 [NASA ADS] [Google Scholar]
  52. González Delgado, R. M., Cid Fernandes, R., Pérez, E., et al. 2016, A&A, 590, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 197, 35 [NASA ADS] [CrossRef] [Google Scholar]
  54. Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332 [NASA ADS] [CrossRef] [Google Scholar]
  55. Guo, Y., Ferguson, H. C., Giavalisco, M., et al. 2013, ApJS, 207, 24 [NASA ADS] [CrossRef] [Google Scholar]
  56. Hopkins, P. F., Quataert, E., & Murray, N. 2011, MNRAS, 417, 950 [NASA ADS] [CrossRef] [Google Scholar]
  57. Hopkins, P. F., Kereš, D., Oñorbe, J., et al. 2014, MNRAS, 445, 581 [NASA ADS] [CrossRef] [Google Scholar]
  58. Ibarra-Medel, H. J., Sánchez, S. F., Avila-Reese, V., et al. 2016, MNRAS, 463, 2799 [NASA ADS] [CrossRef] [Google Scholar]
  59. Ilbert, O., Capak, P., Salvato, M., et al. 2009, ApJ, 690, 1236 [NASA ADS] [CrossRef] [Google Scholar]
  60. Ilbert, O., Arnouts, S., Le Floc’h, E., et al. 2015, A&A, 579, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Karim, A., Schinnerer, E., Martínez-Sansigre, A., et al. 2011, ApJ, 730, 61 [Google Scholar]
  62. Kennicutt, R. C. 1989, ApJ, 344, 685 [Google Scholar]
  63. Kereš, D., Katz, N., Weinberg, D. H., & Davé, R. 2005, MNRAS, 363, 2 [Google Scholar]
  64. Kuntschner, H., Lucey, J. R., Smith, R. J., Hudson, M. J., & Davies, R. L. 2001, MNRAS, 323, 615 [NASA ADS] [CrossRef] [Google Scholar]
  65. Lacerda, E. A. D., Sánchez, S. F., Cid Fernandes, R., et al. 2020, MNRAS, 492, 3073 [Google Scholar]
  66. Lacerda, E. A. D., Sánchez, S. F., Mejía-Narváez, A., et al. 2022, New Astron., 97, 101895 [NASA ADS] [CrossRef] [Google Scholar]
  67. Law, D. R., Cherinka, B., Yan, R., et al. 2016, AJ, 152, 83 [Google Scholar]
  68. Lee, J. C., Kennicutt, R. C., Funes, S. J. J. G., Sakai, S., & Akiyama, S. 2009, ApJ, 692, 1305 [NASA ADS] [CrossRef] [Google Scholar]
  69. Lin, L., Pan, H.-A., Ellison, S. L., et al. 2019, ApJ, 884, L33 [NASA ADS] [CrossRef] [Google Scholar]
  70. Lockman, F. J., Jahoda, K., & McCammon, D. 1986, ApJ, 302, 432 [Google Scholar]
  71. Lotz, J. M., Davis, M., Faber, S. M., et al. 2008, ApJ, 672, 177 [NASA ADS] [CrossRef] [Google Scholar]
  72. Lotz, J. M., Koekemoer, A., Coe, D., et al. 2017, ApJ, 837, 97 [Google Scholar]
  73. Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 [Google Scholar]
  74. Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Martig, M., Bournaud, F., Teyssier, R., & Dekel, A. 2009, ApJ, 707, 250 [Google Scholar]
  76. Matthee, J., & Schaye, J. 2019, MNRAS, 484, 915 [CrossRef] [Google Scholar]
  77. Mérida, R. M., Pérez-González, P. G., Sánchez-Blázquez, P., et al. 2023, ApJ, 950, 125 [CrossRef] [Google Scholar]
  78. Mihos, J. C., & Hernquist, L. 1996, ApJ, 464, 641 [Google Scholar]
  79. Mitchell, P. D., Schaye, J., Bower, R. G., & Crain, R. A. 2020, MNRAS, 494, 3971 [Google Scholar]
  80. Nayyeri, H., Hemmati, S., Mobasher, B., et al. 2017, ApJS, 228, 7 [NASA ADS] [CrossRef] [Google Scholar]
  81. Noeske, K. G., Weiner, B. J., Faber, S. M., et al. 2007, ApJ, 660, L43 [CrossRef] [Google Scholar]
  82. Pérez-González, P. G., Rieke, G. H., Villar, V., et al. 2008a, ApJ, 675, 234 [Google Scholar]
  83. Pérez-González, P. G., Trujillo, I., Barro, G., et al. 2008b, ApJ, 687, 50 [CrossRef] [Google Scholar]
  84. Pérez-González, P. G., Cava, A., Barro, G., et al. 2013, ApJ, 762, 46 [Google Scholar]
  85. Pozzetti, L., Bolzonella, M., Zucca, E., et al. 2010, A&A, 523, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Rodríguez-Puebla, A., Primack, J. R., Avila-Reese, V., & Faber, S. M. 2017, MNRAS, 470, 651 [Google Scholar]
  87. Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267 [NASA ADS] [CrossRef] [Google Scholar]
  88. Salmon, B., Papovich, C., Finkelstein, S. L., et al. 2015, ApJ, 799, 183 [NASA ADS] [CrossRef] [Google Scholar]
  89. Sánchez Almeida, J. 2017, Astrophys. Space Sci. Lib., 430, 67 [CrossRef] [Google Scholar]
  90. Sánchez Almeida, J., Elmegreen, B. G., Muñoz-Tuñón, C., & Elmegreen, D. M. 2014, A&ARv, 22, 71 [Google Scholar]
  91. Sánchez, S. F. 2006, Astron. Nachr., 327, 850 [CrossRef] [Google Scholar]
  92. Sánchez, S. F. 2020, ARA&A, 58, 99 [Google Scholar]
  93. Sánchez, S. F., Kennicutt, R. C., Gil de Paz, A., et al. 2012, A&A, 538, A8 [Google Scholar]
  94. Sánchez, S. F., Pérez, E., Sánchez-Blázquez, P., et al. 2016a, Rev. Mex. Astron. Astrofis., 52, 21 [NASA ADS] [Google Scholar]
  95. Sánchez, S. F., Pérez, E., Sánchez-Blázquez, P., et al. 2016b, Rev. Mex. Astron. Astrofis., 52, 171 [Google Scholar]
  96. Sánchez, S. F., Avila-Reese, V., Hernandez-Toledo, H., et al. 2018, Rev. Mex. Astron. Astrofis., 54, 217 [Google Scholar]
  97. Sánchez, S. F., Avila-Reese, V., Rodríguez-Puebla, A., et al. 2019, MNRAS, 482, 1557 [Google Scholar]
  98. Sánchez, S. F., Walcher, C. J., Lopez-Cobá, C., et al. 2021a, Rev. Mex. Astron. Astrofis., 57, 3 [Google Scholar]
  99. Sánchez, S. F., Barrera-Ballesteros, J. K., Colombo, D., et al. 2021b, MNRAS, 503, 1615 [CrossRef] [Google Scholar]
  100. Sánchez, S. F., Barrera-Ballesteros, J. K., Lacerda, E., et al. 2022, ApJS, 262, 36 [CrossRef] [Google Scholar]
  101. Sánchez-Blázquez, P., Gorgas, J., Cardiel, N., & González, J. J. 2006a, A&A, 457, 787 [CrossRef] [EDP Sciences] [Google Scholar]
  102. Sánchez-Blázquez, P., Gorgas, J., Cardiel, N., & González, J. J. 2006b, A&A, 457, 809 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Sánchez-Blázquez, P., Ocvirk, P., Gibson, B. K., Pérez, I., & Peletier, R. F. 2011, MNRAS, 415, 709 [Google Scholar]
  104. Sánchez-Blázquez, P., Rosales-Ortega, F. F., Méndez-Abreu, J., et al. 2014, A&A, 570, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Santini, P., Fontana, A., Castellano, M., et al. 2017, ApJ, 847, 76 [NASA ADS] [CrossRef] [Google Scholar]
  106. Savaglio, S., Glazebrook, K., Le Borgne, D., et al. 2005, ApJ, 635, 260 [NASA ADS] [CrossRef] [Google Scholar]
  107. Schawinski, K., Urry, C. M., Simmons, B. D., et al. 2014, MNRAS, 440, 889 [Google Scholar]
  108. Schaye, J., Dalla Vecchia, C., Booth, C. M., et al. 2010, MNRAS, 402, 1536 [Google Scholar]
  109. Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521 [Google Scholar]
  110. Schmidt, M. 1959, ApJ, 129, 243 [NASA ADS] [CrossRef] [Google Scholar]
  111. Schreiber, C., Pannella, M., Elbaz, D., et al. 2015, A&A, 575, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Scoville, N., Aussel, H., Brusa, M., et al. 2007, ApJS, 172, 1 [Google Scholar]
  113. Silk, J. 1997, ApJ, 481, 703 [CrossRef] [Google Scholar]
  114. Smee, S. A., Gunn, J. E., Uomoto, A., et al. 2013, AJ, 146, 32 [Google Scholar]
  115. Smethurst, R. J., Lintott, C. J., Simmons, B. D., et al. 2015, MNRAS, 450, 435 [NASA ADS] [CrossRef] [Google Scholar]
  116. Sormani, M. C., Binney, J., & Magorrian, J. 2015, MNRAS, 449, 2421 [NASA ADS] [CrossRef] [Google Scholar]
  117. Speagle, J. S., Steinhardt, C. L., Capak, P. L., & Silverman, J. D. 2014, ApJS, 214, 15 [Google Scholar]
  118. Stark, D. P., Ellis, R. S., Bunker, A., et al. 2009, ApJ, 697, 1493 [NASA ADS] [CrossRef] [Google Scholar]
  119. Stefanon, M., Yan, H., Mobasher, B., et al. 2017, ApJS, 229, 32 [NASA ADS] [CrossRef] [Google Scholar]
  120. Straatman, C. M. S., Spitler, L. R., Quadri, R. F., et al. 2016, ApJ, 830, 51 [NASA ADS] [CrossRef] [Google Scholar]
  121. Sturm, E., González-Alfonso, E., Veilleux, S., et al. 2011, ApJ, 733, L16 [NASA ADS] [CrossRef] [Google Scholar]
  122. Szokoly, G. P., Bergeron, J., Hasinger, G., et al. 2004, ApJS, 155, 271 [NASA ADS] [CrossRef] [Google Scholar]
  123. Tomczak, A. R., Quadri, R. F., Tran, K.-V. H., et al. 2016, ApJ, 817, 118 [Google Scholar]
  124. Trager, S. C., Faber, S. M., Worthey, G., & González, J. J. 2000, AJ, 120, 165 [NASA ADS] [CrossRef] [Google Scholar]
  125. Trussler, J., Maiolino, R., Maraston, C., et al. 2020, MNRAS, 491, 5406 [Google Scholar]
  126. van de Voort, F., & Schaye, J. 2012, MNRAS, 423, 2991 [CrossRef] [Google Scholar]
  127. Vázquez-Mata, J. A., Hernández-Toledo, H. M., Avila-Reese, V., et al. 2022, MNRAS, 512, 2222 [CrossRef] [Google Scholar]
  128. Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 769 [NASA ADS] [CrossRef] [Google Scholar]
  129. Vogelsberger, M., Genel, S., Springel, V., et al. 2014, MNRAS, 444, 1518 [Google Scholar]
  130. Wake, D. A., Bundy, K., Diamond-Stanic, A. M., et al. 2017, AJ, 154, 86 [Google Scholar]
  131. Walcher, J., Groves, B., Budavári, T., & Dale, D. 2011, Ap&SS, 331, 1 [NASA ADS] [CrossRef] [Google Scholar]
  132. Wetzel, A. R., Tinker, J. L., & Conroy, C. 2012, MNRAS, 424, 232 [NASA ADS] [CrossRef] [Google Scholar]
  133. White, S. D. M., & Frenk, C. S. 1991, ApJ, 379, 52 [Google Scholar]
  134. White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 [Google Scholar]
  135. Williams, R. E., Blacker, B., Dickinson, M., et al. 1996, AJ, 112, 1335 [NASA ADS] [CrossRef] [Google Scholar]
  136. Wong, T., Xue, R., Bolatto, A. D., et al. 2013, ApJ, 777, L4 [NASA ADS] [CrossRef] [Google Scholar]
  137. Worthey, G. 1999, ASP Conf. Ser., 192, 283 [NASA ADS] [Google Scholar]
  138. Yan, R., Tremonti, C., Bershady, M. A., et al. 2016, AJ, 151, 8 [Google Scholar]
  139. York, D. G., Adelman, J., Anderson, J. E., et al. 2000, AJ, 120, 1579 [Google Scholar]

Appendix A: Comparison between current SFR measured in Hα and SSPs

As a base check for the reliability of the measurements using full spectral fitting and for comparison purposes we show here the current sSFR-M relation using the SFR measured from the Hα emission line, which can be found in Fig. A.1. The general shape of the relation is very similar to that of Fig. 1, though the sSFR values reach significantly lower values for retired galaxies.

thumbnail Fig. A.1.

Currently observed sSFR of MaNGA galaxies compared to their M* (bottom left), with the SFR measured using the Hα emission line. The black contours enclose 95%, 65% and 35% of the sample, respectively. The green line is the median sSFR within 0.25 dex wide M* bins. In the top and right panels, the distributions in M* and sSFR, respectively, are shown as histograms, with the three morphological bins shown in lines of the corresponding colour.

In Fig. A.2 we show a direct comparison between the SFR measured using the Hα emission line and the value from the stellar populations used in Fig. 1. At relatively high values the two match very well, but for very low SFR the SSP values appear to have a floor. This is not a critical difference as it becomes relevant at around 0.01 M yr−1 which is a very low amount of star formation. The higher floor of the SSP values is likely due to method related reasons, such as the fact that, while the Hα line can be measured in absorption, the minimum light fraction of a stellar population is 0, thus noise can only increase the SFR measurement when the actual value is 0.

thumbnail Fig. A.2.

Comparison between the current SFR in our sample measured using the Hα emission line (X-axis) and full spectral fitting selecting the last 30 Myr (Y axis). The black contours show the general distribution and the green line the 1:1 relation.

Appendix B: Impact of the FoV coverage in MaNGA

Our measurements of the M* and the SFR of the galaxies are intrinsically limited by the coverage that the MaNGA fibre bundles achieve for each galaxy. The MaNGA full sample is defined such that the primary sample, which makes up 2/3 of the galaxies, has ∼80% of the galaxies covered out to 1.5 Re. The secondary sample, with 1/3 of the galaxies, conversely has 80% of them covered out to 2.5 Re. Naturally, these are only selection targets and most of the galaxies will be covered to further galactocentric distances than these.

The resulting loss of light will impact our results, lowering the values of the SFR and M* but the sSFR will be less impacted by this. Our results are mainly based on differences between bins of galaxies selected using M*, morphology and present-day star formation and, therefore, the main concern is that galaxies in different bins will have significantly different loss of light due to the FoV.

To assess this, we have used the Sérsic index measurements from the NASA-Sloan Atlas Catalog and the FoV coverage measurements from Sánchez et al. (2022) to measure the average fraction of light that falls within the FoV by assuming the galaxies can be generally approximated using a Sérsic profile. The resulting average fraction for each of the bins considered is shown in Fig. B.1. It bears mention that the M* range does not reach values as low as those seen in Figs. 2 and 3, due to the fact that the FoV is measured only in the present and, thus, the M* range corresponds better to that of Figs. 1 or A.1.

thumbnail Fig. B.1.

Average fraction of light contained within the FoV for each of the M*, morphology (top), and star formation (bottom) bins used in this article.

The figure shows that, in general, the average loss of light due to the FoV is small (0-20%) and that it is fairly independent of M*, though there is a general trend to lower coverage in more massive galaxies. There are, however, offsets between the morphology and star formation bins we use, such that earlier-type and retired galaxies have lower coverage. Galaxies typically have age gradients such that the stellar populations are older in their centres and therefore the SFR values at high LBT will be less affected by the FoV if the FoV covers most of the galaxy. In any case, the higher loss of light for retired and earlier-type galaxies would mean a relative underestimation of the SFR at high LBT, but in Fig. 2 we show currently retired galaxies as being more star-forming, on average, than those that form stars in the present. Since the sSFR-M* relation is quite flat the offset in M* would also not significantly impact the results.

In Fig. 3 we show that earlier-type galaxies are located at slightly higher sSFR values and significantly higher M* values comparatively. Similarly to Fig. 2, the relative loss of flux in earlier-type galaxies means that these should be shifted to even higher M* values and possibly to higher sSFR values. As such, we conclude that the FoV coverage in MaNGA can have some effect on the measurements, but these effects cannot be the origin of the results we report.

All Figures

thumbnail Fig. 1.

Currently observed sSFR of MaNGA galaxies compared to their M* (bottom left). The black contours enclose 95%, 65% and 35% of the sample respectively. The green line is the median sSFR within 0.25 dex wide M* bins. In the top and right panels, the distributions in M* and sSFR, respectively, are shown as histograms, with the three morphological bins shown in lines of the corresponding colour.

In the text
thumbnail Fig. 2.

Evolution of the sSFR over cosmic time for the galaxies in the MaNGA sample. In each panel the M* and sSFR of the galaxies measured at different LBT are shown. The colour of the data points corresponds to the currently observed morphology of the galaxies. The black contours enclose 95%, 65% and 35% of the sample, respectively. The cyan and violet dots and lines show the median sSFR within 0.5 dex wide M* bins, selecting galaxies above and below 6 Å in EWHα, respectively. The error bars for the dots show the 25th and 75th percentile of the distribution within each M* bin.

In the text
thumbnail Fig. 3.

Evolution of the sSFR over cosmic time for the galaxies in the MaNGA sample. In each panel the M* and sSFR of the galaxies measured at different LBT are shown. There are three sets of contours whose colours correspond to the current morphology of the galaxies and whose levels enclose 95% and 50% of each sub-sample.

In the text
thumbnail Fig. 4.

Two-sample Kolmogorov-Smirnov test for each pair of galaxy populations considered in the article. In each panel we show the p-value for each pair of populations performed individually at each of the 0.5 dex wide M* bins. The dotted line shows the limit of p-value = 0.05, below which the samples are considered to be statistically significantly different. The colours identify the pairs of populations used.

In the text
thumbnail Fig. 5.

sSFR–M* relation for eight cosmic surveys. The data have been adapted into the M*–sSFR plane in some cases. The colour of each line corresponds to the average redshift of its bin and the redshift colour scale is the same for all panels. The dotted lines correspond to M* intervals where sample completeness is not guaranteed as reported in each work. The lower right corner of each field indicates the median scatter of the data.

In the text
thumbnail Fig. A.1.

Currently observed sSFR of MaNGA galaxies compared to their M* (bottom left), with the SFR measured using the Hα emission line. The black contours enclose 95%, 65% and 35% of the sample, respectively. The green line is the median sSFR within 0.25 dex wide M* bins. In the top and right panels, the distributions in M* and sSFR, respectively, are shown as histograms, with the three morphological bins shown in lines of the corresponding colour.

In the text
thumbnail Fig. A.2.

Comparison between the current SFR in our sample measured using the Hα emission line (X-axis) and full spectral fitting selecting the last 30 Myr (Y axis). The black contours show the general distribution and the green line the 1:1 relation.

In the text
thumbnail Fig. B.1.

Average fraction of light contained within the FoV for each of the M*, morphology (top), and star formation (bottom) bins used in this article.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.