X-H (X = C, N, O, Si, P, S) 键离解能计算的密度泛函方法研究

使用了不同密度泛函方法计算X-H (X = C, N, O, Si, P, S) 键离解能，并分析不同密度泛函方法的计算精度。研究发现大多数密度泛函方法包括B3LYP, B3P86, B3PW91, G96LYP, PBE1PBE，和BH&HLYP都明显低估键离解能13-25 kJ/mol。该现象与是否使用无限基组无关，因为即使使用无限基组键离解能仍然被低估。因此密度泛函方法不适合用于键离解能的估算。其中B3P86方法的偏差最小。进一步分析表明，使用限制性开壳层计算并无任何优势，在大多数情况下非限制性开壳层计算实际上比限制性开壳层计算要好。最后，我们发现了密度泛函方法对键离解能的低估是系统的，因此建议利用校准后的UDFT/6-311++G(d, p)方法计算化学键离解能。     

Bond Dissociation Energies and Electronic Structures in a Series of Peroxy Radicals: A Theoretical Study

Quantum chemical calculations are performed to estimate the bond dissociation energies (BDEs) for 18 peroxy radicals. Since DFT methods are researched to have low basis sets sensitivity, these radicals are studied by utilizing the hybrid density functional theory (DFT) (B3LYP, B3P86, B3PW91 and PBE1PBE) in conjunction with the 6‐311G** basis set and the complete basis set (CBS‐Q) method. On the basis of comparisons of the computational results and the experimental values, we evaluate the effectiveness of above methods. It is demonstrated that CBS‐Q method is the best method for computing the reliable BDEs of C—OO bond, with the average absolute errors of 2.1 kcal/mol. So CBS‐Q method is suitable to predict accurate BDEs of C‐OO bond for peroxy compounds. The computational energy gaps between the HOMO and LUMO of studied compounds are almost identical from the point of view of stability and substantial HOMO‐LUMO gaps for all molecules suggest their electronic stability. In addition, substituent effect on the C—OO BDE of peroxy radicals is analyzed. It is noted that the effects of substitution on the C—OO BDE of peroxy radicals are significant. Our results will shed lights on future theoretical and experimental work.     

Intramolecular hydrogen bond in the hydroxycyclohexadienyl peroxy radicals

The hydroxycyclohexadienyl peroxy radicals (HO? C₆H₆? O₂) produced from the reaction of OH‐benzene adduct with O₂ were studied with density functional theory (DFT) calculations to determine their characteristics. The optimized geometries, vibrational frequencies, and total energies of 2‐hydroxycyclohexadienyl peroxy radical II_s and 4‐hydroxycyclohexadienyl peroxy radical III_s were calculated at the following theoretical levels, B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), and B3LYP/6‐311+G(d,p). Both were shown to contain a red‐shifted intramolecular hydrogen bond (O? H … O? H bond). According to atoms‐in‐molecules (AIM) analysis, the intramolecular hydrogen bond in the 2‐hydroxycyclohexadienyl peroxy radical II_s is stronger than that one in 4‐hydroxycyclohexadienyl peroxy radical III_s, and the former is the most stable conformation among its isomers. Generally speaking, hydrogen bonding in these radicals plays an important role to make them more stable. Based on natural bond orbital (NBO) analysis, the stabilization energy between orbitals is the main factor to produce red‐shifted intramolecular hydrogen bond within these peroxy radicals. The hyperconjugative interactions can promote the transfer of some electron density to the O? H antibonding orbital, while the increased electron density in the O? H antibonding orbital leads to the elongation of the O? H bond and the red shift of the O? H stretching frequency. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical Study of Remote Substituent Effects on X-H （X=CH2, NH,O） Bond Dissociation Energies of Azulene

In the study, the X-H （X=CH2, NH, O） bond dissociation energies （BDE） of para-substituted azulene （Y-C10H8X-H） were predicted theoretically for the first time using Density Functronal Theory （DFT） methods at UB3LYP/6-311 ＋ ＋g（2df,2p）//UB3LYP/6-31 ＋g（d） level. It was found that the substituents exerted similar effects on the X-H BDE of azulene as those on benzene, except for 6-substituted 2-methylazulene. Owing to the substituent-dipole interaction, the reaction constants （ρ^＋） of 2- and 6-Y-CIoHsX-H （X=NH and O only） varied violently. The origin of the substituent effects on the X-H BDE of azulene was found, by both GE/RE and SIE theory, to be directly associated with variation of the radical effects, although the ground effects also played a modest role in determining the net. substituent effects.     

First Total Synthesis of (±)‐Latifolin and Its Antioxidant Mechanism

The first total synthesis of (±)‐latifolin has been accomplished in six steps and 47.8% overall yield. To understand the relative importance of phenolic O? H and benzhydryl C? H hydrogen on the antioxidant activity of latifolin, 2,2‐diphenyl‐1‐picrylhydrazyl (DPPH) radical scavenging assay and density functional theory (DFT) studies were carried out. On scavenging DPPH radical in ethanol, the activity of latifolin ( 1 ) bearing phenolic hydrogen is remarkably higher than analogue 10 bearing no phenolic hydrogen. Therefore, Phenolic hydrogen is responsible for latifolin's antioxidant activity rather than benzhydryl C? H hydrogen. Furthermore, the 5‐OH BDE is lower than 2′‐OH and 7‐CH BDEs by a DFT calculation, respectively. Based on theoretical results it is definitely concluded that the phenolic 5‐OH plays a major role in the antioxidant activity of latifolin.     

A theoretical investigation on the structures,densities, detonation properties,and pyrolysis mechanism of the nitro derivatives of phenols

The nitro derivatives of phenols are optimized to obtain their molecular geometries and electronic structures at the DFT‐B3LYP/6‐31G* level. Detonation properties are evaluated using the modified Kamlet–Jacobs equations based on the calculated densities and heats of formation. It is found that there are good linear relationships between density, detonation velocity, detonation pressure, and the number of nitro and hydroxy groups. Thermal stability and pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies (BDEs) at the unrestricted B3LYP/6‐31G* level. The activation energies of H‐transfer reaction is smaller than the BDEs of all bonds and this illustrates that the pyrolysis of the title compounds may be started from breaking O? H bond followed by the isomerization reaction of H transfer. Moreover, the C? NO₂ bond with the smaller bond overlap population and the smaller BDE will also overlap may be before homolysis. According to the quantitative standard of energetics and stability as a high‐energy density compound, pentanitrophenol essentially satisfies this requirement. In addition, we have discussed the effect of the nitro and hydroxy groups on the static electronic structural parameters and the kinetic parameter. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Two Phenolic Diterpenes as Antioxidants

1 INTRODUCTION Experimental studies in recent years show thatphenolic diterpenes are potent antioxidants with highbio-activities of anti-tumor, anti-HIV and anti-mi-crobes[1～4], and they have been developed and ap-plied in medicine, cosmetics, and keeping foods fresh.However, the exact mechanism in which phenolicditerpenes act as those bio-activities are still unex-plored. A series of compounds with similar struc-tures, including carnosol, rosmanol, isorosmanol andso on, were isolated f…     

Theoretical Elucidation on Different Lipid-Oxidation Potentials of Aminoxyl Antioxidants

To elucidate the different lipid-oxidation potentials of aminoxyl antioxidants, a kind of combined density functional theory (DFT) method was employed to calculate C-H bond dissociation enthalpies (BDEs) of a model linoleic acid (LH) and O-H BDEs of hydrogenated aminoxyls. The higher the O-H BDE is, the more potent the amlnoxyl to abstract the H-atom from LH and the stronger the LH-oxidation potential. Accordingly,the prooxidant activity differences of amlnoxyls were elucidated by the different O--H BDEs of hydrogenated amlnoxyls, which were further clarified in terms of distinct electronic effects of the substituents.     

Evaluation of a Combined Quantum Chemical Method Used in Calculating O?H Bond Dissociation Enthalpy

O? H bond dissociation enthalpies (BDE) for a variety of substituted phenols were calculated using a combined quantum chemical method. It is found that the calculated O? H BDE correlated well with the recommended values, except for ortho‐tert‐butyl substituted phenols. For the electron‐donating group substituted phenols the calculated O? H BDE are slightly higher than the recommended values, however, for the electron‐withdrawing group substituted phenols the calculated O? H BDE are slightly lower than the recommended values.     

Theoretical investigation of the structures and properties of fluoromethyl peroxyl radicals

The equilibrium geometries, harmonic vibrational frequencies, charge distributions, spin density distributions, dipole moments, electron affinities (EAs), and C? O bond dissociation energies (BDEs) of HO, CH₃O, CH₂FO, CHF₂O, and CF₃O peroxyl radicals have been calculated using ab initio molecular orbital theory and density functional theory (DFT) at the B3LYP level. The C? H bond dissociation energies of the parent fluoromethanes have been calculated using the same levels of theory. Both the MP2(full) and B3LYP methods, using the 6‐31G(d,p) basis set, are found to be capable of accurately predicting the geometries of peroxyl radicals. Electron correlation accounts for ～25% of the C? H BDE of fluoromethanes and for ～50% of the C? O BDE of the corresponding peroxyl radicals. The B3LYP/6‐31G(d,p) method is found to be comparable to high ab initio levels in predicting C? O BDEs of studied peroxyl radicals and C? H BDEs of the parent alkanes. The progressive fluorine substitution of hydrogen atoms in methyl peroxyl radicals results in shortening of the C? O bond, lengthening of the O? O bond, an increase (decrease) of the spin density on the terminal (inner) oxygen, a decrease in the dipole moments, and an increase in electron affinities. Both C? O BDEs and EAs of peroxyl radicals (RO) correlate well with Taft σ* substituent constants for the R group in peroxyl radicals. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Analysis of pair interaction in a bipolar mesogen 4‐(4‐hydroxylbutyloxy)‐4′‐cyano‐biphenyl: A comparative study based on semiempirical and DFT methods

Structure of 4‐(4‐hydroxylbutyloxy)‐4′‐cyano‐biphenyl (H4CBP) molecule has been optimized using density functional B3LYP with 6‐31G (d) basis set taking crystallographic geometry as input. Using the optimized geometry, electronic structure of the H4CBP molecule has been evaluated on the basis of semiempirical methods and DFT calculations. Intermolecular interaction energy between a pair of H4CBP molecules has been evaluated by using Rayleigh–Schrodinger perturbation theory modified with multicentered multipole expansion method for the electrostatic part while dispersion and repulsion terms have been calculated using Kitaigorodskii formula. The results obtained through semiempirical and DFT calculations have been compared for various interacting conditions, viz.: (a) stacking, (b) in‐plane, and (c) terminal interactions. A comparative analysis of the results has been carried out with a view to examine suitability of different methods to study molecular aggregations in moderately large organic systems. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Thermochemical properties for n‐propyl,iso‐propyl,and tert‐butyl nitroalkanes,alkyl nitrites,and their carbon‐centered radicals

Density functional theory (DFT) based calculations are performed on a series of alkyl nitrites and nitroalkanes representing large‐scale primary, secondary, and tertiary nitro compounds and their radicals resulting from the loss of their skeletal hydrogen atoms. Geometries, vibration frequencies, and thermochemical properties [S°(T) and C°_p(T) (10 K ? T ? 5000 K)] are calculated at the B3LYP/6‐31G(d,p) DFT level. Δ_fH°₂₉₈ values are from B3LYP/6‐31G(d,p), B3LYP/6‐31+G(2d,2p), and the composite CBS‐QB3 levels. Potential energy barriers for the internal rotations have been computed at the B3LYP/6‐31G(d,p) level of theory, and the lower barrier contributions are incorporated into entropy and heat capacity data. The standard enthalpies of formation at 298 K are evaluated using isodesmic reaction schemes with several work reactions for each species. Recommended values derived from the most stable conformers of respective nitro‐ and nitrite isomers include ?30.57 and ?28.44 kcal mol^?1 for n‐propane‐, ?33.89 and ?32.32 kcal mol^?1 for iso‐propane‐, ?42.78 and ?41.36 kcal mol^?1 for tert‐butane‐nitro compounds and nitrites, respectively. Entropy and heat capacity values are also reported for the lower homologues: nitromethane, nitroethane, and corresponding nitrites. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 181–199, 2010     

A Theoretical Study of Ene Reactions in Solution: A Solution‐Phase Translational Entropy Model

Several density functional theory (DFT) methods, such as CAM‐B3LYP, M06, ωB97x, and ωB97xD, are used to characterize a range of ene reactions. The Gibbs free energy, activation enthalpy, and entropy are calculated with both the gas‐ and solution‐phase translational entropy; the results obtained from the solution‐phase translational entropies are quite close to the experimental measurements, whereas the gas‐phase translational entropies do not perform well. For ene reactions between the enophile propanedioic acid (2‐oxo‐1,3‐dimethyl ester) and π donors, the two‐solvent‐involved explicit+implicit model can be employed to obtain accurate activation entropies and free‐energy barriers, because the interaction between the carbonyl oxygen atom and the solvent in the transition state is strengthened with the formation of C?C and O?H bonds. In contrast, an implicit solvent model is adequate to calculate activation entropies and free‐energy barriers for the corresponding reactions of the enophile 4‐phenyl‐1,2,4‐triazoline‐3,5‐dione.     

PCM study of the solvent and substituent effects on bond dissociation energies of the O?NO Bond—A DFT study

A PCM continuum model, at the B3LYP, B3P86, and B3PW91 three‐parameter hybrid DFT methods with 6‐311G** basis set, is used to study the bond dissociation energies (BDEs) of benzyl nitrites. Compared the computed results with the experimental values, it is noted that B3PW91 functional is the best method to compute the BDEs of benzyl nitrites. The solvent and substituent effects on the BDEs of the O? NO bond are analyzed, and it is shown that the BDE of the O? NO bond decreases with the increment of the Hammett constants of substituent groups on benzene for benzyl nitrites except C₆H₅CH₂O? NO. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Quantum chemical computations of 1,3-phenylenediacetic acid

Molecular geometry, vibrational wavenumbers and gauge including atomic orbital (GIAO) 13C NMR and 1H NMR chemical shift values of 1,3-phenylenediacetic acid (C10H10O4), in the ground state have been calculated by using ab initio Hartree-Fock (HF) and density functional theory (DFT/B3LYP) methods with 6-311++G(d,p) basis set for the first time. Comparison of the observed fundamental vibrational modes of 1,3-phenylenediacetic acid and calculated results by density functional B3LYP and Hartree-Fock methods indicate that B3LYP is superior to the scaled Hartree-Fock approach for quantum chemical studies. Geometric parameters (bond lengths and bond angles) and vibrational wavenumbers obtained by the HF and DFT/B3LYP methods are in good agreement with the experimental data. Furthermore, this is the first time the results of the calculated JCH and JCC coupling constants of the C10H10O4 molecule are presented in this study.     

Crystal and geometry‐optimized structure of an anthracene‐based Diels–Alder adduct

Computational calculations of an anthracene‐based Diels–Alder adduct, namely, 17‐ethyl‐1‐hydroxymethyl‐17‐azapentacyclo[6.6.5.0^2,7.0^9,14.0^15,19]nonadeca‐2,4,6,9,11,13‐hexaene‐16,18‐dione, C₂₁H₁₉NO₃, predicting density functional theory (DFT) optimized geometries in the gas phase are compared in terms of accuracy relative to the solid‐state crystal structure and computational cost. Crystal structure determination and Hirshfeld surface analysis of the racemic product reveal that the molecules are linked by O—H…O=C hydrogen bonds between the hydroxy and carbonyl groups, accounting for 9.5% of the intermolecular contacts, while H…H contacts represent 56.9% of the total. Boltzmann population analysis of computed relative rotamer abundances in the gas phase are based on lower‐level geometry optimization and thermochemical corrections coupled with higher‐level electronic energy calculations using the B2PLYP double‐hybrid functional. As expected, the choice of density functional has a greater effect than the basis set on accuracy for all levels of theory. For any given functional, increasing the basis set size did not always correlate with increasingly accurate structures. The hybrid functional B3LYP without dispersion correction routinely gave the closest approximations to the crystal structure where the B3LYP/aug‐cc‐pVDZ combination afforded the best structure (r.m.s. deviation = 0.1314 Å). However, the B3LYP/6‐31+G(d,p) level of theory represents the best compromise between accuracy (r.m.s. deviation = 0.1388 Å) and cost as it yielded appreciably accurate results in a fraction of the time compared to other method combinations.     

Novel silylphenol antioxidants by density functional theory

We have scrutinized five novel silylphenol antioxidants, including 2-silylphenol ( 1 ), 4-silylphenol ( 2 ), 2,6-disilylphenol ( 3 ), 2,4-disilylphenol ( 4 ), and 2,4,6-trisilylphenol ( 5 ), at M06/6–311++G** level of theory. To evaluate the antioxidant efficiency, the electronic effects on O─H bond dissociation energy (BDE) and vertical ionization potential (IP_v) of 1 – 5 are investigated, which are mainly governed by electronic effects. The conductor-like polarized continuum model (CPCM) is applied to measure the antioxidant capacity in the solution phase. The results show that antioxidants with the lowest BDE and IP_v values can efficiently act via hydrogen atom transfer (HAT) and single electron transfer (SET) mechanisms, respectively. The stability of resulting radicals is measured by nucleus independent chemical shift (NICS) index, natural bond orbital (NBO) analysis, and nucleophilicity (N) index. The BDE shows lower values in the gas phase with respect to water, while water exhibits lower IP_v values than gas. Structure 5 turns out as the most efficient antioxidant. The overall order of antioxidant efficiency in both gas and water phases is 5 > 2 > 3 > 4 > 1 .     

Theoretical studies on the structures, thermodynamic properties, detonation properties, and pyrolysis mechanisms of spiro nitramines

Density function theory (DFT) has been employed to study the geometric and electronic structures of a series of spiro nitramines at the B3LYP/6-31G level. The calculated results agree reasonably with available experimental data. Thermodynamic properties derived from the infrared spectra on the basis of statistical thermodynamic principles are linearly correlated with the number of nitramine groups as well as the temperature. Detonation performances were evaluated by the Kamlet-Jacobs equations based on the calculated densities and heats of formation. It is found that some compounds with the predicted densities of ca. 1.9 g/cm3, detonation velocities over 9 km/s, and detonation pressures of about 39 GPa (some even over 40 GPa) may be novel potential candidates of high energy density materials (HEDMs). Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6-31G level and the activation energies (E(a)) with the selected PM3 semiempirical molecular orbital (MO) based on the unrestricted Hartree-Fock model. The relationships between BDE, E(a), and the electronic structures of the spiro nitramines were discussed in detail. Thermal stabilities and decomposition mechanisms of the title compounds derived from the B3LYP/6-31G BDE and the UHF-PM3 E(a) are basically consistent. Considering the thermal stability, TNSHe (tetranitrotetraazaspirohexane), TNSH (tetranitrotetraazaspiroheptane), and TNSO (tetranitrotetraazaspirooctane) are recommended as the preferred candidates of HEDMs. These results may provide basic information for the molecular design of HEDMs.     

The structure of cinnamic acid and cinnamoyl azides,a unique localized π system: The electronic spectra and DFT‐treatment

The electronic absorption spectra of cinnamic acid and some cinnamoyl azides have been recorded in absolute methanol and investigated to explore the structure of the titled compounds. Cinnamic acid and its derivatives have a double bond, ? C?C? , between the aromatic ring and the carboxyl group which disturbs the π electron system of the molecule and inhibits electron delocalization as compared with styrene or benzoic acid. The azide group is neither a strong electron donor nor a strong electron acceptor but it increases conjugation in the molecule. The observed spectra confirm that each of the cinnamic acid and cinnamoyl azide molecules is one of a kind of unique disturbed π‐system and not of different independent π systems, each on a fragment of the molecule as predicted by the quantum theory of atom in molecule calculations. The spectra of cinnamic acid and its derivatives are not the additive spectra of the different fragments of the molecule. The spectra are characterized by few number, low intensity, and high‐energy electronic transitions (absorption bands) in the UV‐vis region. Molecular orbital calculations confirmed the spectral observations. The optimized geometry of the ground state of the studied compounds is calculated using the DFT/B3LYP/6‐31G** level of theory and an explicit molecular orbital analysis is carried out. Excited states are calculated using the TD/DFT procedure as implemented by the Gamess 2009 package of programs. The correspondence between calculated and the observed transition energies is adequate. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

The C-H and alpha(C-X) bond dissociation enthalpies of toluene, C6H5-CH2X (X = F, Cl), and their substituted derivatives: a DFT study

The homolytic C-H bond dissociation enthalpies (BDEs) of toluene and its para- and meta-substituted derivatives have been estimated by using the (RO)B3LYP/6-311++G(2df,2p)//(U)B3LYP/6-311G(d,p) procedure. The performance of two other hybrid functionals of DFT, namely, B3PWP91 and O3LYP, has also been evaluated using the same basis sets and molecules. Our computed results are compared with the available experimental values and are found to be in good agreement. The (RO)B3LYP and (RO)O3LYP procedures are found to produce reliable BDEs for the C-H bonds in toluene and the C-X (X = F, Cl) bond in alpha-substituted toluene (C6H5-CH2X) and their substituted derivatives. The substituent effect on the BDE values has been analyzed in terms of the ground-state effect and the radical effect. The effect of polarization of the C-H bond on the substituent effect is also analyzed. The BDE(C-H) and BDE(C-X) values for alpha-substituted (X = F and Cl) toluenes with a set of para substituents are presented for the first time.