Full Papers
doi.org/10.1002/cssc.202100275
ChemSusChem
with ultrapure water and dried under Ar atmosphere. Finally, the
was carried out using (4×4) grid. The cutoff energy was set to
500 eV and 5000 eV for the orbitals, and the charge density,
respectively, and a Fermi level smearing of 0.1 eV was used. Spin
polarization and dipole correction have been involved for all the
1
2
3
4
5
6
7
8
9
0
1
2
3
4
5
6
7
8
9
0
1
2
3
4
5
6
7
8
9
0
1
2
3
4
5
6
7
8
9
0
1
2
3
4
5
6
7
8
9
0
1
2
3
4
5
6
7
electrode was electrochemically reduced in CO -bubbled 0.2 M
2
NaHCO at À 0.9 V (vs. SCE) for 2 h to obtain the Cu O-derived Cu
3
x
electrode. Cu-based electrodes prepared in the basic precursor
without any additions, and with different amount of KBr or CTAB
are designated as CuÀ Blank, CuÀ KBr(x), CuÀ CTAB(x), respectively,
where x represents the additive concentration.
calculations. To ensure computationally tractable, we use
a
simplified model of CTAB, C6H16N, with a shorter carbon chain
considered since the carbon chain has a negligible influence on the
ERC reaction. Potential asymmetric adsorption sites of each
intermediate were calculated and only the most favorable was
used.
Surface roughness determination
Electrochemical cyclic voltammetry was employed to test the
The Gibbs free energies of typical proton-coupled electron transfer
1
1
1
1
1
1
1
1
1
1
2
2
2
2
2
2
2
2
2
2
3
3
3
3
3
3
3
3
3
3
4
4
4
4
4
4
4
4
4
4
5
5
5
5
5
5
5
5
[18]
double layer charging capacitance of Cu O electrode to characterize
steps are estimated using the CHE approach, where the free
x
the electrochemical reaction surface area participating in the ERC
energy of the solvated proton and electron pair can be by
definition equal to one-half the chemical potential of molecular H2
at 101325 Pa, further as a function of applied potential at all
temperature and PH values. Within the CHE approach, the free
energy of intermediate is calculated by G=E+ZPE–TΔS, in which E
is the DFT-optimized total energy, ZPE is the zero-point vibrational
energy, T is the absolute temperature, and ΔS is the entropy
change. For the adsorbed intermediates, the ZPE and entropy were
calculated using the harmonic approximation with the ASE module
reaction. The test solution is N saturated 0.1 M HClO and the
2
4
potential window is À 0.24 V to À 0.34 V(vs. SCE).
Characterization
Powder X-ray diffraction (XRD) spectra were measured on a Rigaku
Ultima IV diffractometer with CuKα radiation (l=0.15418 nm)
operating at 40 kV and 40 mA over the two theta range of 5–100°
[19]
À 1
vibrations. The value of entropy of gas phase was taken from
[20]
with a scan speed of 10° min . Scanning electron microscopy
NIST,
while the ZPE of gas phase were also obtained from
(SEM) was recorded on JEOL JSM-6360 operating at an acceleration
voltage of 20 kV. Br element quantities were detected by EDX. The
IR spectra were recorded with a Fourier Transform Infra spectrom-
vibrational frequencies from ASE.
To simulate the chemical reaction in electrocatalysis, the effect of
solvent stabilization should be taken into account for adsorbed
species, as well as gas-phase error corrections for the standard free
energy deduced from the limitation of exchange-correlation func-
tional. It is noticed that the gas-phase error corrections are
functional dependence (BEEF-vdW functional in our case), where
the values are taken from reference [21], in which they implement
systematic investigation gas-phase error on few various exchange-
correlation functional used prevalently (PBE, PW91, RPBE, and BEEF-
vdW). Regarding the independent solvation correction of adsorbed
species involved in ERC, we use the result of previous work
obtained from the sensitivity analysis, which is 0 eV, À 0.1 eV,
À 1
eter (Nicolet iS50, Thermo Fisher, USA) in the 2000–400 cm
wavenumber region with ATR technique. H À TPR was carried out
2
on Micromeritics Chemi 2920 II. 10%H /90%Ar stream with flow
2
À 1
rate of 30 mLmin was used to determine the Cu O species and
x
the number of the surface active sites, the temperature rise rate of
À 1
10
(
°
Cmin . The effluent of the H À TPR was detected by the MS
2
MKS104-J0316005) to detect the reduction product.
Electrochemical testing
The Electrochemical measurements were performed on a potentio-
stat (EG&G 2273, Princeton Applied Research) in a self-made H cell
[22]
À 0.25 eV for OCHO*, *CO and *COOH.
(
(
separated by Nafion115). Pt sheet and saturated calomel electrode
SCE) were used as the counter electrode and the reference
electrode, respectively. The prepared Cu-based electrode was used
as the working electrode. CO -saturated 0.2 M NaHCO aqueous
Acknowledgements
2
3
solution was used as the electrolyte. The generated gas products
were analyzed online by a gas chromatograph (GC, Shimadzu GC-
The authors thank the financial grants from the National Natural
Science Foundation of China (No. 21577141).
2014).The hydrogen and CO were quantified by a thermal
conductivity detector (TCD), and hydrocarbons were detected and
quantified by flame ionization detector (FID), respectively. The
À
liquid products including Br were quantified by ion chromatog-
Conflict of Interest
raphy (ICS-1100, Dionex Corporation).
The authors declare no conflict of interest.
Computational details
All DFT calculations were performed with the Quantum ESPRESSO
Keywords: carbon dioxide
·
copper
·
electrocatalysis
·
[15]
code employing the state-of-art Bayesian error estimation func-
hydrothermal synthesis · reduction
[16]
tional with van der Waals correlation (BEEF-vdW), which has been
successfully used for surface chemistry studies. We used soft
pseudopotentials and a vacuum region of 30 Å perpendicular to
the surface. Pseudopotentials for C, O, H, Cu, and N were generated
by using the “atomic” code reported by Dal Corso (v.5.0.2 svn rev.
[
[17]
9
415). In terms of our experimental observation, Cu O-derived
2
the typical (110) facet of Cu with and without CTAB involvement
was modeled using periodic slabs, where (3×3) unit cell and 4
layers were used. The bottom-most layer was fixed in their bulk
truncated positions, while the layers above were relaxed until
reaching a maximum force threshold of 0.03 eV/Å. k-point sampling
ChemSusChem 2021, 14, 1962–1969
1968
© 2021 Wiley-VCH GmbH