High magnetic exchange coupling constants: a density functional theory based study of substituted schlenk diradicals

The Schlenk diradical has been known since 1915. After a detailed experimental work by Rajca, its magnetic nature has remained more or less unexplored. We have investigated by quantum chemical calculations the nature of magnetic coupling in 11 substituted Schlenk diradicals. Substitution has been considered at the fifth carbon atom of the meta-phenylene moiety. The UB3LYP method has been used to study 12 diradicals including the original one. The 6-311G(d,p) basis set has been employed for optimization of molecular geometry in both singlet and triplet states for each species. The singlet optimization has led to the optimization of the broken-symmetry structure for 10 species including the unsubstituted one. This development makes it possible to carry out further broken symmetry calculations in two ways. The triplet calculation has been done using 6-311++G(d,p) basis set and the optimized triplet geometry in both procedures. The broken symmetry calculations have used the optimized geometries of either the triplet states or the broken symmetry solutions. The first method leads to the prediction of electron paramagnetic resonance (EPR) compatible magnetic exchange coupling constant (J) in the range 517-617 cm(-1). A direct optimization of the broken symmetry geometry gives rise to a lower estimate of J, in the range of 411-525 cm(-1) and compatible with macroscopic Curie studies. The calculated J for the unsubstituted Schlenk diradical is 512 cm(-1) that can be compared with 455 cm(-1) estimated by Rajca. In both cases, introduction of groups with +M and +I effects (Ingold's notation) decreases the J value from that for the unsubstituted Schlenk diradical while -I and -M groups at the same position increases J. These trends have been explained in terms of Hammett constants, atomic spin densities, and dihedral angles.     

FC(O)O自由基与NO2反应的量子化学研究

用量子化学DFT, MP2, G3和G3MP2方法对FC(O)O自由基与NO₂的反应机理进行了理论研究. 优化了反应势能面上各驻点的几何结构, 通过内禀反应坐标(IRC)计算和振动分析, 确认了反应中的过渡态, 并用过渡态理论(TST)计算了相关反应的速率常数.     

HPO3异构体结构与稳定性

使用密度泛函方法(B3LYP)和6-311 G(3df,3pd)基组计算得到了HPO3体系的9个异构体与11个过渡态,并用内禀反应坐标理论验证了异构体的异构化过程．结果表明,在HPO3体系中,只有热力学最稳定的平面型异构体HOFO2(E1),(cis,cis)构象的HOOPO(E3)和立体的具有Cs对称性的异构体HP(O)O2(E7)具有较高的动力学稳定性,理论预测结果与实验一致．另外2个立体的HOFO2连接方式的异构体(E2和E4)由于解离为HO(Ⅱ) OPO(^2A1)的解离能较低,因此不能稳定存在。     

Microwave spectrum of 3-butyne-1-thiol: evidence for intramolecular S-H...pi hydrogen bonding

The microwave spectrum of 3-butyne-1-thiol has been studied by means of Stark-modulation microwave spectroscopy and quantum-chemical calculations employing the B3LYP/6-311++G(3df,2pd), MP2/aug-cc-pVTZ, MP2/6-311++G(3df,2pd), and G3 methods. Rotational transitions attributable to two conformers of this molecule were assigned. One of these conformers possesses an antiperiplanar arrangement of the atoms S-C1-C2-C3, while the other is synclinal and stabilized by the formation of an intramolecular hydrogen bond between the H-atom of the thiol group and the pi-electrons of the C[triple bond]C triple bond. The energy difference between these conformers was estimated to be 1.7(4) kJ mol(-1) by relative intensity measurements, with the hydrogen-bonded conformer being lower in energy. The spectra of five vibrationally excited states of the synclinal conformer were observed, and an assignment of these states to particular vibrational modes was made with the aid of a density functional theory (DFT) calculation of the vibrational frequencies at the B3LYP/6-311++G(3df,2pd) level of theory.     

Basis set and method dependence in quantum theory of atoms in molecules calculations for covalent bonds

The influence of various small- and medium-size basis sets used in Hartree-Fock (HF) and density functional theory (DFT)/B3LYP calculations on results of quantum theory of atoms in molecules based (QTAIM-based) analysis of bond parameters is investigated for several single, double, and triple covalent bonds. It is shown that, in general, HF and DFT/B3LYP methods give very similar QTAIM results with respect to the basis set. The smallest 6-31G basis set and DZ-quality basis sets of Dunning type lead to poor results in comparison to those obtained by the most reliable aug-cc-pVTZ. On the contrary, 6-311++G(2df,2pd) and in a somewhat lesser extent 6-311++G(3df,3pd) basis sets give satisfactory values of QTAIM parameters. It is also demonstrated that QTAIM calculations may be sensitive for the method and basis set in the case of multiple and more polarized bonds.     

The molecular structure of fluoromalononitrile as determined by gas-phase electron diffraction and ab initio calculations

The molecular structure of fluoromalononitrile was studied by means of gas-phase electron diffraction and quantum mechanical methods using HF/6-31G(d), MP2/6-311++G(2df,2pd) and DFT/B3LYP/6-31G(d), B3PW91/6-31G(d), B3LYP/6-311++G(2df,2pd) and B3PW91/6-311++G(2df,2pd). The r(g) and angle(alpha) structural parameters we obtained from the present analysis are: CC=1.487(5) A, CN=1.157(3) A, CF=1.386(5) A, CH=1.096 A (ass.), angleCCC=106.7(1.0) degrees , angleCCF=108.0(0.7) degrees , angleCCN=177.6(2.0) degrees . Uncertainties in parenthesis are 3sigma.     

Theoretical investigations on analytical potential energy function and spectroscopic parameters for the state b3Πu of dimer 7Li2

The SAC‐CI (symmetry‐adapted‐cluster configuration‐interaction) method presented in Gaussian 03 program package is applied to investigate the adiabatic potential energy curves (PECs) of ⁷Li₂(b³Π_u). These calculations are performed at numbers of basis sets, such as 6‐311++G(3df,3pd), 6‐311++G(2df,2pd), 6‐311++G(df,pd), D95V++, D95(3df,3pd), D95(d,p), cc‐PVTZ, 6‐311++G and 6‐311++G(d,p). All the ab initio calculated points are fitted to the analytic Murrell‐Sorbie functions and then used to compute the spectroscopic parameters. The analytic potential energy function (APEF) for this b³Π_u state is reported. By comparison, the spectroscopic parameters reproduced by the APEF attained at 6‐311++G(2df,2pd) are found to be very close to the latest experimental findings. With the APEF obtained at the SAC‐CI/6‐311++G(2df,2pd) level of theory, a total of 62 vibrational states is found when J = 0. The complete vibrational levels, classical turning points, inertial rotation and centrifugal distortion constants for these vibrational states are also reported. The reasonable dissociation limit for this state is deduced using the calculated results at present. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Formic acid catalyzed hydrolysis of SO3 in the gas phase: a barrierless mechanism for sulfuric acid production of potential atmospheric importance

Computational studies at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels are performed to explore the changes in reaction barrier height for the gas phase hydrolysis of SO(3) to form H(2)SO(4) in the presence of a single formic acid (FA) molecule. For comparison, we have also performed calculations for the reference reaction involving water assisted hydrolysis of SO(3) at the same level. Our results show that the FA assisted hydrolysis of SO(3) to form H(2)SO(4) is effectively a barrierless process. The barrier heights for the isomerization of the SO(3)···H(2)O···FA prereactive collision complex, which is the rate limiting step in the FA assisted hydrolysis, are found to be respectively 0.59 and 0.08 kcal/mol at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels. This is substantially lower than the ~7 kcal/mol barrier for the corresponding step in the hydrolysis of SO(3) by two water molecules--which is currently the accepted mechanism for atmospheric sulfuric acid production. Simple kinetic analysis of the relative rates suggests that the reduction in barrier height facilitated by FA, combined with the greater stability of the prereactive SO(3)···H(2)O···FA collision complex compared to SO(3)···H(2)O···H(2)O and the rather plentiful atmospheric abundance of FA, makes the formic acid mediated hydrolysis reaction a potentially important pathway for atmospheric sulfuric acid production.     

CH2=CHCl与O(3P)反应的理论研究

用量子化学密度泛函理论和QCISD(Quadratic configuration interaction calculation)方法,对0(^3P)与CH2CHCl的反应进行了理论研究．在UB3LYP／6—311 G(d,p),UB3LYP／6—31 (3df,3pd)计算水平上,优化了反应物、产物、中间体和过渡态的几何构型,并在UQCISD(T)／6—311 G(2df,2pO)水平上计算了单点能量．为了确证过渡态的真实性,在UB3LYP／6—311 G(3df,3pd)水平上进行了内禀坐标(IRC)计算和频率分析,并确定了反应机理．研究结果表明,反应主要产物为CH2CHO和Cl．     

Ar3H＋中多体势能拆分的从头计算分析

使用从头算的方法计算了Ar3H＋的稳定构型及红外振动基频，并与Ar2H＋的稳定结构及红外振动频率作了对比，讨论了二者之间的联系与差别.对Ar3H＋一个特殊构型下的势能线作了扫描，在此势能线下将Ar3H＋的四体势能拆解为两体，半三体及三体势能的加和，初步探讨了高阶势能分解为低阶势能的最佳途径.     

Two Novel Isomers of HPS3 System

Two new isomers of HPS3 system, FIP(S)S2 and HSSPS, are predicted by means of B3LYP method with 6-311 G(3df,3pd) basis set. The two isomers can isomerize into thermodynamically the most stable species HSPS2, which have been experimentally identified,with relatively higher reaction barriers. In view of their higher thermodynamical and kinetic stability and the experimental observation for I-IP(O)O2 and HOOPO in previous study, we can reasonably believe that the two species can be spectrosymmetrically characterized in future experiments.     

Molecular structure of the trans and cis isomers of metal-free phthalocyanine studied by gas-phase electron diffraction and high-level quantum chemical calculations: NH tautomerization and calculated vibrational frequencies

The molecular structure of the trans isomer of metal-free phthalocyanine (H2Pc) is determined using the gas electron diffraction (GED) method and high-level quantum chemical calculations. B3LYP calculations employing the basis sets 6-31G**, 6-311++G**, and cc-pVTZ give two tautomeric isomers for the inner H atoms, a trans isomer having D2h symmetry and a cis isomer having C2v symmetry. The trans isomer is calculated to be 41.6 (B3LYP/6-311++G**, zero-point corrected) and 37.3 kJ/mol (B3LYP/cc-pVTZ, not zero-point corrected) more stable than the cis isomer. However, Hartree-Fock (HF) calculations using different basis sets predict that cis is preferred and that trans does not exist as a stable form of the molecule. The equilibrium composition in the gas phase at 471 degrees C (the temperature of the GED experiment) calculated at the B3LYP/6-311++G** level is 99.8% trans and 0.2% cis. This is in very good agreement with the GED data, which indicate that the mole fraction of the cis isomer is close to zero. The transition states for two mechanisms of the NH tautomerization have been characterized. A concerted mechanism where the two H atoms move simultaneously yields a transition state of D2h symmetry and an energy barrier of 95.8 kJ/mol. A two-step mechanism where a trans isomer is converted to a cis isomer, which is converted into another trans isomer, proceeds via two transition states of C(s) symmetry and an energy barrier of 64.2 kJ/mol according to the B3LYP/6-311++G** calculation. The molecular geometry determined from GED is in very good agreement with the geometry obtained from the quantum chemical calculations. Vibrational frequencies, IR, and Raman intensities have been calculated using B3LYP/6-311++G**. These calculations indicate that the molecule is rather flexible with six vibrational frequencies in the range of 20-84 cm(-1) for the trans isomer. The cis isomer might be detected by infrared matrix spectroscopy since the N-H stretching frequencies are very different for the two isomers.     

One-photon mass-analyzed threshold ionization spectroscopy of 2-chloropropene (2-C3H5Cl) and its vibrational assignment based on the density-functional theory calculations

A high-quality mass-analyzed threshold ionization (MATI) spectrum of 2-chloropropene, 2-C3H5Cl, is reported. Its ionization energy determined for the first time from the 0-0 band position was 9.5395+/-0.0006 eV. Almost all the peaks in the MATI spectrum could be vibrationally assigned utilizing the frequencies calculated at the B3LYP6-311++G(3df,3pd) level and the Franck-Condon factors calculated with the molecular parameters obtained at the same level. In particular, the observed methyl torsional progression could be reproduced very well through quantum-mechanical calculations using the molecular parameters obtained at this level. Dramatic lowering of the torsional barrier inferred from the experimental data was entirely compatible with the B3LYP6-311++G(3df,3pd) results. The torsional barrier and the internal rotational constant determined by fits to six torsional peaks were 53.6 and 5.20 cm(-1), respectively. A brief discussion at the level of molecular orbital is presented to account for the dramatic lowering of the torsional barrier upon ionization.     

Protonated acetylene revisited

The potential energy surface for protonated acetylene has been re-examined with large basis sets and highly correlated methods. The energy difference of 3.6–3.8 kcal/mol between the classical structure and non-classical (bridged) structure computed with CCSD (T)/cc-pVQZ, CCSD(T)/6-311+G(3df,2pd), BD(T)/cc- pVQZ, BD(T)/6-311+G(3df,2pd) and CBS-APNO methods is in very good agreement with the best previous calculations, 3.7–4.0 kcal/mol. In contrast, BLYP, B3LYP, PW91, PBE and TPSS density functional methods do rather poorly, yielding −0.52. 0.29, 1.81, 2.16 and 0.62 kcal/mol, respectively, with the 6-311+G(3df,2pd) basis. MP2 calculations predict the classical structure to be a transition state; however, frequency calculations at the CCSD/6-311+G(3df,2pd) level of theory show that the classical structure is a local minimum. CCSD(T), BD(T) and CBS-APNO energy calculations along the MP2/6-311+G(3df,2pd) reaction path indicate that the classical structure is a shallow local minimum separated from the non-classical structure by a very small barrier of 0.11–0.13 kcal/mol. Because the barrier for proton exchange between the non-classical isomers via the classical structure is broad and nearly flat at the top, the tunneling splitting should be reduced, possibly accounting for the 15% difference between the calculated and experimental barrier heights. Contribution to the Fernando Bernardi Memorial Issue.     

Isomerization of B_6,B_6~- and B_6~+ clusters

The interconversions between isomers with the same spin multiplicity of neutral B6 and charged B6-and B6+ clusters have been investigated at the B3LYP/6-311+G level of theory,including determination of the minimum energy pathways with transition states connecting the corresponding reactants and products.In dynamic calculations,26 isomers were optimized,including 11 novel isomers.In order to further refine the energies,single-point B3LYP/6-311+G(3df) calculations were carried out on the corresponding B3LYP/6...     

Theoretical Study on the Structures and Isomerization of Germylenoid GeH2LiCl

Using geometrical optimization and DFT method at the B3LYP/6-311++G (3df,3pd) level, four equilibrium geometries and one transition state of GeH2LiCl were identified, and the structures at the MP2/6-311++G(3df,3pd) level were calculated simultaneously. We also studied the solvent effects on the structures of Germylenoid GeH2LiCl at the B3LYP/6-311++G (3df,3pd) level. The two more stable forms are suggested to be the p-complex and three-membered ring. The vibrational frequencies and infrared intensities were computed at the B3LYP/6-311++G (3df,3pd) level.     

Role of singlet diradicals in reactions of 2-carbenabicyclo[3.2.1]octa-3,6-diene

The generation of 2-carbenabicyclo[3.2.1]octa-3,6-diene (1) results in the formation of C(8)H(8) hydrocarbons endo-6-ethynylbicyclo[3.1.0]hex-2-ene (4), semibullvalene (5), and 5-ethynyl-1,3-cyclohexadiene (6), and C(8)H(10) hydrocarbons bicyclo[3.2.1]octa-2,6-diene (7), tricyclo[3.2.1.0(4,6)]oct-2-ene (8), and tetracyclo[3.3.0.0(2,8)0(4,6)]octane (9). Focus is placed on three mechanistic pathways for the formation of the C(8)H(10) hydrocarbon fraction: (a) abstraction of hydrogen by triplet carbene 1T to produce an equilibrating set of monoradicals, (b) interconversion of triplet carbene 1T into tricyclic triplet diradical 19T and tetracyclic triplet diradical 20T, and (c) interconversion of singlet 1S with analogous singlet diradical 19S and 20S. Ab initio calculations at the (U)B3LYP/6-311+G(3df,2p)//(U)B3LYP/6-31G(d,p) and broken spin symmetry UBS B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d,p) levels rule out choices (a) and (b) and are consistent with the singlet diradical process.     

CH2FCF3与O(1D)反应机理的理论研究

用量子化学密度泛函理论(DFT)和G3(MP2)B3方法,对O(1↑D)与CH2FCF3的反应进行了研究．在UB3LYP/6-31G(d)计算水平上,优化了反应势能面上各驻点的几何结构,在G3(MP2)B3水平上进行了单点计算,并利用UB3LYP/6-311 G(3df,3pd)计算的波函数进行了电荷密度分析．通过内禀坐标(IRC)计算和振动分析,对反应过渡态进行了确认,并确定了反应机理．     

Inverse hydrogen bonds between SiH_4 and hydrides of Na,Mg and Be

The optimized geometries of the three complexes between MeHn (Me=Na,Mg,Be;n=1 or 2) and SiH4 have been calculated at the B3LYP/6-311++g**,MP2/6-311++g(3df,3pd) and MP2/aug-cc-pvtz levels,respectively.The red-shift inverse hydrogen bonds (IHBs) based on Si-H,an electron donor,were reported.The calculated binding energies with basis set super-position error (BSSE) correction of the three complexes are-5.98,-8.65 and-3.96 kJ mol-1 (MP2/6-311++g(3df,3pd)),respectively,which agree with the results obtained via M...     

甲醛光催化降解反应机理的量子化学研究

分别采用B3LYP,MP2方法在6-311++G(2df,pd)水平研究了甲醛光催化降解反应的微观机理,找到了可能的反应通道,预测反应产物为HCOOH与H2O.并得到了各反应通道的反应物、中间体、过渡态和产物的优化构型、谐振频率.成功地解释了实验结论.从键长和能量的变化角度,讨论了化学反应过程中化学键的变化规律,整个反应通道中各势能面均较低,从理论角度分析该反应室温下能够进行,为空气中的甲醛降解反应的实验研究提供理论依据.