Source author record

Denis Burgarella

Denis Burgarella appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

A hyperluminous obscured quasar at a redshift of z ~ 4.3

In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of $3.91^{+1.69}_{-0.55} \times 10^{13} L_\odot$ and an active galactic nucleus (AGN) fraction of ~89%. We also fit HELP_J100156.75+022344.7 with the Code Investigating GALaxy Emission (CIGALE) code and find a similar result. This is only the second z > 4 hyperluminous obscured quasar discovered to date. The discovery of HELP_J100156.75+022344.7 in the ~ 2deg^2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields which cover ~ 1300deg^2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP_J100156.75+022344.7 at z ~ 4.5 to be $\sim 1.8 \times 10^{-8}$Mpc$^{-3}$. This is slightly higher than the space density of coeval hyperluminous optically selected quasars suggesting that the obscuring torus in z > 4 quasars may have a covering factor $\gtrsim 50\%$.

preprint2022arXiv

Fitting AGN/galaxy X-ray-to-radio SEDs with CIGALE and improvement of the code

Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(θ) \propto 1+\cos θ$, where $θ$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr

preprint2022arXiv

Tempered, Anti-trunctated, Multiple Importance Sampling

Importance sampling is a Monte Carlo method that introduces a proposal distribution to sample the space according to the target distribution. Yet calibration of the proposal distribution is essential to achieving efficiency, thus the resort to adaptive algorithms to tune this distribution. In the paper, we propose a new adpative importance sampling scheme, named Tempered Anti-truncated Adaptive Multiple Importance Sampling (TAMIS) algorithm. We combine a tempering scheme and a new nonlinear transformation of the weights we named anti-truncation. For efficiency, we were also concerned not to increase the number of evaluations of the target density. As a result, our proposal is an automatically tuned sequential algorithm that is robust to poor initial proposals, does not require gradient computations and scales well with the dimension.

preprint2022arXiv

The ALPINE-ALMA [CII] survey: dust attenuation curves at z=4.4-5.5

There is now ample evidence that dust is already present in abundance at high z. However, given the faintness of distant galaxies in the optical and the NIR, datasets are still limited and how the dust affects the emerging radiation of galaxies at very high redshift is not yet fully understood. Using the ALPINE survey, our objective is to quantify the dust attenuation properties in galaxies at z=4.4-5.5, and in particular the shape of their attenuation curve. Using the CIGALE code, we model the stellar populations and their interaction with the dust in order to measure some of the physical properties of a subsample of 23 main-sequence ALPINE galaxies. We find that the attenuation curves span a broad range of properties, from curves that are much steeper than the SMC extinction curve, to shallower than the starburst attenuation curve. The shape of the attenuation curves strongly depends on the V-band attenuation. Galaxies with the lowest attenuation also present the steepest curves. The steepness of such curves is probably the consequence of the combination of the intrinsic physical properties of the dust, the relative distribution of stars and dust in the interstellar medium, and the differential reddening. The broad range of attenuation curves found at z~5 shows that no single attenuation curve is appropriate for main sequence galaxies and that assuming a fixed curve can lead to large errors, for instance in the interpretation and use of the IRX-beta diagram, if SED modeling is not feasible. Great caution should be exercised when correcting high redshift galaxies for the presence of dust using the UV slope beta as it can affect the estimation of both SFR and stellar mass even at low V-band attenuation due to the steepness of the attenuation curve. However, when SED modeling can be used, the impact of the choice of the attenuation curve on the SFR and the stellar mass is limited.

preprint2020arXiv

JWST/MIRI Simulated Imaging: Insights into Obscured Star-Formation and AGN for Distant Galaxies in Deep Surveys

The JWST MIRI instrument will revolutionize extragalactic astronomy with unprecedented sensitivity and angular resolution in mid-IR. Here, we assess the potential of MIRI photometry to constrain galaxy properties in the Cosmic Evolution Early Release Science (CEERS) survey. We derive estimated MIRI fluxes from the spectral energy distributions (SEDs) of real sources that fall in a planned MIRI pointing. We also obtain MIRI fluxes for hypothetical AGN-galaxy mixed models varying the AGN fractional contribution to the total IR luminosity ($\rm frac_{AGN}$). Based on these model fluxes, we simulate CEERS imaging (3.6-hour exposure) in 6 bands from F770W to F2100W using MIRISIM, and reduce these data using JWST PIPELINE. We perform PSF-matched photometry with TPHOT, and fit the source SEDs with X-CIGALE, simultaneously modeling photometric redshift and other physical properties. Adding the MIRI data, the accuracy of both redshift and $\rm frac_{AGN}$ is generally improved by factors of $\gtrsim 2$ for all sources at $z\lesssim 3$. Notably, for pure-galaxy inputs ($\rm frac_{AGN}=0$), the accuracy of $\rm frac_{AGN}$ is improved by $\sim 100$ times thanks to MIRI. The simulated CEERS MIRI data are slightly more sensitive to AGN detections than the deepest X-ray survey, based on the empirical $L_{\rm X}$-$L_{\rm 6μm}$ relation. Like X-ray observations, MIRI can also be used to constrain the AGN accretion power (accuracy $\approx 0.3$ dex). Our work demonstrates that MIRI will be able to place strong constraints on the mid-IR luminosities from star formation and AGN, and thereby facilitate studies of the galaxy/AGN co-evolution.

preprint2020arXiv

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo.

preprint2020arXiv

Search for Optically Dark Infrared Galaxies without Counterparts of Subaru Hyper Suprime-Cam in the AKARI North Ecliptic Pole Wide Survey Field

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $μ$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $μ$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

preprint2020arXiv

UV dust attenuation as a function of stellar mass and its evolution with redshift

Studying the UV dust attenuation, as well as its relation to other galaxy parameters such as the stellar mass, plays an important role in multi-wavelength research. This work relates the dust attenuation to the stellar mass of star forming galaxies, and its evolution with redshift. A sample of galaxies with an estimate of the dust attenuation computed from the infrared excess was used. The dust attenuation vs. stellar mass data, separated in redshift bins, was modelled by a single parameter linear function, assuming a nonzero constant apparent dust attenuation for low mass galaxies. But the origin of this effect is still to be determined and several possibilities are explored (actual high dust content, variation of the dust-to-metal ratio, variation of the stars-dust geometry). The best-fitting parameter of this model is then used to study the redshift evolution of the cosmic dust attenuation and is found to be in agreement with results from the literature. This work also gives evidence to a redshift evolution of the dust attenuation - stellar mass relationship, as is suggested by recent works in the highest redshift range.

preprint2019arXiv

HELP project - a dreamed-of multiwavelength dataset for SED fitting: the influence of used models for the main physical properties of galaxies

The Herschel Extragalactic Legacy Project (HELP) focuses to publish an astronomical multiwavelength catalogue of millions of objects over 1300~deg$^2$ of the Herschel Space Observatory survey fields. Millions of galaxies with ultraviolet--far infrared photometry {make} HELP a perfect sample for testing spectral energy distribution fitting models, and to prepare tools for next-generation data. In the frame of HELP collaboration we estimated the main physical properties of all galaxies from the HELP database and we checked a new procedure to select peculiar galaxies from large galaxy sample and we investigated the influence of used modules for stellar mass estimation.

preprint2019arXiv

Properties of LBGs with [OIII] detection at z $\sim$ 3.5: The importance of including nebular emission data in SED fitting

Nebular emission lines are critical to measure physical properties in the ionized gas (e.g., metallicity, the SFR, or dust attenuation). They also account for a significant fraction of broadband fluxes, in particular at the highest redshifts, and therefore can strongly affect the determination of other physical properties, such as the stellar mass, which are crucial in shaping our understanding of galaxy formation and evolution. We investigate a sample of 51 LBGs at 3.0<z<3.8 with detected [OIII] line emissions and estimated the physical properties of these galaxies to examine the impact of including nebular emission data in the SED fitting. We used CIGALE to fit the rest-frame UV-to-NIR SEDs of these galaxies and their emission line data simultaneously. We ran CIGALE with and without the nebular model or the emission line data, and compare the results to show the importance of including the nebular emission line data in the SED fitting. We find that without the nebular model, the fitting overestimates the stellar mass due to the strong [OIII] lines that are redshifted to the Ks-band, consistent with previous results. The emission line data are necessary to constrain the nebular model in the SED fitting. Examining the Ks-band excess, which is mostly used to estimate the emissions of the [OIII]+Hb lines when there is no spectral data, we find that the estimation and observation are statistically consistent. However, the difference can reach up to more than 1 dex in some catastrophic cases, which shows the importance of obtaining spectroscopic measurements for these lines. We also estimate the equivalent width of the Hb absorption and find it negligible compared to the Hb emission. We conclude that line emission is important to constrain the nebular models and to obtain reliable estimates of the physical properties of galaxies. These data should be taken into account in the SED fitting.

preprint2016arXiv

Cosmic Dawn Intensity Mapper

Cosmic Dawn Intensity Mapper is a "Probe Class" mission concept for reionization studies of the universe. It will be capable of spectroscopic imaging observations between 0.7 to 6-7 microns in the near-Infrared. The primary observational objective is pioneering observations of spectral emission lines of interest throughout the cosmic history, but especially from the first generation of distant, faint galaxies when the universe was less than 800 million years old. With spectro-imaging capabilities, using a set of linear variable filters (LVFs), CDIM will produce a three-dimensional tomographic view of the epoch of reionization (EoR). CDIM will also study galaxy formation over more than 90% of the cosmic history and will move the astronomical community from broad-band astronomical imaging to low-resolution (R=200-300) spectro-imaging of the universe.

preprint2016arXiv

GALEX-SDSS-WISE Legacy Catalog (GSWLC): Star Formation Rates, Stellar Masses and Dust Attenuations of 700,000 Low-redshift Galaxies

In this paper, we present GALEX-SDSS-WISE Legacy Catalog (GSWLC), a catalog of physical properties (stellar masses, dust attenuations and star formation rates (SFRs)) of ~700,000 galaxies with SDSS redshifts below 0.3. GSWLC contains galaxies within the GALEX footprint, regardless of a UV detection, covering 90% of SDSS. The physical properties were obtained from UV/optical SED fitting following Bayesian methodology of Salim et al. (2007), with improvements such as blending corrections for low-resolution UV photometry, flexible dust attenuation laws, and emission line corrections. GSWLC includes mid-IR SFRs derived from IR templates based upon 22 micron WISE observations. These estimates are independent of UV/optical SED fitting, in order to separate possible systematics. The paper argues that the comparison of specific SFRs (SSFRs) is more informative and physically motivated than the comparison of SFRs. SSFRs resulting from the UV/optical SED fitting are compared to the mid-IR SSFRs, and to SSFRs from three published catalogs. For "main sequence" galaxies with no AGN contribution all SSFRs are in very good agreement (within 0.1 dex on average). In particular, the widely-used aperture-corrected SFRs from MPA/JHU catalog show no systematic offsets, in contrast to some integral-field spectroscopy results. For galaxies below the main sequence (log SSFR<-11), mid-IR (S)SFRs based on fixed luminosity-SFR conversion are severely biased (up to 2 dex) because the dust is primarily heated by old stars. Furthermore, mid-IR (S)SFRs are overestimated by up to 0.6 dex for galaxies with AGN, presumably due to non-stellar dust heating. UV/optical (S)SFRs are thus preferred to IR-based (S)SFRs for quenched galaxies and those which host AGN.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2015arXiv

Evolution of mid-infrared galaxy luminosity functions from the entire AKARI NEP-Deep field with new CFHT photometry

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole (NEP) deep field using recently-obtained, wider CFHT optical/near-IR images. AKARI has obtained deep images in the mid-infrared (IR), covering 0.6 deg$^2$ of the NEP deep field. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, we recently obtained CFHT optical and near-IR images over the entire AKARI NEP deep field. These new CFHT images are used to derive accurate photometric redshifts, allowing us to fully exploit the whole AKARI NEP deep field. AKARI's deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$μ$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$μ$m and 12$μ$m luminosities without using a large extrapolation based on spectral energy distribution (SED) fitting, which was the largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained more reliably due to the superior filter coverage. The resulting restframe 8$μ$m, 12$μ$m, and TIR LFs at $0.15<z<2.2$ are consistent with previous works, but with reduced uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. In terms of cosmic infrared luminosity density ($Ω_{\mathrm{IR}}$), we found that the $Ω_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

preprint2013arXiv

A new HST/Herschel deep field at the North Ecliptic Pole: preparing the way for JWST, SPICA and Euclid

We propose a co-ordinated multi-observatory survey at the North Ecliptic Pole. This field is the natural extragalactic deep field location for most space observatories (e.g. containing the deepest Planck, WISE and eROSITA data), is in the continuous viewing zones for e.g. Herschel, HST, JWST, and is a natural high-visibility field for the L2 halo orbit of SPICA with deep and wide-field legacy surveys already planned. The field is also a likely deep survey location for the forthcoming Euclid mission. It is already a multi-wavelength legacy field in its own right (e.g. AKARI, LOFAR, SCUBA-2): the outstanding and unparalleled continuous mid-IR photometric coverage in this field and nowhere else enables a wide range of galaxy evolution diagnostics unachievable in any other survey field, by spanning the wavelengths of redshifted PAH and silicate features and the peak energy output of AGN hot dust. We argue from the science needs of Euclid and JWST, and from the comparative multiwavelength depths, that the logical approach is (1) a deep (H-UDF) UV/optical tile in the NEP over ~10 square arcminutes, and (2) an overlapping wide-field UV/optical HST survey tier covering >100 square arcminutes, with co-ordinated submm SPIRE mapping up to or beyond the submm point source confusion limit over a wider area and PACS data over the shallower HST tier.

preprint2013arXiv

Herschel observations: contraints on dust attenuation and star formation histories at high redshift

SPICA is one of the key projects for the future. Not only its instrument suite will open up a discovery window but they will also allow to physically understand some of the phenomena that we still do not understand in the high-redshift universe. Using new homogeneous luminosity functions (LFs) in the Far-Ultraviolet (FUV) from VVDS and in the Far-Infrared (FIR) from Herschel/PEP and Herschel/HerMES, we studied the evolution of the dust attenuation with redshift. With this information, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density (SFRD_TOT). Our main conclusions are that: 1) the dust attenuation A_FUV is found to increase from z = 0 to $z \sim 1.2 and then starts to decrease until our last data point at z = 3.6; 2) the estimated SFRD confirms published results to z \sim 2. At z > 2, we observe either a plateau or a small increase up to z \sim 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT but the origin of this delay is not understood yet, and SPICA instruments will provide clues to move further in the physical understanding of this delay but also on the detection and redshift measurements of galaxies at higher redshifts. This work is further detailed in Burgarella et al. (2013).

preprint2013arXiv

pcigale: porting Code Investigating Galaxy Emission to Python

We present pcigale, the port to Python of CIGALE (Code Investigating Galaxy Emission) a Fortran spectral energy distribution (SED) fitting code developed at the Laboratoire d'Astrophysique de Marseille. After recalling the specifics of the SED fitting method, we show the gains in modularity and versatility offered by Python, as well as the drawbacks compared to the compiled code.

preprint2012arXiv

Far-Ultraviolet and Far-Infrared Bivariate Luminosity Function of Galaxies: Complex Relation between Stellar and Dust Emission

Far-ultraviolet (FUV) and far-infrared (FIR) luminosity functions (LFs) of galaxies show a strong evolution from $z = 0$ to $z = 1$, but the FIR LF evolves much stronger than the FUV one. The FUV is dominantly radiated from newly formed short-lived OB stars, while the FIR is emitted by dust grains heated by the FUV radiation field. It is known that dust is always associated with star formation activity. Thus, both FUV and FIR are tightly related to the star formation in galaxies, but in a very complicated manner. In order to disentangle the relation between FUV and FIR emissions, we estimate the UV-IR bivariate LF (BLF) of galaxies with {\sl GALEX} and {\sl AKARI} All-Sky Survey datasets. Recently we invented a new mathematical method to construct the BLF with given marginals and prescribed correlation coefficient. This method makes use of a tool from mathematical statistics, so called "copula". The copula enables us to construct a bivariate distribution function from given marginal distributions with prescribed correlation and/or dependence structure. With this new formulation and FUV and FIR univariate LFs, we analyze various FUV and FIR data with {\sl GALEX}, {\sl Spitzer}, and {\sl AKARI} to estimate the UV-IR BLF. The obtained BLFs naturally explain the nonlinear complicated relation between FUV and FIR emission from star-forming galaxies. Though the faint-end of the BLF was not well constrained for high-$z$ samples, the estimated linear correlation coefficient $ρ$ was found to be very high, and is remarkably stable with redshifts (from 0.95 at $z = 0$ to 0.85 at $z = 1.0$). This implies the evolution of the UV-IR BLF is mainly due to the different evolution of the univariate LFs, and may not be controlled by the dependence structure.

preprint2011arXiv

CIGALEMC: Galaxy Parameter Estimation using a Markov Chain Monte Carlo Approach with Cigale

We introduce a fast Markov Chain Monte Carlo (MCMC) exploration of the astrophysical parameter space using a modified version of the publicly available code CIGALE (Code Investigating GALaxy emission). The original CIGALE builds a grid of theoretical Spectral Energy Distribution (SED) models and fits to photometric fluxes from Ultraviolet (UV) to Infrared (IR) to put contraints on parameters related to both formation and evolution of galaxies. Such a grid-based method can lead to a long and challenging parameter extraction since the computation time increases exponentially with the number of parameters considered and results can be dependent on the density of sampling points, which must be chosen in advance for each parameter. Markov Chain Monte Carlo methods, on the other hand, scale approximately linearly with the number of parameters, allowing a faster and more accurate exploration of the parameter space by using a smaller number of efficiently chosen samples. We test our MCMC version of the code CIGALE (called CIGALEMC) with simulated data. After checking the ability of the code to retrieve the input parameters used to build the mock sample, we fit theoretical SEDs to real data from the well known and studied SINGS sample. We discuss constraints on the parameters and show the advantages of our MCMC sampling method in terms of accuracy of the results and optimization of CPU time.

preprint2011arXiv

The Recent Star Formation in NGC 6822: an Ultraviolet Study

We characterize the star formation in the low-metallicity galaxy NGC 6822 over the past few hundred million years, using GALEX far-UV (FUV, 1344-1786 A) and near-UV (NUV, 1771-2831 A) imaging, and ground-based Ha imaging. From GALEX FUV image, we define 77 star-forming (SF) regions with area >860 pc^2, and surface brightness <=26.8 mag(AB)arcsec^-2, within 0.2deg (1.7kpc) of the center of the galaxy. We estimate the extinction by interstellar dust in each SF region from resolved photometry of the hot stars it contains: E(B-V) ranges from the minimum foreground value of 0.22mag up to 0.66+-0.21mag. The integrated FUV and NUV photometry, compared with stellar population models, yields ages of the SF complexes up to a few hundred Myr, and masses from 2x10^2 Msun to 1.5x10^6 Msun. The derived ages and masses strongly depend on the assumed type of interstellar selective extinction, which we find to vary across the galaxy. The total mass of the FUV-defined SF regions translates into an average star formation rate (SFR) of 1.4x10^-2 Msun/yr over the past 100 Myr, and SFR=1.0x10^-2 Msun/yr in the most recent 10 Myr. The latter is in agreement with the value that we derive from the Ha luminosity, SFR=0.008 Msun/yr. The SFR in the most recent epoch becomes higher if we add the SFR=0.02 Msun/yr inferred from far-IR measurements, which trace star formation still embedded in dust (age <= a few Myr).

preprint2010arXiv

Population synthesis modelling of luminous infrared galaxies at intermediate redshift

Luminous InfraRed Galaxies (LIRGs) are particularly important for studying the build-up of the stellar mass from z=1 to z=0. We perform a multiwavelength study of an LIRGs sample in the Extended Chandra Deep Field South at z=0.7, selected at 24 μ\m by MIPS onboard Spitzer Space Telescope and detected in 17 filters. Data go from the near-ultraviolet to the mid-infrared. This multiwavelengths dataset allows us to place strong constraints on the spectral energy distributions (SEDs) of galaxies, and thus to efficiently derive physical parameters such as the SFR, the total infrared luminosity, attenuation parameters, and star formation history. An important part of this work is elaboration of a mock catalogue that allows us to have a reliability criterion for the derived parameters. We studied LIRGs by means of an SED-fitting code CIGALE. At first, this code creates synthetic spectra from the Maraston stellar population models. The stellar population spectra are attenuated by using a synthetic Calzetti-based attenuation law before adding the dust emission as given by the infrared SED library. The originality of CIGALE is that it allows us to perform consistent fits of the dust-affected ultaviolet-to-infrared wavelength range. This technique appears to be a very powerful tool in the case where we can have access to a dataset that is well-sampled over a wide range of wavelengths. We are able to derive a star formation history and to estimate the fraction of infrared luminosity reprocessed by an active galactic nucleus. We study the dust temperatures of our galaxies detected at 70 μ\m and find them colder than predicted by models. We also study the relation between the SFR and the stellar mass and do not find a tight correlation between either of them, but instead a flat distribution and a large scatter, which is interpreted in terms of variations in star formation history.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2008arXiv

Deep 15um AKARI observations in the CDFS: estimating dust luminosities for a MIR-selected sample and for Lyman Break Galaxies and the evolution of L(dust)/L(UV) with the redshift

Deep observations of the CDFS have been secured at 15um with AKARI/IRC infrared space telescope (ESA open time). From these observations, we define a sample of MIR-selected galaxies at 15um and we also obtain 15um flux densities for a sample of LBGs at z=1 already observed at 24um with Spitzer/MIPS. Number counts for the MIR-selected sample show a bump around a 15um flux density of 0.2mJy that can be attributed to galaxies at z>0.4 and at z>0.8 for the fainter part of the bump. This bump seems to be shifted as compared to other works and a possible origin can be the Cosmic variance. Thanks to this dataset, we have tested, on the two above samples at z=1, the validity of the conversions from monochromatic luminosities nu.f(nu) at a rest-frame wavelength of 8um by a comparison with total dust luminosities estimated from Spitzer rest-frame 12um data that we use as a reference. We find that the 8um dust luminosities are not all consistent and that some of them are better when compared to L(dust) evaluated from longer wavelength luminosities. We also find that the rest-frame 8um luminosities provide globally good estimates of L(dust). By comparing our data for the two samples to several libraries of SEDs, we find that models can explain the diversity of the observed f(24)/f(15) ratio quite reasonably for the MIR-selected sample and better for the LBG sample which are less dispersed than the MIR selection. However, when we analyse the luminosity dependence of this ratio, we find important discrepancies. Finally, we revisit the evolution of L(dust)/L(UV) ratio with the redshift z by re-calibrating previous L(dust) at z=2 based on our results and added new data points at higher redshifts. The decreasing trend is amplified as compared to the previous estimate.

preprint2005arXiv

Evolution of morphology in the Chandra Deep Field South

We studied the morphology of galaxies in the Chandra Deep Field South using ACS multi-wavelength data from the Great Observatories Origin Deep Survey and 524 spectroscopic redshifts from the VIMOS VLT Deep Survey completed with 2874 photometric redshift computed from COMBO-17 multi-color data. The rest-frame B-band makes it possible to discriminate two morphological types in an asymmetry-concentration diagram: bulge- and disk-dominated galaxies. The rest-frame color index $B-I$ is found to be very correlated with the morphological classification: wholly bulge-dominated galaxies are redder than disk-dominated galaxies. However color allowed us to distinguish a population of faint blue bulge-dominated galaxies ($B-I<0.9$), whose nature is still unclear. Using the rest-frame B-band classification from $z\sim0.15$ up to $z\sim1.1$, we quantified the evolution of the proportion of morphological types as a function of the redshift. Our large sample allowed us to compute luminosity functions per morphological type in rest-frame B-band. The bulge-dominated population is found to be composite: on the one hand the red ($(B-I)\_{AB}>0.9$), bright galaxies, which seem to increase in density toward low redshifts. On the other hand the blue, compact, faint bulge-dominated galaxies, strongly evolving with the redshift.