Search

Polycyclic aromatic hydrocarbon (PAH) dust emission has been proposed as an effective extinction-independent star formation rate (SFR) indicator in the mid-infrared (MIR), but this may depend on conditions in the interstellar medium. The coverage of the AKARI/Infrared Camera (IRC) allows us to study the effects of metallicity, starburst intensity, and active galactic nuclei on PAH emission in galaxies with $f_{\nu}(L18W)\lesssim 19$ AB mag. Observations include follow-up, rest-frame optical spectra of 443 galaxies within the AKARI North Ecliptic Pole survey that have IRC detections from 7-24 $\mu$m. We use optical emission line diagnostics to infer SFR based on H$\alpha$ and [O II]$\lambda\lambda 3726,3729$ emission line luminosities. The PAH 6.2 $\mu$m and PAH 7.7 $\mu$m luminosities ($L(PAH\ 6.2\ \mu m)$ and $L(PAH\ 7.7\ \mu m)$, respectively) derived using multi-wavelength model fits are consistent with those derived from slitless spectroscopy within 0.2 dex. $L(PAH\ 6.2\ \mu m)$ and $L(PAH\ 7.7\ \mu m)$ correlate linearly with the 24 $\mu$m-dust corrected H$\alpha$ luminosity only for normal, star-forming ``main-sequence" galaxies. Assuming multi-linear correlations, we quantify the additional dependencies on metallicity and starburst intensity, which we use to correct our PAH SFR calibrations at $0

The James Webb Space Telescope (JWST) has discovered a surprising population of bright galaxies in the very early universe (< 500 Myrs after the Big Bang) that is hard to explain with conventional galaxy formation models and whose physical properties remain to be fully understood. Insight into the internal physics of galaxies is captured best via observations of excited-state atomic transitions of ionized gas, but beyond z~7-9, the brightest spectral signatures are redshifted into the mid-infrared regime, where observations are increasingly more difficult. Here, we present the first detection of a hydrogen recombination line (H{\alpha}) and doubly-ionized oxygen ([OIII]4959,5007{\AA}) at z>10 using the JWST Mid-Infrared Instrument, MIRI. These detections place the bright galaxy GHZ2/GLASS-z12 at z=12.33+/-0.02, making it the most distant astronomical object with direct spectroscopic detection of these lines and the brightest confirmed object at this epoch. These observations provide key insights into the conditions of this primeval galaxy, which shows hard ionizing conditions rarely seen in the local Universe and likely driven by compact, young (<30 Myr) star formation. Its oxygen-to-hydrogen abundance is close to a tenth of the solar value, indicating a rapid metal enrichment during the earliest phases of galaxy formation. This study confirms the unique conditions of the brightest and most distant galaxies recently discovered by JWST and the huge potential of mid-IR observations to characterize these systems, opening a range of new possibilities in the study of the very early Universe., Comment: 14 pages, 6 figures. Submitted for publication. Comments are most welcome

It is well known that supermassive black holes (SMBHs) and their host galaxies co-evolve. AGN feedback plays an important role on this symbiosis. To study the effect of the AGN feedback on the host galaxy, a popular method is to study the star-formation rate (SFR) as a function of the X-ray luminosity (L$_X$). However, hydrodynamical simulations suggest that the cumulative impact of AGN feedback on a galaxy is encapsulated in the mass of the SMBH, M$_{BH}$, rather than the L$_X$. In this study, we compare the SFR of AGN and non-AGN galaxies as a function of L$_X$, M$_{BH}$, Eddington ratio (n$_{Edd}$) and specific black hole accretion rate ($\lambda _{sBHAR}$). For that purpose, we use 122 X-ray AGN in the XMM-XXL field and 3371 galaxies from the VIPERS survey and calculate the SFR$_{norm}$ parameter, defined as the ratio of the SFR of AGN to the SFR of non-AGN galaxies with similar stellar mass, M$_*$, and redshift. Our datasets span a redshift range of $\rm 0.5\leq z\leq 1.2$. The results show that the correlation between SFR$_{norm}$ and M$_{BH}$ is stronger compared to that between SFR$_{norm}$ and L$_X$. A weaker correlation is found between SFR$_{norm}$ and $\lambda _{sBHAR}$. No correlation is detected between SFR$_{norm}$ and n$_{Edd}$. These results corroborate the idea that the M$_{BH}$ is a more robust tracer of the cumulative impact of the AGN feedback compared to the instantaneous accretion rate (L$_X$) and, thus, a better predictive parameter of the changes of the SFR of the host galaxy., Comment: A&A accepted, 10 pages, 4 figures, 6 tables

Optical rest-frame spectroscopic diagnostics are usually employed to distinguish between star formation and AGN-powered emission. However, this method is biased against dusty sources, hampering a complete census of the AGN population across cosmic epochs. To mitigate this effect, it is crucial to observe at longer wavelengths in the rest-frame near-infrared (near-IR), which is less affected by dust attenuation and can thus provide a better description of the intrinsic properties of galaxies. AGN diagnostics in this regime have not been fully exploited so far, due to the scarcity of near-IR observations of both AGNs and star-forming galaxies, especially at redshifts higher than 0.5. Using Cloudy photoionization models, we identify new AGN - star formation diagnostics based on the ratio of bright near-infrared emission lines, namely [SIII] 9530 Angstrom, [CI] 9850 Angstrom, [PII] 1.188 $\mu m$, [FeII] $1.257 \mu m$, and [FeII] $1.64 \mu m$ to Paschen lines (either Pa$\gamma$ or Pa$\beta$), providing simple, analytical classification criteria. We apply these diagnostics to a sample of 64 star-forming galaxies and AGNs at 0 < z < 1, and 65 sources at 1 < z < 3 recently observed with JWST-NIRSpec in CEERS. We find that the classification inferred from the near-infrared is broadly consistent with the optical one based on the BPT and the [SII]/H$\alpha$ ratio. However, in the near-infrared, we find $\sim 60 \%$ more AGNs than in the optical (13 instead of 8), with 5 sources classified as 'hidden' AGNs, showing a larger AGN contribution at longer wavelengths, possibly due to the presence of optically thick dust. The diagnostics we present provide a promising tool to find and characterize AGNs from z=0 to z=3 with low and medium-resolution near-IR spectrographs in future surveys., Comment: Paper accepted for publication in A&A on 05/09/2023. Three public Github repositories include: (1) a table with emission line measurements for the paper sample : https://github.com/Anthony96/Line_measurements_nearIR , Cloudy emission line predictions for star-forming galaxies and AGN models : https://github.com/Anthony96/star-forming_models , https://github.com/Anthony96/AGN_models

[Abridged] We combined HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 69 star-forming galaxies (SFGs) at $0.175$; 35 galaxies). At lower $S/N$, simulations show that simultaneous measurements of the size and S\'ersic index become unreliable. We extended our study to fainter sources ($S/N>10$; 69 galaxies) by fixing their S\'ersic index to unity. The S\'ersic index of bright galaxies ($S/N>75$) has a median value of 0.7, which, together with their axis ratio distribution, suggests a disk-like morphology in the rest-MIR. Galaxies above the main sequence (MS; i.e., starbursts) have rest-MIR sizes that are a factor 2 smaller than their rest-optical sizes. The median rest-optical to rest-MIR size ratio of MS galaxies increases with stellar mass, from 1.1 at $10^{9.8}M_\odot$ to 1.6 at $10^{11}M_\odot$. This mass-dependent trend resembles the one found in the literature between the rest-optical and rest-near-infrared sizes of SFGs, suggesting that it is due to radial color gradients affecting rest-optical sizes and that the sizes of the stellar and star-forming components of SFGs are, on average, consistent at all masses. There is, however, a small population of SFGs (15%) with a compact star-forming component embedded in a larger stellar structure. This could be the missing link between galaxies with an extended stellar component and those with a compact stellar component, the so-called blue nuggets., Comment: Published in A&A. 18 pages, 14 figures

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs., Comment: Accepted by A&A. 13 pages, 9 figures

During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far beyond pre-JWST limits. While generally robust, such photometric redshifts can suffer from degeneracies and occasionally catastrophic errors. Spectroscopic measurement is required to validate these sources and to reliably quantify physical properties that can constrain galaxy formation models and cosmology. Here we present JWST spectroscopy that confirms redshifts for two very luminous galaxies with $z > 11$, but also demonstrates that another candidate with suggested $z\approx 16$ instead has $z = 4.9$, with an unusual combination of nebular line emission and dust reddening that mimics the colors expected for much more distant objects. These results reinforce evidence for the early, rapid formation of remarkably luminous galaxies, while also highlighting the necessity of spectroscopic verification. The large abundance of bright, early galaxies may indicate shortcomings in current galaxy formation models, or deviation from physical properties (such as the stellar initial mass function) that are generally believed to hold at later times., Comment: Accelerated Article Preview (AAP) version published in Nature

We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

Infrared (IR), sub-millimetre (sub-mm) and millimetre (mm) databases contain a huge quantity of high quality data. However, a large part of these data are photometric, and are thought not to be useful to derive a quantitative information on the nebular emission of galaxies. The aim of this project is first to identify galaxies at z > 4-6, and in the epoch of reionization from their sub-mm colours. We also aim at showing that the colours can be used to try and derive physical constraints from photometric bands, when accounting for the contribution from the IR fine structure lines to these photometric bands. We model the flux of IR fine structure lines with CLOUDY, and add them to the dust continuum emission with CIGALE. Including or not emission lines in the simulated spectral energy distribution (SED) modifies the broad band emission and colours. The introduction of the lines allows to identify strong star forming galaxies at z > 4 - 6 from the log10 (PSW_250um/PMW_350um) versus log10 (LABOCA_870um/PLW_500um) colour-colour diagramme. By comparing the relevant models to each observed galaxy colour, we are able to roughly estimate the fluxes of the lines, and the associated nebular parameters. This method allows to identify a double sequence in a plot built from the ionization parameter and the gas metallicity. The HII and photodissociation region (PDR) fine structure lines are an essential part of the SEDs. It is important to add them when modelling the spectra, especially at z > 4 - 6 where their equivalent widths can be large. Conversely, we show that we can extract some information on strong IR fine structure lines and on the physical parameters related to the nebular emission from IR colour-colour diagrams., Comment: Paper accepted in Astronomy and Astrophysics on 10 November 2022

We present ALMA Band~7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866~$\mu$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr$^{-1}$. We detect a 5.1$\sigma$ line feature at $338.726\pm0.007$~GHz exactly coinciding with the JWST source position, with a 2\% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$\mu$m at $z=16.01$ or [CII]158$\mu$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11--13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11--17 is favored in every candidate as a very blue (UV continuum slope of $\simeq-2.3$) and luminous ($M_{\rm UV}\simeq[-$24:$-21]$) system. Still, we find in several candidates that reasonable SED fits ($\Delta$ $\chi^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$--5, where such populations predicted from their luminosity functions and EW([OIII]+H$\beta$) distributions are abundant in survey volumes used for the identification of the $z\sim$11--17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in several of the $z\sim$11--17 candidates, indicating the need to consider the relative surface densities of the lower-$z$ contaminants in the ultra high-$z$ galaxy search., Comment: 25 pages, 8 figures, 7 tables. ApJ in press

The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3

We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state., Comment: Submitted to ApJ; 12 Pages, 6 Figures

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain \textit{absolute} spectrophotometry of the current version of the reductions. Compared to $z\sim3$ galaxies in the literature, the $z>5$ galaxies have similar [OIII]$\lambda$5008/H$\beta$ ratios, similar [OIII]$\lambda$4364/H$\gamma$ ratios, and higher ($\sim$0.5 dex) [NeIII]$\lambda$3870/[OII]$\lambda$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$\lambda$3870/[OII]$\lambda$3728, [OIII]$\lambda$4364/H$\gamma$, and [OIII]$\lambda$5008/H$\beta$ emission-line ratios are consistent with an interstellar medium that has very high ionization ($\log(Q) \simeq 8-9$, units of cm~s$^{-1}$), low metallicity ($Z/Z_\odot \lesssim 0.2$), and very high pressure ($\log(P/k) \simeq 8-9$, units of cm$^{-3}$). The combination of [OIII]$\lambda$4364/H$\gamma$ and [OIII]$\lambda$(4960+5008)/H$\beta$ line ratios indicate very high electron temperatures of $4.1<\log(T_e/{\rm K})<4.4$, further implying metallicities of $Z/Z_\odot \lesssim 0.2$ with the application of low-redshift calibrations for ``$T_e$-based'' metallicities. These observations represent a tantalizing new view of the physical conditions of the interstellar medium in galaxies at cosmic dawn., Comment: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

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., Comment: 20 pages, 14 figures. Accepted for publication in A&A

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(\theta) \propto 1+\cos \theta$, where $\theta$ 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, Comment: ApJ accepted. The code is available at https://cigale.lam.fr/news/

In this work we study how to employ the upcoming Legacy Survey of Space and Time (LSST) data to constrain physical properties of normal, star forming galaxies. We use simulated LSST data and existing real observations to test the estimations of the physical properties of galaxies, such as star formation rate (SFR), stellar mass ($M_{star}$), and dust luminosity ($L_{dust}$). We focus on normal star-forming galaxies, as they form the majority of the galaxy population in the universe and therefore are more likely to be observed by the LSST. We perform a simulation of LSST observations and uncertainties of 50,385 real galaxies within redshift range $0

We investigate the timescale with which the IR luminosity decreases after a complete and rapid quenching of star formation using observations of local and high-redshift galaxies. From SED modelling, we derive the time since quenching of a subsample of 14 galaxies from the Herschel Reference Survey suffering from ram-pressure stripping due to the environment of the Virgo cluster and of a subsample of 7 rapidly quenched COSMOS galaxies selected through a state-of-the-art statistical method already tested on the determination of galaxies' star formation history. Three out of the 7 COSMOS galaxies have an optical spectra with no emission line, confirming their quenched nature. Present physical properties of the two samples are obtained as well as the past L$_{IR}$ of these galaxies, just before their quenching, from the long-term SFH properties. This past L$_{IR}$ is shown to be consistent with the L$_{IR}$ of reference samples of normally star-forming galaxies with same $M_*$ and $z$ than each of our quenched galaxies. We put constraints on the present to past L$_{IR}$ ratio as a function of quenching time. The two samples probe different dynamical ranges in terms of quenching age with the HRS galaxies exhibiting longer timescales (0.2-3\,Gyr) compared to the COSMOS one ($<100$\,Myr). Assuming an exponential decrease of the L$_{IR}$ after quenching, the COSMOS quenched galaxies are consistent with short e-folding times less than a couple of hundreds of Myr while the properties of the HRS quenched galaxies are compatible with timescales of several hundreds of Myr. For the HRS sample, this result is consistent with ram pressure stripping due to the environment. For the COSMOS sample, different quenching processes are acting on short to intermediate timescales. Processes such as galaxy mergers, disk instabilities or environmental effects can produce such strong star formation variability., Comment: 14 pages, 14 figures, accepted for publication in A&A

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\%$., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters. The CYGNUS models used in this work are available at https://arc.euc.ac.cy/cygnus-project-arc/

We revisit the dependence of covering factor (CF) of dust torus on physical properties of active galactic nuclei (AGNs) by taking into account an AGN polar dust emission. The CF is converted from a ratio of infrared (IR) luminosity contributed from AGN dust torus ($L_{\rm IR}^{\rm torus}$) and AGN bolometric luminosity ($L_{\rm bol}$), by assuming a non-linear relation between luminosity ratio and intrinsic CF. We select 37,181 type 1 quasars at $z < 0.7$ from the Sloan Digital Sky Survey Data Release 16 quasar catalog. Their $L_{\rm bol}$, black hole mass ($M_{\rm BH}$), and Eddington ratio ($\lambda_{\rm Edd}$) are derived by spectral fitting with QSFit. We conduct spectral energy distribution decomposition by using X-CIGALE with clumpy torus and polar dust model to estimate $L_{\rm IR}^{\rm torus}$ without being affected by the contribution of stellar and AGN polar dust to IR emission. For 5720 quasars whose physical quantities are securely determined, we perform a correlation analysis on CF and (i) $L_{\rm bol}$, (ii) $M_{\rm BH}$, and (iii) $\lambda_{\rm Edd}$. As a result, anti-correlations for CF-$L_{\rm bol}$, CF-$M_{\rm BH}$, and CF-$\lambda_{\rm Edd}$ are confirmed. We find that incorporating the AGN polar dust emission makes those anti-correlations stronger which are compared to those without considering it. This indicates that polar dust wind provably driven by AGN radiative pressure is one of the key components to regulate obscuring material of AGNs., Comment: 20 pages, 14 figures, and 4 tables, accepted for publication in ApJ. Physical properties (e.g., BH mass, Eddington ratio, IR luminosity) of 37,181 quasars at z < 0.7 selected from the SDSS DR16 will be available as a Machine-readable table

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\mu 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., Comment: 31 pages, 24 figures, 3 tables. Re-submitted to ApJS after revision

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

We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170 million objects is constructed by a positional cross match of 51 public surveys. This high resolution optical, near-infrared, and mid-infrared catalogue is designed for photometric redshift estimation, extraction of fluxes in lower resolution far-infrared maps, and spectral energy distribution modelling. It collates, standardises, and provides value added derived quantities including corrected aperture magnitudes and astrometry correction over the Herschel extragalactic wide fields for the first time. $grizy$ fluxes are available on all fields with $g$ band data reaching $5\sigma$ point-source depths in a 2 arcsec aperture of 23.5, 24.4, and 24.6 (AB) mag at the 25th, 50th, and 75th percentiles, by area covered, across all HELP fields. It has $K$ or $K_s$ coverage over 1146 deg.$^2$ with depth percentiles of 20.2, 20.4, and 21.0 mag respectively. The IRAC Ch 1 band is available over 273 deg.$^2$ with depth percentiles of 17.7, 21.4, and 22.2 mag respectively. This paper defines the "masterlist" objects for the first data release (DR1) of HELP. This large sample of standardised total and corrected aperture fluxes, uniform quality flags, and completeness measures provides large well understood statistical samples over the full Herschel extragalactic sky.

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., Comment: International Astronomical Union Proceedings: IAU Symposium 341 - PanModel2018: Challenges in Panchromatic Galaxy Modelling with next generation facilities

Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2., Comment: 10 pages, 6 figures, accepted for publication in PASK AKARI special issue. A summary talk is at https://youtu.be/ZIHEAYj0bqQ

We present a study of 7 star-forming galaxies from the Cosmic Evolution Survey (COSMOS) observed with the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope (HST). The galaxies are located at relatively low redshifts, $z\sim$0.3, with morphologies ranging from extended and disturbed to compact and smooth. To complement the HST observations we also analyze observations taken with the VIMOS spectrograph on the Very Large Telescope (VLT). In our galaxy sample we identify three objects with double peak Lyman-$\alpha$ profiles similar to those seen in Green Pea compact galaxies and measure peak separations of 655, 374, and 275 km s$^{-1}$. We measure Lyman-$\alpha$ escape fractions with values ranging between 5-13\%. Given the low flux levels in the individual COS exposures we apply a weighted stacking approach to obtain a single spectrum. From this COS combined spectrum we infer upper limits for the absolute and relative Lyman continuum escape fractions of $f_{\rm abs}(\rm LyC)$ = 0.4$^{+10.1}_{-0.4}$\% and $f_{\rm res}(\rm LyC)$ = 1.7$^{+15.2}_{-1.7}$\%, respectively. Finally, we find that most of these galaxies have moderate UV and optical SFRs (SFRs $\lesssim$ 10 M$_{\odot}$ yr$^{-1}$)., Comment: 15 pages, 10 figures, MNRAS accepted

We present a study of photometric redshift performance for galaxies and active galactic nuclei detected in deep radio continuum surveys. Using two multi-wavelength datasets, over the NOAO Deep Wide Field Survey Bo\"otes and COSMOS fields, we assess photometric redshift (photo-z) performance for a sample of $\sim 4,500$ radio continuum sources with spectroscopic redshifts relative to those of $\sim 63,000$ non radio-detected sources in the same fields. We investigate the performance of three photometric redshift template sets as a function of redshift, radio luminosity and infrared/X-ray properties. We find that no single template library is able to provide the best performance across all subsets of the radio detected population, with variation in the optimum template set both between subsets and between fields. Through a hierarchical Bayesian combination of the photo-z estimates from all three template sets, we are able to produce a consensus photo-z estimate which equals or improves upon the performance of any individual template set., Comment: 20 pages, 17 figures. Accepted for publication in MNRAS

In most studies of dust in galaxies, dust is only detected from its emission to approximately the optical radius of the galaxy. By combining the signal of 110 spiral galaxies observed as part of the Herschel Reference Survey, we are able to improve our sensitivity by an order-of-magnitude over that for a single object. Here we report the direct detection of dust from its emission that extends out to at least twice the optical radius. We find that the distribution of dust is consistent with an exponential at all radii with a gradient of ~-1.7 dex R$_{25}^{-1}$. Our dust temperature declines linearly from ~25 K in the centre to 15 K at R$_{25}$ from where it remains constant out to ~2.0 R$_{25}$. The surface-density of dust declines with radius at a similar rate to the surface-density of stars but more slowly than the surface-density of the star-formation rate. Studies based on dust extinction and reddening of high-redshift quasars have concluded that there are substantial amounts of dust in intergalactic space. By combining our results with the number counts and angular correlation function from the SDSS, we show that with Milky Way type dust we can explain the reddening of the quasars by the dust within galactic disks alone. Given the uncertainties in the properties of any intergalactic dust, we cannot rule out its existence, but our results show that statistical investigations of the dust in galactic halos that use the reddening of high-redshift objects must take account of the dust in galactic disks., Comment: 15 Pages, 7 Figures. Accepted to MNRAS 2016 July 4. Received 2016 July 4; in original form 2015 December 6

Gas flows play a fundamental role in galaxy formation and evolution, providing the fuel for the star formation process. These mechanisms leave an imprint in the amount of heavy elements. Thus, the analysis of this metallicity signature provides additional constraint on the galaxy formation scenario. We aim to discriminate between four different galaxy formation models based on two accretion scenarios and two different star formation recipes. We address the impact of a bimodal accretion scenario and a strongly regulated star formation recipe. We present a new extension of the eGalICS model, which allows us to track the metal enrichment process. Our new chemodynamical model is applicable for situations ranging from metal-free primordial accretion to very enriched interstellar gas contents. We use this new tool to predict the metallicity evolution of both the stellar populations and gas phase. We also address the evolution of the gas metallicity with the star formation rate (SFR). We then focus on a sub-sample of Milky Way-like galaxies. We compare both the cosmic stellar mass assembly and the metal enrichment process of such galaxies with observations and detailed chemical evolution models. Our models, based on a strong star formation regulation, allow us to reproduce well the stellar mass to gas-phase metallicity relation observed in the local universe. However, we observe a systematic shift towards high masses. Our $Mstar-Zg-SFR relation is in good agreement with recent measurements: our best model predicts a clear dependence with the SFR. Both SFR and metal enrichment histories of our Milky Way-like galaxies are consistent with observational measurements and detailed chemical evolution models. We finally show that Milky Way progenitors start their evolution below the observed main sequence and progressively reach this observed relation at z = 0., Comment: 22 pages, 11 figures

We present Herschel PACS $\&$ SPIRE far-infrared (FIR) and sub-mm imaging observations for a large K-band selected sample of 88 close major-merger pairs of galaxies (H-KPAIRs) in 6 photometric bands (70, 100, 160, 250, 350, and 500 $\mu m$). Among 132 spiral galaxies in the 44 spiral-spiral (S$+$S) pairs and 44 spiral-elliptical (S$+$E) pairs, 113 are detected in at least one Herschel band. Star formation rate (SFR) and dust mass ($M_{\rm dust}$) are derived from the IR SED fitting. Mass of total gas ($M_{\rm gas}$) is estimated by assuming a constant dust-to-gas mass ratio of 0.01. Star forming spiral galaxies (SFGs) in S$+$S pairs show significant enhancements in both specific star formation rate (sSFR) and star formation efficiency (SFE), while having nearly the same gas mass, compared to control galaxies. On the other hand, for SFGs in S$+$E pairs, there is no significant sSFR enhancement and the mean SFE enhancement is significantly lower than that of SFGs in S$+$S pairs. This suggests an important role for the disc-disc collision in the interaction induced star formation. The $M_{\rm gas}$ of SFGs in S$+$E pairs is marginally lower than that of their counterparts in both S$+$S pairs and the control sample. Paired galaxies with and without interaction signs do not differ significantly in their mean sSFR and SFE. As found in previous works, this much larger sample confirms the primary and secondary spirals in S+S pairs follow a Holmberg effect correlation on sSFR., Comment: 21 pages, 14 figures, 5 tables. Accepted for publication in ApJS

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$\mu$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$\mu$m and 12$\mu$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$\mu$m, 12$\mu$m, and TIR LFs at $0.15

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes., Comment: 17 pages, 14 figures, 2 tables, accepted by A&A

Context. Measuring star formation across the Universe is key to constrain models of galaxy formation and evolution. Yet, determining the SFR (star formation rate) of galaxies remains a challenge. Aims. In this paper we investigate in isolation the impact of a variable star formation history on the measurement of the SFR. Methods. We combine 23 state-of-the-art hydrodynamical simulations of 1

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr., Comment: 30 pages, 26 figures, accepted for publication in A&A

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts., Comment: 19 pages, 12 figure, accepted for publication in ApJ

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)., Comment: arXiv admin note: substantial text overlap with arXiv:1304.7000

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.

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies., Comment: Accepted for publication in ApJ. 24 pages, 11 figures. High resolution version of the figures, and composite SEDs presented in the paper, are available at http://www.astro.umass.edu/~pope/Kirkpatrick2012

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., Comment: White paper submitted to the HST / Herschel Deep Fields working group. CANDELS reference added

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content., Comment: 9 pages, 9 figures. Accepted by ApJ

We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 {\mu}Jy. Material of the model (predictions, SED library, mock catalogs...) is available online at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537., Comment: 7 pages, 4 figures, 1 table, ApJL in press

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 $\rho$ 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., Comment: 10 pages, 7 figures, Earth, Planets and Space, in press

We present CIGALE (Burgarella et al. 2005; Noll et al. 2009), a software developed at the Laboratoire d'Astrophysique de Marseille to fit galaxy spectral energy distributions from the rest-frame far-UV to far-IR wavelength range, and to derive some of their physical parameters. We also give some examples of scientific results obtained with CIGALE., Comment: 4 pages, 2 figures, proceeding of a poster to be presented at ADASS XXI in Paris

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., Comment: 12 pages, 8 figures, 4 tables, updated to match the version accepted for publication in ApJ; code available at http://www.oamp.fr/cigale/

We present a new analysis of the dust obscuration in starburst galaxies at low and high redshift. This study is motivated by our unique sample of the most extreme UV-selected starburst galaxies in the nearby universe (z<0.3), found to be good analogs of high-redshift Lyman Break Galaxies (LBGs) in most of their physical properties. We find that the dust properties of the Lyman Break Analogs (LBAs) are consistent with the relation derived previously by Meurer et al. (M99) that is commonly used to dust-correct star formation rate measurements at a very wide range of redshifts. We directly compare our results with high redshift samples (LBGs, BzK, and sub-mm galaxies at z=2-3) having IR data either from Spitzer or Herschel. The attenuation in typical LBGs at z=2-3 and LBAs is very similar. Because LBAs are much better analogs to LBGs compared to previous local star-forming samples, including M99, the practice of dust-correcting the SFRs of high redshift galaxies based on the local calibration is now placed on a much more solid ground. We illustrate the importance of this result by showing how the locally calibrated relation between UV measurements and extinction is used to estimate the integrated, dust-corrected star formation rate density at z=2-6., Comment: Accepted for publication in ApJ Letters (6 pages, 4 figures, 1 table)

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., Comment: White paper supplement to the proposal submitted by the HSLS science team to ESA for Herschel open-time programs

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 \mu\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 \mu\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.