基于CASSCF参考函数的块相关耦合簇方法对烷烃中单键解离势能面研究

用以完全活化空间自洽场(CASSCF)波函数为参考波函数的块相关耦合簇(BCCC)方法(简称CAS-BCCC)研究了烷烃(甲烷和乙烷)中的单键解离过程的势能面(PES). 与其它理论方法比较的结果表明, 该方法可以对所研究的整个解离势能面给出定量准确的描述.     

Alternative single-reference coupled cluster approaches for multireference problems: the simpler, the better

We report a general implementation of alternative formulations of single-reference coupled cluster theory (extended, unitary, and variational) with arbitrary-order truncation of the cluster operator. These methods are applied to compute the energy of Ne and the equilibrium properties of HF and C(2). Potential energy curves for the dissociation of HF and the BeH(2) model computed with the extended, variational, and unitary coupled cluster approaches are compared to those obtained from the multireference coupled cluster approach of Mukherjee et al. [J. Chem. Phys. 110, 6171 (1999)] and the internally contracted multireference coupled cluster approach [F. A. Evangelista and J. Gauss, J. Chem. Phys. 134, 114102 (2011)]. In the case of Ne, HF, and C(2), the alternative coupled cluster approaches yield almost identical bond length, harmonic vibrational frequency, and anharmonic constant, which are more accurate than those from traditional coupled cluster theory. For potential energy curves, the alternative coupled cluster methods are found to be more accurate than traditional coupled cluster theory, but are three to ten times less accurate than multireference coupled cluster approaches. The most challenging benchmark, the BeH(2) model, highlights the strong dependence of the alternative coupled cluster theories on the choice of the Fermi vacuum. When evaluated by the accuracy to cost ratio, the alternative coupled cluster methods are not competitive with respect to traditional CC theory, in other words, the simplest theory is found to be the most effective one.     

Photoluminescence of oxygen-containing surface defects in germanium oxides: a theoretical study

Photoabsorption and photoluminescence properties of nonbridging oxygen -O-Ge[triple bond](NBO), -OO-Ge[triple bond] (peroxy radical), O=Ge=, and (O2)Ge= defects in germanium oxides have been investigated by high-level ab initio calculations. Geometry optimization for excited electronic states of model clusters simulating these defects was carried out at the complete-active-space self-consistent-field level, and relative energies were calculated by various methods including time-dependent density-functional theory, outer-valence Green's functions, equation-of-motion coupled cluster theory with single and double excitations, symmetry-adapted cluster configuration interaction, multireference second-order perturbation theory, and multireference configuration interaction. The results demonstrate that the considered excited states of the aforementioned defects normally exhibit large Stokes shifts and that, with few exceptions, UV photoabsorption is accompanied by red or IR photoluminescence.     

Detailed dynamics of the photodissociation of cyclobutane

Semiclassical electron-radiation-ion dynamics simulations are reported for the photodissociation of cyclobutane into two molecules of ethylene. The results clearly show the formation of the tetramethylene intermediate diradical, with dissociation completed in approximately 400 fs. In addition, the potential energy surfaces of the electronic ground state and lowest excited-state were calculated at the complete-active-space self-consistent-field/multireference second-order perturbation theory (CASSCF/MRPT2) level with 6-31G* basis sets, along the reaction path determined by the dynamics simulations. There are well-defined energy minima and maxima in the intermediate state region. It is found that both C-C-C bond bending and rotation of the molecule (around the central C-C bond) have important roles in determining the features of the potential energy surfaces for the intermediate species. Finally, the simulations and potential energy surface calculations are applied together in a discussion of the full mechanism for cyclobutane photodissociation.     

N[bond]C(alpha) bond dissociation energies and kinetics in amide and peptide radicals. Is the dissociation a non-ergodic process?

Dissociations of aminoketyl radicals and cation radicals derived from beta-alanine N-methylamide, N-acetyl-1,2-diaminoethane, N(alpha)-acetyl lysine amide, and N(alpha)-glycyl glycine amide are investigated by combined density functional theory and M?ller-Plesset perturbational calculations with the goal of elucidating the mechanism of electron capture dissociation (ECD) of larger peptide and protein ions. The activation energies for dissociations of N[bond]C bonds in aminoketyl radicals decrease in the series N[bond]CH(3) > N-CH(2)CH(2)NH(2) > N[bond]CH(2)CONH(2) approximately N[bond]CH(CONH(2))(CH(2))(4)NH(2). Transition state theory rate constants for dissociations of N[bond]C(alpha) bonds in aminoketyl radicals and cation-radicals indicate an extremely facile reaction that occurs with unimolecular rate constants >10(5) s(-1) in species thermalized at 298 K in the gas phase. In neutral aminoketyl radicals the N[bond]C(alpha) bond cleavage results in fast dissociation. In contrast, N[bond]C(alpha) bond cleavage in aminoketyl cation-radicals results in isomerization to ion-molecule complexes that are held together by strong hydrogen bonds. The facile N[bond]C(alpha) bond dissociation in thermalized ions indicates that it is unnecessary to invoke the hypothesis of non-ergodic behavior for ECD intermediates.     

Density Functional Theory for the Investigation of Catalytic Activity of X@Cu12 (X =Cu,Ni, Co,or Fe) for N2O Decomposition

We have investigated the reaction mechanism for N2 O decomposition on Cu13 via density functional theory. It is found that N2 O decomposition on the cluster is more prone to be along the Eley-Rideal(ER) pathway in comparison with the Langmuir-Hinshelwood(LH) channel. There exists structural relaxation for Cu13 cluster in the reaction, which may influence the catalytic activity of cluster for the subsequent N2 O decomposition. The core atom in the Cu13 cluster is substituted with the Fe, Co, or Ni to enhance structural stability and prevent from the obvious configuration relaxation in the reaction. Note that these bimetallic clusters are of icosahedra as the Cu13. They have activities for N2 O dissociation along ER pathway and the heteroatom in the cluster can prevent configuration from relaxation. Finally, the Ni@Cu12 cluster can be as a superior catalyst in a complete catalytic cycle via comparison in this study.     

Comparison of some multireference electronic structure methods in illustrative applications

The performances of several multireference electronic structure methods including complete active space self-consistent field (CASSCF)-based second-order perturbation theory (CASPT2), multireference configuration interaction with single and double excitations (MR-CISD), MR-CISD with the Davidson correction (MR-CISD+Q), and the CASSCF-based block-correlated coupled cluster method (CAS-BCCC4) we developed recently are compared by applying them to study several different chemical problems involving computation...     

Heteronuclear Trimetallic MFe2 and M2Fe (M=V,Nb, and Ta) Clusters for Dinitrogen Activation

Catalysts with heteronuclear metal active sites may have high performance in the nitrogen reduction reaction (NRR), and the in-depth understanding of the reaction mechanisms is crucial for the design of related catalysts. In this work, the dissociative adsorption of N₂ on heteronuclear trimetallic MFe₂ and M₂Fe (M=V, Nb, and Ta) clusters was studied with density functional theory calculations. For each cluster, two reaction paths were studied with N₂ initially on M and Fe atoms, respectively. Mayer bond order analysis provides more information on the activation of N−N bonds. M₂Fe is generally more reactive than MFe₂. The coordination mode of N₂ on three metal atoms can be end-on: end-on: side-on (EES) for both MFe₂ and M₂Fe. In addition, a unique end-on: side-on: side-on (ESS) coordination mode was found for M₂Fe, which leads to a higher degree of N−N bond activation. Nb₂Fe has the highest reactivity towards N₂ when both the transfer of N₂ and the dissociation of N−N bonds are taken into account, while Ta-containing clusters have a superior ability to activate the N−N bond. These results indicate that it is possible to improve the performance of iron-based catalysts by doping with vanadium group metals.     

Modern valence-bond description of the electronic structure of benzocyclobutadiene

The ground state of benzocyclobutadiene, a bicyclic molecule with 8π electrons containing a benzene and a cyclobutadiene ring, is studied by means of modern valence bond (VB) theory in its spin-coupled (SC) form and the complete-active-space self-consistent field (CAS SCF) approach. The CAS SCF wave function is used to optimize the geometry, and SC theory—to obtain a well-correlated and easy to visualize and understand model of the active space hosting the π electrons. It is shown that the π system of the ground state of benzocyclobutadiene can be described with sufficient accuracy as a combination of the π systems of a distorted benzene ring and an isolated double bond. Each of the eight SC orbitals is found to be well-localized about one carbon atom only, with small distortions toward its nearest neighbors. The analysis of the optimal active-space spin-coupling pattern within the SC wave function for benzocyclobutadiene shows unambiguously that this molecule inherits neither the aromatic nor the antiaromatic character of its cyclic components, and it is most appropriate to regard it as nonaromatic. © 1996 John Wiley & Sons, Inc.     

Weak metal-metal bonding in small manganese cluster ions, Mn(N)+ (N < or = 7)

The binding energies of manganese cluster ions Mn(N)+ (N = 5-7) were determined by the photodissociation experiments in the near-infrared and visible-photon-energy ranges. The bond dissociation energies of Mn(N)+, D0(Mn(N-1)+...Mn), were obtained to be 1.70+/-0.08, 1.04+/-0.10, and 1.46+/-0.11 eV, respectively, for N = 5, 6, and 7 from the threshold energies for the two-atom loss processes and the bond dissociation energies of Mn3(+) and Mn4(+) reported previously [A. Terasaki et al., J. Chem. Phys. 117, 7520 (2002)]. Correspondingly, binding energies per atom are obtained to be 0.99+/-0.03, 1.00+/-0.03, and 1.06+/-0.03 eV/at. for N = 5, 6, and 7, respectively. A gradual increase in the binding energy from N = 2 to N = 7 shows an increasing contribution of nonbonding 3d orbitals to the bonding via weak hybridization with valence 4s orbitals as the cluster size increases. These binding energies per atom are still much smaller than the bulk cohesive energy of manganese (2.92 eV/at.), and this finding indicates exceptionally weak metal-metal bonds in this size range.     

H+、NH+4对HMX的N—NO2键解离能的影响

采用密度泛函理论的B3P86/6-31G**方法, 优化了β-HMX及其与H⁺、NH^+₄分别形成的复合物的稳定结构, 计算了β-HMX以及复合物中最弱的N—NO₂键解离能. 结果发现, HMX与H⁺、NH^+₄形成复合物后, 使HMX的构型产生较大变化; 与H⁺结合后, HMX的一个N—NO₂键显著伸长, 键级变小; 但与NH^+₄形成复合物后, HMX中键级最小的N—NO₂键长变化不大. 键解离能计算表明, 同β-HMX相比, 与H⁺形成的两种复合物中N—NO₂键解离能分别降低了近20 和82 kJ·mol^-1, 而HMX与NH^+₄形成的复合物中N—NO₂键解离能仅降低了约8 kJ·mol^-1, 表明H⁺对β-HMX的N—NO₂键的初始热裂解反应有促进作用, 而NH^+₄影响不明显.     

1,2,4,5-四嗪含氮取代基衍生物结构和性质的理论研究

在B3LYP/aug-cc-pvDZ理论水平上研究了CN,NO2,NH2,N3,N2H,NHNH2,N4H和N4H3含氮取代基取代1,2,4,5-四嗪环上的两个氢原子生成的衍乍物,预测了它们的分f构犁、分解能及含能性质.对衍生物分解能的研究结果表明.CN取代的衍生物的分解能比未取代时更高,而其余基团的取代使分解能降低.生成热的研究显示取代基化合物的生成热越大,取代1,2,4,5-四嗪中的氢原子后生成衍生物的牛成热也越大;CN,N3和N4H取代的1,2,4,5-四嗪衍生物的单位原子生成热在83.1～95.2 kJ,比文献报道的三叠氮基-均三嗪的(70.2 kJ)更高;N4H,N3,N4H3,N2H和CN取代的1,2,4,5-四嗪衍生物,生成热在904.9～1496.6 kJ·mol-1,但N4H和N4H3取代的衍生物分解能较小,稳定性较差.     

The vibrational frequencies of CaO2, ScO2, and TiO2: a comparison of theoretical methods

The vibrational frequencies of several states of␣CaO₂, ScO₂, and TiO₂ are computed using density functional theory (DFT), the Hartree-Fock approach, second-order M?ller-Plesset perturbation theory (MP2), and the complete-active-space self-consistent-field theory. Three different functionals are used in the DFT calculations, including two hybrid functionals. The coupled cluster singles and doubles approach including the effect of connected triples, determined using perturbation theory, is applied to selected states. The Becke-Perdew 86 functional appears to be the most cost-effective method of choice, although even this functional does not perform well for one state of CaO₂. The MP2 approach is significantly inferior to the DFT approaches. Received: 3 September 1997 / Accepted: 8 December 1997     

Hydrogen-atom abstraction from the adenine-uracil base pair

The hydrogen-abstracted radicals from the adenine-uracil (AU) base pair have been studied at the B3LYP/DZP++ level of theory. The A(N9)-U and A-U(N1) radicals, which correspond to hydrogen-atom abstraction at the adenine N9 and uracil N1 atoms, respectively, were predicted to be the two lowest-lying among the nine (AU-H) radicals studied in this study. The removal of the amino hydrogen of the adenine moiety that forms a hydrogen bond with the uracil O4 atom in the AU pair resulted in radical A(N6a)-U, which has the smallest base-pair dissociation energy, 5.9 kcal mol(-1). This radical is more likely to dissociate into the two isolated bases than to recover the hydrogen bond with the O4 atom through N6-H bond rotation along the C6-N6 bond. In general, the radicals generated by C-H bond breaking were higher in energy than those arising from N-H bond cleavage, because the unpaired electrons in the carbon-centered radicals were mainly localized on the carbon atom from which the hydrogen atom was removed. However, the highest-lying radical was found to arise from removal of the N3 hydrogen of uracil. The most remarkable structural feature of this radical is a very short C-H...O distance of 2.094 A, consistent with a substantial hydrogen bond. Although this radical lost the N1...H-N3 hydrogen bond between the two bases, its dissociation energy was predicted to be 12.9 kcal mol(-1), similar to that of the intact AU base pair. This is due to the transfer of electron density from the adenine N1 atom to the uracil N3 atom.     

NO在MgO(001)缺陷表面吸附解离的量子化学研究

在密度泛函理论框架下,采用嵌入点电荷簇模型研究了NO在MgO(001)完整和缺陷表面上的吸附。研究结果表明:具有氧缺陷结构表面的催化活性较高,有利于NO键的削弱;当另一个NO分子进攻已吸附的NO分子时,NO键将进一步削弱,直致断裂,并伴有N2O产生,这与UPS和MIES实验观察到的现象一致。Mulliken布居分析指出,底物电子向NO转移,并填充到NO的*反键轨道上,从而导致NO键的削弱,并形成NO-。这也是可能导致形成NO-的原因。研究还表明,具有镁缺陷的MgO(001)表面对NO的解离没有催化活性。     

Ab initio molecular orbital study of Fe(CO) n (n = 1–3)

Geometry optimization was performed for the ground states of FeCO, Fe(CO)₂, and Fe(CO)₃ at various levels of ab initio calculations, and the bond lengths and dissociation energies obtained were in reasonable agreement with experimental results. The nature of bonding was studied for these molecules using a complete-active-space self-consistent-field method. From the Mulliken population analysis, it was found that the traditional donation and back donation mechanism is valid for these molecules, including Fe(CO)₃, which has a pyramidal structure. Received: 27 September 1999 / Accepted: 13 January 2000 / Published online: 19 April 2000