AGN Contribution Research Articles

A GMRT 610 MHz radio survey of the North Ecliptic Pole (NEP, ADF-N) / Euclid Deep Field North

ABSTRACT This paper presents a 610 MHz radio survey covering 1.94 deg$^2$ around the North Ecliptic Pole (NEP), which includes parts of the ${\it AKARI}$ (ADF-N) and Euclid, Deep Fields North. The median 5$\sigma$ sensitivity is 28 $\mu$Jy beam$^{-1}$, reaching as low as 19 $\mu$Jy beam$^{-1}$, with a synthesized beam of 3.6 arcsec $\times$ 4.1 arcsec. The catalogue contains 1675 radio components, with 339 grouped into multicomponent sources and 284 ‘isolated’ components likely part of double radio sources. Imaging, cataloguing, and source identification are presented, along with preliminary scientific results. From a non-statistical sub-set of 169 objects with multiwavelength AKARI and other detections, luminous infrared galaxies (LIRGs) represent 66 ${{\ \rm per\ cent}}$ of the sample, ultra-luminous infrared galaxies (ULIRGs) 4 ${{\ \rm per\ cent}}$, and sources with L$_{IR}$ < 10$^{11}$ $\mathrm{ L}_{\odot }\,$ 30${{\ \rm per\ cent}}$. In total, 56 ${{\ \rm per\ cent}}$ of sources show some AGN presence, though only seven are AGN-dominated. ULIRGs require three times higher AGN contribution to produce high-quality SED fits compared to lower luminosity galaxies, and AGN presence increases with AGN fraction. The PAH mass fraction is not significant, although ULIRGs have about half the PAH strength of lower IR-luminosity galaxies. Higher luminosity galaxies show gas and stellar masses an order of magnitude larger, suggesting higher star formation rates. For LIRGs, AGN presence increases with redshift, indicating that part of the total luminosity could be contributed by AGN activity rather than star formation. Simple cross-matching revealed 13 ROSAT QSOs, 45 X-ray sources, and 61 sub-mm galaxies coincident with GMRT radio sources.

Evolution of gas flows along the starburst to post-starburst to quiescent galaxy sequence

ABSTRACT We measure velocity offsets in the $\rm{Na {}\rm {\small I}}$ λλ5890, 5896 ($\rm{Na {}\rm {\small D}}$) interstellar medium absorption lines to track how neutral galactic winds change as their host galaxies evolve. Our sample of ∼80 000 SDSS spectra at 0.010 < z < 0.325 includes starburst, post-starburst, and quiescent galaxies, forming an evolutionary sequence of declining star formation rate (SFR). We detect bulk flows across this sequence, mostly at higher host stellar masses (log(M*/M⊙)) > 10). Along this sequence, the fraction of outflows decreases (76 ± 2 per cent to 65 ± 4 per cent to a 3σ upper limit of 34 per cent), and the mean velocity offset changes from outflowing to inflowing (−84.6 ± 5.9 to −71.6 ± 11.4 to $76.6\pm 2.3\, \rm km s^{-1}$). Even within the post-starburst sample, wind speed decreases with time elapsed since the starburst ended. These results reveal that outflows diminish as galaxies age. For post-starbursts, there is evidence for an AGN contribution, especially to the speediest outflows: (1) SFR declines faster in time than outflow velocity, a decoupling arguing against massive stellar feedback; (2) of the few outflows strong enough to escape the interstellar medium (9/105), three of the four hosts with measured emission lines are Seyfert galaxies. For discy starburst galaxies, however, the trends suggest flows out of the stellar disc plane (with outflow 1/2-opening angle > 45°) instead of from the nucleus: the wind velocity decreases as the disc becomes more edge-on, and the outflow fraction, constant at $\sim 90~{{\ \rm per\ cent}}$ for disc inclinations i < 45°, steadily decreases from $\sim 90~{{\ \rm per\ cent}}$ to 20 per cent for i > 45°.

Correction to: Deep Herschel observations of the 2 Jy sample: assessing the non-thermal and AGN contributions to the far-IR continuum

Correction to: Deep Herschel observations of the 2 Jy sample: assessing the non-thermal and AGN contributions to the far-IR continuum

Testing the Radio-selection Method of Dual Active Galactic Nuclei in the Stripe 82 Field

We test the merger-induced dual active galactic nuclei (dAGNs) paradigm using a sample of 35 radio galaxy pairs from the Sloan Digital Sky Survey Stripe 82 field. Using Keck optical spectroscopy, we confirm 21 pairs have consistent redshifts, constituting kinematic pairs; the remaining 14 pairs are line-of-sight projections. We classify the optical spectral signatures via emission line ratios, equivalent widths, and excess of radio power above star formation predicted outputs. We find six galaxies are classified as LINERs and seven are AGN/starburst composites. Most of the LINERs are retired galaxies, while the composites likely have AGN contribution. All of the kinematic pairs exhibit radio power more than 10× above the level expected from just star formation, suggestive of a radio AGN contribution. We also analyze high-resolution (0.″3) imaging at 6 GHz from the NSF’s Karl G. Jansky Very Large Array for 17 of the kinematic pairs. We find six pairs (two new, four previously known) host two separate radio cores, confirming their status as dAGNs. The remaining 11 pairs contain single AGNs, with most exhibiting prominent jets/lobes overlapping their companion. Our final census indicates a dAGN duty cycle slightly higher than predictions of purely stochastic fueling, although a larger sample (potentially culled from VLASS) is needed to fully address the dAGN fraction. We conclude that while dAGNs in the Stripe 82 field are rare, the merger process plays some role in their triggering and it facilitates low to moderate levels of accretion.

Enhanced Star Formation Efficiency in the Central Regions of Nearby Quasar Hosts

We combine Atacama Large Millimeter/submillimeter Array and Multi Unit Spectroscopic Explorer observations tracing the molecular gas, millimeter continuum, and ionized gas emission in six low-redshift (z ≲ 0.06) Palomar–Green (PG) quasar host galaxies to investigate their ongoing star formation at roughly kiloparsec-scale resolution. The AGN contribution to the cold dust emission and the optical emission-line flux is carefully removed to derive spatial distributions of the star formation rate (SFR), which, complemented with the molecular gas data, enables the mapping of the depletion time (t dep). We report ubiquitous star formation activity within the quasar host galaxies, with the majority of the ongoing star formation occurring in the galactic center. The rise of the SFR surface density (ΣSFR) toward the nucleus is steeper than that observed for the cold molecular gas surface density, reaching values up to ΣSFR ≈ 0.15–0.80 M ⊙ yr−1 kpc−2. The gas in the nuclear regions is converted into stars at a shortened depletion time (t dep ≈ 0.2–2.0 Gyr), suggesting that those zones can be deemed as starbursts. At large galactocentric radius, we find that the ongoing star formation takes place within spiral arms or H ii region complexes, with an efficiency comparable to that reported for nearby inactive spirals (t dep ≈ 1.8 Gyr). We find no evidence of star formation activity shutoff in the PG quasar host galaxies. On the contrary, these observations shed light on how the central environments of galaxies hosting actively accreting supermassive black holes build up stellar mass.

Mock HUBS observations of hot gas with IllustrisTNG

The lack of adequate X-ray observing capability is seriously impeding the progress in understanding the hot phase of circumgalactic medium (CGM), which is predicted to extend to the virial radius of a galaxy or beyond, and thus in acquiring key boundary conditions for studying galaxy evolution. To this end, the Hot Universe Baryon Surveyor (HUBS) is proposed. HUBS is designed to probe hot CGM by detecting its emission or absorption lines with a non-dispersive X-ray spectrometer of high resolution and high throughput. The spectrometer consists of a 60 × 60 array of microcalorimeters, with each detector providing an energy resolution of 2 eV, and is placed in the focal plane of an X-ray telescope of 1∘ field-of-view. With such a design, the spectrometer is highly optimized for detecting X-ray-emitting hot gas in the CGM of local galaxies, as well as in intra-group medium (IGrM), intra-cluster medium (ICM), or intergalactic medium (IGM). To assess the scientific potential of HUBS, in this work, we created mock observations of galaxies, groups, and clusters at different redshifts with the IllustrisTNG simulation. Focusing exclusively on emission studies, we took into account the effects of light cone, Galactic foreground emission, and background AGN contribution in the mock observations. From the observations, we made mock X-ray images and spectra, analyzed them to derive the properties of the emitting gas in each case, and compared the results with the input parameters from the simulation. The results show that HUBS is well suited for studying hot CGM at low redshifts. The redshift range is significantly extended for measuring IGrM and ICM. The sensitivity limits are also presented for detecting extended emission of low surface brightness.

Galaxy morphologies revealed with Subaru HSC and super-resolution techniques. I. Major merger fractions of LUV ∼ 3–15 L*UV dropout galaxies at z ∼ 4–7

Abstract We perform a super-resolution analysis of the Subaru Hyper Suprime-Cam (HSC) images to estimate the major merger fractions of z ∼ 4–7 dropout galaxies at the bright end of galaxy UV luminosity functions (LFs). Our super-resolution technique improves the spatial resolution of the ground-based HSC images, from ∼1″ to $\lesssim \!\!{0{^{\prime \prime }_{.}}1}$, which is comparable to that of the Hubble Space Telescope, allowing us to identify z ∼ 4–7 bright major mergers at a high completeness value of $\gtrsim \!\!90\%$. We apply the super-resolution technique to 6412, 16, 94, and 13 very bright dropout galaxies at z ∼ 4, 5, 6, and 7, respectively, in a UV luminosity range of LUV ∼ 3–$15\, L_{\rm UV}^*$ corresponding to −24 ≲ MUV ≲ −22. The major merger fractions are estimated to be $f_{\rm merger}\sim 10\%$–$20\%$ at z ∼ 4 and $\sim 50\%$–$70\%$ at z ∼ 5–7, which shows no fmerger difference compared to those of a control faint galaxy sample. Based on the fmerger estimates, we verify contributions of source blending effects and major mergers to the bright-end of double power-law (DPL) shape of z ∼ 4–7 galaxy UV LFs. While these two effects partly explain the DPL shape at LUV ∼ 3–$10\, L_{\rm UV}^*$, the DPL shape cannot be explained at the very bright end of $L_{\rm UV}\gtrsim 10\, L_{\rm UV}^*$, even after the AGN contribution is subtracted. The results support scenarios in which other additional mechanisms such as insignificant mass quenching and low dust obscuration contribute to the DPL shape of galaxy UV LFs.

Multi-Wavelength Study of a Proto-BCG at z = 1.7

In this work we performed a spectral energy distribution (SED) analysis in the optical/infrared band of the host galaxy of a proto-brightest bluster galaxy (BCG, NVSS J103023 + 052426) in a proto-cluster at z = 1.7. We found that it features a vigorous star formation rate (SFR) of ∼570 M⊙/yr and a stellar mass of M*∼3.7×1011M⊙; the high corresponding specific SFR = 1.5±0.5Gyr−1 classifies this object as a starburst galaxy that will deplete its molecular gas reservoir in ∼3.5×108 yr. Thus, this system represents a rare example of a proto-BCG caught during the short phase of its major stellar mass assembly. Moreover, we investigated the nature of the host galaxy emission at 3.3 mm. We found that it originates from the cold dust in the interstellar medium, even though a minor non-thermal AGN contribution cannot be completely ruled out. Finally, we studied the polarized emission of the lobes at 1.4 GHz. We unveiled a patchy structure where the polarization fraction increases in the regions in which the total intensity shows a bending morphology; in addition, the magnetic field orientation follows the direction of the bendings. We interpret these features as possible indications of an interaction with the intracluster medium. This strengthens the hypothesis of positive AGN feedback, as inferred in previous studies of this object on the basis of X-ray/mm/radio analysis. In this scenario, the proto-BCG heats the surrounding medium and possibly enhances the SFR in nearby galaxies.

The Detection of Ionized Carbon Emission at z ∼ 8*

Abstract We present deep Keck/MOSFIRE H-band spectroscopic observations covering the [C iii],C iii]λλ1907, 1909 doublet for three z ∼ 8 galaxy candidates in the AEGIS field. Along with nondetections in two galaxies, we obtain one of the highest-redshift detections to date of [C iii]λ1907 for the galaxy AEGIS-33376, finding z spec = 7.945 ± 0.001. We measure a [C iii]λ1907 flux of 2.24 ± 0.71 × 10−18erg s−1cm−2, corresponding to a rest-frame equivalent width of 20.3 ± 6.5 Å for the single line. Given the not very constraining upper limit for C iii]λ1909 based on strong sky-line contamination, we assume a [C iii]λ1907/C iii]λ1909 doublet ratio of 1.5 and infer a total [C iii][C iii]λλ1907, 1909 equivalent width of 33.7 ± 10.8 Å. We repeat the same reductions and analysis on multiple subsets of our raw data divided on the basis of time and observing conditions, verifying that the [C iii]λ1907 emission is present for AEGIS-33376 throughout our observations. We also confirm that the significance of the [C iii]λ1907 detection in different subsets of our data tracks that of brighter emission features detected on the same multi-slit mask. These multiple tests suggest that the observed emission line is real and associated with the z ∼ 8 target. The strong observed [C iii],C iii]λλ1907, 1909 in AEGIS-33376 likely indicates ISM conditions of low metallicity, high ionization parameter, and a hard ionizing spectrum, although AGN contributions are possible. This single detection represents a sizable increase in the current sample [C iii],C iii]λλ1907, 1909 detections at z > 7, while JWST will provide the first statistical samples of such measurements at these redshifts.

Deep Observations of CO and Free–Free Emission in Ultraluminous Infrared QSO IRAS F07599+6508

Abstract Infrared quasi-stellar objects (IR QSOs) are a rare subpopulation selected from ultraluminous infrared galaxies (ULIRGs) and have been regarded as promising candidates of ULIRG-to-optical QSO transition objects. Here we present NOEMA observations of the CO (1−0) line and 3 mm continuum emission in IR QSO IRAS F07599+6508 at z = 0.1486, which has many properties in common with Mrk 231. The CO emission is found to be resolved with a major axis of ∼6.1 kpc that is larger than the size of ∼4.0 kpc derived for 3 mm continuum. We identify two faint CO features located at a projected distance of ∼11.4 and 19.1 kpc from the galaxy nucleus, both of which are found to have counterparts in the optical and radio bands and may have a merger origin. A systematic velocity gradient is found in the CO main component, suggesting that the bulk of molecular gas is likely rotationally supported. Based on the radio-to-millimeter spectral energy distribution and IR data, we estimate that about 30% of the flux at 3 mm arises from free–free emission and infer a free–free-derived star formation rate of 77 M ⊙ yr−1, close to the IR estimate corrected for the AGN contribution. We find a high-velocity CO emission feature at the velocity range of about −1300 to −2000 km s−1. Additional deep CO observations are needed to confirm the presence of a possible very high velocity CO extension of the OH outflow in this IR QSO.

The XXL Survey

We studied a sample of 274 radio and X-ray selected quasars (XQSOs) detected in the COSMOS and XXL-S radio surveys at 3 GHz and 2.1 GHz, respectively. This sample was identified by adopting a conservative threshold in X-ray luminosity, LX [2−10 keV] ≥ 1044 erg s−1, selecting only the most powerful quasars. A number of previous studies on the origin of radio emission in type-1 quasars have focused on the radio loudness distributions, some claiming to have found evidence for bimodality, pointing toward the existence of two physically different mechanisms for the radio emission. Using available multiwavelength data, we examined various criteria for the selection of radio-loud (RL) and radio-quiet (RQ) XQSOs and found that the number of RL/RQ XQSOs changes significantly depending on the chosen criterion. This discrepancy arises due to the different criteria tracing different physical processes and due to the fact that our sample was selected from flux-limited radio and X-ray surveys. Another approach to study the origin of radio emission in XQSOs is via their radio luminosity functions (RLF). We constructed the XQSO 1.4 GHz RLFs in six redshift bins at 0.5 ≤ z ≤ 3.75. The lower-1.4 GHz luminosity end shows a higher normalization than expected only from AGN contribution in all studied redshift bins. We found that the so-called “bump” is mostly dominated by emission due to star-forming processes within the host galaxies of XQSOs. As expected, AGN-related radio emission is the dominant contribution at the higher-luminosity end of RLF. To study the evolution of the XQSO RLF, we used a combination of analytic forms from the literature to constrain the “bump” due to star formation and the higher-luminosity AGN part of the RLF. We defined two 1.4 GHz luminosity thresholds, Lth, SF and Lth, AGN, below and above which more than 80% of sources contributing to the RLF are dominated by star formation and AGN-related activity, respectively. The two thresholds evolve with redshift, which is most likely driven by the strong evolution of star formation rates of the XQSO host galaxies. We found that both the lower and higher luminosity ends evolve significantly in density, while their luminosity evolution parameters are consistent with being constant. We found that the lower-luminosity end evolves both in density and luminosity, while the higher-luminosity end evolves significantly only in density. Our results expose the dichotomy of the origin of radio emission: while the higher-luminosity end of the XQSO RLF is dominated by AGN activity, the lower-luminosity end is dominated by the star formation-related processes.

AGN jets versus accretion as reionization sources

Context. Cosmic reionization put an end to the dark ages that came after the recombination era. Observations seem to favor the scenario of massive-star photons generated in low-mass galaxies being responsible for the bulk of reionization and, whereas a possible contribution from AGN accretion disks has been widely considered, they are currently thought to have had a minor role in reionization. Aims. We aim to study the possibility of AGN having contributed to reionization not only through their accretion disks, but also through ionizing photons coming from the AGN jets interacting with the intergalactic medium. Methods. We adopt an empirically derived AGN luminosity function at z ≃ 6, use X-ray observations to correct it for the presence of obscured sources, and estimate the density of jetted AGN. We then use analytical calculations to derive the fraction of jet energy that goes into ionizing photons. Finally, we compute the contribution of AGN jets to the H II volume filling factor at redshifts z ≃ 15−5. Results. We show that the contribution of the AGN jet lobes to the reionization of the Universe at z ∼ 6 might have been as high as ≳10% of that of star-forming galaxies, under the most favorable conditions of jetted and obscuration fraction. Conclusions. The contribution of AGN to the reionization, while most likely not dominant, could have been higher than previously assumed, thanks to the radiation originated in the jet lobes.

The OTELO survey

Aims. We take advantage of the capability of the OTELO survey to obtain the Hα luminosity function (LF) at z ∼ 0.40. Because of the deepest coverage of OTELO, we are able to determine the faint end of the LF, and thus better constrain the star formation rate and the number of galaxies at low luminosities. The AGN contribution to this LF is estimated as well. Methods. We make use of the multiwavelength catalogue of objects in the field compiled by the OTELO survey, which is unique in terms of minimum flux and equivalent width. We also take advantage of the pseudo-spectra built for each source, which allow the identification of emission lines and the discrimination of different types of objects. Results. The Hα luminosity function at z ∼ 0.40 is obtained, which extends the current faint end by almost 1 dex, reaching minimal luminosities of log10Llim = 38.5 erg s−1 (or ∼0.002 M⊙ yr−1). The AGN contribution to the total Hα luminosity is estimated. We find that no AGN should be expected below a luminosity of log10L = 38.6 erg s−1. From the sample of non-AGN (presumably, pure SFG) at z ∼ 0.40 we estimated a star formation rate density of ρSFR = 0.012 ± 0.005 M⊙ yr−1 Mpc−3.

Study of AGN contribution on morphological parameters of their host galaxies

AbstractWe tested how the AGN contribution (5%–75% of the total flux) may affect different morphological parameters commonly used in galaxy classification. We carried out all analysis at z ∼ 0 and at higher redshifts that correspond to the COSMOS field. Using a local training sample of >2000 visually classified galaxies, we carried out all measurements with and without the central source, and quantified how the contribution of a bright nuclear point source could affect different morphological parameters, such as: Abraham and Concelice-Bershady indices, Gini, Asymmetry, M20 moment of light, and Smoothness. We found that concentration indexes are less sensitive to both redshift and brightness in comparison to the other parameters. We also found that all parameters change significantly with AGN contribution. At z ∼ 0, up to a 10% of AGN contribution the morphological classification will not be significantly affect, but for ⩾ 25% of AGN contribution late-type spirals follow the range of parameters of elliptical galaxies and can therefore be misclassified early types.

Galaxy evolution studies in clusters: the case of Cl0024 + 1652 cluster galaxies at z ∼ 0.4

AbstractStudying the transformation of cluster galaxies contributes a lot to have a clear picture of evolution of the universe. Towards that we are studying different properties (morphology, star formation, AGN contribution and metallicity) of galaxies in clusters up to z ∼ 1.0 taking three different clusters: ZwCl0024 + 1652 at z ∼ 0.4, RXJ1257 + 4738 at z ∼ 0.9 and Virgo at z ∼ 0.0038. For ZwCl0024 + 1652 and RXJ1257 + 4738 clusters we used tunable filters data from GLACE survey taken with GTC 10.4 m telescope and other public data, while for Virgo we used public data. We did the morphological classification of 180 galaxies in ZwCl0024 + 1652 using galSVM, where 54 % and 46 % of galaxies were classified as early-type (ET) and late-type (LT) respectively. We did a comparison between the three clusters within the clustercentric distance of 1 Mpc and found that ET proportion (decreasing with redshift) dominates over the LT (increasing with redshift) throughout. We finalized the data reduction for ZwCl0024 + 1652 cluster and identified 46 [OIII] and 73 Hβ emission lines. For this cluster we have classified 22 emission line galaxies (ELGs) using BPT-NII diagnostic diagram resulting with 14 composite, 1 AGN and 7 star forming (SF) galaxies. We are using these results, together with the public data, for further analysis of the variations of properties in relation to redshift within z < 1.0.

No signs of star formation being regulated in the most luminous quasars at z ∼ 2 with ALMA

ABSTRACT We present ALMA Band 7 observations at $850\, \mu$m of 20 luminous ($\log \, L_{\rm bol}\,\gt\, 46.9$ [erg s−1]) unobscured quasars at z ∼ 2. We detect continuum emission for 19/20 quasars. After subtracting an AGN contribution, we measure the total far-IR luminosity for 18 quasars, assuming a modified blackbody model, and attribute the emission as indicative of the star formation rate (SFR). Our sample can be characterized with a lognormal SFR distribution having a mean of 140 M⊙ yr−1 and a dispersion of 0.5 dex. Based on an inference of their stellar masses, the SFRs are similar, in both the mean and dispersion, with star-forming main-sequence galaxies at the equivalent epoch. Thus, there is no evidence for a systematic enhancement or suppression (i.e. regulation or quenching) of star formation in the hosts of the most luminous quasars at z ∼ 2. These results are consistent with the Magneticum cosmological simulation, while in disagreement with a widely recognized phenomenological model that predicts higher SFRs than observed here based on the high bolometric luminosities of this sample. Furthermore, there is only a weak relation between SFR and accretion rate on to their supermassive black holes both for average and individual measurements. We interpret these results as indicative of star formation and quasar accretion being fed from the available gas reservoir(s) in their host with a disconnect due to their different physical sizes, temporal scales, and means of gas processing.

The Composite Nature of Dust-obscured Galaxies (DOGs) at z ∼ 2–3 in the COSMOS Field. II. The AGN Fraction

Abstract We present the X-ray properties of 108 Dust-Obscured Galaxies (DOGs; F 24 μm/F R > 1000) in the COSMOS field, all of which are detected in at least three far-infrared bands with the Herschel Observatory. Out of the entire sample, 22 are individually detected in the hard 2–8 keV X-ray band by the Chandra COSMOS Legacy survey, allowing us to classify them as AGN. Six (27%) of them are Compton-thick AGN candidates with column densities N H > 1024 cm−2, while 15 are moderately obscured AGNs with 1022 < N H < 1024 cm−2. Additionally, we estimate AGN contributions to the IR luminosity (8–1000 μm rest-frame) greater than 20% for 19 DOGs based on SED decomposition using Spitzer/MIPS 24 μm and the five Herschel bands (100–500 μm). Only 7 of these are detected in X-rays individually. We performed an X-ray stacking analysis for the 86 undetected DOGs. We find that the AGN fraction in DOGs increases with 24 μm flux and that it is higher than that of the general 24 μm population. However, no significant difference is found when considering only X-ray detections. This strongly motivates the combined use of X-ray and far-IR surveys to successfully probe a wider population of AGNs, particularly for the most obscured ones.

Galaxy properties derived with spectral energy distribution fitting in the Hawaii-Hubble Deep Field-North

We compile multi-wavelength data from ultraviolet to infrared (IR) bands as well as redshift and source-type information, for a large sample of 178 341 sources in the Hawaii-Hubble Deep Field-North field. A total of 145 635 sources among the full sample are classified/treated as galaxies and have redshift information available. We derive physical properties for these sources utilizing the spectral energy distribution fitting code CIGALE that is based on Bayesian analysis. Through various consistency and robustness checks, we find that our stellar-mass and star-formation rate (SFR) estimates are reliable, which is mainly due to two facts. Firstly, we adopt the most up-to-date and accurate redshifts and point spread function-matched photometry; and secondly, we make sensible parameter choices with the CIGALE code and take into account the influences of mid-IR/far-IR data, star-formation history models, and AGN contribution. We release our catalog of galaxy properties publicly (including, e.g., redshift, stellar mass, SFR, age, metallicity, dust attenuation). It is the largest of its kind in this field and should facilitate future relevant studies on the formation and evolution of galaxies.

Consistent modelling of the meta-galactic UV background and the thermal/ionization history of the intergalactic medium

Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment of the intergalactic medium (IGM) opacity for ionizing photons that is able to consistently capture the transition from a neutral to an ionized IGM. Previous synthesis models of the UVB, when used in simulations, yield reionization and thermal histories that are inconsistent with the assumed ionizing emissivities. With our new treatment, this discrepancy is fully resolved. In our fiducial model, galaxies leaking $\lesssim 18\%$ of their Lyman continuum emission drive HI reionization, while AGN drive HeII reionization (completing at $z \simeq 6.2$ and $2.8$ respectively). Due to the limited time available for cooling between HI and HeII reionization, higher IGM temperatures are predicted for late reionization scenarios. In our fiducial model, the predicted temperatures agree well with observational constraints at $z \lesssim 4$, while being slightly high compared to (somewhat uncertain) data above that. Models with a larger contribution of AGN are instead disfavoured by the temperature data, as well as by measurements of the HI and HeII Lyman-$\alpha$ forest opacities. We also present "equivalent-equilibrium" ionization/heating rates that mimic our fiducial UVB model for use in simulation codes that assume ionization equilibrium.

Comparison of cosmological simulations and deep submillimetre galaxy surveys

Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of infrared (IR)-emitting region, IR luminosity function and IRX--$\beta$ relation. After confirming that our models reproduce the local galaxy properties, we examine the STUDIES sample at $z\approx 1-4$, finding that the simulation reproduces the aforementioned quantities except for the $z\gtrsim 2$ IR luminosity function at the massive end ($\sim 10^{13}$ L$_{\odot}$). This means that the current simulation correctly reproduces the overall scaling between the size and luminosity (or star formation rate) of dusty region, but lacks extreme starburst phenomena at $z\gtrsim 2$. We also discuss extinction curves and possible AGN contribution.