Consistent Inclusion of Continuum Solvation in Energy Decomposition Analysis: Theory and Application to Molecular CO2 Reduction Catalysts
Yuezhi Mao,Matthias Loipersberger,Kareesa J. Kron,Jeffrey S. Derrick,Christopher J. Chang,Shaama Mallikarjun Sharada,Martin Head-Gordon
DOI: https://doi.org/10.1039/d0sc05327a
IF: 8.4
2020-01-01
Chemical Science
Abstract:To facilitate computational investigation of intermolecular interactions in the solution phase, we report the developmentof ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework ofenergy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, allthe quantum mechanical states involved in the variational EDA procedure are computed with the presence of solventenvironment so that solvation effects are incorporated in the evaluation of all its energy components. After validation onseveral model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexesthat are related to electrochemical CO2reduction catalysis. For [FeTPP(CO2−κC)]2−(TPP = tetraphenylporphyrin), wereveal that two ortho substituents which yield most favorable CO2binding, −N(CH3)3+ (TMA) and −OH, stabilize thecomplex via through-structure and through-space mechanisms, respectively. The Coulombic interaction between thepositively charged TMA group and activated CO2is found to be largely attenuated by the polar solvent. Furthermore, wealso provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO2via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible “pockets” in that electrostatics isunscreened. For the reactant and product complexes associated with the electron transfer from the p-terphenyl radicalanion to CO2, we demonstrate that the double terminal substitution of p-terphenyl by electron-withdrawing groupsconsiderably strengthens the binding in the product state while moderately weakens that in the reactant state, which areboth dominated by the substituent tuning of the electrostatics component. These applications illustrate that this newextension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantifytheir impacts on chemical reactivity in solution.