等离子体光谱诊断（完整版）

1、华中科技大学研究生课程考试答题本考生XXX璐考生学号D202177286系、年级2021级类别高电压与绝缘技术系 考试科目等离子光谱诊断 考试日期假设佛克特谱线轮廓逐步计算二价分子和原子的光谱概论 本文讲述的是模拟二价分子和原子的电子跃迁导致光谱线的形成的计算机程序。 这 个程序通过解释在计算中所有的独立自旋量和原子谱线对光谱的作用形成一条光谱。 每一条 谱线的总和在光谱学中被称作近似的 佛克特剖线。 这个程序还能生成一种光学稀薄气体中的 谱线或者同时存在发射和吸收情况下的光谱。 这种方法能够允许计算的范围从一团圆柱状的 低温气体伴随着产生的吸收光谱到一个非等温天然气体的高温自吸收发射光谱。

2、如果需要的 话，计算机能够立即产生光谱图。 这个具有敏感性和广域性的工具能形成光栅摄谱仪或者波 长辐射记。 许多介绍这个程序的多用途性的例子已经被提出。 这些信息被给与那些需要这个 程序的拷贝或者一份包含使用说明的阅读者。说明文章中描述的这个程序能够很好的描述 二价分子和原子的电子跃迁形成的光谱 线。在参考文献 1 中讨论到这个方法能极大的扩展和提高程序的多样性使之能适 用于绝大多数原子和二原子分子跃迁。特别的，这个程序还包含：1. 所有相似跃迁忽略自旋态和两倍 入忽略自旋量和两倍 入指的是多重线强度来自于单一的“有效光线2. 所有垂直跃迁忽略自旋态和两倍 A。3. 2ni2工跃迁，忽略两倍入

3、。4. 原子谱线。5. 当分子由于旋转而分裂的时候谱线循环计算将选择终止。6. 对于同核分子这种选择包含谱线浓度的改变。7. 使用一个近似的 活格模型来描述谱线形态。8. 辐射能量转移通过多分子构造有不同的热力学概率。 从以上的例子可以清楚的了解到任何允许的计算的二原子跃迁都有一定程度的 近似；因此只有121 z跃迁能被准确计算。所有单原子跃迁和 2n-> 2工跃迁都 是近似，因为二倍入的影响实在太小。自旋态和两倍入的限制可以随着模型程序 的重组而改变，但随着而来会导致计算机运行需要更多的时间。依据计算机程序，以下输出选择可用 :1. 作为一种波长效应，光谱浓度能够被制成表格。2. 计算

4、机光谱理论能够 机器绘制。3. 以任意波长间距的完整光强谱也能被绘制出。4. 在计算机光谱中，一个任意波长间隔的谱线可以用生长曲线来计算。5. 计算机光谱结合特定灵敏的工具来产生一个描述构造的量，例如产生电压或 者辐射量。构造带的能传递最大超过 99 条直线局部或者高斯曲线。这个敏感元 件可以认为是一个安装好的波长带通过传递来触发一个辐射计或者通过扫描任 何波长间距的存储单元。 作为一个波长去模拟光谱的效果， 这个敏感局部能被改 变，这样的话这种改变能够被一个具有线性色散的工具所记录下来 光谱仪栅格， 这种模拟光谱也能被描述。COSMIC, Barrow Hall, University of

5、 Georgia, Athens, Georgia 都在使用这个 程序。这个程序包含原始程序卡片组， 一个 18 种二原子系统的光谱学数据的 800 位磁带，一个搭配程序扩展校验的输入和输出。当需要这个程序的时候，参考 ARC-10127 号闪光纸。除此之外，一个使用手册在参考文献 2 中登载。理论一个真实的光谱总是由大量的分子旋转谱线， 原子谱线和连续谱线重叠而形 成。如果要用计算机去准确的模拟自然中的光谱， 每一个辐射源必须考虑到。 现 在的程序试图合理的描述分子旋转谱线和原子谱线。 接下来的讨论被辐射体理论 所局限同时又试图简洁的解释这些。在相关的举例中能找到详细介绍。为了合理的计算光强

6、， 需要考虑气体中各种粒子组成辐射源的内在联系。 这 种内在联系在很多文献中讨论过 3，4 ，和准确的理解辐射传输量纲方程的各种不 同形式 3公式1得出光强梯度，I, W/cm2-y-sr, x是根据气体范围和三大根本 局部；第一局部，考虑自发辐射的影响，是独立的光强。其余的局部包含发射和 吸收光谱,在 x 点时候浓度特殊,这意味着剩下的公式是：公式3说明吸收系数a是指在特定情况下那些用来计算特性光谱浓度的参 数现在指的是真实光谱 。在基尔霍夫等式下,由于自发辐射,这个系数 等于光谱浓度。a=E”B入。现在的程序第一次描述自发发射光谱 E入通过总结在 光谱范围的一些点的所有存在的原子和旋转谱线

7、的光谱分配。 在计算中, 这些点 平均分配在比现在存在的最窄的谱线宽度的小尺寸中。 接着,每一点的吸收系数 被黑体辐射效应简单划分。最后,根据公式 3一个给定的气体深度和入射辐 射源形成真实光谱。 这个程序有很高的实用性是因为一个气体层的真实光谱能够 被存储,然后可以被作为一个新的气体层的入射辐射光谱。 这样就允许从各种复 杂的辐射源来计算真实光谱, 而在这些复杂的辐射源中的每一层有详细的热化学 和热力学性质。如果最后的结果是一个微弱的可见光谱， 它将被自发发射光谱E入直接赋予。 这样公式 3中的约束条件 a?l<< 1。这个过程的起始点是计算自发发射光谱 E胚信息的三个根本局部需

8、要每一 个自旋谱线和原子谱线来完成这个实验， 光线中心波长的位置 心L，由自发发射产 生的总体谱线浓度E,适当的光线属性。这三个局部中的第一局部被程序以一种 数学表达来描述从而产生能使用包括二原子分子和原子的信息， 其他两个局部描 述能任意用于产生自发发生光谱的更深入的信息。 关于一个给定的震荡带的最大 自旋量子数的测定和来自同核二原子分子的光谱线增加也已经给定。1 谱线中心的波长A. 原子。每个原子谱线的中心波长被输入程序使得数学表达变得不怎么需要。 原子谱线波长被制成多种形式的表格 5 。B. 双原子分子。关于自旋谱线的中心波长的一般表达在参考文献中给出。 关于 ?， ?'基点带的

9、转移能够直接输入或者被这个程序来计算自旋常数能量，FJ，指的是一种类似于J的电子水平的函数。关于 F J的表达形式的单次跃迁在参考文献 中被给出。在这个程序中，这种表达形式不仅用于单一跃迁也同样适用于那些自旋分裂 可以被忽略的所有跃迁，在三线态情况下的表达为=0在这些跃迁中，三线态的中心组成局部被作为有效的光谱波长来使用。因此，三线态=0 FJ的表达为，当自旋分裂不能被忽略的时候，就必须使用 FJ的一种更准确的表达来 描述单独水平。例如，频谱从2n-T2的跃迁除非在自旋分裂被考虑的情况下 才能被很好的表达。在程序中，F J在这些跃迁下的表达为，忽略 M咅是： 在参考文献中，增加Dv的表达形式被

10、给出。由于高次项被忽略，因此从 12到 8的独立项被约去。这一点在 781页会更深入讨论。2. 由自发发射产生的谱线强度综合A原子。一个单一原子光线产生的自发发射光谱线总体强度公式可以表达 为：爱因斯坦系数能够在参考文献 5中的表格中找到， 在参考文献9中能找到 N 和O还有其他离子的分析函数。就像在参考文献中讨论的一样，为了防止错误， 公式中的分析函数在使用的时候必须单质计算。B.双原子分子。简单的，一个单一自旋谱线的自发发射光谱的总体强度公 式如下：一般的，在一个电子跃迁的多重谱线中，电子跃迁的瞬时值不能由每个跃 迁产生，需要使用一个平均值。这个重和是包含高能级和低能级的多重谱线水平的所有

11、电子跃迁。当公式 13代人公式 12，结果得出谱线总体强度的数值与尼古拉斯和斯图尔特表 述的一致。这个电子跃迁时刻通常是一个测量量，并且与吸收量一起联系使用。 根据电子跃迁， 吸收量和跃迁时刻能够随着波长而改变。 在文献中， 一些分 子变化量的描述能够找到。公式中的康登因素和光线强度因素都可以理论推出。如果，完整的电子多重谱线结合一个有效能量水平，自旋态和X可以被忽略。 如果自旋量需要考虑，例如2n-T2跃迁，c,?,?必须把du取代。如果自旋态 和糸都要被考虑，那么？?需要被1.0取代。二价原子的局部效应能被公式 16最近似计算。 总和是所有考虑到的温度的电子能级和每个电子能级的震荡水平或者

12、这个 公式描述到达假象的峰值的计算。 为了追求准确性和连续性， 现在的程序中使用 双原子粒子浓度， 这样必须不断校验在等价于那些公式的计算结果所使用的边界 条件。必须注意分别在公式 15和16中电子的特殊性，震荡性，自旋温度， 需要假设每个能量模型有不同的玻尔兹曼奉献构成。光强因素，根据电子跃迁类型和 J.As 在说明中提及的自旋多重度，三种类 型的电子跃迁程序都有考虑，但是自旋多重度仅仅考虑2n-J工跃迁。在每个类型中，光线强度因素的表达的讨论在下文提出。1. 平行跃迁人二。这些被高能级P和R还有其他低能级描述的跃迁 通过增加量子数而快速的消逝。 跃迁中的自旋多旋度在程序中被忽略， 强度因素

13、 按照1工工分配。通过这种分配，低能级被认为强度是 0。然而，流失的强度按照P,R层分 布，因为总强度必须遵守下面的准那么。有意思的是，通过限制高能级 K值，这 种强度分配会准确的忽略自旋。2. 垂直跃迁4A= 1。这些跃迁被高能级Q层和稍微低能级P层和R层。 如果电子自旋结果为 0，他们就是唯一的层。如果计算的结果不是 0，那么说明 其他低能级层存在， 但是随着量子数的增加而快速的消逝。 仅仅程序中应用的唯 一强度因素来自参考文献 。3. 2 n-J工跃迁。这些跃迁说明12中可能的层，从 J=0 ± 1， K=0 ±1，±2 都必须要算。如果 ?没有包含在自旋能

14、量公式，来自于四层的自旋谱 线与其他四层有着一样的波长，那些数目可见的层会降到 8. 在总结两个别离强 度因素的时候，两倍谱线的有效强度被发现。三谱线类型。这种谱线类型被 佛克特剖线近似假设。参考文献 15给出佛克特剖线的闭合形 式的近似，导致其本身快速计算而且一般其准确率能到达 3%或者更低。这种近 似在程序中被使用，近似的公式在下面列出。CDv, CDl，佛克特剖线都是给定，一个来自谱线中心的谱线宽度特殊量被加进 自发发射光谱。这种特殊性能被两种方式中的一种构成。 当真实光谱被过量吸收， 我们感兴趣饿范围增加到任意光谱线的中心， 最近谱线的远端也必须考虑。 由于 这个原因，到这些边缘的计算

15、距离为了合理的模拟，必须在程序中具体的输入。 因此，从一些数量的谱线中心在一个大范围安置目标点会显著增加计算机的计算 时间。在那些包含大量闭合空间曲线或者可见稀薄空气的光谱领域， 谱线的远端就 不重要了。 在这种情况下， 就可以忽略掉目标距离的具体输入数据， 程序将根据 具体的谱线类型设置距离。对于一个完全的高斯剖面，?/ ?被?忽略从光线中心往外的 3 个光线长度。 对于洛仑兹或者沃伊特剖面， 光线的边际更强所以是从 光线中心的 5个长度。在最坏的情况下，接近 5%的自发谱线强度超过目标点， 这样在计算中不会考虑。四自旋量子数的最大值 由于这个程序计算双原子线性光谱， 必须在计算中考虑悬着转

16、动谱线量子 数的准那么。使用者通过规定最大的和最小的自旋量子数来获得每个振荡跃迁从 而做出这个选择， 与此同时， 这个程序能在自旋导致分子分裂之前， 计算出最大 可能的自旋量子数。这种找出最大自旋量子数的方法在参考文献 6， 12中也有相 似的介绍。自旋系统中有效电位能量是由旋转系统和莫尔斯势函数总和。 五双原子分子光谱同核分子光谱带在谱线强度中呈现出一中特性交替。 在一些情况下， 由于核 自旋导致多重性的改变或者自旋水平的统计权重的改变从而使得其他谱线完全 消失，这在参考文献 6中有相关描述。这个规律可以用数学式表达在程序中被使 用。程序校验和计算实例一个计算机程序显然只能和那些根底的和模型

17、的数序表达能转化进计算机 的理论模型做的一样好。 首先的局部是定义模型的范围和讨论在应用中的重要假 设。这个局部的主要目的是按照需要定义计算机程序的光谱和对于它的重要性和 多样性做一个简单的计算解释。在开场程序之前， 需要明白的是计算机的运行时间是由产生自发发生光谱所 需要的时间决定的。这个时间能被公式推出。运行时间包含光线总数光线截至点的线宽量叫/ 分钟 过去工作中使用的计算机是IBM7094，常数是5.8 >10-5分钟/条。一些简单 计算的运行时间将在下面给出。1. 光线中心波长通过实验与计算机程序计算结果想比较，得出公式 5准确的表述波长图一中显示的是计算机计算红外跃迁的谱头波长

18、与参考文献 17中的测量结果比 照的具体结果。 每个跃迁的第四谱头波长的计算数据从可用的数据中祛除， 这个 区别以振动量子数的图标表示出来。量子数带由增长的波长决定包含的范围是 4800到11000艾。参考文献 6， 18提供计算机使用的光谱常数。一般来说，误差 在 0.5 埃或者小于一个波长宽度范围。振动量子数能证明公式5，6，10的推导是正确的。的光谱常数和测量波长都在程序中经过验证。 或者， 如果测得新的波长， 这 个程序也能很好确实定适宜的新光谱常数。2. 弱可见光强准确地表述光强取决于自发发射，公式 12，这个光强通过把 1cm 路径长 度的弱可见光强带与 Q 理论的值相比较来检验。

19、 Q 理论不是为了产生详细的光 带构造，而是对于从一个到另一个带序的弱可见光强中能得到有用的结果。这种比较在图 2中显示，表3显示的是在二氧化碳和氮气各占 50%的混合气 体加热到5000k情况下的主要光谱特征，图2中上曲线表示现在计算的结果，而 下曲线那么是 Q 理论的结果。现在的计算显示得到的数据与艾姆斯常规机器中 得到的数据一样。 用这种技术， 每个计算点与它附近的直线相连。 在谱带系统的 重叠局部， Q 理论通过把现在有用的系统相加来现在的计算结果做比较从而制成 表格。谱带系统认为，紫色B2x2，红色A2n x22:,广域?n?n? 跃迁，具体类型在参考文献 19中提到。无量纲跃迁平方

20、，是对每个谱带系统的连续假设，使用值分别为各自 0.56， 0.52和 3.56。康登重叠因数在参考文献 20中提到。红色的光谱分析在参考文献 6， 18 ，21提到。关于 ?2?的光谱分析在参考文献 22提到。线性计算中的光线类型被指 定为 ?=2.0 埃的高斯剖线。 接近 7774埃的臭氧分子也可以进展线性计算， 但是 在这个情况下， 光线太弱而不能在红谱构造中显示出来。 这个例子解释了在二原 子气体中大量自旋谱线是 68100条独立自旋曲线和三原子谱线。 由于这个原因计 算机计算时间在 50分钟。表三给出典型的弱可见光与三系统变化谱带的比较。每个分支在 150次自 旋谱线后终止。这是因为

21、根据计算说明几乎没有超过了 150次以后没有更多的能 量存在。相比较而言， Q 能够应用于更限制的波长间隔。通过比较而选出的谱带 依靠波长的不同来区别目的是为了取满广泛独立的谱线空间。 积分强度的比较大 体上是好的， 最主要的区别发生在系统。 这个结果是合理的， 因为红谱系统构造 很复杂而 Q 理论得到的结果在简洁性上不如这个系统的结果，为了得到比现在的计算更简洁，有用的结果。在表三显示现在的计算与参考文献 9 中提到的三氧 原子系统产生的弱强度光谱相吻合，同时也可以推出公式 11是合理的。 3计算单独高温气体和分子曲线增长的真实光谱。通过使用已经经过计算的自发发射光谱得到结果， 而辐射跃迁通

22、过线性计算1工来表达。他们的计算不是由辐射引起的，因为这个原因，修改很小而且不 能被忽略。 这个计算仍然需要福伊特剖线进展检验， 而且这个程序能够通过篦光 谱仪产生一个同样的光谱。图三显示这个不需要仪器扩展的真实光谱，在计算中，从参考文献6和参考文献23弗兰克-康登因子得到分光光谱。电子跃迁的无穷小量为3.31 X0-2，宽度 是10.3厘米，气体温度是3600K，同物种浓度为1.38 X016part/cco输入光线宽 度 翊=0.038埃，*0.103埃，福伊特宽度 3=0.116埃。这些输入数据与参考文 献【23】中的数据一样。图三中，0, 0R分支的温度最大值曲线到达 3600K。因

23、为这个原因，计算机运行时间是 2分钟。图三b显示图三a的光谱图，它是通过一个1.2埃宽具有顶点反响集 合的矩形光栅功能的光谱仪描绘出来。 文章中表示可一次可以从参考文献 23中引 用3600K，黑体辐射曲线，参照合成光谱来读取这些点。这种符合度很好，只是 在谱带前端的很远位置才有微小的不同。以上的参照是为了那些没有涉及自旋态的12 -12跃迁。然而在一些跃迁中， 自旋态是存在的， 忽略它会在计算大量自吸收光谱的时候产生大量的错误。 修正 自吸收光谱明显的方法就是在计算程序中做出重要修改， 使得每一个独立谱线符 合1口工。但是不幸的是会增加计算机的运行时间。另一种方法是找到一种有 效的谱线类型，

24、 而这种谱线会导致其强度的最大值与正确的多重谱线近视。 两个 限制性的原因必须考虑到：1当自旋分裂光谱比光线宽度还短的时候， 这些光 线将变成一个单一的谱线。 当自旋态被忽略的时候， 在计算机程序中是一种典型 的情况。2当自旋分裂光谱与光线宽度近似时。在这样的情况下，光线不会明 显的重叠， 通过使用一个自旋多重度的一般宽度叠加来比较一个 “真实光线宽 度从而使得光线中心的强度最大值最逼真。 这样可以有效的模糊光线与真实的多 重谱线的误差从而不需要增加光强。 由于这个原因，选择一个有效谱宽就更困难。从两个限制性的原因中， 评估误差。 成长曲线在 7000K 的计算是为了在0， 0和 1，1谱带前

25、端的强度下降梯度。这个结果于归因于三次分裂的费尔班计 算结果相比较。在图四a中，这个比较用完全重叠谱线表示出来， 在图四b 中，这个比较用完全不重叠谱线表示出来。 在大量自吸收情况下， 费尔班结果不 赞同计算中使用完全重叠谱线。 然而，完全不重叠光线计算需要三次使用真实的 光线宽度，费尔班理论中的具体光线宽度很好的符合他的结果。 最大误差为 6%。 与高斯和洛伦茨光线谱宽相联系的参数a由图四a中标志性的表格定义。计 算机需要运行大约四分钟来计算C2曲线生长。下面的例子通过新的部件将更深入的介绍这个程序的可能应用。 结果，现在 没有参考的数据。图5第一次显示这些应用，解释 CO4+1 2 1 2

26、发射谱线途径，作为辐 射气体的几何深度像黑体限制一样不断增长。从 1440埃到 1700埃范围内的谱带 都是需要研究的。 CO 放射具体来说就是一定浓度， 密度下?2?和氩气各 50%被加 热到6000K。在这个计算中，所有的光线必须以 =1.0埃和®=0埃计算。由 于这是个单峰跃迁，所以，不需要自旋量子数的修正。更深入的，对于所有的 ?， ?"跃迁的无穷小跃迁量的取值为 0.76，这个 值与电子吸收光谱常数 f(0.15 到 1545 埃)相对应。参考文献 25提到弗兰克 - 康登因 子，参考文献提到光谱内容。值得注意的是，在表中，接近 6000K的黑体辐射 曲线是一个宽

27、度为5cm的真实光谱。同样的，注意1.0cm的路径长度，在一些地 方，弱可见光谱事实上超过黑体辐射的限制。相似的计算在图5中提出，这个计算被用于产生图6CO 4+生长曲线。关于0，0与4，3谱带首段光强梯度计算将做成以 0.0162 埃宽度的高斯 曲线，这个计算是根据初等多普勒理论和光线宽度参数a，1.0和0.01。a的取值因为超过文章讨论范围因此取值是任意的。 对于生长曲线的计算所需要的其他参数 与图5中所示一样。由于CO(4+)谱带空间严密并且属于红外光谱，从谱带尾 部到谱带头部 0， 0自发发射光谱的强度显著。在计算中，这些包括谱带 200 埃到0， 0 的紫色。在一个附加的计算中，确定

28、的一体化空间的自发发射光谱 的下一个谱带有小于 0.05%超过这个点。计算机把每个曲线生长运行出来需要差 不多 12 分钟。图7 说明一个纯洁的高温气体的真实光谱，而这个光谱通过线色散工具 来决定是否需要附加局部。原始气体是一个温度为6000K，物种浓度为1.0 >015个 /升的 1cm 厚度的气体。一个浓度是 2.0 >1021个/升任意高原子氮气被认为能产生两 条氮气谱线在红外区域临近 10880 埃。参考文献 9中，计算机计算所需要的数据 来自于氮气谱线的自发发射光谱的一体化强度。在参考文献1中红色光谱曲线与其他谱线有区别。对于和N跃迁的高斯谱线宽度，被作为温度6000K的多普勒谱 线宽度输入计算机。 两者的谱线宽度参数假设为 1.0埃。图二中