Investigation on stabilities and spectroscopy of C80O2 based on C80 (D5d) using density functional theory

The relative stabilities of the 17 possible isomers for C₈₀O₂ based on C₈₀ (D_5d) were studied using Becke three parameters plus Lee, Yang, and Parr's (B3LYP) method and 6‐31G (d) basis set in density functional theory. The most stable geometry of C₈₀O₂ was predicted to be 23,24,27,28‐C₈₀O₂ (A) with annulene‐like structures, where the additive bonds are those between two hexagons (6/6 bonds) near the equatorial belt of C₈₀ (D_5d). Electronic spectra of C₈₀O₂ isomers were calculated based on the optimized geometries using intermediate neglect of differential overlap (INDO) calculation. Compared with those of C₈₀ (D_5d), the first absorptions in the electronic spectra of C₈₀O₂ are blue‐shifted owing to the wide energy gaps. ¹³C nuclear magnetic resonance spectra and nucleus independent chemical shifts of the C₈₀O₂ isomers were computed at B3LYP/6‐31G level. The chemical shifts of the bridged carbon atoms in the epoxy structures of C₈₀O₂ compared with those of the bridged carbon atoms in the annulene‐like structures are changed upfield. Generally, the isomers with the annulene‐like structures of C₈₀O₂ are more aromatic than those with the epoxy structures. The addition of the oxygen atoms near the pole of C₈₀ (D_5d) is favorable to improving the aromaticities of C₈₀O₂. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

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.     

Ab initio study on the reaction mechanism of ozone with bromine atom

Ab initio calculations of the potential energy surface (PES) for the Br+O₃ reaction have been performed using the MP2, CCSD(T), and QCISD(T) methods with 6‐31G(d), 6‐311G(d), and 6‐311+G(3df). The reaction begins with a transition state (TS) when the Br atom attacks a terminal oxygen of ozone, producing an intermediate, the bromine trioxide (M), which immediately dissociates to BrO+O₂. The geometry optimizations of the reactants, products, and intermediate and transition states are carried out at the MP2/6‐31G(d) level. The reaction potential barrier is 3.09 kcal/mol at the CCSD(T)/6‐311+G(3df)//MP2 level, which shows that the bromine atom trends intensively to react with the ozone. The comparison of the Br+O₃ reaction with the F+O₃ and Cl+O₃ reactions indicates that the reactions of ozone with the halogen atoms have the similar reaction mechanism. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Long‐range π‐type hydrogen bond in dimers CH2?HF,CH2O?H2O,and CH2O?NH3

Using four basis bets, (6‐311G(d,p), 6‐31+G(d,p), 6‐31++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for the dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. In contrast with the above three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD (T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Steric and Electronic Structure of Selenoanisole: A Quantum-Chemical Study

The potential functions of internal rotation around the C_sp2-Se bond in selenoanisole were ob- tained by quantum-chemical calculations in the approximations HF/3-21G(d), HF/6-31G(d), MP2(f)/6-31G(d), and B3LYP/6-31G(d). The calculations were performed in the range of variation of the torsion angle (between the planes of the benzene ring and C_sp2-Se-C_sp3 bonds) from 0° to 90° with 15° step. The energy minimum is in the region of the orthogonal conformation ( 90°), and the energy maximum, in the region of the planar form ( 0°). The rotation barriers (kJ mol^{- 1}) are as follows: HF/3-21G(d), 9.20; HF/6-31G(d), 13.13; MP2(f)/6-31G(d), 10.25; and B3LYP/6-31G(d), 6.41. The geometric parameters, Koopmans ionization potentials, and dipole moments are given. The energies, degrees of hybridization, populations of the lone electron pairs of Se, energies of their interaction with the antibonding * orbitals of the benzene ring, and electron density distributions were determined in terms of the natural bond orbital approach.     

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     

Rearrangements in Model Peptide‐Type Radicals via Intramolecular Hydrogen‐Atom Transfer

Intramolecular H‐atom transfer in model peptide‐type radicals was investigated with high‐level quantum‐chemistry calculations. Examination of 1,2‐, 1,3‐, 1,5‐, and 1,6[C ? N]‐H shifts, 1,4‐ and 1,7[C ? C]‐H shifts, and 1,4[N ? N]‐H shifts (Scheme 1), was carried out with a number of theoretical methods. In the first place, the performance of UB3‐LYP (with the 6‐31G(d), 6‐31G(2df,p), and 6‐311+G(d,p) basis sets) and UMP2 (with the 6‐31G(d) basis set) was assessed for the determination of radical geometries. We found that there is only a small basis‐set dependence for the UB3‐LYP structures, and geometries optimized with UB3‐LYP/6‐31G(d) are generally sufficient for use in conjunction with high‐level composite methods in the determination of improved H‐transfer thermochemistry. Methods assessed in this regard include the high‐level composite methods, G3(MP2)‐RAD, CBS‐QB3, and G3//B3‐LYP, as well as the density‐functional methods B3‐LYP, MPWB1K, and BMK in association with the 6‐31+G(d,p) and 6‐311++G(3df,3pd) basis sets. The high‐level methods give results that are close to one another, while the recently developed functionals MPWB1K and BMK provide cost‐effective alternatives. For the systems considered, the transformation of an N‐centered radical to a C‐centered radical is always exothermic (by 25 kJ ? mol^?1 or more), and this can lead to quite modest barrier heights of less than 60 kJ ? mol^?1 (specifically for 1,5[C ? N]‐H and 1,6[C ? N]‐H shifts). H‐Migration barriers appear to decrease as the ring size in the transition structure (TS) increases, with a lowering of the barrier being found, for example when moving from a rearrangement proceeding via a four‐membered‐ring TS (e.g., the 1,3[C ? N]‐H shift, CH₃? C(O)? NH^. → ^.CH₂? C(O)? NH₂) to a rearrangement proceeding via a six‐membered‐ring TS (e.g., the 1,5[C ? N]‐H shift, ^.NH? CH₂? C(O)? NH? CH₃ → NH₂? CH₂? C(O)? NH? CH₂^.).     

Quantum-Chemical Study on Positional Selectivity in the Trimethylsilylation and Sulfonation of Pyrrole and <Emphasis Type="Italic">N</Emphasis>-Alkylpyrroles

Positional selectivity (α:β ratio) of electrophilic substitution in pyrrole, N-methylpyrrole, and N-tert-butylpyrrole was analyzed by ab initio [RHF/6-31G(d), MP2/6-31G(d)//RHF/6-31G(d)] and DFT [B3LYP/6-31G(d)] calculations of some σ-complexes derived from the substrates. The results of calculations with the use as model electrophilic species of trimethylsilyl cation [MP2/6-31G(d)//RHF/6-31G(d) and B3LYP/6-31G(d)] and SO₃ molecule [B3LYP/6-31G(d)] instead of proton are fairly consistent with the experimental data, according to which trimethylsilylation of pyrrole and its N-substituted derivatives with trimethylsilyl trifluoromethanesulfonate, as well as sulfonation with pyridine-sulfur trioxide complex, gives the corresponding β-substituted products.     

13C and 15N NMR spectra of aminobenzimidazoles in solution and in the solid state

The ¹³C [hexadeutero‐dimethylsulfoxide (DMSO‐d₆), hexamethyl‐phosphoramide (HMPA)‐d₁₈and solid‐state] and ¹⁵N (solid‐state) NMR spectra of six C‐aminobenzimidazoles have been recorded. The tautomerism of 4(7)‐aminobenzimidazoles and 5(6)‐aminobenzimidazoles has been determined and compared with B3LYP/6‐311 + + G(d,p) calculations confirming the clear predominance of the 4‐amino tautomer and the slight preference for the 6‐amino tautomer. GIAO‐calculated absolute shieldings compare well with experimental chemical shifts. Copyright © 2008 John Wiley & Sons, Ltd.     

An ab initio molecular orbital study of the structures and energies of neutral and charged bimolecular complexes of NH3 with the hydrides AHn (A = N,O, F,P, S,and Cl)

Hartree-Fock 6-31G(d) structures for the neutral, positive ion, and negative ion bimolecular complexes of NH₃ with the first- and second-row hydrides AH_n (AH_n = NH₃, OH₂, FH, PH₃, SH₂, and ClH) have been determined. All of the stable neutral complexes except (NH₃)₂, the positive ion complexes with NH₃ as the proton acceptor, and the negative ion complexes containing first-row anions exhibit conventional hydrogen bonded structures with essentially linear hydrogen bonds and directed lone pairs of electrons. The positive ion complex NH₄⁺ …? OH₂ has the dipole moment vector of H₂O instead of a lone pair directed along the intermolecular line, while the complexes of NH₄⁺ with SH₂, FH, and ClH have structures intermediate between the lone-pair directed and dipole directed forms. The negative ion complexes containing second-row anions have nonlinear hydrogen bonds. The addition of diffuse functions on nonhydrogen atoms to the valence double-split plus polarization 6-31G(d,p) basis set usually decreases the computed stabilization energies of these complexes. Splitting d polarization functions usually destabilizes these complexes, whereas splitting p polarization functions either has no effect or leads to stabilization. The overall effect of augmenting the 6-31G(d,p) basis set with diffuse functions on nonhydrogen atoms and two sets of polarization functions is to lower computed stabilization energies. Electron correlation stabilizes all of these complexes. The second-order Møller–Plesset correlation term is the largest term and always has a stabilizing effect, whereas the third and fourth-order terms are smaller and often of opposite sign. The recommended level of theory for computing the stabilization energies of these complexes is MP2/6-31+G(2d,2p), although MP2/6-31+G(d,p) is appropriate for the negative ion complexes.     

Theoretical studies of the gas‐phase pyrolysis kinetics of ω‐bromonitriles,ZCH2CH2Br [Z = NC,NCCH2, NCCH2CH2]

The gas‐phase elimination of several ω‐bromonitriles (ZCH₂CH₂Br, Z = NC, NCCH₂, NCCH₂CH₂) has been examined at the MP2/6‐31G(d,p), MP2/6‐31G(2d,2p), B3PW91/6‐31G(2d,2p), and MPW91PW91/6‐31(2p,2d) levels of theory. The bromonitriles yield the corresponding cyano‐olefin and HBr gas in a rate‐determining step. The MPW91PW91/6‐31G(2p,2d) results suggest a concerted mechanism, with a polar, four‐membered transition state. The calculated kinetic and thermodynamic parameters were found to be within reasonable agreement with the experimental determinations. Substituent effects are discussed in terms of electronic structure. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 41: 168–175, 2009     

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.     

Highly symmetrical phthalocyanines and perfluorophthalocyanines: The quantum-chemical study

By the quantum-chemical method (U)B3LYP/6-31G(d ₅,p) are determined point symmetry group (D _4h ) and equilibrium structure of phthalocyanine (PcH₂), phthalocyaninates PcBe, PcMg, PcCo, PcNi, PcCu, PcZn, perfluorophthalocyaninates ^FPcNi, ^FPcCu, ^FPcZn, cations Pc⁺Mt, ^FPc⁺Mt and anion PcCo⁻. In the approximation (U)B3LYP/6-311++G(3d ₅ f ₇,p)//6-31G(d ₅,p) is achieved the satisfactory accuracy of the calculation of ionization potentials of the studied molecules. Effect of nuclear relaxation at the ionization is 0.07±0.06 eV; correction for “zero” vibrations does not exceed 0.01 eV. Perfluorination increases ionization potentials by 0.7–0.8 eV.     

Ethane cation decomposition characterization by EBMI spectroscopy: gas‐phase dissociative recombination as a source of secondary products

The decomposition products of the d₆‐ethane cation following charge‐transfer ionization with Ar⁺, under conditions of varying ionization electron current, have been isolated in solid argon matrices at 18 K and examined using Fourier transform infrared spectroscopy. Gas samples containing 1 : 1600 d₆‐ethane : Ar were subjected to electron bombardment by using either a high (pin) or a low (plate) ionization density anode configuration with ionization currents between 20 and 150 μA. Under high ionization density conditions, the observed major products were d₄‐ethene (C₂D₄) and d₂‐acetylene (C₂D₂), with smaller yields of C₂D₅, C₂D₃, and C₂D. The yield of each dehydrogenation product was enhanced with increased current. Analogous experiments employing the low ionization density plate anode resulted in reduced C₂D₆ destruction and the formation of only C₂D₄ and C₂D₂. The results suggest the onset of dissociative recombination processes under high ion density conditions. In this context, the results can be interpreted as a dissociative recombination of primary ion products, which gives rise to further dehydrogenation, and appearance of additional neutral radical products. Copyright © 2012 John Wiley & Sons, Ltd.     

Intra- and intermolecular hydrogen bonds in ethylene glycol,monoethanolamine, and ethylenediamine

The structures of ethylene glycol, aminoethanol, ethylenediamine, and their dimers with the formation of hydrogen bonds of different types are optimized by density functional theory (DFT) using hybrid functional B3LYP in the basis of 6-31++G(d,p), 6-311++G(2d,2p) and aug-CC-pVTZ. Energies of interactions, hydrogen bond parameters, and oscillation frequency are calculated, and NBO analysis is performed. The types of hydrogen bonds formed in dimers of 1,2-disubstituted ethanes X-CH₂-CH₂-Y (X, Y = OH, NH₂) are established.     

Quantum‐chemical,IR, NMR,and X‐ray diffraction studies on 2‐(4‐chlorophenyl)‐1‐methyl‐ 1H‐benzo[d]imidazole

The title molecule, 2‐(4‐chlorophenyl)‐1‐methyl‐1H‐benzo[d]imidazole (C₁₄H₁₁ClN₂), was prepared and characterized by ¹H NMR, ¹³C NMR, IR, and single‐crystal X‐ray diffraction. The molecular geometry, vibrational frequencies, and gauge including atomic orbital (GIAO) ¹H and ¹³C NMR chemical shift values of the title compound in the ground state have been calculated by using the Hartree‐Fock (HF) and density functional theory (DFT/B3LYP) method with 6‐31G(d) basis sets, and compared with the experimental data. The calculated results show that the optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies and GIAO ¹H and ¹³C NMR chemical shifts show good agreement with experimental values. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6‐31G(d) basis set by applying the Onsager and the polarizable continuum model (PCM). Besides, molecular electrostatic potential (MEP), frontier molecular orbitals (FMO) analysis, and nonlinear optical (NLO) properties of the title compound were investigated by theoretical calculations. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

Comparative studies on the structures,infrared spectrum,and thermodynamic properties of phthalocyanine using ab initio Hartree–Fock and density functional theory methods

The molecular structure, vibrational spectrum, standard thermodynamic functions, and enthalpy of formation of free base phthalocyanine (Pc) have been studied using the density functional theory B3LYP procedure, as well as the ab initio Hartree–Fock method. Various basis sets 3‐21G, 6‐31G*, and LANL2DZ have been employed. The results obtained at various levels are discussed and compared with each other and with the available experimental data. It is shown that calculations performed at the Hartree–Fock level cannot produce a reliable geometry and related properties such as the dipole moment of Pc and similar porphyrin‐based systems. Electron correlation must be included in the calculations. The basis set has comparatively less effect on the calculated results. The results derived at the B3LYP level using the smaller 3‐21G and LANL2DZ basis sets are very close to those produced using the medium 6‐31G* basis set. The geometry of Pc obtained at the B3LYP level has D_2h symmetry and the diameter of the central macrocycle is about 4 Å. The enthalpy of formation of Pc in the gas phase has been predicted to be 1518.50 kJ/mol at the B3LYP/6‐311G(2d,2p)//B3LYP/6‐31G* level via an isodesmic reaction. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Correlated one‐electron wave functions

Using four basis sets, 6‐311G(d,p), 6‐31+G(d,p), 6‐311++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the acidic H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. By contrast with above the three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen‐bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is a noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD(T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005