Context. Measuring how the physical properties of galaxies change across cosmic times is essential to understand galaxy formation and evolution. With the advent of numerous ground-based and ...space-borne instruments launched over the past few decades we now have exquisite multi-wavelength observations of galaxies from the far-ultraviolet (FUV) to the radio domain. To tap into this mine of data and obtain new insight into the formation and evolution of galaxies, it is essential that we are able to extract information from their spectral energy distribution (SED). Aims. We present a completely new implementation of Code Investigating GALaxy Emission (CIGALE). Written in python, its main aims are to easily and efficiently model the FUV to radio spectrum of galaxies and estimate their physical properties such as star formation rate, attenuation, dust luminosity, stellar mass, and many other physical quantities. Methods. To compute the spectral models, CIGALE builds composite stellar populations from simple stellar populations combined with highly flexible star formation histories, calculates the emission from gas ionised by massive stars, and attenuates both the stars and the ionised gas with a highly flexible attenuation curve. Based on an energy balance principle, the absorbed energy is then re-emitted by the dust in the mid- and far-infrared domains while thermal and non-thermal components are also included, extending the spectrum far into the radio range. A large grid of models is then fitted to the data and the physical properties are estimated through the analysis of the likelihood distribution. Results. CIGALE is a versatile and easy-to-use tool that makes full use of the architecture of multi-core computers, building grids of millions of models and analysing samples of thousands of galaxies, both at high speed. Beyond fitting the SEDs of galaxies and parameter estimations, it can also be used as a model-generation tool or serve as a library to build new applications.
We present a new, publicly available library of dust spectral energy distributions (SEDs). These SEDs are characterized by only three parameters: the dust mass (Mdust), the dust temperature (Tdust), ...and the mid-to-total infrared color (IR8 ≡ LIR/L8). The latter measures the relative contribution of polycyclic aromatic hydrocarbon (PAH) molecules to the total infrared luminosity. We used this library to model star-forming galaxies at 0.5 < z < 4 in the deep CANDELS fields, using both individual detections and stacks of Herschel and ALMA imaging, and extending this sample to z = 0 using the Herschel Reference Survey. At first order, the dust SED of a galaxy was observed to be independent of stellar mass, but evolving with redshift. We found trends of increasing Tdust and IR8 with redshift and distance from the SFR–M∗ main sequence, and quantified for the first time their intrinsic scatter. Half of the observed variations of these parameters was captured by the above empirical relations, and after subtracting the measurement errors we found residual scatters of ΔTdust/Tdust = 12% and Δlog IR8 = 0.18 dex. We observed second order variations with stellar mass: massive galaxies (M∗ > 1011M⊙) at z ≤ 1 have slightly lower temperatures indicative of a reduced star formation efficiency, while low mass galaxies (M∗ < 1010M⊙) at z ≥ 1 showed reduced PAH emission, possibly linked to their lower metallicities. Building on these results, we constructed high-fidelity mock galaxy catalogs to predict the accuracy of infrared luminosities and dust masses determined using a single broadband measurement. Using a single James Webb Space Telescope (JWST) MIRI band, we found that LIR is typically uncertain by 0.15 dex, with a maximum of 0.25 dex when probing the rest-frame 8 μm, and this is not significantly impacted by typical redshift uncertainties. On the other hand, we found that ALMA bands 8 to 7 and 6 to 3 measured the dust mass at better than 0.2 and 0.15 dex, respectively, and independently of redshift, while bands 9 to 6 only measured LIR at better than 0.2 dex at z > 1, 3.2, 3.8, and 5.7, respectively. Starburst galaxies had their LIR significantly underestimated when measured by a single JWST or ALMA band, while their dust mass from a single ALMA band were moderately overestimated. This dust library and the results of this paper can be used immediately to improve the design of observing proposals, and interpret more accurately the large amount of archival data from Spitzer, Herschel and ALMA.
ABSTRACT
Code Investigating GALaxy Emission (cigale) is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of cigale is ...not able to fit X-ray data, which often provide unique insights into active galactic nucleus (AGN) intrinsic power. We develop a new X-ray module for cigale, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of cigale from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, skirtor. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named ‘x-cigale’). We test x-cigale with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced χ2) is good in general, indicating that x-cigale is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using x-cigale, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies.
The star-formation history (SFH) of galaxies is a key assumption to derive their physical properties and can lead to strong biases. In this work, we derive the SFH of main sequence (MS) galaxies and ...show how the peak SFH of a galaxy depends on its seed mass at, for example, z = 5. This seed mass reflects the galaxy’s underlying dark matter (DM) halo environment. We show that, following the MS, galaxies undergo a drastic slow down of their stellar mass growth after reaching the peak of their SFH. According to abundance matching, these masses correspond to hot and massive DM halos which state could result in less efficient gas inflows on the galaxies and thus could be the origin of limited stellar mass growth. As a result, we show that galaxies, still on the MS, can enter the passive region of the UVJ diagram while still forming stars. The best fit to the MS SFH is provided by a right skew peak function for which we provide parameters depending on the seed mass of the galaxy. The ability of the classical analytical SFHs to retrieve the star-formation rate (SFR) of galaxies from spectral energy distribution (SED) fitting is studied. Due to mathematical limitations, the exponentially declining and delayed SFH struggle to model high SFR, which starts to be problematic at z > 2. The exponentially rising and log-normal SFHs exhibit the opposite behavior with the ability to reach very high SFR, and thus model starburst galaxies, but they are not able to model low values such as those expected at low redshift for massive galaxies. By simulating galaxies SED from the MS SFH, we show that these four analytical forms recover the SFR of MS galaxies with an error dependent on the model and the redshift. They are, however, sensitive enough to probe small variations of SFR within the MS, with an error ranging from 5 to 40% depending on the SFH assumption and redshift; but all the four fail to recover the SFR of rapidly quenched galaxies. However, these SFHs lead to an artificial gradient of age, parallel to the MS, which is not exhibited by the simulated sample. This gradient is also produced on real data as we show using a sample of real galaxies with redshifts between 1.5 and 2.5. Here, we propose an SFH composed of a delayed form to model the bulk of stellar population with the addition of a flexibility in the recent SFH. This SFH provides very good estimates of the SFR of MS, starbursts, and rapidly quenched galaxies at all redshift. Furthermore, when used on the real sample, the age gradient disappears which show its dependency on the SFH assumption made to perform the SED fitting.
Aims. We use high-resolution continuum images obtained with the Atacama Large Millimeter Array (ALMA) to probe the surface density of star formation in z ~ 2 galaxies and study the different physical ...properties between galaxies within and above the star-formation main sequence of galaxies. Methods. We use ALMA images at 870 μm with 0.2 arcsec resolution in order to resolve star formation in a sample of eight star-forming galaxies at z ~ 2 selected among the most massive Herschel galaxies in the GOODS-South field. This sample is supplemented with eleven galaxies from the public data of the 1.3 mm survey of the Hubble Ultra-Deep Field, HUDF. We derive dust and gas masses for the galaxies, compute their depletion times and gas fractions, and study the relative distributions of rest-frame ultraviolet (UV) and far-infrared (FIR) light. Results. ALMA reveals systematically dense concentrations of dusty star formation close to the center of the stellar component of the galaxies. We identify two different starburst regimes: (i) the classical population of starbursts located above the SFR-M⋆ main sequence, with enhanced gas fractions and short depletion times and (ii) a sub-population of galaxies located within the scatter of the main sequence that experience compact star formation with depletion timescales typical of starbursts of ~150 Myr. In both starburst populations, the FIR and UV are distributed in distinct regions and dust-corrected star formation rates (SFRs) estimated using UV-optical-near-infrared data alone underestimate the total SFR. Starbursts hidden in the main sequence show instead the lowest gas fractions of our sample and could represent the last stage of star formation prior to passivization. Being Herschel-selected, these main sequence galaxies are located in the high-mass end of the main sequence, hence we do not know whether these “starbursts hidden in the main sequence” also exist below 1011 M⊙. Active galactic nuclei (AGNs) are found to be ubiquitous in these compact starbursts, suggesting that the triggering mechanism also feeds the central black hole or that the active nucleus triggers star formation.
We study the star formation quenching mechanism in cluster galaxies by fitting the spectral energy distribution (SED) of the Herschel Reference Survey, a complete volume-limited K-band-selected ...sample of nearby galaxies including objects in different density regions, from the core of the Virgo cluster to the general field. The SEDs of the target galaxies were fitted using the CIGALE SED modelling code. The truncated activity of cluster galaxies was parametrised using a specific star formation history with two free parameters, the quenching age QA and the quenching factor QF. These two parameters are crucial for the identification of the quenching mechanism, which acts on long timescales when starvation processes are at work, but is rapid and efficient when ram pressure occurs. To be sensitive to an abrupt and recent variation of the star formation activity, we combined twenty photometric bands in the UV to far-infrared in a new way with three age-sensitive Balmer line absorption indices extracted from available medium-resolution (R ~ 1000) integrated spectroscopy and with Hα narrow-band imaging data. The use of a truncated star formation history significantly increases the quality of the fit in HI-deficient galaxies of the sample, that is to say, in those objects whose atomic gas content has been removed during the interaction with the hostile cluster environment. The typical quenching age of the perturbed late-type galaxies is QA ≲ 300 Myr whenever the activity of star formation is reduced by 50% < QF ≤ 80% and QA ≲ 500 Myr for QF > 80%, while that of the quiescent early-type objects is QA ≃ 1−3 Gyr. The fraction of late-type galaxies with a star formation activity reduced by QF > 80% and with an HI-deficiency parameter HI−def > 0.4 drops by a factor of ~5 from the inner half virial radius of the Virgo cluster (R/Rvir < 0.5), where the hot diffuse X-ray emitting gas of the cluster is located, to the outer regions (R/Rvir > 4). The efficient quenching of the star formation activity observed in Virgo suggests that the dominant stripping process is ram pressure. We discuss the implication of this result in the cosmological context of galaxy evolution.
Context. Over the past few years the Atacama Large Millimeter Array (ALMA) has detected dust-rich galaxies whose cold dust emission is spatially disconnected from the ultraviolet (UV) rest-frame ...emission. This represents a challenge for modeling their spectral energy distributions (SED) with codes based on an energy budget between the stellar and dust components. This could potentially weaken the solidity of the physical parameters measured with these modeling tools. Aims. We want to verify the validity of energy balance modeling on a sample of galaxies observed from the UV to the sub-millimeter rest frame with ALMA and decipher what information can be reliably retrieved from the analysis of the full SED and from subsets of wavelengths. Methods. We select 17 sources at z ≃ 2 in the Hubble Ultra-Deep Field (HUDF) and in the GOODS-South field detected with ALMA and Herschel and for which UV to near-infrared rest-frame ancillary data are available. We fit the data with CIGALE exploring different configurations for dust attenuation and star formation histories, considering either the full dataset or one that is reduced to the stellar and dust emission. We compare estimates of the dust luminosities, star formation rates, and stellar masses. Results. The fit of the stellar continuum alone with the starburst attenuation law can only reproduce up to 50% of the total dust luminosity observed by Herschel and ALMA. This deficit is found to be marginally consistent with similar quantities estimated in the COSMOS field and is found to increase with the specific star formation rate. The combined stellar and dust SEDs are well fitted when different attenuation laws are introduced. Shallow attenuation curves are needed for the galaxies whose cold dust distribution is very compact compared to starlight. The stellar mass estimates are affected by the choice of the attenuation law. The star formation rates are robustly estimated as long as dust luminosities are available. The large majority of the galaxies are above the average main sequence of star forming galaxies and one source is a strong starburst.
Aims.
Dust is observed in the polar regions of nearby active galactic nuclei (AGN) and it is known to contribute substantially to their mid-IR emission and to the obscuration of their UV to optical ...emission. We aim to carry out a statistical test to check whether this component is a common feature based on an analysis of the integrated spectral energy distributions of these composite sources.
Methods.
We selected a sample of 1275 broad-line AGN in the XMM-XXL field, with optical to infrared photometric data. These AGN are seen along their polar direction and we expect a maximal impact of dust located around the poles when it is present. We used X-CIGALE, which introduces a dust component to account for obscuration along the polar directions, modeled as a foreground screen, and an extinction curve that is chosen as it steepens significantly at short wavelengths or is much grayer. By comparing the results of different fits, we are able to define subsamples of sources with positive statistical evidence in favor of or against polar obscuration (if present) and described using the gray or steep extinction curve.
Results.
We find a similar fraction of sources with positive evidence for and against polar dust. Applying statistical corrections, we estimate that half of our sample could contain polar dust and among them, 60% exhibit a steep extinction curve and 40% a flat extinction curve; although these latter percentages are found to depend on the adopted extinction curves. The obscuration in the
V
-band is not found to correlate with the X-ray column density, while
A
V
/
N
H
ratios span a large range of values and higher dust temperatures are found with the flat, rather than with the steep extinction curve. Ignoring this polar dust component in the fit of the spectral energy distribution of these composite systems leads to an overestimation of the stellar contribution. A single fit with a polar dust component described with an SMC extinction curve efficiently overcomes this issue but it fails at identifying all the AGN with polar dust obscuration.
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.
A semi-empirical model is presented that describes the distribution of active galactic nuclei (AGNs) on the cosmic web. It populates dark-matter haloes in N-body simulations (MultiDark) with galaxy ...stellar masses using empirical relations based on abundance matching techniques, and then paints accretion events on these galaxies using state-of-the-art measurements of the AGN occupation of galaxies. The explicit assumption is that the large-scale distribution of AGN is independent of the physics of black hole fuelling. The model is shown to be consistent with current measurements of the two-point correlation function of AGN samples. It is then used to make inferences on the halo occupation of the AGN population. Mock AGNs are found in haloes with a broad distribution of masses with a mode of |${\approx } 10^{12}\, h^{-1} \, M_{\odot }$| and a tail extending to cluster-size haloes. The clustering properties of the model AGN depend only weakly on accretion luminosity and redshift. The fraction of satellite AGN in the model increases steeply toward more massive haloes, in contrast with some recent observational results. This discrepancy, if confirmed, could point to a dependence of the halo occupation of AGN on the physics of black hole fuelling.