DFT calculations on the catalytic oxidation of CO over Si-doped (6,0) boron nitride nanotubes

The oxidation of carbon monoxide (CO) is important for a series of technological and environmental applications. In this work, the catalytic oxidation of CO on Si-doped (6,0) boron nitride nanotubes (BNNTs) is investigated by using density functional theory calculations. Reaction barriers and corresponding thermodynamic parameters were calculated using the M06-2X, B3LYP and wB97XD density functionals with 6-31G* basis set. Our results indicate that a vacancy defect in BNNT strongly stabilizes the Si adatom and makes it more positively charged. This charging enhances the adsorption of reaction gases (O₂ and CO) and results in the change of the electronic structure properties of the tube. The calculated barrier of the reaction CO + O₂ → CO₂ + O_ads on Si-doped BNNTs following the Langmuir–Hinshelwood is lower than that on the traditional noble metal catalysts. The second step of the oxidation would be the Eley–Rideal reaction (CO + O_ads → CO₂) with an energy barrier of about 1.8 and 10.1 kcal/mol at M06-2X/6-31G* level. This suggests that the CO oxidation catalyzed by the Si-doped BNNTs is likely to occur at the room temperature. The results also demonstrate that the activation energies and thermodynamic quantities calculated by M06-2X, B3LYP and wB97XD functionals are consistent with each other.     

Basis set convergence of indirect spin-spin coupling constants in the Kohn-Sham limit for several small molecules

The performance of more than 40 density functionals in predicting indirect spin-spin coupling constants (SSCCs) in the Kohn-Sham basis set limit was tested. For comparison, similar calculations were performed using the RHF, SOPPA, SOPPA(CC2), and SOPPA(CCSD) methods, and the results were estimated toward the complete basis set (CBS) limit. The SSCCs of nine small molecules (N(2), CO, CO(2), NH(3), CH(4), C(2)H(2), C(2)H(4), C(2)H(6), and C(6)H(6)) were calculated using the dedicated Jensen pcJ-n polarization-consistent basis sets and used for the CBS limit estimations within the Kohn-Sham limit. These CBS results were compared with calculations using the aug-cc-pVTZ-J basis set. Among the 41 studied DFT methods, the tHCTHhyb, HSEh1PBE, HSE2PBE, wB97XD, wB97, and wB97X functionals reproduced accurately the experimental (1)J(XH) SSCCs and (3)J(HH60) and (2)J(HH(gem)) in ethane. Similarly, the functionals HSEh1PBE, HSE2PBE, wB97XD, wB97, and wB97X predicted accurately (1)J(CC), and B98, B97-1, B97-2, PBE1PBE, B1LYP, and O3LYP provided accurate (1)J(CO) results in the CO molecule. A very good performance for the calculation of the SSCCs based on the use of the relatively small basis set aug-cc-pVTZ-J was observed.     

Photoinduced Charge Transfer Processes of Zinc Porphyrin Derivatives for Dye-sensitized Solar Cells

We studied the excitation energies of zinc porphyrin(ZnP) and zinc porphyrin analogues(ZnP-R) with long-range corrected(LC) density functional theory. We compared three recent LC functionals, wB97XD, CAM-B3LYP, LC-wPBE and functionals B3LYP with the experiments. It was found that the low energy excitation is well predicted for ZnP with the LC-wPBE functional based on 6-31G(d) basis set. Excited-state geometry optimizations for all the compounds were carried out. It was shown that upon the photo-excitation nearly 1 e transferred to the acceptor unit[dodecafluorosubphthalocyanine, SubPc(F)₁₂] and a highly polarized state formed in compound 5. But after the relaxation of the first excited state, the electron on the acceptor went back to the ZnP a little and then shifted from ZnP to the donor segment[bis(4-hexylphenyl)amino] partially. The computed excited-state radiative lifetime(τ) for compound 5 is 943 μs, which shows a reasonable agreement with the experimental observation. According to the long-lived exciton in strong push-pull compound 5, we proposed that electron injection to semiconductor might occur after vibrational relaxation of excited state.     

Benchmarking of density functionals for the kinetics and thermodynamics of the hydrolysis of glycosidic bonds catalyzed by glycosidases

In recent years, there has been an increased interest in understanding the enzymatic mechanism of glycosidases resorting mostly to DFT and DFT/MM calculations. However, the performance of density functionals (DFs) is well known to be system and property dependent. Trends drawn from general studies, despite important to evaluate the quality of the DFs and to pave the way for the development of new DFs, may be misleading when applied to a single specific system/property. To overcome this issue, we carried out a benchmarking study of 40 DFs applied to the geometry optimization and to the electronic barrier height (E _Barrier) and electronic energy of reaction (E _R) of prototypical glycosidase‐catalyzed reactions. Additionally, we report calculations with SCC‐DFTB and four semiempirical MO methods applied to the same problem. We have used a universal molecular model for retaining glycosidases, comprising only a 22‐atoms system that mimics the active site and substrate. High accuracy reference geometries and energies were calculated at the CCSD(T)/CBS//MP2/aug‐cc‐pVTZ level of theory. Most DFs reproduce the reference geometries extremely well, with mean unsigned errors (MUE) smaller than 0.05 Å for bond lengths and 3° for bond angles. Among the DFs, wB97X‐D, CAM‐B3LYP, B3P86, and PBE1PBE have the best performance in geometry optimizations (MUE = 0.02 Å). Conversely, semiempirical MO and SCC‐DFTB methods yielded less accurate geometries (MUE between 0.09 and 0.17 Å). The inclusion of D3 correction has a small, but still relevant, influence in the geometry predicted by some DFs. Regarding E _Barrier, 11 DFs (MPW1B95, CAM‐B3LYP, M06 ‐ 2X, PBE1PBE, wB97X ‐ D, B1B95, BMK, MN12 – SX, M05, M06, and M11) presented errors below 1 kcal.mol^?1, in relation to the reference energy. Most of these functionals belong to the family of hybrid functionals (H‐GGA, HH‐GGA, and HM‐GGA), which shows a positive influence of HF exchange in the determination of E _Barrier. The inclusion of D3 correction has not affected significantly the E _Barrier and E _R. The use of geometries at the accurate but expensive MP2/aug‐cc‐pVTZ level of theory has a small, albeit not insignificant, influence in the E _Barrier when compared with energies calculated with geometries determined with the DFs (usually a few tenths of kcal.mol^?1, with exceptions). In general, semiempirical MO methods and DFTB are associated with larger errors in the determination of E _Barrier, with unsigned errors from 6.9 to 24.7 kcal.mol^?1.     

A mixed experimental and theoretical study on chelidamate copper(II) complex with 4-methylpyrimidine

A new chelidamate complex, [Cu(chel)(H₂O)₂(mpd)] (chel = chelidamate; mpd = 4-methylpyrimidine), has been synthesized and characterized through a combination of single crystal X-ray analysis, electron paramagnetic resonance (EPR), ultraviolet-visible (UV-vis), and fourier transform infrared spectroscopy (FT-IR). The complex has six-coordinate distorted octahedral geometry around Cu(II). The theoretical vibrational frequencies and optimized geometric parameters (bond lengths and angles) have been calculated using Density Functional Theory (DFT)/B3LYP and Hartree Fock quantum chemical methods with 6-31G(d, p) basis set by Gaussian 09W software. The EPR spectrum of the compound showed that the paramagnetic center has rhombic symmetry. The EPR studies were carried out using the following unrestricted hybrid density functionals: B3LYP, CAM-B3LYP, HSEH1PBE, WB97XD, MPW1PW91, and BPV86. The UV–vis absorption spectra have been examined in different media and compared with the calculated one using TD-DFT method by applying the polarizable continuum model. Natural bond orbital property of complex has been performed by DFT/B3LYP with 6-31G (d, p) basis set.     

Model calculations for the base-pairing specificity of mutagenic exocyclic DNA adduct 1,N 6-ethenoadenine

1,N ⁶-ethenoadenine (εA) has been studied more extensively due to highly fluorescent and mutagenic nature. We report the model calculations for the base-pairing specificity of mutagenic 1,N ⁶-ethenoadenine adduct. To shed insight into mutagenic process of DNA damage based on geometrical characteristics, electronic properties, B3LYP, M06, B97D, and wB97XD methods have been employed for these model calculations. Single point energy calculation at MP2/6-311++G** levels on corresponding B3LYP, M06, B97D, wB97XD-optimized geometries have also been carried out to better estimate the hydrogen-bonding strengths. The polarisable conductor calculation model (CPCM-UAKS) that accounts for the overall polarizability of the solvent has also been employed. The computed reaction enthalpy values concluded that both guanine and thymine are most favorable candidates to be misincorporated to 1,N ⁶-ethenoadenine adduct, which also in good agreement to experimental report (Leonard, Biochemistry 33: 4755–4761, 79).     

Ferric complexes of 3-hydroxy-4-pyridinones characterized by density functional theory and Raman and UV-vis spectroscopies

Deferiprone and other 3-hydroxy-4-pyridinones are used in metal chelation therapy of iron overload. To investigate the structure and stability of these compounds in the natural aqueous environment, ferric complexes of deferiprone and amino acid maltol conjugates were synthesized and studied by computational and optical spectroscopic methods. The complexation caused characteristic intensity changes, a 300× overall enhancement of the Raman spectrum, and minor changes in UV-vis absorption. The spectra were interpreted on the basis of density functional theory (DFT) calculations. The CAM-B3LYP and ωB97XD functionals with CPCM solvent model were found to be the most suitable for simulations of the UV-vis spectra, whereas B3LYP, B3LYPD, B3PW91, M05-2X, M06, LC-BLYP, ωB97XD, and CAM-B3LYP functionals were all useful for simulation of the Raman scattering. Characteristic Raman band frequencies for 3-hydroxy-4-pyridinones were assigned to molecular vibrations. The computed conformer energies consistently suggest the presence of another isomer of the deferiprone-ferric complex in solution, in addition to that found previously by X-ray crystallography. However, the UV-vis and Raman spectra of the two species are similar and could not be resolved. In comparison to UV-vis, the Raman spectra and their combination with calculations appear more promising for future studies of iron sequestrating drugs and artificial metalloproteins as they are more sensitive to structural details.     

Determination of Tryptophan using a p-Sulfonated Calix [4,6,8]arene Modified Gold Electrode

The ability of p-sulfonated calix[n]arenes (n = 4, 6, 8) to form complexes with tryptophan was studied. The electrochemical properties of these complexes immobilized on gold surfaces were examined by cyclic voltammetry. Parameters affecting the performance of the modified electrodes including the arene concentration, scan rate, applied potential, and pH were optimized. Under the optimal conditions, the method had a linear response to tryptophan between 1 × 10^?7 M and 1 × 10^?5 M with a detection limit of 3 × 10^?8 M. The interaction of the tryptophan–p-sulfonated calix[4]arene complex was more stable than the tryptophan–p-sulfonated calix[6]arene and p-sulfonated calix[8]arene complexes. Molecular modeling calculations indicated that electrostatic interactions and structural matching effects were predominant stabilizing factors. The modified electrodes demonstrated good reproducibility and high selectivity, illustrating their effectiveness for analytical measurements.     

Divergent Pathways and Competitive Mechanisms of Metathesis Reactions between 3‐Arylprop‐2‐ynyl Esters and Aldehydes: An Experimental and Theoretical Study

Mechanistic studies of the reaction between 3‐arylprop‐2‐ynyl esters and aldehydes catalyzed by BF₃ ? Et₂O were performed by isotopic labeling experiments and quantum chemical calculations. The reactions are shown to proceed by either a classical alkyne–carbonyl metathesis route or an unprecedented addition–rearrangement cascade. Depending on the structure of the starting materials and the reaction conditions, the products of these reactions can be Morita–Baylis–Hillman (MBH) adducts that are unavailable by traditional MBH reactions or E‐ and Z‐α,β‐unsaturated ketones. ¹⁸O‐Labeling studies suggested the existence of two different reaction pathways to the products. These pathways were further examined by quantum chemical calculations that employed the DFT(wB97XD)/6‐311+G(2d,p) method, together with the conductor‐like screening model for realistic solvation (COSMO‐RS). By using the wB97XD functional, the accuracy of the computed data is estimated to be 1–2 kcal mol^?1, shown by the careful benchmarking of various DFT functionals against coupled cluster calculations at the CCSD(T)/aug‐cc‐pVTZ level of theory. Indeed, most of the experimental data were reproduced and explained by theory and it was convincingly shown that the branching point between the two distinct mechanisms is the formation of the first intermediate on the reaction pathway: either the four‐membered oxete or the six‐membered zwitterion. The deep mechanistic understanding of these reactions opens new synthetic avenues to chemically and biologically important α,β‐unsaturated ketones.     

Assessment of density functionals and force field methods on anion–π interaction in heterocyclic calix complexes

The performance of ten density functionals and four force field methods in describing non-covalent interactions have been assessed by studying the interaction energies and structures of the four anion–π complexes involving tetraoxacalix[2]arene[2]triazine and various anions. Their structures are optimized at MP2/6-311++G(d,p) level, and interaction energies are obtained at DF-MP2-F12/aug-cc-pVDZ level. The result shows that the functional M06-2X predicts the most reliable interaction energy, followed by wB97XD and BHandH. B97D slightly overestimates the interaction energy. Other functionals and force field methods seriously overestimate the interaction energy. For the structures, three functionals M06-2X, wB97XD and BH and H predict the most reliable results, followed by B97D. The force field methods predict the largest deviations. The present work suggests that the functional M06-2X is a reliable method to describe energies and structural properties of the large molecules involving the anion–π interactions.     

Theoretical Modelling of Photoswitching of Hyperpolarisabilities in Ruthenium Complexes

Static excited‐state polarisabilities and hyperpolarisabilities of three Ru^II ammine complexes are computed at the density functional theory (DFT) and several correlated ab initio levels. Most accurate modelling of the low energy electronic absorption spectrum is obtained with the hybrid functionals B3LYP, B3P86 or M06 for the complex [Ru^II(NH₃)₅(MeQ⁺)]³⁺ (MeQ⁺=N‐methyl‐4,4′‐bipyridinium, 3 ) in acetonitrile. The match with experimental data is less good for [Ru^II(NH₃)₅L]³⁺ (L=N‐methylpyrazinium, 2 ; N‐methyl‐4‐{E,E‐4‐(4‐pyridyl)buta‐1,3‐dienyl}pyridinium, 4 ). These calculations confirm that the first dipole‐ allowed excited state (FDAES) has metal‐to‐ligand charge‐transfer (MLCT) character. Both the solution and gas‐phase results obtained for 3 by using B3LYP, B3P86 or M06 are very similar to those from restricted active‐space SCF second‐order perturbation theory (RASPT2) with a very large basis set and large active space. However, the time‐dependent DFT λ_max predictions from the long‐range corrected functionals CAM‐B3LYP, LC‐ωPBE and wB97XB and also the fully ab initio resolution of identity approximate coupled‐cluster method (gas‐phase only) are less accurate for all three complexes. The ground state (GS) two‐state approximation first hyperpolarisability β_2SA for 3 from RASPT2 is very close to that derived experimentally via hyper‐Rayleigh scattering, whereas the corresponding DFT‐based values are considerably larger. The β responses calculated by using B3LYP, B3P86 or M06 increase markedly as the π‐conjugation extends on moving along the series 2 → 4 , for both the GS and FDAES species. All three functionals predict substantial FDAES β enhancements for each complex, increasing with the π‐conjugation, up to about sevenfold for 4 . Also, the computed second hyperpolarisabilities γ generally increase in the FDAES, but the results vary between the different functionals.     

The accuracy of DFT-optimized geometries of functional transition metal compounds: a validation study of catalysts for olefin metathesis and other reactions in the homogeneous phase

We have investigated the performance of eight popular density functionals, four of which are "standard" functionals not including dispersion (B3LYP, BP86, PBE, and TPSS) and four of which have been constructed to account for dispersion (B97D, wB97XD, M06, and M06L), in reproducing 18 molecular structures derived from single-crystal X-ray diffraction experiments on ruthenium-based olefin metathesis catalyst precursors. Our analysis of all the internuclear distances reveals that standard DFT predicts systematically expanded structures. In contrast, all the methods accounting for dispersion give rise to more compact structures, removing the systematic overestimation of internuclear distances. The contracting effect of dispersion is general and also affects chemical bonds, thus reducing the general overestimation of bond lengths. The best overall performance is observed for wB97XD, which offers relatively small statistical errors when considering the overall structure as well as selected distances. Only for the coordination center geometry is the accuracy of wB97XD matched by standard functionals such as PBE and TPSS, whereas M06 and M06L are associated with larger errors. At the other end of the scale, B3LYP is seen to give the largest statistical errors in general, both when considering the complete structures and the geometries of the coordination centers alone. For the organic ligands, however, B3LYP performs clearly better than the other standard functionals although not as well as the functionals accounting for dispersion. Extending the basis sets is seen to improve the structures in particular of the coordination center, thus underlining the importance of using sufficiently flexible basis sets if highly accurate geometries are to be obtained. Similar conclusions to those obtained for the ruthenium catalysts can be drawn from comparisons of the X-ray crystal structures of 10 other organometallic complexes of relevance to homogeneous catalysis, covering first (Ti, Fe, Co, Ni), second (Zr, Mo, Rh, Pd) and third (W, Ir) row transition metals, with those of DFT. The latter analyses thus offer a first indication that the picture obtained for the ruthenium alkylidene complexes may be extended to other classes of relatively large transition metal complexes.     

The relative energies of polypeptide conformers predicted by linear scaling second-order Møller-Plesset perturbation theory

We describe an implementation of the cluster-in-molecule(CIM) resolution of the identity(RI) approximation second-order M?ller–Plesset perturbation theory(CIM-RI-MP2), with the purpose of extending RI-MP2 calculations to very large systems. For typical conformers of several large polypeptides, we calculated their conformational energy differences with the CIM-RI-MP2 and the generalized energy-based fragmentation MP2(GEBF-MP2) methods, and compared these results with the density functional theory(DFT) results obtained with several popular functionals. Our calculations show that the conformational energy differences obtained with CIM-RI-MP2 and GEBF-MP2 are very close to each other. In comparison with the GEBF-MP2 and CIM-RI-MP2 relative energies, we found that the DFT functionals(CAM-B3LYP-D3, LC-?PBE-D3, M05-2X, M06-2X and ?B97XD) can give quite accurate conformational energy differences for structurally similar conformers, but provide less-accurate results for structurally very different conformers.     

Quantum chemical investigations on hydrogen bonding interactions established in the inclusion complex β-cyclodextrin/benzocaine through the DFT,AIM and NBO approaches

Structure and stability of an inclusion complex formed by Benzocaine (BZC) and β-cyclodextrin (β-CD) were investigated computationally using different levels of theory. The conformational research based on PM6 method allowed reach two minimum-energy structures: model A and model B. The lowest conformers have been exposed to fully geometry optimization employing four DFT functionals: B3LYP, CAM-B3LYP, M05-2X and M06-2X. The performed DFT calculations have identified the model B, in which the amino group is located at the primary face of β-CD, as the most stable complex by an amount up to ?40 kcal/mol. Further, the greater stabilization of model B in respect to model A, has been ascertained through AIM and NBO analyses which clarified the main hydrogen bonds HBs interactions governing the reactivity of BZC inside the hydrophobic cavity of β-CD. Finally, the estimated isotropic ¹H nuclear magnetic shielding constants generated from the gauge-including-atomic-orbital calculation have been analyzed and then compared with the available experimental data.     

Density Functional Study of the Reaction Mechanism of Two Oxiimine Alcohol Formations and Their Novel Rearrangements

In this study, the reaction mechanisms of isonitrosoacetophenone (inapH) with ethanolamine (ea) and 1‐phenylethanolamine (pea) have been investigated theoretically using B3LYP/6‐311G(d,p) method to explain why the formation and unexpected rearrangement products occur or not occur. While the reaction between isonitrosoacetophenone (inapH) with ethanolamine gives oximine alcohol ( Ib ), the reaction of 1‐phenylethanolamine with inapH results in the formation of oximine alcohol with a different substituent ( Ia ) and amido alcohol ( IIa ), which is the unexpected rearrangement product. The rearrangement driving forces of compounds from Ia to IIa are calculated as ca. 28 and 23 kJ/mol in the gas and EtOH phases, respectively. These driving forces have been calculated ca. 46 and 45 kJ/mol for the rearrangement of compound Ib to obtain IIb in the same phases, respectively. This high driving force shows that the compound IIb cannot be obtained from rearrangement of compound Ib as described experimentally in the literature. In addition, as the DFT functionals poorly describe dispersion effects, dispersion correction for reaction heat and free‐energy barrier was estimated using the wB97X‐D/6‐311G(d,p). In general, the relative free energies of all molecules calculated from wB97XD method are lower than performed from B3LYP level. The changes of thermodynamic properties for all molecules with temperature ranging from 100 to 500 K have been obtained using the statistical thermodynamic method.     

Theoretical study on nonlinear optical properties of the Li+[calix[4]pyrrole]Li−dimer,trimer and its polymer with diffuse excess electrons

The static (hyper)polarizabilities of the dimer and trimer with diffuse excess electrons, [Li⁺[calix[4]pyrrole]Li^?]_n, are firstly investigated by the DFT(B3LYP) method in detail. For the dimer and trimer, a Li atom inside each calix[4]pyrrole unit is ionized to form a diffuse excess electron. The results show that the dimer and trimer containing two and three excess electrons, respectively, have very large first hyperpolarizablities as 2.3 × 10⁴ and 4.0 × 10⁴ au, which are 30 and 40 times larger than that of the corresponding [calix[4]pyrrole]_n (n = 2, 3) without Li atom. Also, β values of dimer and trimer are twice and four times as large as that of monomer containing one excess electron. Obviously, not only excess electron but also the number of excess electron plays an important role in increasing the first hyperpolarizability. Moreover, the (hyper)polarizabilities of the [Li⁺[calix[4]pyrrole]Li^?]_n polymer are investigated at ab initio level by using the elongation finite‐field (elongation FF) method. All the oligomers of the [Li⁺[calix[4]pyrrole]Li^?]_n with many excess electrons exhibit very large first hyperpolarizability and large second hyperpolarizability. The present investigation shows that by introducing several and more excess electrons into the nonlinear optical (NLO) materials will be an important strategy for improving their NLO properties, which will be helpful for design of NLO materials. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Coupled cluster,MP2, and DFT study of structures,stabilities, vibrations,and bonding properties of XXeOH (X = F,Cl, Br,and I)

The optimized geometries, molecular properties, and stabilities of new noble gas molecules, XXeOH (X = F, Cl, Br, and I), were studied using CCSD, MP2, CAM-B3LYP, and WB97XD methods and large basis sets. All XXeOH molecules showed equilibrium structures with Cs symmetry. The results also showed that some bonds in XXeOH could be presented as a typical ionic bond. An alteration in ion-pair character was observed for IXeOH, showing two OH ^? and IXe ⁺ parts, while in other molecules, they could be presented as XeOH ⁺ and X ^? . Two decomposition routes were proposed for these molecules that showed high exothermic reactions. However, despite their low thermodynamic stabilities, their decomposition rate constants were small and all molecules (except BrXeOH) had high kinetic stabilities, indicating the possibility for identification and characterization of these molecules. However, in addition to the calculation of their vibrational frequencies, NBO atomic charges, and hybridizations, the bonding properties of XXeOH molecules were studied by AIM calculations (to calculate electron densities, bond elipticities, and Laplacian of electron densities) and second-order intramolecular perturbation energies using NBO calculations. Moreover, the ease of formations and relative stabilities of XXeOH molecules were compared using heats of formations, Gibbs free energies of formations and isodesmic reactions. These calculations showed that the stability of XXeOH molecules was decreased from F to I.     

Selective complexation of alkali metal ions and nanotubular cyclopeptides: a DFT study

A density functional theory based on interaction of alkali metal cations (Li⁺, Na⁺, K⁺, Rb⁺ and Cs⁺) with cyclic peptides constructed from 3 or 4 alanine molecule (CyAla3 and CyAla4), has been investigated using mixed basis set (C, H, O, Li⁺, Na⁺ and K⁺ using 6-31+G(d), and the heavier cations: Rb⁺ and Cs⁺ using LANL2DZ). The minimum energy structures, binding energies, and various thermodynamic parameters of free ligands and their metal cations complexes have been determined with B3LYP and CAM-B3LYP functionals. The order of interaction energies were found to be Li⁺ > K⁺ > Na⁺ > Rb⁺ > Cs⁺ and Li⁺ > Na⁺ > K⁺ ? Rb⁺ > Cs⁺, calculated at CAM-B3LYP level for the M/CyAla3 and M/CyAla4 complexes, respectively. Their selectivity trend shows that the highest cation selectivity for Li⁺ over other alkali metal ions has been achieved on the basis of thermodynamic analysis. The main types of driving force host–guest interactions are investigated, the electron-donating O offers lone pair electrons to the contacting LP* of alkali metal cations.     

How accurate are TD-DFT excited-state geometries compared to DFT ground-state geometries?

In this work, we take a different angle to the benchmarking of time-dependent density functional theory (TD-DFT) for the calculation of excited-state geometries by extensively assessing how accurate such geometries are compared to ground-state geometries calculated with ordinary DFT. To this end, we consider 20 medium-sized aromatic organic compounds whose lowest singlet excited states are ideally suited for TD-DFT modeling and are very well described by the approximate coupled-cluster singles and doubles (CC2) method, and then use this method and six different density functionals (BP86, B3LYP, PBE0, M06-2X, CAM-B3LYP, and ωB97XD) to optimize the corresponding ground- and excited-state geometries. The results show that although each hybrid functional reproduces the CC2 excited-state bond lengths very satisfactorily, achieving an overall root mean square error of 0.011 Å for all 336 bonds in the 20 molecules, these errors are distinctly larger than those of only 0.004–0.006 Å with which the hybrid functionals reproduce the CC2 ground-state bond lengths. Furthermore, for each functional employed, the variation in the error relative to CC2 between different molecules is found to be much larger (by at least a factor of 3) for the excited-state geometries than for the ground-state geometries, despite the fact that the molecules/states under investigation have rather uniform chemical and spectroscopic character. Overall, the study finds that even in favorable circumstances, TD-DFT excited-state geometries appear intrinsically and comparatively less accurate than DFT ground-state ones.