scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Infrared luminosity functions based on 18 mid-infrared bands: revealing cosmic star formation history with AKARI and Hyper Suprime-Cam

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

Summary (3 min read)

1 Introduction

  • Mid-infrared (mid-IR) is one of the less explored wavelengths due to the earth's atmosphere, and difficulties in developing sensitive detectors.
  • To overcome these problems, the authors have newly obtained deeper optical data over the entire AKARI NEP wide field, using the Hyper-Suprime Cam on the Subaru telescope.
  • Using the deeper optical data, in this paper, the authors measure mid-infrared galaxy LFs, and estimate total IR LFs (based on the mid-IR SED fitting) from the entire AKARI NEP field.

2 Data

  • To rectify the situation and to fully exploit the AKARI's spacebased data, the authors carried out an optical survey of the AKARI NEP wide field (PI:Goto) using Subaru's new Hyper Suprime-Cam (HSC; Miyazaki et al. 2018) in five optical bands (g, r, i, z, and y, Oi et al. 2018 submitted) .
  • The HSC has a field-of-view (FoV) of 1.5 deg in diameter, covered with 104 red-sensitive CCDs.
  • It has the largest FoV among optical cameras on 8m-class telescopes, and can cover the AKARI NEP wide field (5.4 deg 2 ) with only 4 FoV (Fig. 1 ).
  • See Oi et al. (2018, submitted) for more details of the observation and data reduction.
  • Subaru telescope does not have u * -band capability, while it is critically important to accurately estimate photometric redshifts (photo-z) of low-z galaxies.

3 Analysis

  • Uncertainties of the LF values include fluctuations in the number of sources in each luminosity bin, the photometric redshift uncertainties, the k-correction uncertainties, and the flux errors.
  • To estimate errors, the authors used Monte Carlo simulations from 1000 simulated catalogs.
  • A new flux is also assigned following a Gaussian distribution with the width of flux error.
  • The smaller data points at the faint ends are adopted from the NEP deep field, where AKARI data are deeper (Goto et al. 2015) , and are included in the fit.
  • Vertical arrows show the 8µm luminosity corresponding to the flux limit at the central redshift in each redshift bin.

4.1 The 8µm LF

  • The authors first present monochromatic 8µm LFs, because the 8µm luminosity (L8µm) has been known as a good indicator of the TIR luminosity (Babbedge et al.
  • Often in previous work, SED based extrapolation was needed to estimate the 8µm luminosity.
  • This is not the case for the analysis present in this paper.
  • The smaller data points at the faint ends are adopted from the NEP deep field, where AKARI data are deeper (Goto et al. 2015) , and are included in the fit.
  • Vertical arrows show the 12µm luminosity corresponding to the flux limit at the central redshift in each redshift bin.

templates (2% from the sample).

  • The authors corrected for the completeness using Kim et al. (2012) (25% correction at maximum, with their selection to the 80% completeness limits).
  • Then, the 1/Vmax method was used to compensate for the flux limit.
  • Various previous studies are shown with dashed lines for comparison.
  • Interestingly, the 8µm LFs peaks in the 3rd bin (z∼1), then declines toward z∼2.

4.2 12µm LF

  • The 12µm luminosity L12µm) is also known to correlate well with the TIR luminosity (Spinoglio et al.
  • AKARI's advantage still holds in not needing extrapolation based on SED models.
  • Various previous studies are shown in dash-dotted lines.
  • Similar to the 8µm LF, the evolution becomes less evident between the two higher redshift bins.

4.3 Total IR LFs estimated from mid-IR SED fit

  • The authors caution readers that estimation of the LTIR involves extrapolation to the far-IR wavelength range based on the SED models, and thus invites associated uncertainty, as they further discuss in Section 5.
  • The L18W flux (Matsuhara et al. 2006 ) are used to apply the 1/Vmax method, because it is a wide, sensitive filter (but using the L15 flux limit does not change their main results).
  • For clarity, the authors separated LFs in four different panels at each redshift bin.
  • The TIR LFs show a strong evolution compared to local LFs, but again turns over at z > 1.2.

4.4.1 Total IR Luminosity Density from L8µm

  • LFs First, the authors estimate Total IR Luminosity Density from L8µm LFs.
  • Possible SED evolution, and the presence of AGN will induce further uncertainty.
  • Murata et al. (2014) also reported that L8µm/LTIR is constant at below the main sequence, while it decreases with starburstiness at above the main sequence, concluding that starburst galaxies have deficient PAH emission compared with main-sequence galaxies.
  • Overplotted previous studies are taken from Le Floc'h et al. ( 2005) in the dark-green, dash-dotted line, Magnelli et al. (2013) in the dark-red, dash-dotted line, Huynh et al. (2007) in the dark-yellow, dash-dotted line, Gruppioni et al. (2013) in the pink, dash-dotted line at several redshifts as marked in the figure.

4.4.2 Total IR Luminosity Density from L12µm LFs

  • Due to the same reasons as L8µm (improved statistics, and availability of 140 and 160µm), the authors use the following conversion (Goto et al. 2011b) .
  • The authors caution readers again here for the use of a single conversion for varieties of galaxies with different SFR at different redshifts.
  • Results should be interpreted with this uncertainty in mind.

4.4.3 Integration to TIR density

  • The derived total LFs are multiplied by LTIR and integrated to measure the TIR density ( ΩTIR).
  • Following their previous work, the authors use a double-power law.
  • With the lowest redshift LF, the authors first fit the normalization (Φ * ) and slopes (α, β).
  • The authors also note that ΩIR from 12µm is sensitive to the faint-end slope of 12µm LFs.
  • Much deeper observations are awaited to clarify the issue.

5 Discussion

  • The conversions are based on local star-forming galaxies.
  • The pink dashed line shows the total estimate of IR (TIR LF) and UV (Schiminovich et al. 2005) .
  • Following the results in the literature discussed in Section 4.3, in this section, the authors compare LT IR estimated from L8µm and L12µm from equations 1 and 2 in three overlapping redshift ranges in Fig. 6 using their data.
  • One can immediately notice that the relation deviates at logLTIR >12 (or equivalently at z > 1).
  • 4 and 5 between midand far-IR measurements could be the result of the change in the SED, rather than incorrect measurements on either.

6 Summary

  • Previously AKARI NEP wide field lacked deep optical photometry, and thereby, accurate photo-z, despite the presence of space-based 9-band mid-IR photometry from AKARI.
  • To rectify the situation, the authors have obtained deep optical 5-band imaging covering the entire 5.4 deg 2 of the NEP wide field, using the new Hyper Suprime-Cam mounted on the Subaru 8m telescope.
  • Thanks to the large area coverage, the brightends are better-determined.

Did you find this useful? Give us your feedback

Figures (8)

Content maybe subject to copyright    Report

Open Research Online
The Open University’s repository of research publications
and other research outputs
Infrared luminosity functions based on 18 mid-infrared
bands: revealing cosmic star formation history with
AKARI and Hyper Suprime-Cam*
Journal Item
How to cite:
Goto, Tomotsugu; Oi, Nagisa; Utsumi, Yousuke; Momose, Rieko; Matsuhara, Hideo; Hashimoto, Tetsuya;
Toba, Yoshiki; Ohyama, Youichi; Takagi, Toshinobu; Chiang, Chia-Ying; Kim, Seong Jin; Kilerci Eser, Ece; Malkan,
Matthew; Kim, Helen; Miyaji, Takamitsu; Im, Myungshin; Nakagawa, Takao; Jeong, Woong-Seob; Pearson, Chris;
Barrufet, Laia; Sedgwick, Chris; Burgarella, Denis; Buat, Veronique and Ikeda, Hiroyuki (2019). Infrared luminosity
functions based on 18 mid-infrared bands: revealing cosmic star formation history with AKARI and Hyper Suprime-
Cam*. Publications of the Astronomical Society of Japan, 71(2), article no. 30.
For guidance on citations see FAQs.
c
2019 The Authors
https://creativecommons.org/licenses/by-nc-nd/4.0/
Version: Accepted Manuscript
Link(s) to article on publisher’s website:
http://dx.doi.org/doi:10.1093/pasj/psz009
Copyright and Moral Rights for the articles on this site are retained by the individual authors and/or other copyright
owners. For more information on Open Research Online’s data policy on reuse of materials please consult the policies
page.
oro.open.ac.uk

arXiv:1902.02801v1 [astro-ph.GA] 7 Feb 2019
Publ. Astron. Soc. Japan (2014) 00(0), 1–10
doi: 10.1093/pasj/xxx000
1
Infrared luminosity functions based on 18
mid-infrared bands: revealing cosmic star
formation history with AKARI and Hyper
Suprime-Cam
Tomotsugu GOTO
1
, Nagisa OI
2
, Yousuke UTSUMI
3
, Rieko MOMOSE
1,4
,
Hideo MATSUHARA
5
, Tetsuya HASHIMOTO
1
, Yoshiki TOBA
6
, Youichi
OHYAMA
6
, Toshinobu TAKAGI
7
, Chia Ying CHIANG
1
, Seong Jin KIM
1
, Ece
KILERCI ESER
1
, Matthew MALKAN
8
, Helen KIM
8
, Takamitsu MIYAJI
9
,
Myungshin IM
10
, Takao NAKAGAWA
5
, Woong-Seob JEONG
11,12
, Chris
PEARSON
13,14
, Laia BARRUFET
15
, Chris SEDGWICK
14
, Denis BURGARELLA
16
,
Veronique BUAT
16
and Hiroyuki IKEDA
17
1
National Tsing Hua University, No. 101, Section 2, Kuang-Fu Road, Hsinchu, Taiwan 30013
2
Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
3
Kavli Institute for Particle Astrophysics and Cosmology (KIPAC), SLAC National Accelerator
Laboratory, Stanford University, SLAC, 2575 Sand Hill Road, Menlo Park, CA 94025, USA
4
Department of Astronomy, School of Science, The University of Tokyo 7-3-1 Hongo,
Bunkyo-ku, Tokyo 113-0033, JAPAN
5
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 3-1-1
Yoshinodai, Chuo, Sagamihara, Kanagawa 252-5210, Japan
6
Academia Sinica Institute of Astronomy and Astrophysics, P.O. Box 23-141, Taipei 10617,
Taiwan
7
Japan Space Forum, 3-2-1, Kandasurugadai, Chiyoda-k u, Tokyo 101-0062 Japan
8
Department of Physics and Astronomy, UCLA, Los Angeles, CA, 90095-1547, USA
9
Insitituto de Astronom´ıa, U niversidad Nacional Aut
´
onoma de M
´
exico
10
Astronomy Program, Department of Physics & Astr onomy, FPRD, Seoul National University,
Shillim-Dong, Kwanak-Gu, Seoul 151-742, Korea
11
Korea Astronomy and Space Science Institute (KASI), 776 Daedeok-daero, Yuseong-gu,
Daejeon 34055, Korea
12
Korea University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113,
Korea
13
RAL Space, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire OX11 0QX, UK
14
The Open University, Milton Keynes, MK7 6AA, UK
15
European Space Astronomy Centre, 28691 Villanueva de la Canada, Spain
16
Aix-Marseille Universit, CNRS LAM (Laboratoire dAstrophysique de Marseille) UMR 7326,
13388 Mars eille, France
c
2014. Astronomical Society of Japan.

2 Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0
17
National Astronomical Observatory, 2-21-1 Osawa, Mitaka, Tokyo, Japan
E-mail: tomo@gapp.nthu.edu.tw
Received 2018 June 30; Accepted 2019 January 21
Abstract
Much of the star formation is obscured by dust. For the complete understanding of the cosmic
star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep
mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP)
field (5.4 deg
2
). This took significant amount of satellite’s lifetime, 10% of the entire pointed
observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry,
we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric
coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical
imaging (25.9ABmag) in a small area of 1.0 deg
2
. As a result, there remained thousands of
AKARI’s infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru
telescope, we obtained deep enough optical images of the entire AKARI NEP eld in 5 broad
bands (g 27.5mag). These provided photometric redshift, and thereby IR luminosity for the
previously undetected faint AKARI IR sour ces. Combined with the accurate mid-IR luminosity
measurement, we constr ucted mid-IR LFs, and thereby performed a census of dust-obs cured
CSFH in the entire AKARI NEP eld. We have measured restframe 8µm, 12µm luminosity
functions (LFs), and estimated total infrared LFs at 0.35<z<2.2. Our results are consistent with
our previous work, but with much reduced statistical err ors thanks to the large area coverage
of the new data. We have possibly witnessed the turnover of CSFH at z 2.
Key words: AKARI, infrared galaxies, cosmic star formation history
1 Introduction
Mid-infrared ( mid-IR) is one of the less explored wavelengths
due to the earth’s atmosphere, and difficulties in developing sen-
sitive detectors. NASA’s S pitzer and WISE space telescopes
only had four filters in the mid-IR wavelength range, hamper-
ing studies of distant galaxies.
AKARI space telescope has a potential to revolutionize the
field. Using its 9 continuous mid-IR filters (2-24µm), AKARI
performed a deep imaging survey in the North Ecliptic Pole
(NEP) field over 5.4 deg
2
. Using AKARI’s 9 mid-IR band pho-
tometry, mid-IR SED diagnosis can be performed for thousands
of galaxies, for the first time, over the large enough area to over-
come cosmic variance. Environmental effects on galaxy evolu-
tion can be also investigated with the large volume coverage
(Koyama et al. 2008; Goto et al. 2010a).
However, previously, we were limited by a poor optical cov-
erage both in area and depths. Over this wide area, only shal-
low optical/NIR imaging data have been available (Hwang et al.
2007; Jeon et al. 2010, 2014). Deep optical images are limited
to the central 0.25 deg
2
.
To overcome these problems, we have newly obtained
deeper optical data over the entire AKARI NEP wide field, us-
ing the Hyper-Suprime Cam on the Subaru telescope. Using
the deeper optical data, in this paper, we measure mid-infrared
galaxy LFs, and estimate total IR LFs (based on the mid-IR
SED fitting) from the entire AKARI NEP field. Unless oth-
erwise stated, we assume a cosmology with (h,
m
,
Λ
) =
(0.7,0.3,0.7).
2 Data
To rectify the situation and to fully exploit the AKARI’s space-
based data, we carr ied out an optical survey of the AKARI NEP
wide field (PI:Goto) using Subaru’s new Hyper Suprime-Cam
(HSC; Miyazaki et al. 2018) in five optical bands (g, r,i, z, and
y, Oi et al. 2018 submitted). The HSC has a field-of-view (FoV)
of 1.5 deg in diameter, covered with 104 red-sensitive CCDs. It
has the largest FoV among optical cameras on 8m-class tele-
scopes, and can cover the AKARI NEP wide field (5.4 deg
2
)
with only 4 FoV (Fig.1). The 5 sigma limiting magnitudes are
27.18, 26.71, 26.10, 25.26, and 24.78 mag [AB] in g,r,i,z, and
y-bands, respectively. See Oi et al. (2018, submitted) for more
details of the observation and data reduction.
Subaru telescope does not have u
-band capability, while it
is critically important to accurately estimate photometric red-
shifts (photo-z) of low-z galaxies. Therefore, we obtained u
-
band image of the AKARI NEP wide eld using the Megaprime
camera of Canada France Hawaii Telescope (PI :Goto, Goto

Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 3
45:00.0
17:50:00.0
55:00.0
18:00:00.0
05:00.0
10:00.0
15:00.0
65:00:00.0
30:00.0
66:00:00.0
30:00.0
67:00:00.0
30:00.0
Right ascension
Declination
Fig. 1. HSC three color (g, r, i) composite image of the NEP wide field (5.4 deg
2
). The AKARI NEP wide data exist within the white circle.
et al. 2017). Combining the optical six bands, we have obtained
accurate photo-z in the AKARI NEP field (Oi et al. 2018, sub-
mitted). To the detection limit in L18W filter (18.3 ABmag,
Kim et al. 2012), we have 5078 infrared sources.
In addition to the AKARI’s 9 mid-IR bands, in the AKARI
NEP field, there exit archival deep Spitzer (IRAC1,2,3,4 and
MIPS24, Nayyeri et al. 2018) and WISE (W 1, W 2, W 3, and
W 4) images as well. By combining all available mid-IR bands,
in total we used 18 mid-IR bands, which are one of the most
comprehensive mid-IR data sets for thousands of galaxies.
3 Analysis
To compute LFs, we use the 1/V
max
method, following Goto
et al. (2010b, 2015). Uncertainties of the LF values include
fluctuations in the number of sources in each luminosity bin,
the photometric redshift uncertainties, the k-correction uncer-
tainties, and the flux errors. To estimate errors, we used Monte
Carlo simulations from 1000 simulated catalogs. Each sim-
ulated catalog contains the same number of sources. These
sources are assigned with a new redshift, to follow a Gaussian
distribution centered at the photo-z with the width of z/(1 +
z) ( 0.060, Oi et al. in preparation). A new flux is also as-
signed following a Gaussian distribution with the width of flux
error. For total infrared (TIR) LF errors, we re-performed the
SED fit for the 1000 simulated catalogs. Note that total in-
frared luminosity is estimated based on mid-IR SED tting al-
though we have intensive 18-band filter coverage in mid-IR, as
explained in Section 4.3. We ignored the cosmic variance due
to our much improved volume coverage. All the other err ors

4 Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0
Fig. 2. Restframe 8µm LFs based on the AKARI NEP wide field. The blue
diamonds, the purple triangles, the red squares, and the orange crosses
show the 8µm LFs at 0.28 < z < 0.47, 0.65 < z < 0.90, 1.09 < z < 1.41,
and 1.78 < z < 2.22, respectively. AKARI’s MIR filters can observe rest-
frame 8µm at these redshifts in a corresponding filter. Error bars are esti-
mated from the Monte Carlo simulations (§3). The dotted l ines show analytic
fits with a double-power law. The smaller data points at the faint ends are
adopted from the NEP deep field, where AKARI data are deeper (Goto et al.
2015), and are included in the fit. Vertical arrows show the 8µm luminos-
ity corresponding to the flux limit at the central redshift in each redshift bin.
Overplotted LFs are B abbedge et al. (2006) in the pink dash-dotted lines,
Caputi et al. (2007) in the cyan dash-dotted lines, Huang et al. (2007) in
the dark-yellow dash-dotted lines, F u et al. (2010), in the dark green dash-
dotted line, and Kim et al. (2015) in the bright green dash-dotted line. Best-fit
parameters are presented in Table 1.
described above are added to the Poisson errors for each LF bin
in quadrature.
4 Results
4.1 The 8µm LF
We firs t present monochromatic 8µm LFs, because the 8µm lu-
minosity (L
8µm
) has been known as a good indicator of the TIR
luminosity (Babbedge et al. 2006; Huang et al. 2007; Goto et al.
2011a).
An advantage of AKARI is that we do not need k-correction
because one of the continuous filters always convert the rest-
frame 8µm at our redshift range of 0.28 < z < 2.22. Often in
previous work, SED based extrapolation was needed to estimate
the 8µm luminosity. This was often the largest uncertainty. This
is not the case for the analysis present in this paper.
To estimate the restframe 8µm LFs, we followed our previ-
ous method in Goto et al. (2015) as we briefly summarize below.
We used sources down to 80% completeness limits ( Kim et al.
2012). Galaxies are excluded when SEDs were better fit to QSO
Fig. 3. Restframe 12µm LFs based on the AKARI NEP wide field.
Luminosity unit is logarithmic solar luminosity (L
). The blue diamonds,
the purple triangles, and the red squares show the 12µm LFs at 0.15 < z <
0.35, 0.38 < z < 0.62, and 0.84 < z < 1.16, respectively. The smaller data
points at the faint ends are adopted from the NEP deep field, where AKARI
data are deeper (Goto et al. 2015), and are included in the fit. Vertical arrows
show the 12µm luminosity corresponding to the flux limit at the central red-
shift in each redshift bin. Overplotted LFs are P
´
erez-Gonz
´
alez et al. (2005)
at z=0.3, 0.5 and 0.9 in the dark-cyan dash-dotted l ines, Toba et al. (2014)
at 0< z <0.05 based on WISE in the dark green dash-dotted lines, Rush
et al. (1993) at 0< z <0.3 in the light green dash-dotted lines, and Kim et al.
(2015) at 0< z <0.3 in the pink dash-dotted line. Note Rush et al. (1993) is
at higher redshifts than Toba et al. (2014). Best-fit parameters are presented
in Table 1.
templates (2% from the sample).
We corrected for the completeness using Kim et al. (2012)
(25% correction at maximum, with our selection to the 80%
completeness limits). Four redshift bins of 0.28< z <0.47,
0.65< z <0.90, 1.09< z <1.41, and 1.78< z <2.22, were used,
following our previous work. Then, the 1/V
max
method was
used to compensate for the flux limit.
The resulting restframe 8µm LFs are shown in Fig. 2. The
arrows show the flux limit at the m edian redshift in bin. We
performed the Monte Carlo simulation to obtain errors. They
are smaller than in our previous work (Goto et al. 2010b, 2015)
thanks to the improved area coverage. The faint end marked
with smaller data points are adopted from the NEP deep field,
where AKARI data are deeper (Goto et al. 2015).
Various previous studies are shown with dashed lines for
comparison. Compared to the local LF, our 8µm LFs show
strong evolution in luminosity up to z 0.9. Interestingly, the
8µm LFs peaks in the 3rd bin ( z1), then declines toward z2.

Citations
More filters
Journal ArticleDOI
TL;DR: In this paper, a set of infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies was gathered and the output parameters for all the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M{dust}$ and the $L_{IR}$.
Abstract: We gather infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies. Draine & Li (2007) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the models parameters and that a strong constraint on the 20-60 microns range is mandatory to estimate the relative contribution of the photo-dissociation regions to the IR SED. The SED models tend to systematically under-estimate the observed 500 microns flux densities, especially for low mass systems. We provide the output parameters for all of the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M_{dust}$ and the $L_{IR}$. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the integrated properties of galaxies, such as $M_*$, SFR, metallicity, H$\alpha$ and H-band surface brightness, and the FUV attenuation. A good correlation between the fraction of PAH and the metallicity is found implying a weakening of the PAH emission in galaxies with low metallicities. The intensity of the IRSF and the H-band and H$\alpha$ surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star forming regions and the diffuse evolved population. We use these results to provide a new set of IR templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z=0 reference for cosmological studies. For the same purpose, we put our sample on the SFR-$M_*$ diagram. The templates are compared to the most popular IR SED libraries, enlightening a large discrepancy between all of them in the 20-100 microns range.

134 citations

Journal ArticleDOI
TL;DR: In this article, a sample of 56 sources serendipitously detected in ALMA band 7 as part of the ALMA Large Program to Investigate CII at Early Times (ALPINE) has been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z ǫ ≥ 6.
Abstract: Aims. We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7 as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE). These sources, detected in COSMOS and ECDFS, have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z ≃ 6.Methods. We looked for counterparts of the ALMA sources in all the available multi-wavelength (from HST to VLA) and photometric redshift catalogues. We also made use of deeper UltraVISTA and Spitzer source lists and maps to identify optically dark sources with no matches in the public catalogues. We used the sources with estimated redshifts to derive the 250 μ m rest-frame and total infrared (8–1000 μ m) LFs from z ≃ 0.5 to 6.Results. Our ALMA blind survey (860 μ m flux density range: ∼0.3–12.5 mJy) allows us to further push the study of the nature and evolution of dusty galaxies at high-z , identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint end of the infrared LF, showing little evolution from z ≃ 2.5 to z ≃ 6, and a “flat” slope up to the highest redshifts (i.e. 4.5 ≃ 2 and z ≃ 6, and significantly higher than the optical or ultra-violet derivations, showing a significant contribution of dusty galaxies and obscured star formation at high-z . About 14% of all the ALPINE serendipitous continuum sources are found to be optically and near-infrared (near-IR) dark (to a depth K s ∼ 24.9 mag). Six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [C II] emitters at z ≃ 5. The six HST+near-IR dark galaxies with mid-IR counterparts are found to contribute about 17% of the total SFRD at z ≃ 5 and to dominate the high-mass end of the stellar mass function at z > 3.

50 citations

Journal ArticleDOI
TL;DR: In this article, the authors constrain the rest-frame FUV, NUV and U-band luminosity functions with unprecedented precision from z ∼ 0.2 to z ∼ 3.
Abstract: We constrain the rest-frame FUV (1546 A), NUV (2345 A), and U-band (3690 A) luminosity functions (LFs) and luminosity densities (LDs) with unprecedented precision from z ∼ 0.2 to z ∼ 3 (FUV, NUV) and z ∼ 2 (U band). Our sample of over 4.3 million galaxies, selected from the CFHT Large Area U-band Deep Survey (CLAUDS) and HyperSuprime-Cam Subaru Strategic Program (HSC-SSP) data lets us probe the very faint regime (down to MFUV, MNUV, MU ≃ −15 at low redshift), while simultaneously detecting very rare galaxies at the bright end down to comoving densities 1 it is due to the evolution of both $M^\star _{\rm UV}$ and the characteristic number density $\phi ^\star _{\rm UV}$. In contrast, the U-band LF has an excess of faint galaxies and is fitted with a double-Schechter form; $M^\star _{U}$, both $\phi ^\star _{U}$ components, and the bright-end slope evolve throughout 0.2 < z < 2, while the faint-end slope is constant over at least the measurable 0.05 < z < 0.6. We present tables of our Schechter parameters and LD measurements that can be used for testing theoretical galaxy evolution models and forecasting future observations.

34 citations


Cites background from "Infrared luminosity functions based..."

  • ...…have been possible for some time for high-z galaxies (e.g., Hughes et al. 1998; Chapman et al. 2005; Magnelli et al. 2013; Gruppioni et al. 2013; Goto et al. 2019), they do not yet provide significant insights at very high redshifts (z>∼ 6), nor for low-mass galaxies which have low SFRs and low…...

    [...]

Journal ArticleDOI
TL;DR: In this article, the authors estimate the detections of non-repeating and repeating fast radio bursts separately, based on latest observational constraints on their physical properties including the spectral indices, FRB luminosity functions, and their redshift evolutions.
Abstract: Fast radio bursts (FRBs) are mysterious extragalactic radio signals. Revealing their origin is one of the central foci in modern astronomy. Previous studies suggest that occurrence rates of non-repeating and repeating FRBs could be controlled by the cosmic stellar-mass density (CSMD) and star formation-rate density (CSFRD), respectively. The Square Kilometre Array (SKA) is one of the best future instruments to address this subject due to its high sensitivity and high-angular resolution. Here, we predict the number of FRBs to be detected with the SKA. In contrast to previous predictions, we estimate the detections of non-repeating and repeating FRBs separately, based on latest observational constraints on their physical properties including the spectral indices, FRB luminosity functions, and their redshift evolutions. We consider two cases of redshift evolution of FRB luminosity functions following either the CSMD or CSFRD. At $z\gtrsim2$, $z\gtrsim6$ and $z\gtrsim10$, non-repeating FRBs will be detected with the SKA at a rate of $\sim10^{4}$, $\sim10^{2}$, and $\sim10$ (sky$^{-1}$ day$^{-1}$), respectively, if their luminosity function follows the CSMD evolution. At $z\gtrsim1$, $z\gtrsim2$, and $z\gtrsim4$, sources of repeating FRBs will be detected at a rate of $\sim10^{3}$, $\sim10^{2}$, and $\lesssim10$ (sky$^{-1}$ day$^{-1}$), respectively, assuming that the redshift evolution of their luminosity function is scaled with the CSFRD. These numbers could change by about one order of magnitude depending on the assumptions on the CSMD and CSFRD. In all cases, abundant FRBs will be detected by the SKA, which will further constrain the luminosity functions and number density evolutions.

30 citations

Journal ArticleDOI
TL;DR: In this article, a star formation change parameter, SFRF$(5Myr/SFR$(800Myr), is defined, which is the ratio of the SFR averaged within the last 5 Myr to the Sfr averaged within 800 Myr.
Abstract: To investigate the variability of the star formation rate (SFR) of galaxies, we define a star formation change parameter, SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ which is the ratio of the SFR averaged within the last 5 Myr to the SFR averaged within the last 800 Myr. We show that this parameter can be determined from a combination of H$\alpha$ emission and H$\delta$ absorption, plus the 4000 A break, with an uncertainty of $\sim$0.07 dex for star-forming galaxies. We then apply this estimator to MaNGA galaxies, both globally within Re and within radial annuli. We find that galaxies with higher global SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ appear to have higher SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ at all galactic radii, i.e. that galaxies with a recent temporal enhancement in overall SFR have enhanced star formation at all galactic radii. The dispersion of the SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ at a given relative galactic radius and a given stellar mass decreases with the (indirectly inferred) gas depletion time: locations with short gas depletion time appear to undergo bigger variations in their star-formation rates on Gyr or less timescales. In Wang et al. (2019) we showed that the dispersion in star-formation rate surface densities $\Sigma_{\rm SFR}$ in the galaxy population appears to be inversely correlated with the inferred gas depletion timescale and interpreted this in terms of the dynamical response of a gas-regulator system to changes in the gas inflow rate. In this paper, we can now prove directly with SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ that these effects are indeed due to genuine temporal variations in the SFR of individual galaxies on timescales between $10^7$ and $10^9$ years rather than possibly reflecting intrinsic, non-temporal, differences between different galaxies.

17 citations


Cites background from "Infrared luminosity functions based..."

  • ...…galaxy population also known as the cosmic evolution of the star formation rate density (SFRD), is well established up to redshift of ∼9 (e.g. Lilly et al. 1996; Schiminovich et al. 2005; Bouwens et al. 2011, 2014; Madau & Dickinson 2014; Hagen et al. 2015; Alavi et al. 2016; Goto et al. 2019)....

    [...]

References
More filters
Journal ArticleDOI
TL;DR: In this article, the authors used a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent λ_(eff) = 24, 15, 12, and 8 μm monochromatic and bolometric IR luminosity functions (LF) over the full 0 < z < 2.5 redshift range.
Abstract: Context. Studies of the infrared (IR) emission of cosmic sources have proven essential to constraining the evolutionary history of cosmic star formation and the gravitational accretion of nuclear black holes, because many of these events occur inside heavily dust-extinguished environments. Aims. The Spitzer Space Telescope has provided a large amount of data to constrain the nature and cosmological evolution of infrared source populations. In the present paper we exploit a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent λ_(eff) = 24, 15, 12, and 8 μm monochromatic and bolometric IR luminosity functions (LF) over the full 0 < z < 2.5 redshift range. Methods. Our present analysis is based on a combination of data from deep Spitzer surveys of the VIMOS VLT Deep Survey (VVDS-SWIRE) and GOODS fields. To our limiting flux of S_(24) = 400 μJy, our sample derived from VVDS-SWIRE includes 1494 sources, and 666 and 904 sources brighter than S_(24) = 80 μJy are catalogued in GOODS-S and GOODS-N, respectively, for a total area of ~0.9 square degrees. Apart from a few galaxies, we obtain reliable optical identifications and redshifts for all these sources, providing a rich and robust dataset for our luminosity function determination. The final combined reliable sample includes 3029 sources, the fraction with photometric redshifts being 72% over all redshifts and almost all galaxies at z > 1.5. Based on the multiwavelength information available in these areas, we constrain the LFs at 8, 12, 15, and 24 μm. We also infer the total IR luminosities from our best-fit model of the observed SEDs of each source, and use this to derive the bolometric (8–1000 μm) LF and comoving volume emissivity to z ~ 2.5. Results. In the redshift interval 0 1. The mean redshift of the peak in the source number density shifts with luminosity: the brightest IR galaxies appear to form stars at earlier cosmic times (z > 1.5), while star formation in the less luminous galaxies continues until more recent epochs (z ~ 1 for L_(IR) < 10^(11)_☉), in overall agreement with similar analyses in the literature. Conclusions. Our results are indicative of a rapid increase in the galaxy IR comoving volume emissivity up to z ~ 1 and a constant average emissivity at z > 1. We also appear to measure a difference in the evolutionary rate of the source number densities as a function of luminosity, which is consistent with the downsizing evolutionary patterns reported for other samples of cosmic sources.

207 citations

Journal ArticleDOI
TL;DR: In this paper, a complete view of star formation from the local Universe to z − 4 was provided, using assumptions on earlier star formation history, and compared this evolution with previously published data in an attempt to draw a homogeneous picture of the global evolution of star creation in galaxies.
Abstract: Using new homogeneous luminosity functions (LFs) in the far-ultraviolet (FUV) from VVDS and in the far-infrared (FIR) from Herschel/PEP and Herschel/HerMES, we studied the evolution of the dust attenuation with redshift. With this information, we were able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density (SFRDTOT). By integrating SFRDTOT, we followed the mass building and analyzed the redshift evolution of the stellar mass density (SMD). This article aims at providing a complete view of star formation from the local Universe to z ~ 4 and, using assumptions on earlier star formation history, compares this evolution with previously published data in an attempt to draw a homogeneous picture of the global evolution of star formation in galaxies. Our main conclusions are that: 1) the dust attenuation AFUV is found to increase from z = 0 to z ~ 1.2 and then starts to decrease until our last data point at z = 3.6; 2) the estimated SFRD confirms published results to z ~ 2. At z > 2, we observe either a plateau or a small increase up to z ~ 3 and then a likely decrease up to z = 3.6; 3) the peak of AFUV is delayed with respect to the plateau of SFRDTOT and a probable origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) using assumptions (exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to zform = 10, we integrated SFRDTOT and obtained a good agreement with the published SMDs.

198 citations

Journal ArticleDOI
TL;DR: In this article, a complete view of star formation from the local universe to z = 4 was provided, using assumptions on earlier star formation history, comparing this evolution to what was known before in an attempt to draw a homogeneous picture of the global evolution of star creation in galaxies.
Abstract: Using new homogeneous LFs in the FUV and in the FIR Herschel/PEP and Herschel/HerMES, we study the evolution of the dust attenuation with redshift. With this information in hand, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density SFRD_TOT. By integrating SFRD_TOT, we follow the mass building and analyze the redshift evolution of the stellar mass density (SMD). This letter aims at providing a complete view of star formation from the local universe to z = 4 and, using assumptions on earlier star formation history, compares this evolution to what was known before in an attempt to draw a homogeneous picture of the global evolution of star formation in galaxies. The main conclusions of this letter are: 1) the dust attenuation A_FUV is found to increase from z = 0 to z \sim 1.2 and then starts to decrease up to our last data point at z = 3.6; 2) the estimated SFRD confirms published results up to z = 2. At z > 2, we observe either a plateau or a small increase up to z = 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT and a likely origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) using assumptions (namely exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to z_form = 10, we integrate SFRD_TOT and find a good agreement with the published SMDs.

187 citations


"Infrared luminosity functions based..." refers result in this paper

  • ...This may be qualitatively consistent with previous reports by Herschel that the dust attenuation peaks and declines at z>1.2 (Gruppioni et al. 2013; Burgarella et al. 2013)....

    [...]

Journal ArticleDOI
TL;DR: AKARI (formerly ASTRO-F) is an infrared space telescope designed for an all-sky survey at 10-180 µm and deep pointed surveys of selected areas at 2- 180 µm as discussed by the authors, which will significantly advance our understanding of galaxy evolution, the structure formation of the Universe, the nature of the buried AGNs, and the cosmic infrared background.
Abstract: AKARI (formerly ASTRO-F) is an infrared space telescope designed for an all-sky survey at 10-180 µm, and deep pointed surveys of selected areas at 2-180 µm. The deep pointed surveys with AKARI will significantly advance our understanding of galaxy evolution, the structure formation of the Universe, the nature of the buried AGNs, and the cosmic infrared background. Here we describe the important characteristics of the AKARI mission: the orbit, and the attitude control system, and investigate the optimum survey area based on the updated pre-flight sensitivities of

158 citations


"Infrared luminosity functions based..." refers methods in this paper

  • ...The L18W flux (Matsuhara et al. 2006) are used to apply the 1/Vmax method, because it is a wide, sensitive filter (but using the L15 flux limit does not change our main results)....

    [...]

Journal ArticleDOI
TL;DR: In this article, rest-frame luminosity functions (LFs) at 3.6, 4.5, 5.8, 8 and 24 μm were constructed using optical and infrared (IR) data from the Spitzer Wide-area Infrared Extragalactic (SWIRE) Survey.
Abstract: We construct rest-frame luminosity functions (LFs) at 3.6, 4.5, 5.8, 8 and 24 μm over the redshift range 0 < z < 2 for galaxies and 0 < z < 4 for optical quasi-stellar objects (QSOs), using optical and infrared (IR) data from the Spitzer Wide-area Infrared Extragalactic (SWIRE) Survey. The 3.6- and 4.5-μm galaxy LFs show evidence for moderate positive luminosity evolution up to z~ 1.5, consistent with the passive ageing of evolved stellar populations. Their comoving luminosity density was found to evolve passively, gradually increasing out to z~ 0.5–1 but flattening, or even declining, at higher redshift. Conversely, the 24-μm galaxy LF, which is more sensitive to obscured star formation and/or active galactic nuclei (AGN) activity, undergoes strong positive evolution, with the derived IR energy density and star formation rate (SFR) density ∝ (1 +z)^γ with γ= 4.5^(+0.7)_(−0.6) and the majority of this evolution occurring since z~ 1. Optical QSOs, however, show positive luminosity evolution in all bands, out to the highest redshifts (3 < z < 4). Modelling as L^* ∝ (1 +z)^γ gave γ= 1.3^(+0.1)_(−0.1) at 3.6 μm, γ= 1.0^(+0.1)_(−0.1) at 4.5 μm and stronger evolution at the longer wavelengths (5.8, 8 and 24 μm), of γ~ 3. Comparison of the galaxy LFs to predictions from a semi-analytic model based on cold dark matter (CDM) indicates that an initial mass function (IMF) skewed towards higher mass star formation in bursts compared to locally be preferred. As a result, the currently inferred massive SFRs in distant submm sources may require substantial downwards revision.

145 citations


"Infrared luminosity functions based..." refers background in this paper

  • ...4.1 The 8µm LF We first present monochromatic 8µm LFs, because the 8µm luminosity (L8µm) has been known as a good indicator of the TIR luminosity (Babbedge et al. 2006; Huang et al. 2007; Goto et al. 2011a)....

    [...]

  • ...This is a common exercise with the limited depths of the current IR data (Babbedge et al. 2006; Caputi et al. 2007)....

    [...]

Related Papers (5)
Frequently Asked Questions (6)
Q1. What are the contributions mentioned in the paper "Infrared luminosity functions based on 18 mid-infrared bands: revealing cosmic star formation history with akari and hyper suprime-cam*" ?

Combined with the accurate mid-IR luminosity measurement, the authors constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. 

The L18W flux (Matsuhara et al. 2006) are used to apply the 1/Vmax method, because it is a wide, sensitive filter (but using the L15 flux limit does not change their main results). 

Mid-infrared (mid-IR) is one of the less explored wavelengths due to the earth’s atmosphere, and difficulties in developing sensitive detectors. 

Uncertainties of the LF values includefluctuations in the number of sources in each luminosity bin, the photometric redshift uncertainties, the k-correction uncertainties, and the flux errors. 

It has the largest FoV among optical cameras on 8m-class telescopes, and can cover the AKARI NEP wide field (5.4 deg2) with only 4 FoV (Fig.1). 

even with AKARI’s sensitivity, the observation might not be deep enough to reliably measure the faint-end slope of 12µm LFs, possibly because 12µm does not contain as luminous emission lines as in the case of 8µm.