NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

Theoretical study of geometrical effect on the deoxygenation of epoxide by singlet carbenes

The mechanism of the abstraction of oxygen from epoxide by carbenes has been investigated using the B3LYP/6-31G* and MP2/6-311++G(d,p)//B3LYP/6-31G* levels of theory. According to our model calculations, the reactivity of carbene decreases in the order: four-membered ring carbene > six-membered ring carbene > three-membered ring carbene > five-membered ring carbene vinylidene. Moreover, the present work shows that the singlet–triplet splitting of a carbene can be used as a guide for the prediction of its reactivity.     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Hydrogen peroxide clusters: the role of open book motif in cage and helical structures

Hartree-Fock (HF) calculations using 6-31G*, 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets show that hydrogen peroxide molecular clusters tend to form hydrogen-bonded cyclic and cage structures along the lines expected of a molecule which can act as a proton donor as well as an acceptor. These results are reiterated by density functional theoretic (DFT) calculations with B3LYP parametrization and also by second-order M?ller-Plesset perturbation (MP2) theory using 6-31G* and 6-311++G(d,p) basis sets. Trends in stabilization energies and geometrical parameters obtained at the HF level using 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets are similar to those obtained from HF/6-31G* calculation. In addition, the HF calculations suggest the formation of stable helical structures for larger clusters, provided the neighbors form an open book structure.     

Structural elucidation of nitro-substituted five-membered aromatic heterocycles utilizing GIAO DFT calculations

The GIAO (Gauge Including Atomic Orbitals) DFT (Density Functional Theory) method is applied at the B3LYP/6-31+G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311+G (2d,p)//B3LYP/6-31+G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants for 25 nitro-substituted five-membered heterocycles. Difference (1D NOE) spectra in combination with long-range gHMBC experiments were used as tools for the structural elucidation of nitro-substituted five-membered heterocycles. The assigned NMR data (chemical shifts and coupling constants) for all compounds were found to be in good agreement with theoretical calculations using the GIAO DFT method. The magnitudes of one-bond (1JCH) and long-range (nJCH, n>1) coupling constants were utilized for unambiguous differentiation between regioisomers of nitro-substituted five-membered heterocycles.     

DFT study of the gas‐phase thermal decomposition kinetics of 2‐ethoxypyridine into 2‐pyridone

The mechanism of the gas‐phase elimination kinetics of 2‐ethoxypyridine has been studied through the electronic structure calculations using density functional methods: B3LYP/6‐31G(d,p), B3LYP/6‐31++G(d,p), B3PW91/6‐31G(d,p), B3PW91/6‐31++G(d,p), MPW1PW91/6‐31G(d,p), MPW1PW91/6‐31++G(d,p), PBEPBE/6‐31G(d,p), PBEPBE/6‐31++G(d,p), PBE1PBE1/6‐31G(d,p), and PBE1PBE1/6‐31++G(d,p). The elimination reaction of 2‐ethoxypyridine occurs through a six‐centered transition state geometry involving the pyridine nitrogen, the substituted carbon of the aromatic ring, the ethoxy oxygen, two carbons of the ethoxy group, and a hydrogen atom, which migrates from the ethoxy group to the nitrogen to give 2‐pyridone and ethylene. The reaction mechanism appears to occur with the participation of π‐electrons, similar to alkyl vinyl ether elimination reaction, with simultaneous ethylene formation and hydrogen migration to the pyridine nitrogen producing 2‐pyridone. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

On the Importance of CP-corrected Gradient Optimization in the Study of Hydrogen Bonded Systems

Geometries, harmonic vibrational frequencies and interaction energies of the water-hydrogen sulfide dimer, hydrogen fluoride dimer and glycine zwitterion-water dimer were determined by the counterpoise-corrected (CP-corrected) gradient optimization that explicitly corrects for the basis set superpusition error (BSSE) and CP-uncorrected (normal) gradient opfimization respectively at the B3LYP and MP2 levels of theory, employing the popular Pople‘s standard 6-31G(d), 6-31G(d,p) and 6-311 G(d,p) basis sets in order to assess the importance of CP-corrected gradient optimiTation in the study of hydrogen bonded systems. The normal optimization of these three H-bonded systems obtained using these popular basis sets all yielded erratic results, whereas use of CP-corrected gradient optimization led to consistent results with those from larger basis sets. So this CP receipt becomes useful and necessary to correctly describe large systems, where the use of small basis sets may be necessary.     

Kinetics of the multichannel reaction of methanethiyl radical (CH3S*) with 3O2

The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.     

A theoretical study of nonlinearity in heterocyclic hydrogen-bonded complexes

Ab initio calculations are performed at the MP2/6-311++G(d,p) and DFT/B3LYP/6-311++G(d,p) theoretical levels to obtain geometries, H-bond energies and harmonic infrared vibrational properties for the Cs symmetry structures of heterocyclic hydrogen-bonded complexes, CnHmY-HX. The H-bond lengths in DFT/B3LYP calculation level are in better agreement with the experimental values than the MP2 results. The geometry optimization are interpreted in terms of hydrogen bond nonlinearity represented by theta; and phi angles, once the hydrogen bond is formed among n-electrons pairs of the heteroatom in heterocyclic and the hydrogen atom in HX. The hydrogen bond energy after of the zero-point vibrational energy (ZPE) and basis set superposition error (BSSE) corrections are overestimated at DFT/B3LYP, whereas the MP2 BSSE corrections are very large than corresponding DFT/B3LYP. For example, the BSSE corrections for the C2H4S-HNC complex are 7.60 and 0.09 kJ mol(-1) in MP2 and DFT/B3LYP calculations levels, respectively. The new vibrational modes in infrared harmonic spectrum arising from complexation show several interesting features, especially the intermolecular stretching mode.     

TDDFT investigation on nucleic acid bases: comparison with experiments and standard approach

A comprehensive theoretical study of electronic transitions of canonical nucleic acid bases, namely guanine, adenine, cytosine, uracil, and thymine, was performed. Ground state geometries were optimized at the MP2/6-311G(d,p) level. The nature of respective potential energy surfaces was determined using the harmonic vibrational frequency analysis. The MP2 optimized geometries were used to compute electronic vertical singlet transition energies at the time-dependent density functional theory (TDDFT) level using the B3LYP functional. The 6-311++G(d,p), 6-311(2+,2+)G(d,p), 6-311(3+,3+)G(df,pd), and 6-311(5+,5+)G(df,pd) basis sets were used for the transition energy calculations. Computed transition energies were found in good agreement with the corresponding experimental data. However, in higher transitions, the Rydberg contaminations were also obtained. The existence of pisigma* type Rydberg transition was found near the lowest singlet pipi* state of all bases, which may be responsible for the ultrafast deactivation process in nucleic acid bases.     

The shape of gaseous n-butylbenzene: assessment of computational methods and comparison with experiments

Conformational landscape of neutral and ionized n-butylbenzene has been examined. Geometries have been optimized at the B3LYP/6-31G(d), B3LYP/6-31+G(d,p), B3LYP-D/6-31+G(d,p), B2PLYP/6-31+G(d,p), B2PLYP-D/6-31+G(d,p), B97-D/6-31+G(d,p), and M06-2X/6-31+G(d,p) levels. This study is complemented by energy computations using 6-311++G(3df,2p) basis set and CBS-QB3 and G3MP2B3 composite methods to obtain accurate relative enthalpies. Five distinguishable conformers have been identified for both the neutral and ionized systems. Comparison with experimentally determined rotational constants shows that the best geometrical parameters are provided by B3LYP-D and M06-2X functionals, which include an explicit treatment of dispersion effects. Composite G3MP2B3 and CBS-QB3 methods, and B2PLYP-D, B3LYP-D, B97-D, and M06-2X functionals, provide comparable relative energies for the two sets of neutral and ionized conformers of butyl benzene. An exception is noted however for conformer V(+) the stability of which being overestimated by the B3LYP-D and B97-D functionals. The better stability of neutral conformers I, III, and IV, and of cation I(+) , demonstrated by our computations, is in perfect agreement with conclusions based on micro wave, fluorescence, and multiphoton ionization experiments.     

FT-IR and FT-Raman investigation, computed vibrational intensity analysis and computed vibrational frequency analysis on m-Xylol using ab-initio HF and DFT calculations

The FT-IR and FT-Raman spectra of m-Xylol molecule have been recorded using Bruker IFS 66V spectrometer in the range 4000-100cm(-1). The molecular geometry and vibrational frequencies in the ground state are evaluated using the Hartree-fock (HF) and B3LYP with 6-31+G (d, p), 6-31++G (d, p) and 6-311++G (d, p) basis sets. The computed frequencies are scaled using a suitable scale factors to yield good agreement with the observed values. The HF and DFT analysis agree well with experimental observations. Comparison of the fundamental vibrational frequencies with calculated results by HF and B3LYP methods indicate that B3LYP/6-311++G (d, p) is superior to HF/6-31+G (d, p) for molecular vibrational problems. The complete data of this title compound provide some useful information for the study of substituted benzenes. The influences of Methyl groups on the geometry of benzene and its normal modes of vibrations have also been discussed.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

Quantum vs. classical models of the nitro group for proton chemical shift calculations and conformational analysis

A model based on classical concepts is derived to describe the effect of the nitro group on proton chemical shifts. The calculated chemical shifts are then compared to ab initio (GIAO) calculated chemical shifts. The accuracy of the two models is assessed using proton chemical shifts of a set of rigid organic nitro compounds that are fully assigned in CDCl3 at 700 MHz. The two methods are then used to evaluate the accuracy of different popular post-SCF methods (B3LYP and MP2) and molecular mechanics methods (MMX and MMFF94) in calculating the molecular structure of a set of sterically crowded nitro aromatic compounds. Both models perform well on the rigid molecules used as a test set, although when using the GIAO method a general overestimation of the deshielding of protons near the nitro group is observed. The analysis of the sterically crowded molecules shows that the very popular B3LYP/6-31G(d,p) method produces very poor twist angles for these, and that using a larger basis set [6-311++G(2d,p)] gives much more reasonable results. The MP2 calculations, on the other hand, overestimate the twist angles, which for these compounds compensates for the deshielding effect generally observed for protons near electronegative atoms when using the GIAO method at the B3LYP/6-311++G(2d,p) level. The most accurate results are found when the structures are calculated using B3LYP/6-311++G(2d,p) level of theory, and the chemical shifts are calculated using the CHARGE program based on classical models.     

Conformational stability, vibrational assignmenents, barriers to internal rotations and ab initio calculations of 2-aminophenol (d 0 and d3)

The Raman (3700-100 cm(-1)) and infrared (4000-400 cm(-1)) spectra of solid 2-aminophenol (2AP) have been recorded. The internal rotation of both OH and NH2 moieties produce ten conformers with either Cs or C1 symmetry. However, the calculated energies as well as the imaginary vibrational frequencies reduce rotational isomerism to five isomers. The molecular geometry has been optimized without any constraints using RHF, MP2 and B3LYP levels of theory at 6-31G(d), 6-311+G(d) and 6-31++G(d,p) basis sets. All calculations predict 1 (cis; OH is directed towards NH2) to be the most stable conformation except RHF/6-31++G(d,p) basis set. The 1 (cis) isomer is found to be more stable than 8 (trans; OH is away from the NH2 moiety and the NH bonds are out-of-plane) by 1.7 kcal/mol (598 cm(-1)) as obtained from MP2/6-31G(d) calculations. Aided by experimental and theoretical vibrational spectra, cis and trans 2AP are coexist in solution but cis isomer is more likely present in the crystalline state. Aided by MP2 and B3LYP frequency calculations, molecular force fields, simulated vibrational spectra utilizing 6-31G(d) basis set as well as normal coordinate analysis, complete vibrational assignments for HOC6H4NH2 and DOC6H4ND2 have been proposed. Furthermore, we carried out potential surface scan, to determine the barriers to internal rotations of NH2 and OH groups. All results are reported herein and compared with similar molecules when appropriate.     

Comparison of some representative density functional theory and wave function theory methods for the studies of amino acids

Energies of different conformers of 22 amino acid molecules and their protonated and deprotonated species were calculated by some density functional theory (DFT; SVWN, B3LYP, B3PW91, MPWB1K, BHandHLYP) and wave function theory (WFT; HF, MP2) methods with the 6-311++G(d,p) basis set to obtain the relative conformer energies, vertical electron detachment energies, deprotonation energies, and proton affinities. Taking the CCSD/6-311++G(d,p) results as the references, the performances of the tested DFT and WFT methods for amino acids with various intramolecular hydrogen bonds were determined. The BHandHLYP method was the best overall performer among the tested DFT methods, and its accuracy was even better than that of the more expensive MP2 method. The computational dependencies of the five DFT methods and the HF and MP2 methods on the basis sets were further examined with the 6-31G(d,p), 6-311++G(d,p), aug-cc-pVDZ, 6-311++G(2df,p), and aug-cc-pVTZ basis sets. The differences between the small and large basis set results have decreased quickly for the hybrid generalized gradient approximation (GGA) methods. The basis set convergence of the MP2 results has been, however, very slow. Considering both the cost and the accuracy, the BHandHLYP functional with the 6-311++G(d,p) basis set is the best choice for the amino acid systems that are rich in hydrogen bonds.