Model of the Evolution of Infrared Galaxies
G. Lagache, H. Dole, J-L. Puget




http://lully.as.arizona.edu/Model



Page Content
1- Background of the Model Lagache, Dole, Puget, 2003.
2- The Spitzer Update: the Updated Model Lagache, Dole, Puget, et al, 2004.
3- IDL programs and README for making Predictions
4- IDL save sets containing data like dN/(dlnLdz), dS/dz
5- SED (Spectral Energy Distributions) used in the model, in ASCII tables
6- Source Counts from 7 microns to 2.1 mm in ASCII tables
7- Plots of Source Counts in png
8- Getting the Papers
9- Links to other models on the market

1- Background

Modelling the infrared galaxy evolution using a phenomenological approach

G. Lagache, H. Dole, J-L. Puget, MNRAS, 2003, 338, 555

To characterise the cosmological evolution of the sources contributing to the infrared extragalactic background, we have developped a phenomenological model that constrains in a simple way the galaxy luminosity function evolution with the redshift, and fits all the existing source counts and redshift distributions, Cosmic Infrared Background intensity and fluctuations observations from the mid-infrared to the submillimetre range. The model is based on template spectra of starburst and normal galaxies, and on the local infrared luminosity function. Although the Cosmic Infrared Background can be modeled with very different luminosity functions as long as the radiation production with redshift is the right one, the number counts, and the anisotropies of the unresolved background, imply that the luminosity function must change dramatically with redshift, with a rapid evolution of the high-luminosity sources (Log(L) > 11.5 Lo) from z=0 to z=1 which then stay rather constant up to redshift 5. The derived evolution of the IR luminosity function may be linked to a bimodal star formation process, one associated with the quiescent and passive phase of the galaxy evolution and one associated with the starburst phase, trigerred by merging and interactions. The latter dominates the infrared and submillimetre ouput energy of the Universe.
The model is intended as a convenient tool to plan further observations, as illustrated through predictions for Herschel, Planck and ALMA observations. Our model predictions for given wavelengths, together with some usefull routines, are available for general use hereafter.

The paper is published in MNRAS, 2003, vol 338, Number 3, 555.

In all the work, we use h=0.65, Omega_Lambda = 0.7, Omega_0 = 0.3


2- The Spitzer Update

PAHs contribution to the infrared output energy of the Universe at z~2

G. Lagache, H. Dole, J-L. Puget, P. G. Perez-Gonzalez, E. Le Floc'h, G. H. Rieke, C. Papovich, E. Egami, A. Alonso-Herrero, C. W. Engelbracht, K.D. Gordon, K. A. Misselt, and J. E. Morrison, ApJSS, 2004, in press

We present an updated phenomenological galaxy evolution model to fit the Spitzer 24, 70 and 160 micron number counts as well as all the previous mid and far infrared observations. Only a minor change of the co-moving luminosity density distribution in the previous model, combined with a slight modification of the starburst template spectra mainly between 12 and 30 microns, are required to fit all the data available. We show that the peak in the MIPS 24 microns counts is dominated by galaxies with redshift between 1 and 2, with a non negligible contribution from the z>2 galaxies (~30% at S=0.2 mJy). The very close agreement between the model and number counts at 15 and 24 microns strikingly implies that
(1) the PAHs (Policyclic Aromatic Hydrocarbons) features remain prominent in the redshift band 0.5 to 2.5 and
(2) the IR energy output has to be dominated by ~3 10^11 Lo to ~10^12 Lo galaxies from redshift 0.5 to 2.5.
Combining Spitzer with the Infrared Space Observatory (ISO) deep cosmological surveys gives for the first time an unbiased view of the infrared Universe from z=0 to z=2.5.

The paper is in press in ApJS (Spitzer Special Issue), 2004


3- IDL Programs

The README file

make_integral_counts.pro: Compute the integral and differential number counts

make_z_distrib.pro: Compute the redshift distribution for sources above a given flux

compute_fluc.pro: Compute the fluctuations level made by sources below a given flux

example.scr: Some examples

loadct_plot.pro: A useful program for displaying colors

allpro.tgz : All these files in a single gzipped tar file (extraction: tar -xvzf allpro.tgz)


4- IDL Save Sets

The IDL save files contain the following variables:

- dNdLNLdZ: dN / (dlnLdz) in function of L and z
- z : redshift
- lum_array: array of luminosities in solar luminosity
- DSDZ : dS/dz in function of L and z
- Slz: flux in Jy in function of L and z

The "generic filters" are top-hat filters with delta_lambda / lambda = 1/3. Otherwise, we used the real transmission filters.

ISOCAM 15 micron starburst, ISOCAM 15 micron filter
ISOCAM 15 micron normal galaxies, ISOCAM 15 micron filter

ISOPHOT and PACS 170 micron starburst, ISOPHOT micron 170 filter
ISOPHOT and PACS 170 micron normal galaxies, ISOPHOT 170 micron filter
ISOPHOT 90 micron starburst, ISOPHOT 90 micron filter
ISOPHOT 90 micron normal galaxies, ISOPHOT 90 micron filter

SPITZER/MIPS 24 micron starburst, MIPS 24 micron filter
SPITZER/MIPS 24 micron normal galaxies, MIPS 24 micron filter

SPITZER/MIPS 70 micron starburst, MIPS 70 micron filter
SPITZER/MIPS 70 micron normal galaxies, MIPS 70 micron filter

SPITZER/MIPS 160 micron starburst, MIPS 160 micron filter
SPITZER/MIPS 160 micron normal galaxies, MIPS 160 micron filter

IRAS 60 micron starburst, COBE 60 micron filter
IRAS 60 micron normal galaxies,COBE 60 micron filter

IRAS 100 micron starburst, COBE 100 micron filter
IRAS 100 micron normal galaxies, COBE 100 micron filter

PACS 75 micron starburst, 75 micron generic filter
PACS 75 micron normal galaxies, 75 micron generic filter

PACS 110 micron starburst, 110 micron generic filter
PACS 110 micron normal galaxies, 110 micron generic filter

SPIRE 250 micron starburst 250 micron generic filter
SPIRE 250 micron normal galaxies, 250 micron generic filter

SPIRE and PLANCK 350 micron starburst, 350 micron generic filter
SPIRE and PLANCK 350 micron normal galaxies, 350 micron generic filter

SPIRE and PLANCK 550 micron starburst, 550 micron generic filter
SPIRE and PLANCK 550 micron normal galaxies,550 micron generic filter

PLANCK 850 micron starburst, 850 micron generic filter
PLANCK 850 micron normal galaxies, 850 micron generic filter

SCUBA 850 micron starburst, SCUBA 850 micron filter
SCUBA 850 micron normal galaxies,SCUBA 850 micron filter

MAMBO 1300 micron starburst, 1300 micron generic filter
MAMBO 1300 micron normal galaxies, 1300 micron generic filter

PLANCK 1380 micron starburst, 1380 micron generic filter
PLANCK 1380 micron normal galaxies,1380 micron generic filter

PLANCK 2097 micron starburst, 2097 micron generic filter
PLANCK 2097 micron normal galaxies,2097 generic filter


Luminosity function of the starburst galaxies
Luminosity function of the normal galaxies


5- Spectral Energy Distributions in ASCII


Note that these SED are valid to wavelengths from 4 microns to 3 mm ONLY !

Starburst SED ASCII table
1st column: wavelength in m
2nd column: SED in W/Hz of a 10^9 Lsun galaxy
3rd column: SED in W/Hz of a 10^10 Lsun galaxy
4th column: SED in W/Hz of a 10^11 Lsun galaxy
5th column: SED in W/Hz of a 10^12 Lsun galaxy
6th column: SED in W/Hz of a 10^13 Lsun galaxy

Normal/Cold SED ASCII Table
1st column: wavelength in m
2nd column: SED in W/Hz of a 10^11 Lsun galaxy


Many more Starburst SED (with LogL steps of 0.1 Lsun) in ascii gzipped
Many more Normal/Cold SED (with LogL steps of 0.1 Lsun) in ascii gzipped


6- Source Counts in ASCII

Each table contains 4 columns, in ascii:

- Log(flux[Jy])
- integral source counts N>S in sr^{-1}
- differential source counts dN/dS in sr^{-1} Jy^{-1}
- the contribution to the CIB intensity at that particular wavelength in %

Needless to say, the source counts and the contribution to the CIB are the TOTAL contribution of our components (ie cold AND starbust galaxies).

3.6 microns Spitzer-IRAC
4.5 microns Spitzer-IRAC
5.8 microns Spitzer-IRAC
7 microns ISOCAM
8.0 microns Spitzer-IRAC
15 microns ISOCAM
24 microns Spitzer-MIPS
60 microns IRAS COBE-DIRBE
70 microns Spitzer-MIPS
75 microns HERSCHEL-PACS (generic filter)
90 microns ISOPHOT
100 microns IRAS COBE-DIRBE
110 microns HERSCHEL-PACS (generic filter)
130 microns (generic filter)
140 microns COBE-DIRBE
160 microns Spitzer-MIPS
170 microns ISOPHOT
180 microns (generic filter)
200 microns (generic filter)
250 microns (generic filter)
350 microns HERSCHEL-SPIRE & PLANCK-HFI (generic filter)
550 microns HERSCHEL-SPIRE & PLANCK-HFI (generic filter)
750 microns (generic filter)
849 microns HERSCHEL-SPIRE & PLANCK-HFI (generic filter)
850 microns SCUBA(generic filter)
1380 microns PLANCK-HFI (generic filter)
2097 microns PLANCK-HFI (generic filter)


7- Plots of Source Counts

Integral Source Counts

Differential Source Counts


8- Getting the Papers

- Lagache, Dole, Puget, 2003, MNRAS, 338, 555 in PDF

- Lagache, Dole, Puget et al, 2004, ApJS, in PDF


9- Links to Other Models

- Bruno Guiderdoni et al (e.g astro-ph/9710340)
- Alex King and Michael Rowan-Robinson (e.g. MNRAS 2001)
- Chris Pearson (e.g astro-ph/0011335)
- Kevin Xu, Carol Lonsdale et al (e.g ApJ 2001, 2003)
- Matthew Malkan & Floyd Stecker (e.g astro-ph/0009500)

- unfortunately, other popular models like Chary and Elbaz (2001), Devriendt et al (2000), Balland et al (2001), Franceschini et al (2001) have (to our knowledge) no web page to allow a public download of the models. It is a shame, but fortunately most of authors kindly share their models on request !



Guilaine Lagache and Hervé Dole.
Last modified: Fri Oct 29 13:54:28 CEST 2004

Valid XHTML 1.0!