Computational Modelling — CO2 Reduction to Methane and Ethylene on a Single-Atom Catalyst: A Grand Canonical Quantum Mechanics Study

Measurement evidence

Computational Modelling

CO2 Reduction to Methane and Ethylene on a Single-Atom Catalyst: A Grand Canonical Quantum Mechanics Study · Osella S., Goddard III W.A. · Journal of the American Chemical Society · 2023 · 21319-21329

10 measurement groups · 52 results

Reported values remain attached to the sample, method, conditions, extraction quality and source location that produced them.

GCP-K adsorption-energy calculations and PDOS d-band-centre analysis for CuPc and CuO4 sites.

PcCu 2D MOF reaction model with adsorbates and water cluster · Model

Adsorption energies at zero bias and applied-potential trends; d-band centres from Figure S9.

Atmosphere
Implicit solvent for adsorption calculations.
Geometry
Periodic monolayer reaction model.
Context
Adsorbate-bearing model system.
Measurement source
main p.8, article p.21326 · 2.4 Optimizing the Ethylene Production Yield · Figure 5d; Figure S9
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
CO adsorption energy at first CuO4 site, zero external potential-0.20 eVText
Exact Reported
main p.8, article p.21326 · 2.4 Optimizing the Ethylene Production Yield · Figure 5d
CO adsorption energy at second CuO4 site, zero external potential-0.18 eVText
Exact Reported
main p.8, article p.21326 · 2.4 Optimizing the Ethylene Production Yield · Figure 5d
d-band centre for 2CO-CuO4 site-2.92 eVText
Exact Reported
main p.8, article p.21326 · 2.4 Optimizing the Ethylene Production Yield · Figure S9
d-band centre for CO-CuO4 site-2.85 eVText
Exact Reported
main p.8, article p.21326 · 2.4 Optimizing the Ethylene Production Yield · Figure S9
d-band centre for CO-CuPc siteMarked as a best value within this paper-2.64 eVText
Exact Reported
main p.8, article p.21326 · 2.4 Optimizing the Ethylene Production Yield · Figure S9
CO2 adsorption energy at PcCu catalytic center, zero external potential-0.12 eVText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5d
CO adsorption energy at PcCu catalytic center, zero external potential-0.11 eVText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5d
H adsorption energy at PcCu catalytic center, zero external potential1.85 eVText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5d

QM geometry optimisation of literature PcCu MOF cell.

Reoptimised pristine PcCu 2D MOF monolayer model cell · Model

Optimised monolayer with potassium counterions and S = 3/2 high-spin spin polarisation.

Geometry
Periodic monolayer cell.
Context
Pristine model system.
Measurement source
main p.8, article p.21326 · 3. Experimental Section
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
Optimised cell parameter aa = 18.17536 AText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section
Optimised cell angle alphaalpha = 91.3834 degText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section
Optimised cell parameter bb = 18.17632 AText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section
Optimised cell angle betabeta = 91.0939 degText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section
Optimised cell parameter cc = 20.20189 AText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section
Optimised cell angle gammagamma = 89.8539 degText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section

Spin-polarised DFT in VASP with VASPsol implicit solvation, PBE functional, 500 eV plane-wave cutoff, Grimme D3(BJ), 1 x 1 x 1 Monkhorst-Pack k-point grid; single-point JDFTx at applied potentials with CANDLE implicit solvation; GCP-K/GC-QM Legendre transform.

Reoptimised pristine PcCu 2D MOF monolayer model cell · Model

Geometry optimisations used 1.5 nm z-vacuum, force convergence 0.02 eV/A, electronic convergence 1e-5 eV; vibrational corrections at 298.15 K.

Temperature
298.15
Atmosphere
Implicit water solvent; computational vacuum spacing normal to MOF.
Geometry
Periodic monolayer cell.
Context
Pristine model system.
Measurement source
main p.8, article p.21326 · 3. Experimental Section
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
Geometry optimisation force convergence criterion0.02 eV/AText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section
Plane-wave cutoff energy500 eVText
Exact Reported
main p.8, article p.21326 · 3. Experimental Section

DFT density of states and projected density of states (DOS/PDOS).

Reoptimised pristine PcCu 2D MOF monolayer model cell · Model

Optimised pristine MOF monolayer model; electronic properties computed after geometry optimisation.

Geometry
Periodic monolayer cell.
Context
Pristine model system.
Measurement source
main p.3, article p.21321 · 2.1 Electronic Properties · Figure 1b,c
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
Calculated semiconductor band gapMarked as a best value within this paper~0.3 eVText
Approximate
main p.3, article p.21321 · 2.1 Electronic Properties · Figure 1b
Pristine MOF total spin stateS = 3/2Text
Exact Reported
main p.3, article p.21321 · 2.1 Electronic Properties · Figure 1c

GCP-K free-energy pathway for ethylene formation through *(CH)CO.

PcCu 2D MOF reaction model with adsorbates and water cluster · Model

Computed at 298 K, neutral pH, and applied potential -1.2 V_RHE.

Temperature
298
Atmosphere
Implicit solvent plus protonated four-water cluster.
Geometry
Periodic monolayer reaction model.
Context
Adsorbate-bearing model system.
Measurement source
main p.5, article p.21323 · 2.3 Ethylene Formation Pathways · Figure 3
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
Carbon-protonation branch RDS barrier via *(CH2)COH0.82 eVText
Exact Reported
main p.5, article p.21323 · 2.3 Ethylene Formation Pathways · Figure 3
Relative energy of *(CH)CO branchDeltaG = -4.69 eVText
Exact Reported
main p.5, article p.21323 · 2.3 Ethylene Formation Pathways · Figure 3
Relative energy of second *(CH)CO branchDeltaG = -5.37 eVText
Exact Reported
main p.5, article p.21323 · 2.3 Ethylene Formation Pathways · Figure 3
Water-elimination branch RDS barrier via *(CH)COH1.01 eVText
Exact Reported
main p.5, article p.21323 · 2.3 Ethylene Formation Pathways · Figure 3

GCP-K free-energy pathway for ethylene formation through *(CHO)CO.

PcCu 2D MOF reaction model with adsorbates and water cluster · Model

Computed at 298 K, neutral pH, and applied potential -1.2 V_RHE.

Temperature
298
Atmosphere
Implicit solvent plus protonated four-water cluster.
Geometry
Periodic monolayer reaction model.
Context
Adsorbate-bearing model system.
Measurement source
main p.6, article p.21324 · 2.3 Ethylene Formation Pathways · Figure 4
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
Stabilisation on CO coupling into *CHO to form *(CHO)COMarked as a best value within this paperDeltaG = -0.93 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
Preferred *(CHO)CO pathway RDS barrierMarked as a best value within this paper0.42 eVText
Exact Reported
main p.6, article p.21324 · 2.3 Ethylene Formation Pathways · Figure 4

GCP-K free-energy pathway for methane formation from CO2 on PcCu MOF.

PcCu 2D MOF reaction model with adsorbates and water cluster · Model

Free energies in eV at 298 K; neutral pH; applied potential -1.2 V_RHE unless noted.

Temperature
298
Atmosphere
Implicit solvent plus protonated four-water cluster.
Geometry
Periodic monolayer reaction model.
Context
Adsorbate-bearing model system.
Measurement source
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
*CO to *CHO activation barrier0.22 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
Methane binding energy on surfaceDeltaG = 0.01 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
Initial CO2 adsorption free energyDeltaG = -0.12 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
CO2 bending activation barrier0.32 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
*CO binding energyDeltaG = -0.11 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
Overall stability of *CO relative to free CO2DeltaG = -2.27 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
Final methane formation free energy relative to gas-phase CO2DeltaG = -8.78 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
First proton-electron transfer barrier to *HOCO0.75 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2
Methane-pathway rate-determining barrier to *CHOH0.92 eVText
Exact Reported
main p.4, article p.21322 · 2.2 Methane Formation Pathway · Figure 2

GCP-K potential-dependent free energies for selected intermediates and transition states.

PcCu 2D MOF reaction model with adsorbates and water cluster · Model

Applied-potential series compared with prior experimental FE trends; Figure 5.

Atmosphere
Implicit solvent.
Geometry
Periodic monolayer reaction model.
Context
Adsorbate-bearing model system.
Measurement source
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
*CHOH-formation barrier at -1.0 V_RHE1.00 eVText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5a
*CHOH-formation barrier at -1.6 V_RHE0.74 eVText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5a
*CHOH energy relative to *CHO at -1.0 V_RHE0.40 eV higherText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5a,c
*CHOH energy relative to *CHO at -1.6 V_RHE0.14 eVText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5a,c
Applied potential where *CHOH population exceeds *CHO-1.96 VText
Exact Reported
main p.7, article p.21325 · 2.4 Optimizing the Ethylene Production Yield · Figure 5c

Optimised intermediate geometry; tabulated Cu-C and C-N distances.

PcCu 2D MOF reaction model with adsorbates and water cluster · Model

Different intermediates at -1.2 V.

Geometry
Periodic monolayer reaction model.
Context
Adsorbate-bearing model system.
Measurement source
SI p.S13 · Supporting Information · Table S4
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
Cu-C distance for *CCH intermediateCu-C = 1.91 AText
Exact Reported
main p.6, article p.21324 · 2.3 Ethylene Formation Pathways · Figure 4; Table S4
C-N distance for *CCHOH stabilised structure1.44 AText
Exact Reported
main p.6, article p.21324 · 2.3 Ethylene Formation Pathways · Figure 4; Table S4
*CH C-N distance1.39 ASI Table
Exact Reported
SI p.S13 · Supporting Information · Table S4
*CH Cu-C distance1.94 ASI Table
Exact Reported
SI p.S13 · Supporting Information · Table S4
*(CH)CO Cu-C distances1.98/1.90 ASI Table
Exact Reported
SI p.S13 · Supporting Information · Table S4
*(CHO)CO Cu-C distance1.94 ASI Table
Exact Reported
SI p.S13 · Supporting Information · Table S4

GCP-K quadratic free-energy parameter extraction.

PcCu 2D MOF reaction model with adsorbates and water cluster · Model

GCP parameters a, b, c, differential capacitance, and potential of zero charge for investigated intermediates.

Geometry
Periodic monolayer reaction model.
Context
Adsorbate-bearing model system.
Measurement source
SI p.S15 · Supporting Information · Table S6
PropertyReported valueNormalised valueUncertaintyOrigin and qualitySource
Cu-MOF GCP parameter a0.1612960 V/eSI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6
Cu-MOF differential capacitanceMarked as a best value within this paper3.09989 FSI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6
Cu-MOF potential of zero chargeMarked as a best value within this paper-0.92623 V vs. SHESI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6
*(CH)CO potential of zero charge-0.98136 V vs. SHESI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6
*CHO potential of zero charge-0.90669 V vs. SHESI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6
*(CHO)CO potential of zero charge-1.12833 V vs. SHESI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6
*CHOH potential of zero charge-0.94594 V vs. SHESI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6
*CO2 potential of zero charge-1.18529 V vs. SHESI Table
Exact Reported
SI p.S15 · Supporting Information · Table S6