迭代的多组态相互作用方法（英文）

本文提出了迭代的组态相互作用方法(IMRCI).IMRCI被用来计算H_2O和CH_2(单重态和三重态)在平衡态和远离平衡时的电子能量.IMRCI、MP2、MP3、MP4、CCSD和CCSD(T)还用来计算H_2O、CH_2(单重态和双重态)和N_2的势能曲线.这些计算结果表明IMRCI的结果不依赖初始的多参考态组态函数,并能较快地收敛到完全组态相互作用的计算结果.相比完全组态相互作用的结果,仅需2至4次迭代,IMRCI的误差就能达到10~(-5)hartree数量级.另外,IMRCI还提供了寻找势能面上主要电子组态的一个有效途径.这将有助于单参考态和多参考态理论模型得到准确的计算结果.     

BeCl分子电子激发态的多参考组态相互作用计算

运用含Davidson修正的多参考组态相互作用方法,在aug-cc-pVTZ基组水平上,对BeCl分子基态和相同多重度的几个低电子激发态进行了势能扫描计算.通过群论原理确定各电子态对称性及离解极限.将其中基态(X²Σ⁺)和第一激发态(A²Π})对应的势能曲线拟合到Murrell-Sorbie解析势能函数形式,得到基态(X²Σ⁺)的离解能及主要光谱常数(括号中为文献[6]提供的实验值)为D_e=3.74eV,R_e=0.18173nm(0.17970),w_e=857.4cm¹(847.2),w_ex_e=5.03cm^-1(5.14),B_e=0.7103cm^-1(0.7285),α_e=0.0059cm^-1(0.0069),第一激发态(A²Π)的D_e=3.02eV,R_e=0.18369nm(0.18211),w_e=832.7cm^-1(822.1),w_ex_e=5.93cm^-1(5.24),B_e=0.6953cm^-1(0.7094),α_e=0.0065cm^-1(0.0068),计算结果与实验值符合得较好.另外,通过Level程序求解双原子径向核运动的Schrödinger方程得到J=0时BeCl分子这两个电子态的全部振动能级.     

F−(H2O)+CH3I配体交换反应动力学

本文介绍F^?(H₂O)+CH₃I→[FCH₃I]^?+H₂O在交叉分子束碰撞能量0.3∽2.6 eV的配体交换动力学成像结果. 产物的动能受到弱键结合配合物的稳定性的影响,大量水分子的内部激发不利于中间物有效的能量重新分配,随着碰撞增加,低动能受到抑制. 在0.3 eV时,内部亲核取代非常重要,为形成I^?和I^?(H₂O)的竞争性亲核取代途径提供了依据.     

OH+C2H3F的反应机理及产物支化比反应的理论研究

本文对HOC₂H₃F可能解离通道的势能面进行从头算CCSD(T)/CBS//B3LYP/6-311G(d,p)计算,同时对速率常数进行Rice-Ramsperger-Kassel-Marcus计算. 生成主要产物CH₂CHO+HF最有利的反应途径是OHC₂H₃F→i2→TS14→i6→TS9→i3→TS3→CH₂CHO+HF,其中速率决定步骤是HF通过TS11从CO桥接位置解离,能量比反应物高3.8 kcal/mol. 借助中间态TS14,F原子从C_α迁移到C_β位置生成CH₂O+CH₂F,然后通过中间态TS16,H从O迁移到C_α位置;通过中间态TS5,C-C键断裂生成产物,其能量比反应物低1.8 kcal/mol,比TS11低4.0 kcal/mol.     

O-与CH3CN反应机理的密度泛函理论研究

采用G3MP2B3方法研究了氧负离子与乙腈反应的势能剖面.在(U)B3LYP/6-31+G(d,p)理论水平下分别优化了该反应势能面上反应物、产物、中间体和过渡态的分子结构,采用G3MP2B3方法校正了这些关键点的能量. 势能面上的各个反应路径均通过针对过渡态的內禀反应坐标理论计算加以确定. 分别考察了四个可能的热力学产物通道,即质子转移、氢原子转移、H₂⁺转移和双分子亲核取代反应途径. 其中,经H₂⁺转移生成H₂O的反应通道为该反应的主要产物通道.     

CH3OH在ZnO(0001)表面的光催化解离

本文利用266 nm波长的激光及程序升温脱附的方法研究了甲醇在ZnO(0001)表面的光催化反应. TPD结果显示部分的CH₃OH以分子的形式吸附在ZnO(0001)表面,而另外一部分在表面发生了解离. 实验过程中探测到H₂,CH₃^·,H₂O,CO,CH₂O,CO₂和CH₃OH这些热反应产物. 紫外激光照射实验结果表明光照可以促进CH₃OH/CH₃O^·解离形成CH₂O,在程序升温或光照的过程中它又可以转变为HCOO^-. CH₂OH_Zn与OH_ad反应在Zn位点上形成H₂O分子. 升温或光照都能促进CH₃O^·转变为CH₃^·. 该研究对CH₃OH在ZnO(0001)表面的光催化反应机理提供了一个新的见解.     

对典型多通道反应H+CH3OH的模式特异动力学研究

H+CH₃OH作为典型的多通道反应,在燃烧和星际中起着重要的作用. 本文基于在UCCSD(T)-F12a/AVTZ水平上计算的大量数据点,构建了该体系的全维精确势能面,并基于该势能面,研究了不同产物通道的模式特异动力学. 结果表明,O-H 伸缩、沿C-O轴的扭转以及C$-$H伸缩等模式的振动激发对H₂+CH₃O、H₂+CH₂OH、H₂O+CH₃和H+CH₃OH四个产物通道有着不同的影响. 该研究有助于理解具有多个产物通道的复杂反应的模式特异动力学,进而帮助控制其竞争反应.     

类镍等电子系列离子M1,M2,E2禁戒跃迁特性的理论研究

利用基于全相对论框架下的多组态Dirac-Fock理论方法发展起来的程序包GRASP92和新发展的处理辐射跃迁过程的程序REOS99，计算了类镍等电子系列离子(Z=45—95)的基组态3s²3p⁶3d¹⁰¹S₀以及低激发组态3s²3p⁶3d⁹4l,3s²3p⁵3d¹⁰4l和3s3p⁶3d¹⁰4l(l=s,p,d,f)的能级及其向基态的M1,M2,E2禁戒跃迁概率.通过分析高离化类镍离子在特定的原子序数范围内由于存在能级交叉而产生的强组态相互作用，解释了高离化类镍离子禁戒跃迁概率的反常变化现象，探讨了禁戒跃迁概率受强组态相互作用影响而变化的一般规律.     

硫酸催化乙二醛气体相水化反应机理和速率常数理论计算研究

采用M06-2X和CCSD(T)高阶量化计算和传统过渡态理论研究硫酸催化乙二醛气体相水化反应.对HCOCHO+H₂O, HCOCHO+H₂O+H₂O, HCOCHO+H₂O+H₂O, HCOCHO+H₂O...H₂SO₄和HCOCHO+H₂O+H₂SO₄五个路径的反应机理和速率常数进行了研究.计算结果表明硫酸具有较强的催化能力,能显著减小乙二醛水化反应的能垒,在CCSD(T)/6-311++G(3df,3pd)//M06-2X/6-311++G(3df,3pd)理论水平,当硫酸分子参与乙二醛水化反应时,反应能垒从37.15 kcal/mol减少至7.08 kcal/mol.在室温条件下,硫酸催化乙二醛水化反应的反应速率1.34×10^-11 cm³/(molecule.s),是等量水分子参与乙二醛水化反应的速率的10¹²倍,大于乙二醛与OH自由基反应的反应速率1.10×10^-11 cm³/(molecule.s).这表明大气条件下,硫酸催化乙二醛水化反应可以发生,同乙二醛与OH自由基反应相竞争.     

New ab initio Potential Energy Surfaces for Cl(2P3/2, 2P1/2)+H2 Reaction

采用多参考组态相互作用的计算方法和很大的基组构造了Cl+H₂反应体系的一组新的三维从头算势能面. 该势能面包含Cl+H₂体系的能量最低的三个绝热态,并在非绝热近似下转化为四个非绝热势能面. 另外,旋轨耦合矩阵元也基于Breit-Pauli Hamil-tonian计算得到. 对角化这四个非绝热势能和两个旋轨耦合矩阵元组成的全耦合Hamiltonian,得到了三个考虑旋轨耦合后的绝热势能面.基于这组新势能面的非绝热动力学计算结果与最新的实验符合得很好,很好地解释了     

Highly accurate coupled-cluster singlet and triplet pair energies from explicitly correlated calculations in comparison with extrapolation techniques

Coupled-cluster singles and doubles (CCSD) valence shell correlation energies of the systems CH₂ (<¹A₁ state), H₂O, HF, N₂, CO, Ne, and F₂ are computed by means of standard calculations with correlation-consistent basis sets of the type cc-pVxZ (x = D, T, Q, 5, 6) and by means of explicitly correlated coupled-cluster calculations (CCSD-R12/B) with large uncontracted basis sets of the type 19s14p8d6f4g3h for C, N, O, F, and Ne and 9s6p4d3f for H. These CCSD-R12/B calculations provide reference values for the basis set limit of CCSD theory. The computed correlation energies are decomposed into singlet and triplet pair energies. It is established that the singlet pair energies converge as X^?3 and the triplet pair energies as X^?5 with the cardinal number of the correlation-consistent basis sets, and an extrapolation technique is proposed that takes into account their different convergence behaviour. Applied to the cc-pV5Z and cc-pV6Z results, this new extrapolation yields pair energies with a mean absolute deviation of 0.02 mE_h from the CCSD-R12/B reference values. For the seven systems under study, the extrapolated total valence shell correlation energies agree to within 0.2 mE_h with the CCSD-R12/B benchmark data.     

Photoaddition of water to uracil: A model study of loss of absorbance by pyrimidine adducts using geometry optimization of ground and excited states

Geometries of the ground, first triplet and first singlet excited states of a hydrohydroxide of uracil have been optimized using the SINDO 1 SCF molecular orbital method in conjunction with configuration interaction. The first triplet and first singlet transitions of the molecule are found tc/be of (π-sp*) type which shows their orbitally forbidden character. Thus an explanation as to why the adducts formed by the addition of H₂O, hydroxylamine, bisulphite etc at the CC double bonds of pyrimidines do not absorb in the characteristic 260 nm region has been provided.     

Multireference second-order Brillouin–Wigner perturbation theory§

The multireference, state-specific, second-order, Brillouin–Wigner perturbation theory (MR-BWPT2) is presented. A posteriori corrections are made which, in the case of a single reference function, recover the well-known formula of second-order many-body perturbation theory, i.e. Møller–Plesset ‘MP2’ theory, and in the multireference case suppress terms which scale non-linearly with the number of electrons in the system. Prototype calculations are reported for the (H₂)₄ model in which the degree of quasi-degeneracy can be varied by changing a single geometric parameter. The calculated total energies obtained by a second-order Brillouin–Wigner treatment are compared with those supported by CAS-MP2 (complete active space self-consistent field followed by second-order Møller–Plesset-like treatment of dynamic correlation effects), by MR-BWCC (multireference Brillouin–Wigner coupled cluster expansion), and by full configuration interaction. MR-BWPT2 provides a theoretical apparatus comparable to the widely used MP2 theory, but which can be applied to quantum chemical problems requiring a multireference formalism.     

Isomers of OCS and their reaction with H2O on singlet potential energy surface

The isomers of the carbonyl sulfide (OCS) molecule are investigated in detail at CCSD(T)/cc-pVTZ//MP2/6-311++G(2d,2p) level of theory. One cyclic isomer was identified along with three different linear minima of the OCS molecule. Three interconversion transition states were also located between cyclic and linear forms of OCS. Among these four isomers, the singlet potential energy surface (PES) for the molecule–molecule reaction between the three most energetically favoured isomers of OCS and H₂O has been explored theoretically at the CCSD(T)/cc-pVTZ//MP2/6-311++G(2d,2p) level. This singlet PES comprises of three paths. Path 1 is the reaction of linear OCS molecule with water producing the major product P1 (CO₂?+?H₂S), minor product P2 (S?+?HCOOH) and two isomers via 14 minima and 15 transition states. The Path 2 is an isomerization process in which cyclic isomer of OCS reacts with water molecule via another initial barrierless aduct producing five isomers of the OCS–H₂O system through five interconversion transition states. The reaction of linear COS isomer with water is shown in Path 3. This path produces the radicals SH and COOH from another COS–H₂O complex via a transition state. Among these three products, the product P1 is energetically most favoured. The overall exothermicity of the product channels for the formation of major product P1 on PES is calculated to be about 10.60?kcal/mol possessing initial high entrance barriers of 45.48 and 55.47?kcal/mol in two possible pathways. As the process is favoured thermodynamically but not kinetically, the reaction is expected to be very slow.     

Addition complex and insertion isomers on the potential energy surface of the reaction of indium dimer with water studied with relativistic ECP

Stationary points on the lowest singlet and triplet In₂ + H₂O potential energy surfaces (PESs) have been explored using the coupled cluster method, including single and double excitations with perturbative triples (CCSD(T)), and the density functional theory (DFT), employing the effective core potential (ECP) for indium (In), which accounts for scalar relativistic effects, with the triple-zeta quality basis set. The CCSD(T) calculated binding energy and anharmonic ν₂-bending mode frequency for the triplet ground-state addition complex, In₂^… OH₂(³B₁), are consistent with the complex detected in the matrix isolation infrared (IR) spectroscopic study under the thermal conditions. The two minimum energy crossing points between the triplet and the singlet PESs that have been located between the structures of In₂^…OH₂ and the transition state for the O–H bond breakage are not likely to be thermally accessible under the low-temperature matrix conditions. With the CCSD(T)-calculated In₂ + H₂O reaction profile and anharmonic vibrational frequencies for several In₂(H)(OH) insertion product isomers, we support the IR matrix isolation detection (by two experimental groups) of the lowest energy singlet double-bridged In(μ-H)(μ-OH)In isomer. For the proposed two-step mechanism of H₂ elimination from the In₂(H)(OH) species, the estimated energy barriers are also compatible with experiment.     

Theoretical investigation of the oxidation pathways of the Cl-initiated reaction of 2-methyl-3-buten-2-ol

The mechanism and products of the reaction of 2-methyl-3-buten-2-ol (MBO232) with Cl atoms in the presence of O₂ have been elucidated by performing high-level quantum chemistry calculations. The geometries of the reactants, intermediates, transition states, and products are optimized at the MP2(full)/6-311G(d,?p) level, and their single-point energies are refined at the CCSD(T)/6-311?+?G(d,?p) level. The potential energy surface profiles have been constructed at the CCSD(T)/6-311?+?G(d,?p)//MP2(full)/6-311G(d,?p)?+?0.95?×?ZPE level of theory, and the possible channels involved in the reaction are also discussed. The calculations indicate that the reaction predominantly proceeds via the addition of Cl atoms to the double bond rather than the direct abstraction of the H atoms in MBO232. The nascent adducts (CH₃)₂C(OH)CHCH₂Cl (IM1) and (CH₃)₂C(OH)CHClCH₂ (IM2) do not undergo subsequent isomerization and dissociation reactions, but rather react with O₂. The theoretical results show that the major products are CH₂ClCHO and CH₃C(O)CH₃ for the reaction of MBO232?+?Cl in the presence of O₂, which is in good agreement with the experimental finding.     

Electron paramagnetic resonance of quasi-one-dimensional crystal [Nd2(Cl3COO)6(H2O)3]n · nH2O

A quasi-one-dimensional single crystal of [Nd₂(Cl₃COO)₆(H₂O)₃]_n · nH₂O in which chains are built up of two alternating neodymium-ion dimer fragments is studied by the EPR technique. It is found that anisotropic interactions between neodymium ions in a chain are responsible for the complex shape of the EPR spectrum. Two groups of EPR signals are distinguished in the spectrum. Each group corresponds to one of the chain dimers disturbed by the interaction with neighbors in the chain. The shape of the EPR spectra is interpreted as a superposition of the spectra of chain fragments which have different lengths and are formed by the alternating magnetic triplet and nonmagnetic singlet states in the chain. Consideration is given to two cases when two alternating dimer fragments are either equivalent or nonequivalent to each other. It is shown that the spectral shape is primarily determined by the superposition of the spectra of an isolated triplet state (S=1), two interacting triplet states, and three interacting triplet states whose weighting contributions differ for the above two cases. The tensors of the anisotropic spin-spin interaction are determined, and the contribution from the isotropic component of the interaction is estimated.     

In search of the dark state of 5-methyl-2-hydroxypyrimidine using a numerical DFT/MRCI gradient

In a recent publication, Lobsiger et al. [Phys. Chem. Chem. Phys. 12, 5032 (2010)] presented infrared and electronic absorption spectra of supersonic jet-cooled 5-methyl-2-hydroxypyrimidine (5M2HP), the enol form of deoxythymine. In addition, they reported on the fast nonradiative decay of the S₁ population to a dark state. In the present paper, we have investigated the mechanism and rate constants of this nonradiative decay by means of quantum chemical multi-configuration methods. To this end, minima of the lowest excited singlet and triplet states as well as the minimum-energy crossing point of singlet and triplet potential energy hypersurfaces (PEHs) have been determined employing a numerical DFT/MRCI gradient where DFT/MRCI stands for a combination of density functional theory (DFT) and a semi-empirical multi-reference configuration interaction (MRCI) approach. Rate constants have been calculated in the Condon approximation using a time-dependent approach based on harmonic oscillator functions and electronic spin–orbit coupling matrix elements evaluated at the DFT/MRCI level. It is shown that the first excited triplet state possesses ³(n?→?π*) character in the gas phase. Fast intersystem crossing is mediated by the low-lying ³(π?→?π*) state whose PEH crosses both, the S₁ ¹(n?→?π*) and T₁ ³(n?→?π*) PEHs.     

Size-consistent multireference configuration interaction method through the dressing of the norm of determinants

A new determinant-specific, effective change (‘dressing’) of the norm of the multireference configuration interaction (MRCI) wavefunction is proposed in order to achieve the size-consistency of the MRCI method. The new approach provides a unifying framework for analysis of size-consistent extensions of the MRCI method that are based on the coupled pair functional (CPF) strategy and lead to simplified computations of the analytical gradients. Using the new framework, a generalized multireference full coupled pair functional (MR-FCPF) method is introduced. The MR-FCPF method may be viewed as a functional counterpart of the recently proposed generalized (‘full’) coupled electron pair approximation (CEPA), referred to as the size-consistent self-consistent CI ((SC)²CI) method. A straightforward extension of the MR-FCPF method leads to a pseudo-functional form of the coupled cluster (CC) type formalisms. Therefore, the new approach may be used to introduce a simple alternative to existing CC-type gradient techniques. The new procedure is formally derived and compared with similar methods from the literature. Model systems calculations (H₂O, LiF, CH⁺ ₂) are further used to demonstrate the effect of various approximations and to elucidate the hierarchy of functional MR-CEPA schemes.