03.08.2013 Views

Copyright by Athena Ranice Stacy 2011 - The University of Texas at ...

Copyright by Athena Ranice Stacy 2011 - The University of Texas at ...

Copyright by Athena Ranice Stacy 2011 - The University of Texas at ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

where the sum now extends from the I-front position <strong>at</strong> the previous timestep<br />

t0 to its position <strong>at</strong> the current timestep t0 + ∆t. <strong>The</strong> above-described ray-<br />

tracing is performed separ<strong>at</strong>ely for the H ii and He iii regions.<br />

3.2.5 Photoioniz<strong>at</strong>ion and He<strong>at</strong>ing<br />

Particles determined to be within the extent <strong>of</strong> the H ii and He iii re-<br />

gions are given additional ioniz<strong>at</strong>ion and he<strong>at</strong>ing r<strong>at</strong>es in the chemistry solver:<br />

and<br />

Γ = nn<br />

kion =<br />

∞<br />

νmin<br />

∞<br />

νmin<br />

Fνσν<br />

<br />

Fνσν<br />

dν (3.5)<br />

hν<br />

1 − νmin<br />

ν<br />

<br />

dν , (3.6)<br />

where Fν is the incoming specific flux and σν the photo-ioniz<strong>at</strong>ion cross section.<br />

For a blackbody, we have<br />

Fν =<br />

where r is the distance from the sink.<br />

<strong>by</strong><br />

L∗<br />

4σT 4 eff r2 Bν , (3.7)<br />

Finally, H2 dissoci<strong>at</strong>ion <strong>by</strong> LW radi<strong>at</strong>ion (11.2 to 13.6 eV) is described<br />

kH2 = 1.1 × 10 8 fshield FLW s −1<br />

(3.8)<br />

(Abel et al. 1997), where FLW denotes the radi<strong>at</strong>ion flux, in units <strong>of</strong> erg s −1<br />

cm −2 Hz −1 , <strong>at</strong> hν = 12.87eV, and fshield is the factor <strong>by</strong> which H2 self-shielding<br />

63

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!