He-H2S复合物的六维从头算势能面和束缚态

采用[CCSD(T)]-F12a/aug-cc-pVTZ方法,同时在基组中引入中心键函数(3s3p2d1f1g)构建了He-H$_2$S复合物的高精度六维势能面. 除分子间振动坐标,同时考虑了H₂S分子内的v₁对称伸缩振动Q₁正则模、v₂弯曲振动Q₂正则模和v₃反对称伸缩振动Q₃正则模三种振动模式. 将计算得到的六维势能面在Q₁,Q₂和Q₃方向上分别做积分得到H₂S单体分别处于振动基态、v和v₃激发态下的He-H₂S的三个振动平均势能面. 计算结果表明,每个平均势能面都有一个T形全局极小值、一个平面局部极小值、两个平面内鞍点和一个平面外鞍点. 全局极小值的几何构型位于R=3.46 ?,θ=109.9°和φ=0.0°,势阱深度为35.301 ^-1. 在径向部分采用离散变量表象法和角度部分采用有限基组表象法并结合Lanczos循环算法计算了He-H₂S的振转能级和束缚态. 计算发现He-(para-H₂S)在H₂S的v₂和v₃区域的带心位移分别为0.025 cm^-1和0.031 cm^-1,而He-(ortho-H₂S)的带心位移分别为0.041 cm^-1和0.060 cm^-1,都表现为蓝移.     

He2分子 b3Πg(v=9)态的预解离

利用光外差-浓度调制吸收光谱技术测量了b³Π_g～a³∑_u⁺(9,3)带在12065～12445 cm^-1的光谱,并研究了He₂分子b³Π_g(v=9)态与c³∑_g⁺ 态的预解离相互作用．基于He₂分子c³∑_g⁺ 态的ab initio理论计算势能曲线和b³Π_g态的RKR势能曲线,分析了b³Π_g(v=9)态的预解离机制,并计算了b³Π_g(v=9)态转动能级的预解离线宽,计算结果同实验测量基本一致.     

CH4-Ne体系的三维从头算势能面与光谱的模拟

本文在从头算CCSD(T)/AVXZ(X=T,Q)水平上计算了CH₄-Ne体系的三维势能面,同时在基组中加入了键函数,并且将基组外推到完全基组极限水平. 通过最小二乘法拟合摩斯长程势能函数形式,获得了三维分析势能面,其中对664个从头算格点的拟合的标准偏差仅为0.042 cm^-1. 随后,采用径向离散变量表象和角度有限基组表象相结合的杂化基函数方法,并基于Lanzcos迭代的方法获得了CH₄-Ne体系的振转能级,所预测的红外光谱与实验值非常吻合. 本文首次给出了CH₄-Ne体系的微波光谱数据. CH₄-Ne体系三维摩斯长程势能面分析的表达式为今后CH₄在Ne团簇动力学性质以及其碰撞诱导吸收光谱的研究中奠定了基础.     

O3+与H2碰撞中非解离电荷转移过程的全量子计算

应用全量子的分子轨道强耦合方法，研究了基态的O³⁺(2s²2p ²P)与氢分子碰撞的非解离电荷转移过程，计算了不同方位角(25°，45°，89°)，能量分别为100，500，1000和5000eV/u时的单电子俘获的振动分辨的态选择截面及相应的微分截面.分子轨道强耦合计算中采用了自旋耦合价带理论计算的三原子分子势能面和径向耦合矩阵元.对氢分子的自身转动，采用无限阶的冲量近似方法；对体系的电子运动同H₂或H^+<     

基于激发态运动方程耦合簇的神经网络透热势能矩阵：苯硫酚1πσ*态-介导的光解

本文用最近发展的神经网络拟合方法[Chin. J. Chem. Phys. 34,825 (2021)]构造了一个新的涉及苯硫酚¹πσ^*态-介导光解的¹ππ^*和¹πσ^*态耦合非绝热势能面. 势能面包含了解离过程中的三个关键振动模式,即S-H伸缩、弯曲和扭转振动. 由于单双激发态运动方程耦合簇方法具有简单、效率高、精度高的优点,采用激发态运动方程耦合簇方法计算了苯硫酚激发态¹ππ^*和¹πσ^*的绝热能量. 神经网络拟合绝热激发态S₁和S₂态的均方根误差分别为0.89和1.33 meV,表明神经网络方法具有很高的拟合精度. 在构建非绝热势能面的过程中,仅利用了体系绝热势能,避免了非常耗时的非绝热耦合计算,极大地提高了效率. 为了检测新的非绝热势能面的可靠性,本文进一步展开了苯硫酚光解非绝热动力学模拟. 动力学计算得到的S₁振电态0⁰和3¹的寿命与实验和之前的理论结果均吻合,验证了基于激发态运动方程耦合簇绝热能量构建的非绝热势能面的可靠精确性,并可进一步应用到实际大分子体系中.     

反应N(2D)+H2→NH+H和N(2D)+D2→ND+D的立体动力学理论研究.

基于Ho等人的精确势能面(J. Chem. Phys. 119, 3063(2003))研究,运用准经典轨线方法计算了21.3 kJ/mol碰撞能下反应N(²D)+H₂→NH+H和N(²D)+D₂→ND+D的产物与反应物之间的矢量相关．发现两个反应的产物角分布都是前向和后向呈现峰值分布,产物的转动角动量矢量j′不仅是取向的,而且是在y轴负方向上定向的．两个反应显示出的同位素效应主要归因于同位素质量的差别．     

用神经网络方法拟合的反应体系H+CH4→H2+CH3的一个全域势能 (cited:12)

利用神经网络方法,基于47783个高精度从头算能量点构建了反应体系H+CH₄→H₂+CH₃的一个全域势能面．通过大量的准经典轨线以及量子散射计算测试了势能面的收敛性质．这个势能面对于拟合过程以及从头算点的数目都已经完全收敛,拟合误差很小且比Shepard插值的势能面计算速度更快,代表了此标志性多原子反应体系最好的势能面．     

O(1D)和CD4反应的准经典轨线

在新的全域势能面上, 用准经典轨线方法细致地研究了O(¹D)+CD₄多通道化学反应的动力学．这个势能面是用交换不变多项式方法基于MRC+Q/aug-cc-pVTZ从头算点拟合得到的．通过计算得到了产物OD+CD₃、D+CD₂OD/CD₃O和D₂+DCOD/D₂CO的分支比、平动能分布以及角度分布,结果显示理论与实验吻合得较好, 从而说明了这个反应的同位素取代效应很小． 研究表明,O(¹D)+CD₄反应是经过陷入的抽取机理发生的： 最初主要通过D原子的抽取,并不是之前人们认为的直接C-D键的插入形成CD3OD中间物后再进而解离成各个产物通道．     

反应物NO的转动激发对反应N(4S)+NO(X2Π)→N2(X3Σg-)+O(3P)影响的立体动力学研究

采用准经典轨线方法,在碰撞能为0.6 eV时,研究了反应物NO分子的转动激发对发生在³A"和³A'势能面上的反应N(⁴S)+NO(X²Π)→N₂(X³Σ_g^-)+O(³P)的立体动力学性质的影响. 详细讨论了在反应物分子的不同转动态下发生在两个势能面上反应的矢量性质. 结果表明, 反应物分子NO的转动激发对发生在³A′势能面上的立体动力学性质产生重大影响, 这可能与该势能面上存在一个早期势垒有关.     

BF自由基X1∑+和a3∏态光谱常数和分子常数研究

采用内收缩多参考组态相互作用方法和相关一致基aug-cc-pV6Z, 对BF自由基X¹∑⁺和a³∏ 态的势能曲线进行了研究. 计算是在0.095---1.33 nm的核间距内进行的. 为获得更准确的结果, 计算中还考虑了Davidson修正、相对论修正及核价相关修正对势能曲线的影响. 相对论修正采用的方法是二阶DouglasKroll哈密顿近似, 修正计算是在cc-pV5Z基组水平上进行的. 核价相关修正使用的是cc-pCV5Z基组. 利用得到的势能曲线, 拟合出了各种修正下BF自由基X¹∑⁺和a³∏ 态的光谱常数D_e, R_e, ω_e, ω_ex_e, ω_ey_e, B_e和α_e、并与实验结果进行了比较. 结果表明: 考虑Davidson修正、相对论修正和核价相关修正后得到的光谱常数最接近实验结果. 利用修正后的势能曲线, 通过求解径向振转Schrödinger方程, 找到了转动量子数J = 0时这两个电子态的全部振动态, 并计算了每一电子态前20个振动态的振动能级、惯性转动常数和离心畸变常数, 其值与已有的实验结果较为一致. 本文得到的光谱常数和分子常数达到了很高的精度, 能为进一步的光谱实验提供可靠的参考.     

Microwave spectrum of the hydrogen sulfide molecule H232S in the ground state

Frequencies, line strengths, and intensities of microwave electric dipole rotational transitions of H₂³²S in the ground state have been calculated. The calculation was based on representation of the effective rotational Hamiltonian operator in the form of a Pade operator. As initial information, all known microwave data on frequencies of rotational transitions of H₂³²S in the ground state were used. Some of these data were obtained in the present paper.     

Analysis of ν2 of H233S and H234S

We have extended our analysis of the (010) vibrational state of H₂S, this time using Watson's A-reduced Hamiltonian (through P⁸ terms) in the I′ rotational representation. We have determined separate sets of (010) upper state constants for each isotopomer (H₂³²S, H₂³³S, and H₂³⁴S) by fitting the ν₂ spectral lines, keeping the ground state constants fixed to the values determined by Flaud, Camy-Peyret, and Johns. Determinable coefficients for H₂³²S and a slightly revised set of ν₂ line assignments for H₂³³S and H₂³⁴S are also reported.     

Analysis of ν2 of D2S

The infrared spectrum of ν₂ of D₂S was recorded from 740 to 1100 cm^?1 on the University of Denver 50-cm FTIR spectrometer system. We have assigned 655 transitions from D₂³²S and 129 from D₂³⁴S, and have analyzed them using Watson's A-reduced Hamiltonian evaluated in the I^r representation. We used the recently published D₂³²S and D₂³⁴S ground state Hamiltonian constants [C. Camy-Peyret, J. M. Flaud, L. Lechuga-Fossat and J. W. C. Johns, J. Mol. Spectrosc.109, 300–333 (1985)]. Upper state Hamiltonian constants were obtained from a fit of the ν₂ transitions, keeping the ground state constants fixed while varying the upper state constants. The standard deviation of the D₂³²S ν₂ fit is 0.0025 cm^?1. The standard deviation of the D₂³⁴S ν₂ fit is 0.0041 cm^?1.     

Measurement and analysis of the ν2 band of H2S: Comparison among several reduced forms of the rotational Hamiltonian

The frequencies of 123 lines between 1080 and 1260 cm^?1 in the ν₂ band of H₂S were measured with a tunable diode laser and a grating spectrometer. Ground- and upper-state molecular constants appropriate to Watson's A reduced Hamiltonian in the I^r coordinate representation were obtained from a simultaneous fit of the ν₂ band data and previously measured H₂S rotational microwave frequencies. This fit had a standard deviation comparable to the experimental uncertainties. The A reduced Hamiltonian in the III^r representation, however, was found to fit the data very poorly. The large value of the transformation parameter s₁₁₁ calculated for the A reduced Hamiltonian in the III^r representation confirms that this result is due to the slow convergence of this particular reduced Hamiltonian.     

Centrifugal distortion analysis of pure rotational spectra of H216O,H217O,and H218O

Centrifugal distortion analyses based entirely on high-quality infrared data are carried out for the ground vibrational states of H₂¹⁶O, H₂¹⁷O, and H₂¹⁸O. As a result of the analyses, the values of 27 rotation and distortion constants for each species are determined. By using these constants it was possible to improve considerably the accuracy of the literature values for rotational energy levels at high J_τ, especially for H₂¹⁷O and H₂¹⁸O. The experimental values for the energy levels are deduced from the observed rotational transitions constituting the fitted data.     

The 923-cm−1 band of CF235Cl2 (freon 12), studied by infrared-microwave double resonance

The constants of the 923-cm^?1 band (ν₆) of CF₂³⁵Cl₂ were accurately determined by the application of infrared-microwave double resonance using CO₂ and N₂O lasers. The frequencies of 32 ground-state and 34 (v₆ = 1) state rotational transitions and 14 infrared transitions were measured. The infrared transition frequencies were generally obtained with an accuracy of ± 20 MHz, but in some cases observation of two-photon Lamb-dips allowed the accuracy to be improved to ± 5 MHz. Many double-resonance signals displayed a predominantly “dispersion-type” lineshape and it has been shown that the phase of the observed signals gives information over the relative disposition of the energy levels involved.     

Analysis of ν2 of H2S

The infrared absorption spectrum of ν₂ of H₂S in the region from 1000 to 1500 cm^?1 was obtained with a resolution limit of <0.05 cm^?1 on the University of Denver 50-cm FTIR spectrometer system. We have assigned 387 lines due to H₂³²S, 75 lines due to H₂³⁴S, and 15 lines due to H₂³³S, and have analyzed them using Typke's reduction of Watson's Hamiltonian. Slightly revised ground-state constants for the 32 isotope were obtained from a simultaneous fit of the microwave transitions observed by Helminger, Cook, and De Lucia, combined with weighted averaged ground-state combination differences formed from the infrared bands (010), (020), (100), (001), (110), (011), (210), and (111). The standard deviation for the fit was 0.0018 cm^?1 for the infrared data and 0.000032 cm^?1 for the microwave lines. Upper-state constants for the 32 isotope were obtained from a least-squares fit of the spectral lines of ν₂, keeping the ground-state constants fixed to the values determined by the combination difference fit. The standard deviation of the (010) line fit was 0.0017 cm^?1 for the 32 isotope. Ground-state and upper-state isotopic mass adjustment constants were determined in a simultaneous fit of lines of H₂³³S and H₂³⁴S, keeping the ground-state and upper-state constants for the 32 isotope fixed.     

Linewidth and intensity measurements of SO2 lines at 94 GHz with self-broadening and broadening by H2O and N2

In the frequency range between 91.5 and 95.5 GHz, three rotational lines of the ³²S¹⁶O₂ and two rotational lines of the ³⁴S¹⁶O₂ molecules in the fundamental vibrational state, and also two lines of the ³²S¹⁶O₂ molecule in the v₂ vibrational state, have been investigated. Center frequencies and absolute absorptions have been measured and compared with theoretical values. Furthermore, the self-broadening and broadening by H₂O and N₂ of the transition 23(6,18)–24(5,19) with the line center at 94.064 GHz have been investigated. The following linewidth parameters were found: SO₂-SO₂, 18.2±0.3 MHz/torr; SO₂-N₂, 3.8±0.1 MHz/torr; SO₂-H₂O, 15.2±0.2 MHz/torr. The bridge spectrometer and the measuring method used are also described.     

The infrared spectra of 12C32S, 12C34S, 13C32S,and 12C33S

Using a tunable diode laser spectrometer, the infrared absorption spectra of four isotopic species of carbon monosulfide have been observed in the positive column of a dc discharge of CS₂ and Ar. The wavenumbers of 115 vibration-rotation transitions between 1180.5 and 1266.1 cm^?1 have been measured. These lines were assigned to the 1-0, 2-1, 3-2, and 4-3 bands of ¹²C³²S, the 1-0 and 2-1 bands of ¹²C³⁴S and ¹³C³²S, and the 1-0 band of ¹²C³³S. These new data have been combined with the previous infrared and microwave results to determine Dunham coefficients (Y_ij), the Dunham potential expansion constants (a₀,a₁,a₂,a₃, and a₄), and the classical turning points by the RKR method.     

Dynamic properties of P4O6S and P4O7: P spin–echo and P MAS–NMR investigations

The rotational dynamics of P₄O₆S and P₄O₇ in the solid state were studied by means of ³¹P NMR spectra of spinning and static powder samples in the temperature range of 153–295 K and 295–388 K, respectively. All spectra were simulated to confirm the type of the motion and to extract the time scales as a function of the temperature. Good agreement between experimental and theoretical data was obtained on the basis of a three-site jump model. For P₄O₆S, the activation energy and the pre-exponential factor derived from the lineshape simulations amount to 51(2) kJ/mol and 6(3)·10¹⁵ s⁻¹. For P₄O₇, the spectral analysis yields an activation energy of 67(1) kJ/mol and a pre-exponential factor of 6(2)·10¹⁴ s⁻¹. The dynamic behavior was checked independently by lineshape analyses under both MAS and static conditions. Activation energies are consistent within the errors for the lineshape analyses. Additionally, we have analyzed spin–lattice relaxation measurements, which show the correct trends for the activation energies.