Identification of galaxies that experienced a recent major drop of star formation

[abridged] Building upon a previous study, we define a method to blindly identify galaxies that underwent, and may still be undergoing, a fast downfall of their star-formation activity, that is, a more than 80% drop in star-formation rate (SFR) occurring in less than 500 Myr. Modeling galaxies' SED with a delayed-$\tau$ star formation history (SFH), with and without allowing an instantaneous SFR drop within the last hundreds Myr, we isolate 102 candidates out of a subsample of 6,680 galaxies classified as star-forming from the UVJ criterion in the ZFOURGE catalogues. These galaxies are mostly located in the lower part of the SFR-M$_*$ main sequence (MS) and extend up to a factor 100 below it. They also lie close to the limit between the passive and active regions on the UVJ diagram, indicating that they are in a transition phase. We show that the selected candidates have different physical properties compared to galaxies with similar UVJ colors, namely, lower star-formation rates and different stellar masses. The morphology of the candidates show no preference for a particular type. Among the 102 candidates, only 4 show signs of an AGN activity (from X-ray luminosity or UV-IR SED fitting decomposition). This low fraction of AGNs among the candidates implies that AGN activity may not be the main driver of the recent downfall, although timescale differences and duty cycle must be taken into account. We finally attempt to recover the past position of these galaxies on the SFR-M$_*$ plane, before the downfall of their star-formation and show that some of them were in the starburst region before and are back on the MS. These candidates constitute a promising sample that need more investigation in order to understand the different mechanisms at the origin of the star formation decrease of the Universe since $z$$\sim$2.


Introduction
Ten years ago, the discovery of the tight relation linking the star formation rate (SFR) and stellar mass of star-forming galaxies opened a new window in our understanding of galaxy evolution (Elbaz et al. 2007;Noeske et al. 2007). The main consequence of this relation is that galaxies are forming the bulk of their stars through steady state processes rather than violent episodes of star formation. This galaxy main sequence (MS) is found to hold up to z=4 (Schreiber et al. 2017) and its normalisation and shape to vary with redshift Pannella et al. 2009;Elbaz et al. 2011;Rodighiero et al. 2011;Speagle et al. 2014;Whitaker et al. 2014;Schreiber et al. 2015;Gavazzi et al. 2015;Tomczak et al. 2016). However, what is striking is that the scatter of the main sequence is found to be relatively constant at all masses and over cosmic time (Guo et al. 2013;Ilbert et al. 2015;Schreiber et al. 2015). A possible origin of this scatter is that it could be artificially created by the accumulation of errors in the extraction of photometric measurements and/or in the determination of the SFR and/or stellar mass in relation with model uncertainties. However, several studies have found a coherent variation of physical galaxy properties such as the gas fraction (Magdis et al. 2012), Sersic index and effective radius (Wuyts et al. 2011), U-V color (e.g., Salmi et al. 2012) suggesting that the bulk of the scatter is related to physics and not measurement and model uncertainties. Oscillations of the SFR resulting from a varying infall rate and compaction of star-formation have also been advocated to explain the MS scatter (Tacchella et al. 2016).
The path of galaxies on the SFR-M * plane is still debated. Do galaxies evolve along the MS, growing in mass? Do they reach at some point the starburst region and then go back on the MS or quench? Do they undergo small variations going above and/or below the MS? Some simulations predict that galaxies undergo some fluctuations of star formation activity resulting in variations of their SFR such as compaction or variations of accretion (e.g., Dekel & Burkert 2014;Sargent et al. 2014;Scoville et al. 2017). These variations must be small enough to keep the SFR of the galaxy within the MS scatter. From an observational point of view, Elbaz et al. (2017) showed that some massive compact galaxies exhibiting starburst galaxy properties (short depletion time and high infrared surface density) can be found within the MS. However they have different morphologies and gas fraction compared to "true" starbursts (above the MS), indicating a different origin, possibly being late-stage mergers of gas-rich galaxies. This could be a clue of a possible recent movement of these galaxies from the starburst galaxies region back to the MS. However, recent studies pointed toward alternative interpretations to the in situ scenario implied by the MS based on the young ages found for starburst galaxies (da Cunha et al. 2015;Ma et al. 2015;Mancuso et al. 2016). Thus, there is no consensus on the evolution of galaxies relative to the MS.
To tackle this problem, one has to study the recent star formation history of galaxies. This information is embedded in their spectral energy distribution (SED). However, recovering it through SED modeling is complex and subject to many uncertainties and degeneracies. Although an average SFH of galaxies can be derived assuming that they follow the MS Ciesla et al. 2017), galaxies are expected to exhibit complex SFHs, with short term fluctuations, needing sophisticated SFH parametrizations to model them (e.g., Lee et al. 2010;Pacifici et al. 2013;Behroozi et al. 2013;Pacifici et al. 2016). However, these models are complex to implement and a large library is needed to be able to model all galaxies properties. Instead, numerous studies used simple analytical forms to model galaxies SFH (e.g., Papovich et al. 2001;Maraston et al. 2010;Pforr et al. 2012;Gladders et al. 2013;Simha et al. 2014;Buat et al. 2014;Boquien et al. 2014;Ciesla et al. 2015;Abramson et al. 2016;Ciesla et al. 2016Ciesla et al. , 2017. Furthermore, SFH parameters are known to be difficult to constrain from broadband SED modeling with a general agreement on the difficulty to constrain the age of the galaxy, here defined as the age of the oldest star, from broad-band SED fitting (e.g. Maraston et al. 2010;Pforr et al. 2012;Buat et al. 2014;Ciesla et al. 2015Ciesla et al. , 2017. To understand the origin of the scatter of the MS, we need to use an analytical SFH that must able to recover recent variations of the SFR with a precision better than the scatter of the MS itself, i.e. 0.3 dex. Recently, Ciesla et al. (2017) showed that a delayed SFH to which we add a flexibility on the recent SFH provides SFR that are more accurate than those estimated by other typical analytical SFHs (τ models, delayed, etc...). This SFH was tested in Ciesla et al. (2016). Studying a sample of local galaxies from the Virgo cluster, that is galaxies known to have undergone a fast drop of star formation activity due to ram pressure stripping, we showed that the amplitude of the flexibility can be constrained by broadband SED modeling as long as UV rest frame and near-IR data are available.
Relying on these results, we propose a new approach to probe the recent star formation history of galaxies and identify those that underwent a recent (≤500 Myr) and major (≥80%) drop of their star-formation activity. The goal of our study is to pinpoint recent movements in the SFR-M * plane. In Sect. 2, we present the selected sample of galaxies drawn from the ZFOURGE survey. In Sect. 3, we detail the SED modeling code and procedure used in this work. The definition of the criterion used to select galaxies that underwent a recent and major downfall of their SFR is presented in Sect. 4. The properties of the selected candidates are presented and studied in Sect. 5, Sect. 6, and Sect. 7. The possibility of probing movements inside the star-forming main sequence is discussed in Sect. 8.
Throughout the paper, we use the WMAP7 cosmological parameters (H 0 = 70.4 km Mpc −1 s −1 , Ω 0 = 0.272, Komatsu et al. 2011). SED fitting is performed assuming an IMF of Salpeter (1955), but the results of this work are found to be robust against IMF choice.

The sample
The aim of this work is to identify galaxies experiencing a recent and fast decrease of their star-formation activity. In Ciesla et al. (2016), we showed that the method we use in this study is sensitive to rapid downfall of star-formation that occurred in less than ∼500 Myr. Since the downfall considered here happens in a very short time, we need a large statistics to be able to catch the responsible physical processes in the act. Another key to perform this study is to have a good and large sampling of the UV-NIR SED of galaxies. For these two reasons, we draw our sample from the ZFOURGE/GOODS-South catalogue (Straatman et al. 2016) which contain intermediate bands observations exactly in the key spectral domain where we expect the SED to be sensitive to one star formation history (SFH) or the other.
From the ZFOURGE/GOODS-South parent sample, we select galaxies with a redshift comprised between 0.5 and 3.5 and a stellar mass higher than 10 9 M (as published in Kriek et al. 2009). This corresponds to a final sample of 7,471 galaxies for which we will model the UV to IR SED, and among which 6,680 galaxies are classified as star-forming following a UVJ color selection. The photometric bands used for the fits are listed in Table 1. The photometric catalogues were extensively described in Straatman et al. (2016), and references therein, we refer the reader to these papers for further information on the photometry and catalogue building. We impose of a minimum of three for the signal-to-noise ratio of the flux of each photometric band. In the rest of the paper, we will use the U-V and V-J colors provided by Straatman et al. (2016) as we used them to select our sample. However, the stellar masses (and SFR) used in this work result from the SED fitting presented in Sect. 3.2.

Broadband SED modeling
3.1. The CIGALE code CIGALE 1 (Code Investigating GALaxy Emission) is a SED modeling software package that has two functions: a SED modeling function and a SED fitting function (Roehlly et al. 2014, Boquien et al., in prep). Even though the philosophy of CIGALE, originally presented in Noll et al. (2009), remains, the code has been rewritten in Python and additions have been made in order to optimize its performance and broaden its scientific applications. The SED modeling function of CIGALE allows the building of galaxy SEDs from UV to radio by assuming a combination of modules which model the star formation history (SFH) of the galaxy, the stellar emission from stellar population models (Bruzual & Charlot 2003;Maraston 2005), the nebular lines, the attenuation by dust (e.g., Calzetti et al. 2000), the IR emission from dust (Draine & Li 2007;Casey 2012;Dale et al. 2014), the AGN contribution (Fritz et al. 2006;Ciesla et al. 2015), and the radio emission. CIGALE builds galaxy SEDs taking into account the balance between the energy absorbed by dust in UVoptical and reemitted in the IR.
With the SED fitting function of CIGALE, these modeled SEDs are then integrated into a set of filters to be compared directly to the observations. For each parameter, a probability distribution function (PDF) analysis is made. The output value is the likelihood-weighted mean value of the PDF and the error associated is the likelihood-weighted standard deviation. We use CIGALE to derive the physical properties of galaxies such as stellar masses, instantaneous SFRs, dust attenuation, SFH characteristics, taking into account panchromatic information from the SED. To perform the SED fitting on the ZFOURGE/GOODS-South sample, we rely on photometric redshifts. For the sample used in this work, that is considering the redshift and stellar mass selections, the precision of these redshifts is ∼3%.
In CIGALE, the assumed SFH can be handled in two different ways. The first is to model it using simple analytic functions such as exponential forms or delayed-τ SFH (see for instance, Ciesla et al. 2016Ciesla et al. , 2017. The second is to provide more complex SFHs, such as those provided by semi-analytical models (SAM) and simulations Ciesla et al. 2015), an approach which we will not use in the present study.

Star formation history modeling
To identify galaxies that underwent recent important variations in their SFH from broadband SED fitting, we use the SFH implemented in CIGALE by Ciesla et al. (2016) to model such a rapid decrease of the SFR: where t trunc is the time at which the star formation is instantaneously affected, and r SFR is the ratio between SFR(t > t trunc ) and SFR(t = t trunc ): .
(2) Fig. 1 shows an example of the truncated SFH (see Fig. 3 of Ciesla et al. 2016, to see the impact of each parameter on the shape of the SED). This SFH has been presented and tested intensively in Ciesla et al. (2016) where we showed that the age trunc parameter (age trunc = age − t trunc ) cannot be constrained from broadband SED fitting. However, the r SFR parameter is well-constrained especially in galaxies that were relatively active before the SFH variation. Comprised between 0 et 1, this parameter handles the strength of the rapid star-formation downfall in the SFH. Thus, we empirically showed in Ciesla et al. (2016) that a star-forming galaxy with an output value of r SFR lower than 0.2-0.3 experienced a rapid decrease of its star formation activity in the last ∼500 Myr. However, when allowing r SFR values larger than 1, it is possible to add a recent burst in the star formation activity of the galaxies modeling starburst galaxies. Thus, with this flexible SFH we are able to model a large variety of recent SFH (Ciesla et al. 2017).

SED fitting procedure
To model galaxy SEDs, we used the stellar population models of Bruzual & Charlot (2003) et al. (2000). The energy attenuated in UV-optical is re-emitted in the IR (mostly 8 to 1000 µm) through the models of Dale et al. (2014). However, 71% of the galaxies are not detected byeither Spitzer-MIPS or Herschel-PACS/SPIRE. To apply the same procedure with the same constraints to all galaxies, and to prevent the fit quality to be too strongly driven by the IR, we only fit the UV-to-NIR data of the galaxies, i.e. up to Spitzer/IRAC4. We thus fix the shape of the dust emission since we do no study the dust properties of these galaxies.

Constraints on the output parameters: mock analysis
To understand what are the output parameters of the SED fitting procedure that we can trust and use in our study, i.e. what are the parameters that are constrained with the photometric coverage available for the ZFOURGE sample, we perform a mock analysis. This procedure is thoroughly presented in several works of the literature (e.g., Giovannoli et al. 2011). We briefly describe here the main steps. We first apply our SED fitting procedure on the whole real catalogue in order to retrieve the best fit SED for each individual source. Each best fit SED is then convolved with the set of filters available for the ZFOURGE sample and a noise is added to the mock flux, randomly chosen in a Gaussian distribution with a σ equal to the observed error. The original error of each filter is associated to the perturbed flux. We thus have a mock photometric catalogue for which we know all the associated physical parameters from the output of the first run. A more detailed description on the mock tests implemented in CIGALE can be found in (Giovannoli et al. 2011) and (Boquien et al. 2012). Then, we apply our SED modeling procedure a second time, but on the mock catalogue and compare the output results to the true values of each parameter we want to use in our study. These comparisons are shown in Fig. 2 for the stellar mass, SFR, and r SFR . As expected given the good coverage of the restframe NIR, the stellar mass is well recovered and thus reliable.
Although showing a larger dispersion for values between 0.1 and 10 M /yr, the SFR is relatively well constrained with a relation close to the one-to-one line. The results of the mock analysis for r SFR show that the parameter is, on average, always overestimated. We discuss below that this is not affecting our analysis and the identification of galaxies using this parameter since this overestimation will just make our chosen threshold even more conservative (see Section 4).

Selection criterion
We consider only galaxies defined as star-forming from the UVJ diagram criterion of Whitaker et al. (2011) (see also Tomczak et al. 2014;Spitler et al. 2014;Straatman et al. 2014). Star formation histories of sources that underwent a rapid decrease of their SFR (≤500 Myr) are not well fitted by models assuming a normal delayed-τ SFH, or a exponentially-decreasing one, whereas the truncated SFH proposed manages to reproduce the UV-NIR SED (Ciesla et al. 2016(Ciesla et al. , 2017. We thus expect that the χ 2 obtained with a truncated SFH are significantly lower than the ones from the delayed-τ SFH for galaxies with a recent decrease of SFR as shown in Ciesla et al. (2016). In Fig. 3, we show the distribution of χ 2 obtained using both SFH. However, the two SFH have different number of degree of freedom. To take this into account we compute the Bayesian Information criterion (BIC) for each SFH, using the definition: with χ 2 the non-reduced χ 2 of the fit, k the number of degree of freedom, and n the number of photometric bands fitted. We then calculate ∆BIC the difference between BIC delayed and BIC trunc . If ∆BIC is larger than 2, then there is a significant difference between the two fits (Liddle 2004;Stanley et al. 2015), with the truncated SFH providing a better fit of the data than the delayed-τSFH. However, we impose a ∆BIC≥10 to be conservative. The second criterion is based on the value of r SFR . Indeed, from the study of Virgo cluster galaxies undergoing ram pressure stripping from the intra-cluster medium, we showed in Ciesla et al. (2016) that these galaxies, that underwent a rapid downfall of SFR, have an r SFR of 0.3 and lower. Here, we consider that galaxies with r SFR ≤0.2 are possible candidates, to be more conservative. Although calibrated from a study of galaxies wellknown for undergoing a fast drop of their star formation, the valueof 0.2 is mostly arbitrary. However, the goal of this study is not to select a statistically complete sample to derive properties but rather to use broadband SED fitting to identify rapid phenomena. This threshold translates in a decrease of the SFR of at least 80% in the last few hundreds Myr according to Ciesla et al. (2016) where we tested the sensitivity of the truncated SFH in terms of the strength of the star-formation downfall and its timescale. However, the overestimation of the true r SFR value revealed by the mock analysis (Fig. 2) does not affect our criterion but rather goes in the sense of strengthening it since an output value of 0.2 would be an overestimation of the true one. Finally, we impose a last criterion to reject galaxies with bad fits for both SFH: we reject the 8% galaxies with the largest χ 2 trunc to remove galaxies with catastrophic fits. Quantitatively, we exclude galaxies with χ 2 trunc >88. From these three criteria, we identify galaxies as candidates if: The result of this selection is shown in Fig. 4 where we see that the large majority of the galaxies (98.5%) lie in the region where ∆BIC<10 or r SFR >0.2. We thus identify 102 sources satisfying our criterion. The distribution of the χ 2 obtained with the two SFH are shown in the inset panel of Fig. 3. The best fits of six randomly selected sources, three normal and three candidate galaxies having experienced a rapid downfall of SFR, are presented in Fig. 5. The three normal galaxies are displayed in the top panels where we see that the best fit SEDs are identical whether the normal delayed-τ or the truncated SFH is used. Their best fit SFHs are similar using either of the two SFHs, except for #562 where the truncated SFH prefers the contribution of a recent burst to reproduce the SED whereas the delayed-τ shows a younger age to reach the final SFR. The three other SEDs provide examples of fits results for three candidates. For these galaxies, the delayed-τ SFH does not manage to fit both the UV and the NIR or reproduce the Balmer break. Comparison between the true parameters of the mock SEDs and the results from the SED modeling with CIGALE. A parameter is well-constrained when there is a one-to-one relation between the two. Statistics on the ratio between the output and the input parameters are provided for the stellar mass and SFR: the first quartile of the ratio, the median, and the third quartile. The red dashed line in the r SFR panel indicates the threshold adopted in this study. Models built from a truncated SFH manage to reconcile the UV and NIR emission and provide good fit to the observed fluxes. For the three candidates the best truncated SFH clearly shows a rapid and drastic decrease. For the delayed-τ SFH, in order to provide a equivalent final SFR, the age is shortened or the main peak of SFR is weaken.

The case study of three candidates
To understand the lower quality of delayed-τ SFH fits of compared to the truncated SFH, we choose to study in detail the three candidates, #20946, #20989, and #29927, which SEDs are shown in Fig. 5 and RGB images in Fig. 6. When running CIGALE, we are limited in the number of models to be fitted, and thus on the number of priors, due to computational reasons. The size of the prior grids, if reasonable, does not affect our results on the estimate of the physical parameters since we derive them through their probability distribution function (if the minimum and maximum values are well-chosen and the grid contains a significant number of values). However, it affects the value of the χ 2 and thus can impact our selection criterion as the ∆BIC is estimated using the χ 2 of the fit for each model. Thus, we check if ∆BIC would be the same for these three galaxies if we used an increased grid of τ, age, and E(B-V), and if they would still fulfill the criterion. We run CIGALE using 100 values of τ (from 0.1 to 20 Gyr), of age (from 1 Myr to 13 Gyr), and of E(B-V) (from 0.01 to 1.5), successively. For each test, we provide in Table 3 the resulting χ 2 and ∆BIC for both the truncated and normal delayed-τ SFHs. In each of these tests, the χ 2 obtained with the truncated SFH is the same with variation of less than 5%. For the normal delayed-τ SFH, the χ 2 is improved up to 64% in the case of an increased age grid for #29927. However, despite these improvements in the normal delayed-τ χ 2 , these galaxies still fulfill the selection criterion with a ∆BIC larger than 10. Thus the difference we observe in the quality of the fits using the two SFH is not due to the prior lists that we use for our main run. The UV should be the first spectral range to be impacted by recent variations in the SFH. Thus we test if fitting the SED of our three test galaxies without the UV rest-frame (<300 nm) fluxes would lead to a better agreement between the two SFHs. The runs are made with the increased range of age priors as we just showed that it improves the χ 2 obtained with the normal delayed-τ SFH. The ∆BIC obtained for these runs are negative with an absolute value lower than 10 indicating that both SFH provide an equivalent fit of the data (Table 3). The UV is thus the key domain that the delayed-τ SFH struggles to reconcile with  the optical-NIR part of the SED. However, as shown in Fig. 5, both the delayed-τ and truncated SFHs provide a good fit of the UV emission of #20989. Here, the discrepancy comes from the optical/NIR emission. One can then imagine that with the addition of an older population, that is considering two stellar populations, it would be possible to fit the NIR emission without the need for a recent fast decrease of the SFR. To test this possibility, we fit the SED of #20989 using a double exponentially-declining SFH, available in CIGALE, that allows that mimics two populations combined through an extra free parameter. The χ 2 obtained with this SFH is 82.4 versus 65.2 for the truncated SFH although the two exponentially-declining SFH has more free parameters. Thus a more complicated SFH is not preferred in this case. Interestingly, as shown in Fig. 5 the double exponentiallydeclining SFH favors a main stellar population with quasi constant evolution plus a recent burst of SF followed by a rapid decline. The same test was applied to all the candidates and none of them was significantly better fit by a two population SFH. Thus, we conclude that the candidates underwent rapid and recent variations in their SF, probed by the UV, difficult to be modeled by a smooth SFH. We show RGB color images of the three test candidates in Fig. 6. The sources #20946 (z=0.61) and #20989 (z=0.66) have a red disk and exhibit some UV emission as lanes. Their morphology suggest a merger scenario. In Fig. 7, we place them on the SFR-M * plane along with the MS at their respective redshift. Both #20946 and #20989 lied above the MS prior to theirdecrease of SFR. Using the information provided by the SED fitting on the strength of the SFR decrease and knowing that it happened approximately in less than 500 Myr (the exact time is not constrained from broad-band SED fitting, see Ciesla et al. (2016) and Boselli et al. (2016)), we can recover their position on this plane before the downfall of SFR. They were thus previously even higher above the MS by a factor of 6.7 and 33.3 for #20946 and #20989, respectively, comforting the merger scenario. The question is then to know whether they will later remain on the MS or continue to fall and become passive. Fast decreases of SFR of these orders are consistent with the simulated SFH of mergers computed by Di Matteo et al. (2005). Interestingly, even the morphologies of #20946 and #20989 with a red disk and UV lane across the disk are consistent with the snapshots of their simulation. Since #20946 and #20989 are relatively low mass galaxies, below 1.2×10 10 M , we do not expect these mergers to have a strong dust enshrouded phase, allowing us to identify them from UV-to-NIR broadband SED modeling. The case of #29927 is slightly different as the galaxy is now below the MS. Previously it was indeed above, but, went down in a few hundreds Myr. This source is more massive with logM * =10.51 and at higher redshift (z=1.15). Its morphology (Fig. 6, right panel) shows a large red disk with a central bulge and a UV ring in the middle of the disk. There is also another source lying close to #29927 that could be an indication of possible interactions.

Physical properties of the candidate galaxies
To investigate the properties of these 102 sources, we place them on the main sequence diagram (Fig. 8). They are mostly located in the lower part of the MS indicating a possible transition in their star formation activity. This is consistent with the hypothesis that these sources, selected from their SED, are undergoing a diminution of their SFR. Moreover, we observe that the candidates span a large range of stellar masses, from ∼10 9 M to masses higher than 10 11 M . Adding the information on the r SFR parameter allows us to see that the candidates located in the lower part of the MS still have some residual star formation showing that they are not totally quenched. Indeed, our criterion imposes a drastic diminution of about 80% of their star formation activity in the last few hundred Myr and not necessarily a total quenching. We then display the candidate sources on the UVJ diagram. We recall that our candidates arepre-selected to lie in the star-forming region of this diagram, whichexplains why none of them lies in the quiescent region. However, they are located in a region just below the passive region, close to the line separating the two zones, supporting the hypothesis that these galaxies are undergoing or just underwent a decrease of star formation activity.
Recent studies have shown that galaxies displayed in color diagrams, such as UVJ or NUVrK, show a gradient of sSFR decreasing from the active region toward the passive region (Arnouts et al. 2013). We thus check if our selection is equivalent to a sSFR selection. To test this point, we separate our selected sample in three slices in the UVJ diagram, parallel to the line separating the active and passive region (Fig. 9, top panel). If they were selected from their sSFR value, we would expect that these three slices will be recovered in the MS diagram. But, as shown in Fig. 9 (bottom panel), we recover no gradient of sSFR in the MS diagram, the galaxies are dispersed. This means that our criterion is not just a selection in sSFR and that we are able to select galaxies with different SFH than the galaxies located in the same zone in the UVJ diagram. Indeed, as discussed for the three test sources in Sect. 4.2, the UV is a key domain to probe recent variations in the SFH and "normal" analytical SFHs fail to reconcile the UV and NIR colors of our candidates. But the UV is not used in the UVJ diagram explaining why we can find galaxies with identical UVJ colors but different SFHs.

Reference sample
The position of the candidate galaxies on the UVJ diagram show that they are at the limit between the passive and the active galaxies regions that could imply that these candidates may be undergoing some SFR downfall. Yet, other galaxies are located in the same regions of the UVJ diagram and do not show any specific SED shape. To understand the difference between the candidates and the others, we define a reference sample. For each candidate, we select the galaxies that are located in the UVJ diagram at the same position and in the same redshift bin, within a circle defined as: (5) where (V−J) and (U−V) are the colours of the reference galaxies and (V − J) q and (U −V) q those of the candidate. The 0.2 value is chosen arbitrarily after some tuning to ensure enough statistics to compute average properties of the reference sample. Tests have been made to ensure that the results of this study do not depend on reasonable variation of the adopted radius.

Physical parameters
We compare the physical properties of the candidates to those of the reference sample. We average the properties of the galaxies of each reference sample and compare the results to those of the corresponding candidate. In Fig. 10    from its average properties, is systematically different than the position of the corresponding candidate galaxy. As expected, the SFR of the candidates is systematically lower than the average SFR of normal galaxies located in the same place of the UVJ diagram.

Morphology
With our method, we have adopted an empirical criterion to select galaxies that recently underwent a drastic decrease of their SFR and showed that their physical properties are different from those of galaxies with identical (U-V) and (V-J) colors. We now look at the morphology of the candidates. The morphological classification of Huertas-Company et al. (2016) using deep machine learning is used to separate the candidates between spheroids (24 galaxies), pure disks (31 galaxies), disk plus spheroids (4 galaxies), irregular disks (20 galaxies), and mergers or irregular systems (12 galaxies).
We show in Fig. 11 the candidates in the MS and the UVJ plots, color coded as a function of their morphological classification. There is no segregation of one particular type in a specific stellar mass range, spheroids candidates are found at both low and high masses as well as disks. However, at z>2, the candidates are almost all spheroids or mergers, with only three disks out of 24 galaxies. Comparing the morphological types of the candidates with those of a reference sample is not straightforward. Instead, we look at their effective radius, R e computed by van der Wel et al. (2014). As shown in Fig. 12, the candidates and the galaxies of their reference samples cover the same range of effective radius (see inset panel of Fig. 12). The median effective radius are 2.3 kpc and 2.5 kpc for the reference samples galaxies and the candidates, respectively. The first and third quartiles of both distribution are also consistent: 1.3 and 3.3 kpc for the reference samples galaxies and 1.8 and 3.4 kpc for the candidates. Therefore, we conclude that they have effective radius that are consistent with those of galaxies of similar stellar masses.

AGN activity
A possible source of rapid downfall of star-formation could come from AGN activity as some simulations predict the AGNdriven outflows to be responsible for the downfall of star formation activity (e.g., Baldry et al. 2004;Hopkins et al. 2006) although other results suggest than these outflows hardly affect the gas content of galaxies (Gabor & Bournaud 2014). Out of the 102 candidate galaxies, 5 are detected in X-ray (Cappelluti et al. 2016) and only 4 of them have an X-ray luminosity above 10 42 erg/s: #1817, #3463, #10032, and #30250. The 4 sources with L X >10 42 erg/s are massive with a stellar mass equal or higher than 1.7×10 10 M .
Furthermore, we test whether some of these sources could host an enshrouded AGN emitting in IR. This can only beexplored for sources detected in the MIR/FIR, we thus select the 47 candidate sources detected at MIPS 24 µm. Using the AGN emission modeling of the CIGALE code, we estimate for each source the need for an AGN template to fit the UV to IR SED. Out of the 47 sources, 23 have an output f rac AGN parameter (contribution of the AGN to the total IR luminosity) below 10%. Based on the results of Ciesla et al. (2015), we can conclude that these sources do not likely host a strong IR AGN. Twentythree sources have a trustable f rac AGN parameter above 15%. However, among these 23 likely-AGN-galaxies, 18 of them have only one IR flux on which the decomposition can rely, which  (Schreiber et al. 2015). Bottom panels: Galaxies displayed on the UVJ diagram. Grey dots are the galaxies of the sample well fitted by both SFH models. Colored dots are the sources selected, color-coded according to the r SFR parameter. The black solid lines indicate the UVJ selection according to Whitaker et al. (2011).
is not enough to provide a fiducial conclusion on the presence of an enshrouded AGN in these sources. Five candidates, however, have a high f rac AGN between 30 and 50%: #1415, #10534, #15220, #17144, and #18949. Out of these five sources, the SED of four of them indicates an unusually high IR 24 µm or 100 µm flux compared to the expected ones from the energy absorbed in UV-NIR (CIGALE performs energy balance), suggesting a possible blending problem in the IR photometry. This is confirmed by the fact that all of these four sources have a close companion. Only #18949 shows a clear sign of possible AGN contribution in MIR. Out of the 102 candidates, only 5 of them show a sign of AGN activity either in X-ray or in IR. If it was to be interpreted in terms of duty cycle, considering that our method is sensitive to variations occurring less than 500 Myr prior to observation, it would imply an AGN duty cycle of 5%. Assuming that all the candidates went through an AGN-on phase, it would lead to an AGN visibility timescale of 25 Myr. However, given the poor statistic of AGN detected candidates, we conclude that AGN might not be the main cause for the recent star-formation downfall underwent by the candidate sources.

Probing recent movements relative to the main sequence
The epoch when the fast decrease of star formation occurred remains unconstrained from broadband photometry alone. However, the SFR that the candidate galaxies experienced before the decrease is constrained by the r SFR parameter. This information allows to locate the position of the galaxies with respect to the star-formation MS prior to the downfall of SFR, i.e., to deter-  Fig. 13). The candidates fall in three categories discussed below. The first group is composed of sources that were lying in the starburst region, i.e. above the MS, before undergoing the decrease of star formation activity and are now on the MS (Fig. 14). These sources are thus post-starburst candidates. They have stellar masses lower than 2×10 10 M . Only one source is massive and is the most massive of the sample, #3463. According to the morphological classification of Huertas-Company et al. (2016), 56% of the galaxies of this group are disks. This is almost a factor three more than what is found for galaxies located in the same region of the SFR-M * (see Table 4). The remaining 44% are classified as irregular disks or systems versus 57% found for "normal galaxies". We notice the absence of spheroids in this candidate group when we expect about 23% of the "normal" galaxies to have this morphological type. With the data in hand, we cannot state if these candidates are still in transition and thus just passing through the MS and will reach the passive region afterwards or if, after their star formation burst, they are back on the MS and will stay in this position for a given time. Going further in this analysis would require, at least, spectral information on Hα for instance, to compare different star formation indicators sensitive to different timescales. The second group gathers candidates that were in the starburst region and that are now below the MS (Fig. 15). These sources, which have stellar masses spanning a larger range, passed through the MS in a relative short time which is of the order of several hundred Myr. The distribution in terms of morphological types in this group is more balanced as we find 30% of spheroid galaxies, 30% disks or disk+spheroids, and 40% of irregular disks or system, versus 24%, 21%, and 55% for the reference galaxies, respectively (Table 4).
The third group is composed of galaxies that were on the MS before and are now below it (Fig. 16). These galaxies do not span a specific stellar mass range (Fig. 13). It is composed of a large fraction of galaxies classified as spheroids (36%) as well as pure disks (39%). The median stellar mass of the galaxies classified as spheroids is 2.3×10 10 M whereas the median mass of disks is 6.8×10 9 M . This seems to indicate that the massive galaxies experienced a morphological transformation within the MS, prior to their rapid downfall. This morphological transformation may have been associated with the slow downfall associated to the bending of the MS (Schreiber et al. 2015). For the low mass galaxies, the fact that they are classified as disks is an indication that they have some rotation, despite the fact that being well below the MS. The fraction of irregular systems and disks is less important, 25% versus 44% and 40% for group 1 and 2, respectively, and is almost a factor of two below the proportion of irregular galaxies in the reference sample of this group (Table 4). We find galaxies with low stellar mass (≤10 10 M ) that departed the MS recently and lie now below it.
We find no significant behavior that could be different from one group to another based on the properties derived by our broad-band SED fitting procedure. We showed that the statistical properties of our candidates point toward a recent decrease of star-formation activity in this galaxies. However, we need spectroscopic data to definitely validate the method and its use to quickly identify such galaxies in the act. Furthermore, to understand the different mechanisms that are affecting the candidate galaxies of each group we need additional data in order to look for different star formation indicators, AGN activity, and gas content. Our method can be used to pinpoint galaxies with recent variations of SFH and try to reconstruct these variations, first step to understand the evolution of galaxies in the SFR-M * plane. Although this study focuses on galaxies undergoing a decrease of star formation activity, it could also be used, in principle, to identify movements on the MS plane due to an increase of star formation. We find 233 candidates with r SFR >1 and ∆BIC≥10 that will be explored in a future work.

Conclusions
To probe and study phenomena that occur in a very short amount of time, we need a large statistic to be able to catch them in the act. The goal of this work is to identify galaxies that just underwent a rapid (<500 Myr) and drastic downfall of their starformation rate (more than 80%) from broadband SED modeling since large photometric sample can provide the statistic needed to pinpoint these objects. In this purpose, we choose to apply our method on the ZFOURGE sample to provide both statistics and large photometric coverage of the SED from UV to NIR.
-We fit the whole ZFOURGE sample using the CIGALE code with two different SFHs: one normal delayed-τ SFH and one truncated SFH, and using Bruzual & Charlot (2003) models. We show that a handful of sources are better fit using the truncated SFH, that assume a recent instantaneous break in the SFH, compared to the more commonly used delayed-τ SFH. -Based on the results of Ciesla et al. (2016), we define a criterion based on the quality of the fits and the output strength of the decrease of star formation activity obtained from the SED modeling. From galaxies classified as star-forming with the UVJ diagram, we select a sample of 102 candidates. -We study in details three candidates and show that the difference of fit quality between the two SFH is not due to the priors used for the SED modeling and that the normal delayed-τ struggles to model consistently with the UV rest-frame and the optical-NIR emission. -Placed on the main sequence diagram, these sources lie between the lower part of the MS and the passive galaxies region. In the UVJ diagram, they lie in a region in between the passive and active regions. This confirms the fact that these sources must be in transition.  -To understand their differences, we define for each candidate a reference sample composed of galaxies lying in the same region in the UVJ diagram and compare their physical properties. -We compute the average properties of the reference samples associated to each candidate source. We show that despite their common location on the UVJ diagram, the candidate galaxies are displaced in the MS diagram compared to their reference sample. They have, on average, a lower SFR and an offset in stellar mass that seems to depend on the stellar mass of the candidate. -There is no specific morphology attributed to the candidates.
All types are found. Their effective radius is consistent with what is expected for galaxies of the same stellar masses. -Among the 102 candidates, 4 of them are detected in X-ray with a luminosity larger than 10 42 erg/s. Another source is classified as AGN from an UV-to-NIR SED fitting decompo-sition. This low number imply that AGN activity is not the main responsible for the fast decrease of the star-formation activity of the candidates. -Using the output of the SED fitting, we reconstruct the recent variation of SFR and show that the candidate galaxies are separated in three groups. The first one is composed of galaxies that were in the starburst regions and are now on the MS, which are thus post-starburst candidates. The second group is composed of galaxies that went from above to below the MS. The last group is composed of galaxies that were on the MS and are now below.
In this study, we provided and used a method to blindly identify galaxies undergoing a fast diminution of their star-formation activity from their observed broadband SED. We showed that these galaxies have indeed physical properties different from the main population and characteristic of sources in transition. However, additional work is needed to confirm the variations in the recent Norm. counts Fig. 12: Effective radius of galaxies with redshift between 0.5 and 1 as a function of their stellar mass. Grey dots are all the galaxies in this redshift range. Pink triangles are the candidates of this redshift range, and the blue circle are all the galaxies lying in the same part of the SFR-M * plane than the candidates, that is the galaxies used to build the reference samples. The inset panel shows the distributions of R e in kpc for both the candidates and the galaxies associated to the reference samples. SFH of these galaxies, especially using spectroscopy to look for strong Balmer absorption lines. Ancillary data are also needed to investigate the mechanisms at play and understand the possible role of mergers, AGN, and gas content.