Accurate Thermochemistry for Large Molecules with Modern Density Functionals

Environmental Control of the Magnetic Behavior of Transition Metal Complexes: Density Functional Theory Study of Zeolite Y Embedded Complexes [M(bpy)3]2+@Y (M = Fe2+, Co2+)

Using the supramolecular approach developed for the study of the guest-host interactions in the zeolite Y encapsulated [Fe(bpy)₃]²⁺ compound: [Fe(bpy)₃]²⁺@Y (bpy = 2,2'-bipyridine) [Vargas et al., J. Chem. Theory Comput. 2009, 5, 97-115], we apply density functional theory (DFT) to the study of the influence of zeolite Y encapsulation on the structural and energetic properties of [Co(bpy)₃]²⁺ in the low-spin (LS) and high-spin (HS) states, while revisiting [Fe(bpy)₃]²⁺@Y. Although the accurate prediction of the HS-LS energy difference ΔE_HL^el remains challenging for current DFT methods, they give accurate estimates of its variation Δ(ΔE_HL^el) in a series of complexes of a given transition metal ion. Therefore, denoting [M(bpy)₃]²⁺@Y_SM as the supramolecular model of the inclusion compounds, the values of ΔE_HL^el for the bpy complexes in the gas phase and in the supercage of zeolite Y were determined by combining the DFT estimates of Δ(ΔE_HL^el) in the series {[M(NCH)₆]²⁺, [M(bpy)₃]²⁺, and [M(bpy)₃]²⁺@Y_SM}, with accurate CCSD(T) estimates of ΔE_HL^el in the benchmark complexes [M(NCH)₆]²⁺ (M = Fe, Co) [Lawson Daku et al., J. Chem. Theory Comput., 2012, 8, 4216-4231]. Generalized gradient approximations as well as global and range-separated hybrids were employed. In order to better account for the key role of dispersion, they were also augmented with the semiempirical D2, D3BJ, and D3BJM dispersion corrections when available. The use of the D3BJ and D3BJM corrections led to similar results, and this is only with the use of the D2 scheme that (i) the free and encapsulated [Fe(bpy)₃]²⁺ are correctly predicted as LS species and that (ii) the encapsulation of both complexes translates into a destabilization of their HS state with respect to their LS state. The increase of the HS-LS energy difference is smaller for [Co(bpy)₃]²⁺ than [Fe(bpy)₃]²⁺ because the HS-LS molecular volume difference ΔV_HL in [Co(bpy)₃]²⁺ is ∼50% smaller than in [Fe(bpy)₃]²⁺. Periodic DFT calculations performed on crystalline [M(bpy)₃]²⁺@Y show that the employed [M(bpy)₃]²⁺@Y_SM supramolecular model allows the influence of encapsulation on the geometry and the spin-state energetics of [M(bpy)₃]²⁺ (M = Fe, Co) to be quantitatively captured.

Light-driven reduction of CO2: thermodynamics and kinetics of hydride transfer reactions in benzimidazoline derivatives

CO₂ capture, conversion and storage belong to the holy grail of environmental science. We therefore explore an important photochemical hydride transfer reaction of benzimidazoline derivatives with CO₂ in a polar solvent (dimethylsulfoxide) by quantum-chemical methods. While the excited electronic state undergoing hydride transfer to formate (HCOO^-) shows a higher reaction path barrier compared to the ground state, a charge-transfer can occur in the near-UV region with nearly barrierless access to the products involving a conical intersection between both electronic states. Such radiationless decay through the hydride transfer reaction and formation of HCCO^-via excited electronic states in suitable organic compounds opens the way for future photochemical CO₂ reduction. We provide a detailed analysis for the chemical CO₂ reduction to the formate anion for 15 different benzimidazoline derivatives in terms of thermodynamic hydricities (ΔG_H^-), activation free energies (ΔG‡HT), and reaction free energies (ΔG_rxn) for the chosen solvent dimethylsulfoxide at the level of density functional theory. The calculated hydricities are in the range from 35.0 to 42.0 kcal mol^-1i.e. the species possess strong hydride donor abilities required for the CO₂ reduction to formate, characterized by relatively low activation free energies between 18.5 and 22.2 kcal mol^-1. The regeneration of the benzimidazoline can be achieved electrochemically.

Aryl-Cl vs heteroatom-Si bond cleavage on the route to the photochemical generation of σ,π-heterodiradicals

The photochemistry of aryl chlorides having a X-SiMe₃ group (X = O, NR, S, SiMe₂) tethered to the aromatic ring has been investigated in detail, with the aim to generate valuable ϭ,π-heterodiradicals. Two competitive pathways arising from the excited triplet state of the aromatics have been observed, namely heterolysis of the aryl-chlorine bond and homolysis of the X-silicon bond. The former path is found in chlorinated phenols and anilines, whereas the latter is exclusive in the case of silylated thiophenols and aryl silanes. A combined experimental/computational approach was pursued to explain such a photochemical behavior.Graphical abstract.

Interplay of Protecting Groups and Side Chain Conformation in Glycopyranosides. Modulation of the Influence of Remote Substituents on Glycosylation?

The synthesis and conformational analysis of a series of phenyl 2,3,6-tri-O-benzyl-β-d-thio galacto- and glucopyranosides and their 6S-deuterio isotopomers, with systematic variation of the protecting group at the 4-position, are described. For the galactopyranosides, replacement of a 4-O-benzyl ether by a 4-O-alkanoyl or aroyl ester results in a small but measurable shift in side chain population away from the trans,gauche conformation and in favor of the gauche,trans conformer. In the glucopyranoside series on the other hand, replacement of a 4-O-benzyl ether by a 4-O-alkanoyl or aroyl ester results in a small but measurable increase in the population of the trans,gauche conformer at the expense of the gauche,gauche conformer. The possible modulating effect of these conformational changes on the well-known changes in the anomeric reactivity of glycosyl donors as a function of protecting group is discussed, raising the possibility that larger changes may be observed at the transition state for glycosylation. A comparable study with a series of ethyl 2,3,4-tri-O-benzyl-β-d-thioglucopyranosides reveals that no significant influence in side chain population is observed on changing the O6 protecting group.

MN15-L: A New Local Exchange-Correlation Functional for Kohn-Sham Density Functional Theory with Broad Accuracy for Atoms, Molecules, and Solids

Kohn-Sham density functional theory is widely used for applications of electronic structure theory in chemistry, materials science, and condensed-matter physics, but the accuracy depends on the quality of the exchange-correlation functional. Here, we present a new local exchange-correlation functional called MN15-L that predicts accurate results for a broad range of molecular and solid-state properties including main-group bond energies, transition metal bond energies, reaction barrier heights, noncovalent interactions, atomic excitation energies, ionization potentials, electron affinities, total atomic energies, hydrocarbon thermochemistry, and lattice constants of solids. The MN15-L functional has the same mathematical form as a previous meta-nonseparable gradient approximation exchange-correlation functional, MN12-L, but it is improved because we optimized it against a larger database, designated 2015A, and included smoothness restraints; the optimization has a much better representation of transition metals. The mean unsigned error on 422 chemical energies is 2.32 kcal/mol, which is the best among all tested functionals, with or without nonlocal exchange. The MN15-L functional also provides good results for test sets that are outside the training set. A key issue is that the functional is local (no nonlocal exchange or nonlocal correlation), which makes it relatively economical for treating large and complex systems and solids. Another key advantage is that medium-range correlation energy is built in so that one does not need to add damped dispersion by molecular mechanics in order to predict accurate noncovalent binding energies. We believe that the MN15-L functional should be useful for a wide variety of applications in chemistry, physics, materials science, and molecular biology.

Quantum chemical approach to estimating the thermodynamics of metabolic reactions

Thermodynamics plays an increasingly important role in modeling and engineering metabolism. We present the first nonempirical computational method for estimating standard Gibbs reaction energies of metabolic reactions based on quantum chemistry, which can help fill in the gaps in the existing thermodynamic data. When applied to a test set of reactions from core metabolism, the quantum chemical approach is comparable in accuracy to group contribution methods for isomerization and group transfer reactions and for reactions not including multiply charged anions. The errors in standard Gibbs reaction energy estimates are correlated with the charges of the participating molecules. The quantum chemical approach is amenable to systematic improvements and holds potential for providing thermodynamic data for all of metabolism.