An ab initio study of the structure,dissociation energy,and heat of formation of Na2S

The equilibrium geometries and fundamental frequencies of Na₂S are calculated at HF, MP2(FC, FU), and MP3 with the 6–31G(d) basis set and at HF and MP2(FC, FU) with the 6–31G(d) basis set, respectively. The total energy at MP2(FU)/6–31G(d)-optimized geometry is computed at MP4 with 6–311G(d, p), 6–311 + G(d, p), and 6–311G(2df, p), at QCISD(T)/6–311G(d, p), and at MP2/6–311G(3df, 2p) levels, respectively. The dissociation energy, the atomization energy, and the heat of formation for Na₂S are evaluated using the G1 and G2 models. The calculated results indicated that Na₂S in its ground state was a bent structure (C_2v). Electron correlation corrections on the bending angle are very significant. The equilibrium geometrical parameters are R_e(Na-S) = 2.45 Å and ∠Na-S-Na = 111.13° at the MP2(FU)/6–31G(d) level. The theoretically estimated dissociation energy, total atomization energy, and heat of formation are 67.07, 117.55, and 0.35 kcal mol⁻¹, respectively, at 298.15 K. © 1997 John Wiley & Sons, Inc.     

Structure and electronic properties of hollow‐caged C60 fullerene‐derived (MN4)nC6(10 − n) (M = Zn,Mg, Fe,n = 1−6) complexes

Unique hollow‐caged (MN₄)_nC_{6(10 ? n)} (M = Zn, Mg, Fe, n = 1?6) complexes designed by introduction of n porphyrinoid fragments in C₆₀ fullerene structure were proposed and the atomic and electronic structures were calculated using LC‐DFT MPWB95 and M06 potentials and 6‐311G(d)/6‐31G(d) basis sets. The complexes were optimized using various symmetric configurations from the highest O_h to the lowest C₁ point groups in different spin states from S = 0 (singlet) to S = 7 (quindectet) for M = Fe to define energetically preferable atomic and electronic structures. Several metastable complexes were determined and the key role of the metal ions in stabilization of the atomic structure of the complexes was revealed. For Fe₆N₂₄C₂₄, the minimum energy was reported for C_2h, D_2h, and D_4h symmetry of pentet state S = 2, so the complex can be regarded as unique molecular magnet. It was found that the metal partial density of states determine the nature of HOMO and LUMO levels making the clusters promising catalysts. © 2014 Wiley Periodicals, Inc.     

Theoretical studies on a possible synthesis reaction pathway on N8 (CS) clusters

The potential energy surfaces of N₈ clusters were investigated by density functional theory (DFT) and a possible synthesis reaction pathway for N₈ (C_S) was suggested. The species involved were fully optimized up to the B3LYP/6‐311+G* level of theory. Relative energies were further calculated at the QCISD/6‐311+G*//B3LYP/6‐311+G* level. The reaction rate constants of these steps from the 1 (N₅⁺?N₃^?, complex, C_S) to 2 (N₈, C_S), 2 (N₈, C_S) to 3 (N₈, C_S), 3 (N₈, C_S) to 4 (N₈, D_2d), and 4 (N₈, D_2d) to 5 (N₈, C_S) reactions were predicted by the VTST theory. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1334–1339, 2001     

The N4 molecule has an open-chain triplet C2h structure

The lowest-energy N₄ is computed ab initio to be the planar C_2h(³B_u) open-chain structure 13 . The open-chain N₄ singlet-state structures dissociate on geometry optimization. The tetraazatetrahedrane T_d structure 1 and the tetrazete D_2h structure 2 are minima at MP 2/6-31G *. However, both are higher in energy than 13 (24.1 and 21.2 Kcal/mol [UQCISD ) (T )(full)/6-311+G *//MP 2/6-31G * + ZPE (MP 2/6-31G )*, respectively]. The energy of 13 is 157.5 kcal/mol higher than that of two N₂(¹∑ molecules [UQCISD (T )(Full)/6-311+G *//MP 2/6-31G *] © 1993 John Wiley & Sons, Inc.     

Why betaine crystallizes in high local Cs symmetry. An ab initio MO and DFT study of anhydrous betaine and betaine monohydrate

A theoretical study of the structure, charge distribution, rotational barrier and fundamental vibrations of anhydrous betaine (CH₃)₃NCH₂COO (trimethylglycine) was carried out and compared with available experimental data. Calculations were carried out at HF, MP2 and B3LYP levels using a 6-31+G(d,p) basis set. The calculated rotational barrier of the betaine carboxylic group is 40.5 kJ/mol at the MP4(SDQ)/6-311G(d,p)//HF/6-31+G(d,p) level of theory. The rotation of the carboxylic group changes the molecule from a highly symmetric (C _s ) conformation into a twisted conformation resulting in shortening of the molecule by about 50 pm. Natural population analysis (NPA) indicates intramolecular interaction between the carboxylic oxygen and the nearest methyl hydrogens resulting in internal hydrogen bonding. MP4(SDQ)/6-311G(d,p) single-point NPA calculations on a betaine monohydrate model taken from the X-ray geometry show an expected weakening in the internal hydrogen bond. Calculations explain why betaine preferentially crystallizes in high local C _s symmetry. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 7 December 1998     

Is tetrahedral H42+ a minimum? Anomalous behavior of popular basis sets with the standard p exponents on hydrogen

The nature of the tetrahedral H₄²⁺ stationary point (minimum or triply degenerate saddle) depends remarkably upon the theoretical level employed. Harmonic vibrational analyses with, e.g., the 6-31G** (and 6-31 + +G**) and Dunning's [4s2p1d;2s1p] [D95(d,p)] basis sets using the standard p exponent suggest (erroneously) that the T_d geometry is a minimum at both the HF and MP2 levels. This is not the case at definitive higher levels. The C₃H₄²⁺ structure with an apical H is another example of the failure of the calculations with the 6-31G**, 6-311G**, and D95(d,p) basis sets. Even at MP2/6-31G** and MP2/ cc-pVDZ levels, the C_3v structure has no negative eigenvalues of the Hessian. Actually, this form is a second-order saddle point as shown by the MP2/6-31G** calculation with the optimized exponent. The D_4h methane dication structure is also an example of the misleading performance of the 6-31G** basis set. In all these cases, energy-optimized hydrogen p exponents give the correct results, i.e., those found with more extended treatments. Optimized values of the hydrogen polarization function exponents eliminate these defects in 6-31G** calculations. Species with higher coordinate hydrogens may also be calculated reliably by using more than one set of p functions on hydrogen [e.g., the 6-31G(d,2p) basis set]. Not all cases are critical. A survey of examples, also including some boron compounds, provides calibration. © 1993 John Wiley & Sons, Inc.     

On the additivity of basis set effects in some simple fluorine containing systems

The reactions F + H₂ → HF + H, HF → H + F, F → F⁺ + e^? and F + e^? → F^? were used as simple test cases to assess the additivity of basis set effects on reaction energetics computed at the MP4 level. The 6-31G and 6-311G basis sets were augmented with 1, 2, and 3 sets of polarization functions, higher angular momentum polarization functions, and diffuse functions (27 basis sets from 6-31Gd, p) to 6-31 ++ G(3df, 3pd) and likewise for the 6-311G series). For both series substantial nonadditivity was found between diffuse functions on the heavy atom and multiple polarization functions (e.g., 6-31 + G(3d, 3p) vs. 6-31 + G(d, p) and 6-31G(3d, 3p)). For the 6-311G series there is an extra nonadditivity between d functions on hydrogen and multiple polarization functions. Provided that these interactions are taken into account, the remaining basis set effects are additive to within ±0.5 kcal/mol for the reactions considered. Large basis set MP4 calculations can also be estimated to within ±0.5 kcal/mol using MP2 calculations, est. E_MP4(6-31 ++ G(3df, 3pd)) ≈ E_MP4(6-31G(d, p)) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31G(d, p)) or E_MP4(6-31 + G(d, p) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31 + G(d, p)) and likewise for the 6-311G series.     

Theoretical Study of the Structures,Properties and Spectroscopies on Fullerene Hydrides C<Subscript>26</Subscript>H<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> (<Emphasis Type="Italic">n</Emphasis> = 2, 4, 6, 8)

The structures, stability patterns of C₂₆H _n (n = 2) formed from the initial D _3h C₂₆ fullerene were investigated by use of second-order-Moller–Plesset perturbation theory. The study of the stability patterns of hydrogenation reaction on C₂₆ cage revealed that type (β) carbons were the active site and the analyses of π-orbital axis vector indicated that the reactivity of C₂₆ was the result of the high strain and the hydrogenation reaction on C₂₆ cage was highly exothermic. The calculated ¹³C NMR spectra of C₂₆H _n (n = 2) predicted that the two sp ³ hybridization carbons in C₂₆H _n (n = 2) obviously moved to high field compare with that in D _3h C₂₆. Hence, the C₂₆H₂ should be obtained and detected experimentally. Similarly, the structures and reaction energies of C₂₆H _n (n = 4, 6, 8) were further studied at HF/6-31G*, B3LPY/6-31G* and MP2/6-31G* level. The results suggested the hydrogenation products of C₂₆, C₂₆H _n (n = 4, 6, 8), were more stable than the C₂₆ cage.     

A Theoretical Study of LiHe+ n , NaHe+ n , and MgHe+ n Complexes, with n=1, 2, 3, 4

The LiHe⁺ _n , the NaHe⁺ _n , and the MgHe⁺ _n complexes with n=1, 2, 3, 4 were studied using ab initio calculations with the MP2/6-311+G(3df, 3pd) method. The complexes are found to be stable. For the n=1 complexes, previous results were available and the calculations performed are in good agreement with those results. This lends credibility to the results obtained for the complexes with higher n.     

Exploration on the structure,stability, and isomerization of planar CnB5 (n = 1−7) clusters

The structures, stabilities, nature of bonding, and potential energy surfaces of low‐energy isomers of planar C_nB₅ (n = 1?7) have been systematically explored at the CCSD(T)/6‐311+G(d)//B3LYP/6‐311+G(d) level. Incremental binding energy (IBE) and second order energy difference (Δ²E) analyses demonstrate that C_nB₅ clusters with even n have relatively higher stability. The nature of bonding in these clusters is discussed based on valence molecular orbital (VMO), and Mayer bond order (MBO). Hückel (4n + 2) rule and nucleus‐independent chemical shift (NICS) values suggest that the ground states of C₃B₅, C₄B₅, and C₇B₅ have π aromaticity. VMO, electron localization function (ELF), adaptive natural density partitioning (AdNDP), and NICS analyses reveal the double aromaticity of C₃B₅ cation. CB₅ and C₃B₅ are stable both thermodynamically and kinetically based on isomerization analysis. In addition, the simulated IR spectra are expected to be helpful for future experimental studies of these clusters.     

HNO(H<Subscript>2</Subscript>O)<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> (<Emphasis Type="Italic">n</Emphasis> = 1–4) clusters: a theoretical study

The geometrical structure, binding energy, and vibrational spectra of small clusters of nitrosyl hydride (HNO) and water molecules, HNO(H₂O) _n , where (n = 1–4), have been investigated at the MP2 level of theory, using 6-311++G(2d,2p) basis set. We located three dimers, six trimers, nine tetramers, and three pentamers at the MP2/6-311++G(2d,2p) computational level. Particular attention is given to existence and magnitude of NH···O blue-shifting hydrogen bonds. Blue shifts of the NH stretching frequency upon complex formation in the ranges between 28 and 151 cm⁻¹ is predicted. Cooperative effect in terms of stabilization energy along with the many-body interaction energies analysis was performed for the studied clusters. The Atoms in Molecules (AIM) theory was also applied to explain the nature of the complexes.     

Density functional calculations of the pseudorotational flexibility of tetrahydrofuran

The flexibility of the five-membered ring in tetrahydrofuran was investigated using quantum mechanical methods involving density functional, Hartree-Fock, and many-body perturbation theory (MP2, MP4) calculations. We found that motion along the pseudorotational path of tetrahydrofuran is nearly free. The 0.1 kcal/mol energy barrier for pseudorotation, calculated at the highest MP4(SDQ)/6-311++G(2d,p)//MP2/6-311++G(2d,p) level of theory, agrees well with the experimental value of 0.16±0.03 kcal/mol. Similar results were obtained with the S-VWN, B3-LYP and B-LYP density functional calculations using the 6-31G(d) set of atomic orbitals. Also the density functional dipole moments and geometries were in good agreement with both the MP2 and experimental benchmarks. However, all density functional methods that utilized the default integration grid in the Gaussian 94 program were found to provide false stationary points of the C ₁ symmetry near the pseudorotational angle of 100°. These stationary points disappeared when a denser spherical-product grid was used. Overall, the hybrid B3-LYP functional was found to be the most promising quantum mechanical method for the modeling of biomolecules containing the furanose ring. Received: 17 June 1997 / Accepted: 20 November 1997     

Does CH prefer a C2v rather than a Cs structure?

The closely related C_s ( 1 ) and C_2v ( 3 ) structures of CH have been reinvestigated at many ab initio levels using MP2/6-31G** and MP2/6-311 + + G(2df, 2pd) geometries. The largest basis sets employed were 6-311G(3df, 2p), 6-311 + + G(3df, 3pd), and the Dunning “correlation consistent” polarized triple-split valence basis set (cc-pVTZ). Electron correlation was probed at the MP4 level, but the QCISD method was also used with the largest basis sets. While electron correlation favors 3 over 1 by about 2 kcal/mol, the correlated relative energies with all basis sets employed range from 0.36–1.03 kcal/mol in favor of 1 . The best estimate of this difference, 0.86 kcal/mol, is essentially identical with the (scaled) zero-point energy difference, 0.84 kcal/mol, favoring 3 over 1 . These results indicate that 1 and 3 have almost exactly the same energy at 0 K. Our best value for the dissociation energy of CH is 42.0 kcal/mol [QCISD(T)/6-311 + + G(3df, 3pd)//MP2(fu)/6-311 + + G(2df, 2pd), corrected to 298 K], which agrees very well with the experimental value. © 1992 by John Wiley & Sons, Inc.     

The reaction of nitroso oxides with olefins: Concerted or nonconcerted addition?

The mechanism of the interaction of nitroso oxides (RNOO) with olefins was studied at MCQDPT2/6-311+G(3df, 2p)//CASSCF(10; 9)/6-311G(d) level of theory. The following reaction channels were considered: (1) (3 + 2)-cycloaddition and nonconcerted biradical addition of nitroso oxide (2) through the terminal oxygen atom and (3) through the nitrogen atom to the C=C multiple bond. It was shown for the cases of (A) cis/trans-HNOO + C₂H₄, (B) cis/trans-HNOO + C₂F₄, (C) cis/trans-PhNOO + C₂H₄, and (D) cis/trans-PhNOO + C₂H₃CH₃ model systems that the typical reaction of nitroso oxides with alkenes was cycloaddition. For olefins with a decreased electron density at the multiple bond, as in system B, a substantial contribution of the one-center mechanism with the formation of biradical intermediates is possible.     

Post-Hartree-Fock and DFT level studies on the Cl2CO …︁ Cl2 complex: Accurate molecular parameters,harmonic vibrational frequencies,and interaction energies

The Cl₂CO …︁ Cl₂ complex was studied using ab initio post-Hartree-Fock theory at the MP2 and MP4 levels and, for comparison, the DFT method with 6-311G(2d), 6-311 + G(2d), and Sadlej's medium-size polarized (MSPBS) basis sets. A potential energy search recovered a planar minimum-energy structure characterized by a bent conformation. For this weakly bound complex, the interaction energy corrected for the basis set superposition error amounted to − 0.88, − 1.09, − 1.43, and − 0.38 kcal/mol at the MP4(SDTQ)/6-311G(2d), MP4(SDTQ)/6-311 + G(2d), MP4(SDTQ)/MSPBS, and DFT(Becke3LYP)/6-311 + G(2d) levels of theory, respectively. Two highly symmetrical forms, linear and T-shaped, correspond to transition-state conformers. The analysis of harmonic vibrational frequencies and potential energy distribution was performed at the MP2 and DFT levels with the 6-311 + G(2d) basis set. © 1996 John Wiley & Sons, Inc.     

Theoretical study of structure and vibrational properties of MgnOn (n = 3–10) clusters

The structure and harmonic vibrations of Mg_nO_n (n = 3–10) clusters have been investigated using density functional theory. All structures are found to be cumulenic D_nh rings (equal bonds, alternating angles), with one intense out‐of‐plane mode and three infrared (IR)‐active degenerate modes, of which the highest one is extremely intense and increases asymptotically to 1000 cm^?1 for n = 10 at the B3LYP/6‐311++G(2d,2p) level. Comparisons with C_2n clusters show that B_nN_n and Be_nO_n clusters, the structure and bonding type for the Mg_nO_n clusters are consistent with those of the C_2n (n = 3, 5, 7,…) clusters B_nN_n(n = 3–10) and Be_nO_n(n = 3–10) clusters. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Isomerization study of C5H5NO molecules

The complex potential energy surface (PES) for the isomerization of C₅H₅NO species, including 18 isomers and 23 interconversion transition states, is probed theoretically at the B3LYP/6‐311++G(d,p) and MP2//B3LYP/6‐311++G(d,p) levels of theory. The geometries and relative energies for various stationary points were determined. The zero‐point vibrational energy (ZPVE) corrections have been made to calculate the reliable energy. We predicted a six‐membered ring structure as a global minima isomer I, which is 118.49 and 131.48 kcal · mol^?1 more stable than the least stable, four‐ and three‐membered ring isomer VIII at B3LYP and MP2//B3LYP levels of theory, respectively. The isomers and interconversion transition states have verified by frequency calculation. The intrinsic reaction coordinates (IRC) calculations have been performed to confirm that each transition state is linked by the desired reactants and products. The isomer stability has been studied using relative energies, chemical hardness, and chemical potential. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical design of new sandwich compounds of boron,carbon, nitrogen,and oxygen

Energetic and structural characteristics of new sandwich compounds of the second row elements X(Y₃)₂ and X(Y₄)₂ (X = B^5?, C^4?, N^3?, O^2?, Y = C, Si) stabilized with lithium cations were investigated using the methods of density functional theory B3LYP/6-311+G(d,p) and ab initio MP2(full)/6-311+G(d,p). Carbon sandwich systems were found to be thermodynamically more stable than their silicon analogs. In all systems the formation of hypercoordinated central atom is observed. For the most representative structures the Bader topological analysis of electron density was performed. The key types of distortions in the silicon-containing systems were found.     

HCN exchange on [Cu(HCN)4]+: a quantum chemical investigation

Density functional (B3LYP, B3PW91, X3LYP, BP86, PBEPBE, PW91PW91, and M06) and ab initio (MP2, MP4sdq, CCSD, and CCSD(T)) calculations with extended basis sets (6-311+G**, TZVP, LANL2DZ+p, and SDD+p, the latter including extra polarization and diffuse functions) indicate that HCN exchange on [Cu(HCN)₄]⁺ proceeds via an associative interchange (I_a) mechanism and a D_3h transition structure {[Cu(HCN)₅]⁺}^?. The activation barrier, relative to the model complex [Cu(HCN)₄]⁺·HCN, varies modestly, depending on the computational level. Typical values are 8.0?kcal?M^?1 (B3LYP/6-311+G**), 6.0?kcal?M^?1 (M06/6-311+G**), and 4.8?kcal?M^?1 (CCSD(T)/6-311+G**//MP2(full)/6-311+G**). Inclusion of an implicit solvent model (B3LYP(CPCM)/6-311+G**) leads to an activation barrier of 5.8?kcal?mol^?1. Comparison of the HCN exchange mechanisms on [Li(HCN)₄]⁺ (limiting associative, A) and [Cu(HCN)₄]⁺ (associative interchange, I_a) reveals that π back donation in the equatorial Cu–N bonds in the transition state determines the mechanism.     

Tetra- and hexaatomic cyclic clusters of main-group elements (XY)2 and (XY)3: an ab initio and density functional study

The electronic and molecular structure of planar (cyclic and linear) tetra- and hexaatomic clusters (XY) _n (XY = CC, BN, BeO, LiF; n = 2, 3) was studied using the ab initio CCD(full)/6-311+G** method and density functional approach (B3LYP/6-311+G**). The stability of cyclic clusters C₆, B₃N₃, and Be₃O₃ with D_3h symmetry is mainly determined by the aromaticity of their -electron systems.