Rebecca Sure

Corrected small basis set Hartree‐Fock method for large systems

A quantum chemical method based on a Hartree‐Fock calculation with a small Gaussian AO basis set is presented. Its main area of application is the computation of structures, vibrational frequencies, and noncovalent interaction energies in huge molecular systems. The method is suggested as a partial replacement of semiempirical approaches or density functional theory (DFT) in particular when self‐interaction errors are acute. In order to get accurate results three physically plausible atom pair‐wise correction terms are applied for London dispersion interactions (D3 scheme), basis set superposition error (gCP scheme), and short‐ranged basis set incompleteness effects. In total nine global empirical parameters are used. This so‐called Hartee‐Fock‐3c (HF‐3c) method is tested for geometries of small organic molecules, interaction energies and geometries of noncovalently bound complexes, for supramolecular systems, and protein structures. In the majority of realistic test cases good results approaching large basis set DFT quality are obtained at a tiny fraction of computational cost.

Blind Prediction of Binding Affinities for Charged Supramolecular Host–Guest Systems: Achievements and Shortcomings of DFT-D3

Association free energies ΔGa are calculated for two different types of host-guest systems, the rigid cucurbit[7]uril (CB7) and the basket shaped octa-acid (OA), and a number of charged guest molecules each by quantum chemical methods from first principles in the context of a recent blind test challenge (SAMPL4). For CB7, the overall agreement between theory and experiment is excellent. In comparison with all other submitted calculated relative ΔGa,rel values for this part of the blind test, our results ranked on top. Modeling the binding free energy in the case of the OA host mainly suffers from the problem that the binding situation is undefined with respect to the charge state and due to its intrinsic flexibility the host-guest complex is not represented well by a single configuration, but qualitative features of the binding process such as the proper binding orientation and the order of magnitude of ΔGa are represented in accord with the experimental expectations even though an accurate ranking is not possible.

Substituent effects and supramolecular interactions of titanocene(III) chloride: implications for catalysis in single electron steps.

The electrochemical properties of titanocene(III) complexes and their stability in THF in the presence and absence of chloride additives were studied by cyclic voltammetry (CV) and computational methods. The anodic peak potentials of the titanocenes can be decreased by as much as 0.47 V through the addition of an electron-withdrawing substituent (CO2Me or CN) to the cyclopentadienyl ring when compared with Cp2TiCl. For the first time, it is demonstrated that under the conditions of catalytic applications low-valent titanocenes can decompose by loss of the substituted ligand. The recently discovered effect of stabilizing titanocene(III) catalysts by chloride additives was analyzed by CV, kinetic, and computational studies. An unprecedented supramolecular interaction between [(C5H4R)2TiCl2](-) and hydrochloride cations through reversible hydrogen bonding is proposed as a mechanism for the action of the additives. This study provides the critical information required for the rational design of titanocene-catalyzed reactions in single electron steps.

Mechanistic study of the titanocene(III)-catalyzed radical arylation of epoxides.

An atom-economical and catalytic arylation of epoxide-derived radicals is described. The key step of the catalytic system is a sequential electron and proton transfer for the rearomatization of the radical σ-complex and catalyst regeneration. Kinetic, computational, spectroscopic, and cyclovoltammetric investigations highlight the key issues of the reaction mechanism and catalyst stabilization by collidine hydrochloride. Studies employing radicophiles rule out the participation of cations as reactive intermediates.

Small Atomic Orbital Basis Set First-Principles Quantum Chemical Methods for Large Molecular and Periodic Systems: A Critical Analysis of Error Sources

Abstract In quantum chemical computations the combination of Hartree–Fock or a density functional theory (DFT) approximation with relatively small atomic orbital basis sets of double‐zeta quality is still widely used, for example, in the popular B3LYP/6‐31G* approach. In this Review, we critically analyze the two main sources of error in such computations, that is, the basis set superposition error on the one hand and the missing London dispersion interactions on the other. We review various strategies to correct those errors and present exemplary calculations on mainly noncovalently bound systems of widely varying size. Energies and geometries of small dimers, large supramolecular complexes, and molecular crystals are covered. We conclude that it is not justified to rely on fortunate error compensation, as the main inconsistencies can be cured by modern correction schemes which clearly outperform the plain mean‐field methods.

A systematic study of rare gas atoms encapsulated in small fullerenes using dispersion corrected density functional theory

The most stable fullerene structures from C20 to C60 are chosen to study the energetics and geometrical consequences of encapsulating the rare gas elements He, Ne, or Ar inside the fullerene cage using dispersion corrected density functional theory. An exponential increase in stability is found with increasing number of carbon atoms. A similar exponential law is found for the volume expansion of the cage due to rare gas encapsulation with decreasing number of carbon atoms. We show that dispersion interactions become important with increasing size of the fullerene cage, where Van der Waals forces between the rare gas atom and the fullerene cage start to dominate over repulsive interactions. The smallest fullerenes where encapsulation of a rare gas element is energetically still favorable are He@C48, Ne@C52, and Ar@C58. While dispersion interactions follow the trend Ar > Ne > He inside C60 due to the trend in the rare gas dipole polarizabilities, repulsive forces become soon dominant with smaller cage size and we have a complete reversal for the energetics of rare gas encapsulation at C50.

Computational study of the rate constants and free energies of intramolecular radical addition to substituted anilines.

Summary The intramolecular radical addition to aniline derivatives was investigated by DFT calculations. The computational methods were benchmarked by comparing the calculated values of the rate constant for the 5-exo cyclization of the hexenyl radical with the experimental values. The dispersion-corrected PW6B95-D3 functional provided very good results with deviations for the free activation barrier compared to the experimental values of only about 0.5 kcal mol−1 and was therefore employed in further calculations. Corrections for intramolecular London dispersion and solvation effects in the quantum chemical treatment are essential to obtain consistent and accurate theoretical data. For the investigated radical addition reaction it turned out that the polarity of the molecules is important and that a combination of electrophilic radicals with preferably nucleophilic arenes results in the highest rate constants. This is opposite to the Minisci reaction where the radical acts as nucleophile and the arene as electrophile. The substitution at the N-atom of the aniline is crucial. Methyl substitution leads to slower addition than phenyl substitution. Carbamates as substituents are suitable only when the radical center is not too electrophilic. No correlations between free reaction barriers and energies (ΔG ‡ and ΔG R) are found. Addition reactions leading to indanes or dihydrobenzofurans are too slow to be useful synthetically.

Highly Active Titanocene Catalysts for Epoxide Hydrosilylation: Synthesis, Theory, Kinetics, EPR Spectroscopy

A catalytic system for titanocene-catalyzed epoxide hydrosilylation is described. It features a straightforward preparation of titanocene hydrides that leads to a reaction with low catalyst loading, high yields, and high selectivity of radical reduction. The mechanism was studied by a suite of methods, including kinetic studies, EPR spectroscopy, and computational methods. An unusual resting state leads to the observation of an inverse rate order with respect to the epoxide.

Synthesis, Chiral Resolution, and Absolute Configuration of Dissymmetric 4,15-Difunctionalized [2.2]Paracyclophanes

Despite the fact that functionalized planar chiral [2.2]paracyclophanes have received a lot of attention, the chemistry of pseudo-meta 4,15-distubstituted [2.2]paracyclophanes is largely unexplored. This is mainly due to the fact that the 4,5-dibromo-functionalized [2.2]paracyclophane is much less prone to halogen-metal exchange reactions than its constitutional pseudo-ortho or pseudo-para isomers. Here, we give an account of an efficient protocol to achieve this, which allows the synthesis of a broad variety of 4,15-disubstituted [2.2]paracyclophanes. Furthermore, we were able to resolve several of the racemic compounds via chiral HPLC and assign the absolute configurations of the isolated enantiomers by X-ray diffraction and/or by the comparison of calculated and measured CD-spectra.

From small fullerenes to the graphene limit: A harmonic force-field method for fullerenes and a comparison to density functional calculations for Goldberg-Coxeter fullerenes up to C980

We introduce a simple but computationally very efficient harmonic force field, which works for all fullerene structures and includes bond stretching, bending, and torsional motions as implemented into our open‐source code Fullerene. This gives accurate geometries and reasonably accurate vibrational frequencies with root mean square deviations of up to 0.05 Å for bond distances and 45.5 cm−1 for vibrational frequencies compared with more elaborate density functional calculations. The structures obtained were used for density functional calculations of Goldberg–Coxeter fullerenes up to C980. This gives a rather large range of fullerenes making it possible to extrapolate to the graphene limit. Periodic boundary condition calculations using density functional theory (DFT) within the projector augmented wave method gave an energy difference between −8.6 and −8.8 kcal/mol at various levels of DFT for the reaction C60→graphene (per carbon atom) in excellent agreement with the linear extrapolation to the graphene limit (−8.6 kcal/mol at the Perdew–Burke–Ernzerhof level of theory).