-
The LOFAR Two-metre Sky Survey: The nature of the faint source population and SFR-radio luminosity relation using Prospector
Authors:
Soumyadeep Das,
Daniel J. B. Smith,
Paul Haskell,
Martin J. Hardcastle,
Philip N. Best,
Kenneth J. Duncan,
Marina I. Arnaudova,
Shravya Shenoy,
Rohit Kondapally,
Rachel K. Cochrane,
Alyssa B. Drake,
Gülay Gürkan,
Katarzyna Małek,
Leah K. Morabito,
Isabella Prandoni
Abstract:
Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospecto…
▽ More
Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospector, to model up to 26 bands of UV$-$far-infrared aperture-matched photometry for $\sim$31,000 sources in the ELAIS-N1 field from the LOFAR Two-Metre Sky Survey (LoTSS) Deep fields first data release. One of a new generation of SED-fitting codes, Prospector accounts for potential contributions from radiative active galactic nuclei (AGN) when estimating galaxy properties, including star formation rates (SFRs) derived using non-parametric star formation histories. Combining this information with radio luminosities, we classify 92 per cent of the radio sources as a star-forming galaxy, high/low-excitation radio galaxy, or radio-quiet AGN and study the population demographics as a function of 150 MHz flux density, luminosity, SFR, stellar mass, redshift and apparent $r$-band magnitude. Finally, we use Prospector SED fits to investigate the SFR$-$150 MHz luminosity relation for a sample of $\sim$$133,000~3.6~μ$m-selected $z<1$ sources, finding that the stellar mass dependence is significantly weaker than previously reported, and may disappear altogether at $\log_{10} (\mathrm{SFR}/M_\odot~\mathrm{yr}^{-1}) > 0.5$. This approach makes it significantly easier to classify radio sources from LoTSS and elsewhere, and may have important implications for future studies of star-forming galaxies at radio wavelengths.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Beware the recent past: a bias in spectral energy distribution modelling due to bursty star formation
Authors:
P. Haskell,
S. Das,
D. J. B. Smith,
R. K. Cochrane,
C. C. Hayward,
D. Anglés-Alcázar
Abstract:
We investigate how the recovery of galaxy star formation rates (SFRs) using energy-balance spectral energy distribution (SED) fitting codes depends on their recent star formation histories (SFHs). We use the Magphys and Prospector codes to fit 6,706 synthetic spectral energy distributions of simulated massive galaxies at $1 < z < 8$ from the Feedback in Realistic Environments (FIRE) project. We id…
▽ More
We investigate how the recovery of galaxy star formation rates (SFRs) using energy-balance spectral energy distribution (SED) fitting codes depends on their recent star formation histories (SFHs). We use the Magphys and Prospector codes to fit 6,706 synthetic spectral energy distributions of simulated massive galaxies at $1 < z < 8$ from the Feedback in Realistic Environments (FIRE) project. We identify a previously-unknown systematic error in the Magphys results due to bursty star formation: the derived SFRs can differ from the truth by as much as 1 dex, at large statistical significance ($>5σ$), depending on the details of their recent SFH. SFRs inferred using Prospector with non-parametric SFHs do not exhibit this trend. We show that using parametric SFHs (pSFHs) causes SFR uncertainties to be underestimated by a factor of up to $5\times$. Although this undoubtedly contributes to the significance of the systematic, it cannot explain the largest biases in the SFRs of the starbursting galaxies, which could be caused by details of the stochastic prior sampling or the burst implementation in the Magphys libraries. We advise against using pSFHs and urge careful consideration of starbursts when SED modelling galaxies where the SFR may have changed significantly over the last ~100 Myr, such as recently quenched galaxies, or those experiencing a burst. This concern is especially relevant, e.g. when fitting JWST observations of very high-redshift galaxies.
△ Less
Submitted 8 March, 2024; v1 submitted 24 October, 2023;
originally announced October 2023.
-
Energy balance SED modelling can be effective at high redshifts regardless of UV-FIR offsets
Authors:
P. Haskell,
D. J. B. Smith,
R. K. Cochrane,
C. C. Hayward,
D. Anglés-Alcázar
Abstract:
Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially `decoupled'). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such…
▽ More
Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially `decoupled'). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such galaxies' true properties. In this work, we test how well the SED fitting code Magphys can recover the stellar mass, star formation rate (SFR), specific SFR, dust mass and luminosity by fitting 6,706 synthetic SEDs generated from four zoom-in simulations of dusty, high-redshift galaxies from the FIRE project via dust continuum radiative transfer. Comparing our panchromatic results (using wavelengths 0.4-500$μ$m, and spanning $1<z<8$) with fits based on either the starlight ($λ_\mathrm{eff} \le 2.2\,μ$m) or dust ($\ge 100\,μ$m) alone, we highlight the power of considering the full range of multi-wavelength data alongside an energy balance criterion. Overall, we obtain acceptable fits for 83 per cent of the synthetic SEDs, though the success rate falls rapidly beyond $z \approx 4$, in part due to the sparser sampling of the priors at earlier times since SFHs must be physically plausible (i.e. shorter than the age of the Universe). We use the ground truth from the simulations to show that when the quality of fit is acceptable, the fidelity of Magphys estimates is independent of the degree of UV\FIR offset, with performance very similar to that previously reported for local galaxies.
△ Less
Submitted 14 September, 2023;
originally announced September 2023.
-
The LOFAR Two-metre Sky Survey: the radio view of the cosmic star formation history
Authors:
R. K. Cochrane,
R. Kondapally,
P. N. Best,
J. Sabater,
K. J. Duncan,
D. J. B. Smith,
M. J. Hardcastle,
H. J. A. Röttgering,
I. Prandoni,
P. Haskell,
G. Gürkan,
G. K. Miley
Abstract:
We present a detailed study of the cosmic star formation history over $90$ per cent of cosmic time ($0\lesssim z\lesssim4$), using deep, radio continuum observations that probe star formation activity independent of dust. The Low Frequency Array Two Metre Sky Survey has imaged three well-studied extragalactic fields, Elais-N1, Boötes and the Lockman Hole, reaching $\sim20\,μ\rm{Jy/beam}$ rms sensi…
▽ More
We present a detailed study of the cosmic star formation history over $90$ per cent of cosmic time ($0\lesssim z\lesssim4$), using deep, radio continuum observations that probe star formation activity independent of dust. The Low Frequency Array Two Metre Sky Survey has imaged three well-studied extragalactic fields, Elais-N1, Boötes and the Lockman Hole, reaching $\sim20\,μ\rm{Jy/beam}$ rms sensitivity at $150\,\rm{MHz}$. The availability of high-quality ancillary data from ultraviolet to far-infrared wavelengths has enabled accurate photometric redshifts and the robust separation of radio-bright AGN from their star-forming counterparts. We capitalise on this unique combination of deep, wide fields and robustly-selected star-forming galaxies to construct radio luminosity functions and derive the cosmic star formation rate density. We carefully constrain and correct for scatter in the $L_{150\,\rm{MHz}}-\rm{SFR}$ relation, which we find to be $\sim0.3\,\rm{dex}$. Our derived star formation rate density lies between previous measurements at all redshifts studied. We derive higher star formation rate densities between $z\sim0$ and $z\sim3$ than are typically inferred from short wavelength emission; at earlier times, this discrepancy is reduced. Our measurements are generally in good agreement with far-infrared and radio-based studies, with small offsets resulting from differing star formation rate calibrations.
△ Less
Submitted 24 May, 2023;
originally announced May 2023.
-
The LOFAR Two-metre Sky Survey Deep Fields Data Release 1: V. Survey description, source classifications and host galaxy properties
Authors:
P. N. Best,
R. Kondapally,
W. L. Williams,
R. K. Cochrane,
K. J. Duncan,
C. L. Hale,
P. Haskell,
K. Malek,
I. McCheyne,
D. J. B. Smith,
L. Wang,
A. Botteon,
M. Bonato,
M. Bondi,
G. Calistro Rivera,
F. Gao,
G. Gurkan,
M. J. Hardcastle,
M. J. Jarvis,
B. Mingo,
H. Miraghaei,
L. K. Morabito,
D. Nisbet,
I. Prandoni,
H. J. A. Rottgering
, et al. (4 additional authors not shown)
Abstract:
Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried…
▽ More
Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Radio spectral properties of star-forming galaxies between 150-5000MHz in the ELAIS-N1 field
Authors:
Fangxia An,
M. Vaccari,
P. N. Best,
E. F. Ocran,
C. H. Ishwara-Chandra,
A. R. Taylor,
S. K. Leslie,
H. J. A. Röttgering,
R. Kondapally,
Paul Haskell,
J. D. Collier,
M. Bonato
Abstract:
By combining high-sensitivity LOFAR 150MHz, uGMRT 400MHz and 1,250MHz, GMRT 610MHz, and VLA 5GHz data in the ELAIS-N1 field, we study the radio spectral properties of radio-detected star-forming galaxies (SFGs) at observer-frame frequencies of 150-5,000MHz. We select ~3,500 SFGs that have both LOFAR 150MHz and GMRT 610MHz detections, and obtain a median two-point spectral index of…
▽ More
By combining high-sensitivity LOFAR 150MHz, uGMRT 400MHz and 1,250MHz, GMRT 610MHz, and VLA 5GHz data in the ELAIS-N1 field, we study the radio spectral properties of radio-detected star-forming galaxies (SFGs) at observer-frame frequencies of 150-5,000MHz. We select ~3,500 SFGs that have both LOFAR 150MHz and GMRT 610MHz detections, and obtain a median two-point spectral index of $α_{150}^{610}=-0.51\pm0.01$. The photometric redshift of these SFGs spans $z=0.01-6.21$. We also measure the two-point radio spectral indices at 150-400-610-1,250MHz and 150-610-5,000MHz respectively for the GMRT 610-MHz-detected SFGs, and find that, on average, the radio spectrum of SFGs is flatter at low frequency than at high frequency. At observer-frame 150-5,000MHz, we find that the radio spectrum slightly steepens with increasing stellar mass. However, we only find that the radio spectrum flattens with increasing optical depth at $V$-band at $ν<1$GHz. We suggest that spectral ageing due to the energy loss of CR electrons and thermal free-free absorption could be among the possible main physical mechanisms that drive the above two correlations respectively. In addition, both of these mechanisms could physically explain why the radio spectrum is flatter at low frequency than at high frequency.
△ Less
Submitted 15 February, 2024; v1 submitted 13 March, 2023;
originally announced March 2023.
-
Cosmic evolution of low-excitation radio galaxies in the LOFAR Two-meter Sky Survey Deep Fields
Authors:
R. Kondapally,
P. N. Best,
R. K. Cochrane,
J. Sabater,
K. J. Duncan,
M. J. Hardcastle,
P. Haskell,
B. Mingo,
H. J. A. Röttgering,
D. J. B. Smith,
W. L. Williams,
M. Bonato,
G. Calistro Rivera,
F. Gao,
C. L. Hale,
K. Małek,
G. K. Miley,
I. Prandoni,
L. Wang
Abstract:
Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present…
▽ More
Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present the first measurement of the evolution of the radio luminosity function (LF) of LERGs out to $z\sim2.5$; this shows relatively mild evolution. We split the LERGs into those hosted by quiescent and star-forming galaxies, finding a new dominant population of LERGs hosted by star-forming galaxies at high redshifts. The incidence of LERGs in quiescent galaxies shows a steep dependence on stellar-mass out to $z \sim1.5$, consistent with local Universe measurements of accretion occurring from cooling of hot gas haloes. The quiescent-LERGs dominate the LFs at $z<1$, showing a strong decline in space density with redshift, tracing that of the available host galaxies, while there is an increase in the characteristic luminosity. The star-forming LERG LF increases with redshift, such that this population dominates the space densities at most radio-luminosities by $z \sim 1$. The incidence of LERGs in star-forming galaxies shows a much weaker stellar-mass dependence, and increases with redshift, suggesting a different fuelling mechanism compared to their quiescent counterparts, potentially associated with the cold gas supply present in the star-forming galaxies.
△ Less
Submitted 22 April, 2022; v1 submitted 15 April, 2022;
originally announced April 2022.
-
The LOFAR Two-metre Sky Survey Deep fields: A new analysis of low-frequency radio luminosity as a star-formation tracer in the Lockman Hole region
Authors:
M. Bonato,
I. Prandoni,
G. De Zotti,
P. N. Best,
M. Bondi,
G. Calistro Rivera,
R. K. Cochrane,
G. Gürkan,
P. Haskell,
R. Kondapally,
M. Magliocchetti,
S. K. Leslie,
K. Malek,
H. J. A. Röttgering,
D. J. B. Smith,
C. Tasse,
L. Wang
Abstract:
We have exploited LOFAR deep observations of the Lockman Hole field at 150 MHz to investigate the relation between the radio luminosity of star-forming galaxies (SFGs) and their star formation rates (SFRs), as well as its dependence on stellar mass and redshift. The adopted source classification, SFRs and stellar masses are consensus estimates based on a combination of four different SED fitting m…
▽ More
We have exploited LOFAR deep observations of the Lockman Hole field at 150 MHz to investigate the relation between the radio luminosity of star-forming galaxies (SFGs) and their star formation rates (SFRs), as well as its dependence on stellar mass and redshift. The adopted source classification, SFRs and stellar masses are consensus estimates based on a combination of four different SED fitting methods. We note a flattening of radio spectra of a substantial minority of sources below $\sim 1.4 $ GHz. Such sources have thus a lower "radio-loudness" level at 150 MHz than expected from extrapolations from 1.4 GHz using the average spectral index. We found a weak trend towards a lower $\hbox{SFR}/L_{150 \rm MHz}$ ratio for higher stellar mass, $M_\star$. We argue that such a trend may account for most of the apparent redshift evolution of the $L_{150 \rm MHz}/\hbox{SFR}$ ratio, in line with previous work. Our data indicate a weaker evolution than found by some previous analyses. We also find a weaker evolution with redshift of the specific star formation rate than found by several (but not all) previous studies. Our radio selection provides a view of the distribution of galaxies in the $\hbox{SFR}$-$M_\star$ plane complementary to that of optical/near-IR selection. It suggests a higher uniformity of the star formation history of galaxies than implied by some analyses of optical/near-IR data. We have derived luminosity functions at 150 MHz of both SFGs and radio-quiet (RQ) AGN at various redshifts. Our results are in very good agreement with the T-RECS simulations and with literature estimates. We also present explicit estimates of SFR functions of SFGs and RQ AGN at several redshifts derived from our radio survey data.
△ Less
Submitted 14 September, 2021;
originally announced September 2021.
-
The bright end of the infrared luminosity functions and the abundance of hyperluminous infrared galaxies
Authors:
L. Wang,
F. Gao,
P. N. Best,
K. Duncan,
M. J. Hardcastle,
R. Kondapally,
K. Malek,
I. McCheyne,
J. Sabater,
T. Shimwell,
C. Tasse,
M. Bonato,
M. Bondi,
R. K. Cochrane,
D. Farrah,
G. Gurkan,
P. Haskell,
W. J. Pearson,
I. Prandoni,
H. J. A. Rottgering,
D. J. B. Smith,
M. Vaccari,
W. L. Williams
Abstract:
We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities > 10^13 L_solar, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume. We cross-match H…
▽ More
We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities > 10^13 L_solar, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume. We cross-match Herschel sources and LOFAR sources in Bootes (8.63 deg^2), Lockman Hole (10.28 deg^2), and ELAIS-N1 (6.74 deg^2) with rms sensitivities of around 32, 22, and 20 mJy per beam, respectively. We divide the matched samples into unique and multiple categories. For the multiple matches, we de-blend the Herschel fluxes using the LOFAR positions and the 150-MHz flux densities as priors. We perform spectral energy distribution (SED) fitting, combined with multi-wavelength counterpart identifications and photometric redshift estimates, to derive IR luminosities. The depth of the LOFAR data allows us to identify highly complete (around 92% completeness) samples of bright Herschel sources with a simple selection based on the 250 micron flux density (45, 40, and 35 mJy in Bootes, Lockman Hole, and ELAIS-N1, respectively). Most of the bright Herschel sources fall into the unique category (i.e. a single LOFAR counterpart). For the multiple matches, there is excellent correspondence between the radio emission and the far-IR emission. We find a good agreement in the IR LFs with a previous study out to z around 6 which used de-blended Herschel data. Our sample gives the strongest and cleanest indication to date that the population of HLIRGs has surface densities of around 5 to 18 / deg^2 (with variations due to a combination of the applied flux limit and cosmic variance) and an uncertainty of a factor of 2. In comparison, the GALFORM semi-analytic model significantly under-predicts the abundance of HLIRGs.
△ Less
Submitted 17 November, 2020;
originally announced November 2020.
-
The LOFAR Two-metre Sky Survey Deep fields: The star formation rate - radio luminosity relation at low frequencies
Authors:
D. J. B. Smith,
P. Haskell,
G. Gürkan,
P. N. Best,
M. J. Hardcastle,
R. Kondapally,
W. Williams,
K. J. Duncan,
R. K. Cochrane,
I. McCheyne,
H. J. A. Röttgering,
J. Sabater,
T. W. Shimwell,
C. Tasse,
M. Bonato,
M. Bondi,
M. J. Jarvis,
S. K. Leslie,
I. Prandoni,
L. Wang
Abstract:
In this paper, we investigate the relationship between 150MHz luminosity and star formation rate (the SFR-L150 relation) using 150MHz measurements for a near-infrared selected sample of 118,517 $z<1$ galaxies. New radio survey data offer compelling advantages for studying star formation in galaxies, with huge increases in sensitivity, survey speed and resolution over previous generation surveys, a…
▽ More
In this paper, we investigate the relationship between 150MHz luminosity and star formation rate (the SFR-L150 relation) using 150MHz measurements for a near-infrared selected sample of 118,517 $z<1$ galaxies. New radio survey data offer compelling advantages for studying star formation in galaxies, with huge increases in sensitivity, survey speed and resolution over previous generation surveys, and remaining impervious to extinction. The LOFAR Surveys Key Science Project is transforming our understanding of the low-frequency radio sky, with the 150MHz data over the ELAIS-N1 field reaching an RMS sensitivity of 20uJy/beam over 10 deg$^2$ at 6" resolution. All of the galaxies studied have SFR and stellar mass estimates derived from energy balance SED fitting, using redshifts and aperture-matched forced photometry from the LOFAR Two-metre Sky Survey (LoTSS) deep fields data release. The impact of active galactic nuclei is minimised by leveraging the deep ancillary data alongside outlier-resistant median-likelihood methods. We find a linear and non-evolving SFR-L150 relation, apparently consistent with expectations based on calorimetric arguments, down to the lowest SFRs. However, we also recover compelling evidence for stellar mass dependence in line with previous work on this topic, in the sense that higher mass galaxies have a larger 150MHz luminosity at a given SFR, suggesting that the overall agreement with calorimetric arguments may be a coincidence. We conclude that in the absence of AGN, 150MHz observations can be used to measure accurate galaxy SFRs out to $z=1$ at least, but it is necessary to account for stellar mass in order to obtain 150MHz-derived SFRs accurate to <0.5 dex. Our best-fit relation is $\log_{10} (L_\mathrm{150 MHz} / W\,Hz^{-1}) = (0.90\pm 0.01) \log_{10}(ψ/M_\odot\,\mathrm{yr}^{-1}) + (0.33 \pm 0.04) \log_{10} (M/10^{10}M_\odot) + 22.22 \pm 0.02$. (Abridged)
△ Less
Submitted 16 November, 2020;
originally announced November 2020.