scholarly journals Constraining the recent star formation history of galaxies: an approximate Bayesian computation approach

2020 ◽  
Vol 635 ◽  
pp. A136 ◽  
Author(s):  
G. Aufort ◽  
L. Ciesla ◽  
P. Pudlo ◽  
V. Buat

Although galaxies are found to follow a tight relation between their star formation rate and stellar mass, they are expected to exhibit complex star formation histories (SFH) with short-term fluctuations. The goal of this pilot study is to present a method that identifies galaxies that undergo strong variation in star formation activity in the last ten to some hundred million years. In other words, the proposed method determines whether a variation in the last few hundred million years of the SFH is needed to properly model the spectral energy distribution (SED) rather than a smooth normal SFH. To do so, we analyzed a sample of COSMOS galaxies with 0.5 <  z <  1 and log M* >  8.5 using high signal-to-noise ratio broadband photometry. We applied approximate Bayesian computation, a custom statistical method for performing model choice, which is associated with machine-learning algorithms to provide the probability that a flexible SFH is preferred based on the observed flux density ratios of galaxies. We present the method and test it on a sample of simulated SEDs. The input information fed to the algorithm is a set of broadband UV to NIR (rest-frame) flux ratios for each galaxy. The choice of using colors is made to remove any difficulty linked to normalization when classification algorithms are used. The method has an error rate of 21% in recovering the correct SFH and is sensitive to SFR variations larger than 1 dex. A more traditional SED-fitting method using CIGALE is tested to achieve the same goal, based on fit comparisons through the Bayesian information criterion, but the best error rate we obtained is higher, 28%. We applied our new method to the COSMOS galaxies sample. The stellar mass distribution of galaxies with a strong to decisive evidence against the smooth delayed-τ SFH peaks at lower M* than for galaxies where the smooth delayed-τ SFH is preferred. We discuss the fact that this result does not come from any bias due to our training. Finally, we argue that flexible SFHs are needed to be able to cover the largest possible SFR-M* parameter space.

2020 ◽  
Vol 498 (4) ◽  
pp. 5581-5603
Author(s):  
Sabine Bellstedt ◽  
Aaron S G Robotham ◽  
Simon P Driver ◽  
Jessica E Thorne ◽  
Luke J M Davies ◽  
...  

ABSTRACT We apply the spectral energy distribution (SED) fitting code ProSpect to multiwavelength imaging for ∼7000 galaxies from the GAMA survey at z &lt; 0.06, in order to extract their star formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We show with this approach that we are able to recover the observationally determined cosmic star formation history (CSFH), an indication that stars are being formed in the correct epoch of the Universe, on average, for the manner in which we are conducting SED fitting. We also show the contribution to the CSFH of galaxies of different present-day visual morphologies and stellar masses. Our analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago. In other morphological types, the stellar mass formed later, up to 6 Gyr ago for present-day irregular galaxies. Similarly, the most massive galaxies in our sample were shown to have formed half their stellar mass by 11 Gyr ago, whereas the least massive galaxies reached this stage as late as 4 Gyr ago (the well-known effect of ‘galaxy downsizing’). Finally, our metallicity approach allows us to follow the average evolution in gas-phase metallicity for populations of galaxies and extract the evolution of the cosmic metal mass density in stars and in gas, producing results in broad agreement with independent, higher redshift observations of metal densities in the Universe.


2020 ◽  
Vol 497 (4) ◽  
pp. 4753-4772
Author(s):  
Shuang Zhou ◽  
H J Mo ◽  
Cheng Li ◽  
Médéric Boquien ◽  
Graziano Rossi

ABSTRACT We measure the star formation histories (SFHs) of a sample of low-mass galaxies with M* &lt; 109 M⊙ from the SDSS-IV MaNGA survey. The large number of IFU spectra for each galaxy are either combined to reach a high signal to noise ratio or used to investigate spatial variations. We use Bayesian inferences based on full spectrum fitting. Our analysis based on Bayesian evidence ratio indicates a strong preference for a model that allows the presence of an old stellar population, and that an improper model for the SFH can significantly underestimate the old population in these galaxies. The addition of near-infrared photometry to the constraining data can further distinguish between different SFH model families and significantly tighten the constraints on the mass fraction in the old population. On average more than half of the stellar mass in present-day low-mass galaxies formed at least 8 Gyr ago, while about 30 per cent within the past 4 Gyr. Satellite galaxies on average have formed their stellar mass earlier than central galaxies. The radial dependence of the SFH is quite weak. Our results suggest that most of the low-mass galaxies have an early episode of active star formation that produces a large fraction of their present stellar mass.


2018 ◽  
Vol 615 ◽  
pp. A55 ◽  
Author(s):  
Eduardo Telles ◽  
Jorge Melnick

Aims. We present a UV to mid-IR spectral energy distribution (SED) study of a large sample of SDSS DR13 HII galaxies. These galaxies are selected as starbursts (EW(Hα) > 50Å) and for their high-excitation locus in the upper left region of the BPT diagram. Their photometry was derived from the cross-matched GALEX, SDSS, UKDISS, and WISE catalogs. Methods. We used CIGALE modeling and a SED fitting routine with the parametrization of a three-burst star formation history, and a comprehensive analysis of all other model parameters. We were able to estimate the contribution of the underlying old stellar population to the observed equivalent width of Hβ, and allow for more accurate burst age determination. Results. We found that the star formation histories of HII Galaxies can be reproduced remarkably well by three major eras of star formation. In addition, the SED fitting results indicate that in all cases the current burst produces a small percent of the total stellar mass, i.e., the bulk of stellar mass in HII galaxies has been produced by the past episodes of star formation, and also indicate that at a given age the Hβ luminosity depends only on the mass of young stars favoring a universal IMF for massive stars. Most importantly, the current star formation episodes are maximum starbursts that produce stars at the highest possible rate.


Genetics ◽  
2021 ◽  
Author(s):  
Guillaume Laval ◽  
Etienne Patin ◽  
Pierre Boutillier ◽  
Lluis Quintana-Murci

Abstract During their dispersals over the last 100,000 years, modern humans have been exposed to a large variety of environments, resulting in genetic adaptation. While genome-wide scans for the footprints of positive Darwinian selection have increased knowledge of genes and functions potentially involved in human local adaptation, they have globally produced evidence of a limited contribution of selective sweeps in humans. Conversely, studies based on machine learning algorithms suggest that recent sweeps from standing variation are widespread in humans, an observation that has been recently questioned. Here, we sought to formally quantify the number of recent selective sweeps in humans, by leveraging approximate Bayesian computation and whole-genome sequence data. Our computer simulations revealed suitable ABC estimations, regardless of the frequency of the selected alleles at the onset of selection and the completion of sweeps. Under a model of recent selection from standing variation, we inferred that an average of 68 (from 56 to 79) and 140 (from 94 to 198) sweeps occurred over the last 100,000 years of human history, in African and Eurasian populations, respectively. The former estimation is compatible with human adaptation rates estimated since divergence with chimps, and reveal numbers of sweeps per generation per site in the range of values estimated in Drosophila. Our results confirm the rarity of selective sweeps in humans and show a low contribution of sweeps from standing variation to recent human adaptation.


Author(s):  
Cecilia Viscardi ◽  
Michele Boreale ◽  
Fabio Corradi

AbstractWe consider the problem of sample degeneracy in Approximate Bayesian Computation. It arises when proposed values of the parameters, once given as input to the generative model, rarely lead to simulations resembling the observed data and are hence discarded. Such “poor” parameter proposals do not contribute at all to the representation of the parameter’s posterior distribution. This leads to a very large number of required simulations and/or a waste of computational resources, as well as to distortions in the computed posterior distribution. To mitigate this problem, we propose an algorithm, referred to as the Large Deviations Weighted Approximate Bayesian Computation algorithm, where, via Sanov’s Theorem, strictly positive weights are computed for all proposed parameters, thus avoiding the rejection step altogether. In order to derive a computable asymptotic approximation from Sanov’s result, we adopt the information theoretic “method of types” formulation of the method of Large Deviations, thus restricting our attention to models for i.i.d. discrete random variables. Finally, we experimentally evaluate our method through a proof-of-concept implementation.


2021 ◽  
Vol 62 (2) ◽  
Author(s):  
Jason D. Christopher ◽  
Olga A. Doronina ◽  
Dan Petrykowski ◽  
Torrey R. S. Hayden ◽  
Caelan Lapointe ◽  
...  

2020 ◽  
Vol 501 (2) ◽  
pp. 1568-1590
Author(s):  
Lukas J Furtak ◽  
Hakim Atek ◽  
Matthew D Lehnert ◽  
Jacopo Chevallard ◽  
Stéphane Charlot

ABSTRACT We present new measurements of the very low mass end of the galaxy stellar mass function (GSMF) at z ∼ 6−7 computed from a rest-frame ultraviolet selected sample of dropout galaxies. These galaxies lie behind the six Hubble Frontier Field clusters and are all gravitationally magnified. Using deep Spitzer/IRAC and Hubble Space Telescope imaging, we derive stellar masses by fitting galaxy spectral energy distributions and explore the impact of different model assumptions and parameter degeneracies on the resulting GSMF. Our sample probes stellar masses down to $M_{\star }\gt 10^{6}\, \text{M}_{\odot}$ and we find the z ∼ 6−7 GSMF to be best parametrized by a modified Schechter function that allows for a turnover at very low masses. Using a Monte Carlo Markov chain analysis of the GSMF, including accurate treatment of lensing uncertainties, we obtain a relatively steep low-mass end slope $\alpha \simeq -1.96_{-0.08}^{+0.09}$ and a turnover at $\log (M_T/\text{M}_{\odot})\simeq 7.10_{-0.56}^{+0.17}$ with a curvature of $\beta \simeq 1.00_{-0.73}^{+0.87}$ for our minimum assumption model with constant star formation history (SFH) and low dust attenuation, AV ≤ 0.2. We find that the z ∼ 6−7 GSMF, in particular its very low mass end, is significantly affected by the assumed functional form of the star formation history and the degeneracy between stellar mass and dust attenuation. For example, the low-mass end slope ranges from $\alpha \simeq -1.82_{-0.07}^{+0.08}$ for an exponentially rising SFH to $\alpha \simeq -2.34_{-0.10}^{+0.11}$ when allowing AV of up to 3.25. Future observations at longer wavelengths and higher angular resolution with the James Webb Space Telescope are required to break these degeneracies and to robustly constrain the stellar mass of galaxies on the extreme low-mass end of the GSMF.


Sign in / Sign up

Export Citation Format

Share Document