Department of Applied Biology and Chemical Technology, The Hong Kong Polytechnic University, Hung Hom, Hong Kong.
J Chem Phys. 2011 Sep 28;135(12):124312. doi: 10.1063/1.3640037.
Geometry optimization and harmonic vibrational frequency calculations have been carried out on the X̃(2)A(') state of P(2)H and the X̃(1)A(') state of P(2)H(-) using the restricted-spin coupled-cluster single-double plus perturbative triple excitation [RCCSD(T)] and explicitly correlated unrestricted-spin coupled-cluster single-double plus perturbative triple excitation [UCCSD(T)-F12x] methods. For RCCSD(T) calculations, basis sets of up to the augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality were employed, and contributions from extrapolation to the complete basis set limit and from core correlation of the P 2s(2)2p(6) electrons were also included. For UCCSD(T)-F12x calculations, different atomic orbital basis sets of triple-zeta quality with different associated complementary auxiliary basis sets and different geminal Slater exponents were used. When the P 2s(2)2p(6) core electrons were correlated in these F12x calculations, appropriate core-valence basis sets were employed. In addition, potential energy functions (PEFs) of the X̃(2)A(') state of P(2)H and the X̃(1)A(') state of P(2)H(-) were computed at different RCCSD(T) and UCCSD(T)-F12x levels, and were used in variational calculations of anharmonic vibrational wavefunctions, which were then utilized to calculate Franck-Condon factors (FCFs) between these two states, employing a method which includes allowance for anharmonicity and Duschinsky rotation. The photodetachment spectrum of P(2)H(-) was then simulated using the computed FCFs. Simulated spectra obtained using the RCCSD(T)/aug-cc-pV5Z and UCCSD(T)-F12x(x = a or b)/aug-cc-pCVTZ PEFs are compared and found to be essentially identical. Based on the computed FCFs, a more detailed assignment of the observed vibrational structure than previously reported, which includes "hot bands," has been proposed. Comparison between simulated and available experimental spectra has been made, and the currently most reliable sets of equilibrium geometrical parameters for P(2)H and its anion have been derived. The photodetachment spectrum of P(2)D, yet to be recorded, has also been simulated.
已经使用受限自旋耦合簇单双加微扰三激发 [RCCSD(T)] 和显式关联无限制自旋耦合簇单双加微扰三激发 [UCCSD(T)-F12x] 方法对 P(2)H 的 X̃(2)A(') 态和 P(2)H(-) 的 X̃(1)A(') 态进行了几何优化和简谐振动频率计算。对于 RCCSD(T) 计算,使用了多达扩充相关极化价五重zeta(aug-cc-pV5Z)质量的基组,并包括了对完全基组极限的外推和 P 2s(2)2p(6) 电子的核心相关的贡献。对于 UCCSD(T)-F12x 计算,使用了不同质量的三重 zeta 原子轨道基组,以及不同的关联辅助基组和不同的双子 Slater 指数。当在这些 F12x 计算中相关 P 2s(2)2p(6) 核心电子时,使用了适当的核价基组。此外,还在不同的 RCCSD(T) 和 UCCSD(T)-F12x 水平上计算了 P(2)H 的 X̃(2)A(') 态和 P(2)H(-) 的 X̃(1)A(') 态的势能函数(PEF),并将其用于简谐振动波函数的变分计算,然后利用该方法计算这两个态之间的 Franck-Condon 因子(FCF),其中包括允许非谐性和 Duschinsky 旋转。然后使用计算出的 FCF 模拟了 P(2)H(-) 的光电子谱。使用 RCCSD(T)/aug-cc-pV5Z 和 UCCSD(T)-F12x(x = a 或 b)/aug-cc-pCVTZ PEF 获得的模拟光谱进行了比较,发现它们基本相同。基于计算出的 FCF,提出了比以前报道的更详细的观察到的振动结构的分配,其中包括“热带”。已经进行了模拟光谱与现有实验光谱之间的比较,并推导出了目前最可靠的 P(2)H 和其阴离子的平衡几何参数集。尚未记录的 P(2)D 的光电子谱也已被模拟。