CO+分子离子X2Σ+, A2Π和B2Σ+态的不透明度

上一篇

下一篇

安斯腰力吐, 王佟, 肖利丹, 刘迪, 张夏, 闫冰. CO+分子离子X2Σ+, A2Π和B2Σ+态的不透明度[J]. 物理学报, 2025, 74(12): 123101-1. doi: 10.7498/aps.74.20250380
引用本文: 安斯腰力吐, 王佟, 肖利丹, 刘迪, 张夏, 闫冰. CO+分子离子X2Σ+, A2Π和B2Σ+态的不透明度[J]. 物理学报, 2025, 74(12): 123101-1. doi: 10.7498/aps.74.20250380
Siyaolitu AN, Tong WANG, Lidan XIAO, Di LIU, Xia ZHANG, Bing YAN. Opacities of X2Σ+, A2Π, and B2Σ+ states of CO+ molecule ion[J]. Acta Physica Sinica, 2025, 74(12): 123101-1. doi: 10.7498/aps.74.20250380
Citation: Siyaolitu AN, Tong WANG, Lidan XIAO, Di LIU, Xia ZHANG, Bing YAN. Opacities of X2Σ+, A2Π, and B2Σ+ states of CO+ molecule ion[J]. Acta Physica Sinica, 2025, 74(12): 123101-1. doi: 10.7498/aps.74.20250380

CO+分子离子X2Σ+, A2Π和B2Σ+态的不透明度

    作者简介: 安斯腰力吐:3564414330@qq.com .
    通讯作者: E-mail: xiazhang@jlu.edu.cn.;  E-mail: yanbing@jlu.edu.cn.
  • 中图分类号: 31.50.Df, 31.15.ag, 31.15.aj

Opacities of X2Σ+, A2Π, and B2Σ+ states of CO+ molecule ion

    Corresponding authors: E-mail: xiazhang@jlu.edu.cn.;  E-mail: yanbing@jlu.edu.cn.
  • MSC: 31.50.Df, 31.15.ag, 31.15.aj

  • 摘要: 一氧化碳分子离子(CO+)在大气及天体物理环境中起着关键作用, 其不透明度的理论研究对辐射输运建模具有重要意义. 本文基于实验观测数据, 采用Modified-Morse (MMorse)势函数改进并构建了CO+分子离子X2Σ+, A2Π和B2Σ+电子态的势能曲线, 进一步提取了振动能级和光谱常数. 同时, 利用考虑Davidson修正的多参考组态相互作用(MRCI+Q)方法计算了势能曲线和电偶极跃迁矩. 改进获得的MMorse势与计算得到的势能曲线非常吻合, 且光谱常数和振动能级与其他理论和实验数据符合较好. 结合MMorse势函数和从头计算获得的电偶极跃迁矩, 计算了CO+分子离子在100 atm (1 atm=1.01×105 Pa)压强下, 298—15000 K温度范围内的不透明度, 并探究了不同温度对高温谱的影响. 研究结果表明, 高温环境(T > 5000 K)会导致不同能带系统的谱线展宽与边界模糊, 这种混合效应在T > 10000 K时尤为突出, 揭示了高温下分子离子光谱退化的微观机制. 本研究可以为天体物理领域提供一些理论依据和数据支持. 本文数据集可在https://doi.org/10.57760/sciencedb.j00213.00136中访问获取.
  • 加载中
  • 图 1  CO+的X2Σ+, A2Π和B2Σ+态的势能曲线

    Figure 1.  Potential energy curves for the X2Σ+, A2Π and B2Σ+ states of CO+.

    图 2  CO+分子离子X2Σ+, A2Π和B2Σ+态的改进与实验的振动能级之间误差的绝对值

    Figure 2.  Absolute value of errors of the vibrational energy levels between our refined results and those given by experiment for the X2Σ+, A2Π and B2Σ+ electronic states.

    图 3  一氧化碳阳离子的电偶极跃迁矩随核间距R的变化

    Figure 3.  Transition dipole moments for different states of CO+ cation as a function of internuclear distance R.

    图 4  CO+分子离子的配分函数

    Figure 4.  Partition functions of CO+.

    图 5  在100 atm压强时, CO+分子离子在不同温度下, 各个电子态之间的跃迁对整体不透明度的贡献 (a) 298 K; (b) 2000 K

    Figure 5.  Contributions of electronic state transitions to the total opacity of CO+ cation at different temperatures under pressure of 100 atm: (a) 298 K; (b) 2000 K.

    图 6  在100 atm 压强下, CO+分子离子X2Σ+, A2Π和B2Σ+态的不透明度随温度的变化(T = 298, 2000, 5000, 10000, 15000 K)

    Figure 6.  Opacities of CO+ cation for temperatures of 298 K, 2000 K, 5000 K, 10000 K and 15000 K at a pressure of 100 atm including transitions of B2Σ+-X2Σ+, B2Σ+-A2Π and A2Π-X2Σ+.

    表 1  CO+分子离子X2Σ+, A2Π和B2Σ+态的振动能级间隔(单位cm–1)

    Table 1.  Vibrational energy level intervals for X2Σ+, A2Π和B2Σ+ state of CO+ (in cm–1).

    v X2Σ+ A2Π B2Σ+
    Refined Exp[21] Refined Exp[21] Refined Exp[39]
    1 2183.96 2183.89 1535.08 1535.05 1679.49 1679.55
    2 2153.53 2153.56 1508.02 1508.03 1626.18 1626.61
    3 2123.18 2123.24 1481.08 1481.11 1575.66 1575.77
    4 2092.87 2092.92 1454.26 1454.29 1527.62 1527.02
    5 2062.56 2062.58 1427.55 1427.56 1481.78 1480.36
    6 2032.21 2032.20 1400.93 1400.92 1437.90 1435.78
    7 2001.79 2001.79 1374.40 1374.40 1395.74
    8 1971.26 1971.33 1347.94 1347.95 1355.11
    9 1940.61 1940.82 1321.56 1321.32 1315.81
    10 1909.80 1910.26 1295.25 1295.09 1277.70
    11 1878.80 1879.68 1269.00 1268.82 1240.61
    12 1847.60 1849.06 1242.80 1242.60 1204.42
    13 1816.18 1818.42 1216.65 1217.00 1169.02
    14 1784.50 1787.79 1190.54 1190.00 1134.28
    15 1752.56 1757.19 1164.47 1100.11
    16 1720.33 1726.60 1138.43 1066.44
    17 1687.79 1695.30 1112.43 1033.17
    18 1654.93 1666.00 1086.45 1000.24
    19 1621.72 1632.00 1060.50 967.57
    20 1588.14 1601.00 1034.56 935.11
    下载: 导出CSV

    表 2  CO+分子离子的光谱常数

    Table 2.  Spectroscopic constants of the doublet Λ-S states for CO+.

    Λ-S statesTe/cm–1Reωe/cm–1ωeχe/cm–1Be/cm–1De/cm–1
    X2Σ+Refined01.11522214.7215.23121.976767533.4167
    Theory[30]01.1192214.614.7566731
    Theory[39]2212.0517.741.98056167813.483
    Theory[40]01.11532214.1515.151.9769
    Exp[30]01.11512215.115.2767535
    Exp[39]2214.3904(10)15.1304(31)1.976964(12)
    Exp[41]2214.2219(79)15.1509(31)1.976949(26)
    A2ΠRefined20733.00151.24371561.1713.32361.589546186.90
    Theory[30]205941.2461570.012.8646473
    Theory[40]19628.21.24381562.7913.931.5894
    Exp[30]20733.31.24371562.013.5347128
    B2Σ+Refined45876.69891.16891712.8422.32201.798937988.4876
    Theory[30]459791.1701742.725.9937123
    Theory[39]45868.201742.9419.611.80518437854.825
    Exp[30]45876.71.16871734.127.9237692
    Exp[39]45878.204(95)1734.480(98)28.033(90)1.799491(21)
    Exp[41]45876.724(48)1734.626(86)28.272(38)1.799526(20)
    下载: 导出CSV
  • [1] Liang G Y, Peng Y G, Li R, Wu Y, Wang J G 2020 Chin. Phys. B 29 23101 doi: 10.1088/1674-1056/ab5fb6
    [2] Li R, Sang J Q, Lin X H, Li J J, Liang G Y, Wu Y 2022 Chin. Phys. B 31 103101 doi: 10.1088/1674-1056/ac834f
    [3] Liang G Y, Peng Y G, Li R, Wu Y, Wang J G 2020 Chin. Phys. Lett. 37 123101 doi: 10.1088/0256-307X/37/12/123101
    [4] 陈晨, 赵国鹏, 祁月盈, 吴勇, 王建国 2022 物理学报 71 143102 doi: 10.7498/aps.71.20220043 Chen C, Zhao G P, Qi Y Y, Wu Y, Wang J G 2022 Acta Phys. Sin. 71 143102 doi: 10.7498/aps.71.20220043
    [5] Xiao L, An S, Liu D, Minaev B F, Ågren H, Yan B 2025 Spectrochim. Acta, Part A 330 125704 doi: 10.1016/j.saa.2025.125704
    [6] Herzberg G 1950 Molecular Spectra and Molecular Structure: I (New York: Van Nostrand Reinhold Company) pp240–491
    [7] Davies P B, Rothwell W J 1985 J. Chem. Phys. 83 5450 doi: 10.1063/1.449715
    [8] Krishna Swamy K S 1986 Earth Moon Planets 34 281 doi: 10.1007/BF00145086
    [9] Haridass C, Prasad C V V, Reddy S P 2000 J. Mol. Spectrosc. 199 180 doi: 10.1006/jmsp.1999.7995
    [10] Krupenie P H 1966 The Band Spectrum of Carbon Monoxide (Gaithersburg, MD: National Institute of Standards and Technology) pp1–37
    [11] Huebner W F 1990 Physics and Chemistry of Comets (Berlin Heidelberg: Springer) pp245–303
    [12] Haridass C, Prasad C V V, Reddy S P 1992 ApJ 388 669 doi: 10.1086/171183
    [13] Bembenek Z, Domin U, Kepa R, Porada K, Rytel M, Zachwieja M, Jakubek Z, Janjic J D 1994 J. Mol. Spectrosc. 165 205 doi: 10.1006/jmsp.1994.1123
    [14] Zhuang H, Yang X H, Wu S H, Bi Z Y, Ma L S, Liu Y Y, Chen Y Q 2001 Mol. Phys. 99 1447 doi: 10.1080/00268970110062288
    [15] Yang X H, Wu Y D, Chen Y Q 2007 J. Mol. Spectrosc. 245 84 doi: 10.1016/j.jms.2007.06.007
    [16] Baltzer P, Lundqvist M, Wannberg B, Karlsson L, Larsson M, Hayes M A, West J B, Siggel M R F, Parr A C, Dehmer J L 1994 J. Phys. B: At. Mol. Opt. Phys. 27 4915 doi: 10.1088/0953-4075/27/20/011
    [17] Huber K P, Herzberg G 1979 Molecular Spectra and Molecular Structure 4: Constants of Diatomic Molecules (New York: van Nostrand Reinhold Comp
    [18] Reddy R R, Viswanath R 1990 J. Astrophys. Astron. 11 67 doi: 10.1007/BF02728023
    [19] Kępa R, Malak Z, Szajna W, Zachwieja M 2003 J. Mol. Spectrosc. 220 58 doi: 10.1016/S0022-2852(03)00059-6
    [20] Wu Y D, Yang X H, Guo Y C, Chen Y Q 2008 J. Mol. Spectrosc. 248 81 doi: 10.1016/j.jms.2007.11.013
    [21] Coxon J A, Kępa R, Piotrowska I 2010 J. Mol. Spectrosc. 262 107 doi: 10.1016/j.jms.2010.06.001
    [22] Hamilton P A, Hughes A N, Sales K D 1993 J. Chem. Phys. 99 436 doi: 10.1063/1.465767
    [23] Araújo J P, Ballester M Y, Lugão I G, Silva R P, Martins M P 2024 J. Mol. Model 30 352 doi: 10.1007/s00894-024-06131-8
    [24] Tobias I, Fallon R J, Vanderslice J T 1960 J. Chem. Phys. 33 1638 doi: 10.1063/1.1731475
    [25] Krupenie P H, Weissman S 1965 J. Chem. Phys. 43 1529 doi: 10.1063/1.1696966
    [26] Singh R B, Rai D K 1966 J. Mol. Spectrosc. 19 424 doi: 10.1016/0022-2852(66)90265-7
    [27] Locht R 1977 Chem. Phys. 22 13 doi: 10.1016/0301-0104(77)85204-X
    [28] Honjou N, Sasaki F 1979 Mol. Phys. 37 1593 doi: 10.1080/00268977900101141
    [29] Rosmus P, Werner H J 1982 Mol. Phys. 47 661 doi: 10.1080/00268978200100482
    [30] Okada K, Iwata S 2000 J. Chem. Phys. 112 1804 doi: 10.1063/1.480743
    [31] Mezei J Z, Backodissa-Kiminou R D, Tudorache D E, Morel V, Chakrabarti K, Motapon O, Dulieu O, Robert J, Tchang-Brillet W Ü L, Bultel A, Urbain X, Tennyson J, Hassouni K, Schneider I F 2015 Plasma Sources Sci. Technol. 24 035005 doi: 10.1088/0963-0252/24/3/035005
    [32] Weck P F, Schweitzer A, Kirby K, Hauschildt P H, Stancil P C 2004 Astrophys. J. 613 567 doi: 10.1086/423029
    [33] Mourik T V, Wilson A K, Dunning Jr T H 1999 Mol. Phys. 99 529 doi: 10.1080/00268979909482990
    [34] Werner H, Knowles P J, Knizia G, Manby F R, Schütz M 2012 WIREs Comput. Mol. Sci. 2 242 doi: 10.1002/wcms.82
    [35] Le Roy R J 2017 J. Quant. Spectrosc. Ra. Transfer 186 167 doi: 10.1016/j.jqsrt.2016.05.028
    [36] Yurchenko S N, Lodi L, Tennyson J, Stolyarov A V 2016 Comput. Phys. Commun. 202 262 doi: 10.1016/j.cpc.2015.12.021
    [37] Yurchenko S N, Al-Refaie A F, Tennyson J 2018 Astron. Astrophys. 614 A131 doi: 10.1051/0004-6361/201732531
    [38] https://physics.nist.gov/cgi-bin/ASD/levels_pt.pl [2025-3-11]
    [39] Ventura L R, Da Silva R S, Ballester M Y, Fellows C E 2020 J. Quant. Spectrosc. Ra. Transfer 256 107312 doi: 10.1016/j.jqsrt.2020.107312
    [40] Reddy R R, Nazeer Ahammed Y, Rama Gopal K, Baba Basha D 2004 J. Quant. Spectrosc. Ra. Transfer 85 105 doi: 10.1016/S0022-4073(03)00195-X
    [41] Szajna W, Kepa R, Zachwieja M 2004 Eur. Phys. J. D 30 49 doi: 10.1140/epjd/e2004-00060-0
    [42] Hakalla R, Szajna W, Piotrowska I, Malicka M I, Zachwieja M, Kępa R 2019 J. Quant. Spectrosc. Ra. Transfer 234 159 doi: 10.1016/j.jqsrt.2019.05.030
    [43] Szajna W, Kȩpa R, Field R W, Hakalla R 2024 J. Quant. Spectrosc. Ra. Transfer 324 109059 doi: 10.1016/j.jqsrt.2024.109059
  • 加载中
图( 6) 表( 2)
计量
  • 文章访问数:  73
  • HTML全文浏览数:  73
  • PDF下载数:  1
  • 施引文献:  0
出版历程
  • 收稿日期:  2025-03-24
  • 刊出日期:  2025-06-20

CO+分子离子X2Σ+, A2Π和B2Σ+态的不透明度

    通讯作者: E-mail: xiazhang@jlu.edu.cn.; 
    通讯作者: E-mail: yanbing@jlu.edu.cn.
    作者简介: 安斯腰力吐:3564414330@qq.com
  • 吉林大学原子与分子物理研究所, 长春 130012

摘要: 一氧化碳分子离子(CO+)在大气及天体物理环境中起着关键作用, 其不透明度的理论研究对辐射输运建模具有重要意义. 本文基于实验观测数据, 采用Modified-Morse (MMorse)势函数改进并构建了CO+分子离子X2Σ+, A2Π和B2Σ+电子态的势能曲线, 进一步提取了振动能级和光谱常数. 同时, 利用考虑Davidson修正的多参考组态相互作用(MRCI+Q)方法计算了势能曲线和电偶极跃迁矩. 改进获得的MMorse势与计算得到的势能曲线非常吻合, 且光谱常数和振动能级与其他理论和实验数据符合较好. 结合MMorse势函数和从头计算获得的电偶极跃迁矩, 计算了CO+分子离子在100 atm (1 atm=1.01×105 Pa)压强下, 298—15000 K温度范围内的不透明度, 并探究了不同温度对高温谱的影响. 研究结果表明, 高温环境(T > 5000 K)会导致不同能带系统的谱线展宽与边界模糊, 这种混合效应在T > 10000 K时尤为突出, 揭示了高温下分子离子光谱退化的微观机制. 本研究可以为天体物理领域提供一些理论依据和数据支持. 本文数据集可在https://doi.org/10.57760/sciencedb.j00213.00136中访问获取.

English Abstract

    • 天体物理学中, 辐射不透明度是一种描述物质对辐射吸收能力强弱的量, 通常用来描述光在物质中的传播特性, 其数值主要依靠理论计算得到[15]. 一氧化碳阳离子(CO+)在双原子分子光谱学研究领域中具有重要的科学意义, 其光谱特征为揭示彗星尾部、星际分子云、恒星大气、太阳外层以及行星空间环境中的光化学过程、分子动力学机制及能量传递规律提供了关键的理论和实验依据[68]. 此外, CO+分子离子在环境监测和燃烧化学动力学研究中的应用, 进一步突显了其在科学研究中的广泛意义[9]. 为此, 研究者们通过实验与理论相结合的方法, 对该阳离子的结构、光谱特征及动力学行为进行了系统性研究, 从而为其提供了更为精确和全面的分子数据.

      大量的研究成果表明, CO+分子离子的光谱特征已被深入探讨[10]. CO+分子离子拥有4个重要的带系, 这些带系分别对应于X2Σ+, B2Σ+, A2Πi和C2Δr电子态之间的跃迁. 具体而言, 研究者们已经观测到并分析了多个跃迁带系, 如彗尾带系(A2Πi-X2Σ+)、第一负带系(B2Σ+-X2Σ+)、Baldet-Johnson带系(B2Σ+-A2Πi)以及Marchand-D'Incan-Janin带系(C2Δr-A2Πi)等, 这些为研究CO+分子离子在宇宙中的辐射行为和物理特性提供了宝贵数据[9]. 1990年, 通过微波光谱学等技术, Huebner[11]已对CO+的彗尾带系进行了详细的分析. CO+的彗尾带系出现在光谱区域 11765—32468 cm–1 中, 可用于研究超音速宇宙飞船以逃逸速度飞行时的辐射光谱及热量传输和金星的大气边缘的气体状态, 其中含有大量CO2, CO+为其中物理和化学过程的产物之一. Haridass等回顾了早期对CO+的4个电子跃迁谱带及其红外和微波光谱的研究工作[12], 并于2000年重新研究了12C16O+的彗尾带系(A2Πi-X2Σ+)[9].

      在这些研究中, 使用了各种实验技术, 如红外激光光谱[7]、发射光谱[13]、光外差探测[14], 速度调制光谱技术[14,15]和光电子能谱[16]. 1966 年, Krupenie[10]回顾了当时的实验工作. 1979 年, Huber 和 Herzberg[17]编评了CO+实验光谱参数和分子常数. 1990 年, Reddy和Viswanath[18]通过改进RKR势能重构算法, 计算了CO+分子离子X2Σ+和B2Σ+电子态的高精度势能曲线和基态解离能D0. 1994年, Bembenek等[13]识别和记录了0-2, 1-3和2-1波段并分析了B2Σ+-A2Πi跃迁的1-2和2-0波段. 2003年, Kępa等[19]总结了B2Σ+-X2Σ+跃迁. Wu等[20]在2008年和Coxon等[21]在2010年总结了当时CO+阳离子的各种同位素异构体的A2Πi-X2Σ+和B2Σ+-X2Σ+跃迁.

      在以前的理论工作中, 已经对CO+分子离子进行了许多电子结构光谱计算, 已经确定了许多电子态的势能曲线(PECs)[22,23]. 初步研究始于1960年, 当时Tobias等[24]使用 Rydberg-Klein-Rees (RKR)方法推导出了CO分子的 X1Σ+, a3Π, a'3Σ+, d3∆, e3Σ, A1Π 和 B1Σ+电子态的势能曲线. 1965年, Krupenie和Weissman[25]使用类似的方法, 除了Tobias等[24]计算的曲线外, 还确定了CO+双原子阳离子(X2Σ+, A2Π和B2Σ+态)的势能曲线. Singh和Rai[26]对CO+阳离子的相关研究中, 结合 RKR 方法和原子分子轨道线性组合自洽场(LCAO-MO-SCF)方法获得了比 Krupenie和Weissman[25]之前的研究更准确的势能曲线. Locht[27]在1977 年计算了其他电子态, 如C2Σ+, E2Σ+, F2Σ+和D2Π态. 1979年, Honjou和Sasaki[28]使用MO-SCF-CI方法确定了CO+分子离子的5条解离限. Rosmus 和 Werner[29]在1982年使用多组态自洽场方法(multiconfiguration self-consistent field, MCSCF)计算了X2Σ+和A2Π电子态之间的电偶极跃迁矩. 2000年, Okada和Iwata[30]使用MRCI方法与AVQZ 基组计算了CO+的基态和低能级的几个二重态和四重态的势能曲线和电偶极跃迁矩. 2015年, Mezei等[31]研究了CO+振动激发的解离复合.

      尽管CO+分子离子的光谱特性已在大量实验中被广泛观测, 现有理论研究仍主要集中在CO+分子离子的光谱常数和分子结构计算上. 值得注意的是, 关于大气光谱诊断中至关重要的辐射特性, 目前仍缺乏系统性的研究. 本研究选定典型的冷星大气参考压强100 atm (1 atm = 1.01×105 Pa)[32], 系统地研究了298—15000 K温度范围内CO+分子离子的不透明度. 通过这一研究, 我们期望能够为未来的天文观测和空间探测任务提供更加可靠的分子数据, 进一步推动天体物理学和分子光谱学的交叉发展.

    • 本文的计算中使用Hartree-Fock方法计算了CO+分子离子的分子轨道, 然后采用态平均完全活性空间自洽场(CASSCF)和考虑Davidson修正的多参考组态相互作用(MRCI+Q)方法, 在aug-cc-pV6Z基组[33]下计算了CO+分子离子的X2Σ+, A2Π和B2Σ+电子态的电子结构. 在相同的理论方法下, 分别使用软件包MOLPRO 2012[34]和LEVEL程序[35]计算了CO+离子的电子结构和光谱常数.

      CO+分子离子的对称点群是Cv, 但在Molpro 2012 程序包中使用的是阿贝尔点群, 所以计算中采用了C2v点群对称性. C2v点群包含4个不可约表示( A1, B1, B2, A2), 与分子轨道的对应为: Σ+ = A1, Π = B1 + B2, Σ = A2. 在CASSCF计算中, 选取最外层的8个分子轨道作为活性空间, 包括4个A1、2个B1和2个B2分子轨道. 因此, C+原子在2s2p轨道上的3个电子, O原子2s2p轨道上的6个电子被放在活性空间, 即这9个电子分布在8个分子轨道上. 这些态的从头算核间距都选择在0.8—6 Å的范围之间.

      在Born-Oppenheimer近似条件下, 双原子分子的电子Schrödinger方程可以写成:

      其中$ {\widehat{H}}_{{\mathrm{e}}} $是双原子分子N个电子的哈密顿量, $ {\varPsi }_{{\mathrm{e}}}(r, R) $$ {E}_{{\mathrm{e}}}\left(R\right) $分别是特征函数和特征能量.

      (2)式用a, b和i, j分别对原子核和电子进行标识; Rr分别代表原子核和电子之间的距离; ZaZb是原子核a和b的电荷数; $ {r}_{{\mathrm{a}}i}{\mathrm{和}}{r}_{{\mathrm{b}}i} $分别是原子核a和第i个电子之间的距离、原子核b和第i个电子之间的距离, $ {r}_{ij} $则是第i个电子和第j个电子之间的距离.

      根据计算的电子态能量, 可以通过求解核的Schrödinger方程来获得振转态的特征值和波函数:

      其中$ \mu $是约化质量, $ {E}_{vJ}\left(R\right) $是振转能级, $ {\chi }_{vJ} $是核的波函数.

      结合势能曲线(PECs)和电偶极跃迁矩(TDM), 可以获得分子电子态的跃迁性质. Franck-Condon因子和爱因斯坦系数的描述如下:

      其中$ {v}'{J}' $$ {v}''{J}'' $代表上能级和下能级, $ {\varPsi }_{{v}'{J}'} $$ {\varPsi }_{{v}''{J}''} $分别是上能级和下能级的波函数, $ \nu $是以cm–1为单位的发射频率, M(r)为电偶极跃迁矩, $ S({J}', {J}'') $是能级的Hönl-London因子.

    • 基于实验观测的振动能级数据, 采用最小二乘法拟合[36]结合部分实验振动能级数据改进(refined)分子势能曲线. 随后, 基于分子势能曲线确定了光谱参数和振动能级. 我们所使用的Modified-Morse势函数的解析表达式如下:

      其中, $ \xi =(r-{r}_{{\mathrm{e}}})/(r+{r}_{{\mathrm{e}}}) $, $ {r}_{{\mathrm{e}}} $是平衡核间距, $ {a}_{i} $是系数. $ {A}_{{\mathrm{e}}} $是指势函数$ V\left(r\right) $r趋于无穷大时的渐近线, 与电子最低能态$ \left({T}_{{\mathrm{e}}}\right) $有关. 渐近线$ {A}_{{\mathrm{e}}} $可以用来推算出给定电子态的解离能$ D_{\rm e}={A}_{{\mathrm{e}}}-{T}_{{\mathrm{e}}} $.

    • 结合改进的势能曲线和计算的电偶极跃迁矩, 通过分子光吸收截面计算获得了CO+分子离子的配分函数和不透明度[37]. 用于跃迁 λfλi 的吸收线强度(单位: cm/molecule)计算公式为

      其中, $ {A}_{{\mathrm{f}}{\mathrm{i}}} $是爱因斯坦系数(单位: s–1 ), gns是核统计权重因子, $ {J}_{{\mathrm{f}}} $是f电子态的转动量子数, $ {c_2} = hc/{k_{\text{B}}} $是二阶辐射常数, $ {\tilde E_i} = {E_i}/hc $是项值, $ {\tilde{v}}_{{\mathrm{f}}{\mathrm{i}}} $是跃迁$ {\mathrm{f}}\leftarrow {\mathrm{i}} $的波数, T是环境的温度(单位: K). 对于异核分子, gns是核自旋组合的总数, 计算公式为$ {g_{{\text{ns}}}} = {\text{ }}\left( {2{I_{\text{a}}} + 1} \right)\left( {2{I_{\text{b}}} + 1} \right) $, 其中 IaIb分别是相应的核自旋. $ Q $(T )是配分函数, 其表达式是

      跃迁的吸收截面$ {\sigma }_{{\mathrm{f}}{\mathrm{i}}}\left(\tilde{v}\right) $与相应的积分吸收系数Ifi之间的关系为

      式中$ \tilde{v} $是波数. 引入线型谱函数$ {f}_{{\tilde{v }}_{{\mathrm{f}}{\mathrm{i}}}}\left(\tilde{v }\right) $后, 吸收截面(单位: cm2/molecule)定义为

      其中$ {f}_{{\tilde{v }}_{{\mathrm{f}}{\mathrm{i}}}}\left(\tilde{v}\right) $是面积归一化为1的可积函数,

    • 我们使用上述理论方法计算了两组CO+的3个电子态X2Σ+, A2Π和B2Σ+的势能曲线, 一组为使用Modified-Morse势函数结合实验光谱数据拟合的改进的势能曲线, 另一组为使用MRCI+Q结合aug-cc-pV6Z基组计算的从头计算势能曲线, 如图1所示(从头计算数据在图1中用点代表, 改进的数据是以实线表示). 结果表明, 改进后的势能曲线在整体形态上呈现良好的物理光滑性, 并与理论计算值展现出高度一致性. 图1中显示, X2Σ+, A2Π和B2Σ+三个电子态都是典型的束缚态, 并且CO+分子离子的第一和第二条解离限分别为C+(2Pu) + O(3Pg)和C+(2Pu) + O(1Dg). 这两条解离限的能量差值为16015.52 cm–1, 与实验值[38]仅有147.66 cm–1的偏差(0.92%). 这为后续光谱分析和分子常数研究提供了理论基础.

      本文给出了CO+分子离子改进的结果与实验的X2Σ+, A2Π和B2Σ+电子态的振动能级之间误差的绝对值, 如图2所示. 由图2可知, 与实验结果[23]相比, 最大误差仅有0.15 cm–1. 此外, CO+分子离子的X2Σ+, A2Π和B2Σ+三个电子态的振动能级间隔数据列于表1, 与实验结果[21,39]相比, 最大误差分别为0.81%, 0.05%和0.15%. CO+离子分子的3个电子态的振动能级间隔随振动量子数增大而逐渐减小, 这是由于势能曲线的非谐性而导致的.

      CO+分子离子的X2Σ+, A2Π和B2Σ+态相应的光谱常数是通过对一维核薛定谔方程进行数值求解而获得[34]. 如表2所列, 计算值涵盖了绝热激发能Te、平衡核间距Re、振动常数ωeωeχe、转动常数Be和解离能De等关键参数, 并与实验值及理论结果进行系统性对比. 对于X2Σ+电子态, Re的值为1.1152 Å, 与其他理论值[40]和实验[30]符合得很好, ωe, ωeχeBe与最新的实验值[39]相比仅有微小差异, 偏差分别为0.12%, 0.53%和0.01%, 显著优于采用aug-cc-pV6Z基组的近期理论[39]. 激发态A2Π的ωe, ωeχeDe与实验值[30]的偏差分别为0.05%, 1.55%和2.04%, 其中TeRe值相同, 揭示芯-价电子关联效应的高精度描述能力. 激发态B2Σ+ωeχe与实验值[39]存在显著差异(Δ = 7.2485 cm–1), 但其与aug-cc-pV6Z理论预测值仅偏离2.7120 cm–1; 其余参数Te, ωe, BeDe与实验值的相对偏差均低于1%, 这验证了我们改进的势能曲线模型的可靠性.

      作为本工作的一部分, 计算了CO+分子离子X2Σ+, A2Π和B2Σ+三个电子态之间的电偶极跃迁矩 (B2Σ+-X2Σ+, B2Σ+-A2Π和A2Π-X2Σ+ ), 如图3所示. 由图3可知, 各电偶极跃迁矩在平衡核间距附近(1.05—1.30 Å)呈现显著极值特征, 其中尤为明显的是B2Σ+←X2Σ+跃迁在r = 1.15 Å处达到峰值0.31 a.u. 当核间距r > 3 Å时, 所有电偶极跃迁矩均趋近于零, 这与解离产物C+(³P) + O(³P)的自旋-轨道对称性约束严格相符, 并与过去的理论结果符合较好[30]. 这3个电偶极跃迁矩中B2Σ+←X2Σ+跃迁对总吸收截面的贡献更大.

    • 基于CO+分子离子的势能曲线和所有振-转参数, 本研究精确确定了其配分函数和不透明度. 通过拟合实验数据获得的势能曲线和MRCI+Q方法计算的电偶极跃迁矩, 在局域热力学平衡(LTE)条件下, 依据(8)式系统计算了15—15000 K温度范围内CO+分子离子的配分函数Q(T)随温度的变化, 如图4所示. 作为光谱分析的重要参数, 配分函数由振转能级、简并度和温度共同决定, 并直接关联吸收或发射截面的计算精度. 计算结果表明: Q(T )随温度的升高呈单调递增趋势, 在T = 15—15000 K范围内增长跨越4个数量级, 其物理本质源于高温环境下激发态的布居数增长.

      为了证明CO+分子离子的光谱模拟的准确性, 与实验测量结果进行对比. 图5给出了CO+分子离子在100 atm压强下, 温度为298 K和2000 K时, 各个电子态之间的跃迁对整体不透明度的贡献. B2Σ+-X2Σ+是该体系中最强的电子谱带, 也是实验上观测到的最多、最容易探测的谱带. 此外, A2Π-X2Σ+带系的贡献整体强于B2Σ+-A2Π, 但根据我们的模型, 后者仍在15000—22000 cm–1区域可探测到, 这与实验[42]对这些振动带的观测结果一致. 此外, 在40000 cm–1以上的区域, B2Σ+-X2Σ+带占主导地位, 这也与Szajna等[43]在这个区域观测的振动带相符.

      CO+分子离子不透明度是在上述分子结构数据, 以及电偶极跃迁矩和配分函数的基础上进行计算的. 图6显示了在100 atm 的压强条件下, 不同温度(298 K, 2000 K, 5000 K, 10000 K和15000 K)对CO+分子离子X2Σ+, A2Π和B2Σ+电子态的不透明度的影响. 其中最大的截面出现在T = 298 K且波数($ \tilde{v} $)为 45440 cm–1的范围, 随着温度的升高, 不透明度的最大峰值逐渐减小. CO+分子离子的不透明度分布在红外波段到紫外波段, 随着温度的升高, 激发态A2Π和B2Σ+上的布居数增大, 从而使长波截面增大. 这主要来自于激发态之间的跃迁, 即B2Σ+-A2Π. 在不透明度中可以看到3个明显的主峰, 波长从大到小, 分别对应B2Σ+-X2Σ+, B2Σ+-A2Π和A2Π-X2Σ+跃迁的贡献. 由图6还可以发现, 随着温度的升高, 不同的谱带相互融合. 温度达到10000 K时, 3个主峰对的界限也不再清晰. 这些现象也是由于电子激发态和振转激发态的布局随着温度升高而加强导致的. 正如预期, 光谱在较高温度下变得更加明显和更宽. 该模型成功复现了极端天体物理环境中分子光谱的退分辨特征, 为系外行星大气及恒星风相互作用的辐射传输模拟提供了关键参数集.

    • 本文基于实验数据使用MMorse解析式改进了CO+分子离子X2Σ+, A2Π和B2Σ+电子态的势能曲线, 并给出了振动能级、光谱常数、配分函数和不透明度. 振动能级间隔和光谱常数计算结果与理论结果和现有的实验值符合得很好, 可以保证振转能级的准确性; 所计算的配分函数与不透明度包含了文中所提到的X2Σ+, A2Π和B2Σ+三个电子态振动能级的贡献; 根据改进的势能曲线和理论计算的电偶极跃迁矩, 给出了100 atm压强条件下, 温度为298 K, 2000 K, 5000 K, 10000 K和15000 K时的不透明度. 结果表明, CO+分子离子的不透明度分布在红外波段到紫外波段; 随着温度的升高, 激发态的布居数增加, 不透明度的波长分布范围变宽. 当温度大于10000 K时, 导致不同能带系统的谱线展宽与边界模糊化, 揭示了高温中分子离子光谱退化的微观机制. 本文所得数据可应用于大气混合物的不透明度及辐射输运过程研究, 为天体物理学领域提供一些理论依据和数据支持.

    • 支撑本研究成果的数据集可在科学数据银行https://doi.org/10.57760/sciencedb.j00213.00136中访问获取.

    参考文献 (43)

目录

/

返回文章
返回