Ab initio calculation of inner-sphere reorganization energies of arenediazonium ion couples

The geometries of a series of substituted arenediazonium cations (p-NO2, p-CN, p-Cl, p-F, p-H, m-CH3, p-CH3, p-OH, p-OCH3, p-NH2) and the corresponding diazenyl radicals were optimized at the HF/6-31G, MP2/6-31G, B3LYP/6-31G, B3LYP/TZP, B3PW91/TZP, and CASSCF/6-31G levels of theory. Inner-sphere reorganization energies for the single electron-transfer reaction between the species were computed from the optimized geometries according to the NCG method and compared to experimental values determined by Doyle et al. All levels of theory predicted a CNN bond angle of 180 degrees in the cation. A bent neutral diazenyl radical was predicted at all levels of theory excepting B3LYP/TZP and B3PW91/TZP for the p-Cl-substituted compound. Inner-sphere reorganization energies determined at the HF, MP2, and CASSCF levels of theory correlated poorly with both experimental results and calculated geometries. Density functional methods correlated best with the experimental values, with B3LYP/6-31G yielding the most promising results, although the ROHF/6-31G survey also showed some promise. B3LYP/6-31G calculations correctly predicted the order of the inner-sphere reorganization energies for the series, excluding the halogen-substituted compounds, with values ranging from 42.8 kcal x mol(-1) for the p-NO2-substituted species to 55.1 kcal x mol(-1) for NH2. The magnitudes of these energies were lower than the experimental by a factor of 2. For the specific cases examined, the closed-shell cation geometries showed the expected geometry about the CNN bond, with variations in the CN and NN bond lengths correlating with the electron-donating/withdrawing capacity of the substituent. As predicted by Doyle et al., a large geometry change was observed upon reduction. The neutral diazenyl radicals showed a nominal CNN bond angle of 120 degrees and variations in the CN and NN bond lengths also correlated with the electron-donating/withdrawing capacity of the substituent. Changes in theta(CNN) and r(CN) both correlated well with calculated lambda(inner). The key parameters influencing inner-sphere reorganization energy were the CN and NN bond lengths and the CNN bond angle. This influence is explained qualitatively via resonance models produced from NRT analysis and is related to the amount of CN double bond character. Based on these observations, B3LYP/6-31G calculations are clearly the most amenable for calculating inner-sphere reorganization energies for the single electron-transfer reaction between cation/neutral arenediazonium ion couples.     

Electronic structures of 3d-metal mononitrides

Bond distances, vibrational frequencies, electron affinities, ionization potentials, and dissociation energies of the title molecules in neutral, positively, and negatively charged ions were studied by use of density functional methods B3LYP, BLYP, BHLYP, BPW91, and B3PW91. The calculated results are compared with experiments and previous theoretical studies. It was found that the calculated properties are highly dependent on the functionals employed, in particular for the dissociation energy and vibrational frequency. For neutral species, pure density functional methods BLYP and BPW91 have relatively good performance in reproducing the experimental bond distance and vibrational frequency. For cations, hybrid exchange functional methods B3LYP and B3PW91 are good in predicting the dissociation energy. For both neutral and charged species, BHLYP tends to give smaller dissociation energy.     

Density functional and spin-orbit ab initio study of CF3Br: molecular properties and electronic curve crossing

Quantum chemical calculations of CF(3)Br and the CF(3) radical are performed using density functional theory (DFT) and time-dependent DFT (TDDFT). Molecular structures, vibrational frequencies, dipole moment, bond dissociation energy, and vertical excitation energies of CF(3)Br are calculated and compared with available experimental results. The performance of six hybrid and five hybrid meta functionals in DFT and TDDFT calculations are evaluated. The ωB97X, B3PW91, and M05-2X functionals give very good results for molecular structures, vibrational frequencies, and vertical excitation energies, respectively. The ωB97X functional calculates well the dipole moment of CF(3)Br. B3LYP, one of the most widely used functionals, does not perform well for calculations of the C-Br bond length, bond dissociation energy, and vertical excitation energies. Potential energy curves of the low-lying excited states of CF(3)Br are obtained using the multiconfigurational spin-orbit ab initio method. The crossing point between 2A(1) and 3E states is located near the C-Br bond length of 2.45 ?. Comparison with CH(3)Br shows that fluorination does not alter the location of the crossing point. The relation between the calculated potential energy curves and recent experimental result is briefly discussed.     

An experimental and theoretical study of molecular structure and vibrational spectra of 2-chloronicotinic acid by density functional theory and ab initio Hartree–Fock calculations

In this work, the Fourier transform Raman and Fourier transform infrared spectra of 2-chloronicotinic acid (2-CNA) are recorded in the solid phase. The molecular geometry, vibrational frequencies, infrared intensities and Raman scattering activities of 2-CNA in ground state have been calculated by using ab initio Hartree–Fock (HF) and density functional (B3LYP and B3PW91) methods with 6-31G(d) and 6-311G(d) basis sets level. On the basis of the comparison between calculated and experimental results and the comparison with related molecule, assignments of fundamental vibrational modes are examined. The optimized geometric parameters (bond lengths and bond angles) obtained by using HF show the best agreement with the experimental values of 2-CNA. Comparison of the observed fundamental vibrational frequencies of 2-CNA and calculated results by density functional (B3LYP and B3PW91) and Hartree–Fock methods indicates that B3LYP is superior to the scaled Hartree–Fock and B3PW91 approach for molecular vibrational problems.     

Electronic structures of 4d transition metal monoxides by density functional theory

Bond distances, dissociation energies, ionization potentials and electron affinities of 4d transition metal monoxides from YO to CdO and their positive and negative ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, SVWN, MPW1PW91 and PBE1PBE. It was found that calculated properties are highly dependent on the functionals employed, especially for dissociation energy. For most neutral species, pure density functionals BLYP, BPW91 and BP86 have good performance in predicting dissociation energy than hybrid density functionals B3LYP, B3PW91 and B3P86. In addition, BLYP gives the largest bond distance compared with other density functional methods, while SVWN gives shortest bond distance, largest dissociation energy and electron affinity. For the ground state, the spin multiplicity of the charged species can be obtained by ± 1 of their corresponding neutral species.     

The performance of different density functional methods in the calculation of molecular structures and vibrational spectra of platinum(II) antitumor drugs: cisplatin and carboplatin

A comparison of eight density functional models for predicting the molecular structures, vibrational frequencies, infrared intensities, and Raman scattering activities of platinum(II) antitumor drugs, cisplatin and carboplatin, is reported. Methods examined include the pure density functional protocols (G96LYP, G96PW91, modified mPWPW and original PW91PW91), one‐parameter hybrid approaches (mPW1PW and mPW1LYP), and three‐parameter hybrid models (B3LYP and B3PW91), as well as the HF and MP2 levels of theory. Different effective core potentials (ECPs) and several basis sets are considered. The theoretical results are discussed and compared with the experimental data. It is remarkable that the mPW1PW protocol introduced by Adamo and Barone [J Chem Phys 1998, 108, 664], is clearly superior to all the remaining density functional methods (including B3LYP). The geometry and vibrational frequencies of cisplatin and carboplatin calculated with the mPW1PW method, and the ECP of Hay and Wadt (LanL2DZ basis set) are in better agreement with experiment than those obtained with the MP2 method. The use of more elaborated ECP and the enlargements of basis sets do not significantly improve the results. A clear‐cut assignments of the platinum‐ligand vibrations in cisplatin and carboplatin are presented. It is concluded that mPW1PW is the new reliable method, which can be used in predicting molecular structures and vibrational spectra of large coordination compounds containing platinum(II). © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 901–912, 2001     

Density functional theory studies on the electronic and vibrational spectra of octaethylporphyrin diacid

The ground-state structure and electronic and vibrational spectra of octaethylporphyrin diacid (H4OEP2+) have been studied with the density functional theory. The geometrical parameters computed with B3LYP, PBE1PBE and mPW1PW91 functionals and 6-31G* basis sets are well consistent with the experimental values. Electronic absorption spectrum of H4OEP2+ has been studied with the time-dependent DFT method, and the calculated excitation energies and oscillator strengths are compared with the experimental results. The Raman and IR spectra of H4OEP2+ and the Raman spectrum of its N-deuterated analogue (D4OEP2+) were measured. The observed Raman and IR bands have been assigned based on the frequency calculations at the B3LYP/6-31G* level of theory.     

Density functional study of AuXq (X = O, S, Se, Te, q = +1, 0, -1) molecules

Bond distances, vibrational frequencies, electron affinity, ionization potential, and dissociation energies of the title molecules were studied by use of density functional methods B3LYP, B3P86, B3PW91, BHLYP, BLYP, BP86, mPW1PW91, and PBE1PBE. It was found that the ground electronic state is doublet for neutral species, singlet for the anion, and triplet for the cation, in agreement with experiments and previous theoretical studies. The calculated properties are highly dependent on the functionals employed, in particular for the dissociation energy. The predicted bond distances and vibrational frequencies are in agreement with experiments and previous theoretical results. BP86 and BLYP have relatively good performance in reproducing the experimental results, while BHLYP is the worst functional method compared with the other density functional methods used for the title molecules.     

Comparison between optimized geometries and vibrational frequencies calculated by the DFT methods for the Anderson-type heteropolyanion: hexamolybdoaluminate(III), [AlIII(OH)6Mo6O18]3-

Optimized geometries and vibrational frequencies were calculated for the hexamolybdoaluminate(III), [AlIII(OH)6Mo6O18]3-, Anderson-type heteropolyanion with the HF, B3LYP, B3PW91, B3P86 and B1LYP methods of theory using the LanL2DZ, SDD and combination of LanL2DZ with 6-31G (d, p) basis sets. The agreement between the optimized and experimental geometries was in the decreasing order: HF, B3P86, B3PW91, B1LYP and B3LYP. The calculated frequencies by the B3LYP have the smallest mean root mean square (RMS) error. The effect of the basis set on the calculated bond lengths and frequencies by the density functional calculations (DFT) methods was minor. The agreement between the previously reported IR and Raman spectra and the calculated values is, in general, good.     

Electronic Structures of 5d Transition Metal Monoxides by Density Functional Theory

Bond distances, vibrational frequencies, electron affinities, ionization potentials, and dissociation energies of the diatomic 5d transition metal (except La) monoxides and their positively and negatively charged ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, MPW1PW91, PBE1PBE, and SVWN. Our calculation shows that for each individual species, the calculated properties are quite sensitive to the method used. Compared with hybrid density functional method B3PW91 (B3P86), pure density functional method BPW91 (BP86) gives longer bond distance (lower vibrational frequency) from HfO to PtO for neutral species, HfO⁺ to IrO⁺ for cationic species, and HfO⁻ to AuO⁻ for anionic species. While for B3LYP and BLYP, the trend was observed for cationic species from HfO⁺ to IrO⁺ and anionic species from HfO⁻ to AuO⁻ (except TaO⁻), but not for neutrals. Pure density function methods BLYP, BPW91, and BP86 give larger dissociation energy compared with hybrid density functional methods B3LYP, B3PW91, and B3P86. SVWN in most cases gives the smallest bond distance, while BLYP gives the largest value. MPW1PW91 and PBE1PBE show the same performance in predicting the spectroscopic constants. In addition, useful empirical criteria that one has obtained the ground states of a species and its ions are the spin multiplicities of a neutral and its single charged ions which differs by ±1.     

Thermochemistry of the higher chlorine oxides ClOx (x=3, 4) and Cl2Ox (x=3–7)†

Heats of formation for ClO₃, ClO₄, Cl₂O₃, Cl₂O₄, Cl₂O₅, Cl₂O₆ and Cl₂O₇ molecules are determined at the B3LYP, B3PW91, mPW1PW91 and B1LYP levels of the density functional theory employing a series of extended basis sets, and using Gaussian-3 model chemistries. Modified Gaussian-3 calculations, which employ accurate B3LYP/6-311+G(3d2f) molecular geometries and vibrational frequencies, were also performed. Heats of formation were calculated from both total atomization energies and isodesmic reaction schemes. The latter method in conjunction with Gaussian-3 models leads to the most reliable results. The best values at 298 K for ClO₃, ClO₄, Cl₂O₃ and Cl₂O₄ as derived from an average of G3//B3LYP and G3//B3LYP/6-311+G(3d2f) calculations are 43.1, 54.8, 31.7 and 37.4 kcal mol⁻¹. From calculations carried out at the G3(MP2)//B3LYP and G3(MP2)//B3LYP/6-311+G(3d2f) levels, heats of formation for Cl₂O₅, Cl₂O₆ and Cl₂O₇ are predicted to be 53.2, 52.2 and 61.5 kcal mol⁻¹. All best values are reproduced within 1 kcal mol⁻¹ by using mPW1PW91/6-311+G(3d2f) isodesmic energies. Enthalpy changes for relevant Cl–O bond fission reactions are reported. Comparisons with previous thermodynamics data are made.     

Computational studies of bond dissociation energies, ionization potentials, and heat of formation for NH and NH+. Are hybrid density functional theory methods as accurate as quadratic complete basis set and Gaussian-2 ab initio methods?

Ionization potentials, bond dissociation energies, and heat of formation for NH and NH⁺ molecular species as well as for their elements were computed with highly reliable quadratic complete basis set and Gaussian-2 ab initio methods. The results are compared with experimental results and the assurance of these ab initio approaches is assessed. The same studies were also performed with three hybrid density functional methods (B3LYP, B3P86, and B3PW91) in combination with variously sized basis sets. The computational results are discussed in light of density functional theory reliability for exploring the potential energy of small polar molecular systems. Received: 21 July 1997 / Accepted: 8 December 1997     

Theoretical study on infrared vibrational spectra of boric-acid in gas-phase using density functional methods

Density functional theory (DFT) methods with various exchange-correlation functionals such as SVWN, BVWN, BVWN5, BLYP, B1LYP, B3LYP, B3PW91, and BH and H are employed in a theoretical study of molecular boric-acid in gas-phase. In the calculations, the split valence 6-311++G** and 6-31G* basis sets were used. The geometry, zero-point vibrational energies (ZPVEs), and harmonic infrared vibrational (IR) frequencies are predicted. The calculated C_3h-symmetry geometrical parameters are compared with Hartree–Fock (HF) calculation results and experimental data. IR frequencies predicted by the BLYP, B3LYP, and B3PW91 calculations are in good agreement with experimental data. The frequency calculations presented here also suggest that the C_3h-symmetrical structure corresponds to a minimum in the potential energy surface, but neither is D_3h- or C₃-symmetrical structure.     

Electronic structures and spectroscopic properties of rhenium (I) tricarbonyl photosensitizer: [Re(4,4'-(COOEt)2-2,2'-bpy)(CO)3py]PF6

The ground state and lowest triplet-state structures of [Re(4,4'-(COOEt)(2)-2,2'-bpy)(CO)(3)py]PF(6) photosensitizer (bpy=bipyridine, py=pyridine) have been studied with density functional theory (DFT). Time-dependent density functional theory (TD-DFT) was carried out to predict the photophysical properties of the photosensitizer. The effects of the solvents were evaluated using the conductor-like polarizable continuum (CPCM) method in dichloromethane, chloroform, acetonitrile, acetone, ethanol and dimethylsulfoxide. The electronic transition energies computed with BLYP, MPWPW91, B3LYP and MPW1PW91 functionals are compared with the experimental spectra. Based on the calculated excited energies, the experimental absorption maximum is assigned as metal-to-ligand charge transfer (MLCT) and ligand-to-ligand charge transfer (LLCT) mixed transition, and the luminescence originates from the lowest triplet state that is ascribed as the mixed transition of MLCT/LLCT.     

Comparison of DFT methods for molecular structure and vibration spectra of ofloxacin calculations

Comparison of the performance of different density functional theory (DFT) methods at various basis sets in predicting molecular and vibration spectra of ofloxacin was reported. The methods employed in this study comprise six functionals, namely, mPW1PW91, HCTH, LSDA, PBEPBE, B3PW91 and B3LYP. Different basis sets including LANL2DZ, SDD, LANL2MB, 6-31g, 6-311g and 3-21g were also examined. Comparison between the calculated and experimental data indicates that the mPW1PW91/6-311g level afford the best quality to predict the structure of ofloxacin. The results also indicate that B3LYP/LANL2DZ level show better performance in the vibration spectra prediction of ofloxacin than other DFT methods.     

Comparative Studies of 1,4-Bis[2-（4-Pyridyl）ethenyl]-benzene Using Density Functional Theory

The optimized molecular structure and harmonic vibrational frequencies of a 1,4-bis[2-(4-pyridyl)ethenyl]-benzene(BPENB)molecule were calculated via five popular density functional theory(DFT)methods.On the basis of the comparison between calculated and experimental results,it is concluded that the B3PW91 and B3LYP methods are superior to the others in optimizing structures,and the BPW91 method reproduces the observed fundamental frequencies most satisfactorily.     

Performance of density functionals for calculation of reductive ring-opening reaction energies of Li<Superscript>+</Superscript>-EC and Li<Superscript>+</Superscript>-VC

Reaction energies were determined for reductive ring-opening reactions of Li⁺-coordinated ethylene carbonate (EC) and vinylene carbonate (VC) by using various density functional theory (DFT) and ab-initio methods applying the basis sets up to Dunnings aug-cc-pVQZ. The methods examined include the local density functional (SVWN), the pure gradient-corrected density functionals (BLYP and BPW91), and the hybrid density functionals (B3LYP, B1LYP, B3PW91, and mPW1PW91). Comparison of the DFT results with ab-initio results indicates that the mPW1PW91 approach introduced by Adamo and Barone, is superior to all the other DFT methods (including B3LYP). The performance of more cost-effective Pople-type basis sets ranging from 6-31G(d,p) to 6-311++G(3df,3pd) was assessed at DFT levels of theory by calibrating them with the aug-cc-pVQZ results