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.     

Theoretical isotropic hyperfine coupling constants of third-row nuclei (29Si, 31P, and 33S)

In a previous paper (Hermosilla, L.; Calle, P.; Garcia de la Vega, J. M.; Sieiro, C. J. Phys. Chem. A 2005, 109, 1114), an adequate computational protocol for the calculation of isotropic hyperfine coupling constants (hfcc's) was proposed. The main conclusion concerns the reliability of the scheme B3LYP/TZVP//B3LYP/6-31G* in the predictions of hfcc's with low computational cost. In the present study, we gain insight into the behavior of the above functional/basis set scheme on nuclei of the third row, for which few systematic studies have been carried out up to the present date. The systems studied are neutral, cationic, anionic, localized, and conjugated radicals, containing (29)Si, (31)P, and (33)S nuclei. After carrying out a regression analysis, we conclude that density functional theory (DFT) predictions on the hfcc's of the third-row nuclei are reliable for B3LYP/TZVP by using an optimized geometry with B3LYP/6-31G* combination. By comparison with other much more computationally demanding schemes, namely, B3LYP/cc-pVTZ and B3LYP/cc-pVQZ, we conclude that the B3LYP functional in conjunction with the TZVP basis set is the most useful computational protocol for the assignment of experimental hfcc's, not only for nuclei of first and second rows, but also for those of the third row.     

DFT study of the EPR spectral pattern of propagating methacrylic radicals

The EPR spectral pattern observed in the bulk polymerization of methacrylic monomers was theoretically investigated by DFT methodology. The conformational analysis of the propagating radical by the rotation around the C–C_β bond, was performed using the B3LYP/6-31G^* computational scheme. To obtain accurate values of the isotropic hyperfine coupling constants (hfccs) a higher level protocol, B3LYP/TZVP//B3LYP/6-31G^*, was applied. The experimental 13-line spectrum registered at the first stage of the polymerization was assumed to correspond to a free rotating radical in a fluid medium and it was simulated just considering the most stable conformation. The 9-line spectrum registered at high conversions was interpreted in terms of highly hindered rotational conformers frozen in the very viscous medium. This spectrum was well reproduced by a model which considers the sum of the individual spectra of the conformations spread around the most probable. Each of these contributing spectra was obtained on the basis of the computed hfccs for the considered conformations and weighted by his relative Boltzmann population according to the DFT analysis. Besides, the calculated hfccs showed an excellent agreement with those predicted by the Heller–McConnell relationship, which confirms the coherence of the DFT methodology for this kind of calculations.     

Assessment of density functional theory methods for the computation of heats of formation and ionization potentials of systems containing third row transition metals

The performance of several different density functional theory (DFT) methods, including GGA, hybrid-GGA, meta-GGA, and hybrid-meta-GGA methods, have been assessed in terms of their ability to accurately compute both heats of formation and ionization potentials of systems containing third row transition metals. Two different basis sets were used in this study: 6-31G** and TZVP. It is found that the triple-zeta quality TZVP basis set generally produces the best results for both heats of formation and ionization potentials. One important observation made in this study is that the inclusion of exact exchange terms in DFT methods generally results in more consistently accurate results for both heats of formation and ionization potentials of transition metal systems. In general, DFT methods do not yield good ionization potential results for systems containing titanium or zinc. For heats of formation, it is found that the hybrid-meta-GGA functional, TPSS1KCIS, yields the best overall results when combined with the TZVP basis set, while PBE1PBE (hybrid-GGA) gives the best results for the 6-31G** basis. The hybrid-GGA functional, B3LYP, is found to produce the lowest overall errors for ionization potentials when combined with both 6-31G** and TZVP.     

Interaction of guanine, its anions, and radicals with lysine in different charge states

Modification in DNA or protein structure can severely affect DNA-protein interactions and the functioning of biological systems. Some new insights into radiation-induced effects of guanine-lysine interactions have been obtained here by theoretical investigations. Geometries of zwitterionic and non-zwitterionic lysine in different charge states (neutral, radical cation, and protonated cation) were optimized employing the B3LYP/6-31G** and B3LYP/AUG-cc-pVDZ levels of hybrid density functional theory (DFT) and using the second-order M?ller-Plesset perturbation theory along with the 6-31G** basis set. In the case of neutral lysine in the gas phase, no zwitterionic structure was obtained. The non-zwitterionic structures of lysine in radical and protonated cationic forms are appreciably more stable than the corresponding zwitterionic structures in the gas phase as obtained at all levels of theory employed here. Binding of guanine and different dehydrogenated guanine radicals with lysine in different charge states was studied at the B3LYP/6-31G** level of DFT. When guanine makes a complex with the lysine radical cation, large amounts of spin and positive charge densities are transferred from the lysine radical cation to guanine and the guanine is thus converted from its normal form to the radical cationic form. Complexation of the lysine radical cation with the H1-hydrogen-abstracted guanine radical leads to CO2 liberation and proton transfer from lysine. These results are compared with the available experimental ones.     

平面五元水分子簇的量子化学研究

选用Gaussian03的B3LYP/6-31G(d,p)、DMol3的BLYP/DNP和deMon的BLYP/TZVP等方法计算了甲烷水合物(结构-1)中平面五元水分子簇的结合能和氢键能,作了基组重叠误差(BSSE)和色散能(dispersion)的修正,估算了次级相互作用的贡献.在DMol3程序中使用了大型数值基组DNP,将基组重叠误差降至最低.在Gaussi-an03的B3LYP/6-31G(d,p)计算中,采用平衡法(Counterpoise)校正基组重叠误差.两种计算方法给出了一致的结果,证实了在使用6-31G(d,p)基组时,一对水分子在平衡距离的基组重叠误差高达8 kJ/mol.为估算色散能的贡献,使用了新近发展的包含色散能的密度泛函的DFT程序deMon计算了五元水分子簇.用多种方法计算出了经基组重叠误差和色散能修正的五元水分子簇的分子间结合能和氢键能的较为精确的势能超曲面,为甲烷和其他气体水合物的分子动力学模拟提供了依据.     

MP2, DFT and ab initio calculations on thioxanthone

Density functional theory (DFT), HF and MP2 calculations have been carried out to investigate thioxanthone molecule using the standard 6-31+G(d,p) basis set. The results of MP2 calculations show a butterfly structure for thioxanthone. The calculated results show that the predicted geometry can well reproduce the structural parameters. The predicted vibrational frequencies were assigned and compared with experimental IR spectra. A good harmony between theory and experiment is found. The theoretical electronic absorption spectra have been calculated using CIS method. ¹³C and ¹H NMR of the title compound have been calculated by means of B3LYP density functional method with 6-31+G(d,p) basis set. The comparison of the experimental and the theoretical results indicate that density functional B3LYP method is able to provide satisfactory results for predicting NMR properties.     

Electronic effects on O-H proton dissociation energies of phenolic cation radicals: a DFT study

The electronic effects on O-H proton dissociation energies (PDEs) of para- and meta-substituted phenolic cation radicals have been investigated by density functional theory (DFT) using B3LYP function on a 6-31G(d, p) basis set. The calculation results indicate that electron-donating groups raise the O-H PDE and electron-withdrawing groups reduce the parameter, which are opposite to the electronic effects on O-H bond dissociation energies (BDEs). In addition, the electronic effects on O-H PDE are much stronger than those on O-H BDE. The differences result from the distinct electronic effects on stabilities of phenolic cation radicals and parent phenols. The finding also implies the proton-transfer process is unlikely a rate-controlling step for phenolic antioxidants to scavenge free radicals. Moreover, like O-H BDE, O-H PDE correlate better with the resonance parameter R+ than with field/inductive parameter F. Therefore, O-H PDEs of para-substituted phenolic cation radicals are mainly governed by the resonance effect.     

基于噻吩和亚苯基合成的新物质的结构和电子性质

Theoretical study on the geometries and electronic properties of new conjugated compounds based on thiophene and phenylene was carried out. The theoretical ground-state geometries and electronic structures of the studied molecules were obtained using the density functional theory (DFT) method at B3LYP level with 6-31G(d) basis set. The electronic properties were determined by ZINDO/s, CIS/3-21G(d), and TD//B3LYP/3-21G(d) calculations performed on the B3LYP/6-31(d) optimized geometries. The effects of the ring structure and the substituents on the geometries and electronic properties of these materials were discussed. The results of this study indicate how the electronic properties can be tuned by the backbone ring or side group and suggest these compounds as good candidates for opro-electronic applications.     

Structural and Electronic Properties of New Materials Based on Thiophene and Phenylene

Theoretical study on the geometries and electronic properties of new conjugated compounds based on thiophene and phenylene was carried out. The theoretical ground-state geometries and electronic structures of the studied molecules were obtained using the density functional theory (DFT) method at B3LYP level with 6-31G(d) basis set. The electronic properties were determined by ZINDO/s, CIS/3-21G(d), and TD//B3LYP/3-21G(d) calculations performed on the B3LYP/6-31(d) optimized geometries. The effects of the ring structure and the substituents on the geometries and electronic properties of these materials were discussed. The results of this study indicate how the electronic properties can be tuned by the backbone ring or side group and suggest these compounds as good candidates for opto-electronic applications.     

Density functional theory study on the structure and vibrational spectra for cyanuric chloride

The molecular structure and vibrational spectra of cyanuric chloride have been investigated by density functional theory (DFT) using standard B3LYP/6-31G* and B3LYP/6-311+G** method and basis set combinations. The DFT force field transformed to natural internal coordinates was corrected by a well-established set of scale factors that were found to be transferable to the title compound. Both the calculated structural parameters and vibrational frequencies are in good agreement with the available experimental data.     

Magnetic properties of single crystal alpha-benzoin oxime: An EPR study

The electron paramagnetic resonance (EPR) spectra of gamma irradiated single crystals of alpha-benzoinoxime (ABO) have been examined between 120 and 440 K. Considering the dependence on temperature and the orientation of the spectra of single crystals in the magnetic field, we identified two different radicals formed in irradiated ABO single crystals. To theoretically determine the types of radicals, the most stable structure of ABO was obtained by molecular mechanic and B3LYP/6-31G(d,p) calculations. Four possible radicals were modeled and EPR parameters were calculated for the modeled radicals using the B3LYP method and the TZVP basis set. Calculated values of two modeled radicals were in strong agreement with experimental EPR parameters determined from the spectra. Additional simulated spectra of the modeled radicals, where calculated hyperfine coupling constants were used as starting points for simulations, were well matched with experimental spectra.     

Theoretical calculations of hyperfine coupling constants for muoniated butyl radicals

The hyperfine coupling constants (HFCCs) of all the butyl radicals that can be produced by muonium (Mu) addition to butene isomers (1- and 2-butene and isobutene) have been calculated, to compare with the experimental results for the muon and proton HFFCs for these radicals reported in paper II (Fleming, D. G.; et al. J. Phys. Chem. A 2011, 10.1021/jp109676b) that follows. The equilibrium geometries and HFCCs of these muoniated butyl radicals as well as their unsubstituted isotopomers were treated at both the spin-unrestricted MP2/EPR-III and B3LYP/EPR-III levels of theory. Comparisons with calculations carried out for the EPR-II basis set have also been made. All calculations were carried out in vacuo at 0 K only. A C-Mu bond elongation scheme that lengthens the equilibrium C-H bond by a factor of 1.076, on the basis of recent quantum calculations of the muon HFCCs of the ethyl radical, has been exploited to determine the vibrationally corrected muon HFCCs. The sensitivity of the results to small variations around this scale factor was also investigated. The computational methodology employed was "benchmarked" in comparisons with the ethyl radical, both with higher level calculations and with experiment. For the β-HFCCs of interest, compared to B3LYP, the MP2 calculations agree better with higher level theories and with experiment in the case of the eclipsed C-Mu bond and are generally deemed to be more reliable in predicting the equilibrium conformations and muon HFCCs near 0 K, in the absence of environmental effects. In some cases though, the experimental results in paper II demonstrate that environmental effects enhance the muon HFCC in the solid phase, where much better agreement with the experimental muon HFCCs near 0 K is found from B3LYP than from MP2. This seemingly better level of agreement is probably fortuitous, due to error cancellations in the DFT calculations, which appear to mimic these environmental effects. For the staggered proton HFCCs of the butyl radicals exhibiting no environmental effect in paper II, the best agreement with experiment is consistently found from the B3LYP calculations, in agreement also with benchmark calculations carried out for the ethyl radical.     

Theoretical Studies on Electronic Spectrum Property of 2—（2—Hydroxyphenyl）pyridine Via Time—Dependence Density Functional Theory Method

The geometrical structures of 2-(2-hydorxyphenyl) pyridine(PP) and its protonation states were optimized by means of the B3LYP/6-31G(d) method.For all the selected systems,the existence of H-bond is in favor of the stability of the systems.On the basis of the optimized geometrical structures,their electronic spectrum properties were studied by time-dependent density functional theory(TD-DFT)methosd via a hybrid runction of B3LYP and 6-31G(d) basis set.The TD-DFT calculation result predicts the absorption spectrum of PP at 324nm (3.82eV),which is in very good agreement with the experimental value of 322nm (3.85eV)determined in solvent chloroform.The absorption spectra of the two protoation states both exert a red shift in various pH media.     

Proton affinities of ketones, vicinal diketones and α-keto esters: a computational study

The proton affinities of seven different ketones, vicinal diketones, and α-keto esters (acetophenone, 2,2,2-trifluoroacetophenone, 2,3-butanedione, 1-phenyl-1,2-propanedione, methyl pyruvate, ethyl benzoylformate, and ketopantolactone) have been evaluated theoretically using the conventional ab initio HF and several post-HF methods (MP2, MP4, CCSD), density functional methods with the B3LYP hybrid functional, as well as some ab initio model chemistries [CBS-4M, G2(MP2), and G3(MP2)//B3LYP]. The chemical compounds studied are frequently used substrates in the asymmetric hydrogenation over chirally modified platinum catalysts where the protonation properties of the chiral modifier and the substrates are of great interest. In most cases, the proton affinities (PAs) evaluated with the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) methods are in good agreement with the existing experimental ones. However, the previously reported PA of 2,3-butanedione seems to be too high by 10-15 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP model chemistries predict proton affinities that are systematically higher and lower than the experimental PAs, respectively. If proton affinities are evaluated as the average of the PAs calculated with these two theoretical methods a very good agreement with the experimental results is obtained. The mean absolute deviation (MAD) from experiment of this combination method for the PAs of 13 test molecules is 4.0 kJ mol⁻¹. For 9 molecules composed only of first-row atoms the MAD is 2.5 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP methods provide significant savings in computational time and disk space compared to the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) models. Therefore, it is suggested that if no experimental or highly accurate theoretical data is available (due to computational cost), the proton affinities of similar compounds as investigated in this paper, can be evaluated with the combination method. For the studied molecules, this method gives the following PAs (in kJ mol⁻¹): 788 (2,3-butanedione, exptl 802); 798 (2,2,2-trifluoroacetophenone, exptl 799); 811 (ketopantolactone); 813 (methyl pyruvate); 825 (1-phenyl-1,2-propanedione); 862 (acetophenone, exptl 861); 865 (ethyl benzoylformate).     

Spectroscopic, quantum chemical DFT/HF study and synthesis of [2.2.1] hept-2'-en-2'-amino-N-azatricyclo [3.2.1.0(2,4)] octane

The compound 4-N-bicyclo [2.2.1] hept-2'-en-2'-amino-N-azatricyclo [3.2.1.0(2,4)] octane (2) has been synthesized and characterized by elemental analysis, IR, UV-vis, mass and NMR. Density functional theory (DFT) and Hartree-Fock (HF) calculations have been carried out for the title compound by using the standard 6-31G* basis set. The calculated results show that the predicted geometry can well reproduce the structural parameters. Predicted vibrational frequencies have been assigned and compared with experimental IR spectra and they complement each other. The theoretical electronic absorption spectra have been calculated by using CIS, TD-DFT and ZINDO methods. The (13)C NMR and (1)H NMR of compound (2) have been calculated by means of Becke 3-Lee-Yang-Parr (B3LYP) density functional method with 6-31G* basis set. Comparison between the experimental and the theoretical results indicates that density functional B3LYP method is able to provide satisfactory results for predicting NMR properties. On the basis of vibrational analyses, the thermodynamic properties of the title compound at different temperatures have been calculated.     

Effect of substituents on spin density in benzimidazole nitronyl nitroxide radicals studied by electron spin resonance

Several novel benzimidazole‐3‐oxide‐1‐oxyl radicals with substituents at 5 and/or 6 position were synthesized. The ESR analysis of nitrogen hyperfine coupling constants (hfccs) revealed that substituents at 5 and 6‐position affect the spin density to greater extent than substituents on the phenyl ring at 2‐position. Density functional theory calculations of nitrogen hfccs were performed using several different Pople type basis sets, as well as double and triple zeta quality individual gauge for localized orbital (IGLO‐II, IGLO‐III) and electron paramagnetic resonance (EPR‐II, EPR‐II) basis sets. Experimental and theoretical hfccs are compared. Copyright © 2009 John Wiley & Sons, Ltd.     

The O-H stretching mode of a prototypical photoacid as a local dielectric probe

We investigate the OH stretch vibrational frequency shifts of a prototype photoacid, 2-naphthol (2N), when dissolved in solvents of low polarity. We combine femtosecond mid-infrared spectroscopy and a theoretical model based on the Pullin-van der Zwan-Hynes perturbative approach to explore vibrational solvatochromic effects in the ground S(0) and the first electronically excited (1)L(b) states. The model is parametrized using density functional theory (DFT), at the B3LYP/TZVP and TD-B3LYP/TZVP levels for the 2N chromophore in the S(0) and (1)L(b) states, respectively. From the agreement between experiment and theory we conclude that vibrational solvatochromic effects are dominated by the instantaneous dielectric response of the solvent, while time-dependent nuclear rearrangements are of secondary importance.     

Comparison of density functionals for the study of the high spin low spin gap in Fe(III) spin crossover complexes

A detailed investigation of the accuracy of different quantum mechanical methods for the study of iron(III) spin crossover complexes is presented. The energy spin state gap between the high and low spin states; ΔE (HS‐LS) of nine iron(III) quinolylsalicylaldiminate complexes were calculated with nine different DFT functionals, then compared. DFT functionals: B3LYP, B3LYP‐D3, B3LYP*, BH&HLYP, BP86, OLYP, OPBE, M06L, and TPSSh were tested with six basis sets: 3‐21G*, dgdzvp, 6‐31G**, cc‐pVDZ, Def2TZVP, and cc‐pVTZ. The cations from the X‐ray crystal structures of [Fe(qsal‐OMe)₂]Cl·MeCN·H₂O, [Fe(qsal‐OMe)₂]Cl·2MeOH·0.5H₂O, [Fe(qsal‐OMe)₂]BF₄·MeOH, [Fe(qsal‐OMe)₂]NCS·CH₂Cl₂, [Fe(qsal‐F)₂]NCS, [Fe(qsal‐Cl)₂]NCS·MeOH, [Fe(qsal‐Br)₂]NCS·MeOH, [Fe(qsal‐I)₂]OTf·MeOH, and [Fe(qsal)₂]NCS?CH₂Cl₂ were used as starting structures. The results show that B3LYP, B3LYP‐D3, OLYP, and OPBE with a 6‐31G**, Def2TZVP, and cc‐pVTZ basis set give reasonable results of ΔE (HS‐LS) compared with the experimental data. The enthalpy of [Fe(qsal‐I)₂]⁺ calculated with an OLYP functional and cc‐pVTZ basis set (1.48 kcal/mol) most closely matches the experimental data (1.34 kcal/mol). B3LYP* yields an enthalpy of 5.92 kcal/mol suggesting it may be unsuitable for these Fe(III) complexes, mirroring recent results by Kepp (Inorg . Chem ., 2016, 55 , 2717–2727).