Gas-phase acidities of disubstituted methanes and of enols of carboxamides substituted by electron-withdrawing groups

The gas-phase acidities DeltaG degrees (acid) of some 20 amides/enols of amides RNHCOCHYY'/RNHC(OH)=CYY' [R = Ph, i-Pr; Y, Y' = CO(2)R', CO(2)R' ', or CN, CO(2)R', R', R' ' = Me, CH(2)CF(3), CH(CF(3))(2)], the N-Ph and N-Pr-i amides of Meldrum's acid, 1,3-cyclopentanedione, dimedone, and 1,3-indanedione, and some N-p-BrC(6)H(4) derivatives and of nine CH(2)YY' (Y, Y' = CN, CO(2)R', CO(2)R' '), including the cyclic carbon acids listed above, were determined by ICR. The acidities were calculated at the B3LYP/6-31+G//B3LYP/6-31+G level for both the enol and the amide species or for the carbon acid and the enol on the CO in the CH(2)YY' series. For 12 of the compounds, calculations were also conducted with the larger base sets 6-311+G and G-311+G. The DeltaG degrees (acid) values changed from 341.3 kcal/mol for CH(2)(CO(2)Me)(2) to 301.0 kcal/mol for PhNHC(OH)=C(CN)CH(CF(3))(2). The acidities increased for combinations of Y and Y' based on the order CO(2)Me < CO(2)CH(2)CF(3) < CN, CO(2)CH(CF(3))(2) for a single group and reflect the increased electron-withdrawal ability of Y,Y' coupled with the ability to achieve planarity of the crowded anion. The acidities of corresponding YY'-substituted systems follow the order N-Ph enols > N-Pr-i enols > CH(2)YY'. Better linear relationships between DeltaG degrees (acid) values calculated for the enols and the observed values than those for the values calculated for the amides suggest that the ionization site is the enolic O-H of most of the noncyclic trisubstituted methanes. The experimental DeltaG degrees (acid) value for Meldrum's acid matches the recently reported calculated value. The calculated structures and natural charges of all species are given, and the changes occurring in them on ionization are discussed. Correlations between the DeltaG degrees (acid) values and the pK(enol) values, which are linear for the trisubstituted methanes, excluding YY' = (CN)(2) and nonlinear for the CH(2)YY' systems, are discussed.     

The determination of the equilibrium structures of oxygen,ozone, and hydrogen peroxide using the ab initio and density functional theory methods

The geometries and energies of small oxygen containing molecules are studied by both the ab initio and density functional theory (DFT) methods. The RHF, MP2, and QCISD(T) ab initio methods, BHandH, BHandHLYP, BeckeSLYP, Becke3P86 DFT hybrid methods, BLYP, and the BP86 non-local DFT methods with the 3-21G¹, 6-31G(d,p), 6-311 + G(2d,2p) and 6-311 + + G(3df,3pd) basis sets were used for the computational study. The obtained results from the different methods were compared to the experimental values. The suitability of the DFT methods for reproducing experimental data were discussed.     

Quantum Chemical Model of a Lysozyme+Substrate Complex and a Cooperative Process in It

This paper presents an ab initio (RHF/6-31G** and MP2(full)6-31G**) and density functional (DFT) study of the structure and energetics of formation of an intermolecular complex which is the simplest model of an active center lysozyme with a substrate. The calculated energy of complex formation is 41.4 (RHF), 53.4 (MP2), and 52.7 kcal/mole (DFT). The proton transfer reaction is a concerted reaction having an energy barrier of 41.1 (RHF), 31.6 (MP2), and 25.3 (DFT) kcal/mole.     

Vitamin E models. Conformational analysis and stereochemistry of tetralin, chroman, thiochroman and selenochroman

Tetralin, chroman as well as its' S and Se containing congeners were subjected to ab initio (RHF/3-21G and RHF/6-31G(d)) and DFT (B3LYP/6-31G(d)) computation. Molecular geometries and the activation energies for ring inversions were determined with full geometry optimizations.     

Computational study on hydroxybenzotriazoles as reagents for ester hydrolysis

1-Hydroxybenzotriazole (1) and several of its derivatives (2-5) demonstrate potent esterolytic activity toward activated esters such as p-nitrophenyl diphenyl phosphate (PNPDPP) and p-nitrophenyl hexanoate (PNPH) in cationic micelles at pH 8.2 and 25 degrees C. The deprotonated anionic forms of such reagents act as reactive species in the hydrolysis of ester. To rationalize the origin of their nucleophilic character, a detailed ab initio/DFT computational study has been performed on 1-5 along with additional hydroxybenzotriazole derivatives (6-13). The geometries of 1-hydroxybenzotriazoles (1-13) and their corresponding bases are discussed in detail. All calculations were carried out using different methods, i.e., restricted Hartree-Fock (RHF) and hybrid ab initio/DFT (B3LYP) using 6-31G and 6-31+G basis sets. Free energy of protonation ("fep") of the 1-hydroxybenzotriazoles (1-13), free energy of solvation DeltaG(aq), and the corresponding pK(a) values have been calculated. Solvation-free energies were calculated using density functional theory and the polarizable continuum model. In addition, to examine the reliability of calculated fep, benzaldehyde oxime (14) and 2-methyl propionaldehyde oxime (15) have been computed as reference systems using different methods and basis sets, the experimental feps of which are known. Our experimental finding shows that the compound 4 is the most effective catalyst for the hydrolytic cleavages of PNPDPP and PNPH. This has been predicted from our calculated fep, pK(a), and natural charge analysis results as well. In general, the introduction of electron-withdrawing substituents on 1-hydroxybenzotriazoles facilitates the lowering of pK(a) and fep. As the pK(a) values are lowered, a greater percentage of such hydroxybenzotriazoles remain in their deprotonated, anionic forms at pH 8.2. Since the anionic forms are nucleophilic, pK(a) lowering should enhance their ester cleaving capacity. However, such substitution also decreases the charge density on the catalytically active oxido atom (O(7)). Taking these two factors together, the derivatives are only modestly better nucleophiles in comparison to the parent 1-hydroxybenzotriazole. Interestingly, the introduction of electron-donating groups does not significantly enhance the charge accumulation on the oxido atom (O(7)) of 1-hydroxybenzotriazoles.     

乙烯醇钠的从头算研究(英文)

在ＲＨＦ／３－２１Ｇ和ＲＨＦ／６－３１＋Ｇ水平上优化了乙烯醇钠的的平衡构型及过渡态，比较了不同构型的稳定性。在ＲＨＦ／６－３１＋Ｇ水平上探讨了气相中由乙醛和氢氧化钠生成乙烯醇钠的机理。并在ＭＰ２（ｆｕｌ）／６－３１＋Ｇ／／ＭＰ２（ｆｕｌ）／６－３１＋Ｇ水平下计算了反应热     

An exhaustive conformational analysis of N-acetyl-L-cysteine-N-methylamide. Identification of the complete set of interconversion pathways on the ab initio and DFT potential energy hypersurface

The full conformational space of N-acetyl-l-cysteine-N-methylamide was explored by ab initio (RHF/ 6-31G(d)) and DFT (B3LYP/6-31G(d)) computations. Multidimensional conformational analysis predicts 81 structures in N-acetyl-l-cysteine-N-methylamide, but only 47 relaxed structures were previously determined at the RHF/3-21G level of theory. These structures were now optimized using RHF/6-31G(d) and B3LYP/6-31G(d) approaches. Seven conformational migrations were observed when recalculated at higher level of theory. Besides these major changes, only smaller conformational shifts were operative for the remaining stationary points. The exploration of the whole conformational space of N-acetyl-l-cysteine-N-methylamide, including the transition-state structures allowing the conformational interconversion among the low-energy forms, was analyzed in this study. Our results offer new insights into the influence of polar side chains on the conformational preferences of peptide structures.     

Vibrational spectra and structure of benzil and its 18O- and d10-labelled derivatives: a quantum chemical and experimental study

Geometry and vibrational spectroscopic data of benzil-d0 benzil-d10 and benzil-18O calculated at various levels of theory (RHF/6-31G*, B3LYP/6-31G*, BLYP/6-31G*) are reported. The theoretical results are discussed mainly in terms of the comparisons with infrared (4000-100 cm(-1)) and Raman (4000-50 cm(-1)) spectral data. The calculated isotopic frequency shifts, induced by the 18O- and d10-labeling, are in a good agreement with the measured values. A complete vibrational assignment was made with the help of ab initio force field calculations. The data thus obtained were used for reassigning some vibrational frequencies. The results of the optimized molecular structure obtained on the basis of RHF and the DFT calculations are presented and compared with the experimental X-ray diffraction for the benzil-d0 single crystal. It turns out that the best structural parameters are predicted by the B3LYP/6-31G* method.     

Cyano-, nitro-, and alkoxycarbonyl-activated observable stable enols of carboxylic acid amides

A search for the enol structures of several amides YY'CHCONHPh with Y,Y' = electron-withdrawing groups (EWGs) was conducted. When Y = CN, Y' = CO(2)Me the solid structure is that of the enol (8b) MeO(2)CC(CN)=C(OH)NHPh, whereas in solution the NMR spectrum indicate the presence of both the amide MeO(2)CCH(CN)CONHPh (8a) and 8b. When Y = NO(2), Y' = CO(2)Et the main compound in CDCl(3) is the amide, but <10% of enol(s), presumably EtO(2)CC(NO(2))=C(OH)NHPh (9b), are also present. When Y = COEt, Y' = CO(2)Me or Y = COMe, Y' = CO(2)Et (10 and 11) enolization in solution and of 11 also in the solid state occurs at the carbonyl rather than at the ester site. With Y = Y' = CN a rapid exchange between the amide (NC)(2)CHCONHPh (12a) and a tautomer, presumably the enol, take place in several solvents on the NMR time scale. With YY' = barbituric acid moiety the species in DMSO-d(6) is an enol of an amide although which CONH group enolizes is unknown. B3LYP/6-31G calculations showed that the enol (NC)(2)C=C(OH)NH(2) (13b) is more stable by DeltaG of 0.4 kcal/mol than (NC)(2)CHCONH(2) (13a) due to a combination of stabilization of 13b and destabilization of 13a and both are much more stable than the hydroxyimine and ketene imine tautomers. The effect of Y,Y' and the solvent on the relative stabilization of enols of amides is discussed.     

Theoretical Study on Electronic Gain-and-loss Properties of TEMPO and Its Derivates in Charge/Discharge Processes

Theoretical study on the electronic structures and related properties of 2,2,6,6-tetramethyl-1-piperidinyloxy (TEMPO) and its cationic lipid derivates in the charge/discharge processes has been carried out using the density functional theory (DFT) at the (U)B3LYP/6-31G(d,p) or 6-31+G(d,p) level. The changes and regularities of geometric and electronic properties of these compounds in the charge/discharge processes were revealed in detail. The compu-tational results show that the substitute group plays a very important role in the electronic structures and related properties of TEMPOs during the charge/discharge processes. It isvery interesting to find that after getting an electron, TEMPO is more stable in singlet state but the lipid is more stable in triplet state. For TEMPO, both the charge and the discharge processes greatly influence the electronic properties of N and O atoms of the radical part. For the cationic lipid, the discharge process mainly influences the pyridinium head and the charge process mainly influences the free radical head. Moreover, the solvent effect plays an important role in some bond lengths and the charge population of the free radical head. In addition, the UV-Vis absorption spectra of TEMPO and the lipid were calculated and simulated using TDDFT at the 6-31G(d,p) or 6-31+G(d,p) level, in satisfying agreement with the experimental ones.     

Conformational stability, optimized geometries, vibrational and electronic spectra of methacryloyl bromide in ground and excited electronic states

In order to understand conformational isomerism in methacryloyl bromide (MABR) in the ground (S(0)) and the first excited (S(1)) electronic states and to interpret the vibrational and electronic spectra of its conformers in the S(0) state, quantum mechanical calculations using Density Functional Theory (DFT) and RHF methods with extended basis sets 6-31G, 6-31G** and 6-311+G(d,p) have been conducted. In RHF calculations, electron correlation effects have been included at the M?ller-Plesset MP2 level. It is inferred that in both the electronic states the molecule may exist in two isomeric forms-s-trans and s-cis; the former being more stable than the later by about 1.629 kcal mol(-1) in the S(0) state and by about 2.218 kcal mol(-1) in the S(1) state. Electronic transition tends to increase the s-trans/s-cis and s-cis/s-trans, rotational barriers from 7.059 kcal mol(-1) (2468.1 cm(-1)) and 5.428 kcal mol(-1) (1897.8 cm(-1)) in S(0) state to 23.594 kcal mol(-1) (8249.4 cm(-1)) and 21.376 kcal mol(-1) (7473.9 cm(-1)) in the S(1) state. Completely optimized geometries of the two conformers in S(0) state reveal that while there is no significant difference in their bond lengths, some of the bond angles associated with COBr group are appreciably different. Electronic excitation tends to change both the bond lengths and bond angles. Based on suitably scaled DFT and RHF results obtained from the use of 6-31G** and 6-311+G(d,p) basis sets, a complete assignment is provided to the fundamental vibrational bands of both the s-trans and s-cis conformers in terms of frequency, form and intensity of vibrations and potential distribution across the symmetry coordinates in the S(0) state and a comparison has been made with experimental assignments. A theoretical prediction of the electronic transitions in the near UV-region in the two conformers and their tentative assignment has been provided on the basis of CI level calculations using 6-31G basis set.     

Density functional theory predicts the barriers for radical fragmentation in solution

[reaction: see text] N-Methoxypyridyl radicals formed by one-electron reduction of the corresponding cationic heterocycles undergo N-O bond cleavage. Experimental activation free energies for a series of these bond fragmentations are compared to corresponding barriers determined from electronic structure calculations. The DFT barriers agree well with those from experiment, being smaller than the latter values by an average value of ca. 1 kcal/mol, for rate constants varying over almost 3 orders of magnitude, or within ca. 3 kcal/mol over 8 orders of magnitude of rate constant. For a model compound, the B3PW91/6-31+G hybrid density functional method is also found to be in good agreement with the MCSCF-MRMP2 method. One of the reactions is found by DFT to have no minimum for the reactant radical, consistent with a truly barrierless reaction.     

Double Bond Migration Mechanism in Allyl Systems Involving the Hydroxide Ion. 2. Polarizable Continuum Model (PCM)

Isomerization processes of a double bond site in propene and methylthiopropene molecules with the hydroxide ion were studied in the framework of the RHF/6-31+G*, MP2/6-31+G*, and B3LYP/6-31+G* (density functional) ab initio methods. The solvent effect was taken into account using PCM in its IEFPCM and SCIPCM versions. It is shown that to construct a reaction profile for propene rearrangement, it suffices to perform geometry optimization of stationary points within the Born–Onsager model with further refinement of the energy using IEFPCM. The reaction profiles obtained display that the multiple bond migration mechanism involving the hydroxide ion proton is energetically preferable to the two-stage mechanism forming a solvated carbanion for the propene molecule and for the methylthiopropene molecule that forms a much more stable carbanion.     

Density functional theory study of 14N isotropic hyperfine coupling constants of organic radicals

Nitrogen hyperfine coupling constants (hfccs) of organic radicals have been calculated by density functional theory (DFT) methodology. The capability of the B3LYP functional, combined with 6-31G*, TZVP and EPR-III basis sets, to reproduce experimental nitrogen coupling constant data has been analyzed for 109 neutral, cationic and anionic radicals, all of them containing at least one nitrogen atom. The results indicate that the selection of the basis set plays an important role in the accuracy of DFT calculations of hfccs, mainly in relation with the composition of the primitive functions and the quantum number of those functions. The main conclusion obtained is the high reliability of the scheme B3LYP/6-31G* for the prediction of nitrogen hfccs with very low computational cost.     

Effect of heteroatom insertion at the side chain of 5-alkyl-1H-tetrazoles on their properties as catalysts for ester hydrolysis at neutral pH

[reaction: see text] Herein we introduce tetrazole and its suitably designed derivatives as powerful ester-cleaving reagents. By first performing a detailed ab initio computational study, we found that, in the side chain of 5-alkyl-1H-tetrazoles, introduction of a heteroatom (e.g., N, O, or S at the alpha-position of the tetrazole ring) raises the charge on the tetrazole nucleus significantly. All calculations have been performed using restricted Hartree-Fock (RHF) and hybrid ab initio/DFT (B3LYP) methods employing 6-31G* and 6-31+G* basis sets. To estimate the nucleophilicity of these reagents, the charges on conjugate bases of various tetrazole derivatives have been calculated using natural population (NBO) analysis in gas phase and in water. Free energy of protonation (fep) of the 1H-tetrazole derivatives (1-7), free energy of solvation, deltaG(aq), and the corresponding pKa values have been calculated by self-consistent reaction field (SCRF) methods applying the polarized continuum model (PCM). Since the calculation indicates that incorporation of heteroatom leads to enhanced nucleophilicity in their deprotonated anionic tetrazole forms, a series of 5-substituted 1H-tetrazole derivatives have been synthesized. These compounds indeed catalyze the hydrolysis of p-nitrophenyl diphenyl phosphate (PNPDPP) and p-nitrophenyl hexanoate (PNPH) efficiently in cationic cetyl trimethylammonium bromide (CTABr) micelles at pH 7.0 and 25 degrees C. The pseudo-first-order rate constants (k(obs)) were determined for each catalyst against both substrates. The experimental and theoretical results show that, to achieve better k(obs) values for the cleavage of PNPDPP and PNPH under micellar conditions, charge on the N- atom (nucleophile) of conjugate base is important. Replacing the alpha-CH2 in alkyl substituent with S (3), NH (4), or O (5) enhances the accumulation of charge on N- in conjugate bases of tetrazoles and subsequently increases their intrinsic nucleophilic reactivity toward hydrolytic reactions. Significantly large rate enhancements were observed for the cleavage of PNPDPP and PNPH at pH 7.0 in the presence of catalytic system 5/CTABr over background (only CTABr). Tetrazole 4 (alpha-isomer) showed 4-5-fold superior reactivity over 6 (beta-isomer) under identical conditions. Natural charges obtained from NBO analysis (B3LYP/6-31+G*) are -0.94 and -0.852 on N- in the conjugate bases of 4 and 6, respectively. This also predicts that 4 is a better nucleophile than 6. All the newly synthesized tetrazole derivatives in micellar media display true catalytic properties by cleaving several fold excess of substrates.     

Double Bond Migration Mechanism in Allyl Systems Involving the Hydroxide Ion. 3. Supermolecular Approach

The mechanism of a proton transfer in an allyl system involving the hydroxide ion was investigated by the RHF/6-31+G* and MP2/6-31+G* methods with one, two, and four water molecules (water regarded as a solvent) included in the calculation; in the fourth model, the reactant system and two H₂O molecules were placed into continuum. A possible mechanism of intramolecular proton transfer involving easy exchange of water molecules between the first and second coordination spheres of the propenide ion is considered. The results are compared with those obtained with gas-phase and polarizable continuum models.     

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.     

Conformational analysis of 2,4-dialkyl-1,3,2-dioxaborinanes

Computer modeling using the quantum-chemical empirical MM+ method and nonempirical RHF//6-31G(d), MP2//6-31G(d), and, in individual cases, RHF//3-21G methods was employed to study the potential energy surface of 2,4-dialkyl-1,3,2-dioxaborinanes. The optimal geometry and ¹H NMR spectral data gave values for ∆G ⁰ of the methyl and hexyl substituents at C-4 of the heterocyclic ring equal to 0.6 and 1.6 kcal/mol, respectively.     

Geminal interaction and effect of “positive charge”

Quantum-chemical calculations of the molecules containing Cl-C-M group (M=C, N, O, F) were carried out using RHF/6-31G(d) and MP2/6-31+G(d) methods. The calculations are aimed at verifying concept on “positive charge” effect of the central atom of the group suggested as understanding of non-inductive effect of heteroatom M on the Cl (Y) atom in Y-Z-M moieties (geminal interaction). It is shown that C atom in such systems bears negative or small positive charge, and hence “positive charge” effect is excluded. The non-inductive effect originates from the polarization of Cl-C (Y-Z) bonds due to the influence of the charge on any M atom directly through space on any Y and Z.