羰基镍簇Ni(CO)n(n=1-4)的结构和Ni-CO键解离性质的密度泛函理论研究

在密度泛函理论框架下, 应用不同泛函计算了配合物Ni(CO)n(n=1~4)的平衡几何构型和振动频率. 考察了泛函和基组重叠误差对预测Ni—CO键解离能的影响. 计算结果表明, 用杂化泛函能得到与实验一致的优化几何构型和较合理的振动频率. 对Ni(CO)n(n=2~4)体系, 用“纯”泛函, 如BP86和BPW91, 可得到与CCSD(T)更符合、 并与实验值接近的解离能. 当解离产物出现单个金属原子或离子(如金属羰基配合物的完全解离)时, BSSE校正项的计算中应保持金属部分的电子结构一致. 只有考虑配体基组和不考虑配体基组两种情况下金属的电子构型与配合物中金属的构型一致时, 才能得到合理的BSSE校正, 从而预测合理的解离能.     

钛,锆和铪羰基化合物M(CO)_n(M=Ti,Zr,Hf;n=4～7)的理论研究

利用密度泛函方法对标题化合物的平衡几何、热化学及振动频率进行了理论预测,发现这3种金属原子都有相似的M(CO)n(n=4～7)结构.全局最低构型对M(CO)7都是单态C3v戴帽八面体7S-1,对M(CO)6都是三重态D3d畸变八面体6T(而对应的单重态M(CO)5仅比它低不到21 kJ·mol-1).对M(CO)n(n=5,4)都是三重态6S-1,其构型分别为从6T中移去1个或2个CO基的衍生物5T和4T.此外,五重态的D3h的三角双锥M(CO)5和单态的Td四面体M(CO)4以及能量更高的含有C和O同时与金属成键的独特配位CO基的M(CO)6和M(CO)3也被发现.最后,给出M(CO)7→M(CO)6+CO反应的离解能.并讨论了金属18价电子的Ti(CO)7存在的可能性.     

ML42-几何和电子结构的理论研究(M=Ni2+,Pd2+,Pt2+;L=CN-,Cl-)

对具有16电子的镍族二价阴离子系列的几何和电子结构进行了研究．在从头算和密度泛函（DFT）的水平上计算得到的几何结构与实验结果具有较好的可比性．SVWN，BP86和B3LYP计算的键长、键角一般具有相同的趋势，但SVWN计算的键长有时偏短．NP2个别情况下表现出异于DFT的趋势．电荷分布计算结果表明：在高自旋态的（扭曲）四面体NiCl4^2-中，Ni具有较高的正电荷，NP2的计算值甚至接近两个正电荷；相对而言，低自旋态的平方四边形中金属原子拥有明显较低的正电荷．计算结果证明了传统观点，即Cl^-是比CN^-更弱的电子受体，这点在NiCl4^2-的情形下体现得特别明显，Cl^-几乎没有为Ni^2＋提供电子或负电荷．     

在高氯酸溶液中Rh(Ⅳ)和Ir(Ⅳ)离子交换行为的研究

本文证实了在高氯酸溶液中Rh(Ⅳ)主要以Rh(H_2O)_6~(4+)形式存在。研究了Rh(H_2O)_6~(4+)和Ir(H_2O)_6~(4+)在阳离子交换树脂上的行为。测定了两者在聚苯乙烯磺酸型阳离子交换树脂—HCIO_4间的分配比。根据Rh(H_2O)_6~(4+)和Ir(H_2O)_6~(4+)在阳离子交换柱上的洗脱曲线和回收率,讨论了在高氯酸溶液中Rh和Ir分离的可能性。     

AuX(X=O,S)低电子态的理论研究

通常要用多参考态方法才能合理处理需考虑旋轨耦合(SOC)效应的开壳层分子如AuO和AuS的低电子态.事实上,通过选取合适的参考态,采用运动方程耦合簇方法(EOM-CC)也能计算这些分子的一些低电子态,而且EOM-CC方法是单参考态方法,使用起来比多参考态方法更加简单.本文采用最近发展的含旋轨耦合的EOM-CC计算电离能的方法(EOMIP-CC),选取对应的负离子为参考态,在CCSD级别上计算了AuO和AuS低电子态的性质.在不考虑旋轨耦合时,通过比较EOMIP-CCSD和EOMIP-CCSDT的结果考察EOMIPCCSD的精度.此外,与EOMIP-CCSDT结果相比,如果自旋污染较为显著而且T1的模较大时,UCCSD(T)方法对能量最低的某一特定对称性的电子态的所对应的电离能误差约为0.1-0.15 eV.在考虑了旋轨耦合效应后,我们的方法得到的键长和振动频率与实验值吻合较好.另一方面,虽然EOMIP-SOC-CCSD高估了能量较高的2Δ3/2态、2Σ+1/2态和2Π1/2态的能量,但是对于其它能量更低的电子态,它们的能量与已有实验值误差在0.2 eV左右.这显示我们所用的含SOC的EOMIP-CCSD方法对原本需要用多参考态方法才能处理的AuO和AuS低电子态能给出可靠的结果.     

乙烯酮自由基(HCCO·)与二氧化氮(NO2)反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//MP2/6-311G(d,p)+ZPE水平上对反应HCCO+NO2进行了计算, 建立了反应势能面. 此反应由反应物通过三步反应到达产物. 首先, NO2的O原子进攻HCCO自由基中与H相邻的C原子, 形成异构体1[ONOC(H)CO]或2[H(CONOC)O]. 然后, 异构体1和2通过N-O键的断裂形成产物NO和OC(H)CO. 最后, 产物中的OC(H)CO可以通过C-C键的断裂进一步分解为HCO和CO. 由HCCO+NO2反应得到产物NO+HCO+CO.     

乙炔基自由基(C2H·)与二氧化氮(NO2)气体反应势能面的理论研究

在 CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPE 水平上对反应C2H+NO2 进行了计算, 建立了反应势能面并得到了3种产物. 利用RRKM理论估算了反应的总速率和分支比. 总速率为1.427×10-12×T0.556×exp(190.547/T) cm3*molecule-1*s-1, 其中主要产物P1(HCCO+NO)比例大于96%, 次要产物P2(HCNO+CO)和P3(HCN+CO2)小于4%.     

硫族铅化物阳离子低电子态的理论研究

采用我们最近发展的含旋轨耦合的运动方程耦合簇计算电离能(EOMIP-CC)方法,在CCSD级别上计算了硫族铅化物PbS、PbSe、PbTe阳离子低电子态的平衡键长和谐振频率以及绝热和垂直电离能,得到的结果与已有的实验值吻合较好.不考虑旋轨耦合(SOC)的情况下通过与CCSD(T)的计算结果比较,考察了三重激发对计算结果的影响,结果显示考虑三重激发的贡献后得到的键长和频率结果与实验值吻合更好.计算结果表明PbTe+中2Π态的能量分裂明显大于PbS+和PbSe+中2Π态的能量分裂,但是PbTe+中2Π1/2和2Σ1/2态之间的相互耦合则明显弱于PbS+和PbSe+中这两个态之间的耦合.PbTe+中2Π1/2和2Σ1/2态之间耦合很弱,一方面是因为2Σ+态和2Π态的能量差比PbS+和PbSe+中2Σ+态和2Π态的能量差大,另一方面还由于PbTe+中2Π1/2和2Σ1/2态之间的旋轨耦合矩阵元只是PbS+和PbSe+中2Π1/2和2Σ1/2态之间的旋轨耦合矩阵元的一半.这些计算结果为PbS+、PbSe+、PbTe+阳离子的低电子态性质提供了新的理论数据,可以为将来的实验数据提供参考.     

(Mo_3S_4)_2M~(8+)(M=Mo,Sn,Cu_2)型夹心簇合物的电子结构研究

经对一系列夹心簇合物(H_2O)_9Mo_3S_4MS_4Mo_3(H_2O)_9~(8+)(M=Mo,Sn,Cu_2)的量子化学计算和分析,发现该簇合物是通过配体Mo_3S_4~(4+)中的多中心d-pπ轨道与中心原子键合成簇,其成簇机理与经典的以有机共轭环为配体的夹心配合物十分相似,研究了夹心簇的构型、磁性及电子光谱等性质的分子轨道本质,揭示了多中心d-pπ键的变化与簇合物性质的内在联系。     

在高氯酸溶液中Rh(Ⅲ)和Ir(Ⅲ)的离子交换行为及分离的研究

本文研究了在高氯酸溶液中,Rh(H_2O)_6~(3+)和Ir(H_2O)_6~(3+)在强酸性阳离子交换树脂上的行为,测得在不同酸度的高氯酸溶液中Rh(H_2O)_6~(3+)和Ir(H_2O)_6~(3+)的分配比。分别测得Rh(H_2O)_6~(3+)和Ir(H_2O)_6~(3+)在柱上的洗脱曲线和分离曲线,讨论了在高氯酸中Rh(Ⅲ)和Ir(Ⅲ)分离的可能性。     

A truncated version of reduced multireference coupled-cluster method with singles and doubles and noniterative triples: application to F2 and Ni(CO)n (n=1, 2, and 4)

A perturbatively truncated version of the reduced multireference coupled-cluster method with singles and doubles and noniterative triples RMR CCSD(T) is described. In the standard RMR CCSD method, the effect of all triples and quadruples that are singles or doubles relative to references spanning a chosen multireference (MR) model space is accounted for via the external corrections based on the MR CISD wave function. In the full version of RMR CCSD(T), the remaining triples are then handled via perturbative corrections as in the standard, single-reference (SR) CCSD(T) method. By using a perturbative threshold in the selection of MR CISD configuration space, we arrive at the truncated version of RMR CCSD(T), in which the dimension of the MR CISD problem is significantly reduced, thus leaving more triples to be treated perturbatively. This significantly reduces the computational cost. We illustrate this approach on the F2 molecule, in which case the computational cost of the truncated version of RMR CCSD(T) is only about 10%-20% higher than that of the standard CCSD(T), while still eliminating the failure of CCSD(T) in the bond breaking region of geometries. To demonstrate the capabilities of the method, we have also used it to examine the structure and binding energy of transition metal complexes Ni(CO)n with n=1, 2, and 4. In particular, Ni(CO)2 is shown to be bent rather than linear, as implied by some earlier studies. The RMR CCSD(T) binding energy differs from the SR CCSD(T) one by 1-2 kcal/mol, while the energy barrier separating the linear and bent structures of Ni(CO)2 is smaller than 1 kcal/mol.     

Benchmark calculations on the electron detachment energies of MO3* and M2O6* (M = Cr, Mo, W)

Neutral and anionic molecules of the monomers and dimers of the group VIB transition metal oxides (MO3 and M2O6) were studied with density functional theory (DFT) and coupled cluster CCSD(T) theory. Franck-Condon simulations of the photoelectron spectra were carried out for the transition from the ground state of the anion to that of the neutral molecule. Molecular structures from the DFT and CCSD(T) methods are compared. Electron detachment energies reported in the literature were evaluated. The calculated adiabatic and vertical electron detachment energies (ADEs and VDEs) were compared with the experimental results. CCSD(T) gives results within 0.12 eV for the ADEs. CCSD(T) predicts VDEs that are in error by as much as 0.3 eV for M = Cr. DFT hybrid functionals were found to give poor results for the ADEs and VDEs for M = Cr due to the substantial amount of multireference character in the wavefunction, whereas the pure DFT functionals give superior results. For M = Mo and W, excellent agreement was found for both CCSD(T) and many DFT fucntionals. The BP86 functional yields the best overall results for the VDEs of all the metal oxide clusters considered. Heats of formation calculated at the CCSD(T) level extrapolated to the complete basis set limit are also in good agreement with available experimental data.     

Density functional theory and ab initio studies of the structure and energetics of digallium tetraoxide, Ga2O4, in the gas phase

A systematic investigation on the neutral and anionic digallium tetraoxide, Ga(2)O(4) has been carried out by using density functional theory (DFT), second-order M?ller-Plesset perturbation theory (MP2), and the coupled cluster approach with single and double substitutions and a perturbative treatment of the triple excitations [CCSD(T)]. The geometry of neutral Ga(2)O(4) has been proposed earlier, from an experimental study, to adopt a D(2d) symmetry (J. Phys. Chem. 1979, 83, 656). However, the current research reveals that, out of the several isomers considered for neutral and anionic digallium tetraoxide, the (3)B(1u) and (2)B(3g) of the planar D(2h) geometry (7a-D(2h)) are the lowest-energy states for Ga(2)O(4) and Ga(2)O(4)(-). Our computations rule out the D(2d) geometry (3-D(2d)) as a viable contender for neutral Ga(2)O(4). The (3)B(2) (3-D(2d)) state is located above the (3)B(1u) (7-D(2h)) state by at least 4.26 eV. The energies of the low lying states, geometrical parameters, and energetic features (VEDE, AEDE, and AEA) are reported. The AEA of Ga(2)O(4) is calculated to be 3.94 eV (B3LYP), 3.24 eV (MP2), 3.42 eV [CCSD(T)//B3LYP], and 3.38 eV [CCSD(T)//MP2], respectively. In addition, the dissociation energy, D(e), for the process Ga(2)O(4) ((3)B(1u)) → 2GaO(2) ((2)A(2)) is 3.59 eV (B3LYP), 5.08 eV (MP2), 4.82 eV [CCSD(T)//B3LYP], and 4.80 eV [CCSD(T)//MP2]. The results obtained in this work are critically analyzed, discussed, and compared with those of the analogous metal oxides.     

紫外激光诱导C60与M(CO)6(M=Cr,Mo, W)的配位反应

研究了用紫外激光(355 nm)诱导C_(60)与金属羰基化合物M(CO)_6(M＝Cr；Mo，W)的配位反应，合成了具有C_(4υ)对称性的配位络合物(η~2－C_(60))M(CO)_5，并初步讨论了C_(60)与M(CO)_6反应的动力学过程．     

Theoretical study on the reaction mechanism of the gas-phase H2/CO2/Ni(3D) system

The ground-state potential energy surface (PES) in the gas-phase H2/CO2/Ni(3D) system is investigated at the CCSD(T)//B3LYP/6-311+G(2d,2p) levels in order to explore the possible reaction mechanism of the reverse water gas shift reaction catalyzed by Ni(3D). The calculations predict that the C-O bond cleavage of CO2 assisted by co-interacted H2 is prior to the dissociation of the H2, and the most feasible reaction path for Ni(3D) + H2 + CO2 --> Ni(3D) + H2O + CO is endothermic by 12.5 kJ mol(-1) with an energy barrier of 103.9 kJ mol(-1). The rate-determining step for the overall reaction is predicted to be the hydrogen migration with water formation. The promotion effect of H2 on the cleavage of C-O bond in CO2 is also discussed and compared with the analogous reaction of Ni(3D) + CO2 --> NiO + CO, and the difference between triplet and singlet H2/CO2/Ni systems is also discussed.     

Low-lying electronic states of M(3)O(9)(-) and M(3)O(9)(2-) (M = Mo, W)

Multiple low-lying electronic states of M(3)O(9)(-) and M(3)O(9)(2-) (M = Mo, W) arise from the occupation of the near-degenerate low-lying virtual orbitals in the neutral clusters. We used density functional theory (DFT) and coupled cluster theory (CCSD(T)) with correlation consistent basis sets to study the structures and energetics of the electronic states of these anions. The adiabatic and vertical electron detachment energies (ADEs and VDEs) of the anionic clusters were calculated with 27 exchange-correlation functionals including one local spin density approximation functional, 13 generalized gradient approximation (GGA) functionals, and 13 hybrid GGA functionals, as well as the CCSD(T) method. For M(3)O(9)(-), CCSD(T) and nearly all of the DFT exchange-correlation functionals studied predict the (2)A(1) state arising from the Jahn-Teller distortion due to singly occupying the degenerate e' orbital to be lower in energy than the (2)A(1)' state arising from singly occupying the nondegenerate a(1)' orbital. For W(3)O(9)(-), the (2)A(1) state was predicted to have essentially the same energy as the (2)A(1)' state at the CCSD(T) level with core-valence correlation corrections included and to be higher in energy or essentially isoenergetic with most DFT methods. The calculated VDEs from the CCSD(T) method are in reasonable agreement with the experimental values for both electronic states if estimates for the corrections due to basis set incompleteness are included. For M(3)O(9)(2-), the singlet state arising from doubly occupying the nondegenerate a(1)' orbital was predicted to be the most stable state for both M = Mo and W. However, whereas M(3)O(9)(2-) was predicted to be less stable than M(3)O(9)(-), W(3)O(9)(2-) was predicted to be more stable than W(3)O(9)(-).     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

Structural and electronic properties of reduced transition metal oxide clusters, M4O10 and M4O10- (M = Cr, W), from photoelectron spectroscopy and quantum chemical calculations

Anion photoelectron spectroscopy and quantum chemical calculations at the density functional theory (DFT), coupled cluster theory (CCSD(T)), and complete active space self-consistent field (CASSCF) theory levels are employed to study the reduced transition metal oxide clusters M(4)O(10)(-) (M = Cr, W) and their neutrals. Photoelectron spectra are obtained at 193 and 157 nm photon energies, revealing very different electronic structures for the Cr versus W oxide clusters. The electron affinity and HOMO-LUMO gap are measured to be 3.68 ± 0.05 and 0.7 eV, respectively, for the Cr(4)O(10) neutral cluster, as compared to 4.41 ± 0.04 and 1.3 eV for W(4)O(10). A comprehensive search is performed to determine the ground-state structures for M(4)O(10) and M(4)O(10)(-), in terms of geometry and electronic states by carefully examining the calculated relative energies at the DFT, CCSD(T), and CASSCF levels. The ground states of Cr(4)O(10) and Cr(4)O(10)(-) have tetrahedral structures similar to that of P(4)O(10) with the anion having a lower symmetry due to a Jahn-Teller distortion. The ground states of W(4)O(10) and W(4)O(10)(-) have butterfly shape structures, featuring two fused five-member rings with a metal-metal multiple bond between the central metal atoms. The much stronger WW bonding than the CrCr bonding is found to be the primary cause for the different ground state structures of the reduced Cr(4)O(10)(0/-) versus W(4)O(10)(0/-) oxide clusters. The photoelectron spectra are assigned by comparing the experimental and theoretical adiabatic and vertical electron detachment energies, further confirming the determination of the ground electronic states of M(4)O(10) and M(4)O(10)(-). The time-dependent DFT method is used to calculate the excitation energies of M(4)O(10). The TD-DFT results in combination with the self-consistently calculated vertical detachment energies for some of the excited states at the DFT and CCSD(T) levels are used to assign the higher energy bands. Accurate clustering energies and heats of formation of M(4)O(10) are calculated and used to calculate accurate reaction energies for the reduction of M(4)O(12) to M(4)O(10) by CH(3)OH, as well as for the oxidation of M(4)O(10) to M(4)O(12) by O(2). The performance of the DFT method with the B3LYP and BP86 functionals in the calculations of the relative energies, electron detachment energies, and excitation energies are evaluated, and the BP86 functional is found to give superior results for most of these energetic properties.     

Ab initio pseudopotential study of M2As− and M2Br+ (M = Cu,Ag, Au)

The equilibrium geometries and the vibration frequencies of M₂As⁻ and M₂Br⁺ (M = Cu, Ag, Au) are calculated at the Hartree–Fock (HF) and the second‐order Møller–Plesset (MP2) levels with pseudopotentials. The calculated results indicate that the species have a bent structure (C_2v). The electron correlation corrections on the geometrical structure are investigated at the MP2 level, the bond angles are reduced by 10°–20° for considered species. The electron correlation effects on the geometry of the Au₂As⁻ are studied particularly at MP2, MP3, MP4, CCSD and CCSD(T) levels. Comparing the species containing Ag and Au, the relativistic effects slightly short the bond lengths of the species. The bonding possibility of the Au₂As⁻ is predicted. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 80: 38–43, 2000