59 double rec7323 , rec7332, rec3730 , rec3726 , rec2471,
60 reco23tot , reco22tot;
66 fprintf(
ioQQQ,
" lines_lv1_li_ne called\n" );
72 " start level 1 ines" );
75 " cooling due to collisional ionization of heavy elements" );
92 " C 1 1656, collision strength from van Regemoter");
94 linadd(rec,1656,
"REC ",
'i',
95 " C 1 1656 recomb; n.b. coll deexcitation not in" );
99 " C 1 9850, coll excit" );
108 lindst(rec,9850,
"C Ir",ipnt,
'r',
true,
109 " C I 9850 recombination contribution" );
114 " total intensity, all processes, C I 9850");
119 "C 1 8727; equivalent to 4363" );
131 " C 2 158 micron, both e- and H0, H2 in excitation ");
141 " total intensity of C II] 2326, all lines in the multiplet " );
154 " photoproduction, Helfand and Trefftz");
164 " total intensity of C 2 1335");
166 linadd(rec,1335,
"REC ",
'i',
167 " C 2 1335 recombination," );
185 lindst(pump*0.387 * 5.08e-12/(1.+
dense.eden/1e12) ,3920,
"pump",ipnt,
't',
true ,
186 " CII 3918.98/3920.68 is only pumped, no recombination part");
191 pump *= 0.305 * 0.387 * 3.02e-12;
193 " recombination part of C II 6580 line " );
195 " pumped part of line C II 6580" );
199 lindst((rec+pump)/(1.+
dense.eden/1e12),6580,
"TOTL",ipnt,
't',
true ,
200 " total intensity, all processes, C II 6580");
227 " total C 3] 977, recombination + collisional + pumped excitation ");
229 linadd(rec,977,
"C3 R",
'i',
230 " dielectronic recombination contribution to C 3 977 " );
232 linadd(p386,977,
"P386",
'i',
233 " C 3 977 pumped by continuum near 386A" );
238 "C 3 1909 collision, both lines together");
247 " C 3 1908 j-2 to ground" );
250 " the 13C forbidden line of C III " );
257 linadd(3.1e-19*fac*corr,1909,
"C3 R",
'i',
258 " C 3 1909 recombination from Storey" );
261 " C 3 1909 following relax following inner shell photoionization" );
270 " C 3* 1175, excited state line, above 1909 ");
273 " dielectronic recombination contribution to C 3 1175 " );
284 rec,1549,
"TOTL",
'i',
"total intensity of C 4 1549, all processes " );
291 "inward part of C 4 " );
304 " part of line destroyed by photoionization of Balmer continuum " );
306 linadd(rec,1549,
"C4 r",
'i',
307 " recombination C 4 1549 from CV" );
310 " Li seq 2s 3p Li seq transition");
318 for(
int i=0; i <
NI_NDP; ++i )
324 lindst(
nitro.xN3466, 3466,
"N 1", ipnt,
't',
true,
" [N I] 3466.50 only" );
326 lindst(
nitro.xN3467, 3468,
"N 1", ipnt,
't',
true,
" [N I] 3466.54 only" );
330 " [N I] 3466.54, 3466.50 together" );
334 lindst(
nitro.rec5199, 5199,
"TOTr", ipnt,
'r',
true,
335 " estimate of contribution to [N I] 5199 by recombination" );
341 chem =
mole.dissoc_rate(
"N") * 3.82e-12 *
nitro.quench_5200;
342 lindst( chem, 5199,
"chem", ipnt,
'r',
true,
343 " upper limit to [N I] 5199 produced by chemistry" );
348 3.82e-12 *
nitro.quench_5200;
349 lindst( ctRate, 5199,
"H CT", ipnt,
'r',
true,
350 " upper limit to [N I] 5199 produced by charge transfer" );
353 lindst(
nitro.pump5199, 5199,
"pump", ipnt,
'r',
true,
354 " estimate of contribution to [N I] 5199 by FUV pumping" );
357 lindst(
nitro.xN5200 + (ctRate+chem)*0.6, 5200,
"N 1", ipnt,
't',
true,
358 " [N I] 5200 - all processes - stat weight is 6 - total in term is 10" );
360 lindst(
nitro.xN5198 + (ctRate+chem)*0.4, 5198,
"N 1", ipnt,
't',
true,
361 " [N I] 5198 - all processes - stat weight is 4 - total in term is 10" );
364 lindst(
nitro.xN5200+
nitro.xN5198 + (ctRate+chem), 5199,
"TOTL", ipnt,
't',
false,
365 " [N I] 5200 + 5198 together" );
368 lindst(
nitro.xN10397, 10397,
"N 1", ipnt,
't',
true,
" [N I] 10397.7 only");
369 lindst(
nitro.xN10398, 10398,
"N 1", ipnt,
't',
true,
" [N I] 10398.2 only");
370 lindst(
nitro.xN10407, 10407,
"N 1", ipnt,
't',
true,
" [N I] 10407.2 only");
371 lindst(
nitro.xN10408, 10408,
"N 1", ipnt,
't',
true,
" [N I] 10407.6 only");
375 " [N I] 10398.2, 10397.7, 10407.6, 10407.2 together");
387 lindst(
nitro.c6584/(1.+1./2.951),6584,
"N 2",ipnt,
't',
true ,
388 " N 2 6584 alone " );
391 lindst(
nitro.c6584/(1.+2.951),6548,
"N 2",ipnt,
't',
true,
392 " N 2 6548 alone " );
394 efficn2 = 4e-3/(4e-3 + 5.18e-6*
dense.eden/
phycon.sqrte);
397 " N 2 6584 alone, recombination contribution" );
402 0.87*
dense.cdsqte)*3.46e-12;
425 rec =
nitro.xN2_A3_tot * HBeta *
435 " N 2 5755 total, collisions plus charge transfer plus recombination" );
439 " N 2 5755 collisional contribution" );
441 lindst(ctRate,5755,
"C T ",ipnt,
'r',
true,
442 " N 2 5755 charge transfer contribution " );
444 lindst( rec ,5755,
"N 2r",ipnt,
'r',
true,
445 " N 2 5755 recombination contribution" );
448 " N 2 fine structure line ");
451 " N 2 fine structure line " );
454 " N 2 2140 intercombination line " );
463 " N 2 1084, CS guess from g-bar " );
466 " dielectronic recombination contribution to N 2 1085" );
471 linadd(rnii*0.377*0.75*3.02e-12*efficn2,6584,
"N2cn",
'i',
472 " continuum pumped N 2 6584 " );
475 linadd(rnii*0.0117*3.46e-12*efficn2,5755,
"N2cn",
'i',
476 " continuum pumped N 2 5755" );
494 lindst(pump*0.236 * 6.01e-12/(1.+
dense.eden/1e12) ,3311,
"pump",ipnt,
'r',
true,
495 " NII 3311.42 - 3331.31 (6 lines) are only pumped, no recombination part" );
498 lindst(pump*0.186 * 5.18e-12/(1.+
dense.eden/1e12) ,3840,
"pump",ipnt,
'r',
true,
499 " NII 3829.8-3856.06 (6 lines) are only pumped, no recombination part" );
502 lindst(pump*0.025 * 5.52e-12/(1.+
dense.eden/1e12) ,3609,
"pump",ipnt,
'r',
true,
503 " NII 3593.60/3609.1/3615.86 (3 lines) are only pumped, no recombination part" );
507 lindst(pump*0.186*0.595 * 4.31e-12/(1.+
dense.eden/1e12) ,4640,
"pump",ipnt,
'r',
true ,
508 " NII 4601.5-4643.1 (6 lines) are only pumped, no recombination part");
511 lindst(pump*0.025*0.442 * 3.97e-12/(1.+
dense.eden/1e12) ,5010,
"pump",ipnt,
'r',
true,
512 " NII 5002.7/5010.6/5045.1 (3 lines) are only pumped, no recombination part " );
517 pump *= 0.236 * 0.626 * 3.50e-12;
520 " recombination part of N II 5679 line" );
523 " pumped part of line N II 5679 " );
526 lindst((rec+pump)/(1.+
dense.eden/1e12),5679,
"TOTL",ipnt,
't',
true,
527 " total intensity, all processes, N II 5679 " );
530 " [N 3] 57 micron fine structure line");
539 " total intensity of N III] 1750, all lines in the multiplet " );
573 rn3tot = (rb + raten3)*0.439 + rn3mor;
575 sum = raten3*4.29e-12;
577 lindst(sum,4640,
"N3cn",ipnt,
'r',
true ,
578 " continuum pumped \"Bowen\" N 3, optically thin excited line ");
581 sum = rb*4.29e-12*0.834;
583 lindst(sum,4634,
"N3cn",ipnt,
'r',
true ,
584 " continuum pumped \"Bowen\" N 3, optically thin excited line");
586 sum = rb*4.29e-12*(1. - 0.834);
588 lindst(sum,4642,
"N3cn",ipnt,
'r',
true,
589 " continuum pumped \"Bowen\" N 3, optically thin excited line" );
604 " N 3 990, all processes ");
606 linadd(rec+rn3tot*2.01e-11*fac,990,
"extr",
'i',
607 " total N 3 990, both electron excitation and continuum pumping" );
609 linadd(rec,990,
"rec ",
'i',
610 " part of N 3 990 due to recombination " );
612 linadd(rn3tot*2.01e-11,990,
"N 3p",
'r',
613 " N 3 989.8, continuum pumped" );
616 " N 4] 1486, total intensity of both lines" );
623 " the N IV] slow transition by itself " );
634 " N 4 765, collisionally excited");
637 " N 4 765 recombination," );
640 " continuum pumping of NV 1240, N 5 1240, total emission, collisions plus pumping " );
644 linadd(sum,1240,
"Inwd",
'i',
645 " inward part of N 5 " );
652 " N 5 209, 2s-3p Li seq ");
656 " oxygen total Oxygen I 6300, including continuum optical depth ");
668 rate_OH_dissoc =
mole.findrate(
"PHOTON,OH=>O,H");
669 r12 = rate_OH_dissoc * 0.55 * 3.16e-12 *
CoolHeavy.c6300_frac_emit;
672 6300.,
"OH p",ipnt ,
'i' ,
false,
673 " the intensity of [OI] 6300 line due to OH photodestruction");
676 6363.,
"OH p",ipnt ,
'i' ,
false,
677 " the intensity of [OI] 6363 line due to OH photodestruction ");
681 " total Oxygen I 6363, including continuum optical depth " );
687 r13 = rate_OH_dissoc * 0.05 * 3.57e-12 * 0.94*
CoolHeavy.c5577_frac_emit;
688 lindst( r13 , 5577.,
"OH p",ipnt ,
'i',
false,
689 " 94% of excitations to highest level decay via 5577" );
692 " O I fine structure line ");
695 " O I fine structure line ");
698 " total collisional cooling due to 6-level OI atom" );
701 " total collisional heating due to 6-level OI atom " );
709 lindst(sum,8446,
"6lev",ipnt,
'i',
false,
710 " OI 8446 from six level atom" );
714 lindst(sum,1304,
"6lev",ipnt,
'i',
false,
715 " OI 1304 from six level atom " );
719 lindst(sum,1039,
"6lev",ipnt,
'i',
false ,
720 " OI 1039 from six level atom");
724 lindst(sum,4368,
"6lev",ipnt,
'i',
false,
725 " OI 4368 from six level atom" );
729 lindst(sum,13100,
"6lev",ipnt,
'i',
false ,
730 " OI 1.3 micron from six level atom");
734 lindst(sum,11300,
"6lev",ipnt,
'i',
false ,
735 " OI 1.1 micron from six level atom");
739 lindst(sum,29000,
"6lev",ipnt,
'i',
false ,
740 " OI 2.9 micron from six level atom");
744 lindst(sum,46000,
"6lev",ipnt,
'i',
false ,
745 " OI 4.6 micron from six level atom");
763 reco23tot =
CoolHeavy.O2_A3_tot * HBeta *
784 rec2471 = reco23tot *
CoolHeavy.O2471*2471./7325. * 8.05e-12/2.72e-12;
801 rec3726 = reco22tot *
CoolHeavy.O3726 * 5.34e-12;
802 rec3730 = reco22tot *
CoolHeavy.O3730 * 5.34e-12;
806 9.7e-5/(9.7e-5 +
dense.eden*1.15e-6/
phycon.sqrte));
810 lindst(fac ,3727,
"TOTL",ipnt,
't',
true,
811 " O II 3727, all lines of multiplet together " );
815 lindst( fac ,7325,
"TOTL",ipnt,
't',
true ,
816 " O II 7325, all lines of multiplet together");
819 " line produced by photoionization of Oo; already in TOTL" );
824 " line produced by photoionization of Oo; already in TOTL" );
827 " collisional contribution to line " );
830 " collisional contribution to line " );
834 " five level atom calculations; D5/2 - S3/2" );
837 " D3/2 - S3/2 transition" );
841 " both 2P 1/2 and 3/2 to ground " );
845 " P1/2-D5/2 and P3/2-D5/2 together" );
848 " P1/2-D3/2 and P3/2-D3/2 together " );
850 linadd( rec3730 ,3729,
"O 2r",
'i',
851 " recombination contribution refer o2 rec Liu, X-W., Storey, P.J., Barlow, M.J., Danziger, I.J.,refercon Cohen, M., & Bryce, M., 2000, MNRAS, 312, 585 recombination contributions five level atom calculations; D5/2 - S3/2 " );
853 linadd( rec3726 ,3726,
"O 2r",
'i',
854 " D3/2 - S3/2 transition" );
856 linadd(rec2471,2471,
"O 2r",
'i',
857 " both 2P 1/2 and 3/2 to ground " );
858 linadd(rec7323,7323,
"O 2r",
'i',
859 " P1/2-D5/2 and P3/2-D5/2 together " );
861 linadd(rec7332,7332,
"O 2r",
'i',
862 " P1/2-D3/2 and P3/2-D3/2 together " );
865 " O II 833.8 coll excit ");
881 lindst(pump*0.336 * 6.37e-12/(1.+
dense.eden/1e12) ,3120,
"pump",ipnt,
'r',
true,
882 " OII 3113.62 - 3139.68 (8 lines) are only pumped, no recombination part" );
885 lindst(pump*0.147 * 6.03e-12/(1.+
dense.eden/1e12) ,3300,
"pump",ipnt,
'r',
true,
886 " OII 3277.56 - 3306.45 (6 lines) are only pumped, no recombination part" );
889 lindst(pump*0.087 * 5.29e-12/(1.+
dense.eden/1e12) ,3762,
"pump",ipnt,
'r',
true,
890 " OII 3739.76/3762.47/3777.42 (3 lines) are only pumped, no recombination part" );
895 lindst(rec,4651,
"O 2r",ipnt,
'r',
true,
896 " O II 4651 total recombination, 4638.86-4696.35 (8 lines) " );
900 linadd(pump* 0.336 * 0.933 * 4.27e-12/(1.+
dense.eden/1e12),4651,
"O 2p",
'i',
901 " pumped part of line O II 4651 " );
907 " recombination contribution to O II 4341 line " );
909 linadd(pump* 0.147 * 0.661 * 4.58e-12/(1.+
dense.eden/1e12),4341,
"O 2p",
'i',
910 " pumped part of line O II 4341 " );
913 lindst(rec+pump* 0.147 * 0.661 * 4.58e-12/(1.+
dense.eden/1e12),4341,
"TOTL",ipnt,
't',
true,
914 " total intensity, all processes, O II 4341" );
920 linadd(rec/(1.+
dense.eden/1e12),3736,
"O 2r",
'i',
"\n recombination part of O II 3736 line " );
921 linadd(pump* 0.087 * 0.763 * 5.33e-12/(1.+
dense.eden/1e12),3736,
"O 2p",
'i',
922 " pumped part of line O II 3736" );
925 lindst((rec+pump* 0.087 * 0.763 * 5.33e-12)/(1.+
dense.eden/1e12),3736,
"TOTL",ipnt,
't',
true,
926 " total intensity, all processes, O II 3736" );
933 "total intensity of OIII] 1665, all processes " );
939 " contribution to OIII 1665 due to inner shell (2s^2) ionization " );
941 linadd(
oxy.AugerO3*1.20e-11*efac*0.27,1665,
"Augr",
'i',
942 " contribution to OIII 1665 due to K-shell ionization " );
946 " O III 5007 alone, collisions, tot OIII is this times 1.333 fac = c5007/(1.+1./2.887) >>chng 01 may 04, branching ratio had been 2.887, revised to 3 as per refer o3 as Storey, P.J., & Zeippen, C.J., 2000, 312, 813-816 ");
950 " O III 4959 alone, collisions, tot OIII is this times 4" );
954 " O III 4931 alone, collisions >>chng 01 jul 11, added this line >>refer o3 as Nussbaumer, H., & Storey, P., 1981, A&A, 99, 177 >>refer o3 as Mathis, J.S., & Liu, X.-W., 1999, ApJ, 521, 212-216 ");
957 " O III 5007 lost through excited state photo" );
960 effec = 1.6/(1.6 + 0.9*
dense.cdsqte);
974 linadd(fac,4363,
"TOTL",
'i',
975 " O III 4363, sum of rec, coll, ct excitation" );
979 " O III 4363,collisions from five level atom " );
981 lindst(r4363,4363,
"Rec ",ipnt,
'r',
true,
982 " O III 4363 recombination, coefficient from Burgess and Seaton " );
986 " collisional excitation of 2321, 5-level atom");
987 linadd(ct4363,4363,
"C EX",
'i' ,
988 " call linadd( c4363*0.236 , 2321 , 'O 3','c') charge exchange, Dalgarno+Sternberg ApJ Let 257, L87. ");
991 5592,
"C EX",
'i',
" charge exchange rate, D+S " );
994 " O III 88 micron, collisionally excited");
997 " O III 52 micron, collisionally excited ");
1005 " O III 834A, collisions and dielectronic recombination ");
1008 " O III 834A, dielectronic recombination only" );
1011 " O IV 26 micron ");
1020 " total intensity of O IV] 1402, all lines in the multiplet " );
1034 " inner shell photoionization, relaxation " );
1044 " O IV 789A, dielectronic recombination only" );
1053 " O V 630A, dielectronic recombination only" );
1056 " O V 1218], total intensity of both lines " );
1060 " the slow transition by itself" );
1065 " BS O V 5112, recombination " );
1068 " O VI 1035, total of pumping and collisional excitation " );
1072 linadd(sum,1035,
"Inwd",
'i',
1073 " inward part of OVI line" );
1080 "O VI 150, Li seq 2s 3p ");
1083 " neon Neon II 12.8 micron ");
1086 " Ne III fine structure line ");
1089 " Ne III fine structure line ");
1093 " Ne III 3869, of 3968+3869 doublet" );
1097 " Ne III 3968, of 3968+3869 doublet" );
1101 " NeIII auroral line " );
1105 " NeIII auroral line");
1109 " Ne IV 2424, collisional excitation" );
1113 " Ne IV N=3-2 lines, three level atom approx, this is the sum of the 4714.5, 4724.2, 4725.5 lines" );
1117 " Ne IV N=3 lines, three level atom approx");
1121 " Ne V 3426 of 3426, 3346 doublet" );
1125 " Ne V 3346 of 3426, 3346 doublet " );
1133 " collisionally excited" );
1140 lindst(
CoolHeavy.c1134,1141,
"Ne 5",ipnt,
't',
true,
" both components of 5S-3P 1146.1, 1137.0 doublet " );
1145 " Ne VII 895, collisionally excited, both lines " );
1148 " Ne VII 895, only fast transition ");
1151 " Ne VII 890, single line " );
1154 " Ne VIII 774, collisionally excited " );
1158 linadd(sum,774,
"Inwd",
'i',
1159 " inward part of NeVIII 774 line" );
1162 " the NeVIII 770 780 doublet ");
1167 " Ne VIII 88 2s 3p, collisionally excited ");
1171 fprintf(
ioQQQ,
" lines_lv1_li_ne returns\n" );