H_%282%29S和H_%282%29S%5E%28%2B%29的振动光谱和电子能谱计算.pdf

第41卷第4期2024年8 月J.At.Mol.Phys.,2024,41:041005(8pp)原子与分子物理学报JOURNAL OF ATOMIC AND MOLECULAR PHYSICSH,S和H,S+的振动光谱和电子能谱计算Vol.41 No.4Aug.2024郭甜，徐建刚，张佳美，张思雨，陈艳南，张云光（西安邮电大学理学院，西安7 10 12 1）摘要：本文基于MOLPRO软件包使用从头算方法计算了星际分子H,S及其阳离子H,S+的势能面及光电子能谱：首先，在（U)CCSD/cc-pVQZ理论水平下获取了H,S沿法线坐标展开的势能面，势能面直观描述了不同振动模式耦合对分子能量变化的影响,S-H键的非对称伸缩振动和面内弯曲振动共同作用使得系统势能变化范围明显变大：振动多组态相互作用方法被用来计算非谐振动频率和振动光谱，计算结果显示，倍频和组合频之间出现了强烈的费米共振，使得相应波段处的红外强度显著增强最后，使用拉曼波函数和收缩不变Krylov子空间方法首次计算了H,SXA一H,S+XB,的光电子能谱：此项研究有助于进一步理解星际分子的内部结构，并为实验研究及星际观测提供参考，关键词：势能面；VMRCI；振动光谱；光电子能谱中图分类号：0 56 1Vibrational spectra and photoelectron spectra calculations of H,S and H,S文献标识码：AD0I:10.19855/j.1000-0364.2024.041005GUO Tian,XU Jian-Gang,ZHANG Jia-Mei,ZHANG Si-Yu,CHEN Yan-Nan,ZHANG Yun-Guang(School of Science,Xi an University of Posts and Telecommunications,Xi an 710121,China)Abstract:The potential energy surfaces and photoelectron spectra of the interstellar molecule H,S and its cationH,S+were calculated using ab initio methods based on the MOLPRO software package.First,the potential ener-gy surfaces of H,S expanded along the normal coordinates were obtained at the(U)CCSD/cc-pVQZ level oftheory.The potential energy surfaces visually describe the effect of different vibrational modes coupling on themolecular energy change,and the combined effect of the asymmetric stretching vibration of the S-H bond andthe in-plane bending vibration makes the potential energy change range of the system significantly larger.Thevibrational multi-reference configuration interaction method was used to calculate the anharmonic vibrational fre-quencies and vibrational spectrum,and the calculations showed that strong Fermi resonances appear between theovertones and combination bands,resulting in a significant enhancement of the infrared intensity at the corre-sponding bands.Finally,the photoelectron spectrum of H,S XA,-H,S+xB,was calculated for the first timeusing the Raman wave function and the contracted invariant Krylov subspaces method.This study will help to fur-ther understand the internal structure of interstellar molecules and provide a reference for experimental studiesand interstellar observations.Key words:Potential energy surface;VMRCI;Vibrational spectroscopy;Photoelectron spectroscopy中起着重要的作用，它的微波和红外光谱在许多1 引 言振动状态下显示出不寻常的谱带强度 例如，H,S分子的物理和化学性质在星际介质大气H,S的基频振动比类似的三原子分子（如H,O和收稿日期：2 0 2 2-11-17基金项目：国家自然科学基金（2 2 10 30 6 1）；陕西省自然科学基金（2 0 2 2 JQ-569）作者简介：郭甜（1997 一），女，在读硕士，主要研究方向：分子结构与分子光谱：E-mail：g d a t i a n 16 3.c o m通讯作者：张云光.E-mail：z y g s r 2 0 10 16 3.c o m041005-1第41卷H,Se)弱两到三个数量级，不仅如此，它的组合频V+V2、V 2+v 和1+V,也比类似分子弱得多2.3它还是一种轻质、非刚性的不对称转子分子，引起了众多学者的关注197 5年cook等人对H,S旋转光谱的微波研究，提供了精确的旋转常数和四次畸变常数4 1990 年，Kauppi和Halonen运用微扰理论计算了硫化氢分子的势能面5，并用非线性最小二乘法优化了势能参数之后，Poly-ansky提出对运动动能算符的非绝热修正6 Ko-zin和Jensen通过拟合高分辨率光谱学实验中观察到的随机能级7 ，使用莫尔斯振子刚性弯曲内部动力学代码优化了势能面：Martin于19 9 5年通过从头计算得到了电子基态的H,S势能面8 Xie等人采用VSCF-CI方法计算了H,S的高激发振动态的能量9.Cazzoli等人获得了硫化氢的两种同位素的旋转光谱10 .Lodi等人计算了H,S分子从头算偶极矩面，发现核心修正和相对论修正对谱线强度有显著贡献 Carvajal 等人提出了一种基于非谐阶梯算子分析分子振动结构的新模型，并将其应用于H,S 分子I.H,S由于其自身性质引发了人们对其更深层次的研究兴趣；另一方面，电子状态可以用于确定不同分子性质的光离过程，如电离能或结构参数，用于识别分子UV/VIS的电子激发或具有特征跃迁的荧光过程，如有机发光二极管已有较为强大的实验技术来测量和研究这些转变对理论方面而言，简单绝热跃迁的计算较为容易，但这些计算不能描述包含重要信息的光谱的精细结构12 对分子光电子能谱不同跃迁的光谱带的计算和分配变得越来越重要，振动贡献或垂直激励等性质可以用FranckCo n d o n f a c t o r s（FCFs）来解释13,14.因此，本文计算了H,S 及H,S*势能面及振动光谱，并基于初态（H,SXA,）和末态（H2SxB,）势能数据绘制相应光电子能谱.2理论计算本文所有计算在MOLPRO2018量子化学软件包中C2v点群下进行，以HartreeFo c k（H F)方程为起点，对H,S和H,S+分子轨道进行初始猜测，采用(非限制）单双耦合簇理论（U)CCSD结合相关一致基组cc-pVQZ对H,S和H,S+优化几何结构，使用同一层次理论计算了分子振动产生的谐振频率及一维、二维势能面，势能面使用多模展开，并在3D项之后截断.原子与分子物理学报V(q)=ZV(qi)+ZV(qi,i)+1Dk3DV.(qi)=V(q.)-Vo,V,(qi,q)=V(qi,q,)-ZV(qr)-Vo,(3)Vu(qi9i,h)=Vo(qi9j,k)-ZVa(qr,g.)-E(q.)-VoT,sei.j,kreijTs势能面V()沿法坐标q展开，V.表示平衡结构的总能量，高阶耦合项是对低阶贡献的修正一旦势能面确定，就可以通过求解薛定谔方程来确定振动波函数，沃森哈密顿的振动自洽场理论（VSCF）被用于确定并优化初始猜测得到的单模波函数，在这些初始计算中使用了2 0 个分布式高斯基.H=12在后续的计算中，使用Culot和Lievin提出的VMCSCF方法15 修正静态相关效应，并将VCI系数冻结，使用一个向量进行VMRCI161计算修正动态相关效应，得到H,S和H,S*的基频、倍频以及组合频最后使用拉曼波函数(RWF)方法和收缩不变Krylov子空间（CIKS）中基于残差的特征对计算算法(RACE)17,18)模拟了 H,S XA,H,S*xB,的光电子能谱.3结果与讨论3.1H,S及H,S*的空间几何结构H,S与水分子相似，在C2v点群下共享一组共同的对称适应分子轨道，也就是说，中心原子S的原子轨道与H,的1轨道的对称和反对称组合相互作用。表1列出了在(U)CCSD/cc-pVQZ 理论水平下优化得到的H,S、H,S*几何结构参数.与实验值的相比，本次计算得到H,S中S-H的键长略微高于实验值，仅相差0.0 0 0 6 A，0 H s H 键角偏离约0.5，一致性良好相比H,S，H,S*中的S-H键长和QHsH键角分别扩展至1.356 4A，93.3624，当H,S阳离子达到平衡状态时，键长和角度分别增加了约1.5%和1%此外，由于H2S+的平衡几何结构参数实验值的缺乏，我们还提供了其他学者的计算结果以供对比，S-H键长相差0.0 0 0 32 A，键角偏离约0.2 8 041005-2第4期2D(1)(2)TEi.j(4)82+V(g)(5)第41卷表1 H,S及H,S*的平衡几何参数Table 1Equilibrium geometric parameters of H,S and H,s*H,STsH/(A)1.3368CHSH/()92.60823.2势能面3.2.1一维势能面图1中上图和下图分别为面内弯曲振动和S-H键对称拉伸振动产生的一维非对称势能曲线，而S-H键的非对称拉伸振动产生了对称势能曲线，如图2 所示；当分子振动模式所代表的对称性不可约表示不同于A,时，会产生对称势能曲线。在能量大小的变化上，能明显看出S-H键的拉伸振动对分子势能影响较大，是面内弯曲振动的2 倍至3倍；而S-H键的非对称拉伸振动对分子势能变化的影响最大在对称势能曲线产生的过程中，使用分子对称性还会使加快计算速度，使CPU时间减少一倍.0.030.0250.020.0150.010.0050-50-40-30-20-100q,(A1)0.080.070.060.050.040.030.020.010-50-40-30-20-100102030q,(A1)图1H,S振动产生的非对称势能曲线Fig.1Asymmetric potential energy curves generated byH,S vibration3.2.2二维相对势能面H,S三个基频振动模式两两组合共可以产生郭甜，等：H,S和H,S*的振动光谱和电子能谱计算Exp.19H,S+1.33621.356492.1193.36241020304050041005-3第4期0.120.1Cal.20 0.081.3567293.08020.060.040.020-40-30-20-10010203040q;(B2)图2 H,S振动产生的对称势能曲线Fig.2Symmetric potential energy curve generated byH,S vibration三个二维势能面，势能面的展开采用法线坐标，横坐标数值的大小和正负代表分子处于振动状态时各原子空间位置的变化，每个二维势能面沿法坐标方向生成16(2 4）个网格点由于零势能处的选取，使得纵坐标势能产生了负值，让分子在振动时的能量变化更为明显，并且在两个横向坐标同时达到最大值时变化最大在图3中的二维相对势能面中，仍能明显看出在横坐标对称性不同于A,时，所在横轴沿坐标轴方向延伸的势能变化表现出对称性，如图3中(b)和(c)中横坐标方向的33.2.3二维完整势能面三模耦合项用电子密度碎片化来模拟，它包含了对势能的1D贡献使用参考几何RGHF方法对多维势能面展开中高阶项进行建模，该方法基于平衡结构的局域电子密度，其关键步骤是在给定的参考结构下，将电子密度定位到分子的原子上，然后将这些原子中心的部分密度转移到扭曲结构的原子位置图4展示了H,S振动产生的二维完整势能面，其中，（a)为沿S-H键对称拉伸振动(q2)和H-S-H面内弯曲模式(qi）法线坐标展开的2 D势能面的特征，仅在横坐标同时为O时分子处于能量最低状态，S-H键对称拉伸振动较H-S-H面内弯曲振动分子能量变化更为明显；图4（c)中的势能变化特征类似而对图4(b)而言，S-H非对称伸缩模式(qs)和H-S-H面内弯曲模式（q2）振动模式共同作用使得纵坐标分子能量变化范围明显变大，且分子势能在q3横坐标正向和负向最大时达到最大值，也就是说，此时分子振动最为剧烈.第41卷(a)0.0050-0.005-0.01-0.015-0.02-0.025-0.033020100-10原子与分子物理学报(b)0.0050-0.005-0.01-0.005-0.01-0.015-0.0150.02-0.02-0.0230.02503-0.0350.03-0.04-50-400-305-1020第4期0.00500.005-0.0050-0.01-0.015-0.02-0.025-0.03-0.035-0.04-50404030-302010001020-20q,(A.)-30-40-5050403020100q,(Al)-10q,(B.)-20-30-40.5040302010q,(A,)(c)0.250.20.250.20.150.10.050-0.05-0.1403020100-10q,(B.)-20-30-4030图3H,S振动产生的二维相对势能面Fig.3 Two-dimensional relative potential energy surfaces generated by H,S vibration(a)0.0.0.8:820.013020100-10-20q2(A,)-30-40-50500.150.10.050-0.05-0.1-50-40-30-200-10q,(A.)1020(b)0.110.090.10.080.40.070.350.060.30.050.250.040.20.030.150.020.10.010.0500-50-30-40-20-100qi(A.)0302010400.40.350.30.250.20.150.10.050-5040-4030-3020-2010-1000q(A.)q:(B,)-10-20-30-40301020(c)0.120.10.120.10.080.060.040.0200.080.06.0.040.020-50403020100-10q:(B,)-20-30-405040300-40-300102010qi(A.)图4H,S振动产生的二维完整势能面Fig.4 Two-dimensional complete potential energy surfaces generated by H,S vibration在势能面的计算中，对称性可以被用在单个表面的构造中以节约时间和资源，在势能变化中，S-H键的拉伸振动对分子体系的能量影响相对较大，并且不同振动模式的耦合对也会影响分子的041005-4第41卷势能变化，这些依赖于直线法线坐标的多维高级势能面都是通过MOLPRO完全自动生成.3.3振动频率H,S分子三个基频振动模式如图1所示，Vi为S-H键的对称拉伸的振动，V2为OHsH 的面内弯曲振动，V3为S-H键的非对称拉伸振动对比计算得到的谐振频率和VSCF修正的结果，可以得出H,S在VSCF计算后对振动模式间非谐效应的对三个振动模式的修正分别为：6 6.6 7 cml、27.85cml、8 4.8 4c ml；H,S*在VSCF计算后非谐效应的修正大概为：31.13cm=l、6 5.0 7cm-l、8 5.12 c m=1在不同方法下获得的H,S及H,S*的基频如表2 所示：H,S中VMRCI的计算结果与实验结果之间的平均绝对误差约为2 1.6 9cm-，而VSCF的结果与实验值之间的平均绝对误差约39.8 1cm=1，说明引入对相关效应的校正后使得误差减少至二分之一，此外，使用VMRCITable 2Harmonic and and anharmonic vibration frequencies(cm-)of H,S and H,S+H2S+ModesSym.AV2AiV3B23.4振动光谱H,S振动产生的红外光谱中，最明显的吸收峰出现在1195cm处，由?sn的面内弯曲振动产生，如图6 所示；而第二强烈的吸收峰出现在图中的近红外光谱区域，这是由于组合频v2+V3(5203.45cm=）和倍频2 vi（52 0 3.48 c m）之间的耦合，相应能级间的带隙较小（仅相差0.0 3cm-），产生费米共振，使得为对应波段处的红外吸收强度变大在2 6 0 0 cm=1附近有两个吸收峰，高低峰分别由S-H键的非对称和对称拉伸振动产生；2 38 0.37 cm-1处出现吸收峰来源于倍频2 v2在图7 的拉曼光谱中能明显看出三个强度较高的峰值，其中，2 6 0 0 cm-附近显示了最强烈的拉曼活性，对比而言，相同频率处红外光谱显示出的吸收峰较为微弱，高低峰同样分别由分子基频振动中 S-H键的非对称和对称拉伸振动产生，另外一个峰值出现1195cm-处H-S-H键的面内弯曲振动，另有2 36 2.0 8 cm-1、38 7 1.2 8郭甜，等：H,S和H,S*的振动光谱和电子能谱计算表2 H,S及 H,S*的谐振频率和非谐振动频率（cm-）H2SHarmVSCF2742.452675.781224.411196.562756.742671.90第4期引入相关动态效应的校正对S-H键的拉伸振动的结果更为明显最精确的理论结果和实验数据之间显示出一定的偏差很大程度上是因为在计算中忽略了势能展开中的高阶耦合项，这可能导致基本振动产生约10 cm的误差表3汇总了倍频和组合频的计算结果，经过与有限的实验值对比，一致性较好，此外，表中2 v和+V振动带之间发生耦合，相应带隙小于1cm-1，从而产生费米共振对于H,S*而言，失电子使得其与H,S在相同振动模式下时的基频、倍频以及组合频减小，只有在组合频2 v2+2v较H,S振动频率增大.图5H,S的基频振动模式Fig.5 The vibration modes of H,SVMRCIExp.21,22641.362614.41195.951182.62652.572628.4cm1、532 9.40 c m处分别出现极为微弱的活性。强度拉曼光谱作为红外光谱的补充，帮助我们更准确地分析分子的振动结构：3.5光电子能谱图8 中的(a)和(b)分别为 T=300 K时使用RWF和CIKS模拟的H,SXA,H,S+XB,光电子能谱在振动构型基的选择中，最大总激励水平为8，同时激励模式的最大数量为3，单模最大激励水平为8；取近似对应于模态收缩Frank-Condon因子（MCFCFs）分布平均值的模态；根据初始状态的相对基态振动能级做出预筛选，110-3以下的初始状态将被忽略，RWF使用基于非齐次薛定调方程的时间无关本征态无ansatz 评估光谱，方程的解与光谱强度直接相关，在Lanczos基下RWF计算的收敛阅值设置为trwF=210-4；CIKS用Lanczos或RACE方法确定具有最大FranckC o n d o n 因子的本征态，通过RACE方法得到的光电子能谱如图8（b)所示，在该方法下，041005-5Harm2618.921200.832624.63VSCF2553.851169.702539.51VMRCI2518.371169.332520.21第41卷Modes2v3212v23v33v3v24v34v14v2V1+V2V2+V3V1+V3Vi+2v3V2+2v32v1+V3Vi+V2+V32v2+V3V2+212V2+V1Vi+3v3V2+3v12v1+2v3Vi+v2+2v32v2+2v33V+V3Vi+2v2+v32V1+V2+V33v2+V3V2+3v12v2+2vi3V2+V1“实验值来源于文献2 1-2 4定义了FCF之和的阈值定义为0.99.从图中容易看出，两种方法绘制光电子能谱图总体一致，但RWF方法得到的光谱图更为精细，较CIKS而言出现的峰值更多，在横坐标为11-11.5eV以及12.5eV后仍有的微弱峰值，如图8（a)中的放大部分所示图8(b)中12.2 eV至12.3eV处有两个峰值，而对应于(a)中的同等位置处有三个峰原子与分子物理学报表3H,S及H,S*的倍频、组合频（cm)Table 3Overtones and combination bands(cm-)of H,S and H,S+Sym.VMCSCFAi5335.13A5211.97A2381.21B27981.20A7870.59A13548.81A110654.61AI10430.76A4658.53A3884.59B23884.59B25344.37A17719.32A16489.46B27959.68B26301.26B25021.02A16373.77A4986.99B210647.79B29630.28A10157.4A7755.52A7652.14B210812.63B29074.86B27627.97B26169.2A9007.00A8677.58A6105.36第4期H2SH,S*VMRCIExp.5295.315240.85203.485145.12380.3723547867.377844.623541.9910656.7210356.834656.993851.113865.105203.457699.766452.777680.416289.84991.086362.934985.3610362.249001.179886.327871.287606.8210116.488817.807507.086123.538988.148542.166111.55041005-6VSCF5013.345061.792331.207422.677523.763484.059768.969939.844627.373779.23700.743789.13686.315147.365077.587535.296383.56136.787569.136201.454939.234825.116288.36186.074839.689928.138522.5410009.848636.387252.2610014.058670.167317.795955.578625.657302.705970.19值，由于峰值个数的增加，再加之图像绘制时高斯函数中半宽大小的设置，致使RWF中相应的峰略微右移，4结论本文基于MOLPRO量子化学程序包，使用从头算方法在(U)CCSD/cc-pVQZ 理论水平下计算VMCSCF5075.324962.662326.827595.67498.883463.910134.469944.274536.043821.663726.464797.27336.516202.897584.727120.224829.186094.354795.9810134.099189.827388.327373.238434.1210296.438696.07317.695942.158605.428303.644803.29VMRCI5039.584950.812326.097477.287472.063457.3310136.629863.714534.083664.683598.414801.217323.606168.037302.366093.434799.696166.544805.309855.848579.317321.227350.408259.019620.098417.487191.405903.228584.638343.594804.84483.第41卷0.50.40.20.101000 1500 200025003000 3500400045005000 5500图6 H,S的振动产生的红外光谱Fig.6IR spectrum of H,SRAMANSPECTRUM160F1401201008060402001000150020002500300035004000450050005500图7 H,S分子振动产生的拉曼光谱Fig.7Raman spectrum of H,SPHOTOELECTRONSPECTRUM0.60.50.40.30.20.1(b)0.50.40.30.20.1图 8 H,S XA,H,S+B,的光电子能谱Fig.8The XA,-XB,photoelectron spectrum of H,S了H,S 以及阳离子H,S+的势能面及光电子能谱.从势能面的特征中可以看出不同振动模式对分子能量变化的影响，S-H键拉伸振动产生的能量为郭甜，等：H,S和H,S*的振动光谱和电子能谱计算IR SPECTRUMWavenumber(cm-l)Wave number(cm-)0.0250.020.0150.010.0051010.110.210.310.410.510.610.712.510.51111.5第4期面内弯曲振动的2 至3倍，并且S-H键的非对称拉伸振动和面内弯曲振动共同作用会使得势能变化范围明显变大：VMCSCF和VMRCI分别被用来修正静态相关效应和动态相关效应以获取非谐振动频率和振动光谱，结果表明，较H,S而言，失电子使得H,S*在相同模式下的非谐振动频率变小此外，H,S分子在组合频v+Vs（52 0 3.45cm=）和倍频2 vi（52 0 3.48 c m=）处的带隙小于1cm1，产生强烈的费米共振，致使红外光谱中相应波段处的吸收强度明显变大，超过个别基频振动产生的吸收峰最后，基于势能面使用RWF和CIKS方法首次计算了H,S XA,H,S+X?B,的光电子能谱，两种方法所得电子能谱图像总体一致，但RWF方法所得光电子能谱图像显示出更多的峰值，为其实验研究和星际探测提供理论数据支撑参考文献：1Lodi L,Yurchenko S N,Tennyson J.The dipole mo-ment surface for hydrogen sulfide H,S J.J.Quant.Spectrosc.Ra.2015,161:41.2Botschwina P,Zilch A,Werner H J,et al.Ab initiocalculation of potential energy surfaces and spectro-scopic properties of H,S and H,S*J.J.Chem.Phys.,1986,85:5107.3 Kozin I,Jensen P.The potential energy surface and0.015the fourfold clustering of the rovibrational energy levels0.01in H,S molecule J.Mol.Spectrosc.,1993,163:0.00人1313.51212.5eV141313.54Cook R L,Lucia F,Helminger P.Molecular forcefield and structure of water:Recent microwave resultsJ.J.Mol.Spectrosc.,1974,53:62.5 Kauppi E,Halonen L.Potential energy surfaces for theground electronic state of water,hydrogen sulfide,andhydrogen selenide J.J.Phys.Chem.,1990,94:5779.6Ol.P,Tennyson J,Jensen P.The potential energysurface of H,0 J.J.Chem.Phys.,1996,105:6490.7Kozin I N,Jensen P.Fourfold clusters of rovibrational14energy levels for H,S studied with a potential energysurface derived from experiment J.J.Mol.Spec-trosc.,1994,163:483.8Martin J,Faracois J,Gijbels R.The high level ab ini-tio studies on the H,S molecules in the electronicground states J.J.Mol.Spectrosc Pronicial.,1995,169:445.041005-7第41卷9Xie J K,Yan G S.Theoretical study of the potentialenergy surface and the highly excited rovibrationalspectroscopy of H,S molecule J.Chinese Sci.Bull.,2001,46:650.10Cazzoli G,Puzzarini C.The rotational spectrum of hy-drogen sulfide:The H,S and H,S isotopologues re-visited J.J.Mol.Spectrosc.,2014,298:31.11 Carvaja M,Lemus R.Toward a global analysis methodbased on anharmonic ladder operators:application tohydrogen sulide(H,)J.J.Phys.Chem.A,2015,119:12823.12 Chong D P.Computational study of the electron spectraof vapor-phase indole and four azaindoles J.Mole-cules,2021,26:1947.13Chang J L,Kuo W H,Huang Y J,et al.Identificationof the photoelectron spectra of HFCS via computingFranck-Condon factors J.Comput.Theor.Chem.,2021,1204:113393.14 Yarzhemsky V G,Teterin Y A.Satellite excitationsand final state interactions in atomic photoionizationJ.Atoms,2022,10:73.15 Heislbetz S,Pfeiffer F,Rauhut G.Configuration se-lection within vibrational multiconfiguration self-con-sistent field theory:Application to bridged lithiumcompounds J J.C h e m.Ph y s.,2 0 11,134:204108.16 Meier P,Oschetzki D,Pfeiffer F,et al.Towards anautomated and efficient calculation of resonating vibra-tional states based on state-averaged multiconfigura-原子与分子物理学报tional approaches J.J.Chem.Phys.,2015,143:244111.17 Petrenko T,Rauhut G.A general approach for calcu-lating strongly anharmonic vibronic spectra with a highdensity of states:The X?B,X A,photoelectronspectrum of difluoromethane J.J.Chem.TheoryComput.,2017,13:551518Petrenko T,Rauhut G.A new efficient method for thecalculation of interior eigenpairs and its application tovibrational structure problems J.J.Chem.Phys.,2017,146:124101.19 Cook R L,Lucia F,Helminger P.Molecular forcefield and structure of hydrogen sulfide:recent micro-wave results J.J.Mol.Struct.,1975,28:237.20Neufeld D A,Godard B,Gerin,M.Sulphur-bearingmolecules in diffuse molecular clouds:new results fromSOFIA/GREAT and the IRAM 30 m telescope J.Astron.Astrophys.Rev.,2015,577:A49.21Bacic Z,Light J C.Theoretical methods for rovibra-tional states of floppy molecules J.Annu.Rev.Phys.Chem.,1989,40:469.22Allen H C,Plyler E K.Infrared spectrum of hydrogensulide.II.The 5100 cm region J.J.Chem.Phys.,1954,22:1104.23Hutson J M.Intermolecular forces from the spectrosco-py of van der Waals molecules J.Annu.Rev.Phys.Chem.,1990,41:123.24 Lane W C,Edwards T,Gilis J R,et