These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
EDA_PCM_ChemRxiv.pdf (5.71 MB)

Consistent Inclusion of Continuum Solvation in Energy Decomposition Analysis: Theory and Application to Molecular CO2 Reduction Catalysts

submitted on 08.10.2020 and posted on 08.10.2020 by Yuezhi Mao, Matthias Loipersberger, Kareesa Kron, Jeffrey Derrick, Christopher Chang, Shaama Mallikarjun Sharada, Martin Head-Gordon

To facilitate computational investigation of intermolecular interactions in the solution phase, we report the development of ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework of energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, all the quantum mechanical states involved in the variational EDA procedure are computed with the presence of solvent environment so that solvation effects are incorporated in the evaluation of all its energy components. After validation on several model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexes that are related to electrochemical CO2 reduction catalysis. For [FeTPP(CO2−κC)]2− (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO2 binding, −N(CH3)3+ (TMA) and −OH, stabilize the complex via through-structure and through-space mechanisms, respectively. The Coulombic interaction between the positively charged TMA group and activated CO2 is found to be largely attenuated by the polar solvent. Furthermore, we also provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO2 via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible “pockets” in that electrostatics is unscreened. For the reactant and product complexes associated with the electron transfer from the p-terphenyl radical anion to CO2 , we demonstrate that the double terminal substitution of p-terphenyl by electron-withdrawing groups considerably strengthens the binding in the product state while moderately weakens that in the reactant state, which are both dominated by the substituent tuning of the electrostatics component. These applications illustrate that this new extension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantify their impacts on chemical reactivity in solution.


Email Address of Submitting Author


Stanford University


United States

ORCID For Submitting Author


Declaration of Conflict of Interest

MHG is a part owner of Q-Chem Inc.


Logo branding