Carlo Adamo

TOWARD RELIABLE DENSITY FUNCTIONAL METHODS WITHOUT ADJUSTABLE PARAMETERS :THE PBE0 MODEL

We present an analysis of the performances of a parameter free density functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amount of exact exchange. The results obtained for structural, thermodynamic, kinetic and spectroscopic (magnetic, infrared and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chemistry and condensed matter physics.

Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters: The mPW and mPW1PW models

Starting from an analysis of the low-density and large gradient regions which dominate van der Waals interactions, we propose a modification of the exchange functional introduced by Perdew and Wang, which significantly enlarges its field of applications. This is obtained without increasing the number of adjustable parameters and retaining all the asymptotic and scaling properties of the original model. Coupling the new exchange functional to the correlation functional also proposed by Perdew and Wang leads to the mPWPW model, which represents the most accurate generalized gradient approximation available until now. We next introduce an adiabatic connection method in which the ratio between exact and density functional exchange is determined a priori from purely theoretical considerations and no further parameters are present. The resulting mPW1PW model allows to obtain remarkable results both for covalent and noncovalent interactions in a quite satisfactory theoretical framework encompassing the free elec...

Extensive TD-DFT Benchmark: Singlet-Excited States of Organic Molecules

Extensive Time-Dependent Density Functional Theory (TD-DFT) calculations have been carried out in order to obtain a statistically meaningful analysis of the merits of a large number of functionals. To reach this goal, a very extended set of molecules (∼500 compounds, >700 excited states) covering a broad range of (bio)organic molecules and dyes have been investigated. Likewise, 29 functionals including LDA, GGA, meta-GGA, global hybrids, and long-range-corrected hybrids have been considered. Comparisons with both theoretical references and experimental measurements have been carried out. On average, the functionals providing the best match with reference data are, one the one hand, global hybrids containing between 22% and 25% of exact exchange (X3LYP, B98, PBE0, and mPW1PW91) and, on the other hand, a long-range-corrected hybrid with a less-rapidly increasing HF ratio, namely LC-ωPBE(20). Pure functionals tend to be less consistent, whereas functionals incorporating a larger fraction of exact exchange tend to underestimate significantly the transition energies. For most treated cases, the M05 and CAM-B3LYP schemes deliver fairly small deviations but do not outperform standard hybrids such as X3LYP or PBE0, at least within the vertical approximation. With the optimal functionals, one obtains mean absolute deviations smaller than 0.25 eV, though the errors significantly depend on the subset of molecules or states considered. As an illustration, PBE0 and LC-ωPBE(20) provide a mean absolute error of only 0.14 eV for the 228 states related to neutral organic dyes but are completely off target for cyanine-like derivatives. On the basis of comparisons with theoretical estimates, it also turned out that CC2 and TD-DFT errors are of the same order of magnitude, once the above-mentioned hybrids are selected.

TD-DFT Performance for the Visible Absorption Spectra of Organic Dyes: Conventional versus Long-Range Hybrids.

The π → π* transitions of more than 100 organic dyes from the major classes of chromophores (quinones, diazo, ...) have been investigated using a Time-Dependent Density Functional Theory (TD-DFT) procedure relying on large atomic basis sets and the systematic modeling of solvent effects. These calculations have been performed with pure (PBE) as well as conventional (PBE0) and long-range (LR) corrected hybrid functionals (LC-PBE, LC-ωPBE, and CAM-B3LYP). The computed wavelengths are systematically guided by the percentage of exact exchange included at intermediate interelectronic distance, i.e., the λmax value always follows the PBE > PBE0 > CAM-B3LYP > LC-PBE > LC-ωPBE > HF sequence. The functional giving the best estimates of the experimental transition energies may vary, but PBE0 and CAM-B3LYP tend to outperform all other approaches. The latter functional is shown to be especially adequate to treat molecules with delocalized excited states. The mean absolute error provided by PBE0 is 22 nm (0.14 eV) with no deviation exceeding 100 nm (0.50 eV):  PBE0 is able to deliver reasonable estimates of the color of most organic dyes of practical or industrial interest. By using a calibration curve, we found that the LR functionals systematically allow an even more consistent description of the low-lying excited-state energies than the conventional hybrids. Indeed, linearly corrected LR approaches yield an average error of 10 nm for each dye family. Therefore, when such statistical treatments can be designed for given sets of dyes, a simple and rapid theoretical procedure allows both a chemically sound and a numerically accurate description of the absorption wavelengths.

TOWARD RELIABLE ADIABATIC CONNECTION MODELS FREE FROM ADJUSTABLE PARAMETERS

Abstract Methods rooted in the adiabatic connection theorem offer a significant improvement over standard density functional approximations. Among them, the most successful model is the so-called B3LYP functional, which includes three empirical parameters optimized by a best fit to atomization energies. We introduce a related functional, B1LYP, in which the ratio between Hartree-Fock and density functional exchange is determined a priori from purely theoretical considerations and no further parameters are present. The numerical results obtained on a standard molecular data set and on some ‘delicate’ chemical systems indicate that this model has essentially the same performance as the B3LYP model, but within a more satisfactory theoretical framework.

Accurate excitation energies from time-dependent density functional theory: Assessing the PBE0 model

We analyze the performance of a density functional model obtained by combining the Perdew–Burke–Erzenrhof (PBE) generalized gradient functional with a predetermined amount of exact exchange for predicting vertical electronic excitation energies within a time-dependent approach. Four molecules, namely, CO, H2CO, (CH3)2CO, and C2H4, have been chosen as benchmark cases. Our results show that this model (PBE0) provides accurate excitations both to valence and Rydberg states. Furthermore, the results are numerically close to those obtained using asymptotically correct exchange-correlation functionals. The performance of the PBE0 model for predicting excitation energies in larger molecules is assessed for benzene, pyridine, and naphthalene. Here, the PBE0 model provides results which are in fairly good agreement with experimental data and of similar quality to those obtained by more sophisticated (and time-consuming) post-Hartree–Fock methods.

A Qualitative Index of Spatial Extent in Charge-Transfer Excitations

With the aim of defining the spatial extent associated to an electronic transition, of particular relevance in the case of charge-transfer (CT) excitations, a new index, evaluated only from the computed density for the ground and excited state, is here derived and tested on a family of molecules that can be considered as prototypes of push-pull chromophores.The index (DCT) allows to define the spatial extent associated to a given transition as well as the associated fraction of electron transferred. By definition of centroids of charges associated to the density increase and depletion zones upon excitation, a qualitative and easy to visualize measure of the spatial extent of the donor and the acceptor moieties within a given molecular system is also given. Finally, an index (t) allowing to define the presence eventually pathologic CT transitions for time-dependent density functional theory treatment in conjunction with standard generalized gradient approximation or hybrid functional, that is through space CT, is disclosed.

Accurate simulation of optical properties in dyes.

Since Antiquity, humans have produced and commercialized dyes. To this day, extraction of natural dyes often requires lengthy and costly procedures. In the 19th century, global markets and new industrial products drove a significant effort to synthesize artificial dyes, characterized by low production costs, huge quantities, and new optical properties (colors). Dyes that encompass classes of molecules absorbing in the UV-visible part of the electromagnetic spectrum now have a wider range of applications, including coloring (textiles, food, paintings), energy production (photovoltaic cells, OLEDs), or pharmaceuticals (diagnostics, drugs). Parallel to the growth in dye applications, researchers have increased their efforts to design and synthesize new dyes to customize absorption and emission properties. In particular, dyes containing one or more metallic centers allow for the construction of fairly sophisticated systems capable of selectively reacting to light of a given wavelength and behaving as molecular devices (photochemical molecular devices, PMDs).Theoretical tools able to predict and interpret the excited-state properties of organic and inorganic dyes allow for an efficient screening of photochemical centers. In this Account, we report recent developments defining a quantitative ab initio protocol (based on time-dependent density functional theory) for modeling dye spectral properties. In particular, we discuss the importance of several parameters, such as the methods used for electronic structure calculations, solvent effects, and statistical treatments. In addition, we illustrate the performance of such simulation tools through case studies. We also comment on current weak points of these methods and ways to improve them.

On the Performances of the M06 Family of Density Functionals for Electronic Excitation Energies

We assessed the accuracy of the four members of the M06 family of functionals (M06-L, M06, M06-2X, and M06-HF) for the prediction of electronic excitation energies of main-group compounds by time-dependent density functional theory. This is accomplished by comparing the predictions both to high-level theoretical benchmark calculations and some experimental data for gas-phase excitation energies of small molecules and to experimental data for midsize and large chromogens in liquid-phase solutions. The latter comparisons are carried out using implicit solvation models to include the electrostatic effects of solvation. We find that M06-L is one of the most accurate local functionals for evaluating electronic excitation energies, that M06-2X outperforms BHHLYP, and that M06-HF outperforms HF, although in each case, the compared functionals have the same or a similar amount of Hartree-Fock exchange. For the majority of investigated excited states, M06 emerges as the most accurate functional among the four tested, and it provides an accuracy similar to the best of the other global hybrids such as B3LYP, B98, and PBE0. For 190 valence excited states, 20 Rydberg states, and 16 charge transfer states, we try to provide an overall assessment by comparing the quality of the predictions to those of time-dependent Hartree-Fock theory and nine other density functionals. For the valence excited states, M06 yields a mean absolute deviation (MAD) of 0.23 eV, whereas B3LYP, B98, and PBE0 have MADs in the range 0.19-0.22 eV. Of the functionals tested, M05-2X, M06-2X, and BMK are found to perform best for Rydberg states, and M06-HF performs best for charge transfer states, but no single functional performs satisfactorily for all three kinds of excitation. The performance of functionals with no Hartree-Fock exchange is of great practical interest because of their high computational efficiency, and we find that M06-L predicts more accurate excitation energies than other such functionals.

TD-DFT Assessment of Functionals for Optical 0-0 Transitions in Solvated Dyes.

Using TD-DFT, we performed simulations of the adiabatic energies of 40 fluorescent molecules for which the experimental 0-0 energies in condensed phase are available. We used six hybrid functionals (B3LYP, PBE0, M06, M06-2X, CAM-B3LYP, and LC-PBE) that have been shown to provide accurate transition energies in previous TD-DFT assessments, selected two diffuse-containing basis sets, and applied the most recent models for estimating bulk solvation effects. In each case, the correction arising from the difference of zero-point vibrational energies between the ground and the excited states has been consistently determined. Basis set effects have also been carefully studied. It turned out that PBE0 and M06 are the most effective functionals in terms of average deviation (mean absolute error of 0.22-0.23 eV). However, both the M06-2X global hybrid that contains more exact exchange and the CAM-B3LYP range-separated hybrid significantly improve the consistency of the prediction for a relatively negligible degradation of the average error. In addition, we assessed (1) the cross-structure/spectra relationships, (2) the importance of solvent effects, and (3) the differences between adiabatic and vertical energies.