Clémence Corminboeuf

Induced magnetic fields in aromatic [n]-annulenes: interpretation of NICS tensor components

The components of nucleus-independent chemical shift (NICS) tensors for Dnhn-annulenes are discussed as indexes of the aromatic character of electronic π systems. The component corresponding to the principal axis perpendicular to the ring plane, NICSzz, is found to be a good measure for the characterisation of the π system of the ring. Isotropic NICS values at ring centres contain large influences from the σ system and from all three principal components of the NICS tensor. At large distances away from the ring center, NICSzz, which is dominated by contributions from the π system, characterizes NICS well.

Comprehensive Benchmarking of a Density-Dependent Dispersion Correction

Standard density functional approximations cannot accurately describe interactions between nonoverlapping densities. A simple remedy consists in correcting for the missing interactions a posteriori, adding an attractive energy term summed over all atom pairs. The density-dependent energy correction, dDsC, presented herein, is constructed from dispersion coefficients computed on the basis of a generalized gradient approximation to Becke and Johnsons exchange-hole dipole moment formalism. dDsC also relies on an extended Tang and Toennies damping function accounting for charge-overlap effects. The comprehensive benchmarking on 341 diverse reaction energies divided into 18 illustrative test sets validates the robust performance and general accuracy of dDsC for describing various intra- and intermolecular interactions. With a total MAD of 1.3 kcal mol(-1), B97-dDsC slightly improves the results of M06-2X and B2PLYP-D3 (MAD = 1.4 kcal mol(-1) for both) at a lower computational cost. The density dependence of both the dispersion coefficients and the damping function makes the approach especially valuable for modeling redox reactions and charged species in general.

Why are the Interaction Energies of Charge-Transfer Complexes Challenging for DFT?

The description of ground state charge-transfer complexes is highly challenging. Illustrative examples include large overestimations of charge-transfer by local and semilocal density functional approximations as well as inaccurate binding energies. It is demonstrated here that standard density functionals fail to accurately describe interaction energies of charge-transfer complexes not only because of the missing long-range exchange as generally assumed but also as a result of the neglect of weak interactions. Thus, accounting for the missing van der Waals interactions is of key importance. These assertions, based on the evaluation of the extent of stabilization due to dispersion using both DFT coupled with our recent density-dependent dispersion correction (dDsC) and high-level ab initio computations, reflect the imperfect error-cancellation between the overestimation of charge-transfer and the missing long-range interactions. An in-depth energy decomposition analysis of an illustrative series of four small ambidentate molecules (HCN, HNC, HF, and ClF) bound together with NF3 provides the main conclusions, which are validated on a prototypical organic charge-transfer complex (i.e., tetrathiafulvalene-tetracyanoquinodimethane, TTF-TCNQ). We establish that the interaction energies for charge-transfer complexes can only be properly described when using well-balanced functionals such as PBE0-dDsC, M06-2X, and LC-BOP-LRD.

Visualizing and Quantifying Interactions in the Excited State

Determining the location and nature of the electron pairs within a molecule provides an intuitive representation of electronic structures. Yet, most of the available theoretical representations are not suitable for describing excited state phenomena. The density overlap region indicator (DORI) scalar field, which depends only on the density and its derivatives, overcomes previous limitations, while keeping the intuitiveness of popular scalar fields. Here, its usefulness is demonstrated by pinpointing visual and numerical DORI signatures for both intra- and intermolecular excited state situations.

A System-Dependent Density-Based Dispersion Correction.

Density functional approximations fail to provide a consistent description of weak molecular interactions arising from small electron density overlaps. A simple remedy to correct for the missing interactions is to add a posteriori an attractive energy term summed over all atom pairs in the system. The density-dependent energy correction, presented herein, is applicable to all elements of the periodic table and is easily combined with any electronic structure method, which lacks the accurate treatment of weak interactions. Dispersion coefficients are computed according to Becke and Johnsons exchange-hole dipole moment (XDM) formalism, thereby depending on the chemical environment of an atom (density, oxidation state). The long-range ∼R(-6) potential is supplemented with higher-order correction terms (∼R(-8) and ∼R(-10)) through the universal damping function of Tang and Toennies. A genuine damping factor depending on (iterative) Hirshfeld (overlap) populations, atomic ionization energies, and two adjustable parameters specifically fitted to a given DFT functional is also introduced. The proposed correction, dDXDM, dramatically improves the performance of popular density functionals. The analysis of 30 (dispersion corrected) density functionals on 145 systems reveals that dDXDM largely reduces the errors of the parent functionals for both inter- and intramolecular interactions. With mean absolute deviations (MADs) of 0.74-0.84 kcal mol(-1), PBE-dDXDM, PBE0-dDXDM, and B3LYP-dDXDM outperform the computationally more demanding and most recent functionals such as M06-2X and B2PLYP-D (MAD of 1.93 and 1.06 kcal mol(-1), respectively).

A generalized-gradient approximation exchange hole model for dispersion coefficients.

A simple method for computing accurate density-dependent dispersion coefficients is presented. The dispersion coefficients are modeled by a generalized gradient-type approximation to Becke and Johnsons exchange hole dipole moment formalism. Our most cost-effective variant, based on a disjoint description of atoms in a molecule, gives mean absolute errors in the C(6) coefficients for 90 complexes below 10%. The inclusion of the missing long-range van der Waals interactions in density functionals using the derived coefficients in a pair wise correction leads to highly accurate typical noncovalent interaction energies.

H2O2 Generation by Decamethylferrocene at a Liquid|Liquid Interface

Hydrogen peroxide generation at a liquid|liquid interface occurs with a yield of 20 % with respect to the concentration of reducing agent (decamethylferrocene). The liquid|liquid interface supplies electrons from the reducing agent and protons from the aqueous phase to drive the reduction of O2 into H2O2, which is extracted into the aqueous phase during the course of reaction (see picture; DCE=1,2-dichloroethane).

Oxygen Reduction Catalyzed by a Fluorinated Tetraphenylporphyrin Free Base at Liquid/Liquid Interfaces

The diprotonated form of a fluorinated free base porphyrin, namely 5-(p-aminophenyl)-10,15,20-tris(pentafluorophenyl)porphyrin (H(2)FAP), can catalyze the reduction of oxygen by a weak electron donor, namely ferrocene (Fc). At a water/1,2-dichloroethane interface, the interfacial formation of H(4)FAP(2+) is observed by UV-vis spectroscopy and ion-transfer voltammetry, due to the double protonation of H(2)FAP at the imino nitrogen atoms in the tetrapyrrole ring. H(4)FAP(2+) is shown to bind oxygen, and the complex in the organic phase can easily be reduced by Fc to produce hydrogen peroxide as studied by two-phase reactions with the Galvani potential difference between the two phases being controlled by the partition of a common ion. Spectrophotometric measurements performed in 1,2-dichloroethane solutions clearly evidence that reduction of oxygen by Fc catalyzed by H(4)FAP(2+) only occurs in the presence of the tetrakis(pentafluorophenyl)borate (TB(-)) counteranion in the organic phase. Finally, ab initio computations support the catalytic activation of H(4)FAP(2+) on oxygen.

Evaluation of aromaticity: A new dissected NICS model based on canonical orbitals

Several methods to address aromaticity in terms of nucleus-independent chemical shifts (NICS) are compared. These include NICS at the ring centre NICS(0), NICS 1 A above the ring plane NICS(1), aromatic ring current shielding (ARCS), and dissected NICS, i.e. NICS calculated from selected π orbitals NICSπ, again in the ring plane and 1 A above. The methods are tested on the basis of density-functional theory (DFT) and the individual gauge for local orbitals (IGLO) technique. Applications include simple organic rings (C4H4, C4H42+, C6H6, C5H5−, C7H7+) and transition metal carbonyl complexed molecules Fe(CO)3C4H4 and Cr(CO)3C6H6.

Simultaneous Visualization of Covalent and Noncovalent Interactions Using Regions of Density Overlap

We introduce a density-dependent bonding descriptor that enables simultaneous visualization of both covalent and noncovalent interactions. The proposed quantity is tailored to reveal the regions of space, where the total electron density results from a strong overlap of shell, atomic, or molecular densities. We show that this approach is successful in describing a variety of bonding patterns as well as nonbonding contacts. The Density Overlap Regions Indicator (DORI) analysis is also exploited to visualize and quantify the concept of electronic compactness in supramolecular chemistry. In particular, the scalar field is used to compare the compactness in molecular crystals, with a special emphasis on quaterthiophene derivatives with enhanced charge mobilities.