class Continuum and Line Analysis Single-dish Software - IRAM
class Continuum and Line Analysis Single-dish Software - IRAM
class Continuum and Line Analysis Single-dish Software - IRAM
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2.5. FITTING AN HYPERFINE STRUCTURE 7<br />
From these equations, we deduce:<br />
τ(vi + p2) = p4 · ri (assumption A4)<br />
�<br />
τi = p4 · S<br />
i<br />
Tant(vi + p2) = p1 � −p4ri 1 − e �<br />
p4<br />
Tant(vi + p2) ≈ p1 · ri<br />
where the last equality holds in the optically thin regime.<br />
This implies that the physical meaning of p4 depends on the value of S: If the relative<br />
intensities are normalized to unity, then S = 1 <strong>and</strong> p4 equals the sum of all centerline opacities.<br />
The results of the HFS fitting procedure are:<br />
<strong>Line</strong> T ant * Tau V lsr Delta V Tau main<br />
1 1.313 ( 0.018) 3.783 ( 0.001) 0.589 ( 0.003) 0.230 ( 0.006)<br />
where T ant * Tau=p1, V lsr=p2, Delta V=p3 <strong>and</strong> Tau main=p4.<br />
According to assumptions A1 <strong>and</strong> A5, the hyperfine structure fitting procedure allows you to<br />
deduce the excitation temperature, since (assuming the Rayleigh-Jeans regime is valid, which is<br />
not true at λ < 3mm...):<br />
Tant(v) = T ∗ A(v) = Beff<br />
Finally, the excitation temperature is given by:<br />
[Tex − Tbg](1 − e<br />
Feff<br />
−τ(v) ) (2.4)<br />
Tex = Tbg + Feff p1<br />
Beff p4<br />
Note: the main group opacity is limited to the range 0.1 − 30 since outside these limits, the<br />
problem becomes degenerate because the opacity no longer appears in the equations (in the<br />
optically thin limit, line ratio no longer depend on the opacity, <strong>and</strong> in the thick limit, exp(−τ) ≪<br />
1).<br />
2.5.4 Typical <strong>Analysis</strong> Sequence<br />
This routine produces the figure 2.1.<br />
1 set plot histo<br />
2 set format brief<br />
3 clear<br />
4 clear alpha<br />
5 file in prov.30m<br />
6 find<br />
7 get f<br />
8 modify source TOTO<br />
9 set unit v c<br />
10 method hfs hfs-n2hp.lin<br />
11 minimize<br />
(2.5)