An ab initio study of the structures and energetics of the planar ground and 90°‐twisted excited states of substituted ethylenes

Ab initio molecular orbital calculations are performed on the planar ground states (S₀), the 90°‐twisted triplet (T₁), and pyramidalized singlet (S₁) excited states of ethylene, methaniminium cation (MC), monocyano‐ (MCE), 1,1‐dicyano‐ (DCE), 1,1‐dihydroxy‐ (DHE), and 1,1‐dicyano‐2,2‐dihydroxy (DCHE) ethylenes. Equilibrium geometries are optimized at the Hartree–Fock (HF) level with the 6‐31G* basis set. Electron correlation corrections are estimated by optimizing the HF/6‐31G* geometries at the (U)MP2/6‐31G* level and then by carrying out single‐point calculations at the fourth‐order Møller–Plesset perturbation theory ((U)MP4/6‐311G**//MP2/6‐31G*). The effects of various types of perturbations on the structures, energetics, dipole moments, and state ordering of S₀, S₁, and T₁ are carefully investigated. “Positive” S₁–T₁ splittings are estimated at the HF level for all studied molecules, while “negative” S₁–T₁ splittings are obtained at the MP2 level for MC, DHE, and DCHE. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 242–254, 2001     

BSSE‐corrected geometry and harmonic and anharmonic vibrational frequencies of formamide–water and formamide–formamide dimers

The basis set superposition error (BSSE) influence in the geometry structure, interaction energies, and intermolecular harmonic and anharmonic vibrational frequencies of cyclic formamide–formamide and formamide–water dimers have been studied using different basis sets (6‐31G, 6‐31G**, 6‐31++G**, D95V, D95V**, and D95V++**). The a posteriori “counterpoise” (CP) correction scheme has been compared with the a priori “chemical Hamiltonian approach” (CHA) both at the Hartree–Fock (HF) and second‐order Møller–Plesset many‐body perturbation (MP2) levels of theory. The effect of BSSE on geometrical parameters, interaction energies, and intermolecular harmonic vibrational frequencies are discussed and compared with the existing experimental data. As expected, the BSSE‐free CP and CHA interaction energies usually show less deep minima than those obtained from the uncorrected methods at both the HF and MP2 levels. Focusing on the correlated level, the amount of BSSE in the intermolecular interaction energies is much larger than that at the HF level, and this effect is also conserved in the values of the force constants and harmonic vibrational frequencies. All these results clearly indicate the importance of the proper BSSE‐free correlation treatment with the well‐defined basis functions. At the same time, the results show a good agreement between the a priori CHA and a posteriori CP correction scheme; this agreement is remarkable in the case of large and well‐balanced basis sets. The anharmonic frequency correction values also show an important BSSE dependence, especially for hydrogen bond stretching and for low frequencies belonging to the intermolecular normal modes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Structure and stability of X4Y4H4 (XY=CC,BN)

Ab initio calculations have been carried out to study the structures and relative stabilities of the planar eight‐membered ring B₄N₄H₄ and its isoelectronic species C₈H₄ at the HF/6‐31G*, MP2/6‐31G*, MP2/6‐311G**, and MP4SDQ/6‐31G* levels. The analyses of Milliken population, vibration frequencies, π‐molecular orbital components, and orbital energy levels were used to evaluate the relative stabilities of these two similar systems. The homodesmotic reactions were also taken to be a useful index of relative stability for X₄Y₄H₄ (XY=CC, BN) and gave the resonance energies with MP4SDQ/6‐31G* of C₈H₄ (?37.2 kcal/mol) < B₄N₄H₄ (?29.2 kcal/mol). Furthermore, we calculated the thermodynamic functions of these reactions to discuss the influence of temperature. It is concluded that B₄N₄H₄ may exist in theory and could be a little more stable than C₈H₄. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 293–298, 2001     

Density functional theory study of a molecular allosteric switch for 2,2′‐bipyridyl‐3,3′‐15‐crown‐5

A classical model of “molecular machine,” which acts as an ON–OFF switch for 2,2′‐bipyridyl‐3,3′‐15‐crown‐5 ( L ), has been theoretically studied. It is highly important to understand the mechanism of this switch. The alkali‐metal cations (Na⁺ and K⁺) and W(CO)₄ fragment are introduced to coordinate with the different active sites of L , respectively. The density functional theory (DFT) method is used for understanding the stereochemical structural natures and thermodynamic properties of all the target molecules at B3LYP/6‐31G(d) and SDD (Stuttgart–Dresden) level, together with the corresponding effective core potential (ECP) for tungsten (W). The fully optimized geometries have been performed with real frequencies, which indicate the minima states. The nucleophilicity of L has been investigated by the Fukui functions. The natural bond orbital analysis is used to study the intermolecular charge‐transfer interactions and explore the origin of the internal forces of the molecular switch. In addition, the binding energies, enthalpies, Gibbs free energies, and the cation exchange energies have been studied for L , W(CO)₄ L , and their corresponding complexes. The properties of the complexes displayed by in presence or absence of the W(CO)₄ fragment are also analyzed. The calculated results of allosterism displayed by L are in a good agreement with the experimental results. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Effects of different basis sets and donor‐acceptor groups on linear and second‐order nonlinear optical properties and molecular frontier orbital energies

The calculation of molecular hyperpolarizability, molecular frontier orbital energies of some donor‐acceptor oxadiazoles ( 5a – f , 8a – f , and 9a – f ) have been investigated using ab initio methods and different basis sets. Ab initio optimizations were performed at the Hartree–Fock (HF) and density functional (Beckee‐3–Lee–Yang–Parr; B3LYP) levels of theory with 6‐31G basis set. The polarizability (<α>), anisotropy of polarizability (Δα), and ground‐state dipole moment (μ), first hyperpolarizability (β), and molecular frontier orbital (HOMO, highest occupied molecular orbital and LUMO, lowest unoccupied molecular orbital) energies of 5a – f , 8a – f , and 9a – f have been calculated at the HF and B3LYP methods with 6‐31G, 6‐31G(d), 6‐31+G(d), 6‐31++G(d,p), 6‐311G, 6‐311G(d), 6‐311+G(d), and 6‐311++G(d,p) basis sets. Also, the molecular hardness (η) and electronegativity (χ) parameters have been obtained using molecular frontier orbital energies. The <α>, Δα, μ, β, HOMO, LUMO energies, η and χ parameters have been investigated as dependence on the choice of method and basis set. The variation graphics of <α>, Δα, μ, β, η, and χ parameters using HF and B3LYP methods with different basis sets are presented. We have examined the frontier molecular orbital pictures of 5a – f , 8a – f , and 9a – f using B3LYP/6‐31++G(d,p) level. The 5a – f , 8a – f , and 9a – f display significant linear, second‐order molecular nonlinearity, and molecular parameters and provide the basis for future design of efficient nonlinear optical materials having the 1,3,4‐oxadiazole core. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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     

Are water–aromatic complexes always stabilized due to π–H interactions? LMP2 study

Eight complexes of various aromatic molecules with water have been studied theoretically at the local Møller–Plesset 2nd order theory (LMP2)/aug‐cc‐pVTZ(‐f)//LMP2/6‐31+G* level of theory. Two types of complexes can be formed, depending on the electronic structure of aromatic molecules. Donor hydrocarbons form A‐type complexes, while aromatics bearing electron‐withdrawing substituents form B‐type complexes. A‐type complexes are stabilized due to π–H interactions with the OH bond pointing to the aromatic molecule plane, while B‐type complexes have geometry with the oxygen atom pointing to the aromatic molecule plane stabilized by the interaction of highest occupied molecular orbital (HOMO) of water molecule with π* orbitals of the aromatics. It has been found that a (? HOMO–lowest unoccupied molecular orbital (LUMO)/2 value of aromatic molecule, which can be called “molecular electronegativity,” is useful to predict the type of complex formed by aromatic molecule and water. Aromatic hydrocarbons with “molecular electronegativity” of <0.15 tend to form A‐type complexes, while aromatic molecules with “molecular electronegativity” of <0.15 a.u. form B‐type complexes. The binding energy of water–aromatic complexes undergoes a minimum in the area of switching from A‐type to B type complexes, which can be rationalize in terms of frontier orbital interactions. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical study of hydrogen bonds between acetylene and selected proton donor systems

The equilibrium structures, the binding energies, and the second‐order energy components of a series of hydrogen‐bonded complexes involving acetylene are studied. The strength of the binding energy of the selected systems (HF … HCCH, HCl … HCCH, HCN … HCCH, and HCCH … HCCH) was different, ranging from a very weak interaction to a strong interaction. Calculations have been carried out at both the Hartree–Fock and correlated (second‐order Møller–Plesset perturbation theory) levels of theory, using several different basis sets [6‐31G(d,p), 6‐311G(d,p), 6‐31G++(d,p), 6‐311G++(d,p), 6‐31++G(2d,2p) and 6‐311++G(2d,2p)]. The widely used a posteriori Boys–Bernardi counterpoise (CP) correction scheme has been compared with the a priori CHA/CE, CHA–MP2, and CHA–PT2 methods, using the chemical Hamiltonian approach (CHA). The results show that at both levels the CP and the appropriate CHA results are very close to each other. Only the monomer‐based CHA‐PT2 theory gives slightly overcorrected results, reflecting that the charge transfer and polarization effects are not taken into account in this method up to second order. We have also applied our earlier developed energy decomposition scheme in order to decompose the second‐order energy contribution into different physically meaningful components. The results show that at large and intermediate intermolecular distances, the second‐order intermolecular contribution is almost equal to the sum of different physically meaningful components (e.g., polarization, charge transfer, dispersion), while at shorter distances the slightly strong overlap effects fairly disturb this simple additivity. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical studies on heats of formation for cubylnitrates using density functional theory B3LYP method and semiempirical MO methods

The heats of formation (HOF) have been calculated for all the 21 cubylnitrate compounds using the semiemprical molecular orbital (MO) methods (MINDO/3, MNDO, AM1, and PM3) and for 8 of 21 cubylnitrates containing 1–4 ? ONO₂ groups using the density functional theory (DFT) method at the B3LYP/6‐31G* level by means of designed isodesmic reactions. The cubane cage skeletons in cubylnitrate molecules have been kept in setting up isodesmic reactions to produce more accurate and reliable results. It is found that there are good linear relationships between the HOFs of the 8 cubylnitrates calculated using B3LYP/6‐31G* and two semiempirical MO (PM3 and AM1) methods, and the linear correlation coefficients of PM3 and AM1 methods are 0.9901 and 0.9826, respectively. Subsequently, the accurate HOFs at B3LYP/6‐31G* level of other 13 cubylnitrates containing 4–8 ? ONO₂ groups are obtained by systematically correcting their PM3‐calculated HOFs. Compared with noncaged nitrates, all the 21 cubylnitrates have high heats of formation implying that they may be very powerful energetic materials and have highly exploitable value. The relationship between the HOFs and the molecular structures of cubylnitrates has been discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

The effect of electron correlation on the conformational space of melatonin

The pineal gland hormone melatonin regulates several physiological processes including circadian rhythm and also alleviates oxidative stress‐induced degenerative diseases. In spite of its important biological roles, no high level ab initio conformational study has been conducted to reveal its structural features. In this work, the conformational flexibility of melatonin was investigated using correlated ab initio calculations. Conformers, obtained previously at the Hartree‐Fock level (HF/6‐31G*), were fully optimized using second order Møller‐Plesset perturbation theory applying the frozen core approximation (MP2(FC)/6‐31G*). Furthermore, single‐point MP4(SDQ,FC)/6‐31G*//MP2(FC)/6‐31G* computations were performed to investigate the effect of higher order perturbation terms. The HF and MP2 conformational spaces are considerably different: the initial 128 structures converged into 102 different local minima as confirmed by frequency calculations; 28 new minima appeared and 26 previous HF local minima disappeared; no “all‐trans” C3 side chain conformations are seen at the MP2(FC) level. The MP2 global minimum conformation is stabilized by an aromatic‐side chain interaction. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008     

Fast,efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method

The AM1‐BCC method quickly and efficiently generates high‐quality atomic charges for use in condensed‐phase simulations. The underlying features of the electron distribution including formal charge and delocalization are first captured by AM1 atomic charges for the individual molecule. Bond charge corrections (BCCs), which have been parameterized against the HF/6‐31G* electrostatic potential (ESP) of a training set of compounds containing relevant functional groups, are then added using a formalism identical to the consensus BCI (bond charge increment) approach. As a proof of the concept, we fit BCCs simultaneously to 45 compounds including O‐, N‐, and S‐containing functionalities, aromatics, and heteroaromatics, using only 41 BCC parameters. AM1‐BCC yields charge sets of comparable quality to HF/6‐31G* ESP‐derived charges in a fraction of the time while reducing instabilities in the atomic charges compared to direct ESP‐fit methods. We then apply the BCC parameters to a small “test set” consisting of aspirin, d ‐glucose, and eryodictyol; the AM1‐BCC model again provides atomic charges of quality comparable with HF/6‐31G* RESP charges, as judged by an increase of only 0.01 to 0.02 atomic units in the root‐mean‐square (RMS) error in ESP. Based on these encouraging results, we intend to parameterize the AM1‐BCC model to provide a consistent charge model for any organic or biological molecule. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 132–146, 2000     

A BSSE-free SCF algorithm for intermolecular interactions. II. Sample calculations on hydrogen-bonded complexes

The usefulness and reliability of the recent BSSE -free SCF algorithm based on the “chemical Hamiltonian approach” (CHA /F ) is demonstrated by calculating potential curves for several hydrogen-bonded complexes with 4-31G , 6-31G , and 6-31G ** basis sets. It is concluded that the CHA /F scheme gives results that are numerically close to those of the Boys–Bernardi a posteriori correction scheme but are free from the “overcompensation” characteristic of the latter at smaller distances and given basis sets. © 1992 John Wiley & Sons, Inc.     

Ab initio study of reactions of hydroxyl radicals with chloro‐ and fluoro‐substituted methanes

Ab initio calculations at the unrestricted Hartree–Fock (UHF) level have been performed to investigate the hydrogen abstraction reactions of ^? OH radicals with methane and nine halogen‐substituted methanes (F, Cl). Geometry optimization and vibrational frequency calculations have been performed on all reactants, adducts, products, and transition states at the UHF/6‐31G* level. Single‐point energy calculations at the MP2/6‐31++G* level using the UHF/6‐31G* optimized geometries have also been carried out on all species. Pre‐ and postreaction adducts have been detected on the UHF/6‐31G* potential energy surfaces of the studied reactions. Energy barriers, ΔE^?, reaction energies, ΔE_r, reaction enthalpies, ΔH_r, and activation energies, E_a, have been determined for all reactions and corrected for zero‐point energy effects. Both E_a and ΔH_r come into reasonable agreement with the experiment when correlation energy is taken into account and when more polarized and diffuse basis sets are used. The E_a values, estimated at the PMP2/6‐31++G* level, are found to be in good agreement with the experimental ones and correctly reproduce the experimentally observed trends in fluorine and chlorine substitution effects. A linear correlation between E_a and ΔH_r is obtained, suggesting the presence of an Evans–Polanyi type of relationship. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 426–440, 2001     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

Theoretical study on intermolecular interaction of epoxyethane dimer

Ab initio calculations at Hartree–Fock and fourth‐order Mø ller–Plesset (MP4) correlation correction levels with 6‐31G* basis set have been performed on the epoxyethane dimer. Dimer binding energies have been corrected for the basis set superposition error (BSSE) and the zero‐point energy. The greatest corrected dimer binding energy is −8.36 kJ/mol at the MP4/6‐31G*//HF/6‐31G* level. The natural bond orbital analysis has been performed to trace the origin of the weak interactions that stabilize dimer. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 94–98, 2000     

Potential energy surface of thionylimide

The potential energy surface (PES) of thionylimide has been searched using ab initio MO and density functional calculations. The electronic structures of the isomers of HNSO have been studied using the HF/6‐31+G*, MP2(full)/6‐31+G*, and B3LYP/6‐31+G* levels. Final energies of these molecules have been calculated at the high‐accuracy G2 and CBS‐Q levels. The probable pathways of isomerization of thionylimide to its isomers (e.g., thiocyanic acid, HONS, nitrosothiols) have been explored by studying the three‐ or four‐membered transition states. This study identified total eight possible isomers ( 1–8 ) of HNSO, of which four ( 1–4 ) have already been realized experimentally. Of the remaining four ( 5–8 ), at least two ( 5, 7 ) can be generated experimentally. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

C28 (D2): Fullerene growth mechanism

Reactions among the mono‐ to polycyclic carbon clusters have been analyzed using semiempirical AM1 and HF/6‐31G* methods. The C₂₈ (D₂) fullerene cage has been considered. Various precursors are chosen with the appropriate carbon belts. It is observed that the reactions between the precursors and the belts are essentially endoergic in nature, whereas the reactions between the stable intermediates and the final belts are exoergic. Further, the second step of the process is really akin to an annealing mechanism. Also, it has been observed that, in the annealing process, a cascade type of bond‐forming mechanism is in operation. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Alternating,gradient, block,and block–gradient copolymers of ethylene and norbornene by Pd–Diimine‐Catalyzed “living” copolymerization

We demonstrate, in this article, the facile synthesis of a broad class of low‐polydispersity ethylene–norbornene (E–NB) copolymers having various controllable comonomer composition distributions, including gradient, alternating, diblock, triblock, and block–gradient, through “living”/quasiliving E–NB copolymerization facilitated with a single Pd – diimine catalyst ( 1 ). This synthesis benefits from two remarkable features of catalyst 1 , its high capability in NB incorporation and high versatility in rendering E–NB “living” copolymerization at various NB feed concentrations ([NB]₀) while under an ethylene pressure of 1 atm and at 15 °C. At higher [NB]₀ values between 0.42 and 0.64 M, E–NB copolymerization with 1 renders nearly perfect alternating copolymers. At lower [NB]₀ values (0.11–0.22 M), gradient copolymers yield due to gradual reduction in NB concentration, with the starting chain end containing primarily alternating segments and the finishing end being hyperbranched polyethylene segments. Through two‐stage or three‐stage “living” copolymerization with sequential NB feeding, diblock or triblock copolymers containing gradient block(s) have been designed. This work thus greatly expands the family of E–NB copolymers. All the copolymers have controllable molecular weight and relatively low polydispersity (with polydispersity index below 1.20). Most notably, some of the gradient and block–gradient copolymers have been found to exhibit the characteristic broad glass transitions as a result of their possession of broad composition distribution. © 2012 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem, 2013