A comparative study on the nature and strength of O–O, S–S, and Se–Se bond

A computational study on dichalcogenide molecules (R₂X₂; X = O, S, Se; R = H, CH₃, NH₂) has been carried out employing B3LYP and MP2 levels using 6-31+G*, 6-311+G*, 6-311++G**, and PVDZ basis sets. The relative energies have been evaluated at G2MP2 also. The rotational barriers and bond dissociation energies indicate that S–S bond is stronger than Se–Se and O–O bond. NBO analysis at MP2/6-31+G* suggest the presence of partial π character between X–X bond that decreases in the order S–S > Se–Se > O–O. Fuki functions for nucleophilic and electrophilic attack fail to distinguish the reactivity of S and Se. The proton affinities of the O₂H₂, S₂H₂, Se₂H₂ decrease in the order Se > S > O.     

Quantum-chemical Study on the Structures and Properties of Uracil-BX3 （X = F, Cl） Complexes

1 INTRODUCTION The intermolecular interaction of bases in DNA or RNA is of immense interest and significance to che- mists and biologists alike. The interactions of these bases with metal cations, solvent molecules and other small molecules or ions would affect the struc- ture and biological properties or recognition process,which has been investigated widely[1～8]. Boron contained compounds are electron deficient com- pounds and have been extensively used as catalysts in chemical react…     

Investigation of the use of B3LYP zero-point energies and geometries in the calculation of enthalpies of formation

The use of B3LYP/6–31G^* zero-point energies and geometries in the calculation of enthalpies of formation has been investigated for the enlarged G2 test set of 148 molecules [J. Chem. Phys. 106 (1997) 1063]. A scale factor of 0.96 for the B3LYP zero-point energies gives an average absolute deviation nearly the same as scaled HF/6–31G^* zero-point energies for G2, G2(MP2), and B3LYP/6–311 + G(3df,2p) enthalpies. A scale factor of 0.98, which has been recommended in some studies, increases the average absolute deviation by about 0.2 kcal/mol. Geometries from B3LYP/6–31G^* are found to do as well as MP2/6–31G^* geometries in the calculation of the enthalpies of formation.     

Dimers, trimers and oligomers of sulfur oxides: an ab initio and density functional study

Ab initio (HF/3-21G^*), DFT (B3LYP with basis sets 6-31G^*, 6-311+G^* and 6-311+G(2d)) and, in some cases, MP2/6-31G^* calculations, were done on cyclic dimers, trimers, etc. and on acyclic oligomers (with OH and H on the ends) of sulfur monoxide and sulfur dioxide. The four cyclic (SO)_n molecules were (S–O)₂ (1,3,2,4-dioxadithietane, 1a), (S–O)₃ (1,3,5,2,4,6-trioxatrithiane, 2a), (S(=O))₄ (tetrathietane 1,2,3,4-tetraoxide, 1b), and (S(=O))₆ (hexathiane 1,2,3,4,5,6-hexaoxide, 2b). The four cyclic (SO₂)_n molecules were the dioxide of 1a (1,3,2,4-dioxadithietane 2,4-dioxide, 1c), the trioxide of 2a (1,3,5,2,4,6-trioxatrithiane 2,4,6-trioxide, 2c), the tetraoxide of 1b (tetrathietane 1,1,2,2,3,3,4,4-octaoxide, 1d) and the hexaoxide of 2b (hexathiane 1,1,2,2,3,3,4,4,5,5,6,6-dodecaoxide, 2d). The 16 acyclic molecules (oxides of disulfane, trisulfane, etc. and oxides of oxadisulfane, dioxatrisulfane, etc.) were (–S–O–)_n, (–S(=O)–)_n, (–S(=O)O–)_n, and (–S(=O)₂–)_n, with n from 2 to 5 and HO, H at the ends. Most of these species are relative minima on the B3LYP/6-31G^* potential energy surface. In energy content, the SO dimer, etc. lie below, and the SO₂ dimer, etc. above, their SO_x components, at all the electron-correlated levels.     

Implications of hyperconjugative effects on bond lengths of allylic systems. An NBO investigation

The influence of hyperconjugative interactions on bond lengths of some allylic compounds (H₂CCH–CH₂–M(CH₃)₃; M=C, Si, Ge) has been investigated through NBO calculations using ab initio and density functional methods. The optimized structural parameters, at the B3LYP/6-31+G(d,p) and HF/6-31+G(d,p) levels, showed a good agreement with the resonance theory. Partial geometry optimization with orbital interactions removed confirmed the observations and revealed that σ→σ* interactions, together with the more common σ→π* ones, play an important role in determining the variations in bond lengths on going from C to Ge.     

吡啶-BH~3相互作用复合物的理论研究

对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能，得到四种构型，在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性，然后用B3LYP作振动频率分析，计算了各构型的垂直电离势，最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明，N-B直接相连的构型最稳定，其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大，另外两种构型稳定性介于二者之间，解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Binding energies of germanium clusters, Gen (n = 2−5)

Gaussian-2 (G2) theory for third-row non-transition elements is used to calculate energies of germanium clusters, Ge_n (n = 2−5). The G2 energies are used to derive accurate binding energies for the clusters. The results for Ge₂ and Ge₃ are in agreement with experiment while there is some disagreement for Ge₄ and Ge₅. The binding energies are also calculated using the B3LYP density functional method with the 6–311 + G(3df,2p) basis set and compared with the G2 results and experiment.     

Intramolecular hydrogen bond in the push–pull CF3-aminoenones: DFT and FTIR study,NBO analysis

Postulated conformers of trifluoromethylated β-aminoenones stabilized by intramolecular NH?O and N?HO bonds were studied by IR and NMR spectroscopy and evaluated with quantum chemical calculations (B3LYP/6-311+G(d,p), MP2/6-311+G(d,p)//B3LYP/6-311+G(d,p) and MP2/6-31G(d,p)) and NBO analysis. The influence of the nature of EWG, substituents at the nitrogen atom and double bond, and of orbital interactions of heteroatoms and double bonds in these structures on the proton affinity of basic and acid centers, strength of hydrogen bonds, and the energy of tautomeric transfers is discussed. The theoretical results agree satisfactorily with the experimental observations.     

Molecular structure and conformational stability of allylisocyanate—post-Hartree–Fock and density functional theory studies

The molecular structure and conformational stability of allylisocyanate (CH₂CHCH₂NCO) molecule was studied using the ab initio and DFT methods. The geometries of possible conformers, C-gauche (δ=120°, θ=0°) (δ=C=C–C–N and θ=C–C–N=C) and C-cis N-trans (δ=0° and θ=180°) were optimized employing HF/6-31G^*, MP2/6-31G^* levels of theory of ab initio and BLYP, B3LYP, BPW91 and B3PW91 methods of DFT implementing the atomic basis set 6-311+G(d,p). The structural and physical parameters of the above conformers were discussed with the experimental and theoretical values of the related molecules, methylisocyanate and 3-fluoropropene. It has been found that the N=C=O bond angle is not linear as the experimental result for both the conformers and the theoretical bond angle is 173°. The rotational potential energy surfaces have been performed at the HF/6-31G^*, and MP2/6-31G^* levels of theory. The Fourier decomposition potentials were analysed at the HF/6-31G^*, and MP2/6-31G^* levels of theory. The HF/6-31G^* level of theory predicted that the C-gauche conformer is more stable than the C-cis N-trans conformer by 0.41 kJ/mol, but the MP2 and DFT methods predicted the C-cis N-trans conformer is found to be more stable than the C-gauche conformer. The calculated chemical hardness value at the HF/6-31G^* level of theory predicted the C-cis N-trans form is more stable than C-gauche form, whereas the chemical hardness value at the MP2/6-31G^* level of theory favours the slight preference towards the C-gauge conformer.     

DFT study on the reaction mechanisms of polyfluorosulfonate ester with F

Gas-phase reaction of C(1)F₃S(2)O₂O(3)C(4)H₂C(5)F₃ and F⁻(16) is investigated using DFT method. The geometries of various stationary points and their relative energies are obtained from 6-31+G*, 6-311G**, and 6-311++G** levels. In the S_N2(C) reaction leading to the cleavage of the C(4)–O(3) bond, the reaction complex results from attacking of F⁻ at a hydrogen atom H11 attached to carbon atom C(4). Afterward, F⁻ is attacking the atom C(4) from the backside of the atom O(3) with the help of the neighboring effect, and meanwhile a multi-membered ring, F(16)–H(11)–C(4)–C(5)–F(16), is being formed. The S_N2(C) reaction is irreversible. On the contrary, the S_N2(S) reaction leading to the cleavage of the S(2)–O(3) bond is reversible, and it is initiated by attacking of F⁻ at the atom S(2) from the backside of the atom O(3). The products of the reaction CF₃SO₃CH₂CF₃ +F⁻ should be, thermodynamically, controlled due to the reversibility of the S_N2(S) reaction, and those result, chemospecifically, from the cleavage of the C–O bond. At last, the SCRF calculations confirm that the solvent effect is preferable to the S_N2(C) reaction.     

胸腺嘧啶-BH3构型与性质的DFT和MP2研究

Geometries and combination energies are predicated at B3LYP / 6-31G（d）and MP2 / 6-31G（d）level for thymine-BH3 complexes and 5 geometries have been obtained. Then single point energy calculations using larger basis sets（6-311 + G（2df）and aug-cc-pVDZ）and vibrational analysis and natural bond orbital analysis are carried out on the 5 optimized conformers. The outcome indicates that the conformers with the boron atom combined with O directly are relatively stable ones,（a）and（b）,with the combination energies of 90. 4 and 88. 0 kJ / mol （B3LYP / 6-31G（d）,BSSE corrected）. The fact is that the nitrogen atom offers electron to the empty atomic orbital of boron which produces the conformers（c）and（d）. Only one conformer is found which is formed because two carbon atoms offer π electron to the empty orbital of boron. The charge transference exists in all the conformers. The combination energies have a good line relation with their charge transference. The calculated results show that when the complex forms their IR spectrum moved to the red side and the frequency shifts are relative to the stabilities of the complexes.     

The NICS (Nucleus-Independent Chemical Shift) as a probe of the relative stability of β-chalcogenovinylaldehydes stabilized through intramolecular chalcogen–chalcogen interactions

Density functional calculations, at the B3LYP/6-311+G(3df,2p) level, have been carried out for the complete series of β-chalcogenovinylaldehydes, CH(X)–CHCH–YH (X, Y=O, S, Se, Te), to estimate the strength of H–XY or XY–H intramolecular chalcogen–chalcogen interactions, through the use of appropriate homodesmotic reactions. For the same set of compounds the value of the nucleus-independent chemical shift (NICS), on points 1 Å above the corresponding ring critical point, has been obtained at the B3LYP/6-311+G(3df,2p) level. For non-stabilizing chalcogen–chalcogen interactions the NICS value is positive, while the opposite is found when the interaction is stabilizing. In general, there is a good linear correlation between both magnitudes and therefore, we can conclude that NICS value is a reliable probe of the strength of intramolecular chalcogen–chalcogen interactions in this set of compounds.     

An ab initio and density functional theory study of the dimethylzinc-hydrogen selenide adduct: (CH3)2Zn:SeH2

The molecular structure (equilibrium geometry) and binding energy of the dimethylzinc (DMZn)-hydrogen selenide (H₂Se) adduct, (CH₃)₂Zn:SeH₂, have been computed with ab initio molecular orbital and density functional theory (DFT) methods and, where possible, compared with experimental results. The structure of the precursors DMZn and H₂Se are perturbed to only a small extent upon adduct formation. (CH₃)₂Zn:SeH₂ was found to be 3 kcal mol⁻¹ less stable than the precursors at the B3LYP/6-311 + G(2d,p)//B3LYP/6-311 + G(2d,p) level of computation, indicating that the (CH₃)₂Zn:SeH₂ adduct is unlikely to be a stable gas-phase species under chemical vapour deposition conditions. Further calculations at the B3LYP/6-311 + G(2d,p)//B3LYP/6-311 + G(2d,p) level of computation suggest that the 1:2 adduct species, (CH₃)₂Zn:(SeH₂)₂, is much less stable than the 1:1 adduct and consequently the precursors by 19 kcal mol⁻¹.     

The C7-C10 cycloalkanes revisited

The conformers of cycloheptane through cyclodecane have been examined at the B3LYP/6-311+G* and MP2/6-311+G* theoretical levels, with some additional calculations at the CCD/6-311+G* and CCSD(T)/6-311++G** levels. With cyclooctane, B3LYP predicts that the boat-chair and crown conformers have similar energies, whereas MP2 and CCSD(T) predict that the crown conformer is 2 kcal/mol higher in energy. The latter is in agreement with the electron diffraction data. With cyclononane, B3LYP predicts that two of the higher-energy conformers found in molecular mechanics calculations should convert to one of the lower-energy conformers. However, MP2/6-311+G* optimizations find them to be true minima on the potential energy surface. B3LYP systematically predicts larger C-C-C bond angles for these compounds than either MP2 or CCD. The results of molecular mechanics MM4 calculations are generally in good agreement with those obtained using MP2.     

Infrared and Raman spectra, conformational stability, ab initio calculations, and vibrational assignments for methylaminothiophosphoryl difluoride

Infrared spectra (4000–50 cm⁻¹) of the vapor, amorphous and crystalline solids and Raman spectra (3600–10 cm⁻¹) of the liquid with qualitative depolarization data as well as the amorphous and crystalline solids of methylaminothiophosphoryl difluoride, CH₃N(H)P(=S)F₂, and three deuterated species, CD₃N(H)P(=S)F₂, CH₃N(D)P(=S)F₂, and CD₃N(D)P(=S)F₂, have been recorded. The spectra indicate that in the vapor, liquid and amorphous solid a small amount of a second conformer is present, whereas only one conformer remains in the low temperature crystalline phase. The near-infrared spectra of the vapor confirms the existence of two conformers in the gas phase. Asymmetric top contour simulation of the vapor shows that the trans conformer is the predominant vapor phase conformer. From a temperature study of the Raman spectrum of the liquid the enthalpy difference between the trans and near-cis conformers was determined to be 368±15 cm⁻¹ (4.41±0.2 kJ/mol), with the trans conformer being thermodynamically preferred. Ab Initio calculations with structure optimization using the 6-31G(d) and 6-311+G(d,p) basis sets at the restricted Hartree–Fock (RHF) and/or with full electron correlation by the perturbation method to second order (MP2) support the occurrence of near-trans (5° from trans) and near-cis (20° from cis) conformers. From the RHF/6-31G(d) calculation the near-trans conformer is predicted to be the more stable form by 451 cm⁻¹ (5.35 kJ/mol) and from the MP2/6-311+G(d,p) calculation by 387 cm⁻¹ (4.63 kJ/mol). All of the normal modes of the near-trans rotamer have been assigned based on infrared band contours, depolarization values and group frequencies and the assignment is supported by the normal coordinate calculation utilizing harmonic force constants from the MP2/6-31G(d) ab initio calculations.     

Ab initio and DFT energetics of silylenic X–CNSi (X=H, F, Cl, and Br)

Singlet–triplet energy splitting for 24 silylenic reactive intermediates, X–CNSi (where X=H, F, Cl and Br), are compared and contrasted at 11 levels of theory: B1LYP/6-31++G**, B3LYP/6-31++G**, B1LYP/6-311++G**, B3LYP/6-311++G**, MP3/6-31G*, MP3/6-311++G**, MP2/6-31+G**, MP2/6-311++G**, MP4 (SDTQ)/6-311++G**, QCISD(T)/6-311++G** and CCSD(T)/6-311++G**. Each X-substituted silylenic species may either be singlet (s) or triplet (t), with one of the following three structures: 3-X-2-aza-1-silacyclopropenylidene (1_s-X, 1_t-X); [(X-imino)methylene]silylene (2_s-X, 2_t-X); and X-cyanosilylene (3_s-X, 3_t-X). For all X–CNSi species studied, orders of singlet–triplet energy separations (ΔE_s-t,X), appear as a function of electro-negativity (F>Cl>Br>H). For the six H–CNSi isomers (X=H), stability order is: 3_s-H>1_s-H>2_t-H>3_t-H>2_s-H>1_t-H. Likewise, stability order for the six isomers with X=F, is: 3_s-F>3_t-F>1_s-F>1_t-F>2_s-F>2_t-F. For X=Cl, the order of stability is: 3_s-Cl>1_s-Cl>3_t-Cl>2_t-Cl>1_t-Cl>2_t-Cl. Finally, the order of stability for six isomers of Br–CNSi is: 3_s-Br>3_t-Br>1_s-Br>2_s-Br>2_t-Br>1_t-Br. The lowest energy minimum, among all 24 species scrutinized, appears to be the singlet acyclic 3_s-X. Triplet silylene 2_t-H is suggested to be more stable than its corresponding 2_s-H at MP3, MP2 and DFT levels of theory. Comparisons between relative stabilities; multiplicities and geometrical parameters of 1–3 are discussed.     

Density functional theory analysis of CaRgn complexes: (Rg=He, Ne, Ar; n=1–4)

CaRg_n⁺ (Rg=He, Ne, Ar) complexes with n=1–4, are investigated by performing using the B3LYP/6-311+G (3df) density functional theory calculations. The CaHe_n⁺ (n=1–4) complexes are found to be stable. In the case of CaNe_n⁺ and CaAr_n⁺, stable structures and stationary point were found only for n=1 and 2. For n=3 in the C_3V and the D_3h point group as well as for n=4 in the T_d (tetrahedral) point group a saddle point (imaginary frequency) is observed and global minimum could be obtained along the potential energy surface.     

FTIR studies of the CH3CN–BF3 complex in solid Ar, N2, and Xe: Matrix effects on vibrational spectra

FTIR spectra of the four isotopically substituted 1:1 complexes of acetonitrile and boron trifluoride were recorded in Ar, N₂ and Xe matrices. In Ar, previously unreported three vibrational modes of the complex were clearly observed. Several significant vibrational bands were also observed in N₂ and Xe. The observed frequency shifts on complexation, Δν, were qualitatively in good agreement with the computational results, which were calculated at the B3LYP/6-311++G(d,p) level using the optimized geometry of the C_3v eclipsed conformation. The observed magnitudes of Δν for most of the complex bands were larger than the calculated values. The BF₃ symmetric deformation mode is an exception. The observed frequency shits for this mode were smaller than the calculated values, especially in N₂. This suggests that even an inert matrix can significantly affect the vibrational spectrum of the complex.