-
Nitric oxide (NO) is one of atmospheric molecules of interest and has attracted considerable attention due to its important role in the chemical process taking place in a flow field of hypersonic vehicle, in which the thermodynamic properties are required in the calculation of the aerothermodynamic flow field. Moreover, the total internal partition function is the key to calculating the thermodynamic properties of high-temperature gases. For diatomic molecules, according to the product approximation, the total internal partition function is split into three parts: electronic, vibration and rotation partition function. In this paper, by using the quantum statistical ensemble theory based on some classical thermodynamic and statistical formulae, the thermodynamic properties of NO are analyzed and discussed. Firstly, in order to obtain an accurate energy of molecule, the variational algebraic method (VAM) is employed to calculate the full vibrational energy, the resultis in good agreement with the experimental result and thus yielding the realistic predictions of the unobserved higher vibrational energy that converges to the dissociation limit. Secondly, an attempt is to use the full VAM vibrational energy, the Rydberg-Klein-Rees (RKR) vibrational energy, the simple Harmonic oscillator (SHO) model and the quantum-mechanical vibrational energy obtained by the multiconfiguration self-consistent-field (MCSCF) to calculate the vibrational partition function. Then, with the rotational contributions from the Müller-McDowell formula, the internal partition function can be determined by combining the product of electronic, vibration and rotation partition functions. Thirdly, according to the thermodynamic and statistical formulae, it is easy to calculate the internal energy, entropy and heat capacity for the NO molecule in a range of 1000-5000 K. Comparison of different calculated heat capacities with the experimental ones reveals the heat capacity, of which vibrational contributions determined by the full VAM vibrational energy accord better with the experimental ones, with the maximum relative error being no more than 2.4%, whereas it can be seen that those thermodynamic results evaluated from the SHO model attest to a failure for the summation of infinite vibrational energy. The thermodynamic results of NO may have proper applications in areas that can be of great importance in theoretical and (or) experimental aspects. 1. 引 言
在行星大气中扮演重要角色的NO气体是造成酸雨、光化学烟雾和臭氧层破坏等大气污染的主要成份之一[1-3], 且NO的分子结构与性质在生物、原子分子物理和药物学等领域一直受到人们的广泛关注[4-9]. 近年来, NO气体宏观热力学性质也成为了人们感兴趣的研究课题之一. 早在1987年, Jaffe[10]将内能划分为振动、转动和电子贡献计算了NO气体的热容, 并通过这些热容数据计算了气体混合物(含NO)的热力学温度和焓. 2005年, Capitelli等[11]获得了与火星大气相关的NO, NO2, N2O等成份的热力学性质, 为高温混合气体热力学模型的研究提供了重要的数据支撑. 2009年, Babou等[12]在多温度模型框架下, 采用双温度分区计算了NO分子在高温和非平衡介质中的热力学函数、配分函数、平均能量和旋转比热. 2012年, 邓伦华等[13]利用光外差-浓度调制吸收光谱技术测量了NO分子在12530—12850 cm–1波段的吸收光谱, 在该波段内标识出了324条跃迁谱线, 随后采用标准的
$^4{\Sigma ^ - }{ - ^4}{\Pi _i}$ 哈密顿量模型并借助最小二乘法拟合其中的267条谱线, 获得了体系的主要分子常数和精细结构分子常数.基于统计物理学的微观振动信息人们也能得到与双原子分子振动相关的宏观热力学量[14]. 对于NO分子微观振动信息的理论研究主要基于1950年Huber和Herzberg[15]建立的计算双原子分子振动能级表达式. 早在1975年, Billingsley[16]便采用优化价态构型多构型自洽场方法得到了基态NO分子的电子偶极矩函数和势能曲线. 2006年, Reddy等[17]依靠实验技术用改进的RKR (Rydberg-Klein-Rees)方法构建了基态NO和NO+的势能曲线, 并采用五参数(Hulburt-Hirschfelder, H-H)函数拟合获得了这些电子态的离解能. 2018年, Qin等[18]使用RKR方法获得了NO分子低阶电子态的势能曲线, 利用该势能曲线通过求解薛定谔方程确定了振动能级, 计算了NO分子平衡和非平衡状态的热力学性质.
以上研究为了解NO分子电子态的势能曲线和振动信息做出了巨大贡献. 然而, 实验上现阶段的光谱测量技术一般只能获得体系低激发振动量子态的物理信息, 对于大多数NO分子电子激发态而言, 高激发振动量子态尤其是处于离解区的数据却很难得到. 由于这一客观事实的存在, 人们通常采用简单的谐振子模型(simple harmonic oscillatormodel, SHO)来计算体系的微观振动信息, 再借助量子统计系综理论来获得系统的宏观热力学量. 如此以来, 用SHO模型计算得到的一部分实际上不存在的虚高能级也参与了统计计算, 从而导致最终的热力学计算结果往往出现较大误差. 为了获得NO分子电子态可靠的完全振动能级从而进一步研究其微观振动对宏观热力学性质的影响, 本文基于课题组前期建立的变分代数法(variational algebraic method, VAM)和能精确计算分子离解能的新解析表达式[19], 研究获得了基态NO分子的完全振动能级, 随后采用量子统计系综理论[14]计算了NO分子的热力学性质, 并将计算获得的理论结果与实验数据进行了分析和比较.
2. 理论与方法
量子统计系综理论的重要基础是获取系统正确的配分函数, 再将摩尔内能、摩尔熵、摩尔振动热容表示为含有配分函数的表达形式, 从而确定系统的全部热力学性质[14]. 基态的NO分子在不考虑电子、振动和转动的相互作用时可由乘积近似模型得到其配分函数[20]:
$ Q(T) = {Q_{\rm{n}}}{\rm{(}}T){Q_{\rm{e}}}(T){Q_{\rm{\upsilon }}}(T){Q_{\rm{r}}}(T), $
(1) 式中, T为系统的开氏温度,
${Q_{\rm{n}}}(T),\; {Q_{\rm{e}}}(T),\; {Q_{\rm{\upsilon }}}(T), $ Qr(T)分别表示核配分函数、电子配分函数、振动配分函数以及转动配分函数. 电子配分函数为$ {Q_{\rm{e}}}(T) = \sum\limits_i^{} {{g_{\rm{e}}}(i){{\rm{e}}^{[ - {c_2} \cdot {E_{\rm{e}}}(i)/T]}}}, $
(2) 式中,
${g_{\rm{e}}}(i)$ 表示电子能级i的简并度,${c_2}$ 是第二辐射常数,${E_{\rm{e}}}(i)$ 表示第i个电子能级能量, 通常取基态能量${E_{\rm{e}}} = 0$ . 另外, 本文中的NO为异核双原子分子, 其核配分函数等于原子核基态的简并度${d_{\rm{s}}} = 1$ . 故NO的配分函数$Q(T)$ 可表示为$ Q(T) = {Q_{\rm{\upsilon }}}(T){Q_{\rm{r}}}(T). $
(3) 采用正则量子系综, 体系的振动配分函数可定义为
$ {Q_{\rm{\upsilon }}}(T) = \sum\limits_{\upsilon = 0}^{{\upsilon _{{\rm{max}}}}} {{g_{\rm{\upsilon }}}{{\rm{e}}^{ - \beta {\kern 1pt} {E_{\rm{\upsilon }}}}}}, $
(4) 式中,
${g_{\rm{\upsilon }}}$ 表示能级的简并度, 此处${g_{\rm{\upsilon }}} = 1$ ,${E_\upsilon }$ 是分子内部的振动能级,$\beta = 1/(kT)$ , 玻尔兹曼常数k = 1.380062 × 10–23 J/K.由于NO具有两个电子的双重分裂影响, 其转动配分函数是
${}^2{\Pi _{{1 / 2}}}$ 与${}^2{\Pi _{{3 / 2}}}$ 两个电子态的和[21], 可表示为$ {Q_{\rm{r}}}(T){\rm{ }} = \sum\limits_{{\rm{all\ states}}\ i} {{g_i}} {{\rm{e}}^{ - hc{E_i}/(kT)}}{\rm{ }}, $
(5) 式中,
${g_i}$ 表示能级的简并度, h为普朗克常数, c为光速, 再根据McDowell建立的解析公式[22], NO的转动配分函数可表示为$ {Q_{\rm{r}}}{\rm{(}}T) = \sum\limits_{J = 0} ( 2J + 1){{\rm{e}}^{ - \frac{{hc}}{{kT}}{E_{\rm{r}}}}}, $
(6) 其中Er是转动能量, 由转动光谱常数B和转动量子数J决定:
$ {E_{\rm{r}}} = BJ(J + 1). $
(7) 将欧拉-麦克劳林公式应用于(6)式可以等价为如下的渐进的展开式[22]:
$ {Q_{\rm{r}}}(T)\sim \frac{T}{{{\varTheta _{\rm{r}}}}}\left[ {1 + \frac{1}{3}\frac{{{\varTheta _{\rm{r}}}}}{T} + \frac{1}{{15}}{{\left( {\frac{{{\varTheta _{\rm{r}}}}}{T}} \right)}^2} + \frac{4}{{315}}{{\left( {\frac{{{\varTheta _{\rm{r}}}}}{T}} \right)}^3}} \right]. $
(8) (8)式中
${\varTheta _{\rm{r}}} = hcB/k$ , NO的转动光谱常数B = 1.70488847 cm–1 [23].对于常见的双原子分子组成的宏观系统, 内能是系统能量在所有微观状态上的平均值, 摩尔内能可表示为
$ U = {N_{\rm{A}}}k{T^2}\frac{\partial }{{\partial T}}\ln Q(T), $
(9) (9)式中
${N_{\rm{A}}}$ 为阿伏伽德罗常数. 统计系统的摩尔熵表示为$ S = {N_{\rm{A}}}k\left[ {\frac{\partial }{{\partial T}}\left( {T\ln Q(T)} \right)} \right]. $
(10) 统计物理中可观测的物理量是摩尔热容, 可表示为
$ C_{\rm{m}}^{} = {N_{\rm{A}}}k\frac{\partial }{{\partial T}}\left[ {{T^2}\frac{\partial }{{\partial T}}\ln Q(T)} \right]. $
(11) 根据上述公式通过对系统的所有能级求和可计算得到宏观热力学量. 然而, 由于实验条件的限制, 对于大多数双核体系很难获得其电子态包含高激发态在内的完全振动能级, 因此, 人们通常采用简单谐振子模型(SHO)来描述振动运动对无限个振动能级进行求和获得热力学性质, 等同于认为最大振动量子数
${\upsilon _{{\rm{max}}}} \to \infty $ , 得到配分函数的解析形式:$ Q_{{\rm{SHO}}}^{{\rm{cal}}}(T) \!=\! Q_{\rm{\upsilon }}^{{\rm{SHO}}}(T){Q_{\rm{r}}}(T) \!=\! \frac{{{{\rm{e}}^{ - \beta {\kern 1pt} \hbar {\omega _{\rm{e}}}/2}}}}{{1 \!-\! {{\rm{e}}^{ - \beta {\kern 1pt} \hbar {\omega _{\rm{e}}}}}}}{Q_{\rm{r}}}(T), $
(12) 其中
$Q_{\rm{\upsilon }}^{{\rm{SHO}}}(T)$ 表示谐振子模型的振动配分函数,$\hbar = {h / {(2{\text{π}})}}$ ,${\omega _{\rm{e}}}$ 为光谱常数, 通过(12)式可得到相关的热力学量的解析解. 实际上, 真实双原子分子各电子态的最大振动量子数${\upsilon _{{\rm{max}}}}$ 并非是无穷大的, 且(12)式对高于离解能实际上不存在的虚高能级进行统计求和也是不合理的, 因此不应采用简单谐振子振动模型来计算宏观系统的热力学性质.为了获得某分子电子态的完全振动能级, 孙卫国课题组在前期的研究工作中提出了变分代数方法(VAM), 该方法的建立是根据Born-Oppenheimer近似, 利用二阶微扰理论求解薛定谔方程, 获得如下双原子分子体系振动能级的表达式[24]:
$ \begin{split} {E_{\rm{\upsilon }}} =\; & {\omega _0} + ({\omega _{\rm{e}}} + {\omega _{{\rm{e0}}}})\Big( {\upsilon + \frac{1}{2}} \Big) - {\omega _{\rm{e}}}{x_{\rm{e}}}{\Big( {\upsilon + \frac{1}{2}} \Big)^2}\\ & + {\omega _{\rm{e}}}{y_{\rm{e}}}{\Big( {\upsilon + \frac{1}{2}} \Big)^3} + {\omega _{\rm{e}}}{z_{\rm{e}}}{\Big( {\upsilon + \frac{1}{2}} \Big)^4} \\ & + {\omega _{\rm{e}}}{t_{\rm{e}}}{\Big( {\upsilon + \frac{1}{2}} \Big)^5} + {\omega _{\rm{e}}}{s_{\rm{e}}}{\Big( {\upsilon + \frac{1}{2}} \Big)^6}\\ & + {\omega _{\rm{e}}}{r_{\rm{e}}}{\Big( {\upsilon + \frac{1}{2}} \Big)^7} + \cdots, \\[-17pt] \end{split} $
(13) 其中
$({\omega _0}, {\omega _{\rm{e}}}, {\omega _{{\rm{e0}}}}, {\omega _{\rm{e}}}{x_{\rm{e}}}, {\omega _{\rm{e}}}{y_{\rm{e}}}, {\omega _{\rm{e}}}{z_{\rm{e}}}, {\omega _{\rm{e}}}{t_{\rm{e}}}, \cdots )$ 是双原子分子的振动光谱常数, υ是振动量子数,${E_{\rm{\upsilon }}}$ 是振动能级. 为了抵消实验上可能出现的测量误差, 在每个能级节点上添加了一个微小变分项$\text{δ} {E_{{\rm{\upsilon }}i}}$ [19], 并将(13)式改写为$ {{AX}} = {{E}} + \text{δ}{{E}} = {{E}}', $
(14) 其中
${{A}}$ 为系数矩阵,${{X}}$ 为双核分子光谱常数,${{E}}$ 为振动能级, 分别为$ \begin{split} &{{A}_{{\rm{\upsilon }}k}} = {\left( {\upsilon + \frac{1}{2}} \right)^k},\;\; {{X}} = {\left( {\begin{array}{*{20}{c}} {{\omega _0}} \\ {{{\omega '}_{\rm{e}}}} \\ { - {\omega _{\rm{e}}}{x_{\rm{e}}}} \\ {{\omega _{\rm{e}}}{y_{\rm{e}}}} \\ \vdots \end{array}} \right)_{m \times 1}}, \\ &{{E}} = {\left( {\begin{array}{*{20}{c}} {{E_{{\upsilon _1}}} + \text{δ} {E_{{\upsilon _1}}}} \\ {{E_{{\upsilon _2}}} + \text{δ} {E_{{\upsilon _2}}}} \\ {{E_{{\upsilon _3}}} + \text{δ} {E_{{\upsilon _3}}}} \\ \vdots \\ {{E_{{\upsilon _m}}} + \text{δ} {E_{{\upsilon _m}}}} \end{array}} \right)_{m \times 1}}, \end{split} $
(15) 其中
${\omega '_{\rm{e}}} = {\omega _{\rm{e}}} + {\omega _{{\rm{e0}}}}$ . 随着光谱实验装置的不断改进, 实验上产生的误差不会太大, 有约束条件$ \text{δ} {E_{{\upsilon _i}}} \to 0. $
(16) 对于某双原子分子电子态, 只需从该电子态已测量获得的n个有限实验振动能级数据中选取m个能级, 由于不同分子的不同电子态的物理性质不同, 能级的展开形式应该存在一定的差异性, 因此贡献项的个数选择也有所不同, 根据不同电子态的能级性质自由选取贡献项的个数(即
$ m = 3, 4, 5,\cdots $ )代入(14)式中进行数值计算, 就可获得$N = C_n^m$ 组振动光谱常数$({\omega _0}, {\omega '_{\rm{e}}}, {\omega _{\rm{e}}}{x_{\rm{e}}}, {\omega _{\rm{e}}}{y_{\rm{e}}}, {\omega _{\rm{e}}}{z_{\rm{e}}}, {\omega _{\rm{e}}}{t_{\rm{e}}}, \cdots )$ . 再利用以下的物理判据选出一组最佳的振动光谱常数${X'}$ :$ \frac{{{\rm{d}}{E_{\rm{\upsilon }}}}}{{{\rm{d}}\upsilon }}{\Big|_{\upsilon = {\upsilon _{{\rm{max}}}}}} = 0, $
(17) $ \Delta {E_{{\upsilon _{{\rm{max}}}},{\upsilon _{{\rm{max}}}} - 1}} = {E_{{{\rm{\upsilon }}_{{\rm{max}}}}}} - {E_{{\upsilon _{{\rm{max}}}} - 1}} \to \;{\rm{small}}\;{\rm{enough}}, $
(18) $ D_{\rm{e}}^{{\rm{exp}}} - {E_{{\upsilon _{{\rm{max}}}}}} \to \;{\rm{small}}\;{\rm{enough}}, $
(19) $ \Delta E({\rm{exp}},{\rm{cal}}) = \sqrt {\frac{1}{m}{{\sum\limits_{\upsilon = 0}^{m - 1} {\left| {{E_{\upsilon,{\rm{exp}}}} - {E_{\upsilon,{\rm{cal}}}}} \right|}^2 }}} \to 0, $
(20) $ D_{\rm{e}}^{{\rm{cal}}} \cong {E_{{\upsilon _{{\rm{max}}}}}} + \frac{{\Delta E_{{\upsilon _{{\rm{max}}}},{\upsilon _{{\rm{max}}}} - 1}^2}}{{\Delta E_{{\upsilon _{{\rm{max}}}},{\upsilon _{{\rm{max}}}} - 2}^{} - \Delta {E_{{\upsilon _{{\rm{max}}}},{\upsilon _{{\rm{max}}}} - 1}}}}, $
(21) $ {E_{{\upsilon _{{\rm{max}}}}}} < {\rm{D}}_e^{\rm {cal}} \leqslant D_e^{\rm {exp}}, $
(22) $ 0 < {\rm{Error}} = \frac{{D_{\rm{e}}^{{\rm{exp}}} - D_{\rm{e}}^{{\rm{cal}}}}}{{\Delta {E_{{\upsilon _{{\rm{max}}}},{\upsilon _{{\rm{max}}}} - 1}}}} \leqslant 1. $
(23) (17)—(20)式是振动收敛条件, (21)式是VAM方法给出的离解能计算公式, (22)式表示计算出的振动能级正确收敛于离解能, 得到最优的一组振动光谱常数后, 将其代入振动能级的精确表达式(13)式中, 就能够对体系未知的振动能级进行预测. (17)—(23)式中, “exp”表示实验已有的数据, “cal”表示VAM计算获得的数据,
$D_{\rm{e}}^{}$ 表示体系的离解能,$D_{\rm{e}}^{{\rm{cal}}}$ 表示用最高的三个振动能级代入(21)式计算得到的离解能.综上所述, 对于基态NO分子, 加入微小的变分项抵消实验振动能级误差并利用代数方法求解其完全振动能级
$\left\{ {E_{\rm{\upsilon }}^{{\rm{VAM}}}} \right\}$ 和离解能$D_{\rm{e}}^{{\rm{cal}}}$ 的方法称为变分代数法(VAM). 获得的VAM振动配分函数则为$Q_{\rm{\upsilon }}^{{\rm{VAM}}}(T)$ , 因此(11)式中的摩尔热容量表示为$ C_{{\rm{VAM}}}^{{\rm{cal}}} = {N_{\rm{A}}}k\frac{\partial }{{\partial T}}\left( {{T^2}\frac{\partial }{{\partial T}}{\rm{ln}}\left[Q_{\rm{\upsilon }}^{{\rm{VAM}}}(T){Q_{\rm{r}}}(T)\right]} \right). $
(24) 相似地, 将(24)式中的
$Q_{\rm{\upsilon }}^{{\rm{VAM}}}(T)$ 替换成$Q_{\rm{\upsilon }}^{{\rm{RKR}}}(T)$ ,$Q_{\rm{\upsilon }}^{{\rm{MCSCF}}}(T)$ ,$Q_{\rm{\upsilon }}^{{\rm{SHO}}}(T)$ 则可以获得RKR热容值$C_{{\rm{RKR}}}^{{\rm{cal}}}$ , 偶极子的多构型自洽场 (multi- configuration self-consistent-field calculation of the dipole moment function, MCSCF)热容值$C_{{\rm{MCSCF}}}^{{\rm{cal}}}$ 和谐振子热容值$C_{{\rm{SHO}}}^{{\rm{cal}}}$ .如此一来, 用VAM方法求得的热力学量的具体表现形式就不同于谐振子模型无穷项求和的形式, 减小了对实际上不存在的虚高能级进行统计求和带来的误差.
3. 应用与讨论
首先应用VAM方法计算了基态NO分子的完全振动能级. 表1列出了基态NO分子VAM完全振动能级
$\left\{ {E_{\rm{\upsilon }}^{{\rm{VAM}}}} \right\}$ 并计算了相对于实验上获得的有限精确振动能级$\left\{ {E_{\rm{\upsilon }}^{{\rm{exp}}}} \right\}$ [17]的绝对误差. 如表1所列, VAM方法选出了6条($\upsilon = 0, {\rm{1}}, {\rm{3}}, {\rm{5}}, {\rm{6}}, {\rm{7}}$ )最能反映基态NO分子振动行为的实验振动能级(黑体标出). 利用这6条数据计算得到了收敛于离解能的完全振动能级, 并且在低激发振动能级区域与实验振动能级符合得很好. 表1还列出了MCSCF方法得到的振动能级[16]及其与实验振动能级的绝对误差.$\upsilon $ $E_{\rm{\upsilon }}^{{\rm{exp}}}$[17] $E_{\rm{\upsilon }}^{{\rm{MCSCF}}}$[16] $E_{\rm{\upsilon }}^{{\rm{VAM}}}$ $E_{\rm{\upsilon }}^{{\rm{exp}}} - E_{\rm{\upsilon }}^{{\rm{MCSCF}}}$ $E_{\rm{\upsilon }}^{{\rm{exp}}} - E_{\rm{\upsilon }}^{{\rm{VAM}}}$ υ $E_{\rm{\upsilon }}^{{\rm{VAM}}}$ 0 948.50 948.60 948.50 –0.10 0 26 40240.70 1 2824.50 2824.60 2824.50 –0.10 0 27 41310.84 2 4627.30 4672.40 4672.27 –45.10 –44.97 28 42340.94 3 6491.90 6492.10 6491.90 –0.20 0 29 43329.70 4 8283.50 8283.90 8283.43 –0.40 0.07 30 44275.74 5 10046.90 10047.80 10046.90 –0.90 0 31 45177.59 6 11782.30 11783.70 11782.30 –1.40 0 32 46033.69 7 13489.40 13491.70 13489.60 –2.30 –0.20 33 46842.41 8 15171.80 15168.73 34 47602.00 9 16823.90 16819.60 35 48310.62 10 18448.00 18442.06 36 48966.35 11 20044.00 20035.93 37 49567.15 12 21611.80 21600.99 38 50110.90 13 23151.30 23136.97 39 50595.35 14 24662.40 24643.56 40 51018.15 15 26144.80 26120.38 41 51376.86 16 27598.50 27567.04 42 51668.92 17 29023.20 28983.04 43 51891.64 18 30367.88 44 52042.23 19 31720.95 20 33041.62 21 34329.18 22 35582.84 23 36801.78 24 37985.07 25 39131.73 $D_{\rm{e}}^{{\rm{exp}}}$ 52155.68 $D_{\rm{e}}^{{\rm{cal}}}$ 52155.68 注: $E_{\rm{\upsilon }}^{{\rm{VAM}}}$中表示VAM计算所需的已知实验振动能级用黑体标出. 将VAM方法和MCSCF方法计算的振动能级与实验能级进行比较, VAM的误差明显小于MCSCF的误差. 在实验上获得的有限精确振动能级区域, 除振动量子数为2的能级外, VAM振动能级与实验值的误差均在[-0.2, 0.067] cm-1范围内. 振动量子数为2的能级误差很大, 实验上获得的数据为4627.300 cm–1, VAM计算值为4672.274 cm–1, MCSCF的计算值为4672.000 cm–1, 两种方法误差均达到45 cm–1左右. 王建坤等[25]用多组态相互作用理论(multireference configuration interaction, MRCI)以cc-pVDZ为基组得到的理论计算结果也在4672 cm–1附近. 更直观地, 从图1可以看出, 使用(18)式和(19)式来判断VAM振动能级的数值结果在其离解区的收敛性是可靠的.
获得基态NO分子精确的完全振动能级后, 根据量子统计系综理论计算了1000—5000 K温度范围内NO分子气体系统摩尔内能、摩尔熵和摩尔热容. 在高温低压气体条件下, 系统处于理想气体状态, 去掉分子热运动的热容量, 剩余分子振动和转动对应热容量
$C_{{\rm{exp}}}^{} = C_{\rm{p}}^{{\rm{exp}}} - R - 3 R/2 = C_{\rm{p}}^{{\rm{exp}}} -$ 5R/2, 其中普适气体常数R = 8.31441 J·K–1·mol–1. 采用Chase[26]得到的NO分子气体定压摩尔热容量的实验数据与理论计算结果进行比较. 表2和表3列出了量子统计系综理论在考虑不同振动配分函数贡献($Q_{\rm{\upsilon }}^{{\rm{RKR}}}(T)$ ,$Q_{\rm{\upsilon }}^{{\rm{MCSCF}}}(T)$ ,$Q_{\rm{\upsilon }}^{{\rm{SHO}}}(T)$ ,$Q_{\rm{\upsilon }}^{{\rm{VAM}}}(T)$ )条件下计算得到的内能和熵.T/K $U_{{\rm{RKR}}}^{{\rm{cal}}}$ $U_{{\rm{MCSCF}}}^{{\rm{cal}}}$ $U_{{\rm{SHO}}}^{{\rm{cal}}}$ $U_{{\rm{VAM}}}^{{\rm{cal}}}$ 1000 21.28 21.28 21.27 21.28 1100 22.61 22.61 22.59 22.61 1200 23.98 23.99 23.95 23.98 1300 25.39 25.40 25.34 25.39 1400 26.84 26.84 26.77 26.84 1500 28.31 28.31 28.23 28.31 1600 29.80 29.80 29.70 29.80 1700 31.31 31.31 31.20 31.31 1800 32.84 32.84 32.71 32.84 1900 34.38 34.39 34.24 34.38 2000 35.94 35.94 35.77 35.94 2100 37.51 37.51 37.32 37.51 2200 39.09 39.09 38.88 39.09 2300 40.68 40.68 40.45 40.68 2400 42.28 42.28 42.03 42.28 2500 43.89 43.89 43.61 43.89 2600 45.50 45.50 45.19 45.50 2700 47.12 47.12 46.79 47.12 2800 48.75 48.75 48.38 48.75 2900 50.38 50.38 49.99 50.38 3000 52.01 52.01 51.59 52.01 3100 53.66 53.65 53.20 53.66 3200 55.30 55.30 54.81 55.30 3300 56.95 56.95 56.43 56.95 3400 58.60 58.60 58.05 58.60 3500 60.26 60.26 59.67 60.26 3600 61.92 61.92 61.29 61.92 3700 63.58 63.58 62.91 63.59 3800 65.25 65.24 64.54 65.25 3900 66.91 66.91 66.17 66.92 4000 68.58 68.58 67.80 68.60 4100 70.26 70.25 69.43 70.27 4200 71.93 71.93 71.07 71.95 4300 73.61 73.61 72.70 73.63 4400 75.29 75.28 74.34 75.32 4500 76.97 76.96 75.97 77.01 4600 78.65 78.64 77.61 78.69 4700 80.33 80.32 79.25 80.39 4800 82.01 82.01 80.89 82.08 4900 83.69 83.69 82.53 83.78 5000 85.38 85.37 84.17 85.48 T/K $S_{{\rm{RKR}}}^{{\rm{cal}}}$ $S_{{\rm{MCSCF}}}^{{\rm{cal}}}$ $S_{{\rm{SHO}}}^{{\rm{cal}}}$ $S_{{\rm{VAM}}}^{{\rm{cal}}}$ 1000 60.50 60.50 60.42 60.50 1100 61.76 61.76 61.67 61.76 1200 62.96 62.96 62.86 62.96 1300 64.09 64.09 63.97 64.09 1400 65.16 65.16 65.03 65.16 1500 66.17 66.17 66.03 66.17 1600 67.13 67.13 66.99 67.13 1700 68.05 68.05 67.89 68.05 1800 68.92 68.92 68.76 68.92 1900 69.76 69.76 69.58 69.76 2000 70.56 70.56 70.37 70.56 2100 71.32 71.32 71.13 71.32 2200 72.06 72.06 71.85 72.06 2300 72.77 72.77 72.55 72.77 2400 73.45 73.45 73.22 73.45 2500 74.10 74.10 73.87 74.10 2600 74.73 74.73 74.49 74.73 2700 75.35 75.35 75.09 75.35 2800 75.94 75.94 75.67 75.94 2900 76.51 76.51 76.23 76.51 3000 77.06 77.06 76.78 77.06 3100 77.60 77.60 77.30 77.60 3200 78.12 78.12 77.82 78.12 3300 78.63 78.63 78.31 78.63 3400 79.12 79.12 78.80 79.13 3500 79.60 79.60 79.27 79.61 3600 80.07 80.07 79.72 80.07 3700 80.53 80.53 80.17 80.53 3800 80.97 80.97 80.60 80.97 3900 81.41 81.40 81.03 81.41 4000 81.83 81.83 81.44 81.83 4100 82.24 82.24 81.84 82.25 4200 82.65 82.64 82.24 82.65 4300 83.04 83.04 82.62 83.05 4400 83.43 83.42 83.00 83.43 4500 83.80 83.80 83.36 83.81 4600 84.17 84.17 83.72 84.18 4700 84.53 84.53 84.08 84.55 4800 84.89 84.89 84.42 84.90 4900 85.24 85.23 84.76 85.25 5000 85.58 85.57 85.09 85.60 内能、熵、热容均是重要的热力学参数, 但在实验中常常测定的是热容的值
${C_{{\rm{exp}}}}$ . 在表4中将得到的热容$C_{{\rm{VAM}}}^{{\rm{cal}}}$ 与实验测量的热容${C_{{\rm{exp}}}}$ 进行比较, 在1200 K时与实验值的相对误差仅为1.16%. 在整个温度区间范围内, 虽然相对误差随着温度的升高而增大, 但是最大值仅为2.36%, 小于采用其他方法计算的误差.T/K ${C_{{\rm{exp}}}}$ $C_{{\rm{RKR}}}^{{\rm{cal}}}$ $C_{{\rm{MCSCF}}}^{{\rm{cal}}}$ $C_{{\rm{SHO}}}^{{\rm{cal}}}$ $C_{{\rm{VAM}}}^{{\rm{cal}}}$ ${\delta _{{\rm{RKR}}}}$a ${\delta _{{\rm{MCSCF}}}}$b ${\delta _{{\rm{SHO}}}}$c ${\delta _{{\rm{VAM}}}}$d 1000 13.20 13.04 13.04 10.72 13.04 1.22% 1.22% 18.82% 1.22% 1100 13.68 13.52 13.52 11.24 13.52 1.18% 1.18% 17.85% 1.18% 1200 14.09 13.93 13.93 11.73 13.93 1.16% 1.16% 16.77% 1.16% 1300 14.44 14.27 14.27 12.17 14.27 1.17% 1.17% 15.69% 1.17% 1400 14.74 14.56 14.56 12.58 14.56 1.19% 1.19% 14.66% 1.19% 1500 14.99 14.81 14.81 12.94 14.81 1.21% 1.21% 13.70% 1.21% 1600 15.22 15.03 15.03 13.26 15.03 1.23% 1.23% 12.83% 1.23% 1700 15.41 15.22 15.21 13.55 15.22 1.26% 1.26% 12.04% 1.26% 1800 15.58 15.38 15.38 13.81 15.38 1.29% 1.29% 11.35% 1.29% 1900 15.73 15.52 15.52 14.04 15.52 1.33% 1.33% 10.73% 1.33% 2000 15.86 15.65 15.64 14.25 15.65 1.36% 1.37% 10.18% 1.36% 2100 15.98 15.76 15.76 14.43 15.76 1.41% 1.41% 9.71% 1.41% 2200 16.09 15.86 15.86 14.59 15.86 1.45% 1.45% 9.29% 1.45% 2300 16.18 15.95 15.94 14.74 15.95 1.48% 1.49% 8.92% 1.48% 2400 16.27 16.03 16.03 14.88 16.03 1.52% 1.53% 8.59% 1.52% 2500 16.35 16.10 16.10 15.00 16.10 1.56% 1.57% 8.31% 1.56% 2600 16.43 16.17 16.17 15.10 16.17 1.61% 1.61% 8.07% 1.61% 2700 16.50 16.23 16.23 15.20 16.23 1.65% 1.65% 7.85% 1.65% 2800 16.56 16.28 16.28 15.29 16.28 1.69% 1.70% 7.67% 1.69% 2900 16.62 16.34 16.33 15.38 16.34 1.74% 1.74% 7.51% 1.73% 3000 16.68 16.38 16.38 15.45 16.39 1.77% 1.78% 7.37% 1.77% 3100 16.73 16.43 16.43 15.52 16.43 1.82% 1.83% 7.25% 1.80% 3200 16.78 16.47 16.47 15.58 16.47 1.87% 1.88% 7.16% 1.85% 3300 16.83 16.51 16.51 15.64 16.51 1.91% 1.92% 7.08% 1.88% 3400 16.88 16.55 16.54 15.69 16.55 1.96% 1.97% 7.01% 1.92% 3500 16.92 16.58 16.58 15.74 16.59 2.01% 2.02% 6.96% 1.95% 3600 16.96 16.61 16.61 15.79 16.62 2.06% 2.07% 6.91% 1.99% 3700 17.00 16.64 16.64 15.83 16.66 2.12% 2.13% 6.89% 2.03% 3800 17.04 16.67 16.67 15.87 16.69 2.18% 2.19% 6.86% 2.06% 3900 17.08 16.69 16.69 15.91 16.72 2.24% 2.25% 6.85% 2.09% 4000 17.11 16.72 16.71 15.94 16.75 2.31% 2.32% 6.85% 2.13% 4100 17.15 16.74 16.74 15.97 16.78 2.38% 2.40% 6.86% 2.16% 4200 17.18 16.76 16.75 16.00 16.81 2.46% 2.48% 6.86% 2.18% 4300 17.21 16.77 16.77 16.03 16.83 2.55% 2.57% 6.88% 2.21% 4400 17.24 16.79 16.79 16.05 16.86 2.64% 2.66% 6.90% 2.24% 4500 17.28 16.80 16.80 16.08 16.89 2.74% 2.76% 6.93% 2.26% 4600 17.31 16.81 16.81 16.10 16.91 2.85% 2.86% 6.96% 2.29% 4700 17.34 16.82 16.82 16.12 16.94 2.96% 2.99% 7.00% 2.31% 4800 17.36 16.83 16.83 16.14 16.96 3.09% 3.11% 7.04% 2.32% 4900 17.39 16.83 16.83 16.16 16.99 3.22% 3.25% 7.08% 2.35% 5000 17.42 16.84 16.83 16.18 17.01 3.37% 3.38% 7.13% 2.36% 注: a, ${\delta _{{\rm{RKR}}}} = \left| {C_{{\rm{exp}}}^{} - C_{{\rm{RKR}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}}^{} \times 100\% $; b, ${\delta _{{\rm{MCSCF}}}} = \left| {C_{{\rm{exp}}}^{} - C_{{\rm{MCSCF}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}}^{} \times 100\% $; c, $ {\delta _{{\rm{SHO}}}} = \left| {C_{{\rm{exp}}} - C_{{\rm{SHO}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}} \times 100\% $; d,$ {\delta _{{\rm{VAM}}}} = \left| {C_{{\rm{exp}}} - C_{{\rm{VAM}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}} \times 100\% .$ 表4还列出了
$C_{{\rm{RKR}}}^{{\rm{cal}}}$ ,$C_{{\rm{MCSCF}}}^{{\rm{cal}}}$ 和$C_{{\rm{SHO}}}^{{\rm{cal}}}$ 以及这些结果与实验测量的热容$C_{{\rm{exp}}}$ 的相对误差. 由表4知, 仅将实验获得的8个振动能级进行统计求和计算热容$C_{{\rm{RKR}}}^{{\rm{cal}}}$ 时, 由于大量缺乏高阶振动能级对系统的统计贡献, 热容结果总是小于实验测得的结果, 且与实验值的误差随着温度的升高而增大, 在5000 K时与实验值的相对误差达到了3.37%. 将MCSCF振动能级用于统计公式计算热容$C_{{\rm{MCSCF}}}^{{\rm{cal}}}$ 时, 其参与计算的最高振动能级能量为29023.2 cm–1, 仅达到离解能的55.65%, 仍然缺乏高阶振动能级对热容的贡献. 因此其结果与实验值的误差还是很大, 相对误差达到了3.38%. 而简单谐振子模型对实际上不存在的虚高能级进行统计求和, 偏离了物理事实, 计算得到的热容误差与实验值相比误差最大. 其相对误差最小时也有6.85%, 远大于使用VAM完全振动能级计算的热容与实验测得的热容值的误差. 说明简单谐振子模型对无限能级求和方法对求解宏观气体系统的热力学性质有一定的局限性. 本文未考虑激发态对配分函数的影响, 所以计算得到的热容值均小于实验值. 综上所述, 使用VAM完全振动能级作为振动贡献计算获得的热容量$C_{{\rm{VAM}}}^{{\rm{cal}}}$ 与实验测得热容$C_{{\rm{exp}}}$ 相比误差最小.图2可以更直观地体现不同方法获得的热容计算值与热容实验值的相对误差. 由图2可知VAM热容
$C_{{\rm{VAM}}}^{{\rm{cal}}}$ 与实验上测得的热容${C_{{\rm{exp}}}}$ 误差最小, 谐振子热容$C_{{\rm{SHO}}}^{{\rm{cal}}}$ 与实验上测得的热容${C_{{\rm{exp}}}}$ 误差最大. 在1000—2800 K温度范围内,$C_{{\rm{RKR}}}^{{\rm{cal}}}$ ,$C_{\rm VAM}^{\rm cal}$ 和$C_{{\rm{MCSCF}}}^{{\rm{cal}}}$ 三者的值较为吻合, 在2900—5000 K较高温度范围内, 较$C_{{\rm{VAM}}}^{{\rm{cal}}}$ 而言,$C_{{\rm{RKR}}}^{{\rm{cal}}}$ 和$C_{{\rm{MCSCF}}}^{{\rm{cal}}}$ 热容值明显偏离了实验热容${C_{{\rm{exp}}}}$ . 这一结果表明, 计算热容时振动部分采用VAM方法获得完全振动能级的计算结果精度优于其他方法获得的计算结果.4. 结 论
使用前期建立的VAM方法, 基于实验上获得的基态NO分子8个振动能级得到了正确收敛于实验离解能的45个完全振动能级. 发现振动量子数为2的振动能级与实验值有较大的差异, 但其他理论研究者的结果也在本文的研究结果附近. 分别将VAM完全振动能级、实验上得到的有限的振动能级、MCSCF方法得到的振动能级和简谐振子模型得到的振动能级作为振动贡献, 将欧拉-麦克劳林渐进展开公式的计算结果作为转动贡献, 用量子系综理论计算了NO气体系统的摩尔内能、摩尔熵和摩尔热容. 将不同方法计算获得的摩尔热容与实验值进行比较发现: 缺乏高阶振动能级会给求解系统的热力学性质带来较大误差; 对实际上不存在的虚高能级进行统计求和也会给计算结果带来很大误差; 振动部分采用VAM完全振动能级所计算的热容与实验上测得的热容最为吻合. 这说明使用谐振子模型和不完全振动能级对计算热力学性质均有一定的局限性, 应当使用有限的完全振动能级进行统计求和.
[1] Brown W A, King D A 2000 J. Phys. Chem. B 104 2578
[2] Palmer R M J, Ferrige A G, Moncada S 1987 Nature 327 524
[3] Wang Z Z, Zhao Y Y, Sun R, et al. 2019 Fuel 253 1424Google Scholar
[4] Goldstein I, Lue T F, Padma-Nathan H, et al. 1998 N. Engl. J. Med. 338 1397Google Scholar
[5] Meng Q T, Liu X G, Zhang Q G, Han K L 2005 Chem. Phys. 316 93Google Scholar
[6] 王德华, 王雅静, 薛艳丽, 李洪云, 林圣路 2007 物理学报 56 6209Google Scholar
Wang D L, Wang Y J, XueY L, Li H Y, Lin S L 2007 Acta Phys. Sin. 56 6209Google Scholar
[7] Shrestha K P, Seidel L, Zeuch T, Mauss F 2019 Combust. Sci. Technol. 191 1628Google Scholar
[8] Chen L, Wang D, Wang J D, Weng D, Cao L 2019 J. Rare Earths 37 829Google Scholar
[9] Matzkin A, Raoult M, Gauyacq D 2003 Phys. Rev. A 68 061401Google Scholar
[10] Jaffe R L 1987 AIAA 22nd Thermophysics Conference New York, USA, June 8–10, 1999 p1633
[11] Capitelli M, Colonna G, Giordano D, et al. 2005 Tables of Internal Partition Functions and Thermodynamic Properties of High-temperature Mars-atmosphere Species from 50 K to 50000 K (Netherlands: European Space Agency Publications Division) pp3–19
[12] Babou Y, Rivière P, Perrin M Y, Soufiani A 2009 Int. J. Thermophys. 30 416Google Scholar
[13] 邓伦华, 李传亮, 朱圆月, 何文艳, 陈扬骎 2012 物理学报 61 194208Google Scholar
Deng L H, Li C L, Zhu Y Y, He W Y, Chen Y Q 2012 Acta Phys. Sin. 61 194208Google Scholar
[14] Pathria R K1977 Statistical Mechanics (London: Pergamon press) pp100–107
[15] Huber K P, Herzberg G 1950 Molecular Spectra and Molecular Structure: Spectra of Diatomic Molecules (New York: Van Nostrand Reinhold Company) pp9–11
[16] Billingsley F P 1975 J. Chem. Phys. 62 864Google Scholar
[17] Reddy R R, Ahammed Y N, Basha D B, et al. 2006 J. Quant. Spectrosc. Radiat. Transfer 97 344Google Scholar
[18] Qin Z, Zhao J M, Liu L H 2018 J. Quant. Spectrosc.Radiat. Transfer 210 1Google Scholar
[19] Zhang Y, Sun W G, Fu J, Fan Q C, et al. 2014 Spectrochim. Acta Part A 117 442Google Scholar
[20] 周金伟, 李吉成, 石志广, 陈小天, 卢晓卫 2014 光学学报 34 130001Google Scholar
Zhou J W, Li J C, Shi Z G, Chen X T, Lu X W 2014 Acta Optic. Sin. 34 130001Google Scholar
[21] Gamachea R R, Kennedya S, Hawkinsb R, Rothmanb L S 2000 J. Mol. Struct. 517 407
[22] McDowell R S 1988 J. Chem. Phys. 88 356Google Scholar
[23] Amiot C 1982 J. Mol. Spectrosc. 94 150Google Scholar
[24] Sun W G, Hou S L, Feng H, Ren W Y 2002 J. Mol. Spectrosc. 215 93Google Scholar
[25] Wang J K, Yang Z Y, Wu Z S 2010 Acta Photonica Sin. 39 1312Google Scholar
[26] Chase M W 1998 Journal of Physical and Chemical Reference DataMonograph No.9 (New York: National Institute of Standards and Technology Gaithersburg) pp641–643
-
表 1 基态NO分子不同振动能级间的比较 (单位: cm–1)
Table 1. Comparison of different vibrational levels of NO in the ground state (in cm–1).
$\upsilon $ $E_{\rm{\upsilon }}^{{\rm{exp}}}$[17] $E_{\rm{\upsilon }}^{{\rm{MCSCF}}}$[16] $E_{\rm{\upsilon }}^{{\rm{VAM}}}$ $E_{\rm{\upsilon }}^{{\rm{exp}}} - E_{\rm{\upsilon }}^{{\rm{MCSCF}}}$ $E_{\rm{\upsilon }}^{{\rm{exp}}} - E_{\rm{\upsilon }}^{{\rm{VAM}}}$ υ $E_{\rm{\upsilon }}^{{\rm{VAM}}}$ 0 948.50 948.60 948.50 –0.10 0 26 40240.70 1 2824.50 2824.60 2824.50 –0.10 0 27 41310.84 2 4627.30 4672.40 4672.27 –45.10 –44.97 28 42340.94 3 6491.90 6492.10 6491.90 –0.20 0 29 43329.70 4 8283.50 8283.90 8283.43 –0.40 0.07 30 44275.74 5 10046.90 10047.80 10046.90 –0.90 0 31 45177.59 6 11782.30 11783.70 11782.30 –1.40 0 32 46033.69 7 13489.40 13491.70 13489.60 –2.30 –0.20 33 46842.41 8 15171.80 15168.73 34 47602.00 9 16823.90 16819.60 35 48310.62 10 18448.00 18442.06 36 48966.35 11 20044.00 20035.93 37 49567.15 12 21611.80 21600.99 38 50110.90 13 23151.30 23136.97 39 50595.35 14 24662.40 24643.56 40 51018.15 15 26144.80 26120.38 41 51376.86 16 27598.50 27567.04 42 51668.92 17 29023.20 28983.04 43 51891.64 18 30367.88 44 52042.23 19 31720.95 20 33041.62 21 34329.18 22 35582.84 23 36801.78 24 37985.07 25 39131.73 $D_{\rm{e}}^{{\rm{exp}}}$ 52155.68 $D_{\rm{e}}^{{\rm{cal}}}$ 52155.68 注: $E_{\rm{\upsilon }}^{{\rm{VAM}}}$中表示VAM计算所需的已知实验振动能级用黑体标出. 表 2 使用RKR, MCSCF, SHO和VAM振动能级计算的摩尔内能 (单位: J·K–1·mol–1)
Table 2. Calculated molar internal energy using RKR, MCSCF, SHO, and VAM vibrational energies as the vibrational contributions, respectively (in J·K–1·mol–1).
T/K $U_{{\rm{RKR}}}^{{\rm{cal}}}$ $U_{{\rm{MCSCF}}}^{{\rm{cal}}}$ $U_{{\rm{SHO}}}^{{\rm{cal}}}$ $U_{{\rm{VAM}}}^{{\rm{cal}}}$ 1000 21.28 21.28 21.27 21.28 1100 22.61 22.61 22.59 22.61 1200 23.98 23.99 23.95 23.98 1300 25.39 25.40 25.34 25.39 1400 26.84 26.84 26.77 26.84 1500 28.31 28.31 28.23 28.31 1600 29.80 29.80 29.70 29.80 1700 31.31 31.31 31.20 31.31 1800 32.84 32.84 32.71 32.84 1900 34.38 34.39 34.24 34.38 2000 35.94 35.94 35.77 35.94 2100 37.51 37.51 37.32 37.51 2200 39.09 39.09 38.88 39.09 2300 40.68 40.68 40.45 40.68 2400 42.28 42.28 42.03 42.28 2500 43.89 43.89 43.61 43.89 2600 45.50 45.50 45.19 45.50 2700 47.12 47.12 46.79 47.12 2800 48.75 48.75 48.38 48.75 2900 50.38 50.38 49.99 50.38 3000 52.01 52.01 51.59 52.01 3100 53.66 53.65 53.20 53.66 3200 55.30 55.30 54.81 55.30 3300 56.95 56.95 56.43 56.95 3400 58.60 58.60 58.05 58.60 3500 60.26 60.26 59.67 60.26 3600 61.92 61.92 61.29 61.92 3700 63.58 63.58 62.91 63.59 3800 65.25 65.24 64.54 65.25 3900 66.91 66.91 66.17 66.92 4000 68.58 68.58 67.80 68.60 4100 70.26 70.25 69.43 70.27 4200 71.93 71.93 71.07 71.95 4300 73.61 73.61 72.70 73.63 4400 75.29 75.28 74.34 75.32 4500 76.97 76.96 75.97 77.01 4600 78.65 78.64 77.61 78.69 4700 80.33 80.32 79.25 80.39 4800 82.01 82.01 80.89 82.08 4900 83.69 83.69 82.53 83.78 5000 85.38 85.37 84.17 85.48 表 3 使用RKR, MCSCF、SHO和VAM振动能级计算的摩尔熵(单位: J·K–1·mol–1)
Table 3. Calculated molar entropy using RKR, MCSCF, SHO, and VAM vibrational energies as the vibrational contributions, respectively (in J·K–1·mol–1).
T/K $S_{{\rm{RKR}}}^{{\rm{cal}}}$ $S_{{\rm{MCSCF}}}^{{\rm{cal}}}$ $S_{{\rm{SHO}}}^{{\rm{cal}}}$ $S_{{\rm{VAM}}}^{{\rm{cal}}}$ 1000 60.50 60.50 60.42 60.50 1100 61.76 61.76 61.67 61.76 1200 62.96 62.96 62.86 62.96 1300 64.09 64.09 63.97 64.09 1400 65.16 65.16 65.03 65.16 1500 66.17 66.17 66.03 66.17 1600 67.13 67.13 66.99 67.13 1700 68.05 68.05 67.89 68.05 1800 68.92 68.92 68.76 68.92 1900 69.76 69.76 69.58 69.76 2000 70.56 70.56 70.37 70.56 2100 71.32 71.32 71.13 71.32 2200 72.06 72.06 71.85 72.06 2300 72.77 72.77 72.55 72.77 2400 73.45 73.45 73.22 73.45 2500 74.10 74.10 73.87 74.10 2600 74.73 74.73 74.49 74.73 2700 75.35 75.35 75.09 75.35 2800 75.94 75.94 75.67 75.94 2900 76.51 76.51 76.23 76.51 3000 77.06 77.06 76.78 77.06 3100 77.60 77.60 77.30 77.60 3200 78.12 78.12 77.82 78.12 3300 78.63 78.63 78.31 78.63 3400 79.12 79.12 78.80 79.13 3500 79.60 79.60 79.27 79.61 3600 80.07 80.07 79.72 80.07 3700 80.53 80.53 80.17 80.53 3800 80.97 80.97 80.60 80.97 3900 81.41 81.40 81.03 81.41 4000 81.83 81.83 81.44 81.83 4100 82.24 82.24 81.84 82.25 4200 82.65 82.64 82.24 82.65 4300 83.04 83.04 82.62 83.05 4400 83.43 83.42 83.00 83.43 4500 83.80 83.80 83.36 83.81 4600 84.17 84.17 83.72 84.18 4700 84.53 84.53 84.08 84.55 4800 84.89 84.89 84.42 84.90 4900 85.24 85.23 84.76 85.25 5000 85.58 85.57 85.09 85.60 表 4 不同摩尔热容及其与实验值的相对误差 (单位: J·K–1·mol–1)
Table 4. Comparisons of different molar capacities with observed experimental
${C_{{\rm{exp}}}}$ (in J·K–1·mol–1).T/K ${C_{{\rm{exp}}}}$ $C_{{\rm{RKR}}}^{{\rm{cal}}}$ $C_{{\rm{MCSCF}}}^{{\rm{cal}}}$ $C_{{\rm{SHO}}}^{{\rm{cal}}}$ $C_{{\rm{VAM}}}^{{\rm{cal}}}$ ${\delta _{{\rm{RKR}}}}$a ${\delta _{{\rm{MCSCF}}}}$b ${\delta _{{\rm{SHO}}}}$c ${\delta _{{\rm{VAM}}}}$d 1000 13.20 13.04 13.04 10.72 13.04 1.22% 1.22% 18.82% 1.22% 1100 13.68 13.52 13.52 11.24 13.52 1.18% 1.18% 17.85% 1.18% 1200 14.09 13.93 13.93 11.73 13.93 1.16% 1.16% 16.77% 1.16% 1300 14.44 14.27 14.27 12.17 14.27 1.17% 1.17% 15.69% 1.17% 1400 14.74 14.56 14.56 12.58 14.56 1.19% 1.19% 14.66% 1.19% 1500 14.99 14.81 14.81 12.94 14.81 1.21% 1.21% 13.70% 1.21% 1600 15.22 15.03 15.03 13.26 15.03 1.23% 1.23% 12.83% 1.23% 1700 15.41 15.22 15.21 13.55 15.22 1.26% 1.26% 12.04% 1.26% 1800 15.58 15.38 15.38 13.81 15.38 1.29% 1.29% 11.35% 1.29% 1900 15.73 15.52 15.52 14.04 15.52 1.33% 1.33% 10.73% 1.33% 2000 15.86 15.65 15.64 14.25 15.65 1.36% 1.37% 10.18% 1.36% 2100 15.98 15.76 15.76 14.43 15.76 1.41% 1.41% 9.71% 1.41% 2200 16.09 15.86 15.86 14.59 15.86 1.45% 1.45% 9.29% 1.45% 2300 16.18 15.95 15.94 14.74 15.95 1.48% 1.49% 8.92% 1.48% 2400 16.27 16.03 16.03 14.88 16.03 1.52% 1.53% 8.59% 1.52% 2500 16.35 16.10 16.10 15.00 16.10 1.56% 1.57% 8.31% 1.56% 2600 16.43 16.17 16.17 15.10 16.17 1.61% 1.61% 8.07% 1.61% 2700 16.50 16.23 16.23 15.20 16.23 1.65% 1.65% 7.85% 1.65% 2800 16.56 16.28 16.28 15.29 16.28 1.69% 1.70% 7.67% 1.69% 2900 16.62 16.34 16.33 15.38 16.34 1.74% 1.74% 7.51% 1.73% 3000 16.68 16.38 16.38 15.45 16.39 1.77% 1.78% 7.37% 1.77% 3100 16.73 16.43 16.43 15.52 16.43 1.82% 1.83% 7.25% 1.80% 3200 16.78 16.47 16.47 15.58 16.47 1.87% 1.88% 7.16% 1.85% 3300 16.83 16.51 16.51 15.64 16.51 1.91% 1.92% 7.08% 1.88% 3400 16.88 16.55 16.54 15.69 16.55 1.96% 1.97% 7.01% 1.92% 3500 16.92 16.58 16.58 15.74 16.59 2.01% 2.02% 6.96% 1.95% 3600 16.96 16.61 16.61 15.79 16.62 2.06% 2.07% 6.91% 1.99% 3700 17.00 16.64 16.64 15.83 16.66 2.12% 2.13% 6.89% 2.03% 3800 17.04 16.67 16.67 15.87 16.69 2.18% 2.19% 6.86% 2.06% 3900 17.08 16.69 16.69 15.91 16.72 2.24% 2.25% 6.85% 2.09% 4000 17.11 16.72 16.71 15.94 16.75 2.31% 2.32% 6.85% 2.13% 4100 17.15 16.74 16.74 15.97 16.78 2.38% 2.40% 6.86% 2.16% 4200 17.18 16.76 16.75 16.00 16.81 2.46% 2.48% 6.86% 2.18% 4300 17.21 16.77 16.77 16.03 16.83 2.55% 2.57% 6.88% 2.21% 4400 17.24 16.79 16.79 16.05 16.86 2.64% 2.66% 6.90% 2.24% 4500 17.28 16.80 16.80 16.08 16.89 2.74% 2.76% 6.93% 2.26% 4600 17.31 16.81 16.81 16.10 16.91 2.85% 2.86% 6.96% 2.29% 4700 17.34 16.82 16.82 16.12 16.94 2.96% 2.99% 7.00% 2.31% 4800 17.36 16.83 16.83 16.14 16.96 3.09% 3.11% 7.04% 2.32% 4900 17.39 16.83 16.83 16.16 16.99 3.22% 3.25% 7.08% 2.35% 5000 17.42 16.84 16.83 16.18 17.01 3.37% 3.38% 7.13% 2.36% 注: a, ${\delta _{{\rm{RKR}}}} = \left| {C_{{\rm{exp}}}^{} - C_{{\rm{RKR}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}}^{} \times 100\% $; b, ${\delta _{{\rm{MCSCF}}}} = \left| {C_{{\rm{exp}}}^{} - C_{{\rm{MCSCF}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}}^{} \times 100\% $; c, $ {\delta _{{\rm{SHO}}}} = \left| {C_{{\rm{exp}}} - C_{{\rm{SHO}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}} \times 100\% $; d,$ {\delta _{{\rm{VAM}}}} = \left| {C_{{\rm{exp}}} - C_{{\rm{VAM}}}^{{\rm{cal}}}} \right|/C_{{\rm{exp}}} \times 100\% .$ -
[1] Brown W A, King D A 2000 J. Phys. Chem. B 104 2578
[2] Palmer R M J, Ferrige A G, Moncada S 1987 Nature 327 524
[3] Wang Z Z, Zhao Y Y, Sun R, et al. 2019 Fuel 253 1424Google Scholar
[4] Goldstein I, Lue T F, Padma-Nathan H, et al. 1998 N. Engl. J. Med. 338 1397Google Scholar
[5] Meng Q T, Liu X G, Zhang Q G, Han K L 2005 Chem. Phys. 316 93Google Scholar
[6] 王德华, 王雅静, 薛艳丽, 李洪云, 林圣路 2007 物理学报 56 6209Google Scholar
Wang D L, Wang Y J, XueY L, Li H Y, Lin S L 2007 Acta Phys. Sin. 56 6209Google Scholar
[7] Shrestha K P, Seidel L, Zeuch T, Mauss F 2019 Combust. Sci. Technol. 191 1628Google Scholar
[8] Chen L, Wang D, Wang J D, Weng D, Cao L 2019 J. Rare Earths 37 829Google Scholar
[9] Matzkin A, Raoult M, Gauyacq D 2003 Phys. Rev. A 68 061401Google Scholar
[10] Jaffe R L 1987 AIAA 22nd Thermophysics Conference New York, USA, June 8–10, 1999 p1633
[11] Capitelli M, Colonna G, Giordano D, et al. 2005 Tables of Internal Partition Functions and Thermodynamic Properties of High-temperature Mars-atmosphere Species from 50 K to 50000 K (Netherlands: European Space Agency Publications Division) pp3–19
[12] Babou Y, Rivière P, Perrin M Y, Soufiani A 2009 Int. J. Thermophys. 30 416Google Scholar
[13] 邓伦华, 李传亮, 朱圆月, 何文艳, 陈扬骎 2012 物理学报 61 194208Google Scholar
Deng L H, Li C L, Zhu Y Y, He W Y, Chen Y Q 2012 Acta Phys. Sin. 61 194208Google Scholar
[14] Pathria R K1977 Statistical Mechanics (London: Pergamon press) pp100–107
[15] Huber K P, Herzberg G 1950 Molecular Spectra and Molecular Structure: Spectra of Diatomic Molecules (New York: Van Nostrand Reinhold Company) pp9–11
[16] Billingsley F P 1975 J. Chem. Phys. 62 864Google Scholar
[17] Reddy R R, Ahammed Y N, Basha D B, et al. 2006 J. Quant. Spectrosc. Radiat. Transfer 97 344Google Scholar
[18] Qin Z, Zhao J M, Liu L H 2018 J. Quant. Spectrosc.Radiat. Transfer 210 1Google Scholar
[19] Zhang Y, Sun W G, Fu J, Fan Q C, et al. 2014 Spectrochim. Acta Part A 117 442Google Scholar
[20] 周金伟, 李吉成, 石志广, 陈小天, 卢晓卫 2014 光学学报 34 130001Google Scholar
Zhou J W, Li J C, Shi Z G, Chen X T, Lu X W 2014 Acta Optic. Sin. 34 130001Google Scholar
[21] Gamachea R R, Kennedya S, Hawkinsb R, Rothmanb L S 2000 J. Mol. Struct. 517 407
[22] McDowell R S 1988 J. Chem. Phys. 88 356Google Scholar
[23] Amiot C 1982 J. Mol. Spectrosc. 94 150Google Scholar
[24] Sun W G, Hou S L, Feng H, Ren W Y 2002 J. Mol. Spectrosc. 215 93Google Scholar
[25] Wang J K, Yang Z Y, Wu Z S 2010 Acta Photonica Sin. 39 1312Google Scholar
[26] Chase M W 1998 Journal of Physical and Chemical Reference DataMonograph No.9 (New York: National Institute of Standards and Technology Gaithersburg) pp641–643
期刊类型引用(1)
1. 文琳,樊群超,蹇君,范志祥,李会东,付佳,马杰,谢锋. 基于SO分子振转能级计算其宏观气体摩尔热容. 物理学报. 2022(17): 264-271 . 百度学术 其他类型引用(3)
Catalog
Metrics
- Abstract views: 7818
- PDF Downloads: 63
- Cited By: 4