Showing 1–50 of 222 results for author: Burgarella, D

Rapid Response Mode observations of GRB 160203A: Looking for fine-structure line variability at z=3.52

Authors: G. Pugliese, A. Saccardi, V. D Elia, S. D. Vergani, K. E. Heintz, S. Savaglio, L. Kaper, A. de Ugarte Postigo, D. H. Hartmann, A. De Cia, S. Vejlgaard, J. P. U. Fynbo, L. Christensen, S. Campana, D. van Rest, J. Selsing, K. Wiersema, D. B. Malesani, S. Covino, D. Burgarella, M. De Pasquale, P. Jakobsson, J. Japelj, D. A. Kann, C. Kouveliotou , et al. (4 additional authors not shown)

Abstract: Gamma-ray bursts are the most energetic known explosions. Despite fading rapidly, they allow to measure redshift and important properties of their host-galaxies. We report the photometric and spectroscopic study of GRB 160203A and its host-galaxy. Fine-structure absorption lines, detected in the afterglow at different epochs, allow us to investigate variability due to the strong fading background… ▽ More Gamma-ray bursts are the most energetic known explosions. Despite fading rapidly, they allow to measure redshift and important properties of their host-galaxies. We report the photometric and spectroscopic study of GRB 160203A and its host-galaxy. Fine-structure absorption lines, detected in the afterglow at different epochs, allow us to investigate variability due to the strong fading background source. We obtained two optical to near-infrared spectra of the afterglow with X-shooter on ESO/VLT, 18 min and 5.7 hrs after the burst, allowing us to investigate temporal changes of fine-structure absorption lines. We measured HI column density log N(HI/cm-2)=21.75+/-0.10, and several heavy-element ions along the GRB sight-line in the host-galaxy: SiII,AlII,AlIII,CII,NiII,SiIV,CIV,ZnII,FeII, and FeII and SiII fine structure transitions from energetic levels excited by the afterglow, at a redshift z=3.518. We measured [M/H]TOT=-0.78+/-0.13 and [Zn/Fe]FIT=0.69+/-0.15, representing the total(dust-corrected) metallicity and dust depletion, respectively. We detected additional intervening systems along the line of sight at z=1.03,z=1.26,z=1.98,z=1.99,z=2.20 and z=2.83. We could not measure significant variability in the fine-structure lines throughout all the observations and determined an upper limit for the GRB distance from the absorber of d<300 pc, adopting the canonical UV pumping scenario. However, we note that the quality of our data is not sufficient to conclusively rule out collisions as an alternative mechanism. GRB 160203A belongs to a growing sample of GRBs with medium resolution spectroscopy, provided by the Swift/X-shooter legacy program, which enables detailed investigation of the interstellar medium in high-redshift GRB host-galaxies. In particular, this host galaxy shows relatively high metal enrichment and dust depletion already in place when the universe was only 1.8 Gyr old. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 16 pages, 9 figures, 2 appendices, A&A accepted

Physical properties of strong 1 < z < 3 Balmer and Paschen lines emitters observed with JWST

Authors: L. -M. Seillé, V. Buat, V. Fernández, M. Boquien, Y. Roehlly, A. Boselli, A. Calabrò, R. O. Amorín, B. E. Backhaus, D. Burgarella, N. J. Cleri, M. Dickinson, N. P. Hathi, B. W. Holwerda, A. M. Koekemoer, L. Napolitano, F. Pacucci, C. Robertson, L. Y. A. Yung

Abstract: The ultraviolet continuum traces young stars while the near-infrared unveils older stellar populations and dust-obscured regions. Balmer emission lines provide insights on gas properties and young stellar objects but are highly affected by dust attenuation. The near-infrared Paschen lines suffer less dust attenuation and can be used to measure star formation rates (SFRs) in star-forming regions ob… ▽ More The ultraviolet continuum traces young stars while the near-infrared unveils older stellar populations and dust-obscured regions. Balmer emission lines provide insights on gas properties and young stellar objects but are highly affected by dust attenuation. The near-infrared Paschen lines suffer less dust attenuation and can be used to measure star formation rates (SFRs) in star-forming regions obscured by dust clouds. We select 13 sources between redshifts 1 and 3 observed with HST, JWST/NIRCam and NIRSpec based on the availability of at least one Balmer and one Paschen line with S/N > 5. With a newly-developed version of CIGALE, we fit their hydrogen line equivalent widths (EWs) and photometric data. We assess the impacts of the removal of spectroscopic data by comparing the quality of the fits of the spectro-photometric data to those with photometric data only. We compare the single (BC03) vs binary (BPASS) stellar populations models in the fitting process of spectro-photometric data. We derive the differential attenuation and explore different attenuation recipes by fitting spectro-photometric data with BC03. For each stellar model and for each input dataset (with and without EWs), we quantify the deviation on the SFRs and stellar masses from the "standard" choice. On average, the SFRs are overestimated and the stellar masses are underestimated when EWs are not included as input data. We find a major contribution of the H$α$ emission line to the broadband photometric measurements of our sources, and a trend of increasing contribution with specific SFR. Using the BPASS models has a significant impact on the derived SFRs and stellar masses. We show that a flexible attenuation recipe provides more accurate estimates of the dust attenuation parameters, especially the differential attenuation which agrees with the original value of Charlot & Fall (2000). △ Less

Submitted 16 April, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: Submitted to Astronomy & Astrophysics, 19 pages, 15 figures

Overcoming Confusion Noise with Hyperspectral Imaging from PRIMAger

Authors: James M. S. Donnellan, Seb J. Oliver, Matthieu Bethermin, Longji Bing, Alberto Bolatto, Charles M. Bradford, Denis Burgarella, Laure Ciesla, Jason Glenn, Alexandra Pope, Stephen Serjeant, Raphael Shirley, JD T. Smith, Chris Sorrell

Abstract: The PRobe far-Infrared Mission for Astrophysics (PRIMA) concept aims to perform mapping with spectral coverage and sensitivities inaccessible to previous FIR space telescopes. PRIMA's imaging instrument, PRIMAger, provides unique hyperspectral imaging simultaneously covering 25-235 $μ$m. We synthesise images representing a deep, 1500 hr deg$^{-2}$ PRIMAger survey, with realistic instrumental and c… ▽ More The PRobe far-Infrared Mission for Astrophysics (PRIMA) concept aims to perform mapping with spectral coverage and sensitivities inaccessible to previous FIR space telescopes. PRIMA's imaging instrument, PRIMAger, provides unique hyperspectral imaging simultaneously covering 25-235 $μ$m. We synthesise images representing a deep, 1500 hr deg$^{-2}$ PRIMAger survey, with realistic instrumental and confusion noise. We demonstrate that we can construct catalogues of galaxies with a high purity ($>95$ per cent) at a source density of 42k deg$^{-2}$ using PRIMAger data alone. Using the XID+ deblending tool we show that we measure fluxes with an accuracy better than 20 per cent to flux levels of 0.16, 0.80, 9.7 and 15 mJy at 47.4, 79.7, 172, 235 $μ$m respectively. These are a factor of $\sim$2 and $\sim$3 fainter than the classical confusion limits for 72-96 $μ$m and 126-235 $μ$m, respectively. At $1.5 \leq z \leq 2$, we detect and accurately measure fluxes in 8-10 of the 10 channels covering 47-235 $μ$m for sources with $2 \leq$ log(SFR) $\leq 2.5$, a 0.5 dex improvement on what might be expected from the classical confusion limit. Recognising that PRIMager will operate in a context where high quality data will be available at other wavelengths, we investigate the benefits of introducing additional prior information. We show that by introducing even weak prior flux information when employing a higher source density catalogue (more than one source per beam) we can obtain accurate fluxes an order of magnitude below the classical confusion limit for 96-235 $μ$m. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 14 pages, 11 figures

Confusion of extragalactic sources in the far infrared: a baseline assessment of the performance of PRIMAger in intensity and polarization

Authors: Matthieu Béthermin, Alberto D. Bolatto, François Boulanger, Charles M. Bradford, Denis Burgarella, Laure Ciesla, James Donnellan, Brandon S. Hensley, Jason Glenn, Guilaine Lagache, Enrique Lopez-Rodriguez, Seb Oliver, Alexandra Pope, Marc Sauvage

Abstract: Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in in… ▽ More Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in intensity and polarization. We aim to provide predictions of the confusion level and its consequences for future surveys. We produced simulated PRIMAger maps affected only by the confusion noise using the simulated infrared extragalactic sky (SIDES) semi-empirical simulation. We then estimated the confusion limit in these maps and extracted the sources using a basic blind extractor. By comparing the input galaxy catalog and the extracted source catalog, we derived various performance metrics as completeness, purity, and the accuracy of various measurements. In intensity, we predict that the confusion limit increases rapidly with increasing wavelength. The confusion limit in polarization is more than 100x lower. The measured flux density is dominated by the brightest galaxy in the beam, but other objects also contribute at longer wavelength (~30% at 235 $μ$m). We also show that galaxy clustering has a mild impact on confusion in intensity (up to 25%), while it is negligible in polarization. In intensity, a basic blind extraction will be sufficient to detect galaxies at the knee of the luminosity function up to z~3 and 10$^{11}$ M$_\odot$ main-sequence galaxies up to z~5. In polarization for a conservative sensitivity, we expect ~8 000 detections up to z=2.5 opening a totally new window on the high-z dust polarization. Finally, we show that intensity surveys at short wavelength and polarization surveys at long wavelength tend to reach confusion at similar depth. There is thus a strong synergy. △ Less

Submitted 28 April, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Comments: 21 pages, 15 figures, 4 tables, submitted to A&A, minor update on the names of the sensitivity scenarii (for consistency with other papers)

Multi-purpose InSTRument for Astronomy at Low-resolution: MISTRAL@OHP

Authors: J. Schmitt, C. Adami, M. Dennefeld, F. Agneray, S. Basa, J. C. Brunel, V. Buat, D. Burgarella, C. Carvalho, G. Castagnoli, N. Grosso, F. Huppert, C. Moreau, F. Moreau, L. Moreau, E. Muslimov, S. Pascal, S. Perruchot, D. Russeil, J. L. Beuzit, F. Dolon, M. Ferrari, B. Hamelin, A. LevanSuu, K. Aravind , et al. (9 additional authors not shown)

Abstract: MISTRAL is the new Faint Object Spectroscopic Camera mounted at the folded Cassegrain focus of the 1.93m telescope of Haute-Provence Observatory. We describe the design and components of the instrument and give some details about its operation. We emphasise in particular the various observing modes and the performances of the detector. A short description is also given about the working environmen… ▽ More MISTRAL is the new Faint Object Spectroscopic Camera mounted at the folded Cassegrain focus of the 1.93m telescope of Haute-Provence Observatory. We describe the design and components of the instrument and give some details about its operation. We emphasise in particular the various observing modes and the performances of the detector. A short description is also given about the working environment. Various types of objects, including stars, nebulae, comets, novae, galaxies have been observed during various test phases to evaluate the performances of the instrument. The instrument covers the range of 4000 to 8000A with the blue setting, or from 6000 to 10000A with the red setting, at an average spectral resolution of 700. Its peak efficiency is about 22% at 6000A. In spectroscopy, a limiting magnitude of 19.5 can be achieved for a point source in one hour with a signal to noise of 3 in the continuum (and better if emission lines are present). In imaging mode, limiting magnitudes of 20-21 can be obtained in 10-20mn (with average seing conditions of 2.5 arcsec at OHP). The instrument is very users-friendly and can be put into operations in less than 15mn (rapid change-over from the other instrument in use) if required by the science (like for Gamma-Rays Bursts). Some first scientific results are described for various types of objects, and in particular for the follow-up of GRBs. While some further improvements are still under way, in particular to ease the switch from blue to red setting and add more grisms or filters, MISTRAL is ready for the follow-up of transients and other variable objects, in the soon-to-come era of e.g. the SVOM satellite and of the Rubin telescope. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: Accepted in A&A

Evidence of extreme ionization conditions and low metallicity in GHZ2/GLASS-z12 from a combined analysis of NIRSpec and MIRI observations

Authors: Antonello Calabro, Marco Castellano, Jorge A. Zavala, Laura Pentericci, Pablo Arrabal Haro, Tom J. L. C. Bakx, Denis Burgarella, Caitlin M. Casey, Mark Dickinson, Steven L. Finkelstein, Adriano Fontana, Mario Llerena, Sara Mascia, Emiliano Merlin, Ikki Mitsuhashi, Lorenzo Napolitano, Diego Paris, Pablo G. Perez-Gonzalez, Guido Roberts-Borsani, Paola Santini, Tommaso Treu, Eros Vanzella

Abstract: GHZ2/GLASS-z12 has been recently observed by JWST with both NIRSpec and MIRI spectrographs, making it the most distant galaxy ($z_{spec}=12.34$) with complete spectroscopic coverage from rest-frame UV to optical. It is identified as a strong CIV$_{1549}$ emitter with many other detected emission lines (NIV], HeII, OIII], NIII], CIII], [OII], [NeIII], [OIII], and H$α$), including a remarkable OIII… ▽ More GHZ2/GLASS-z12 has been recently observed by JWST with both NIRSpec and MIRI spectrographs, making it the most distant galaxy ($z_{spec}=12.34$) with complete spectroscopic coverage from rest-frame UV to optical. It is identified as a strong CIV$_{1549}$ emitter with many other detected emission lines (NIV], HeII, OIII], NIII], CIII], [OII], [NeIII], [OIII], and H$α$), including a remarkable OIII$_{1333}$ Bowen fluorescence line. We analyze in this paper the joint NIRSpec+MIRI spectral data set. Combining six optical diagnostics (namely R2, R3, R23, O32, Ne3O2, and Ne3O2Hd), we find extreme ionization conditions, with O32 $=1.39 \pm 0.19$ and Ne3O2 $=0.37 \pm 0.18$ in stark excess compared to typical values in the ISM at lower redshifts. These line properties are compatible either with an AGN or with a compact, very dense star-forming environment ($Σ_{\rm SFR}$ $\sim 10^2$-$10^3$ Msun/yr/kpc$^2$), with a high ionization parameter ($\log_{10}$(U) $=-1.75 \pm 0.16$), a high ionizing photon production efficiency $\log(ξ_{\rm ion}) = 25.7_{-0.1}^{+0.2}$, and a low, although not pristine, metal content ranging between $5\%$ and $11\%$ Z$_\odot$, indicating a rapid enrichment of the ISM in the last few Myrs. These properties also suggest that a substantial amount of ionizing photons ($\sim 10\%$) are leaking outside. The general lessons learned from GHZ2 are the following: (i) the UV to optical combined nebular indicators are broadly in agreement with UV-only or optical-only indicators. (ii) UV+optical diagnostics fail to discriminate between an AGN and star-formation in a low metallicity, high density, and extreme ionization environment. (iii) comparing the nebular line ratios with local analogs may be approaching its limits at $z \gtrsim 10$, as this approach is potentially challenged by the unique conditions of star formation experienced by galaxies at these extreme redshifts. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 14 pages, 3 figures, and 1 table. Submitted to ApJ

Detection of ionized hydrogen and oxygen from a very luminous and young galaxy 13.4 billion years ago

Authors: Jorge A. Zavala, Marco Castellano, Hollis B. Akins, Tom J. L. C. Bakx, Denis Burgarella, Caitlin M. Casey, Óscar A. Chávez Ortiz, Mark Dickinson, Steven L. Finkelstein, Ikki Mitsuhashi, Kimihiko Nakajima, Pablo G. Pérez-González, Pablo Arrabal Haro, Veronique Buat, Bren Backhaus, Antonello Calabrò, Nikko J. Cleri, David Fernández-Arenas, Adriano Fontana, Maximilien Franco, Mauro Giavalisco, Norman A. Grogin, Nimish Hathi, Michaela Hirschmann, Ryota Ikeda , et al. (12 additional authors not shown)

Abstract: 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… ▽ More 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α) and doubly-ionized oxygen ([OIII]4959,5007Å) 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. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 14 pages, 6 figures. Submitted for publication. Comments are most welcome

Galaxy main sequence and properties of low-mass Lyman-alpha Emitters towards reionisation viewed by VLT/MUSE and JWST/NIRCam

Authors: I. Goovaerts, R. Pello, D. Burgarella, T. T. Thai, J. Richard, A. Claeyssens, P. Tuan-Anh, R. C. Arango-Toro, L. Boogaard, T. Contini, Y. Guo, I. Langan, N. Laporte, M. Maseda

Abstract: Faint, star-forming galaxies likely play a dominant role in cosmic reionisation. Strides have been made in recent years to characterise these populations at high redshifts ($z>3$). Now for the first time, with JWST photometry beyond 1$\,μm$ in the rest frame, we can derive accurate stellar masses and position these galaxies on the galaxy main sequence. We seek to assess the place of 96 individual… ▽ More Faint, star-forming galaxies likely play a dominant role in cosmic reionisation. Strides have been made in recent years to characterise these populations at high redshifts ($z>3$). Now for the first time, with JWST photometry beyond 1$\,μm$ in the rest frame, we can derive accurate stellar masses and position these galaxies on the galaxy main sequence. We seek to assess the place of 96 individual Lyman-alpha emitters (LAEs) selected behind the A2744 lensing cluster with MUSE spectroscopy on the galaxy main sequence. We also compare derived stellar masses to Lyman-alpha luminosities and equivalent widths to better quantify the relationship between the Lyman-alpha emission and the host galaxy. These 96 LAEs lie in the redshift range $2.9<z<6.7$, and their range of masses extends down to $10^6\,\mathrm{M_{\odot}}$ (over half with $\mathrm{M_{\star}}<10^8\,\mathrm{M_{\odot}}$). We use the JWST/NIRCam and HST photometric catalogs from the UNCOVER project, giving us excellent wavelength coverage from $450\,\mathrm{nm}$ to $4.5\,μm$. We find a main sequence relation for these low mass LAEs of the form: $\mathrm{log\,SFR}=(0.88\pm0.07 - 0.030\pm0.027\times t)\,\mathrm{log\,M_{\star}} - ( 6.31\pm0.41 - 0.08\pm0.37\times t)$. This is in approximate agreement with best-fits of previous collated studies, however, with a steeper slope and a higher normalisation. This indicates that low-mass LAEs towards the epoch of reionisation lie above typical literature main sequence relations derived at lower redshift and higher masses. Additionally, comparing our results to UV-selected samples, we see that while low-mass LAEs lie above these typical main sequence relations, they are likely not singular in this respect at these masses and redshifts. While low-mass galaxies have been shown to play a significant role in cosmic reionisation, our results point to no special position for LAEs in this regard. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 10 pages, 5 figures, accepted for publication in Astronomy and Astrophysics

CEERS: 7.7 $μ$m PAH Star Formation Rate Calibration with JWST MIRI

Authors: Kaila Ronayne, Casey Papovich, Guang Yang, Lu Shen, Mark Dickinson, Robert Kennicutt, Anahita Alavi, Pablo Arrabal Haro, Micaela Bagley, Denis Burgarella, Aurélien Le Bail, Eric Bell, Nikko Cleri, Justin Cole, Luca Costantin, Alexander de la Vega, Emanuele Daddi, David Elbaz, Steven Finkelstein, Norman Grogin, Benne Holwerda, Jeyhan Kartaltepe, Allison Kirkpatrick, Anton Koekemoer, Ray Lucas , et al. (11 additional authors not shown)

Abstract: We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 $μ$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 $μ$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy di… ▽ More We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 $μ$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 $μ$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy distribution (SED) modeling of these data to measure dust-corrected far-UV (FUV) luminosities, $L_{FUV}$, and UV-derived SFRs. We then fit SED models to the JWST/MIRI 7.7 - 21 $μ$m CEERS data to derive rest-frame 7.7 $μ$m luminosities, $L_{770}$, using the average flux density in the rest-frame MIRI F770W bandpass. We observe a correlation between $L_{770}$ and $L_{FUV}$, where log $L_{770}$ is proportional to (1.27+/-0.04) log $L_{FUV}$. $L_{770}$ diverges from this relation for galaxies at lower metallicities, lower dust obscuration, and for galaxies dominated by evolved stellar populations. We derive a "single-wavelength" SFR calibration for $L_{770}$ which has a scatter from model estimated SFRs (${σ_{ΔSFR}}$) of 0.24 dex. We derive a "multi-wavelength" calibration for the linear combination of the observed FUV luminosity (uncorrected for dust) and the rest-frame 7.7 $μ$m luminosity, which has a scatter of ${σ_{ΔSFR}}$ = 0.21 dex. The relatively small decrease in $σ$ suggests this is near the systematic accuracy of the total SFRs using either calibration. These results demonstrate that the rest-frame 7.7 $μ$m emission constrained by JWST/MIRI is a tracer of the SFR for distant galaxies to this accuracy, provided the galaxies are dominated by star-formation with moderate-to-high levels of attenuation and metallicity. △ Less

Submitted 13 October, 2023; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 20 pages, 11 figures, 2 tables, submitted to ApJ

Cosmic Evolution Early Release Science Survey (CEERS): Multi-classing Galactic Dwarf Stars in the deep JWST/NIRCam

Authors: B. W. Holwerda, Chih-Chun Hsu, Nimish Hathi, Laura Bisigello, Alexander de la Vega, Pablo Arrabal Haro, Micaela Bagley, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Casey Papovich, Nor Pirzkal, Kyle Cook, Clayton Robertson, Caitlin M Casey, Christian Aganze, Pablo G. Pérez-González, Ray A. Lucas, Shardha Jogee, Stephen Wilkins, Denis Burgarella, Allison Kirkpatrick

Abstract: Low mass (sub)stellar objects represent the low end of the initial mass function, the transition to free-floating planets and a prominent interloper population in the search for high-redshift galaxies. Without proper motions or spectroscopy, can one identify these objects photometrically? JWST/NIRCam has several advantages over HST/WFC3 NIR: more filters, a greater wavelength range, and greater sp… ▽ More Low mass (sub)stellar objects represent the low end of the initial mass function, the transition to free-floating planets and a prominent interloper population in the search for high-redshift galaxies. Without proper motions or spectroscopy, can one identify these objects photometrically? JWST/NIRCam has several advantages over HST/WFC3 NIR: more filters, a greater wavelength range, and greater spatial resolution. Here, we present a catalogue of (sub)stellar dwarfs identified in the Cosmic Evolution Early Release Science Survey (CEERS). We identify 518 stellar objects down to $m_F200W \sim 28$ using half-light radius, a full three magnitudes deeper than typical HST/WFC3 images. A kNN nearest neighbour algorithm identifies and types these sources, using four HST/WFC3 and four NIRCam filters, trained on SpeX spectra of nearby brown dwarfs. The kNN with four neighbors classifies well within two subtypes: e.g M2$\pm$2 or T4$\pm$2, achieving $\sim$95% precision and recall. More granular typing results in worse metrics. In CEERS, we find 9 M8$\pm$2, 2 L6$\pm$2, 1 T4$\pm$2, and 15 T8$\pm$2. We compare the observed long wavelength NIRCam colours -- not used in the kNN -- to those expected for brown dwarf atmospheric models. The NIRCam F356W-F444W and F410M-F444W colours are redder by a magnitude for the type assigned by the kNN, hinting at a wider variety of atmospheres for these objects. We find a 300-350pc scale-height for M6$\pm$2 dwarfs plus a second structural component and a 150-200pc scale-height for T6$\pm$2 type dwarfs, consistent with literature values. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 18 pages, 24 figures, 3 tables, submitted to MNRAS

The relation of cosmic environment and morphology with the star formation and stellar populations of AGN and non-AGN galaxies

Authors: G. Mountrichas, G. Yang, V. Buat, B. Darvish, M. Boquien, Q. Ni, D. Burgarella, L. Ciesla

Abstract: In this work, we study the relation of cosmic environment and morphology with the star-formation (SF) and the stellar population of galaxies. Most importantly, we examine if this relation differs for systems with active and non-active supermassive black holes. For that purpose, we use 551 X-ray detected active galactic nuclei (AGN) and 16,917 non-AGN galaxies in the COSMOS-Legacy survey, for which… ▽ More In this work, we study the relation of cosmic environment and morphology with the star-formation (SF) and the stellar population of galaxies. Most importantly, we examine if this relation differs for systems with active and non-active supermassive black holes. For that purpose, we use 551 X-ray detected active galactic nuclei (AGN) and 16,917 non-AGN galaxies in the COSMOS-Legacy survey, for which the surface-density field measurements are available. The sources lie at redshift of $\rm 0.3<z<1.2$, probe X-ray luminosities of $\rm 42<log\,[L_{X,2-10keV}(erg\,s^{-1})]<44$ and have stellar masses, $\rm 10.5<log\,[M_*(M_\odot)]<11.5$. Our results show that isolated AGN (field) have lower SFR compared to non AGN, at all L$_X$ spanned by our sample. However, in denser environments (filaments, clusters), moderate L$_X$ AGN ($\rm log\,[L_{X,2-10keV}(erg\,s^{-1})]>43$) and non-AGN galaxies have similar SFR. We, also, examine the stellar populations and the morphology of the sources in different cosmic fields. For the same morphological type, non-AGN galaxies tend to have older stellar populations and are less likely to have undergone a recent burst in denser environments compared to their field counterparts. The differences in the stellar populations with the density field are, mainly, driven by quiescent systems. Moreover, low L$_X$ AGN present negligible variations of their stellar populations, in all cosmic environments, whereas moderate L$_X$ AGN have, on average, younger stellar populations and are more likely to have undergone a recent burst, in high density fields. Finally, in the case of non-AGN galaxies, the fraction of bulge-dominated (BD) systems increases with the density field, while BD AGN are scarce in denser environments. Our results are consistent with a scenario in which a common mechanism, such as mergers, triggers both the SF and the AGN activity. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A. 10 pages, 5 figures, 3 tables

Journal ref: A&A 675, A137 (2023)

JWST CEERS probes the role of stellar mass and morphology in obscuring galaxies

Authors: Carlos Gómez-Guijarro, Benjamin Magnelli, David Elbaz, Stijn Wuyts, Emanuele Daddi, Aurélien Le Bail, Mauro Giavalisco, Mark Dickinson, Pablo G. Pérez-González, Pablo Arrabal Haro, Micaela B. Bagley, Laura Bisigello, Véronique Buat, Denis Burgarella, Antonello Calabrò, Caitlin M. Casey, Yingjie Cheng, Laure Ciesla, Avishai Dekel, Henry C. Ferguson, Steven L. Finkelstein, Maximilien Franco, Norman A. Grogin, Benne W. Holwerda, Shuowen Jin , et al. (16 additional authors not shown)

Abstract: 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 JW… ▽ More 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. △ Less

Submitted 4 September, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: Accepted by A&A. 13 pages, 9 figures

Journal ref: A&A 677, A34 (2023)

Spectroscopic Confirmation of CEERS NIRCam-selected Galaxies at $\boldsymbol{z \simeq 8-10}$

Authors: Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Seiji Fujimoto, Vital Fernández, Jeyhan S. Kartaltepe, Intae Jung, Justin W. Cole, Denis Burgarella, Katherine Chworowsky, Taylor A. Hutchison, Alexa M. Morales, Casey Papovich, Raymond C. Simons, Ricardo O. Amorín, Bren E. Backhaus, Micaela B. Bagley, Laura Bisigello, Antonello Calabrò, Marco Castellano, Nikko J. Cleri, Romeel Davé, Avishai Dekel, Henry C. Ferguson, Adriano Fontana , et al. (23 additional authors not shown)

Abstract: We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a… ▽ More We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a mostly neutral intergalactic medium. The presence (absense) of strong breaks (strong emission lines) give high confidence that these two galaxies are at z>9.6, but the break-derived redshifts have large uncertainties given the low spectral resolution and relatively low signal-to-noise of the CEERS NIRSpec prism data. The two z~10 sources are relatively luminous (M_UV<-20), with blue continua (-2.3<beta<-1.9) and low dust attenuation (A_V=0.15(+0.3,-0.1)); and at least one of them has high stellar mass for a galaxy at that redshift (log(M_*/M_sol)=9.3(+0.2,-0.3)). Considered together with spectroscopic observations of other CEERS NIRCam-selected high-z galaxy candidates in the literature, we find a high rate of redshift confirmation and low rate of confirmed interlopers (8.3%). Ten out of 34 z>8 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z>9.6. We find that z>8 photometric redshifts are generally in agreement (within uncertainties) with the spectroscopic values. However, the photometric redshifts tend to be slightly overestimated (average Delta(z)=0.50+/-0.12), suggesting that current templates do not fully describe the spectra of very high-z sources. Overall, our results solidifies photometric evidence for a high space density of bright galaxies at z>8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z>8. △ Less

Submitted 6 July, 2023; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: Published in ApJL. 19 pages, 9 figures, 7 tables. File with Table 6 included in source .tar file

Journal ref: Astrophys. J. Lett. 951, L22 (2023)

A VLT/VIMOS view of two $Planck$ multiple-cluster systems: structure and galaxy properties

Authors: R. Wicker, N. Aghanim, V. Bonjean, E. Lecoq, M. Douspis, D. Burgarella, E. Pointecouteau

Abstract: We analysed spectroscopic data obtained with VLT-VIMOS for two multiple-cluster systems, PLCKG$214.6+36.9$ and PLCKG$334.8-38.0$, discovered via their thermal Sunyaev-Zel'dovich signal by $Planck$. Combining the Optical spectroscopy, for the redshift determination, and photometric data from galaxy surveys (SDSS, WISE, DESI), we were able to study the structure of the two multiple-cluster systems,… ▽ More We analysed spectroscopic data obtained with VLT-VIMOS for two multiple-cluster systems, PLCKG$214.6+36.9$ and PLCKG$334.8-38.0$, discovered via their thermal Sunyaev-Zel'dovich signal by $Planck$. Combining the Optical spectroscopy, for the redshift determination, and photometric data from galaxy surveys (SDSS, WISE, DESI), we were able to study the structure of the two multiple-cluster systems, to determine their nature and the properties of their member galaxies. We found that the two systems are populated mainly with passive galaxies and that PLCKG$214.6+36.9$ consists of a pair of clusters at redshift $z = 0.445$ and a background isolated cluster at $z = 0.498$, whereas the system PLCKG$334.8-38.0$ is a chance association of three independent clusters at redshifts $z = 0.367$, $z =0.292$, and $z = 0.33$. We also find evidence for remaining star formation activity in the highest-redshift cluster of PLCKG$214.6+36.9$, at $z = 0.498$. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 12 pages, 9 Figures, 5 Tables. Submitted to A&A, comments are welcome

Journal ref: A&A 676, A61 (2023)

Confirmation and refutation of very luminous galaxies in the early universe

Authors: Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Callum T. Donnan, Denis Burgarella, Adam Carnall, Fergus Cullen, James S. Dunlop, Vital Fernández, Seiji Fujimoto, Intae Jung, Melanie Krips, Rebecca L. Larson, Casey Papovich, Pablo G. Pérez-González, Ricardo O. Amorín, Micaela B. Bagley, Véronique Buat, Caitlin M. Casey, Katherine Chworowsky, Seth H. Cohen, Henry C. Ferguson, Mauro Giavalisco, Marc Huertas-Company , et al. (12 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 15 August, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: Accelerated Article Preview (AAP) version published in Nature

CEERS Key Paper VI: JWST/MIRI Uncovers a Large Population of Obscured AGN at High Redshifts

Authors: G. Yang, K. I. Caputi, C. Papovich, P. Arrabal Haro, M. B. Bagley, P. Behroozi, E. F. Bell, L. Bisigello, V. Buat, D. Burgarella, Y. Cheng, N. J. Cleri, R. Dave, M. Dickinson, D. Elbaz, H. C. Ferguson, S. L. Finkelstein, N. A. Grogin, N. P. Hathi, M. Hirschmann, B. W. Holwerda, M. Huertas-Company, T. Hutchison, E. Iani, J. S. Kartaltepe , et al. (13 additional authors not shown)

Abstract: Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN populat… ▽ More Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN population in the distant universe. We estimate the source properties of MIRI-selected objects by utilizing spectral energy distribution (SED) modelling, and classify them into star-forming galaxies (SF), SF-AGN mixed objects, and AGN. The source numbers of these types are 418, 111, and 31, respectively, from 4 MIRI pointings covering $\sim 9$ arcmin$^2$. The sample spans a redshift range of $\approx 0$--5. We derive the median SEDs for all three source types, respectively, and publicly release them. The median MIRI SED of AGN is similar to the typical SEDs of hot dust-obscured galaxies and Seyfert 2s, for which the mid-IR SEDs are dominantly from AGN-heated hot dust. Based on our SED-fit results, we estimate the black-hole accretion density (BHAD; i.e., total BH growth rate per comoving volume) as a function of redshift. At $z<3$, the resulting BHAD agrees with the X-ray measurements in general. At $z>3$, we identify a total of 27 AGN and SF-AGN mixed objects, leading to that our high-$z$ BHAD is substantially higher than the X-ray results ($\sim 0.5$ dex at $z \approx 3$--5). This difference indicates MIRI can identify a large population of heavily obscured AGN missed by X-ray surveys at high redshifts. △ Less

Submitted 15 May, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: ApJL in press, as part of the CEERS Focus Issue. 13 pages and 7 figures (median SEDs in Fig. 5 are available upon request)

A CEERS Discovery of an Accreting Supermassive Black Hole 570 Myr after the Big Bang: Identifying a Progenitor of Massive z > 6 Quasars

Authors: Rebecca L. Larson, Steven L. Finkelstein, Dale D. Kocevski, Taylor A. Hutchison, Jonathan R. Trump, Pablo Arrabal Haro, Volker Bromm, Nikko J. Cleri, Mark Dickinson, Seiji Fujimoto, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Casey Papovich, Nor Pirzkal, Sandro Tacchella, Jorge A. Zavala, Micaela Bagley, Peter Behroozi, Jaclyn B. Champagne, Justin W. Cole, Intae Jung, Alexa M. Morales, Guang Yang, Haowen Zhang, Adi Zitrin , et al. (27 additional authors not shown)

Abstract: We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$α$-break by Hubble and with a Ly$α$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra unc… ▽ More We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$α$-break by Hubble and with a Ly$α$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$α$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$β$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $μ$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch. △ Less

Submitted 29 August, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 26 pages, 14 figures, 6 tables, published in ApJL

Journal ref: The Astrophysical Journal Letters, 2023, Volume 953, Issue 2, id.L29, 26 pp

CEERS Spectroscopic Confirmation of NIRCam-Selected z > 8 Galaxy Candidates with JWST/NIRSpec: Initial Characterization of their Properties

Authors: Seiji Fujimoto, Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Rebecca L. Larson, Denis Burgarella, Micaela B. Bagley, Peter Behroozi, Katherine Chworowsky, Michaela Hirschmann, Jonathan R. Trump, Stephen M. Wilkins, L. Y. Aaron Yung, Anton M. Koekemoer, Casey Papovich, Nor Pirzkal, Henry C. Ferguson, Adriano Fontana, Norman A. Grogin, Andrea Grazian, Lisa J. Kewley, Dale D. Kocevski, Jennifer M. Lotz, Laura Pentericci , et al. (16 additional authors not shown)

Abstract: We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5μ$m either with the NIRSpec prism or its three medium resolution gratings… ▽ More We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5μ$m either with the NIRSpec prism or its three medium resolution gratings. For $z\simeq9$ photometric candidates, we achieve a high confirmation rate of $\simeq$90\%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at $z>10$, where the strong [OIII] and H$β$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$α$ continuum break is not detected with the current sensitivity. Compared with HST-selected bright galaxies ($M_{\rm\,UV}\simeq-22$) that are similarly spectroscopically confirmed at $z\gtrsim8$, these NIRCam-selected galaxies are characterized by lower star formation rates (SFR$\simeq4\,M_{\odot}$~yr$^{-1}$) and lower stellar masses ($\simeq10^{8}\,M_{\odot}$), but with higher [OIII]+H$β$ equivalent widths ($\simeq$1100$Å$), and elevated production efficiency of ionizing photons ($\log(ξ_{\rm\,ion}/{\rm\,Hz\,erg}^{-1})\simeq25.8$) induced by young stellar populations ($<10$~Myrs) accounting for $\simeq20\%$ of the galaxy mass, highlighting the key contribution of faint galaxies to cosmic reionization. Taking advantage of the homogeneous selection and sensitivity, we also investigate metallicity and ISM conditions with empirical calibrations using the [OIII]/H$β$ ratio. We find that galaxies at $z\sim8-9$ have higher SFRs and lower metallicities than galaxies at similar stellar masses at $z\sim2-6$, which is generally consistent with the current galaxy formation and evolution models. △ Less

Submitted 26 May, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

CEERS Key Paper IV: Galaxies at $4 < z < 9$ are Bluer than They Appear -- Characterizing Galaxy Stellar Populations from Rest-Frame $\sim 1$ micron Imaging

Authors: Casey Papovich, Justin Cole, Guang Yang, Steven L. Finkelstein, Guillermo Barro, Véronique Buat, Denis Burgarella, Pablo G. Pérez-González, Paola Santini, Lise-Marie Seillé, Lu Shen, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Laura Bisigello, Antonello Calabrò, Caitlin M. Casey, Marco Castellano, Katherine Chworowsky, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Mark Dickinson, Henry C. Ferguson, Adriano Fontana , et al. (24 additional authors not shown)

Abstract: We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame… ▽ More We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame spectral-energy distribution (SED) to nearly 1 micron for galaxies in this redshift range. By modeling the galaxies' SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å) $-$ $I$-band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower (stellar)-mass-to-light ratios. The MIRI data reduce the stellar masses by $\langle Δ\log M_\ast\rangle=0.25$ dex at $4<z<6$ (a factor of 1.8) and 0.37 dex at $6<z<9$ (a factor of 2.3). This also reduces the star-formation rates (SFRs) by $\langle Δ\log\mathrm{SFR} \rangle=0.14$ dex at $4<z<6$ and 0.27 dex at $6<z<9$. The MIRI data also improve constraints on the allowable stellar mass formed in early star-formation. We model this using a star-formation history that includes both a "burst' at $z_f=100$ and a slowly varying ("delayed-$τ$") model. The MIRI data reduce the allowable stellar mass by 0.6 dex at $4<z< 6$ and by $\approx$1 dex at $6<z<9$. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early universe ($z\approx9$). Therefore, observations of rest-frame $\gtrsim$1 $μ$m are paramount for constraining the stellar-mass build-up in galaxies at very high-redshifts. △ Less

Submitted 25 March, 2023; v1 submitted 30 December, 2022; originally announced January 2023.

Comments: Updated with accepted ApJ version. Part of the CEERS Focus Issue. 27 pages, many figures (4 Figure Sets, available upon reasonable request)

The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques

Authors: Camilla Pacifici, Kartheik G. Iyer, Bahram Mobasher, Elisabete da Cunha, Viviana Acquaviva, Denis Burgarella, Gabriela Calistro Rivera, Adam C. Carnall, Yu-Yen Chang, Nima Chartab, Kevin C. Cooke, Ciaran Fairhurst, Jeyhan Kartaltepe, Joel Leja, Katarzyna Malek, Brett Salmon, Marianna Torelli, Alba Vidal-Garcia, Mederic Boquien, Gabriel G. Brammer, Michael J. I. Brown, Peter L. Capak, Jacopo Chevallard, Chiara Circosta, Darren Croton , et al. (30 additional authors not shown)

Abstract: The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have signifi… ▽ More The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 25 pages, 11 figures. Accepted for publication in ApJ

CEERS Key Paper I: An Early Look into the First 500 Myr of Galaxy Formation with JWST

Authors: Steven L. Finkelstein, Micaela B. Bagley, Henry C. Ferguson, Stephen M. Wilkins, Jeyhan S. Kartaltepe, Casey Papovich, L. Y. Aaron Yung, Pablo Arrabal Haro, Peter Behroozi, Mark Dickinson, Dale D. Kocevski, Anton M. Koekemoer, Rebecca L. Larson, Aurelien Le Bail, Alexa M. Morales, Pablo G. Perez-Gonzalez, Denis Burgarella, Romeel Dave, Michaela Hirschmann, Rachel S. Somerville, Stijn Wuyts, Volker Bromm, Caitlin M. Casey, Adriano Fontana, Seiji Fujimoto , et al. (42 additional authors not shown)

Abstract: 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 det… ▽ More 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. △ Less

Submitted 4 November, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: Replaced with published version

Identification of Large Equivalent Width Dusty Galaxies at 4 $<$ z $<$ 6 from Sub-mm Colours

Authors: Denis Burgarella, Patrice Theulé, Véronique Buat, Lisa Gouiran, Lorie Turco, Médéric Boquien, Tom J. L. C. Bakx, Akio K. Inoue, Yoshinobu Fudamoto, Yuma Sugahara, Jorge Zavala

Abstract: 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.… ▽ More 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. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: Paper accepted in Astronomy and Astrophysics on 10 November 2022

Journal ref: A&A 671, A123 (2023)

ALMA FIR View of Ultra High-redshift Galaxy Candidates at $z\sim$ 11-17: Blue Monsters or Low-$z$ Red Interlopers?

Authors: Seiji Fujimoto, Steven L. Finkelstein, Denis Burgarella, Chris L. Carilli, Véronique Buat, Caitlin M. Casey, Laure Ciesla, Sandro Tacchella, Jorge A. Zavala, Gabriel Brammer, Yoshinobu Fudamoto, Masami Ouchi, Francesco Valentino, M. C. Cooper, Mark Dickinson, Maximilien Franco, Mauro Giavalisco, Taylor A. Hutchison, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Takashi Kojima, Rebecca L. Larson, Eric J. Murphy, Casey Papovich, Pablo G. Pérez-González , et al. (28 additional authors not shown)

Abstract: 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~$μ$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr… ▽ More 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~$μ$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$σ$ 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$μ$m at $z=16.01$ or [CII]158$μ$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 ($Δ$ $χ^{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$β$) 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. △ Less

Submitted 26 July, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 25 pages, 8 figures, 7 tables. ApJ in press

Comparison of star formation histories of AGN and non-AGN galaxies

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, Q. Ni, E. Pouliasis, D. Burgarella, P. Theule, I. Georgantopoulos

Abstract: We use AGN with $\rm L_{X} \sim 10^{42.5-44}\,erg\,s^{-1}$, from the COSMOS-Legacy survey that lie within the UltraVISTA region and cross match them with the LEGA-C catalogue. The latter provides measurements of the calcium break, D$_n$4000, and H$_δ$ Balmer line that allow us to study the stellar populations of AGN and compare them with a galaxy reference catalogue. Our samples consist of 69 AGN… ▽ More We use AGN with $\rm L_{X} \sim 10^{42.5-44}\,erg\,s^{-1}$, from the COSMOS-Legacy survey that lie within the UltraVISTA region and cross match them with the LEGA-C catalogue. The latter provides measurements of the calcium break, D$_n$4000, and H$_δ$ Balmer line that allow us to study the stellar populations of AGN and compare them with a galaxy reference catalogue. Our samples consist of 69 AGN and 2176 non-AGN systems, within $\rm 0.6<z<1.3$, that satisfy the same photometric selection criteria. We construct the SEDs of both population and use the CIGALE code to investigate the effect of the two indices in the SED fitting process. Our analysis shows that the inclusion of D$_n$4000 and H$_δ$ allows CIGALE to constrain better the ages of the stellar populations. Furthermore, we find an increase of the estimated stellar masses by, on average, $\sim 0.2$ dex, in particular for systems with young stars (D$_n$4000$\,<1.5$), when the two indices are included in the SED fitting. We then compare the D$_n$4000 and H$_δ$ of AGN with sources in the reference catalogue, accounting for the different stellar mass of the two populations. Low to moderate L$_X$ AGN tend to reside in galaxies with older stellar populations and are less likely to have experienced a recent star formation burst, compared to galaxies in the control sample. We, also, compare the two populations as a function of their morphology (bulge-dominated, BD, vs. non-BD) and compactness. A similar fraction of AGN and non-AGN systems are classified as non-BD ($\sim 70\%$). Our analysis shows that BD AGN tend to have younger stellar populations compared to BD, non-AGN systems. On the other hand, non-BD AGN have, on average, older stellar populations and are less likely to have experienced a burst compared to non-BD sources in the reference sample. Furthermore, AGN tend to prefer more compact systems compared to non-AGN. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: A&A accepted. 15 pages, 18 figures

Journal ref: A&A 667, A145 (2022)

CEERS Key Paper III: The Resolved Host Properties of AGN at 3 < z < 5 with JWST

Authors: Dale D. Kocevski, Guillermo Barro, Elizabeth J. McGrath, Steven L. Finkelstein, Micaela B. Bagley, Henry C. Ferguson, Shardha Jogee, Guang Yang, Mark Dickinson, Nimish P. Hathi, Bren E. Backhaus, Eric F. Bell, Laura Bisigello, Véronique Buat, Denis Burgarella, Caitlin M. Casey, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Darren Croton, Emanuele Daddi, Adriano Fontana, Seiji Fujimoto, Jonathan P. Gardner, Eric Gawiser , et al. (34 additional authors not shown)

Abstract: 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 redshif… ▽ More 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. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Submitted to ApJ; 12 Pages, 6 Figures

Dusty Starbursts Masquerading as Ultra-high Redshift Galaxies in JWST CEERS Observations

Authors: Jorge A. Zavala, Veronique Buat, Caitlin M. Casey, Denis Burgarella, Steven L. Finkelstein, Micaela B. Bagley, Laure Ciesla, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Maximilien Franco, E. F. Jim'enez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Aurélien Le Bail, E. J. Murphy, Casey Papovich, Sandro Tacchella, Stephen M. Wilkins, Itziar Aretxaga, Peter Behroozi, Jaclyn B. Champagne, Adriano Fontana, Mauro Giavalisco, Andrea Grazian , et al. (99 additional authors not shown)

Abstract: 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 als… ▽ More 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. △ Less

Submitted 30 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Published in The Astrophysical Journal Letters (updated to match the published version)

Journal ref: ApJL 943 L9 (2023)

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: 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 f… ▽ More 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. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 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

The Physical Conditions of Emission-Line Galaxies at Cosmic Dawn from JWST/NIRSpec Spectroscopy in the SMACS 0723 Early Release Observations

Authors: Jonathan R. Trump, Pablo Arrabal Haro, Raymond C. Simons, Bren E. Backhaus, Ricardo O. Amorín, Mark Dickinson, Vital Fernández, Casey Papovich, David C. Nicholls, Lisa J. Kewley, Samantha W. Brunker, John J. Salzer, Stephen M. Wilkins, Omar Almaini, Micaela B. Bagley, Danielle A. Berg, Rachana Bhatawdekar, Laura Bisigello, Véronique Buat, Denis Burgarella, Antonello Calabrò, Caitlin M. Casey, Laure Ciesla, Nikko J. Cleri, Justin W. Cole , et al. (39 additional authors not shown)

Abstract: 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 wav… ▽ More 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]$λ$5008/H$β$ ratios, similar [OIII]$λ$4364/H$γ$ ratios, and higher ($\sim$0.5 dex) [NeIII]$λ$3870/[OII]$λ$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$λ$3870/[OII]$λ$3728, [OIII]$λ$4364/H$γ$, and [OIII]$λ$5008/H$β$ 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]$λ$4364/H$γ$ and [OIII]$λ$(4960+5008)/H$β$ 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. △ Less

Submitted 19 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

Spatial disconnection between stellar and dust emissions: the test of the Antennae Galaxies (Arp 244)

Authors: L. -M. Seillé, V. Buat, W. Haddad, A. Boselli, M. Boquien, L. Ciesla, Y. Roehlly, D. Burgarella

Abstract: The detection with of the Atacama Large Millimeter Array (ALMA) of dust-rich high redshift galaxies whose cold dust emission is spatially disconnected from the ultraviolet emission bears a challenge for modelling their spectral energy distributions (SED) with codes based on an energy budget between the stellar and dust components. We test the validity of energy balance modelling on a nearby resolv… ▽ More The detection with of the Atacama Large Millimeter Array (ALMA) of dust-rich high redshift galaxies whose cold dust emission is spatially disconnected from the ultraviolet emission bears a challenge for modelling their spectral energy distributions (SED) with codes based on an energy budget between the stellar and dust components. We test the validity of energy balance modelling on a nearby resolved galaxy with vastly different ultraviolet and infrared spatial distributions and infer what information can be reliably retrieved from the analysis of the full spectral energy distribution. We use 15 broadband images of the Antennae Galaxies ranging from far-ultraviolet to far-infrared and divide Arp 244 into 58 square ~1 kpc$^2$ regions. We fit the data with CIGALE to determine the star formation rate, stellar mass and dust attenuation of each region. We compare these quantities for the addition of the 58 regions to the ones obtained for Arp 244 as a whole and find that both estimates are consistent within one sigma. We present the spatial distribution of these physical parameters as well as the shape of the attenuation curve across the Antennae Galaxies . We also observe a flattening of the attenuation curves with increasing attenuation and dust surface density in agreement with the predictions of hydrodynamical simulations coupled with radiative transfer modelling. △ Less

Submitted 16 July, 2022; originally announced July 2022.

Comments: 13 pages, 9 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 665, A137 (2022)

Unsupervised classification of CIGALE galaxy spectra

Authors: J Dubois, D Fraix-Burnet, J Moultaka, P Sharma, D Burgarella

Abstract: Aims. The present study aims at providing a deeper insight into the power and limitation of an unsupervised classification algorithm (called Fisher-EM) on spectra of galaxies. This algorithm uses a Gaussian mixture in a discriminative latent subspace. To this end, we investigate the capacity of this algorithm to segregate the physical parameters used to generate mock spectra and the influence of t… ▽ More Aims. The present study aims at providing a deeper insight into the power and limitation of an unsupervised classification algorithm (called Fisher-EM) on spectra of galaxies. This algorithm uses a Gaussian mixture in a discriminative latent subspace. To this end, we investigate the capacity of this algorithm to segregate the physical parameters used to generate mock spectra and the influence of the noise on the classification. Methods. With the code CIGALE and different values for nine input parameters characterising the stellar population, we have simulated a sample of 11 475 optical spectra of galaxies containing 496 monochromatic fluxes. The statistical model and the optimum number of clusters is given in Fisher-EM by the integrated completed likelihood (ICL) criterion. We repeated the analyses several times to assess the robustness of the results. Results. Two distinct classifications can be distinguished in the case of the noiseless spectra. The one above 13 clusters disappears when noise is added, while the classification with 12 clusters is very robust against noise down to a signal to noise ratio (SNR) of 3. At SNR=1, the optimum is 5 clusters, but the classification is still compatible with the previous one. The distribution of the parameters used for the simulation shows an excellent discrimination between classes. A higher dispersion both in the spectra within each class and in the parameter distribution, leads us to conclude that despite a much higher ICL, the classification with more than 13 clusters in the noiseless case is not physically relevant. Conclusions. This study yields two conclusions valid at least for the Fisher-EM algorithm. Firstly, the unsupervised classification of spectra of galaxies is both reliable and robust to noise. Secondly, such analyses are able to extract the useful physical information contained in the spectra and to build highly meaningful classifications. In an epoch of data-driven astrophysics, it is important to trust unsupervised machine learning approaches that do not require training samples which are unavoidably biased. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: Astronomy and Astrophysics - A&A, EDP Sciences, In press

Journal ref: A&A 663, A21 (2022)

Tempered, Anti-trunctated, Multiple Importance Sampling

Authors: Grégoire Aufort, Pierre Pudlo, Denis Burgarella

Abstract: 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 Mu… ▽ More 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. △ Less

Submitted 16 June, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comparison of the star formation of X-ray selected active galactic nuclei in eFEDS with star-forming galaxies

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla, K. Malek, R. Shirley

Abstract: We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with… ▽ More We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with $\sim 17000$ sources. Both samples have the same photometric coverage, from optical to far-infrared. We construct the spectral energy distributions (SEDs) of all sources and fit them using the CIGALE code, utilizing the same templates and parametric grid for both samples. We account for the mass incompleteness and exclude quiescent systems from both datasets. These allow us to compare the SFR of the two populations in a uniform manner, minimising systematic effects. Based on our analysis, AGN at low and moderate L$_X$ ($\rm L_{X,2-10keV} < 10^{44}\,ergs^{-1}$), have SFR that is lower, or at most, equal to that of star-forming galaxies, in agreement with previous studies. The large number of luminous X-ray AGN available in the eFEDS dataset, enable us to expand our investigations at higher L$_X$ to test previous, tentative results. At $\rm L_{X,2-10keV} > 10^{44.2}\,ergs^{-1}$, the SFR of AGN appears enhanced, by $\sim 30\%$, compared to that of star-forming sources, for systems with stellar mass, $\rm 10.5 < log\,[M_*(M_\odot)] < 11.5$, confirming indications found in previous studies. The most massive sources $\rm log\,[M_*(M_\odot)] > 11.5$, present a flat SFR$_{norm}$-L$_X$ relation up to $\rm L_{X,2-10keV} \sim 10^{44.5}\,ergs^{-1}$, with SFR similar to that of star-forming galaxies. However, at higher L$_X$ ($\rm L_{X,2-10keV} \sim 10^{45}\,ergs^{-1}$), we find indications that the SFR of these massive AGN hosts may be enhanced compared to that of non-AGN systems. △ Less

Submitted 5 May, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 14 pages, 7 figures, 2 tables

Journal ref: A&A 663, A130 (2022)

IR characteristic emission and dust properties of star-forming galaxies at 4.5 $<$ z $<$ 6.2

Authors: D. Burgarella, J. Bogdanoska, A. Nanni, S. Bardelli, M. Bethermin, M. Boquien, V. Buat, A. L. Faisst, M. Dessauges-Zavadsky, Y. Fudamoto, S. Fujimoto, M. Giavalisco, M. Ginolfi, C. Gruppioni, N. P. Hathi, E. Ibar, G. C. Jones, A. M. Koekemoer, K. Kohno, B. C. Lemaux, D. Narayanan, P. Oesch, M. Ouchi, D. A. Riechers, F. Pozzi , et al. (9 additional authors not shown)

Abstract: The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amou… ▽ More The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amount of information in a collective way. From a subsample of 27 ALMA-detected galaxies at z > 4.5, we built an infrared spectral energy distribution composite template. It was used to fit, with CIGALE, the 105 galaxies (detections and upper limits) in the sample from the FUV to the FIR. The derived physical parameters provide information to decipher the nature of the dust cycle and of the stellar populations in these galaxies. The derived IR composite template is consistent with the galaxies in the studied sample. A delayed star formation history with tau_main = 500 Myrs is slightly favored by the statistical analysis as compared to a delayed with a final burst or a continuous star formation history. The position of the sample in the star formation rate (SFR)- M_star diagram is consistent with previous papers. The redshift evolution of the log M_star versus A_FUV relation is in agreement with evolution in the redshift of this relation. This evolution is necessary to explain the cosmic evolution of the average dust attenuation of galaxies. Evolution is also observed in the L_dust/ L_FUV (IRX) versus UV slope beta_FUV diagram: younger galaxies have bluer beta_FUV. We modeled the shift of galaxies in the IRX versus the beta_FUV diagram with the mass-weighted age as a free parameter, and we provide an equation to make predictions. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 40 pages, 13 figures

Journal ref: A&A 664, A73 (2022)

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

Authors: Mederic Boquien, Veronique Buat, Denis Burgarella, Sandro Bardelli, Matthieu Bethermin, Andreas Faisst, Michele Ginolfi, Nimish Hathi, Gareth Jones, Anton Koekemoer, Brian Lemaux, Desika Narayanan, Michael Romano, Daniel Schaerer, Daniela Vergani, Giovanni Zamorani, Elena Zucca

Abstract: 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… ▽ More 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. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 20 pages, 14 figures. Accepted for publication in A&A

Journal ref: A&A 663, A50 (2022)

A blast from the infant Universe: the very high-z GRB 210905A

Authors: A. Rossi, D. D. Frederiks, D. A. Kann, M. De Pasquale, E. Pian, G. Lamb, P. D'Avanzo, L. Izzo, A. J. Levan, D. B. Malesani, A. Melandri, A. Nicuesa Guelbenzu, S. Schulze, R. Strausbaugh, N. R. Tanvir, L. Amati, S. Campana, A. Cucchiara, G. Ghirlanda, M. Della Valle, S. Klose, R. Salvaterra, R. Starling, G. Stratta, A. E. Tsvetkova , et al. (30 additional authors not shown)

Abstract: We present a detailed follow-up of the very energetic GRB 210905A at a high redshift of z = 6.312 and its luminous X-ray and optical afterglow. We obtained a photometric and spectroscopic follow-up in the optical and near-infrared (NIR), covering both the prompt and afterglow emission from a few minutes up to 20 Ms after burst. With an isotropic gamma-ray energy release of Eiso = 1.27E54 erg, GRB… ▽ More We present a detailed follow-up of the very energetic GRB 210905A at a high redshift of z = 6.312 and its luminous X-ray and optical afterglow. We obtained a photometric and spectroscopic follow-up in the optical and near-infrared (NIR), covering both the prompt and afterglow emission from a few minutes up to 20 Ms after burst. With an isotropic gamma-ray energy release of Eiso = 1.27E54 erg, GRB 210905A lies in the top ~7% of gamma-ray bursts (GRBs) in terms of energy released. Its afterglow is among the most luminous ever observed. It starts with a shallow evolution that can be explained by energy injection, and it is followed by a steeper decay, while the spectral energy distribution is in agreement with slow cooling in a constant-density environment within the standard fireball theory. A jet break at ~ 46.2+-16.3 d (~6.3 d rest-frame) has been observed in the X-ray light curve; however, it is hidden in the H band due to the contribution from the likely host galaxy, the fourth GRB host at z > 6 known to date. We derived a half-opening angle of 8.4+-1.0 degrees, which is the highest ever measured for a z>6 burst, but within the range covered by closer events. The resulting collimation-corrected gamma-ray energy release of 1E52 erg is also among the highest ever measured. The moderately large half-opening angle argues against recent claims of an inverse dependence of the half-opening angle on the redshift. The total jet energy is likely too large to be sustained by a standard magnetar, and it suggests that the central engine of this burst was a newly formed black hole. Despite the outstanding energetics and luminosity of both GRB 210905A and its afterglow, we demonstrate that they are consistent with those of less distant bursts, indicating that the powering mechanisms and progenitors do not evolve significantly with redshift. △ Less

Submitted 4 August, 2022; v1 submitted 9 February, 2022; originally announced February 2022.

Comments: 19 pages, 12 figures, 5 tables, accepted in Astronomy & Astrophysics

Journal ref: A&A 665, A125 (2022)

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

Authors: Guang Yang, Médéric Boquien, William N. Brandt, Véronique Buat, Denis Burgarella, Laure Ciesla, Bret D. Lehmer, Katarzyna E. Małek, George Mountrichas, Casey Papovich, Estelle Pons, Marko Stalevski, Patrice Theulé, Shifu Zhu

Abstract: 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). Re… ▽ More 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 △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: ApJ accepted. The code is available at https://cigale.lam.fr/news/

Fitting spectral energy distributions of FMOS-COSMOS emission-line galaxies at z$\sim$1.6: Star formation rates, dust attenuation, and [OIII]$λ$5007 emission-line luminosities

Authors: J. A. Villa-Vélez, V. Buat, P. Theulé, M. Boquien, D. Burgarella

Abstract: We perform SED fitting analysis on a COSMOS sample covering UV-to-FIR wavelengths with emission lines from the FMOS survey. The sample of 182 objects with H$α$ and [OIII]$\lambda5007$ emission spans over a range of $1.40<\rm{z}<1.68$. We obtain robust estimates of stellar mass ($10^{9.5}-10^{11.5}~\rm{M_\odot}$) and SFR ($10^1-10^3~\rm{M_\odot}~\rm{yr}^{-1}$) from the Bayesian analysis with CIGALE… ▽ More We perform SED fitting analysis on a COSMOS sample covering UV-to-FIR wavelengths with emission lines from the FMOS survey. The sample of 182 objects with H$α$ and [OIII]$\lambda5007$ emission spans over a range of $1.40<\rm{z}<1.68$. We obtain robust estimates of stellar mass ($10^{9.5}-10^{11.5}~\rm{M_\odot}$) and SFR ($10^1-10^3~\rm{M_\odot}~\rm{yr}^{-1}$) from the Bayesian analysis with CIGALE fitting continuum photometry and H$α$. We obtain a median attenuation of A$_\rm{Hα}=1.16\pm0.19$ mag and A$_\rm{[OIII]}=1.41\pm0.22$ mag. H$α$ and [OIII]$\lambda5007$ attenuations are found to increase with stellar mass, confirming previous findings. A difference of $57$% in the attenuation experienced by emission lines and continuum is found in agreement with the lines being more attenuated than the continuum. New CLOUDY HII-region models in CIGALE enable good fits of H$α$, H$β$, [OIII]$\lambda5007$ emission lines with differences smaller than $0.2$ dex. Fitting [NII]$\lambda6584$ line is challenging due to well-known discrepancies in the locus of galaxies in the BPT diagram at intermediate redshifts. We find a positive correlation for SFR and dust-corrected L$_\rm{[OIII]\lambda5007}$ and we derive the linear relation $\log_{10}\rm{(SFR/\rm{M}_\odot~\rm{yr}^{-1})}=\log_{10} (\rm{L}_{[\rm{OIII]}}/\rm{ergs~s^{-1}})-(41.20\pm0.02)$. Leaving the slope as a free parameter leads to $\log_{10}\rm{(SFR/\rm{M}_\odot~\rm{yr}^{-1})}=(0.83\pm0.06)\log_{10}(\rm{L}_{[\rm{OIII]}}/\rm{ergs~s^{-1}})-(34.01\pm2.63)$. Gas-phase metallicity and ionization parameter variations account for a $0.24$ dex and $1.1$ dex of the dispersion, respectively. An average value of $\log\rm{U}\approx-2.85$ is measured for this sample. Including HII-region models to fit simultaneously photometry and emission line fluxes are paramount to analyze future data from surveys such as MOONS and PFS. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: Accepted for publication in A&A, 17 pages, 19 figures

Journal ref: A&A 654, A153 (2021)

Polar dust obscuration in broad-line active galaxies from the XMM-XXL field

Authors: V. Buat, G. Mountrichas, G. Yang, M. Boquien, Y. Roehlly, D. Burgarella, M. Stalevski, L. Ciesla, P. Theulé

Abstract: Dust is observed in the polar regions of nearby AGN and it is known to contribute substantially to their mid-IR emission and to the obscuration of their UV to optical emission. We selected a sample of 1275 BLAGN in the XMM-XXL field, with optical to infrared photometric data. These AGN are seen along their polar direction and we expect a maximal impact of dust located around the poles when it is p… ▽ More Dust is observed in the polar regions of nearby AGN and it is known to contribute substantially to their mid-IR emission and to the obscuration of their UV to optical emission. We selected a sample of 1275 BLAGN in the XMM-XXL field, with optical to infrared photometric data. These AGN are seen along their polar direction and we expect a maximal impact of dust located around the poles when it is present. We used X-CIGALE, which introduces a dust component to account for obscuration along the polar directions, modeled as a foreground screen, and an extinction curve that is chosen as it steepens significantly at short wavelengths or is much grayer. By comparing the results of different fits, we are able to define subsamples of sources with positive statistical evidence in favor of or against polar obscuration and described using the gray or steep extinction curve. We find a similar fraction of sources with positive evidence for and against polar dust. Applying statistical corrections, we estimate that half of our sample could contain polar dust and among them, 60% exhibit a steep extinction curve and 40% a flat extinction curve; although these latter percentages are found to depend on the adopted extinction curves. The obscuration in the V-band is not found to correlate with the X-ray column density, while A_V/N_H ratios span a large range of values and higher dust temperatures are found with the flat, rather than with the steep extinction curve. Ignoring this polar dust component in the fit of the spectral energy distribution of these composite systems leads to an overestimation of the stellar contribution. A single fit with a polar dust component described with an SMC extinction curve efficiently overcomes this issue but it fails at identifying all the AGN with polar dust obscuration. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: 12 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

The evolution of merger fraction of galaxies at z < 0.6 depending on the star formation mode in the AKARI NEP Wide field

Authors: Eunbin Kim, Ho Seong Hwang, Woong-Seob Jeong, Seong Jin Kim, Denis Burgarella, Tomotsugu Goto, Tetsuya Hashimoto, Young-Soo Jo, Jong Chul Lee, Matthew Malkan, Chris Pearson, Hyunjin Shim, Yoshiki Toba, Simon C. -C. Ho, Daryl Joe Santos, Hiroyuki Ikeda, Helen K. Kim, Takamitsu Miyaji, Hideo Matsuhara, Nagisa Oi, Toshinobu Takagi, Ting-Wen Wang

Abstract: We study the galaxy merger fraction and its dependence on star formation mode in the5.4 square degrees of the North Ecliptic Pole-Wide field. We select 6352 galaxies withAKARI 9μm detections, and identify mergers among them using the Gini coefficientand M20derived from the Subaru/HSC optical images. We obtain the total infraredluminosity and star formation rate of galaxies using the spectral energ… ▽ More We study the galaxy merger fraction and its dependence on star formation mode in the5.4 square degrees of the North Ecliptic Pole-Wide field. We select 6352 galaxies withAKARI 9μm detections, and identify mergers among them using the Gini coefficientand M20derived from the Subaru/HSC optical images. We obtain the total infraredluminosity and star formation rate of galaxies using the spectral energy distributiontemplates based on one band, AKARI 9μm. We classify galaxies into three differentstar formation modes (i.e. starbursts, main sequence, and quiescent galaxies) andcalculate the merger fractions for each. We find that the merger fractions of galaxiesincrease with redshift atz<0.6. The merger fractions of starbursts are higher thanthose of main sequence and quiescent galaxies in all redshift bins. We also examinethe merger fractions of far-infrared detected galaxies which have at least one detectionfromHerschel/SPIRE. We find thatHerscheldetected galaxies have higher mergerfraction compared to non-Herscheldetected galaxies, and bothHerscheldetected andnon-Herscheldetected galaxies show clearly different merger fractions depending onthe star formation modes. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: 13 pages, 11 figures, Accepted for publication by MNRAS

Getting ready for the LSST data -- estimating the physical properties of $z<2.5$ main sequence galaxies

Authors: Gabriele Riccio, Katarzyna Małek, Ambra Nanni, Mederic Boquien, Veronique Buat, Denis Burgarella, Darko Donevski, Mahmoud Hamed, Peter Hurley, Raphael Shirley, Agnieszka Pollo

Abstract: 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 sta… ▽ More 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<z<2.5$. In order to achieve this goal, we used the unique multi-wavelength data from the Herschel Extragalactic Legacy Project (HELP) survey. Our analysis focus on two fields: ELAIS-N1 and COSMOS. To obtain galaxy physical properties we fit their Spectral Energy Distributions (SEDs) using the Code Investigating GALaxy Emission (CIGALE). We compare the main galaxy physical properties obtained from the fit of the observed multi-wavelength photometry of galaxies (from UV to FIR) to the ones obtained from the simulated LSST optical measurements only. The stellar masses estimated based on the LSST measurements are in agreement with the full UV-FIR SED estimations, as they depend mainly on the UV and optical emission, well covered by LSST in the considered redshift range. We obtain a clear overestimation of SFR, $L_{dust}$, $M_{dust}$ estimated with LSST only, highly correlated with redshift. We investigate the cause of this overestimation and we conclude that it is related to an overestimation of the dust attenuation, both UV and NIR. We find that it is necessary to employ auxiliary rest-frame mid-infrared observations, simulated UV observations, or FUV attenuation (AFUV)- Mstar relation, to correct the overestimation. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: Accepted to A&A

Journal ref: A&A 653, A107 (2021)

Galaxy properties of type 1 and 2 X-ray selected AGN and comparison among different classification criteria

Authors: G. Mountrichas, V. Buat, I. Georgantopoulos, G. Yang, V. A. Masoura, M. Boquien, D. Burgarella

Abstract: We present analyses of host galaxy properties of type 1 and type 2 X-ray selected AGNs in the XMM-XXL field, which have available optical spectroscopic classification. We model their optical to far-infrared spectral energy distributions (SEDs) using the X-CIGALE code. X-CIGALE allows the fitting of X-ray flux and accounts for the viewing angle of dusty torus and the attenuation from polar dust. By… ▽ More We present analyses of host galaxy properties of type 1 and type 2 X-ray selected AGNs in the XMM-XXL field, which have available optical spectroscopic classification. We model their optical to far-infrared spectral energy distributions (SEDs) using the X-CIGALE code. X-CIGALE allows the fitting of X-ray flux and accounts for the viewing angle of dusty torus and the attenuation from polar dust. By selecting matched type 1 and 2 subsamples in the X-ray luminosity and redshift parameter space, we find that both types live in galaxies with similar star formation. However, type 2 AGN tend to reside in more massive systems ($10.87^{+0.06}_{-0.12}\,\rm M_\odot$) compared to their type 1 counterparts ($10.57^{+0.20}_{-0.12}\,\rm M_\odot$). In the second part of our analysis, we compare the spectroscopic classification with that from the SED fitting. X-CIGALE successfully identifies all spectroscopic type 2 sources either by estimating an inclination angle that corresponds to edge on viewing of the source or by measuring increased polar dust in these systems. $\sim 85\%$ of spectroscopic type 1 AGN are also identified as such, based on the SED fitting analysis. There is a small number of sources ($\sim 15\%$ of the sample), that present broad lines in their spectra, but show strong indications of obscuration, based on SED analysis. These, could be systems that are viewed face on and have an extended dust component along the polar direction. The performance of X-CIGALE in classifying AGN is similar at low and high redshifts, under the condition that there is sufficient photometric coverage. Finally, usage of optical/mid-IR colour criteria to identify optical red AGN (${\it{u}}-\rm W3$), suggests that these criteria are better suited for IR selected AGN and their efficiency drops for the low to moderate luminosity sources included in X-ray samples. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: A&A accepted

Journal ref: A&A 653, A70 (2021)

The role of AGN and its obscuration on the position of the host galaxy relative to Main Sequence

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla, K. Malek, R. Shirley

Abstract: We use X-ray Active Galactic Nuclei (AGN) observed by the Chandra X-ray Observatory within the 9.3 deg$^2$ Bo$\rm \ddot{o}$tes field of the NDWFS to study whether there is a correlation between X-ray luminosity (L$_X$) and star formation rate (SFR) of the host galaxy, at $\rm 0.5<z<2.0$, with respect to the position of the galaxy to the main sequence (SFR$_{norm}$). About half of the sources in th… ▽ More We use X-ray Active Galactic Nuclei (AGN) observed by the Chandra X-ray Observatory within the 9.3 deg$^2$ Bo$\rm \ddot{o}$tes field of the NDWFS to study whether there is a correlation between X-ray luminosity (L$_X$) and star formation rate (SFR) of the host galaxy, at $\rm 0.5<z<2.0$, with respect to the position of the galaxy to the main sequence (SFR$_{norm}$). About half of the sources in the X-ray sample have spectroscopic redshifts. We also construct a reference galaxy catalogue. For both datasets, we use photometric data from optical to the far infrared, compiled by the HELP project and apply spectral energy distribution (SED) fitting, using the X-CIGALE code. We exclude quiescent sources from both the X-ray and the reference samples. We also account for the mass completeness of our dataset, in different redshifts bins. Our analysis highlights the importance of studying the SFR-L$_X$ relation, in a uniform manner, taking into account the systematics and selection effects. Our results suggest that, in less massive galaxies ($\rm log\,[M_*(M_\odot)] \sim 11$), AGN enhances the SFR of the host galaxy by $\sim 50\%$ compared to non AGN systems. A flat relation is observed for the most massive galaxies. SFR$_{norm}$ does not evolve with redshift. The results, although tentative, are consistent with a scenario in which, in less massive systems, both AGN and star formation (SF) are fed by cold gas, supplied by a merger event. In more massive galaxies, the flat relation could be explained by a different SMBH fuelling mechanism that is decoupled from the star formation of the host galaxy (e.g. hot diffuse gas). Finally, we compare the host galaxy properties of X-ray absorbed and unabsorbed sources. Our results show no difference which suggests that X-ray absorption is not linked with the properties of the galaxy. △ Less

Submitted 20 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 653, A74 (2021)

HELP: The Herschel Extragalactic Legacy Project

Authors: R. Shirley, K. Duncan, M. C. Campos Varillas, P. D. Hurley, K. Malek, Y. Roehlly, M. W. L. Smith, H. Aussel, T. Bakx, V. Buat, D. Burgarella, N. Christopher, S. Duivenvoorden, S. Eales, A. Efstathiou, E. A. Gonzalez Solares, M. Griffin, M. Jarvis, B. Lo Faro, L. Marchetti, I. McCheyne, A. Papadopoulos, K. Penner, E. Pons, M. Prescott , et al. (7 additional authors not shown)

Abstract: We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the H… ▽ More We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or "open" methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenisation of calibration, meta data and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP database. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions. △ Less

Submitted 12 May, 2021; originally announced May 2021.

A hyperluminous obscured quasar at a redshift of z ~ 4.3

Authors: Andreas Efstathiou, Katarzyna Malek, Denis Burgarella, Peter Hurley, Seb Oliver, Veronique Buat, Raphael Shirley, Steven Duivenvoorden, Vicky Papadopoulou Lesta, Duncan Farrah, Kenneth J. Duncan, Maria del Carmen Campos Varillas

Abstract: 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 a… ▽ More 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\%$. △ Less

Submitted 8 April, 2022; v1 submitted 14 March, 2021; originally announced March 2021.

Comments: 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/

Journal ref: 2021MNRAS503L..11E

How does the Polar Dust affect the Correlation between Dust Covering Factor and Eddington Ratio in Type 1 Quasars Selected from the Sloan Digital Sky Survey Data Release 16?

Authors: Yoshiki Toba, Yoshihiro Ueda, Poshak Gandhi, Claudio Ricci, Denis Burgarella, Veronique Buat, Tohru Nagao, Shinki Oyabu, Hideo Matsuhara, Bau-Ching Hsieh

Abstract: 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… ▽ More 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 ($λ_{\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) $λ_{\rm Edd}$. As a result, anti-correlations for CF-$L_{\rm bol}$, CF-$M_{\rm BH}$, and CF-$λ_{\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. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 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

Origins Space Telescope: From First Light to Life -- ESA Voyage 2050 White Paper

Authors: M. C. Wiedner, S. Aalto, J. Birkby, D. Burgarella, P. Caselli, V. Charmandaris, A. Cooray, E. De Beck, J. -M. Desert, M. Gerin, J. Goicoechea, M. Griffin, P. Hartogh, F. Helmich, M. Hogerheijde, L. Hunt, A. Karska, Q. Krall, D. Leisawitz, G. Melnick, M. Meixner, M. Mikako, Ch. Pearson, D. Rigopoulou, T. Roellig , et al. (2 additional authors not shown)

Abstract: The Origins Space Telescope (Origins) is one of four science and technology definition studies selected by National Aeronautics and Space Administration (NASA) in preparation of the 2020 Astronomy and Astrophysics Decadal survey in the US. Origins will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. It is designed… ▽ More The Origins Space Telescope (Origins) is one of four science and technology definition studies selected by National Aeronautics and Space Administration (NASA) in preparation of the 2020 Astronomy and Astrophysics Decadal survey in the US. Origins will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. It is designed to answer three major science questions: How do galaxies form stars, make metals, and grow their central supermassive black holes from reionization? How do the conditions for habitability develop during the process of planet formation? Do planets orbiting M-dwarf stars support life? Origins operates at mid- to far-infrared wavelengths from ~2.8 to 588 μm, is more than 1000 times more sensitive than prior far-IR missions due to its cold (~4.5 K) aperture and state-of-the-art instruments. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 24 pages, 12 figures, 5 tables. Submitted as White paper to Voyage 2050 Long-term planning of the ESA Science Programme. arXiv admin note: substantial text overlap with arXiv:1912.06213

X-ray flux in the SED modelling: An application of X-CIGALE in the XMM-XXL field

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla

Abstract: X-CIGALE, built upon the spectral energy distribution (SED) code of CIGALE, implements important new features: The code accounts for obscuring material in the polars of the AGN and has the ability to fit X-ray fluxes. In this work, we use ~2500 spectroscopic, X-ray AGN from the XMM-XXL-North field and examine the improvements the new features bring in the SED modelling analysis. Based on our resul… ▽ More X-CIGALE, built upon the spectral energy distribution (SED) code of CIGALE, implements important new features: The code accounts for obscuring material in the polars of the AGN and has the ability to fit X-ray fluxes. In this work, we use ~2500 spectroscopic, X-ray AGN from the XMM-XXL-North field and examine the improvements the new features bring in the SED modelling analysis. Based on our results, X-CIGALE successfully connects the X-ray with the UV luminosity in the whole range spanned by our sample (log Lx(2-10 keV) = (42 - 46) erg/s). The addition of the new features globally improves the efficiency of X-CIGALE in the estimation and characterization of the AGN component. The classification into type 1 and type 2 based on their inclination angle is improved, especially at redshifts lower than 1. The statistical significance of the AGN fraction, fracAGN, measurements is increased, in particular for luminous X-ray sources (log LX > 45 erg/s). These conclusions hold under the condition that (mid-) IR photometry is available in the SED fitting process. The addition of polar dust increases the AGN fraction and the efficiency of the SED decomposition to detect AGN among X-ray selected sources. X-CIGALE estimates a strong AGN (fracAGN > 0.3) in more than 90% of the infrared selected AGN and 75% of X-ray detected AGN not selected by IR colour criteria. The latter drops to ~50% when polar dust is not included. The ability of X-CIGALE to include X-ray information in the SED fitting process can be instrumental in the optimal exploitation of the wealth of data that current (eROSITA) and future (ATHENA) missions will provide us. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: A&A accepted, 19 pages, 20 figures, 5 tables

Journal ref: A&A 646, A29 (2021)

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

Authors: Guang Yang, Casey Papovich, Micaela B. Bagley, Veronique Buat, Denis Burgarella, Mark E. Dickinson, David Elbaz, Steven L. Finkelstein, Adriano Fontana, Norman A. Grogin, Intae Jung, Jeyhan S. Kartaltepe, Allison Kirkpatrick, Pablo G. Pérez-González, Norbert Pirzkal, L. Y. Aaron Yung

Abstract: 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 poi… ▽ More 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. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 31 pages, 24 figures, 3 tables. Re-submitted to ApJS after revision

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

Authors: Dominik A. Riechers, Hooshang Nayyeri, Denis Burgarella, Bjorn H. C. Emonts, David L. Clements, Asantha Cooray, Rob J. Ivison, Seb Oliver, Ismael Perez-Fournon, Dimitra Rigopoulou, Douglas Scott

Abstract: 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^1… ▽ More 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. △ Less

Submitted 30 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 20 pages, 11 figures, 6 tables; to appear in the Astrophysical Journal (accepted November 29, 2020)

Extinction-free Census of AGNs in the $AKARI$/IRC North Ecliptic Pole Field from 23-band Infrared Photometry from Space Telescopes

Authors: Ting-Wen Wang, Tomotsugu Goto, Seong Jin Kim, Tetsuya Hashimoto, Denis Burgarella, Yoshiki Toba, Hyunjin Shim, Takamitsu Miyaji, Ho Seong Hwang, Woong-Seob Jeong, Eunbin Kim, Hiroyuki Ikeda, Chris Pearson, Matthew Malkan, Nagisa Oi, Daryl Joe D. Santos, Katarzyna Małek, Agnieszka Pollo, Simon C. -C. Ho, Hideo Matsuhara, Alvina Y. L. On, Helen K. Kim, Tiger Yu-Yang Hsiao, Ting-Chi Huang

Abstract: In order to understand the interaction between the central black hole and the whole galaxy or their co-evolution history along with cosmic time, a complete census of active galactic nuclei (AGN) is crucial. However, AGNs are often missed in optical, UV and soft X-ray observations since they could be obscured by gas and dust. A mid-infrared (mid-IR) survey supported by multiwavelength data is one o… ▽ More In order to understand the interaction between the central black hole and the whole galaxy or their co-evolution history along with cosmic time, a complete census of active galactic nuclei (AGN) is crucial. However, AGNs are often missed in optical, UV and soft X-ray observations since they could be obscured by gas and dust. A mid-infrared (mid-IR) survey supported by multiwavelength data is one of the best ways to find obscured AGN activities because it suffers less from extinction. Previous large IR photometric surveys, e.g., $WISE$ and $Spitzer$, have gaps between the mid-IR filters. Therefore, star forming galaxy (SFG)-AGN diagnostics in the mid-IR were limited. The $AKARI$ satellite has a unique continuous 9-band filter coverage in the near to mid-IR wavelengths. In this work, we take advantage of the state-of-the-art spectral energy distribution (SED) modelling software, CIGALE, to find AGNs in mid-IR. We found 126 AGNs in the NEP-Wide field with this method. We also investigate the energy released from the AGN as a fraction of the total IR luminosity of a galaxy. We found that the AGN contribution is larger at higher redshifts for a given IR luminosity. With the upcoming deep IR surveys, e.g., $JWST$, we expect to find more AGNs with our method. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Comments: 14 pages, 18 figures. Accepted for publication in MNRAS. For associated video, please see https://www.youtube.com/watch?v=B01jL0Bol9Q&feature=emb_logo