This HTML file was created by the program doc_tsuite.pl. doc_tsuite.txt contains a tab delimited list of the files.


feii_blr_n09_p18.in FeII model


title FeII model
table agn
atom feii
print lines column
print lines sort intensity
hden 9
phi(h) 18
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n09_p18.lin" last
punch feii column densities "feii_blr_n09_p18.fe2col" last
punch feii relative populations all "feii_blr_n09_p18.lev" last
punch feii relative populations range 0 200 "feii_blr_n09_p18.lv1" last
punch feii relative populations range 201 371 "feii_blr_n09_p18.lv2" last
punch overview "feii_blr_n09_p18.ovr" last
punch dr "feii_blr_n09_p18.dr" last
punch convergence reason "feii_blr_n09_p18.cnv"

This is one of the 5 models that sample the LOC plane.


feii_blr_n09_p18_Z20.in FeII model


title FeII model
table agn
print lines column
print lines sort intensity
atom feii
abundances starburst 20
hden 9
phi(h) 18
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n09_p18_Z20.lin" last
punch feii colum density "feii_blr_n09_p18_Z20.fe2col" last
punch feii populations all "feii_blr_n09_p18_Z20.lev" last
punch overview "feii_blr_n09_p18_Z20.ovr" last
punch dr "feii_blr_n09_p18_Z20.dr" last
punch convergence reason "feii_blr_n09_p18_Z20.cnv"

This is one of the 5 models that sample the LOC plane.


feii_blr_n11_p20.in FeII model


title FeII model
table agn
print lines column
print lines sort intensity
atom feii
hden 11
phi(h) 20
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n11_p20.lin" last
punch feii column densities "feii_blr_n11_p20.fe2col" last
punch feii populations all "feii_blr_n11_p20.lev" last
punch overview "feii_blr_n11_p20.ovr" last
punch dr "feii_blr_n11_p20.dr" last

This is one of the 5 models that sample the LOC plane.


feii_blr_n11_p20_Z20.in FeII model


title FeII model
table agn
print lines column
print lines sort intensity
atom feii
failures 2
abundances starburst 20
hden 11
phi(h) 20
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch cooling "feii_blr_n11_p20_Z20.col"
punch feii lines "feii_blr_n11_p20_Z20.lin" last
punch feii column densites "feii_blr_n11_p20_Z20.fe2col" last
punch feii populations all "feii_blr_n11_p20_Z20.lev"
punch heating "feii_blr_n11_p20_Z20.het"
punch element calcium "feii_blr_n11_p20_Z20.ca"
punch overview "feii_blr_n11_p20_Z20.ovr" last
punch dr "feii_blr_n11_p20_Z20.dr"

This is one of the 5 models that sample the LOC plane.


feii_blr_n12_p19.in FeII model


title FeII model
table agn
print lines column
print lines sort intensity
atom feii
hden 12
phi(h) 19
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n12_p19.lin" last
punch feii populations all "feii_blr_n12_p19.lev" last
punch feii column densities "feii_blr_n12_p19.fe2col" last
punch overview "feii_blr_n12_p19.ovr" last
punch dr "feii_blr_n12_p19.dr" last

This is one of the 5 models that sample the LOC plane.


feii_blr_n12_p19_Z20.in FeII model


title FeII model
table agn
print lines column
print lines sort intensity
atom feii
abundances starburst 20
hden 12
phi(h) 19
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n12_p19.lin" last
punch feii column densiteis "feii_blr_n12_p19.fe2col" last
punch feii populations all "feii_blr_n12_p19.lev" last
punch convergence reason "feii_blr_n12_p19.cvg"
punch overview "feii_blr_n12_p19.ovr" last
punch dr "feii_blr_n12_p19.dr" last

This is one of the 5 models that sample the LOC plane.


feii_blr_n13_p18.in BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1


title BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1
print lines column
print lines sort intensity
atom feii
c
c this is a very low ionization cloud
c the conditions, and some lines, are surprisingly sensitive
c to the treatment of hydrogen molecules
c
table agn
hden 13
phi(h) 18
stop column density 22
double
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n13_p18.lin" last
punch feii column densites "feii_blr_n13_p18.fe2col" last
punch feii populations all "feii_blr_n13_p18.lev" last
punch overview "feii_blr_n13_p18.ovr" last
punch dr "feii_blr_n13_p18.dr" last
punch convergence reason "feii_blr_n13_p18.cvr"

This is one of the 5 models that sample the LOC plane.


feii_blr_n13_p18_Z20.in BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1


title BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1
print lines column
print lines sort intensity
atom feii
abundances starburst 20
c
c this is a very low ionization cloud
c the conditions, and some lines, are surprisingly sensitive
c to the treatment of hydrogen molecules
c
table agn
hden 13
phi(h) 18
stop column density 19
double
iterate convergence limit 15
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n13_p18_Z20.lin" last
punch feii column densities "feii_blr_n13_p18_Z20.fe2col" last
punch feii populations all "feii_blr_n13_p18_Z20.lev" last
punch overview "feii_blr_n13_p18_Z20.ovr" last
punch dr "feii_blr_n13_p18_Z20.dr" last
punch convergence reason "feii_blr_n13_p18_Z20.cvr"

This is one of the 5 models that sample the LOC plane.


feii_blr_n13_p22.in FeII model


title FeII model
table agn
print lines column
print lines sort intensity
atom feii
hden 13
phi(h) 22
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n13_p22.lin" last
punch feii populations all "feii_blr_n13_p22.lev" last
punch feii column density "feii_blr_n13_p22.col" last
punch overview "feii_blr_n13_p22.ovr" last
punch dr "feii_blr_n13_p22.dr" last

This is one of the 5 models that sample the LOC plane.


feii_blr_n13_p22_Z20.in FeII model


title FeII model
table agn
print lines column
print lines sort intensity
atom feii
abundances starburst 20
hden 13
phi(h) 22
stop column density 23
iterate convergence
normalize to "totl" 1216 100
init "c84.ini"
punch feii lines "feii_blr_n13_p22_Z20.lin" last
punch feii column densities "feii_blr_n13_p22_Z20.fe2col" last
punch feii populations all "feii_blr_n13_p22_Z20.lev" last
punch overview "feii_blr_n13_p22_Z20.ovr" last
punch dr "feii_blr_n13_p22_Z20.dr"
punch heating "feii_blr_n13_p22_Z20.het"
punch cooling "feii_blr_n13_p22_Z20.col"

This is one of the 5 models that sample the LOC plane.


h2_HTT91.in Hollenbach et al. 1991 low-density PDR


title Hollenbach et al. 1991 low-density PDR
iterate
c
c
hden 3
c
c this is hot star continuum
black 30,000
intensity 0.2 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c
turbulence 1.5 km/sec
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
grains orion, abundance log 0.16
grains PAH
atom h2
init file="ism.ini"
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl-7 ar-8 fe -6.60
normalize to 157.6m "C 2"
sphere
case b
cosmic rays, background
c stop when gas is fully neutral
stop efrac -10
c this is to stop at an intended place, since results would be
c very dependent on details if we stop on temperature
stop thickness 18.954
c stop when gas is cold
stop temperature 10 linear
c add this to mimic unmodelled neutral gas
double optical depths
c
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print ages
c
punch overview last "h2_HTT91.ovr"
punch dr "h2_HTT91.dr "
punch molecules last "h2_HTT91.mol"
punch results last "h2_HTT91.rlt"
punch continuum units microns last "h2_HTT91.con"
punch heating last "h2_HTT91.het"
punch cooling last "h2_HTT91.col"
c
c h2_HTT91.in

This is the Hollenbach et al 1991 Low-density PDR The case b command appears because the Lyman lines are vastly optically thick in this environment. If the command is removed the Lyman lines will be optically thin, which will result in fast fluorescent excitation of the (nearly totally neutral) hydrogen atom.

there is very little CO in this model since it is not thick enough for the UV pumping lines to become optically thick


h2_orion_hii_pdr.in constant gas pressure orion into pdr


title constant gas pressure orion into pdr
stop temperature 10 linear
c
c print lots of faint CO lines
print line faint -4
stop AV 1,000 point
constant gas pressure
sphere
c
c the incident continuum is two parts
c star and flux of photons striking it
table star kurucz 39,600K
Q(H) 49
radius 17.4507
c plus hot brems
brems 6
phi(h) 10
c
c add cosmic rays, which are important at depth
cosmic rays, background
c
c we have a spherical geometry but want to simulate observing
c through a spectrometer's slit. use the aperture
c command for this
aperture beam
c
c the observed microturbulence
turbulence 8 km/sec no pressure
c
c density and abundances
hden 4
init file="ism.ini"
abundances hii region no grains
grains orion
grains pah
atom h2
c
set nend 2000
c
punch overview last "h2_orion_hii_pdr.ovr"
punch hydrogen 21 cm last "h2_orion_hii_pdr.21cm"
punch heating "h2_orion_hii_pdr.het"
punch cooling "h2_orion_hii_pdr.col"
punch dr last "h2_orion_hii_pdr.dr"
punch results last "h2_orion_hii_pdr.rlt"
punch continuum last "h2_orion_hii_pdr.con" units microns
punch hydrogen lya last "h2_orion_hii_pdr.lya"
punch grain charge last "h2_orion_hii_pdr.grnchr"
punch grain potential last "h2_orion_hii_pdr.grnpot"
punch H2 lines last "h2_orion_hii_pdr.h2lin"
punch H2 column density last "h2_orion_hii_pdr.h2col"
punch molecules "h2_orion_hii_pdr.mol"
c
c h2_orion_hii_pdr.in

05 dec 19, had stopped at 1 pc, gas/dust got very cold and H2 stopped forming, H2 went to H0 - this below the CO network out of the water. stop at AV of 1000 instead. Is this loss of H2 at low grain temperature physical? Is it ever seen?


h2_pdr_leiden_f2.in low density high flux model 2 as defined in e-mail


title low density high flux model 2 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use
c half the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c hydrogen density
hden 3.
c
c commands controlling geometry =========
c
c other commands for details =========
atom h2
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch overview "h2_pdr_leiden_f2.ovr"
punch leiden lines "h2_pdr_leiden_f2.lin"
punch leiden "h2_pdr_leiden_f2.lei"
punch dr "h2_pdr_leiden_f2.dr"
punch grain physical conditions "h2_pdr_leiden_f2.grn"
punch H2 lines "h2_pdr_leiden_f2.h2lin" all
punch H2 column density "h2_pdr_leiden_f2.h2col"
punch H2 populations matrix zone "h2_pdr_leiden_f2.pop"
c
c
c h2_pdr_leiden_f2.in
c class pdr
c ========================================
c


h2_pdr_leiden_f3.in high density low flux model 3 as defined in e-mail


title high density low flux model 3 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
c ice formation is exteme and will establish the thickness, which would
c make this very sensitive to details - turn off ices
no grain molecules
atom h2
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch overview "h2_pdr_leiden_f3.ovr"
punch leiden lines "h2_pdr_leiden_f3.lin"
punch leiden "h2_pdr_leiden_f3.lei"
punch dr "h2_pdr_leiden_f3.dr"
punch grain physical conditions "h2_pdr_leiden_f3.grn"
punch H2 lines "h2_pdr_leiden_f3.h2lin" all
punch H2 column density "h2_pdr_leiden_f3.h2col"
punch H2 populations matrix zone "h2_pdr_leiden_f3.pop"
punch H2 rates "h2_pdr_leiden_f3.rat"
c
c
c h2_pdr_leiden_f3.in
c class pdr
c ========================================
c


h2_pdr_leiden_f4.in high density high flux model 4 as defined in e-mail


title high density high flux model 4 as defined in e-mail
c
c commands controlling continuum =========
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c
c commands controlling geometry =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
c
c hydrogen density
hden 5.5
c
c other commands for details =========
atom h2
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch overview "h2_pdr_leiden_f4.ovr"
punch leiden lines "h2_pdr_leiden_f4.lin"
punch leiden "h2_pdr_leiden_f4.lei"
punch dr "h2_pdr_leiden_f4.dr"
punch grain physical conditions "h2_pdr_leiden_f4.grn"
punch H2 lines "h2_pdr_leiden_f4.h2lin" all
punch H2 column density "h2_pdr_leiden_f4.h2col"
punch H2 populations matrix zone "h2_pdr_leiden_f4.pop"
punch H2 solomon "h2_pdr_leiden_f4.sol"
c
c
c h2_pdr_leiden_f4.in
c class pdr
c ========================================
c


h2_pdr_leiden_hack_f1.in low density low flux model 1 as defined in e-mail


title low density low flux model 1 as defined in e-mail
atom h2
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
c
c commands for density & abundances =========
c hydrogen density
hden 3.
grains ism 1.16
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch leiden lines "h2_pdr_leiden_hack_f1.lin"
punch leiden "h2_pdr_leiden_hack_f1.lei"
punch dr "h2_pdr_leiden_hack_f1.dr"
punch grain physical conditions "h2_pdr_leiden_hack_f1.grn"
punch overview "h2_pdr_leiden_hack_f1.ovr"
c
punch H2 lines "h2_pdr_leiden_hack_f1.h2lin"
punch H2 column density "h2_pdr_leiden_hack_f1.h2col"
punch H2 populations matrix zone "h2_pdr_leiden_hack_f1.pop"
c

c h2_pdr_leiden_hack_f1.in c class pdr c ======================================== c


h2_pdr_leiden_v1.in model 5 as defined in e-mail


title model 5 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c is half the requested value
c they want for the model is actually twice the value below
table draine 5
c
c commands for density & abundances =========
c hydrogen density
hden 3.
grains ism 1.16 no qheat
c add PAHs
grains PAH no qheat 3 function
atom h2
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c
c commands controlling output =========
punch overview "h2_pdr_leiden_v1.ovr"
punch leiden lines "h2_pdr_leiden_v1.lin"
punch leiden "h2_pdr_leiden_v1.lei"
punch dr "h2_pdr_leiden_v1.dr"
punch grain dust temperature "h2_pdr_leiden_v1.grn"
c
c
c
c
c h2_pdr_leiden_v1.in
c class pdr
c ========================================
c


h2_pdr_leiden_v2.in model 6 as defined in e-mail


title model 6 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
c
c commands for density & abundances =========
c add PAHs and grains
grains PAH no qheat 3 function
grains ism 1.16 no qheat
c hydrogen density
hden 3.
c
c commands controlling geometry =========
c
c other commands for details =========
atom h2
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c
c commands controlling output =========
punch overview "h2_pdr_leiden_v2.ovr"
punch leiden lines "h2_pdr_leiden_v2.lin"
punch leiden "h2_pdr_leiden_v2.lei"
punch dr "h2_pdr_leiden_v2.dr"
punch grain dust temperature "h2_pdr_leiden_v2.grn"
c
c
c h2_pdr_leiden_v2.in
c class pdr
c ========================================
c


h2_pdr_leiden_v3.in model 7 as defined in e-mail


title model 7 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
c
c commands for density & abundances =========
c add PAHs and grains
grains PAH no qheat 3 function
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
atom h2
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c
c commands controlling output =========
punch overview "h2_pdr_leiden_v3.ovr"
punch leiden lines "h2_pdr_leiden_v3.lin"
punch leiden "h2_pdr_leiden_v3.lei"
punch dr "h2_pdr_leiden_v3.dr"
punch grain dust temperature "h2_pdr_leiden_v3.grn"
punch heating "h2_pdr_leiden_v3.het"
punch cooling "h2_pdr_leiden_v3.col"
punch molecules "h2_pdr_leiden_v3.mol"
c
c
c h2_pdr_leiden_v3.in
c class pdr
c ========================================
c


h2_pdr_leiden_v4.in model 8 as defined in e-mail


title model 8 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
c
c commands for density & abundances =========
grains PAH no qheat 3 function
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
atom h2
c
c commands controlling output =========
punch overview "h2_pdr_leiden_v4.ovr"
punch leiden lines "h2_pdr_leiden_v4.lin"
punch leiden "h2_pdr_leiden_v4.lei"
punch dr "h2_pdr_leiden_v4.dr"
punch grain dust temperature "h2_pdr_leiden_v4.grn"
c
c
c h2_pdr_leiden_v4.in
c class pdr
c ========================================
c


h2_th85orion.in Tielens and Hollenbach pdr model with orion grains, Table 2, paper b


title Tielens and Hollenbach pdr model with orion grains, Table 2, paper b
failures 3
c
print line sort intensity
print line column
iterate
c
c cosmic background
background
cosmic rays, background
c
set nend 2000
c simulate effects of gas we do not model
double
c
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
blackbody, t=75K
intensity 2.7 (total)
c
c this is the attenuated hot stellar continuum
black 30,000
intensity 2.204 range 0.4412 to 1 Ryd
c
c this will remove all ionizing radiation
extinguish 24 0
c
hden 5.362
grains orion
grains PAH
atom h2
init file="ism.ini"
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl-7 ar-8 fe -6.60
turbulence 2.7 km/sec
normalize to "C 2" 157.6m
sphere
case b
cosmic rays, background
c this should be the one actually used
stop AV 33.2
c stop when gas is fully neutral
stop efrac -10
c stop when gas is cold
stop temperature 10 linear
c stop at thickness so that would make results very
c detail dependent
stop thickness 19
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print ages
punch overview last "h2_th85orion.ovr" no hash
punch pdr last "h2_th85orion.pdr"
punch molecules last "h2_th85orion.mol"
punch dr "h2_th85orion.dr"
punch hydrogen physical conditions last "h2_th85orion.hyd"
punch results last "h2_th85orion.rlt"
punch continuum units microns last "h2_th85orion.con"
punch fine conti unit last micron "h2_th85orion_con.fin" range 0.7 1. skip 2
* punch fine opaci unit micron "h2_th85orion_opc.fin" range 0.7 1. skip 2
punch heating last "h2_th85orion.het"
punch cooling last "h2_th85orion.col"
punch hydrogen 21 cm last "h2_th85orion.21cm"
c
c
c h2_th85orion.in

This is the Tielens and Hollenbach (1985a, b) standard model of the Orion photodissociation region (PDR). Specifically, this is my attempt at their Table 2 of paper 2, ApJ 291, p749. The case b command appears because the Lyman lines are vastly optically thick in this environment. If the command is removed the Lyman lines will be optically thin, which will result in fast fluorescent excitation of the (nearly totally neutral) hydrogen atom.

// >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985a, ApJ, 291, 722 // >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985b, ApJ, 291, 746