J/A+A/356/1010 Ortho:para-H2 ratio in C and J shocks (Wilgenbus+, 2000)
The ortho:para-H2 ratio in C- and J- type shock
Wilgenbus D., Cabrit S., Pineau des Forets G., Flower D.
<Astron. Astrophys. 356, 1010 (2000)>
=2000A&A...356.1010W 2000A&A...356.1010W
ADC_Keywords: Atomic physics ; Interstellar medium
Keywords: Magnetohydrodynamics (MHD) - molecular processes - shock waves -
ISM: individual objects: HH 54 - ISM: molecules -
infrared: ISM: lines and bands
Description:
Tables 2a-f and 4a-f contain an extensive set of parameters derived
from our grid of C-type (tables 2a-f) and J-type (tables 4a-f) planar
stationary shock models.
Models of C-type shocks were calculated for
* initial density, nH = 103, 104, 105, 106 cm-3
* shock speed, Vs = 10, 20, 30, 40 km/s (except for nH=106 cm-3,
where Vs = 10, 15, 20, 25, 30 km/s)
* initial ortho:para-H2 ratio, o/p=0.01, 1, 2, 3
Models of J-type shocks were calculated for
* initial density, nH = 103, 104, 105, 106 cm-3
* shock speed, Vs = 5, 10, 15, 20, 25 km/s
* initial ortho:para-H2 ratio, o/p=0.01, 1, 2, 3
Tables 2a and 4a list for each model some selected physical parameters
(width, duration, temperature, ortho:para ratios ...), tables 2b and
4b contain H2-line fluxes in 0-0 S (0)...S(11), 1-0 S(0)...S(11), and
2-0 S(1) lines, and total H2 luminosity. Tables 2c-f and 4c-f list
column densities in each energy level of H2 (each table corresponds to
one initial density nH=103, 104, 105 or 106 cm-3).
File Summary:
--------------------------------------------------------------------------------
FileName Lrecl Records Explanations
--------------------------------------------------------------------------------
ReadMe 80 . This file
table2a.dat 121 68 C-type shock model selected physical parameters
table2b.dat 243 68 C-type shock model H2-line fluxes and luminosity
table2c.dat 32 784 C-type shock models column densities (nH=103cm-3)
table2d.dat 32 784 C-type shock models column densities (nH=104cm-3)
table2e.dat 32 784 C-type shock models column densities (nH=105cm-3)
table2f.dat 32 980 C-type shock models column densities (nH=106cm-3)
table4a.dat 121 80 J-type shock model selected physical parameters
table4b.dat 243 80 J-type shock model H2-line fluxes and luminosity
table4c.dat 32 980 J-type shock models column densities (nH=103cm-3)
table4d.dat 32 980 J-type shock models column densities (nH=104cm-3)
table4e.dat 32 980 J-type shock models column densities (nH=105cm-3)
table4f.dat 32 980 J-type shock models column densities (nH=106cm-3)
--------------------------------------------------------------------------------
Byte-by-byte Description of file: table2a.dat table4a.dat
--------------------------------------------------------------------------------
Bytes Format Units Label Explanations
--------------------------------------------------------------------------------
1- 4 F4.2 --- o/p Initial ortho-to-para H2 ratio
6- 11 E6.0 cm-3 nH Initial density, nH=n(H)+2n(H2)
13- 14 I2 km/s Vs Shock velocity
16- 21 F6.1 10-10T B0 Initial magnetic induction (microGauss)
23- 29 E7.1 cm L Shock width (region where Tn>50K)
31- 37 E7.1 yr t Flow time (for the neutrals to cross the
width of the shock)
39- 45 E7.1 cm-3 nH(50K) Post-shock density at Tn=50K
50- 53 F4.1 km/s |Vi-Vn|max Maximum ion-neutral drift velocity
59- 65 E7.1 --- (H/nH)max Maximum value of n(H)/nH where nH=n(H)+2n(H2)
68- 72 I5 K Tnmax Maximum temperature of neutrals
75- 79 I5 K TjCAM Mean rotational excitation temperature for
levels v=0,J=5...8 (1)
82- 86 I5 K Tjv:1 Mean rotational excitation temperature for
levels v=1,J=1...11 (1)
90- 94 F5.2 --- (o/p)max Local ortho:para-H2 ratio at the rear of the
shock (where Tn=50K)
98-102 F5.2 --- N(o)/N(p) Ratio of column densities of ortho and para
H2, at the rear of the shock
108-112 F5.2 --- (o/p)CAM Mean of empirical o/p ratio for levels
v=0,J=5...8 (2)
117-121 F5.2 --- (o/p)v:1 Mean of empirical o/p ratio for levels
v=1,J=1...11 (2)
--------------------------------------------------------------------------------
Note (1): Excitation temperature is defined in section 4.1, equ. (6).
Note (2): Empirical determination of the observed ortho:para ratio is made using
column densities of H2, see section 4.1, equ. (5).
--------------------------------------------------------------------------------
Byte-by-byte Description of file: table2b.dat table4b.dat
--------------------------------------------------------------------------------
Bytes Format Units Label Explanations
--------------------------------------------------------------------------------
1- 4 F4.2 --- o/p Initial ortho-to-para H2 ratio
6- 11 E6.0 cm-3 nH Initial density, nH=n(H)+2n(H2)
13- 14 I2 km/s Vs Shock velocity
16- 21 F6.1 10-10T B0 Initial magnetic induction
23- 29 E7.1 mW/m2 Lmech Initial mechanical energy flux,
Lmech=1/2.ρ.Vs3, where ρ=1.4nH.m(H)
31- 34 F4.1 % pH2 Percentage of Lmech radiated in all H2 lines
included in the model
36- 39 F4.1 % pCAM Percentage of Lmech radiated in H2
transitions 0-0 S(2)...S(7).
42- 48 E7.1 mW/m2/sr 0-0(0) H2 line flux (1)
50- 56 E7.1 mW/m2/sr 0-0(1) H2 line flux (1)
58- 64 E7.1 mW/m2/sr 0-0(2) H2 line flux (1)
66- 72 E7.1 mW/m2/sr 0-0(3) H2 line flux (1)
74- 80 E7.1 mW/m2/sr 0-0(4) H2 line flux (1)
82- 88 E7.1 mW/m2/sr 0-0(5) H2 line flux (1)
90- 96 E7.1 mW/m2/sr 0-0(6) H2 line flux (1)
98-104 E7.1 mW/m2/sr 0-0(7) H2 line flux (1)
106-112 E7.1 mW/m2/sr 0-0(8) H2 line flux (1)
114-120 E7.1 mW/m2/sr 0-0(9) H2 line flux (1)
122-128 E7.1 mW/m2/sr 0-0(10) H2 line flux (1)
131-137 E7.1 mW/m2/sr 0-0(11) H2 line flux (1)
140-146 E7.1 mW/m2/sr 1-0(0) H2 line flux (1)
148-154 E7.1 mW/m2/sr 1-0(1) H2 line flux (1)
156-162 E7.1 mW/m2/sr 1-0(2) H2 line flux (1)
164-170 E7.1 mW/m2/sr 1-0(3) H2 line flux (1)
172-178 E7.1 mW/m2/sr 1-0(4) H2 line flux (1)
180-186 E7.1 mW/m2/sr 1-0(5) H2 line flux (1)
188-194 E7.1 mW/m2/sr 1-0(6) H2 line flux (1)
196-202 E7.1 mW/m2/sr 1-0(7) H2 line flux (1)
204-210 E7.1 mW/m2/sr 1-0(8) H2 line flux (1)
212-218 E7.1 mW/m2/sr 1-0(9) H2 line flux (1)
220-226 E7.1 mW/m2/sr 1-0(10) H2 line flux (1)
229-235 E7.1 mW/m2/sr 1-0(11) H2 line flux (1)
237-243 E7.1 mW/m2/sr 2-1(1) H2 line flux (1)
--------------------------------------------------------------------------------
Note (1): values < 1.0E-10 are set to 1.0E-99
--------------------------------------------------------------------------------
Byte-by-byte Description of file: table2[cdef].dat table4[cdef].dat
--------------------------------------------------------------------------------
Bytes Format Units Label Explanations
--------------------------------------------------------------------------------
1 I1 --- V Vibrational state of the upper level
3- 4 I2 --- J Rotational state of the upper level
6- 10 I5 K Evj/k Energy of the upper level
12- 15 F4.2 --- o/p Ortho-to-para H2 ratio
17- 22 E6.0 cm-3 nH Initial density, nH=n(H)+2n(H2)
24- 25 I2 km/s Vs Shock velocity
27- 32 F6.1 0.43[cm-2] ln(N(V,J)/g) ln (N(V,J)/g) value (1)
--------------------------------------------------------------------------------
Note (1): ln (N(V,J)/g) = logarithm of the column density in the upper level,
normalized by its statistical weight g=(2J+1)(2I+1), with I=0 for J
even, and I=1 for J odd.
--------------------------------------------------------------------------------
Acknowledgements: David Wilgenbus
(End) David Wilgenbus [DEMIRM, Obs. Paris], Patricia Bauer [CDS] 24-Feb-2000