An atomic and molecular database for analysis of submillimetre line observations

1. An atomic and molecular database for analysis of submillimetre line observations Schöier, F. L.; van der Tak, F. F. S.; van Dishoeck, E. F.; Black, J. H. reporter: 李英杰

2. 作者简介 F.L.Schöier: Stockholm Observatory, Alba nova University center, Stockholm univerci ty, Leiden Observatory; Sweden. F.F.S. Van der Tak: Max-Planck-Institut für Radioastronomie. J.H.Black: Onsala Space Observatory; Sweden. E. F. van Dishoeck: Leiden Observatory; The Netherlands.

3. E. F. van Dishoeck • 莱顿天文台的分子天体物理教授，是荷兰皇家科学院和美国国家科学院成员。 • 2000获荷兰最高科学奖(Spinozaprem ie), 2001年获英国the Bourke Award of the Royal Society of Chemistry，and the Gold Medal of the Royal Dutch Chemical Society in 1994。 • She was a Junior Fellow of the Harvard Society of Fellows。

4. Abstract • 分子能级、统计权重、爱因斯坦-A系数； • 各种分子的碰撞速率系数； • 利用量子化学理论计算和实验将碰撞速率系数推至高能； • 数据库网址： http://www.strw.leidenuniv.nl/∼moldata 这一数据库可为non-LTE辐射转移过程提供数据，是当前和未来(亚)毫米和红外的重要分析工具。

5. 1.Introduction • 现已发现各种各样的分子，这些分子是探测星际和环恒星包层的很好探针，如： • CO是探测总气体质量的探针； • CS、HCO+和HCN常用来探测密度结构； • H2CO和CH3OH被用来探知温度与密度结构(mm与亚mm波) • CH3OCH3和CH3CN探测化学构造； • SiO和SO2探测激波； • 为提取天体物理参数，需考虑谱线辐射转移，常用的是局部近似，如LVG().当然更复杂的非局部近似，近来也有发展。

6. LVG模型 LVG即大速度梯度模型是用来近似计算逃逸概率的。其基本思想是：一个局部区域中产生的辐射只与局部的分子、原子发生相互作用，此相互作用的尺度为∆𝑙=∆𝑣/(𝑑𝑣/𝑑𝑙)，其中，∆𝑣是运动介质气体的局部热运动或湍动运动的速度宽度，dv/dt是系统运动介质的速度梯度。 求出能级布局数，从而可求分子谱线的激发温度、发射强度以及谱线亮温等，并且调整模型参数使其与观测相符可得感兴趣源的物理特征和运动学特征。 但如若有脉泽出现，运用LVG就会有问题，因为脉泽辐射强度与光深的绝对值成正相关，而光深则与速度场梯度呈负相关，因而若出现脉泽则表明速度场梯度不能很大，因而用非局域近似可能会更好。

7. 数据库简介 • 网站地址：http://www.strw.leidenuniv.nl/∼moldata. 该网站数据大部分都是基于振动基态的，只有少部分分子考虑了低的振动激发态，如CO。同时该数据还包含了少数多原子分子数据，如HCN， H2CO 等。点击选择后会得到一个.dat的文件，包含能级、统计权重、跃迁频率、自发发射概率、碰撞速率系数、跃迁强度以及计算的来源。 • RADEX：包含在上面的网站中，用逃逸概率的方法来编码统计平衡下的辐射转移方程的。RADEX与LVG模型相互比较，可用来快速分析大量的观测数据，从而获取感兴趣区的物理条件，如：密度、动能温度、旋转温度，以及粗略的估计柱密度等。

8. 2. Energy levels • 1.分子态由电子能级表出：2S+1ΛΛ+Σ.习惯上常用Λ=Σ，Π，Δ，Φ分别指代Λ=0，1，2，3，其中|Λ|=|ML|为沿轴向分量的轨道角动量。 • 2.分子的能级由电子能级、振动能级和转动能级组成。 • 求解电子能级，首先得求出不同电子态的势能曲线，用的是量子化学中的从头计算法。 • 若分子势能曲线有极小值则说明相应的电子态是稳定的。

9. 1. Energy levels 本文关注的是mm与亚mm波段，其中落在mm与亚mm波段的有： 附注：线性分子是这一陀螺模型中竖直方向长度为零的情况。

10. 2.1 General consideration 本文中能级一般仅包含振动基态，能量可高至~1000K。 只有CO、CS分子包含了振动激发态。 正(ortho)或仲(para)分子或是甲醇的A-和E-型将被作为不同分子分别对待。

11. 正、仲氢与A-、E-甲醇 对于正仲之分可以以氢分子为例。对于正态氢，两个氢核的自旋平行，总核自旋为1，只有奇数的转动量子数J存在；对于仲态氢则两个氢核自旋反平行，只有偶数的J存在。因正仲间的辐射跃迁是不可能的，且观测到的偶J和奇J的能级布局在低J态时明显的处于热平衡，因而可以用正仲态的布局数比率可靠地估计分子云的温度。 对于A-和E-，它们是由于甲醇的三重对称造成的。在E型中，-J<k<J，其中k为转动角动量在甲基三个氢原子所在平面上的投影量子数，甲醇是振动双重简并的；在A型中，对于k≥0，0≤k≤J，由于不对称造成doublets split，分别记为A+和A-。

12. 2.2 Linear molecules • 对于电子态为1Σ的线性双原子、多原子分子：E=BJ(J+1)，B=(2I)-1 • 重线性分子统计权重：g=(2J+1).如HC3N • 若考虑离心变形则需要高次项：∝J2(J+1)2. • 对于nΣ，n≠1的则需要考虑总角动量与总电子自旋的耦合效应。如SO、CN。 • 对于基态为2Π，即当沿轴向分量的轨道角动量量子数不为1时则需要考虑Λ-双重能级。如NO，NS和OH。

13. 2.3 Non-linear molecular • 通常用Ray’s参量κ来描述不对称度，决定于三个转动常数(A,B,C)： κ=(2B-A-C)/(A-C) • 长椭圆对称:κ=-1, (B=C, e.g. CH3CN); • 扁椭圆对称：κ=1, (B=A, e.g. NH3); • 不对称：κ<<1, (e.g. H2O); • 对于对称陀螺分子：E=BJ(J+1)+(A-B)K2.其中K为总角动量J在对称轴上的投影。e.g. NH3. • 轻微不对称陀螺分子(e.g.H2CO)：

14. 2.4 Hyperfine splitting 考虑核自旋与转动的耦合将会出现超精细结构。 超精细结构会提供光深的信息。 一个普遍假设：每一超精细分量有着相同的激发温度。 对激发态进行分析时，可以不考虑超精细结构。 例子：HCN和OH分子。

15. 3. Radiative rates • 发射速率Aul(u为上能级，l为下能级)： 其中μ为电子偶极矩，一般取μe，并假定电 子偶极矩对同位素是一样的。S为跃迁强度。 • 对于J±1的线性分子有：Sul=gl； • 对于对称陀螺分子，从能级J、K跃迁至J-1、K：Sul=(J2-K2)/J. • 对于非对称陀螺分子，简单的S的表达式不存在。 • 偶极矩的数据依赖于从头计算法的理论计算。

16. 4. Collosional rate coefficient4.1 General considerations 碰撞速率(式一)： 在上面的式子中ncol是被碰撞粒子(如氢气或是He)的柱密度，γ是向下碰撞速率系数，单位为立方厘米每秒。碰撞速率是碰撞散射截面σ的麦克斯韦平均。μ是系统的约化质量，E是质心碰撞能。 另：除了辐射主导(物质与电子碰撞占主导)情况外，一般物质会与H2碰撞。

17. 4.1 General considerations 其中被碰撞物质是He时，γ的关系为； 若被碰粒子是其他情况时也可以用相同的方法来处理。 另有一点，某些元素是可以相互间替换的，如：HCO+与HCS+， O与S，其碰撞速率大小不变。

18. 4.2 Accuracy of adopted rate coefficients • 碰撞速率系数由两个步骤获得： • 决定碰撞系统的相互作用势V； • 计算碰撞动力学。 • 相互作用势能的模型(需用到高级量子化学模型)： • CI：configuration interaction，最精确但耗费时间。 • SCF；perturbation method； • Electron Gas method，已废弃，但仍被用来对天体物理系统做动力学计算。

19. 分子天体物理学中的两大难题 对于碰撞速率系数，只有少部分可以由实验测出，而且大部分实验反映的是许多碰撞事件的平均,正因为如此，碰撞速率系数的测定是分子天体物理中最困难的问题之一。 另一方面，碰撞截面的测定还只是在特定能级间进行，而且获得的只是相对散射截面而非绝对散射截面，因而有效碰撞截面的测定也是分子天体物理中最困难的问题之一。 不过理论与实验的比较还是能够获得一些关于碰撞速率系数不确定度的信息。

20. 4.2 Accuracy of adopted (examples)

21. 4.3 Adopted collisional rate coefficients • 对于线性分子，用于将向下碰撞速率系数(J=Ju→Jl,Ju>Jl)外推至其它温度T的公式为 其中y=Eul/kT,参数a，b，c是用最小二乘 法由拟合每一个J而得到的。 • 精度通常可达20%。 • 对于温度不在能获取碰撞速率系数的温度范围之内的，RADEX不做进一步外推。

23. 6. Extrapolation of collisional rate coefficents6.1 Linear molecules 为推广至高的转动能级，一般采用IOS近似方法，且以基态的碰撞跃迁系数为起算点 其中为Wigner 3-j symbol。该公式条件为：碰撞动能大于转动能级裂距。并且对于高的转动能级精度低。

24. 6.1 Linear molecules 为改进高转动能级计算时的缺陷，一般乘以下式： 有时为推广到更高能级也会用到高次项： 以下是CO的例子：

25. Winger 3-j symbol

26. 6.2. Non-linear molecules 例：SO2： 获得的数据只包含最低的50个能级。 速率系数只取决于温度，与T1/2成正比。 去激发碰撞截面为常数。 若将collision partner由He改为H2则需要乘以1.4.

27. Collisional de-excitation rate coefficients for the lowest 50 states of SO2 summed over all lower levels, calculated from the data by Green for various temperatures.

28. State-to-state de-excitation rate coefficients for SO2as fractions of the total downward rate coefficient as a function of the number of levels by which the transition is changed. The light (coloured) curves are values from Green (1995) at T=25 K for the 10th, 20th, 30th,40thand 50th state above ground. The thick black curve is the normalized mean of the light (coloured) curves,adoptedhere to extrapolate Green’s rate coefficients tohigher-lying levels. The states are labelled in order of increasing energy.

29. The End！