{"id":45,"date":"2013-05-03T22:20:08","date_gmt":"2013-05-03T22:20:08","guid":{"rendered":"http:\/\/portillo.ca\/nebula\/?page_id=45"},"modified":"2014-10-21T22:10:03","modified_gmt":"2014-10-21T22:10:03","slug":"modelling-nebular-structure","status":"publish","type":"page","link":"http:\/\/portillo.ca\/nebula\/the-details\/modelling-nebular-structure\/","title":{"rendered":"Modelling Nebular Structure"},"content":{"rendered":"<p>The Pillars of Creation is a false colour image with SII 6176\/6731 \u00c5 in red, H\u03b1 6536 \u00c5 in green, and OIII 5007 \u00c5 in blue. To determine the properties of the photoevaporative flow, Hester et al. used a spectral synthesis code (Cloudy 85) to predict the emission from these lines from different models of the nebula. The interactive <a title=\"Nebula Predictor\" href=\"http:\/\/portillo.ca\/nebula\/nebula-predictor\/\">nebula predictor<\/a> shows the output of a grid of <a title=\"Cloudy &amp; Associates\" href=\"http:\/\/www.nublado.org\/\">Cloudy 13<\/a> models.<\/p>\n<p>An important part of predicting the spectrum is the radiative transfer equation, with the incoming radiation determined by the star:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7BdI_%7B%5Cnu%7D%7D%7Bds%7D+%3D+-I_%7B%5Cnu%7D+%5Cfrac%7Bd%5Ctau_%7B%5Cnu%7D%7D%7Bds%7D+%2B+j_%7B%5Cnu%7D++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"&#92;frac{dI_{&#92;nu}}{ds} = -I_{&#92;nu} &#92;frac{d&#92;tau_{&#92;nu}}{ds} + j_{&#92;nu}  \" class=\"latex\" \/><\/p>\n<p>Photoionization absorbs radiation:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7Bd%5Ctau_%7B%5Cnu%7D%7D%7Bds%7D+%3D+%5Csum%5Climits_%7BX%2Ci%7D+n%28X%5E%7B%2Bi%7D%29a_%7B%5Cnu%7D%28X%5E%7B%2Bi%7D%29+&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"&#92;frac{d&#92;tau_{&#92;nu}}{ds} = &#92;sum&#92;limits_{X,i} n(X^{+i})a_{&#92;nu}(X^{+i}) \" class=\"latex\" \/><\/p>\n<p>Recombination provides emissivity:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=j_%7B%5Cnu%7D%28X%5E%7B%2Bi%7D%2CT%29+%3D+%5Cfrac%7B2h%5Cnu%5E%7B3%7D%7D%7Bc%5E%7B2%7D%7D+%28%5Cfrac%7Bh%5E%7B2%7D%7D%7B2%5Cpi+mkT%7D%29%5E%7B3%2F2%7Da_%7B%5Cnu%7D%5Cexp%28-%5Cfrac%7Bh%28%5Cnu-%5Cnu_%7B0%7D%29%7D%7BkT%7D%29n%28X%5E%7B%2Bi%7D%29n%28e%29+&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"j_{&#92;nu}(X^{+i},T) = &#92;frac{2h&#92;nu^{3}}{c^{2}} (&#92;frac{h^{2}}{2&#92;pi mkT})^{3\/2}a_{&#92;nu}&#92;exp(-&#92;frac{h(&#92;nu-&#92;nu_{0})}{kT})n(X^{+i})n(e) \" class=\"latex\" \/><\/p>\n<p>In practise, only hydrogen and helium need to be included in radiative transfer because of their overwhelming abundances.<\/p>\n<p>To determine the abundances of the different ionization states of each species, ionization equilibrium is used. In any place in the nebula for each ionization state, ionizations (left hand side) are balanced by recombinations (right hand side):<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=n%28X%5E%7B%2Bi%7D%29%5Cint%5Climits_%7B%5Cnu_%7B0%7D%7D%5E%7B%5Cinfty%7D%5Cfrac%7B4%5Cpi+J_%7B%5Cnu%7D%7D%7Bh%5Cnu%7Da_%7B%5Cnu%7D%28X%5E%7B%2Bi%7D%29d%5Cnu+%3D+n%28X%5E%7B%2B%28i%2B1%29%7D%29n%28e%29%5Calpha_%7BG%7D%28X%5E%7B%2Bi%7D%2CT%29++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"n(X^{+i})&#92;int&#92;limits_{&#92;nu_{0}}^{&#92;infty}&#92;frac{4&#92;pi J_{&#92;nu}}{h&#92;nu}a_{&#92;nu}(X^{+i})d&#92;nu = n(X^{+(i+1)})n(e)&#92;alpha_{G}(X^{+i},T)  \" class=\"latex\" \/><\/p>\n<p>The temperature is determined by thermal equilibrium for the gas:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=G+-+L_%7BR%7D+%3D+L_%7BFF%7D+%2B+L_%7BC%7D++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"G - L_{R} = L_{FF} + L_{C}  \" class=\"latex\" \/><\/p>\n<p>With <em>G<\/em> being energy input from photoionization, <em>L<sub>R<\/sub><\/em> being energy losses from recombination, <em>L<sub>FF<\/sub><\/em> being energy losses from free-free radiation, and <em>L<sub>R<\/sub><\/em> being energy losses to collisionally excited lines. Each of these terms as a sum over all species and ionization states.<\/p>\n<p>When a photoionization occurs, the excess energy <em>h<\/em>(<em>\u03bd<\/em>&#8211;<em>\u03bd<sub>0<\/sub><\/em>) of the photon comprises the kinetic energy of the freed electron:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=G%28X%5E%7B%2Bi%7D%29+%3D+n%28X%5E%7B%2Bi%7D%29%5Cint%5Climits_%7B%5Cnu_%7B0%7D%7D%5E%7B%5Cinfty%7D%5Cfrac%7B4%5Cpi+J_%7B%5Cnu%7D%7D%7Bh%5Cnu%7D+h%28%5Cnu-%5Cnu_%7B0%7D%29+a_%7B%5Cnu%7D%28X%5E%7B%2Bi%7D%29d%5Cnu++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"G(X^{+i}) = n(X^{+i})&#92;int&#92;limits_{&#92;nu_{0}}^{&#92;infty}&#92;frac{4&#92;pi J_{&#92;nu}}{h&#92;nu} h(&#92;nu-&#92;nu_{0}) a_{&#92;nu}(X^{+i})d&#92;nu  \" class=\"latex\" \/><\/p>\n<p>When a recombination occurs, the kinetic energy of the electron 1\/2 <em>m<\/em> <em>v<\/em><sup>2<\/sup> is lost from the gas:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=L_%7BR%7D%28X%5E%7B%2Bi%7D%29%3Dn%28X%5E%7B%2B%28i%2B1%29%7D%29n%28e%29%5Cint%5Climits_%7B0%7D%5E%7B%5Cinfty%7D+%5Csigma%28X%5E%7B%2Bi%7D%2Cv%29v+%5Cfrac%7B1%7D%7B2%7Dmv%5E%7B2%7Df%28v%29dv++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"L_{R}(X^{+i})=n(X^{+(i+1)})n(e)&#92;int&#92;limits_{0}^{&#92;infty} &#92;sigma(X^{+i},v)v &#92;frac{1}{2}mv^{2}f(v)dv  \" class=\"latex\" \/><\/p>\n<p><em>L<sub>R<\/sub><\/em> &lt; <em>G<\/em> because slower electrons are preferentially captured.<\/p>\n<p>Electrons and ions will scatter off each other, releasing some radiation:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=L_%7BFF%7D%28X%5E%7B%2Bi%7D%29+%3D+%5Cfrac%7B32%5Cpi+e%5E%7B6%7D+i%5E%7B2%7D%7D%7B3%5E%7B3%2F2%7Dhmc%5E%7B3%7D%7D+%28%5Cfrac%7B2%5Cpi+k+T%7D%7Bm%7D%29%5E%7B1%2F2%7Dg_%7BFF%7Dn%28e%29n%28X%5E%7B%2B%28i%2B1%29%7D%29++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"L_{FF}(X^{+i}) = &#92;frac{32&#92;pi e^{6} i^{2}}{3^{3\/2}hmc^{3}} (&#92;frac{2&#92;pi k T}{m})^{1\/2}g_{FF}n(e)n(X^{+(i+1)})  \" class=\"latex\" \/><\/p>\n<p>In the above three mechanisms, only hydrogen and helium are important because of their abundances.<\/p>\n<p>Metals, however, have many collisionally excited lines and thus contribute to <em>L<sub>C<\/sub><\/em>. The atoms and ions in the gas can be excited or de-excited by radiation (<em>A<\/em> terms) and collisions (<em>q<\/em> terms). Transitions from excitation state <em>j<\/em> to <em>k<\/em> (left hand side) must be balanced by transitions from excitation state <em>k<\/em> back to <em>j<\/em> (right hand side), determining the population <em>n<sub>j<\/sub><\/em> of each state of species <em>X<sup>+i<\/sup><\/em>:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=%5Csum%5Climits_%7Bk%7Dn_%7Bk%7Dn%28e%29q_%7Bkj%7D+%2B+%5Csum%5Climits_%7Bk%7Dn_%7Bk%7DA_%7Bkj%7D+%3D+%5Csum%5Climits_%7Bk%7Dn_%7Bj%7Dn%28e%29q_%7Bjk%7D+%2B+%5Csum%5Climits_%7Bk%7Dn_%7Bj%7DA_%7Bjk%7D++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"&#92;sum&#92;limits_{k}n_{k}n(e)q_{kj} + &#92;sum&#92;limits_{k}n_{k}A_{kj} = &#92;sum&#92;limits_{k}n_{j}n(e)q_{jk} + &#92;sum&#92;limits_{k}n_{j}A_{jk}  \" class=\"latex\" \/><\/p>\n<p>Then, spontaneous de-excitations from each state determine the cooling rate:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=L_%7BC%7D%28X%5E%7B%2Bi%7D%29+%3D+%5Csum%5Climits_%7Bj%7Dn_%7Bj%7D%5Csum%5Climits_%7Bk%3Cj%7DA_%7Bjk%7Dh%5Cnu_%7Bjk%7D++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"L_{C}(X^{+i}) = &#92;sum&#92;limits_{j}n_{j}&#92;sum&#92;limits_{k&lt;j}A_{jk}h&#92;nu_{jk}  \" class=\"latex\" \/><\/p>\n<p>These equations must be solved iteratively. To start, the \u201con-the-spot\u201d approximation can be made, which assumes that photons coming from recombination are immediately absorbed. Thus, the radiation intensity at any point is the sum of intensity from the star (first term) and from local recombinations:<\/p>\n<p><img decoding=\"async\" src=\"http:\/\/s0.wp.com\/latex.php?latex=J_%7B%5Cnu%7D+%3D+%5Cfrac%7B%5Cpi+F_%7B%5Cnu%2C%5Cstar%7D+R%5E%7B2%7D%7D%7Br%5E%7B2%7D%7D+%5Cexp%28-%5Ctau_%7B%5Cnu%7D%29+%2B+%5Cfrac%7Bj_%7B%5Cnu%7D%7D%7B%5Csum%5Climits_%7BX%2Ci%7D+n%28X%5E%7B%2Bi%7D%29+a_%7B%5Cnu%7D+%28X%5E%7B%2Bi%7D%29%7D++&#038;bg=0f0f0f&#038;fg=ffffff&#038;s=3&#038;c=20201002\" alt=\"J_{&#92;nu} = &#92;frac{&#92;pi F_{&#92;nu,&#92;star} R^{2}}{r^{2}} &#92;exp(-&#92;tau_{&#92;nu}) + &#92;frac{j_{&#92;nu}}{&#92;sum&#92;limits_{X,i} n(X^{+i}) a_{&#92;nu} (X^{+i})}  \" class=\"latex\" \/><\/p>\n<p>Then the ionization and thermal equilibrium equations can be integrated outwards from the star to get approximate ionization populations <em>n<\/em>(<em>X<sup>+i<\/sup><\/em>) and temperature <em>T<\/em> at each point in the nebula. Then the radiative transfer equation can be integrated using these approximate solutions, and then the ionization and thermal equilibrium conditions re-integrated using the new radiative transfer solution for <em>J<sub>\u03bd<\/sub><\/em>. This procedure can be iterated until the solution converges.<\/p>\n<p>(For more details on modelling nebular structure, consult Osterbrock &amp; Ferland&#8217;s <a title=\"Astrophysics of Gaseous Nebulae and Active Galactic Nuclei\" href=\"http:\/\/uscibooks.com\/ost2.htm\">Astrophysics of Gaseous Nebulae and Active Galactic Nuclei<\/a>.)<\/p>\n<p>Now take <a title=\"A Closer Look\" href=\"http:\/\/portillo.ca\/nebula\/the-details\/a-closer-look\/\">a closer look<\/a> at &#8220;The Pillars of Creation&#8221; to see how these models of nebular structure help us understand the properties of the Eagle Nebula.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>The Pillars of Creation is a false colour image with SII 6176\/6731 \u00c5 in red, H\u03b1 6536 \u00c5 in green, and OIII 5007 \u00c5 in blue. To determine the properties of the photoevaporative flow, Hester et al. used a spectral &hellip; <a href=\"http:\/\/portillo.ca\/nebula\/the-details\/modelling-nebular-structure\/\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"parent":13,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"jetpack_post_was_ever_published":false,"footnotes":""},"class_list":["post-45","page","type-page","status-publish","hentry"],"jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/P3uAvH-J","_links":{"self":[{"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/pages\/45","targetHints":{"allow":["GET"]}}],"collection":[{"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/comments?post=45"}],"version-history":[{"count":28,"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/pages\/45\/revisions"}],"predecessor-version":[{"id":249,"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/pages\/45\/revisions\/249"}],"up":[{"embeddable":true,"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/pages\/13"}],"wp:attachment":[{"href":"http:\/\/portillo.ca\/nebula\/wp-json\/wp\/v2\/media?parent=45"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}