Computational Design of an Electro-Organocatalyst for Conversion of CO2 into Formaldehyde

28 December 2022, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Density functional theory calculations employing a hybrid implicit/explicit solvation method were used to explore a new strategy for electrochemical conversion of CO2 using an electro-organocatalyst. A particular structural motif is identified that consists of an electron-rich >N−C=C−N< (enediamine) backbone which is capable of activating CO2 by formation of a C−C bond while subsequently facilitating the transfer of electrons from a chemically inert cathode to ultimately produce formaldehyde. Unlike transition metal-based electrocatalysts, the electroorganocatalyst is not constrained by scaling relations between the formation energies of activated CO2 and adsorbed CO, nor is it expected to be active for the competing hydrogen evolution reaction. The rate limiting steps are found to occur during two proton-coupled electron transfer (PCET) sequences and are associated with the transfer of a proton from a proton transfer mediator to a carbon atom on the electro-organocatalyst. The difficulty of this step in the second PCET sequence necessitates an electrode potential of −0.94 V vs. RHE to achieve the maximum turnover frequency. In addition, it is postulated that the electro-organocatalyst should also be capable of forming long chain aldehydes by successively carrying out reductive aldol condensation to grow the alkyl chain one carbon at a time. While much effort is still required to bring this conceptual design to reality, we are optimistic that the new directions for CO2 conversion opened up by this initial work will one day result in a practical device for electrochemical conversion of CO2 into multicarbon products.

Supplementary materials

Title
Description
Actions
Title
Computational Design of an Electro-Organocatalyst for Conversion of CO2 into Formaldehyde
Description
Density functional theory calculations were performed using the Vienna Ab-initio Simulation Package (VASP) combined with the VASPsol extension for including an implicit electrolyte in the simulation. The exchange-correlation energy was computed at the generalized gradient approximation level using the Bayesian error estimation functional with van der Waals correlation (BEEF-vdW) 4. The wave functions were constructed using plane waves with an energy up to 400 eV, while the projector augmented wave (PAW) method was used to represent the oscillations in the core region. The Kohn-Sham orbital populations were determined using an error function distribution having a width of 0.2 eV and the energy was extrapolated to zero temperature. All species were placed in a 20×20×20 Å unit cell and only the Γ- point was used to sample the Brillouin zone.
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.