Density functional theory study on molecular structure and vibrational IR spectra of isocytosine

Density functional theory with the combined Becke3-LYP exchange-correlation energy functional [DFT(B3-LYP) method] using the 6-31G(d, p) basis set is applied to predict molecular parameters (geometries, rotational constants, dipole moments) and vibrational IR spectra (harmonic wavenumbers, absolute intensities) of six tautomers of the isocytosine molecule. The results are compared with the corresponding data calculated at the conventional ab initio Hartree-Fock (HF) level using the same basis set and with available experimental data. Calculations show that (a) three amino tautomers are slightly nonplanar species with, evidently, a distorted amino group, (b) the DFT (B3-LYP)/6-31G(d, p) method predicts better molecular parameters, than do the HF calculations, and (c) the DFT(B3-LYP)-calculated vibrational IR spectra of isocytosine agree well with the available recorded IR spectra, and they show marked improvement over the IR spectra predicted at the HF/6-31G(d, p) level. Tautomeric stabilities of isocytosine are discussed on the basis of computed electronic energies by the DFT(B3-LYP) and ab initio approaches [including the MP2 and MP4(SDQ) calculations of electronic energies] and predicted zero-point vibrational energies by DFT(B3-LYP) and HF methods. This relative energies at 0 K of the tautomeric forms of isocytosine predicted by both conventional ab initio and DFT(B3-LYP) methods correlate well with the experimental data, showing the predominance of the aminohydroxy tautomer of isocytosine for an isolated molecule. © 1997 John Wiley & Sons, Inc.     

Ring conserved isodesmic reactions: A new method for estimating the heats of formation of aromatics and PAHs

Density functional theory (DFT) has been used along with isodesmic reaction schemes to estimate heats of formation for aromatics and polynuclear aromatic hydrocarbons (PAHs). Calculations have been performed for 42 molecules, 12 of which have uncertain or unknown experimental values, using the B3-LYP functional with the small 6-31G(d) basis set. Heats of formation for the group of test molecules were estimated using both conventional bond separation (BS) isodesmic reactions as well as a new technique of ring conserved (RC) isodesmic reactions which is able to correct systematic errors in B3-LYP calculations. When a ring conserved isodesmic reaction based on delocalization energies is used, the estimated heat of formation is more accurate than that obtained by the bond separation technique. The methodology for creating and using appropriate ring conserved isodesmic reactions is discussed. The present scheme also compares favorably against a recently developed bond centered group additivity scheme that was tested against a large number of PAH molecules.     

Cyanovinyl radical: an illustration of the poor performance of unrestricted perturbation theory and density functional theory procedures in calculating radical stabilization energies

Stabilization energies for the 1-cyanovinyl radical (CH₂=CCN) have been calculated using a variety of conventional ab initio (M?ller–Plesset, quadratic configuration interaction and coupled-cluster) and density functional theory (B-LYP, B3-LYP) procedures, as well as with a range of compound methods. Compared with a high-level benchmark value (that predicts a stabilization energy of 17.1 kJ mol⁻¹), UMP2 and UMP4 give the wrong sign and magnitude of the stabilization energy (both methods predicting desta- bilization instead of stabilization), while B-LYP and B3-LYP overestimate the degree of stabilization. The RMP2, RMP4, QCISD(T) and CCSD(T) techniques, and several, but not all, variants of G2 and CBS theories give radical stabilization energies in good agreement with the benchmark value. Received: 15 June 1998 / Accepted: 19 August 1998 /  Published online: 15 February 1999     

The performance of B3-LYP density functional theory in describing SN2 reactions at saturated carbon

The performance of the B3-LYP variant of density functional theory when used in conjunction with the 6-31G(d) and 6-311 + G(3df, 2p) basis sets in describing the prototypical gas-phase S_N2 reactions of Cl⁻ + CH₃Cl and CH₃Br has been examined in detail. Reasonable values of the complexation energies (ΔH_comp) for the ion-molecule complexes formed in these reactions are obtained. However, the overall (ΔH_ovr^#) and central (ΔH_cent^#) barriers for these reactions calculated using the B3-LYP functional are significantly underestimated when compared with G2(+) or experimental results. This implies that the B3-LYP energies for the Cl(H₃C)Cl⁻ (D_3h) and Cl(H₃C)Br⁻ (C_3v) transition structures are relatively too low. The B3-LYP errors appear to be systematic, with similar errors being found for corresponding quantities for the two reactions examined.     

Density functional calculations of the pseudorotational flexibility of tetrahydrofuran

The flexibility of the five-membered ring in tetrahydrofuran was investigated using quantum mechanical methods involving density functional, Hartree-Fock, and many-body perturbation theory (MP2, MP4) calculations. We found that motion along the pseudorotational path of tetrahydrofuran is nearly free. The 0.1 kcal/mol energy barrier for pseudorotation, calculated at the highest MP4(SDQ)/6-311++G(2d,p)//MP2/6-311++G(2d,p) level of theory, agrees well with the experimental value of 0.16±0.03 kcal/mol. Similar results were obtained with the S-VWN, B3-LYP and B-LYP density functional calculations using the 6-31G(d) set of atomic orbitals. Also the density functional dipole moments and geometries were in good agreement with both the MP2 and experimental benchmarks. However, all density functional methods that utilized the default integration grid in the Gaussian 94 program were found to provide false stationary points of the C ₁ symmetry near the pseudorotational angle of 100°. These stationary points disappeared when a denser spherical-product grid was used. Overall, the hybrid B3-LYP functional was found to be the most promising quantum mechanical method for the modeling of biomolecules containing the furanose ring. Received: 17 June 1997 / Accepted: 20 November 1997     

Testing the performance of density functionals for the calculation of energetic properties of complex-forming radical-molecule reactions

The performance of some popular and some more recent density functional methods for the calculation of energies of stationary points on the potential surfaces of radical-molecule reactions was examined. The functionals studied are B3-LYP, BH&H, BH&H-LYP, MPW1K, MPWB1K, TPSS, TPSSh, BB1K, M05 and M05-2X, in conjunction with nine different AO basis sets. The reaction energies, barrier heights and the relative energies of the pre-and post-reaction complexes were compared with those obtained at the CCSD(T)/CBS limit for the reactions of OH radicals with HOOH and CH₃OOH. Very poor barrier heights are provided by the B3-LYP, TPSS and TPSSh functionals. The best overall performance was obtained with the BB1K, MPW1K and MPWB1K functionals. In these reactions all of the studied functionals provide converged results only if they are used with large basis sets like aug-cc-pVTZ and def2-TZVP. The data show that before relying on a functional for a specific reaction, it is desirable to make some test calculations on the performance. The same functional can predict some relative energies very well and some others very poorly even in systems including chemically similar reactants.     

A Density Functional Study on the Hydrolysis Process of Non-classical Transplatin（Ⅱ） with Two Same Planar Heterocycle Amines

In the present work, the hydrolysis process of non-classical transplatin（Ⅱ） with two same planar heterocycle amines has been studied using hybrid density functional theory （B3LYP） and IEF-PCM solvation models. Optimizations were performed at the B3LYP level using a combined basis set of （LanL2DZ＋6-31 ＋G（d,p）） with single-point energy evaluations using the B3LYP/6-31 ＋ ＋G（3df,2pd） approach in vacuo and in aqueous solution. For the obtained structures of reactants, intermediates, transition states, and products, both thermodynamic （reaction energies and Gibbs energies） and kinetic （reaction barriers） characteristics were estimated. In comparison with cisplatin, decreased activation energies were obtained. The result implies that the non-classical transplatin with two same planar heterocycle amines increases the equatorial steric effect and lowers reaction barriers, which may assist in designing novel Pt-based anticancer drugs.     

How to compute isomerization energies of organic molecules with quantum chemical methods

The reaction energies for 34 typical organic isomerizations including oxygen and nitrogen heteroatoms are investigated with modern quantum chemical methods that have the perspective of also being applicable to large systems. The experimental reaction enthalpies are corrected for vibrational and thermal effects, and the thus derived "experimental" reaction energies are compared to corresponding theoretical data. A series of standard AO basis sets in combination with second-order perturbation theory (MP2, SCS-MP2), conventional density functionals (e.g., PBE, TPSS, B3-LYP, MPW1K, BMK), and new perturbative functionals (B2-PLYP, mPW2-PLYP) are tested. In three cases, obvious errors of the experimental values could be detected, and accurate coupled-cluster [CCSD(T)] reference values have been used instead. It is found that only triple-zeta quality AO basis sets provide results close enough to the basis set limit and that sets like the popular 6-31G(d) should be avoided in accurate work. Augmentation of small basis sets with diffuse functions has a notable effect in B3-LYP calculations that is attributed to intramolecular basis set superposition error and covers basic deficiencies of the functional. The new methods based on perturbation theory (SCS-MP2, X2-PLYP) are found to be clearly superior to many other approaches; that is, they provide mean absolute deviations of less than 1.2 kcal mol-1 and only a few (<10%) outliers. The best performance in the group of conventional functionals is found for the highly parametrized BMK hybrid meta-GGA. Contrary to accepted opinion, hybrid density functionals offer no real advantage over simple GGAs. For reasonably large AO basis sets, results of poor quality are obtained with the popular B3-LYP functional that cannot be recommended for thermochemical applications in organic chemistry. The results of this study are complementary to often used benchmarks based on atomization energies and should guide chemists in their search for accurate and efficient computational thermochemistry methods.     

Density functional theory rationalization of the substituent effects in trifluoromethyl-pyridinol derivatives

The influence of α-substitution in the structure, bonding and thermochemical properties of trifluoromethyl-pyridinol derivatives and their protonated counterparts has been studied by means of density functional theory. The geometries of the neutral and protonated species were optimized at the B3-LYP/6-311G(d,p) level of theory. Final energies were obtained through single point B3-LYP/6-311+G(3df,2p) calculations.The relative orientation of the different substituents within the heterocycle ring favours the formation of unexpected intramolecular hydrogen bonds (IHB), which have been characterized by means of the Atoms in Molecules theory of Bader. Although weak, these IHB are of great importance for understanding the gas phase structure and the thermodynamical properties of these compounds. Surprisingly, most of the substituted investigated pyridinols present proton affinities below or close to that calculated for the unsubstituted pyridine molecule. Only pyridinols bearing strong σ or π donor activating groups show proton affinities greater than that of pyridine.     

Hammett correlations in the chemistry of 3-phenylpropyl radicals

The energetics and kinetics of the reaction of variously substituted benzyl radicals with a model alkene were calculated at the G3(MP2)-RAD//B3-LYP/6-31G(d) level of theory to determine whether such reactions are amenable to Hammett analysis. The reactions were studied both in the gas phase and in toluene solution in the temperature range 298-353 K; calculations include 1D-hindered rotor corrections for low frequency torsional modes, and the solvation energies were calculated using COSMO-RS at the BP/TZP level of theory. The addition reaction was found to be dominated by radical stabilization effects, but under circumstances where olefin substituent effects were decoupled from aryl substituent effects, a modest polar effect comes into play, which is enhanced by solvation. Reasonable correlations with empirical substituent parameters such as Hammett σ and σ(?) are observed for the enthalpy of activation, but additional entropic factors act to decrease the degree of correlation with respect to free energies and rate coefficients, confirming hypotheses from earlier experimental work. Substituent effects on the reverse β-fragmentation reaction, and potential cyclization of the 3-phenylpropyl radicals formed by addition are also discussed.     

N-甲基b-磺内酰胺氨解反应机理的理论研究

The aminolysis and the effect of water on the aminolysis processes of n-methyl β-sultam have been studied using density functional theory （DFF） method at the B3LYP/6-31G＊ level. The stationary structures and energies have been investigated for both reactions to find two different reaction channels. Specific and general solvent effects have been evaluated and the most favored pathway was found. The presence of solvent disfavors the reaction, whereas the participation of water in the aminolysis reaction plays a positive role and reduces the activation energy greatly. All transition states in the assisted aminolysis are 35-70 kJ/mol lower than those for the non-assisted reaction.     

Thermochemistry and kinetics of hydrogen abstraction by methyl radical from polycyclic aromatic hydrocarbons

Thermodynamic and kinetic properties relating to hydrogen abstraction by methyl radical from various sites in polycyclic aromatic hydrocarbons (PAHs) have been investigated. The reaction enthalpies (298 K), barriers (0 K), and activation energies and pre-exponential factors (700-1100 K), have been calculated by means of density functional theory, specifically with B3-LYP/6-311G(d,p) geometries, followed by BMK/6-311+G(3df,2p) single-point energy calculations. For uncongested sites in the PAHs, a reasonable correlation is obtained between reactivities (as characterized by the reaction barriers) and reaction enthalpies. This is reflected in a Bell-Evans-Polanyi (BEP) relationship. However, for congested sites, abstraction is accompanied both by lower reaction enthalpies (due to relief of steric strain) and also by reduced reactivities (due to significantly increased steric hindrance effects in the transition structures), so that the BEP relationship does not hold. In addition, the reaction enthalpies and kinetic parameters for the series of linear acenes indicate that abstraction is more difficult from the central rings.     

Electronically excited states and visible region photodissociation spectroscopy of Au(m)(+) . Ar(n) clusters (m=7-9): molecular dimensionality transition?

Photodissociation spectra were determined for Au(m)(+) . Ar(n) (m=7; n=0-3 and m=8,9; n=0,1) in the photon energy range of 2.14-3.02 eV. Experimental data were compared with predictions of dipole allowed transitions using time-dependent density functional theory (TDDFT) as applied to cluster structures from both DFT (B3-LYP functional) and ab initio calculations at the MP2 level. Argon adduct formation does not significantly perturb the bare metal cluster core structure, but it does change the metal cluster spectrum for highly symmetric cluster structures. The photodissociation spectra are consistent with a transition from planar to three-dimensional gold cluster core geometries between m=7 and m=8 for both n=0 and 1. TDDFT predictions for favored isomers describe experimental absorption features to within +/-0.25 eV. We also discuss size-dependent trends in TDDFT transition energies for the lowest energy two- and three-dimensional structures of Au(m)(+)(m=3-9).     

An evaluation of harmonic vibrational frequency scale factors

Scale factors for obtaining fundamental vibrational frequencies, low-frequency vibrational frequencies, zero-point vibrational energies (ZPVEs), and thermal contributions to enthalpy and entropy have been derived through a least-squares approach from harmonic frequencies determined at more than 100 levels of theory. Wave function procedures (HF, MP2, QCISD, QCISD(T), CCSD, and CCSD(T)) and a large and representative range of density functional theory (DFT) approaches (B3-LYP, BMK, EDF2, M05-2X, MPWB1K, O3-LYP, PBE, TPSS, etc.) have been examined in conjunction with basis sets such as 6-31G(d), 6-31+G(d,p), 6-31G(2df,p), 6-311+G(d,p), and 6-311+G(2df,p). The vibrational frequency scale factors were determined by a comparison of theoretical harmonic frequencies with the corresponding experimental fundamentals utilizing a standard set of 1066 individual vibrations. ZPVE scale factors were generally obtained from a comparison of the computed ZPVEs with experimental ZPVEs for a smaller standard set of 39 molecules, though the effect of expansion to a 48 molecule data set was also examined. In addition to evaluating the scale factors for a wide range of levels of theory, we have also probed the effect on scale factors of varying the percentage of incorporated exact exchange in hybrid DFT calculations using a modified B3-LYP functional. This has revealed a near-linear relationship between the magnitude of the scale factor and the proportion of exact exchange. Finally, we have investigated the effect of basis set size on HF, MP2, B3-LYP, and BMK scale factors by deriving values with basis sets ranging from 6-31G(d) up to 6-311++G(3df,3pd) as well as with basis sets in the cc-pVnZ and aug-cc-pVnZ series and with the TZV2P basis.     

High-level ab initio versus DFT calculations on (H2O2)2 and H2O2–H2O complexes as prototypes of multiple hydrogen bond systems

The performance of B-LYP, B-P86, B3-LYP, B3-P86, and B3-PW91 density functionals to describe multiple hydrogen bond systems was studied. For this purpose we have chosen the dimers of hydrogen peroxide and the hydrogen peroxide–water complexes. The geometries and vibrational frequencies obtained with a 6-311+G(d,p) basis set were compared with those obtained at the MP2 level using the same basis set expansion. The corresponding dimerization energies were obtained using a 6-311+G(3df,2p) basis set and compared with those obtained using the G2(MP2) theory. Red shiftings of the OH donor stretching frequencies were predicted by all approaches investigated; however, in all cases, the DFT values were sizably larger than the MP2 ones. Similarly, the blue shifting of the torsion of the hydrogen peroxide subunit was larger when evaluated at the DFT level. All functionals reproduced the G2(MP2) relative stabilities of the different local minima quite well. With the exception of the B-LYP and B3-PW91 approaches, all functionals yielded binding energies which deviated from the G2(MP2) values by less than 0.5 kcal/mol, provided that G2-type basis sets were used and that the corresponding BSSE corrections were included. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1124–1135     

线性BC2nB (n＝1～12)的结构特征和电子光谱的理论研究

应用密度泛函理论, 在B3LYP/6-31G^*水平上优化得到了线性簇合物BC_2nB (n＝1～12, D_(h)的平衡几何构型, 并计算了它们的谐振动频率. 在优化平衡几何构型下, 通过TD-B3LYP/cc-pvDZ和TD-B3LYP/cc-pvTZ计算, 分别得到了n＝1～12和n＝1～7的电子跃迁的垂直激发能和对应的振子强度. 在B3LYP/6-311＋G^*水平上计算得到了簇合物BC_2nB (n＝1～12, D_(h)的电离能. 基于计算结果, 导出了BC_2nB体系电子跃迁能以及第一电离能与体系大小n的解析表达式.     

Theoretical calculations of glycine and alanine gas-phase acidities

The gas-phase acidities of glycine and alanine were determined by using a variety of high level theoretical methods to establish which of these would give the best results with accessible computational efforts. MP2, MP4, QCISD, G2 ab initio procedures, hybrid Becke3-LYP (B3LYP) and gradient corrected Becke-Perdew (BP) and Perdew-Wang and Perdew (PWP) nonlocal density functionals were used for the calculations. A maximum deviation of approximately 13 and 18 kJ/mol from experimental data was observed for the computed delta Hacid and delta Gacid values, respectively. The best result was obtained at G2 level, but comparable reliability was reached when the considerably less time consuming B3LYP, BP, and PWP density functional approaches were employed.     

Internal rotation in peroxynitrous acid (ONOOH)

Using higher levels of wave-function-based electronic structure theory than previously applied, as well as density functional theory (B-LYP and B3-LYP functionals), all theoretical models conclude that three ONOOH conformers are stationary point minima, in disagreement with some of the previous studies that we survey. In order of increasing energy, these are the cis-cis, cis-perp, and trans-perp conformers. Basis sets including diffuse functions seem to be needed to obtain a qualitatively correct representation of the internal rotation potential energy surface at higher levels of theory. Internal rotation about the peroxide bond involving the cis-cis, cis-gauche transition structure (TS), cis-perp, and cis-trans TS conformers is studied in detail. To help ascertain the relative stability of the cis-perp conformer, multireference configuration interaction energy calculations are carried out, and rule of thumb estimates of multireference character in the ground-state wave functions of the ONOOH conformers are considered. CCSD(T)/aug-cc-pVTZ physical properties (geometries, rotational constants, electric dipole moments, harmonic vibrational frequencies, and infrared intensities) are compared with the analogous experimental data wherever possible, and also with density functional theory. Where such experimental data are nonexistent, the CCSD(T) and B3-LYP results are useful representations. For example, the electric dipole moment |mu(e)| of the cis-cis conformer is predicted to be 0.97+/-0.03 D. CCSD(T) energies, extrapolated to the aug-cc-pVNZ limit, are employed in isodesmic reaction schemes to derive zero Kelvin heats of formation and bond dissociation energies of the ONOOH stationary point minima. In agreement with recent gas-phase experiments, the peroxide bond dissociation energies of the cis-cis and trans-perp conformers are calculated as 19.3+/-0.4 and 16.0+/-0.4 kcalmol, respectively. The lowest energy cis-cis conformer is less stable than nitric acid by 28.1+/-0.4 kcalmol at 0 K.     

A direct ab initio dynamics study of the water-assisted tautomerization of formamide

Direct ab initio dynamics calculations based on a canonical variational transition-state theory with several multidimensional semiclassical tunneling approximations were carried out to obtain rate constants for the water-assisted tautomerization of formamide. The accuracy of the density functionals, namely, B-LYP, B3-LYP, and BH&H-LYP, were examined. We found that the BH&H-LYP method yields the most accurate transition-state properties when comparing it to ab initio MP2 and QCISD results, whereas B-LYP and B3-LYP methods predict barrier heights too low. Reaction path information was calculated at both the MP2 and nonlocal hybrid BH&H-LYP levels using the 6–31G(d,p) basis set. At the BH&H-LYP level, we found that the zero-point energy motion lowers the barrier to tautomerization in the formamide-water complex by 3.6 kcal/mol. When tunneling is considered, the activation energy at the BH&H-LYP level at 300 K is 17.1 kcal/mol. This is 3.4 kcal/mol below the zero-point-corrected barrier and 7.0 kcal/mol below the classical barrier. Excellent agreement between BH&H-LYP and MP2 rate constants further supports the use of BH&H-LYP for rate calculations of large systems. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 861–874, 1997