Ab initio study of anomeric effect in 2,2-difluoroglycine

The optimized molecular structures of seven conformations of 2,2-difluoroglycine have been obtained from ab initio calculations. For conformers in which the lone pair of electrons on the nitrogen are antiperiplanar to one of the C–F bonds, that C–F bond is longer than the other C–F bond, which is synperiplanar to the lone pair of electrons. Conformers which have these features are the most stable conformers of those examined. This observation is explained in terms of an anomeric effect of the 1p(N)→σ^*(C–F). At the MP2/6-31G^* level of calculation, conformers IV and V are 21.5 and 18.7 kJ/mol, respectively, more stable than the least stable conformer, VI, which does not exhibit an anomeric effect. Conformer VII was found to be exceptionally stable, in addition to an anomeric effect, this conformer also exhibits features of a FH–O hydrogen bond.     

Ab initio studies of the conformers and conformational distribution of the gaseous hydroxyamino acid threonine

Systematic and extensive conformational search has been performed to characterize the gas-phase threonine structures. A total of 1296 unique trial structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were optimized at the B3LYP/6-311G* level of the theory and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 71 conformers were found and their rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies and vertical ionization energies of all the conformers were determined. Single-point energies were also calculated at the MP2/6-311G(2df,p) and B3LYP/6-311G(2df,p) levels. Characteristic H-bonding types were classified and listed for all the conformers. The conformational distributions of gaseous threonine at various temperatures were calculated.     

Ab initio studies of the intramolecular amide hydrolysis in N-methylmaleamic acids

The intramolecular amide hydrolysis reactions of N-methylmaleamic acids (NMMA) are studied at the MP2/6-31G**//RHF/4-31G level of theory as model reactions of peptide bond cleavage by a proteolytic enzyme. In contrast to the previously reported results for a bimolecular reaction model of peptide hydrolysis, the unimolecular reactions studied here proceed via the concerted pathway in which the C–O bond formation and the release of methylamine occur simultaneously in preference to the stepwise one. The determination of an intrinsic reaction coordinate shows that the reaction is facilitated by the intramolecular proton transfer from the undissociated carboxyl group to the nitrogen of the leaving amine group. Mainly because of the increase in activation energy, methyl substitution at the 2-position retards the hydrolysis reaction rate by a factor of 14 compared to the reaction of the unsubstituted molecule. In contrast, additional methyl substitution at 3-position leads to 35-fold increase in the reaction rate. These variations of reactivity are caused by the charge redistributions in the amide group induced by methyl substitutions and the resulting changes in electrophilicity of the aminocarbonyl carbon.     

CBS-4 ab initio results on the formation and stability of N-aminoazetidine radical cation and derivatives

Experimental evidence has recently indicated the role of substituted azetidine radical cations in the electron ionization mass spectrometric fragmentation of N,N-dimethyltydrazone derivatives of 5- and 6-unsaturated straight-chain aldehydes. Although the 2+2 cycloaddition of ethylene and methyleneamine to form neutral azetidine has a very high energy barrier and is only mildly exothermic, CBS-4 ab initio calculations show that formation of the proposed radical cation intermediate is strongly exothermic with a low barrier. Further calculations show that the general energetics remain essentially unchanged even when substituents of the type present in the mass spectrometric studies are included. These calculations demonstrate that the proposed new mass spectral fragmentation mechanism is both thermodynamically and kinetically feasible.     

Ab initio and DFT studies of the vibrational spectra of benzofuran and some of its derivatives

The vibrational spectra of benzofuran and some of its derivatives have been systematically investigated by ab initio and density functional B3LYP methods. The harmonic vibrational wavenumbers and intensity of vibrational bands were calculated at ab initio and DFT levels invoking different basis sets up to 6-311++g**. Vibrational assignments have been made and it has been found that the calculated DFT frequencies agree well in most cases with the observed frequencies for each molecule. Conformational studies have also been carried out and it is evident from ab initio calculations that 2(3H) benzofuranone is more stable than 3(2H) benzofuranone in support to our earlier semiempirical results.     

CBS-4 ab initio results on the formation and stability of N-aminoazetidine radical cation and derivatives

Experimental evidence has recently indicated the role of substituted azetidine radical cations in the electron ionization (EI) mass spectrometric fragmentation of N,N-dimethylhydrazone derivatives of 5- and 6-unsaturated straight-chain aldehydes. Although the 2+2 cycloaddition of ethylene and methyleneamine to form neutral azetidine has a very high energy barrier and is only mildly exothermic, CBS-4 ab initio calculations show that formation of the proposed radical-cation intermediate is strongly exothermic with a low barrier. Further calculations show that the general energetics remain essentially unchanged even when substituents of the type present in the mass spectrometric studies are included. These calculations demonstrate that the proposed new mass spectral fragmentation mechanism is both thermodynamically and kinetically feasible.     

Photophysics of inter- and intra-molecularly hydrogen-bonded systems: Computational studies on the pyrrole–pyridine complex and 2(2′-pyridyl)pyrrole

The role of electron and proton transfer processes in the photophysics of hydrogen-bonded molecular systems has been investigated with ab initio electronic-structure calculations. We discuss generic mechanisms of the photophysics of a hydrogen-bonded aromatic pair (pyrrole–pyridine), as well as an intra-molecularly hydrogen-bonded π system composed of the same molecular sub-units (2(2′-pyridyl)pyrrole). The reaction mechanisms are discussed in terms of excited-state minimum-energy paths, conical intersections and the properties of frontier orbitals. A common feature of the photochemistry of these systems is the electron-driven proton transfer (EDPT) mechanism. In the hydrogen-bonded complex, a highly polar charge transfer state of ¹ππ^* character drives the proton transfer, which leads to a conical intersection of the S₁ and S₀ surfaces and thus ultrafast internal conversion. In 2(2′-pyridyl)pyrrole, out-of-plane torsion is additionally needed for barrierless access to the S₁–S₀ conical intersection. It is pointed out that the EDPT process plays an essential role in the fluorescence quenching in hydrogen-bonded aromatic complexes, the function of organic photostabilizers, and the photostability of biological molecules.     

Ab initio and DFT studies on the mechanism of the 1,3-dipolar cycloaddition reaction between nitrone and sulfonylethene chloride

The molecular mechanism for the 1,3-dipolar cycloaddition of nitrone with sulfonylethene chlorides has been studied using ab initio and DFT methods at the HF, MP2 and B3LYP levels together with the 6-31G* basis set. Relative rates, stereo and regioselectivity, have been analysed and discussed. For this cycloaddition four reactive channels associated with the formation of two pairs of diastereoisomeric regioisomers have been characterized. Analysis of the geometries of the corresponding transition structures shows that the cycloaddition takes place along a concerted but asynchronous mechanism. Activation energies as asynchronicity are dependent on the computation level. Thus, while HF calculations gave large barriers, MP2 calculations tend to underestimate them. DFT calculations gave reasonable values. These 1,3-dipolar cycloadditions present an endo stereoselectivity while the meta regioselectivity depends on the computational level. Thus, while HF and DFT calculations predict meta path, in agreement with the experimental results, MP2 calculation predict ortho regioselectivity. The frontier molecular orbitals analysis shows that the reaction is controlled by the (HOMO_dipole–LUMO_{dipolarophile}) interaction in agreement with the charge transfer analysis carried out at the transition structures. Inclusion of diffuse functions at the B3LYP/6-31+G* level increases the energy barriers about 4 kcal/mol, giving a similar endo/meta selectivity. Solvent effects have been taken into account, by means of self-consistent reaction field.     

Ab initio and density function theory computational studies of the CH4 + H → CH3 + H2 reaction

The activation barrier for the CH₄ + H → CH₃ + H₂ reaction was evaluated with traditional ab initio and Density Functional Theory (DFT) methods. None of the applied ab initio and DFT methods was able to reproduce the experimental activation barrier of 11.0-12.0 kcal/mol. All ab initio methods (HF, MP2, MP3, MP4, QCISD, QCISD(T), G1, G2, and G2MP2) overestimated the activation energy. The best results were obtained with the G2 and G2MP2 ab initio computational approaches. The zero-point corrected energy was 14.4 kcal mol⁻¹. Some of the exchange DFT methods (HFB) computed energies which were similar to the highly accurate ab initio methods, while the B3LYP hybrid DFT methods underestimated the activation barrier by 3 kcal mol⁻¹. Gradient-corrected DFT methods underestimated the barrier even more. The gradient-corrected DFT method that incorporated the PW91 correlational functional even generated a negative reaction barrier. The suitability of some computational methods for accurately predicting the potential energy surface for this hydrogen radical abstraction reaction was discussed.     

Ab initio studies on the mechanism of dimerization reactions of ketene imine and bis(trifluoromethyl)ketene imine

The dimerization reactions of ketene imine and bis(trifluoromethyl)ketene imine were studied theoretically. All the dimerization processes take place in a concerted but asynchronous manner, each proceeding through a four-membered ring transition state. For the ketene imine dimerization reactions, three different processes have almost equal activation barriers, while for the three bis(trifluoromethyl)ketene imine dimerization processes the reaction giving symmetrical a four-membered heterocyclic product has the lowest activation barrier. Received: 15 July 1998 / Accepted 3 September 1998 / Published online: 17 December 1998     

Ab initio study on the thermolysis of β-hydroxyolefins in gas phase

Thermolysis studies of β-hydroxyolefins in gas phase were realized using ab initio MP2 and DFT methods at the 6-31G^* levels to explore the possibility of determining a possible concerted process with a six-membered cyclic transition state (TS). Vibrational frequency calculations were carried out in order to confirm the stationary states, including TS structures. IRC calculations have been performed in all cases in order to verify that localized TS structures connect with the corresponding minimum stationary points associated with the reactant and products. With the aim of corroborating the postulated mechanism in the experimental study, we present a theoretical study in order to calculate the rate constants and the activation parameters. The results obtained are in accordance with the experimental conclusions.     

Investigation of the effects of ionic hydrogen bonds (COH–C and N–HOC) in crystalline dl-proline by ab initio and DFT calculated NQR parameters

In this paper, we have calculated the nuclear quadrupole resonance (NQR) parameters of the quadrupole nuclei involved in the hydrogen bonds (CO⁻H–C and ⁺N–H⁻OC) in the monomer and pentameric cluster of dl-proline by HF and B3LYP methods and basis sets of 6-311+G^* and 6-311++G^**. These computations are performed on the basis of X-ray diffraction structural data of dl-proline. The results indicate that the calculations including hydrogen-bonding (HB) interactions (in pentamer) are in better agreement with the experimental data than those in which these interactions are neglected (in monomer). The quantum chemical calculations show that the intermolecular hydrogen-bonding interactions play an important role in determination of the NQR parameters of ¹⁴N, ²H of group and ¹⁷O.     

Ab initio and density functional studies on the structure and vibrational spectra of 2-hydroxy-1,4-naphthoquinone-1-oxime derivatives

Structure and vibrational frequencies of lawsoneoxime and its C₃-substituted (R=CH₃, NH₂, Cl, NO₂) derivatives in keto and nitrosophenol forms have been obtained employing the Hartree–Fock and density functional methods. Charge distributions in different conformers have been studied using the molecular electrostatic potential topography as a tool. For all these derivatives except for nitrolawsoneoxime the amphi conformer in the keto form is predicted to be of lowest energy, which can partly be attributed to hydrogen bonding through the oximino nitrogen. In the nitro derivative, however, the preference to form a six membered ring owing to O–H–O hydrogen-bonded interactions makes the anti conformer (keto) the stablest. Further one of the nitrosophenol conformers of nitrolawsoneoxime turns out to be very close in energy (0.21 kJ mol^–1 higher) to this anti conformer. The consequences of hydrogen bonding on charge distribution and vibrational spectra are discussed.     

Ab initio studies on the pyrolysis mechanism of 2-haloacetic acids in the gas phase

The dehydrohalogenation mechanism of 2-haloacetic acids (XCH2CO2H, X=F, Cl and Br) has been studied theoretically by HF/3-21G and AM1 methods. The results indicate that these reactions are most probably proceeded in terms of a polar five-membered cyclic transition state in the gas phase. Their microscopic processes are beleived to be a stepwise reaction and the rate-determining step is the first one. By comparing the energy barriers of different 2-haloacetic acids, it can be realized that 2-fluoroacetic acid is easier to react than 2-chloroacetic and 2-bromoacetic acids.     

Ab initio calculation of the electronic coupling element in bimetallic model compounds [M---L---M], M = Be, Mg, Zn; L = O, S, ---CH2---, ---CC---; electron correlation effects and dependence on the bridge nature

Ab initio calculations of the electron-transfer matrix element V_ab in bimetallic model systems [M---M]⁺ and [M---L---M]⁺ with M = Be, Zn and L = O, S, ---CH₂--- and ---CC---, are reported. Electron correlation effects are estimated through a second-order method in which only the differential space contributing to the transfer integral H_ab is considered. V_ab is found to be noticeably larger for bridged systems, showing that the through bond mechanism is pre-eminent in these systems. In most of the calculations, electron correlation contributions are found to increase V_ab, the corrections in some cases being as large as 35%. Analysis of the second-order contributions in these cases shows that intruder states involving ligand or metal ligand excitations strongly mix with the zeroth-order model wave functions.     

Kinetic parameters for the reaction of hydroxyl radical with CH3OCH2F (HFE‐161) in the temperature range of 200–400 K: Transition state theory and Ab initio calculations

Rate coefficients for the reaction of the hydroxyl radical with CH₃OCH₂F (HFE‐161) were computed using transition state theory coupled with ab initio methods, viz., MP2, G3MP2, and G3B3 theories in the temperature range of 200–400 K. Structures of the reactants and transition states (TSs) were optimized at MP2(FULL) and B3LYP level of theories with 6‐31G* and 6‐311++G** basis sets. The potential energy surface was scanned at both the level of theories. Five different TSs were identified for each rotamer. Calculations of Intrinsic reaction coordinates were performed to confirm the existence of all the TSs. The kinetic parameters due to all different TSs are reported in this article. The rate coefficients for the title reaction were computed to be k = (9 ± 1.08) × 10^?13 exp [?(1,713 ± 33)/T] cm³ molecule^?1 s^?1 at MP2, k = (7.36 ± 0.42) × 10^?13 exp [?(198 ± 16)/T] cm³ molecule^?1 s^?1 at G3MP2 and k = (5.36 ± 1.57) × 10^?13 exp [?(412 ± 81)/T] cm³ molecule^?1 s^?1 at G3B3 theories. The atmospheric lifetimes of CH₃OCH₂F at MP2, G3MP2, and G3B3 level of theories were estimated to be 20, 0.1, and 0.3 years, respectively. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Ab initio quantum chemical study on the mechanism of exceptional behavior of lysine for ion yields in MALDI--role of vibrational entropic contribution in thermally averaged proton affinities

To elucidate the mechanism of the exceptional behavior of lysine for the ionization (protonation) yields in matrix-assisted laser desorption/ionization (MALDI) observed by Nishikaze and Takayama [Rapid Commun. Mass Spectrom. 2006, 20, 376], the temperature dependences of proton affinity (PA) and gas phase basicity for 20 amino acids are theoretically analyzed with correlated ab initio molecular orbital method under ideal gas condition. We have found that two different conformations, the linear structure with elongation of the side chain and the folded one having intramolecular hydrogen bonding, play important roles for the exceptional behavior of lysine. At low temperatures of around 298 K, the most stable conformation of the protonated lysine is the folded structure due to the formation of intramolecular hydrogen bonding. Meanwhile, at high temperatures, the Gibbs free energy of linear structure of protonated lysine becomes lower than that of the folded one because of the increment of vibrational entropic contribution. To explicitly take account of the contribution of the free energies, we have proposed the effective PA values thermally averaged using the ratio of Boltzmann distributions for two conformations. Since the effective PA value for lysine drastically decreases as the temperature increases above 1000 K, the linear correlation is clearly obtained between our effective PA values at high temperature and the ion yields in MALDI.     

Ab initio and density functional theory studies on the mechanism of nucleophilic vinylic substitution of 4<Emphasis Type="Italic">H</Emphasis>-pyran-4-one and 2-methyl-4<Emphasis Type="Italic">H</Emphasis>-pyran-4-one with ammonia

 Nucleophilic vinylic substitutions of 4H-pyran-4-one and 2-methyl-4H-pyran-4-one with ammonia were calculated by the B3LYP method using the 6-31G(d,p) basis set. Bulk solvent effects of aqueous solution were estimated by the polarized continuum and Poisson–Boltzmann self-consistent reaction field models using the 6-311+G(d,p) basis set. In the gas phase different mechanisms were found for the two reaction systems calculated. The reaction of 4H-pyran-4-one proceeds through enol, whereas a feasible path for the less reactive 2-methyl-4H-pyran-4-one is the mechanism through a keto intermediate. Addition of ammonia in concert with proton transfer is the rate-determining step ofthe reaction. The mechanism proceeding either by a bimolecular nucleophilic substitution (S_N2) or by one involving a tetrahedral zwitterionic intermediate is shown to be unlikely in the gas phase or nonpolar solution. The effects of bulk solvent not only consist in a reduction of the various activation barriers by about 25–40 kJ mol⁻¹ but also in a change in the reaction mechanism. Received 26 May 2002 / Accepted 26 July 2002 / Published online: 14 February 2003