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.

educes the LW dissoci<strong>at</strong>ion r<strong>at</strong>e. This self-shielding factor depends upon the<br />

H2 column density NH2. With the above ray-tracing scheme, we determine<br />

NH2 along each ray <strong>by</strong> summing the contribution from each bin. We then<br />

use results from Draine and Bertoldi (1996) to determine the value for fshield.<br />

We note a recent upd<strong>at</strong>e to their fshield fitting formula (Wolcott-Green et al.<br />

<strong>2011</strong>; Wolcott-Green and Haiman <strong>2011</strong>), but do not expect this to significantly<br />

affect the results for our particular case. Because <strong>of</strong> the unusually high column<br />

densities within the molecular disk (NH2 > ∼ 10 26 cm −2 ), we calcul<strong>at</strong>e fshield<br />

using equ<strong>at</strong>ion 37 from Draine and Bertoldi (1996), which is more accur<strong>at</strong>e for<br />

large NH2 than their simple power-law expression in their equ<strong>at</strong>ion 36.<br />

<strong>The</strong>se he<strong>at</strong>ing, ioniz<strong>at</strong>ion, and dissoci<strong>at</strong>ion r<strong>at</strong>es are accounted for<br />

<strong>at</strong> every timestep, while the ray-tracing procedure is performed every fifth<br />

timestep. <strong>The</strong> r<strong>at</strong>es are upd<strong>at</strong>ed every 10 yr as they evolve with the proto-<br />

stellar mass and accretion r<strong>at</strong>e, because these determine the values <strong>of</strong> L∗ and<br />

Teff.<br />

3.2.6 Protostellar evolution model<br />

Our ray-tracing algorithm requires an input <strong>of</strong> protostellar effective<br />

temper<strong>at</strong>ure Teff and luminosity L∗. This is the sum <strong>of</strong> Lacc, the accretion<br />

luminosity, and Lphoto, the luminosity gener<strong>at</strong>ed <strong>at</strong> the protostellar surface:<br />

L∗ = Lacc + Lphoto =<br />

GM∗ ˙M<br />

R∗<br />

+ Lphoto, (3.9)<br />

where M∗ is the protostellar mass, ˙M the accretion r<strong>at</strong>e, and R∗ the proto-<br />

stellar radius. <strong>The</strong> photospheric luminosity is either due to Kelvin Helmholtz<br />

64

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

Saved successfully!

Ooh no, something went wrong!