Theoretical Investigations into Defected Graphene ... - ACS Publications

Oct 10, 2017 - examine a variety of different active sites in N-doped graphene to give a ... energy source (such as wind or solar) has attracted much ...
0 downloads 0 Views 2MB Size
Research Article pubs.acs.org/journal/ascecg

Cite This: ACS Sustainable Chem. Eng. 2017, 5, 11080-11085

Theoretical Investigations into Defected Graphene for Electrochemical Reduction of CO2 Samira Siahrostami,*,† Kun Jiang,‡ Mohammadreza Karamad,† Karen Chan,§ Haotian Wang,‡ and Jens Nørskov†,§ †

SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, 443 Via Ortega, Stanford, California 94305, United States ‡ Rowland Institute, Harvard University, Cambridge, Massachusetts 02142, United States § SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States S Supporting Information *

ABSTRACT: Despite numerous experimental efforts that have been dedicated to studying carbon-based materials for electrochemical reduction of CO2, a rationalization of the associated trends in the intrinsic activity of different active motifs has so far been elusive. In the present work, we employ density functional theory calculations to examine a variety of different active sites in N-doped graphene to give a comprehensive outline of the trends in activity. We find that adsorption energies of COOH* and CO* do not follow the linear scaling relationships observed for the pure transition metals, and this unique scaling is rationalized through differences in electronic structure between transition metals and defected graphene. This finding rationalizes most of the experimental observations on the carbon-based materials which present promising catalysts for the two-electron reduction of CO2 to CO. With this simple thermodynamic analysis, we identify several active sites that are expected to exhibit a comparable or even better activity to the state-of-the-art gold catalyst, and several configurations are suggested to be selective for CO2RR over HER. KEYWORDS: Density functional theory (DFT), Calculated limiting potential, Hydrogen evolution reaction (HER), Scaling relation, Free energy diagram



compared to metal surfaces.25 Pristine graphene does not exhibit any activity. However, introducing dopants20,26,27 and defects28 during the synthesis tailors the electronic structure and catalytic properties of nanostructured carbon materials.28 N-doping has been shown to significantly enhance the CO2 reduction activity.29−33 Different types of naturally occurring defects can also modulate the catalytic activity.28,34,35 It is interesting to note that there is an interplay between N-doping and the number of generated defects.36 N-doping is energetically more favorable near the defects especially around the pentagons and represents the important role of defect-induced strain in the stability of the N-dopant.36 Figure 1 displays the performance of several carbon-based materials that have been tested and reported for the electrochemical reduction of CO2.29,30,37−41 The reported overpotentials and Faradaic efficiencies in Figure 1a correspond to the maximum collected CO in each report. The partial current

INTRODUCTION Global average carbon dioxide (CO2) concentration at the surfacethe main driver of recent climate changehas reached 400 ppm for the first time in recorded history. In recent years, several strategies have been proposed to mitigate atmospheric CO2 concentration, including CO2 capture, storage, and conversion into valuable products. Among these, electrochemical conversion and reduction of CO2 coupled with a renewable energy source (such as wind or solar) has attracted much attention as an approach to generate carbon-neutral fuels or industrial chemicals. Many potential electrocatalysts have been suggested, including transition metals (e.g., Cu,1−4 Au, and Ag), transition metal oxides (e.g., RuO2, oxide- derived copper),5−11 transition metal dichalcogenides (e.g., MoS2, MoSe2),12,13 metalfunctionalized porphyrin like structure,s14−17 and metal−organic frameworks.18,19 Carbon-based materials are particularly interesting potential catalysts for the CO2 reduction reaction (CO2RR) due to their low cost and ability to form a wide range of nanostructures.20−24 They are chemically inactive at negative potential ranges and present high overpotentials for the hydrogen evolution reaction © 2017 American Chemical Society

Received: August 31, 2017 Revised: October 6, 2017 Published: October 10, 2017 11080

DOI: 10.1021/acssuschemeng.7b03031 ACS Sustainable Chem. Eng. 2017, 5, 11080−11085

Research Article

ACS Sustainable Chemistry & Engineering

between each successive step at the equilibrium potential. The potential at which the electrochemical reaction takes place is called the limiting potential, which can be obtained by the difference between the equilibrium potential and overpotential. At the limiting potential, all the reaction steps are downhill in free energy.



RESULTS AND DISCUSSION Taking into account the two-dimensional graphene sheets as a model system, we introduce different types of defects as depicted in Figure 1a, including, mono vacancies (MV) as well as several double vacancies (DV), such as 555-777, 555-6-777, and Stone Wales 55-77. We also consider two different line defects including 558 and 555-777. Here, 5, 6, 7, and 8 refer to the pentagon, hexagon, heptagon, and octagon, respectively, in the defect configuration. Further details on the relative stability of different modeled defect configurations are provided in Table S2 in the Supporting Information. Our calculations indicate that MV is too reactive to contribute to the CO2 reduction catalysis, while some of the DVs can catalyze CO2RR. Using DFT calculations and the computational hydrogen electrode approach (SI), the CO2 reduction reaction mechanism has been explored and reported on transition metal surfaces.45 A simple thermochemical analysis shows the pathway for CO2 reduction to CO to go through COOH* and CO* intermediates

Figure 1. Performance map for CO2RR over reported N-doped carbonbased materials as well as N-doped reduced graphene oxide (N-rGO) studied in this work (experimental details are in the Supporting Information). Au and Ag are reported for comparison. Representative CO (a) Faradaic efficiencies, (b) current density, and their corresponding overpotentials are plotted for each reported catalyst. Data for the plots are listed in the Supporting Information, Table S1.

densities at the same potential are displayed in Figure 1b. We also report an in-house made N-doped reduced graphene oxide (NrGO) as a representative CO2RR catalyst (details are in the Supporting Information (SI)). Among them, nitrogen-doped carbon nanotubes (NCNT) have shown great promise in the electrochemical reduction of CO2 to CO at very low overpotentials.29−33 Apart from CO, other products such as acetate,41 methane,30 and formic acid38,42 have also been reported on different N-doped carbon-based materials. The aforementioned experimental reports show that different carbon catalysts present diverse activities, and in most cases, the major product of the CO2 reduction is found to be CO. The interplay between N-doping and naturally occurring defects generates a wide range of possible active motifs, contributing to the overall performance of the material. Different types of Ndoping can be formed including pyridinic, pyrrolic, and graphitic nitrogen. N-doping in different types of defects provides a wide range of diverse graphitic nitrogen with different activities.35 Through a simple thermochemical analysis, Wu et al. have suggested pyridinic-N to have the highest activity for the reduction of CO2 to CO.33 In this work, we examine the activity of a wide range of different possible active sites in the carbonbased materials with and without N-doping in pristine and defective carbon structures for the electrochemical reduction of CO2. We use density functional theory calculations (DFT) in conjunction with the computational hydrogen electrode (CHE)44 model to determine the adsorption energies of key intermediates of CO2 reduction reaction in order to get the trends in activity and draw comparisons with currently known carbon-based materials (details are in the Supporting Information). The CHE model exploits that the chemical potential of a proton−electron pair is equal to that of gas-phase H2, at Uelec = 0.0 V vs the reversible hydrogen electrode (RHE).44 The maximum potential allowed by thermodynamics is obtained by shifting the chemical potential of the electrons at the equilibrium potential of U = −0.12 V. This allows us to obtain the theoretical overpotential, which is defined as the maximum difference

CO2 + * + (H+ + e−) → COOH*

(1)

COOH* + (H+ + e−) → CO* + H 2O

(2)

CO* → CO + *

(3)

Total: CO2 + 2(H+ + e−) → CO + H 2O E0 = −0.12 V

(4)

If the surface adsorbs CO weakly, the CO molecule readily desorbs from the surface, and the major product will be CO gas. Au and Ag, for example, bind CO weakly and are selective for the two-electron reduction of CO2 to CO.45 On the other hand, for the transition metal surfaces that bind CO, further reduced products past CO may be produced, where the activity has been suggested to be limited by the reduction of CO* to CHO*.1 Figure 2b displays the calculated chemisorption energies of COOH* vs CO* on the studied active motifs as well as transition metals (111) surfaces. As can be seen, there is a linear scaling relation between COOH* and CO* adsorption energies in the case of transition metals.45 The orange dashed line shows the chemical potential of CO in the gas phase, separating the region that CO can be bound to the surface and therefore gets further reduced from no adsorption regime. Copper is in the border and weakly binds CO, while gold and silver are in the no adsorption regime. We note that most of the examined active sites do not fall on the known scaling line for the transition metals because they do not bind CO. This is the reason that carbon-based materials represent especially interesting alternatives for the two-electron reduction of CO2 to CO. For the examined active sites in Figure 2a with positive CO* binding energy, the activity toward CO production can be determined exclusively by the COOH* binding strength. An optimal COOH* binding strength, not too weak and not too strong, is desirable in line with the Sabatier principle (Figure 2c). The equilibrium potential for two-electron reduction of CO2 to CO (eq 4) is −0.12 V vs RHE. Using thermodynamic 11081

DOI: 10.1021/acssuschemeng.7b03031 ACS Sustainable Chem. Eng. 2017, 5, 11080−11085

Research Article

ACS Sustainable Chemistry & Engineering

Figure 2. (a) Optimized structures of different examined carbon structures with and without N-doped, where LD stands for line defect. Here, 5, 6, 7, and 8 refer to the pentagon, hexagon, heptagon, and octagon in the defect configuration. Color code: C, gray; N, blue. (b) Calculated chemisorption energies of COOH* vs CO* intermediates for examined carbon-based materials. Numbers refer to the structures shown in panel (a). Transition metals (111) are depicted in black, adapted from ref 45 for comparison. Energies for CO* are relative to gas-phase CO and H2 and for COOH* are relative to gas-phase CO2 and H2. The vertical dashed line corresponds to the chemical potential of CO in the gas phase. (c) Activity volcano for two-electron reduction of CO2 to CO. (d) Free energy diagram of CO2 to CO conversion on different examined active sites displayed in panel (a) at equilibrium potential of −0.12 V vs RHE. (e) The p-projected density of states for pure graphene and N-doped defective carbon labeled as number (5) in panel (a) as an example. 11082

DOI: 10.1021/acssuschemeng.7b03031 ACS Sustainable Chem. Eng. 2017, 5, 11080−11085

Research Article

ACS Sustainable Chemistry & Engineering

obtain an estimate of the selectivity. This analysis does not apply for the catalysts with strong CO* binding energy, where the CO desorption process is slow and the Faradaic efficiency is negligible. Figure 3a and b display the UL(CO2), UL(H2), and (UL(CO2) − UL(H2)) for the configurations shown in the free energy

considerations alone, the ideal catalyst should have a COOH* binding energy of 0.12 eV at zero applied potential. The horizontal red line in Figure 2b displays the adsorption energy of the ideal catalyst that effectively catalyzes the twoelectron reduction of CO2 to CO. The vertical distance between this line and each point in the graph can be used as a metric for activity. Interestingly, many of the examined active sites show higher activity than Au for the two-electron reduction of CO2 to CO. This information can be translated this into a thermodynamic activity volcano (Figure 2c) with COOH binding energy as descriptor and a free energy diagram at the equilibrium potential of −0.12 V vs RHE (Figure 2d). These plots suggest that the COOH* formation is rather difficult and is the bottleneck step for CO production on pyrrolic-N (Figure2a (15)) and N-doped at pristine graphene (Figure 2a (1)). Therefore, the overall CO2 reduction activity is severely limited by formation of COOH*. Pyridinic-N (14), on the other hand, has a surmountable *COOH formation energy at relatively low overpotentials, and CO desorption is facile. This has been previously reported by Wu et al.,33 who assigned the CO2 reduction activity on N-doped carbon nanotube to the pyridinic-N based on a similar thermodynamic analysis. Figure 2b and c show that several defects with and without N-doped are as active as pyridinic-N for the CO2RR. For example, the 558-line defect (Figure2a (12)) and N-doped in the 555-6-777 (Figure2a (7)) defective structure show comparable activity to pyridinic-N. Moreover, there are a number of other sites, such as N-doped 5577 (Figure2a (5)), 555-6-777 (Figure2a (6)), N-doped 558-line defect (Figure2a (13)), and 555-777 defect (Figure2a (2)), that exhibit lower activity than pyridinic-N but represent activities close to or higher than Au and Ag. The difference in scaling behavior between transition metals and defective graphene can be rationalized through their differences in electronic structure. CO binds to transition metals through the interplay between attractive π* and repulsive σ interactions with the metal d-states.46 Pristine graphene is a semiconductor with no states at the Fermi level (Figure 2e),47 which does not provide attractive interactions with CO. Intense electron localization, however, is present around N-dopants and defects in graphene (Figure 2e).47 Such sites introduce the possibility for covalent bonds to COOH, which is a radical, but not CO, which has a closed shell configuration. This bonding picture is similar to that suggested for transition metal dichalcogenide edges48 and p-block dopants in metals,49 both of which show similar deviations from transition metal scaling relations. HER Analysis. In general, the carbon-based materials represent high overpotentials for a hydrogen evolution reaction (HER) compared to metal surfaces.25 Herein, we analyze the competition between HER and CO2RR using the thermodynamic limiting potentials UL. It was shown that the difference between the thermodynamic limiting potentials for CO2RR and HER, UL(CO2) − UL(H2), follows trends in selectivity on transition metals, and the magnitude of this difference for the metals is less than 0.5 V on the metals that reduce CO2 with significant Faradaic efficiency (FE).45 For metals binding CO* weaker than Cu, UL(CO2) is determined by the COOH* formation from CO2. For Cu and metals binding CO* stronger than Cu, protonation of CO* to CHO* determines the limiting potential. Assuming that the corresponding electrochemical barriers scale similarly in both carbon-based and transition metal catalysts, we can apply the same analysis in the present work to

Figure 3. (a) Calculated limiting potentials for HER and CO2RR. (b) Calculated UL(CO2) − UL(H2) on the active sites studied in Figure 2d compared to Au(111).

diagram (Figure 2d) with respect to gold. The more positive UL(CO2) − UL(H2) corresponds to higher selectivity toward CO2 reduction. This analysis therefore suggests that of the configurations with surmountable COOH* formation energies shown in the Figure 2(d) configurations (7), (12), and (14) should be selective for CO2RR over HER. Pyridinic-N (configuration (14) has the most positive UL(CO2) − UL(H2) and therefore is suggested to be the most selective site for the CO2RR over HER.



CONCLUSION In summary, we have studied a diverse range of possible active sites in N-doped carbon-based materials and showed that the calculated adsorption energies do not follow the scaling relations established for transition metal surfaces. This difference in scaling is rationalized by the presence of electron localization in the defected sites of graphene which binds COOH but not CO. This unique scaling rationalizes the numerous experimental reports of low overpotentials for CO2 reduction to CO on carbon-based materials. We used the calculated COOH* adsorption energy as an activity metric to map out the trends in activity of different active site motifs. COOH* adsorption energies vary significantly depending on the type of active site, 11083

DOI: 10.1021/acssuschemeng.7b03031 ACS Sustainable Chem. Eng. 2017, 5, 11080−11085

Research Article

ACS Sustainable Chemistry & Engineering

Nanoscale Copper Catalyst. J. Am. Chem. Soc. 2014, 136 (38), 13319− 13325. (4) Kuhl, K. P.; Cave, E. R.; Abram, D. N.; Jaramillo, T. F. New Insights into the Electrochemical Reduction of Carbon Dioxide on Metallic Copper Surfaces. Energy Environ. Sci. 2012, 5, 7050. (5) Karamad, M.; Hansen, H. A.; Rossmeisl, J.; Nørskov, J. K. Mechanistic Pathway in the Electrochemical Reduction of CO 2 on RuO 2. ACS Catal. 2015, 5 (7), 4075−4081. (6) Kas, R.; Kortlever, R.; Milbrat, A.; Koper, M. T. M.; Mul, G.; Baltrusaitis, J.; D’Alessandro, D. M.; Smit, B.; Long, J. R.; Ritter, S. K.; et al. Electrochemical CO2 Reduction on Cu2O-Derived Copper Nanoparticles: Controlling the Catalytic Selectivity of Hydrocarbons. Phys. Chem. Chem. Phys. 2014, 16 (24), 12194. (7) Le, M.; Ren, M.; Zhang, Z.; Sprunger, P. T.; Kurtz, R. L.; Flake, J. C. Electrochemical Reduction of CO2 to CH3OH at Copper Oxide Surfaces. J. Electrochem. Soc. 2011, 158 (5), E45. (8) Frese, K. W. J. Electrochem. Soc. 1991, 138 (11), 3338−3344. (9) Liu, L.-M.; Zhang, C.; Thornton, G.; Michaelides, A. Structure and Dynamics of Liquid Water on Rutile TiO_{2}(110). Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 82 (16), 161415. (10) Sen, S.; Liu, D.; Palmore, G. T. R. Electrochemical Reduction of CO 2 at Copper Nanofoams. ACS Catal. 2014, 4, 3091−3095. (11) Chang, T. Y.; Liang, R. M.; Wu, P. W.; Chen, J. Y.; Hsieh, Y. C. Electrochemical Reduction of CO2 by Cu2O-Catalyzed Carbon Clothes. Mater. Lett. 2009, 63 (12), 1001−1003. (12) Asadi, M.; Kumar, B.; Behranginia, A.; Rosen, B. A.; Baskin, A.; Repnin, N.; Pisasale, D.; Phillips, P.; Zhu, W.; Haasch, R.; et al. Robust Carbon Dioxide Reduction on Molybdenum Disulphide Edges. Nat. Commun. 2014, 5, 4470. (13) Asadi, M.; Kim, K.; Liu, C.; Addepalli, A. V.; Abbasi, P.; Yasaei, P.; Phillips, P.; Behranginia, A.; Cerrato, J. M.; Haasch, R.; et al. Nanostructured Transition Metal Dichalcogenide Electrocatalysts for CO2 Reduction in Ionic Liquid. Science 2016, 353 (6298), 467−470. (14) Tripkovic, V.; Vanin, M.; Karamad, M.; Björketun, M. E.; Jacobsen, K. W.; Thygesen, K. S.; Rossmeisl, J. Electrochemical CO 2 and CO Reduction on Metal-Functionalized Porphyrin-like Graphene. J. Phys. Chem. C 2013, 117 (18), 9187−9195. (15) Shen, J.; Kortlever, R.; Kas, R.; Birdja, Y. Y.; Diaz-Morales, O.; Kwon, Y.; Ledezma-Yanez, I.; Schouten, K. J. P.; Mul, G.; Koper, M. T. M. Electrocatalytic Reduction of Carbon Dioxide to Carbon Monoxide and Methane at an Immobilized Cobalt Protoporphyrin. Nat. Commun. 2015, 6, 8177. (16) Varela, A. S.; Ranjbar Sahraie, N.; Steinberg, J.; Ju, W.; Oh, H.-S.; Strasser, P. Metal-Doped Nitrogenated Carbon as an Efficient Catalyst for Direct CO 2 Electroreduction to CO and Hydrocarbons. Angew. Chem. 2015, 127 (37), 10908−10912. (17) Cheng, M.-J.; Kwon, Y.; Head-Gordon, M.; Bell, A. T. Tailoring Metal-Porphyrin-Like Active Sites on Graphene to Improve the Efficiency and Selectivity of Electrochemical CO 2 Reduction. J. Phys. Chem. C 2015, 119 (37), 21345−21352. (18) Lin, S.; Diercks, C. S.; Zhang, Y.-B.; Kornienko, N.; Nichols, E. M.; Zhao, Y.; Paris, A. R.; Kim, D.; Yang, P.; Yaghi, O. M.; et al. Covalent Organic Frameworks Comprising Cobalt Porphyrins for Catalytic CO2 Reduction in Water. Science 2015, 349 (6253), 1208−1213. (19) Hod, I.; Sampson, M. D.; Deria, P.; Kubiak, C. P.; Farha, O. K.; Hupp, J. T. Fe-Porphyrin-Based Metal−Organic Framework Films as High-Surface Concentration. ACS Catal. 2015, 5, 6302. (20) Dai, L. Functionalization of Graphene for Efficient Energy Conversion and Storage. Acc. Chem. Res. 2013, 46 (1), 31−42. (21) Wang, H.; Yuan, X.; Zeng, G.; Wu, Y.; Liu, Y.; Jiang, Q.; Gu, S. Three Dimensional Graphene Based Materials: Synthesis and Applications from Energy Storage and Conversion to Electrochemical Sensor and Environmental Remediation. Adv. Colloid Interface Sci. 2015, 221, 41−59. (22) Sun, Y.; Wu, Q.; Shi, G. Graphene Based New Energy Materials. Energy Environ. Sci. 2011, 4 (4), 1113. (23) Brownson, D. A. C.; Kampouris, D. K.; Banks, C. E. An Overview of Graphene in Energy Production and Storage Applications. J. Power Sources 2011, 196 (11), 4873−4885.

allowing for a wide range of possible activities. Based on this simple thermodynamic picture, several configurations are suggested to be selective for CO2RR over HER, and pyridinicN is suggested to be the most selective site for the CO2RR over HER. Future work will investigate activation energies and kinetics in order to fully evaluate the efficiency and selectivity of candidate active site motifs.



COMPUTATIONAL DETAILS



ASSOCIATED CONTENT

Density functional theory calculations for the carbon-based materials were performed using the QUANTUM ESPRESSO program package50 and the Atomic Simulation Environment (ASE).51 The electronic wave functions were expanded in plane waves up to a cutoff energy of 500 eV, while the electron density was represented on a grid with an energy cutoff of 5000 eV. Core electrons are approximated with ultrasoft pseudopotentials.52 We use the BEEF-vdW exchange−correlation functional,53 which has been shown to accurately describe chemisorption as well as physisorption properties on graphene. For adsorption studies, a supercell of lateral size of (4 × 4) for nondefective and (8 × 8) for defective structures are used, and a vacuum region of at least 18 Å is used to decouple the periodic images. The first Brillouin zones are sampled with (4 × 4 × 1) and (2 × 2 × 1) Monkhorst−Pack54 k-points for the (4 × 4) and (8 × 8) supercells, respectively. All the atoms were allowed to relax until the maximum force on each atom was below 0.05 eV/A for each calculation. S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acssuschemeng.7b03031. Information as mentioned in the text. (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Samira Siahrostami: 0000-0002-1192-4634 Kun Jiang: 0000-0003-3148-5058 Karen Chan: 0000-0002-6897-1108 Haotian Wang: 0000-0002-3552-8978 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This material is based upon work performed by the Joint Center for Artificial Photosynthesis, a DOE Energy Innovation Hub, as follows: Density functional theory calculations was supported through the Office of Science, U.S. Department of Energy, under Award No. DE-SC0004993. Analysis of experimental data was supported by the Rowland Junior Fellows Program at Rowland Institute, Harvard University.



REFERENCES

(1) Peterson, A. A.; Abild-Pedersen, F.; Studt, F.; Rossmeisl, J.; Nørskov, J. K. How Copper Catalyzes the Electroreduction of Carbon Dioxide into Hydrocarbon Fuels. Energy Environ. Sci. 2010, 3 (9), 1311. (2) Schouten, K. J. P.; Kwon, Y.; van der Ham, C. J. M.; Qin, Z.; Koper, M. T. M. A New Mechanism for the Selectivity to C1 and C2 Species in the Electrochemical Reduction of Carbon Dioxide on Copper Electrodes. Chem. Sci. 2011, 2 (10), 1902. (3) Manthiram, K.; Beberwyck, B. J.; Alivisatos, A. P. Enhanced Electrochemical Methanation of Carbon Dioxide with a Dispersible 11084

DOI: 10.1021/acssuschemeng.7b03031 ACS Sustainable Chem. Eng. 2017, 5, 11080−11085

Research Article

ACS Sustainable Chemistry & Engineering (24) Pumera, M. Graphene-Based Nanomaterials for Energy Storage. Energy Environ. Sci. 2011, 4 (3), 668−674. (25) Yang, N.; Waldvogel, S. R.; Jiang, X. Electrochemistry of Carbon Dioxide on Carbon Electrodes. ACS Appl. Mater. Interfaces 2016, 8, 28357. (26) Liu, J.; Song, P.; Ning, Z.; Xu, W. Recent Advances in Heteroatom-Doped Metal-Free Electrocatalysts for Highly Efficient Oxygen Reduction Reaction. Electrocatalysis 2015, 6 (2), 132−147. (27) Wong, W. Y.; Daud, W. R. W.; Mohamad, A. B.; Kadhum, A. A. H.; Loh, K. S.; Majlan, E. H. Recent Progress in Nitrogen-Doped Carbon and Its Composites as Electrocatalysts for Fuel Cell Applications. Int. J. Hydrogen Energy 2013, 38 (22), 9370−9386. (28) Terrones, H.; Lv, R.; Terrones, M.; Dresselhaus, M. S. The Role of Defects and Doping in 2D Graphene Sheets and 1D Nanoribbons. Rep. Prog. Phys. 2012, 75 (6), 062501. (29) Sharma, P. P.; Wu, J.; Yadav, R. M.; Liu, M.; Wright, C. J.; Tiwary, C. S.; Yakobson, B. I.; Lou, J.; Ajayan, P. M.; Zhou, X.-D. NitrogenDoped Carbon Nanotube Arrays for High-Efficiency Electrochemical Reduction of CO2: On the Understanding of Defects, Defect Density, and Selectivity. Angew. Chem., Int. Ed. 2015, 54 (46), 13701−13705. (30) Li, W.; Seredych, M.; Rodríguez-Castellón, E.; Bandosz, T. J. Metal-Free Nanoporous Carbon as a Catalyst for Electrochemical Reduction of CO2 to CO and CH4. ChemSusChem 2016, 9 (6), 606− 616. (31) Chai, G.-L.; Guo, Z.-X. Highly Effective Sites and Selectivity of Nitrogen-Doped graphene/CNT Catalysts for CO 2 Electrochemical Reduction. Chem. Sci. 2016, 7 (2), 1268−1275. (32) Liu, Y.; Chen, S.; Quan, X.; Yu, H. Efficient Electrochemical Reduction of Carbon Dioxide to Acetate on Nitrogen-Doped Nanodiamond. J. Am. Chem. Soc. 2015, 137 (36), 11631−11636. (33) Wu, J.; Yadav, R. M.; Liu, M.; Sharma, P. P.; Tiwary, C. S.; Ma, L.; Zou, X.; Zhou, X.-D.; Yakobson, B. I.; Lou, J.; et al. Achieving Highly Efficient, Selective, and Stable CO2 Reduction on Nitrogen-Doped Carbon Nanotubes. ACS Nano 2015, 9 (5), 5364−5371. (34) Jiang, Y.; Yang, L.; Sun, T.; Zhao, J.; Lyu, Z.; Zhuo, O.; Wang, X.; Wu, Q.; Ma, J.; Hu, Z. Significant Contribution of Intrinsic Carbon Defects to Oxygen Reduction Activity. ACS Catal. 2015, 5 (11), 6707− 6712. (35) To, J. W. F.; Ng, J. W. D.; Siahrostami, S.; Koh, A. L.; Lee, Y.; Chen, Z.; Fong, K. D.; Chen, S.; He, J.; Bae, W.-G.; et al. HighPerformance Oxygen Reduction and Evolution Carbon Catalysis: From Mechanistic Studies to Device Integration. Nano Res. 2017, 10, 1163. (36) Hou, Z.; Wang, X.; Ikeda, T.; Terakura, K.; Oshima, M.; Kakimoto, M.; Miyata, S. Interplay between Nitrogen Dopants and Native Point Defects in Graphene. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85 (16), 165439. (37) Jhong, H.-R. M.; Tornow, C. E.; Smid, B.; Gewirth, A. A.; Lyth, S. M.; Kenis, P. J. A. A Nitrogen-Doped Carbon Catalyst for Electrochemical CO 2 Conversion to CO with High Selectivity and Current Density. ChemSusChem 2017, 10 (6), 1094−1099. (38) Zhang, S.; Kang, P.; Ubnoske, S.; Brennaman, M. K.; Song, N.; House, R. L.; Glass, J. T.; Meyer, T. J. Polyethylenimine-Enhanced Electrocatalytic Reduction of CO 2 to Formate at Nitrogen-Doped Carbon Nanomaterials. J. Am. Chem. Soc. 2014, 136 (22), 7845−7848. (39) Wu, J.; Yadav, R. M.; Liu, M.; Sharma, P. P.; Tiwary, C. S.; Ma, L.; Zou, X.; Zhou, X.-D.; Yakobson, B. I.; Lou, J.; et al. Achieving Highly Efficient, Selective, and Stable CO 2 Reduction on Nitrogen-Doped Carbon Nanotubes. ACS Nano 2015, 9 (5), 5364−5371. (40) Rosen, J.; Hutchings, G. S.; Lu, Q.; Rivera, S.; Zhou, Y.; Vlachos, D. G.; Jiao, F. Mechanistic Insights into the Electrochemical Reduction of CO 2 to CO on Nanostructured Ag Surfaces. ACS Catal. 2015, 5 (7), 4293−4299. (41) Liu, Y.; Chen, S.; Quan, X.; Yu, H. Efficient Electrochemical Reduction of Carbon Dioxide to Acetate on Nitrogen-Doped Nanodiamond. J. Am. Chem. Soc. 2015, 137 (36), 11631−11636. (42) Wang, H.; Jia, J.; Song, P.; Wang, Q.; Li, D.; Min, S.; Qian, C.; Wang, L.; Li, Y. F.; Ma, C.; et al. Efficient Electrocatalytic Reduction of CO2 by Nitrogen-Doped Nanoporous Carbon/Carbon Nanotube

Membranes - A Step Towards the Electrochemical CO2 Refinery. Angew. Chem. 2017, 129, 7955. (43) Zhu, W.; Michalsky, R.; Metin, Ö .; Lv, H.; Guo, S.; Wright, C. J.; Sun, X.; Peterson, A. A.; Sun, S. Monodisperse Au Nanoparticles for Selective Electrocatalytic Reduction of CO2 to CO. J. Am. Chem. Soc. 2013, 135 (45), 16833−16836. (44) Nørskov, J. K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J. R.; Bligaard, T.; Jónsson, H. Origin of the Overpotential for Oxygen Reduction at a Fuel-Cell Cathode. J. Phys. Chem. B 2004, 108 (46), 17886−17892. (45) Shi, C.; Hansen, H. A.; Lausche, A. C.; Nørskov, J. K. Trends in Electrochemical CO2 Reduction Activity for Open and Close-Packed Metal Surfaces. Phys. Chem. Chem. Phys. 2014, 16 (10), 4720−4727. (46) Nilsson, A.; Pettersson, L. G. M. Chemical Bonding on Surfaces Probed by X-Ray Emission Spectroscopy and Density Functional Theory. Surf. Sci. Rep. 2004, 55, 49−167. (47) Jing, Y.; Zhou, Z. Computational Insights into Oxygen Reduction Reaction and Initial Li 2 O 2 Nucleation on Pristine and N-Doped Graphene in Li−O 2 Batteries. ACS Catal. 2015, 5 (7), 4309−4317. (48) Hong, X.; Chan, K.; Tsai, C.; Nørskov, J. K. How Doped MoS 2 Breaks Transition-Metal Scaling Relations for CO 2 Electrochemical Reduction. ACS Catal. 2016, 6 (7), 4428−4437. (49) Lim, H.-K.; Shin, H.; Goddard, W. A.; Hwang, Y. J.; Min, B. K.; Kim, H. Embedding Covalency into Metal Catalysts for Efficient Electrochemical Conversion of CO 2. J. Am. Chem. Soc. 2014, 136 (32), 11355−11361. (50) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I.; et al. QUANTUM ESPRESSO: A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys.: Condens. Matter 2009, 21 (39), 395502. (51) Atomic Simulation Environment (ASE). Center for Atomic Scale Material Design (CAMD), Technical University of Denmark, Lyngby. https://wiki.fysik.dtu.dk/ase (accessed October 2017). (52) Adllan, A. A.; Dal Corso, A. Ultrasoft Pseudopotentials and Projector Augmented-Wave Data Sets: Application to Diatomic Molecules. J. Phys.: Condens. Matter 2011, 23 (42), 425501. (53) Wellendorff, J.; Lundgaard, K. T.; Møgelhøj, A.; Petzold, V.; Landis, D. D.; Nørskov, J. K.; Bligaard, T.; Jacobsen, K. W. Density Functionals for Surface Science: Exchange-Correlation Model Development with Bayesian Error Estimation. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85 (23), 235149. (54) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13 (12), 5188−5192.

11085

DOI: 10.1021/acssuschemeng.7b03031 ACS Sustainable Chem. Eng. 2017, 5, 11080−11085