Project

General

Profile

Actions

Feature #2

closed

Add intergalactic extinction in stuff

Added by Emmanuel Bertin over 8 years ago. Updated almost 8 years ago.

Status:
Closed
Priority:
Normal
Start date:
05/27/2013
Due date:
05/28/2013
% Done:

100%

Estimated time:
3.00 h
Spent time:

Files

igmabs.pdf (42.6 KB) igmabs.pdf IGM transmission curves computed for various redshifts with the new code on top of MEGACAM filter responses Emmanuel Bertin, 05/28/2013 07:16 PM
igmabs.png (81.6 KB) igmabs.png IGM transmission curves computed for various redshifts with the new code on top of MEGACAM filter responses Emmanuel Bertin, 05/28/2013 07:16 PM
Actions #2

Updated by Emmanuel Bertin over 8 years ago

We use the average model from Madau 1995 and Madau et al. 1996 for both Lyman absorbants:

(...) For $\lambda_{\beta}(1+z_{\rm em}) < \lambda_{\rm obs} < \lambda_{\alpha} (1+z_{\rm em})$, where $\lambda_{\alpha}=1216\,\mathring{A}$ and $\lambda_{\beta}=1026\,\mathring{A}$, a galaxy's continuum intensity is attenuated by the combined blanketing effect of many Ly$\alpha$ forest absorption lines, with effective opacity
$$ \tau_{\rm eff}=0.0036\left({\lambda_{\rm obs}\over \lambda_{\alpha}}\right) ^{3.46},$$
(Press, Rybicki, & Schneider 1993). Hence, line blanketing from Ly$\alpha$ alone will produce $\ge1\,$mag of attenuation in the continuum spectrum shortward of $6000\,\mathring{A}$ of a galaxy at $z_{\rm em}\ge 4$.

When $\lambda_{\rm obs}\le\lambda_{\beta}(1+z_{\rm em})$, a significant contribution to the blanketing opacity comes from the higher order lines of the Lyman series. A standard curve of growth analysis has been applied to numerically compute the attenuation expected from line blanketing of Ly$\beta$, $\gamma$, $\delta$ plus 13 higher order members. In the wavelength range $\lambda_{i+1}(1+z_{\rm em})<\lambda_{\rm obs}<\lambda_i(1+z_{\rm em})$, the total optical depth can be written as the sum of the contributions from the $j\rightarrow 1$ transitions,
$$ \tau_{\rm eff}=\sum_{j=2,i} A_j({\lambda_{\rm obs}\over \lambda_j})^{3.46}, $$
where $A_j=(1.7\times 10^{-3}, 1.2\times 10^{-3},9.3\times 10^{-4})$, and $\lambda_j=(1026, 973, 950\,\mathring{A})$ for Ly$\beta$, $\gamma$, and $\delta$, respectively. It is worth noting that heavy element absorbers make a negligible contribution to the blanketing optical depth at high redshifts, as this is dominated by those lines which lie at the transition between the linear and the flat part of the curve of growth, i.e., with $N_{\rm HI} \approx 10^{13.6}\,{\rm cm}^{-2}$ in the case of Ly$\alpha$.
Continuum absorption by HI affects photons observed at $\lambda_{\rm obs}<\lambda_L(1+z_{\rm em})$, where $\lambda_L=912\,\mathring{A}$ is the Lyman limit. An approximate (within 5%) integration of equation yields for the effective photoelectric optical depth along the line of sight:
$$\tau_{eff}=0.25x_c^3 (x_{em}^{0.46}-x_c^{0.46})+9.4x_c^{1.5} (x_{em}^{0.18}-x_c^{0.18})-0.7x_c^3(x_c^{-1.32}-x_{em}^{-1.32})-0.023 (x_{em}^{1.68}-x_c^{1.68})$$
(...)

and for metal lines (although Madau 1995 claims that their effect is negligible for $2.5\ge z\ge4.5$:

$$ \tau_{\rm eff}=0.0017\left({\lambda_{\rm obs}\over \lambda_{\alpha}}\right) ^{1.68}$$

The complete list of $\lambda_j$ and $A_j$ is taken from the xspec IGM attenuation module by Martin Still and Frank Marshall:

$\lambda_j$ (Å) 1215.67 1025.72 972.537 949.743 937.803 930.748 926.226 923.150
$A_j$ 0.0036 0.0017 0.0011846 0.0009410 0.0007960 0.0006967 0.0006236 0.0005665
$\lambda_j$ (Å) 920.963 919.352 918.129 917.181 916.429 915.824 915.329 914.919 914.576
$A_j$ 0.0005200 0.0004817 0.0004487 0.0004200 0.0003947 0.000372 0.0003520 0.0003334 0.00031644

We obtain the following transmission curves for sources at $z=1,2,3,4,5,6$, plotted over MEGACAM 's global filter responses:
IGM transmission curves computed for various redshifts with the new code on top of MEGACAM filter responses

Actions #3

Updated by Emmanuel Bertin over 8 years ago

  • % Done changed from 0 to 80
Actions #4

Updated by Emmanuel Bertin over 8 years ago

Actions #5

Updated by Emmanuel Bertin over 8 years ago

  • Status changed from New to In Progress
  • % Done changed from 80 to 90

New IGM_TYPE MADAU_AVERAGE option implemented and tested. Commit on the public SVN pending because of data policy issues with public release (this commit includes new DECam passband data that require authorization from the original authors/providers).

Actions #6

Updated by Emmanuel Bertin over 8 years ago

New "stuff":http://astromatic.net/software/stuff 1.26.0 with IGM extinction option (turned on by default) installed on the morpho cluster.

Actions #7

Updated by Emmanuel Bertin over 8 years ago

  • Subject changed from Add intergalactic extinction to Add intergalactic extinction in stuff
  • Private changed from No to Yes
Actions #8

Updated by Emmanuel Bertin over 8 years ago

  • Private changed from Yes to No
Actions #9

Updated by Emmanuel Bertin over 8 years ago

  • Status changed from In Progress to Resolved
  • % Done changed from 90 to 100

Not committed yet but installed on morpho and running.

Actions #10

Updated by Emmanuel Bertin over 8 years ago

  • Status changed from Resolved to Closed
Actions #11

Updated by Emmanuel Bertin almost 8 years ago

  • Project changed from 24 to Stuff
Actions #12

Updated by Emmanuel Bertin almost 8 years ago

  • Parent task set to #298
Actions

Also available in: Atom PDF