Design of Efficient Catalysts with Double Transition Metal Atoms on

Apr 19, 2016 - Mahmood , J.; Lee , E. K.; Jung , M.; Shin , D.; Jeon , I. P.; Jung , S. M.; Choi , H. J.; Seo , J. M.; Bae , S. Y.; Sohn , S. D. Nitro...
1 downloads 17 Views 2MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Letter 2

Design of Efficient Catalysts with Double Transition Metal Atoms on CN Layer Xiyu Li, Wenhui Zhong, Peng Cui, Jun Li, and Jun Jiang J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b00096 • Publication Date (Web): 19 Apr 2016 Downloaded from http://pubs.acs.org on April 22, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Design of Efficient Catalysts with Double Transition Metal Atoms on C2N Layer Xiyu Li,1,† Wenhui Zhong, 2,† Peng Cui,1 Jun Li,1 Jun Jiang1,* 1

Hefei National Laboratory for Physical Sciences at the Microscale, Collaborative Innovation Center of Chemistry for Energy Materials, Hefei Science Center of CAS, School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, Anhui 230026, P. R. China 2 Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Institute of Applied Physics, Guizhou Synergetic Innovation Cen-ter of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Normal College, Gaoxin Road 115, Guiyang, Guizhou 550018, P. R. China Supporting Information Placeholder

ABSTRACT: Heterogeneous catalysis often involves molecular adsorptions to charged catalyst site, and reactions triggered by catalyst charges. Here we use first-principles simulations to design oxygen reduction reaction (ORR) catalyst based on double transition metal (TM) atoms stably supported by two-dimensional crystal C2N. It not only holds characters of low cost and high durability, but also effectively accumulates surface polarization charges on TMs and later deliveries to adsorbed O2 molecule. The Co-Co, Ni-Ni, and Cu-Cu catalysts, exhibit high adsorption energies and extremely low dissociation barriers for O2, as compared with their single-atom counterparts. Co-Co on C2N presents less than the half value of the reaction barrier of bulk Pt catalysts in the ORR rate-determining steps. These catalytic improvements are well explained by the dependences of charge polarization on various systems, which opens up a new strategy for optimizing TM catalytic performance with the least metal atoms on porous low-dimensional materials.

TOC Graphic

Heterogeneous catalysis utilizing the ensemble of a small number of atoms on certain substrate often exhibits excellent catalytic performance for many important chemical reactions including oxygen reduction reaction, carbon monoxide oxidation, and hydrogen production.1-4 It provides a promising way to minimize the use of noble metals or even substitute them with a few non-noble transition metal (TM) atoms,5,6 so as to meet the ultimate goal of cost-effective catalysis. Recent years have witnessed remarkable catalyst designs of several TM atoms.5,7 However, the search for appropriate support substrate remains

very challenging, and the control of optimal TM atom number or size is extremely hard, both of which hindered the realistic applications. Metal substrates have been widely used to support small TM nano-clusters8,9 to offer highly active reaction sites, but their applications are often limited by the high cost of bulk metal. Nonmetallic supports can utilize carbon-based materials such as C3N4, while electrochemical reactions would be jeopardized by the low efficiency in charge delivery.10 Multi-transition metal active sites stabilized by organic and/or inorganic ligands have shown to be very promising11-14, but their surface-to-volume ratios are not always ideal and the charge injection and transportation in ligands are not very efficient. Graphene owns high conductivity, but TM atoms often drift easily on graphene surface to form aggregated particles, which heavily reduces active catalytic sites. Strategies have been proposed to fix several TM atoms on graphene through strong chemical interactions such as replacing C atom with TM, which however request extra endeavor and sophisticated techniques to decorate the perfect graphene.15 Meanwhile, the optimal numbers of TM atoms for superior performance remain unknown for most cases, prohibiting the rational design or optimization of catalysts. The motif of a single TM atom was widely proposed,4,5 but it is also stated that several atoms being together is necessary to ensure high catalysis performance.2,3 Moreover, the difficulties to identify or control the actual number of TM atoms impede a thorough examination of the size dependence of catalysis performance from the experimental point of view. Very recently, Baek. et. al. reported a novel two-dimensional (2D) layered crystal with uniform hole vacancies named C2N holey 2D (C2N-h2D), easily prepared through wet chemical reaction.16 As seen in Figure 1a, its unique porous structure originated from 2D graphene layer provides ideal support for metallic active center, and the high surface-to-volume ratio property ensures sufficient exposure of TM atoms to interact with reactant molecules.17 Unlike C3N4 layer which breaks a graphene into fragments connected by non-conductive nitrogen nodes, the C2N layer holds continuous network structure so as to maintain high electric conductivity.16 More importantly, it contains wellordered big vacancies to anchor one or two TM atoms tightly (Figure 1b and c). The sp2-bonded nitrogen atoms at the vacancy hole edge can provide the coordination sites for TM atom owing to their lone-pair electrons, which would couple strongly with TM and thereby prohibit TM drifting and aggregating. Spectroscopic tools such as x-ray or Raman are able to measure the TM-TM bond, making it feasible to identify (and later to control) one and

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

double TM configurations on C2N. Taking advantages of high durability and conductivity, the systems of C2N encapsulating one or double TM atoms in the vacancy (namely TM@C2N or TM2@C2N) offer opportunities to design efficient catalysts by optimizing the number of TM atoms. Here, we applied first-principles theoretical simulations to study the system of four types of TM atoms (Pt, Co, Ni and Cu) anchored on a C2N monolayer, based on which the catalytic activity of oxygen were examined. Promoting O2 activation is of fundamental importance for chemical reactions in many biologic applications and electrochemical field, such as CO oxidation18 and ORR19. Tremendous efforts have been devoted to exploring durable, efficient and low-cost catalysts to accelerate the sluggish kinetics of ORR in fuel cell, which critically restricts electrochemical performance. Pt-based catalyst probably exhibits the best performance in acidic and alkaline electrolyte, while its utilization is limited by the scarcity and poor durability. In our designed TM@C2N and TM2@C2N system, merits of low cost and morphological variety can be expected. The strong couplings of several TM atoms with sp2-bonded N and later with reactant molecule would cause d electron redistribution, so as to promote the activities of non-noble Co, Ni and Cu based systems to be better than (or at least as good as) those expensive Pt-based catalysts.

Figure 1. Top and side view of the structure of the monolayer C2N (h2D form) in a 2×2 supercell (a), and the structures of TMN-hole in TM@C2N (b) and TM2@C2N (c). Our simulations were carried at the spin-polarized density functional theory (DFT) level using the VASP package (See details in Support Information).20 The Perdew, Burke, and Ernzerhof (PBE) exchange-correlation functional within a generalized gradient approximation (GGA)13 and the projector augmented-wave (PAW)21 potential were employed. The GGA+U method was applied to describe partially filled d-orbitals by considering coulomb and exchange corrections.22 The method of climbing image nudged elastic band (CI-NEB) was used for transition state search.23 The optimized structure of a monolayer C2N-h2D holds a unit cell of 12 C and 6 N atoms, where four unit cells form a uniform hole (Figure 1a). Our calculated lattice parameter of 8.32 Å matches with experiment.16 The computed energy band structure (Figure S1) agrees with previous reports.24 The deposition of one TM atom of Pt, Co, Ni, Cu (TM@C2N in Figure 1b) causes negligible structural change, while double TM atoms anchoring slightly distorts the 2D plane (TM2@C2N in Figure 1c). Each TM atom bonds strongly to two N atoms with bonding lengths of 1.93 ~ 2.13 Å (Figure S2, S3). The binding energy (Eb) was calculated to examine the structural stability, where Eb = ETMn + EC2N ETMn@C2N, and ETMn@C2N, ETMn, and EC2N denote energies of the hybrid structure and isolated parts. In Table S1, the smallest Eb of TM on C2N is 3.85 eV, suggesting strong couplings and good stability. More importantly, the addition of another TM on TM@C2N to form TM2@C2N effectively reduces system energy (increases Eb) for 1.91 ~ 5.09 eV. Since the computed binding energies of the single or double TM on C2N are slightly larger than the corresponding TM bulk cohesive energies (Table S1),

Page 2 of 9

one can expect good stability for the anchoring of TM on C2N.25 Meanwhile, the computed diffusion barriers for TM atoms to move from one hole vacancy to the next hole site were 2.97 ~ 3.91 eV (Table S2), indicating strong chemical bondings between TM and C2N. This prohibits the drifting of TM atoms to aggregate. In contrast, TM atoms on other 2D layered structures such as graphene and h-BN have very small diffusion barriers (0.02~0.69 eV in Table S2), causing the serious problem of TM aggregation. Meanwhile, since TM atoms can drift on the non-vacancy surface part of C2N before being trapped by vacancies, one can adjust the existence of TM@C2N and TM2@C2N configurations by changing TM deposition density, and later identify them with TM-TM bond measurement. Meanwhile, Bader charge analysis reveals that the C2N part can extract 0.47 ~ 0.96 e- electrons from TM, which are increased to 0.64 ~ 1.50 e- in TM2@C2N. These accumulate polarization charges in C2N, which partially occupied the previous conduction bands (Figure S4) and induced a metallic feature favoring electrochemical reactions. Here the magnetic couplings of TMs are weak (Table S3). We then move to investigate the ORR catalytic properties. It is well-established that the initial adsorption manner of O2 molecule often influence or determine the entire reaction pathway. For the pure monolayer C2N, the physical adsorption of O2 (Eb = 0.06 eV) and nearly zero polarization charge suggest very low activity of O2 on C2N (Figure S5). While for the case of TM@C2N, there are strong positive charges remaining on TM sites, favoring O2 adsorption. Starting with two configurations of O-O bond being parallel or vertical to the C2N plane, we found that O2 is always adsorbed to TM sites. The most stable configurations were achieved by maximizing the adsorption energy (Eads = ETMn@C2N + EO2 - ETMn@C2N-O2), where ETMn@C2N-O2, ETMn@C2N, and EO2 represent energies of the hybrid structure and separated parts. In Figure 2a, O2 is adsorbed on Pt@C2N with a side-on configuration, in which both O atoms connect to Pt with equal distances. In contrast, the adsorption of O2 to Ni/Cu@C2N takes the end-on configuration, where only one O atom bonds to TM and the other one is far from metal. O2 on Co@C2N also chooses the end-on mode except the disconnected O is still close to Co. As a result, the order of O2 adsorption energy is Pt > Co > Ni > Cu (Table 1). Nevertheless, the high values of Eads (0.64 ~ 1.22 eV) confirm strong chemical adsorption. The adsorption also elongates the O-O bond length from 1.22 Å in free O2 to 1.28 ~ 1.36 Å in TM@C2N (Table 1). Electronic structures reflected by partial density of states (PDOS) in Figure 2b, show that electronic coupling mainly exists between TM-3d and O2-2p orbitals. Importantly, TM@C2N always donates electrons to the adsorbed O2, accumulating negative charges (yellow bubbles in Figure 2c) which normally help trigger reduction reaction.

Figure 2. (a) The optimized configurations of O2 adsorbed on TM@C2N (TM= Pt, Co, Ni, Cu). (b) The partial density of states

ACS Paragon Plus Environment

Page 3 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(PDOS) of O2 on TM@C2N. (c)The distribution of exchange/polarized charges between O2 and TM@C2N with the isovalue of 0.003 e Å-3 (computed as the difference charge of systems before and after O2 adsorption). Cyan and yellow bubbles represent positive and negative charges, respectively. The adsorption of O2 would then be followed by its dissociation, which is often the rate determining step for the ORR processes.26 With negative charges in the O2 on Pt@C2N, the O-O bond scission would occur via the transition state (TS) in Figure S6, where the dissociating O-O bond is stretched from 1.36 Å to 2.10 Å. In the TS geometry, a newly formed O atom via the O-O bond scission moves to an atop position of Pt. In the final state (FS) of reaction, two O atoms locate at the atop position via Pt-O bonds. The O-O distance is extended to 2.78 Å in FS, while the Pt-O bonds are decreased from 2.00 Å in IS to 1.80 Å in FS. This process is exothermic by 0.84 eV, and needs to overcome a barrier of 1.29 eV (Figure S6). Substituting the Pt atom by Co, Ni and Cu, we identified similar reaction processes and obtained the O2 dissociation barriers of 1.56 eV, 2.36 eV and 2.79 eV (Figure S6), respectively. Unfortunately, these barriers are still too high to enable efficient ORR. Table 1. The computed adsorption energies (Eads), O-O bond lengths, and charges donated from TM@C2N and TM2@C2N catalysts to O2 obtained by Bader charge analysis. Eads (eV)

Pt

O-O bond (Å)

O2 charge (e-)

TM@ C 2N

TM2 @C2N

TM@ C 2N

TM2 @C2N

TM@ C 2N

TM2 @C2N

1.22

1.77

1.36

1.40

-0.51

-0.62

Co

1.18

2.39

1.33

1.47

-0.51

-0.81

Ni

0.81

1.83

1.29

1.44

-0.37

-0.76

Cu

0.64

1.41

1.28

1.45

-0.33

-0.63

We thus turn to the TM2@C2N systems with double TM atoms anchoring. In Figure 3a, O2 on TM2@C2N all take the side-on adsorption configurations where both O atoms bond to two TM atoms and pull them out the C2N plane. Such strong bindings result in high adsorption energies (1.41 ~ 2.39 eV in Table 1), and cause strong couplings between TM-3d and O2-2p orbitals as reflected by PDOS in Figure 3b. These enable O2 to extract more polarized charges (0.62 ~ 0.81 e- in Table 1) from TM2 and C2N (Figure 3c), in comparing with TM@C2N (Figure 2c). The O-O bond length (1.22 Å in free O2) on TM2@C2N is enlarged to 1.40 ~ 1.47 Å (Table 1), implying high activity. Taking O2 on Pt2@C2N as an example (Figure S7), the O-O bond (1.40 Å) scission would occur via a transition state of 1.43 Å O-O bond, and arrive at the final state of 2.47 Å. The reaction is exothermic by 1.75 eV, and only needs to overcome a barrier of 0.63 eV. Amazingly, the dissociation of O2 on Co2@C2N is exothermic by 1.25 eV with no energy barrier at all (Figure S7). This could be ascribed to the increase of d-electrons in Co-Co which in turn shift up the d-band center position of the single Co on C2N. The barriers of O-O bond scission on the Ni2@C2N and Cu2@C2N are also lowered to 0.11 eV and 0.56 eV, respectively (Figure S7). Therefore, the system of double TMs of Co-Co, Ni-Ni, Cu-Cu anchoring together on C2N is facile to dissociate O2.

Figure 3. (a) The optimized configurations of O2 adsorbed on TM2@C2N (TM= Pt, Co, Ni, Cu). (b) The partial density of states (PDOS) of O2 on TM2@C2N. (c) The distribution of exchange/polarized charges between O2 and TM2@C2N with the isovalue of 0.003 e Å-3 (computed as the difference charge of systems before and after O2 adsorption). Cyan and yellow bubbles represent positive and negative charges, respectively. According to above studies, the O2 activation with double TM atoms are evidently superior to those of the single TM atom, in adsorbing and dissociating O2. It was reported that O2 dissociation by single Pt on N-doped graphene27, MoS228, and single Co/Fe atom embedded in nitrogen-modified graphene29 require activation energy of 2.39 eV, 1.44 eV, and 1.96/1.19 eV, respectively. More importantly, with only a few metal utilizations, our TM2@C2N systems could dissociate O2 with energy barriers as low as those of the corresponding bulk metals.30,31 O2 molecule is the electron-accepting agent in ORR. Those performance improvements could thus be ascribed to the change of charge polarization, as revealed by our recent work on metal32 or semiconductor33 catalysts. It is found that metals can accumulate more positive charges in the double TM atoms (0.64 ~ 1.50 e+) than the single TM (0.47 ~ 0.96 e+), and O2 extracts more electrons from TM2@C2N (0.62 ~ 0.81 e-) than TM@C2N (0.33 ~ 0.51 e-). Putting all systems together, Figure 4a shows that except for the single Pt and double Pt-Pt cases (Pt is known for being excellent in adsorbing O2), the O2 adsorption ability increases with the increase of TM positive polarization charges following the order of Cu < Ni < Co < Cu-Cu < Ni-Ni < Co-Co. More strikingly, the O2 dissociation barrier decreases almost linear with the increase of negative polarization charge accumulated in O2 in Figure 4b. Therefore, it is convenient to predict the catalytic performance of several TM atoms by evaluating the induced polarization charges.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 9

amazingly low O2 dissociation barrier. The C2N material serves as an ideal substrate to trap single or double TM atoms and prohibit TM aggregation, offering high durable structures to adsorb molecules and accommodate charges. C2N also holds good electric conductivity to deliver charges for electrochemical and photocatalysis reactions. The striking improvement of catalytic O2 activity by Co-Co, Ni-Ni, and Cu-Cu catalysts is explained by the accumulation of polarization charges, which highlights the importance of atom number control in designing efficient and economic metal catalysts. With the most stable and inexpensive Co2@C2N structure, we demonstrated that the rate-determining barriers along two full ORR pathways are only 0.39 eV, which is less than the half value of the reported bulk Pt catalysts. Future work with better control over metal size by porous lowdimensional materials would further boost catalytic activities.

ASSOCIATED CONTENT Figure 4. (a) The dependence of O2 adsorption energy on TM polarization charge (positive) in TM@C2N and TM2@C2N. (b) The dependence of O2 dissociation barrier on its polarization charge (negative) induced by adsorption on TM@C2N and TM2@C2N. (c) The O2 and OOH dissociation initiated pathways for ORR catalyzed by Co2@C2N. The reaction energy (∆E) and activation energy (Ea) for all steps are given in parentheses with the form of (∆E, Ea).

AUTHOR INFORMATION Corresponding Author *E-mail:[email protected] Notes The authors declare no competing financial interest.

Author Contributions The effect of external electric field induced by electrochemistry or photocatalysis was modeled by adding one extra electron, in which ~0.25 e- electrons are trapped by O2 (Table S4) to facilitate reduction. The accumulation of polarization charges explains the order of O2 dissociation barrier (Figure S8), suggesting good catalytic O2 activity by double TMs. It is found that neither the O2 adsorption nor the extra charge can cause big structural deviation (Table S5), demonstrating C2N as an ideally stable substrate to support TM and deliver charges. As a demonstration, we have chosen the most stable system of Co2@C2N to investigate the full ORR through two pathways initiated by O2 or OOH dissociation in Figure 4c (Detailed structures and reaction steps are in Figure S9, S10 and S11), which are thermodynamically viable at critical electrode potentials as suggested by free energy diagram calculations (Figure S12). Considering the low dissociation barrier of OOH (0.09 eV), it is unnecessary to study the hydrogenation reaction of OOH and the pathway starting with HOOH dissociation. The results show that a four-electron O2 dissociation pathway would be kinetically favorable, in which the H2O formation through OH hydrogenation reaction is the rate-determining step with an activation barrier of 0.39 eV. This value is less than the half value of that in bulk Pt catalysts (~0.80 eV).34,35 We thus achieve ORR catalytic performance as good as bulk Pt with much less expensive yet highly durable Co-Co structure. Meanwhile, the varying TM concentration was examined by modeling a 2×1 C2N supercell with one or two TM atoms binding to different positions. The variation of geometries (Figure S13) and polarization charges (Table S6 and S7) are found to be very small. The shortest distance between two TM sites at adjacent C2N vacancies is 8.32 Å, based on which the influence of adjacent TM atoms to the neighboring C2N hole should be negligible. The solvent effect is found to be insignificant to the reaction barriers and ORR pathways (Table S8, Figure S14). In summary, using double TM atoms supported by the novel 2D material C2N, we designed low cost, high durable, and highly efficient ORR catalysts. The Cu-Cu, Co-Co and Ni-Ni systems, made of non-noble metals, exhibit both high adsorption and low activation barrier for O2. Specifically, Ni-Ni and Co-Co provide

†These authors contributed equally.

ACKNOWLEDGMENT This work was financially supported by the 973 Program (No. 2014CB848900), NSFC (No. 21473166 and 21303027), CAS Strategic Priority Research Program B (No. XDB01020000), Hefei Science Center CAS (2015HSC-UP011).

Supporting Information Available Optimized geometries and formation energies of TM1,2@C2N, atomic structures of the initial state, transition state, and final state for O2 dissociation absorbed on TM1,2@C2N. The elementary reactions along the O2 and OOH dissociation ORR pathway, and DFT-calculated free energy diagrams following O2 dissociation pathway for Co2@C2N. The water-solvent effect on O2 dissociations for TM2@C2N and ORR along O2 and OOH dissociation pathways for Co2@C2N.

REFERENCES (1) Calle-Vallejo, F.; Tymoczko, J.; Colic, V.; Vu, Q. H.; Pohl, M. D.; Morgenstern, K.; Loffreda, D.; Sautet, P.; Schuhmann, W.; Bandarenka, A. S. Finding Optimal Surface Sites on Heterogeneous Catalysts by Counting Nearest Neighbors. Science 2015, 350, 185189. (2) Ding, K.; Gulec, A.; Johnson, A. M.; Schweitzer, N. M.; Stucky, G. D.; Marks, L. D.; Stair, P. C. Identification of Active Sites in CO Oxidation and Water-Gas Shift over Supported Pt Catalysts. Science 2015, 350, 189-192. (3) Stephens, I. E. L.; Elias, J. S.; Shao-Horn, Y. The Importance of Being Together. Science 2015, 350, 164-165. (4) Yang, M.; Li, S.; Wang, Y.; Herron, J. A.; Xu, Y.; Allard, L. F.; Lee, S.; Huang, J.; Mavrikakis, M.; Flytzani-Stephanopoulos, M. Catalytically Active Au-O(OH)x Species Stabilized by Alkali Ions on Zeolites and Mesoporous Oxides. Science 2014, 346, 1498-1501. (5) Zhang, X. F.; Guo, J. J.; Guan, P. F.; Liu, C. J.; Huang, H.; Xue, F. H.; Dong, X. L.; Pennycook, S. J.; Chisholm, M. F. Catalytically Active Single-Atom Niobium in Graphitic Layers. Nat. Commun. 2012, 4, 1924. (6) Zultanski, S.; Fu, G. C. Nickel-Catalyzed Carbon-Carbon BondForming Reactions of Unactivated Tertiary Alkyl Halides: Suzuki Arylations. J. Am. Chem. Soc. 2015, 135, 624-627.

ACS Paragon Plus Environment

Page 5 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(7) Yan, H.; Cheng, H.; Yi, H.; Lin, Y.; Yao, T.; Wang, C. L.; Li, J. J.; Wei, S. Q.; Lu, J. L. Single-Atom Pd1/Graphene Catalyst Achieved by Atomic Layer Deposition: Remarkable Performance in Selective Hydrogenation of 1,3-Butadiene. J. Am. Chem. Soc. 2015, 137, 10484-10487. (8) Kyriakou, G.; Boucher, M. B.; Jewell, A. D.; Lewis, E. A.; Lawton, T. J.; Baber, A. E.; Tierney, H. L.; Flytzani-Stephanopoulos, M.; Sykes, E. C. H. Isolated Metal Atom Geometries as a Strategy for Selective Heterogeneous Hydrogenations. Science 2012, 335, 1209-1212. (9) Fu, Q.; Luo, Y. Active Sites of Pd-Doped Flat and Stepped Cu(111) Surfaces for H2 Dissociation in Heterogeneous Catalytic Hydrogenation. ACS Catal. 2013, 3, 1245-1252. (10) Ma, X. G.; Lv, Y. H.; Xu, J.; Liu, T. F.; Zhang, R. Q.; Zhu, Y. F. A Strategy of Enhancing the Photoactivity of g-C3N4 via Doping of Nonmetal Elements: A First-Principles Study. J. Phys. Chem. C 2012, 116, 23485-23493. (11) Sartorel, A.; Bonchioa, M.; Campagnab, S.; Scandola, F. Tetrametallic Molecular Catalysts for Photochemical Water Oxidation. Chem. Soc. Rev. 2013, 42, 2262-2280. (12) Hammarstrom, L. Accumulative Charge Separation for Solar Fuels Production: Coupling Light-Induced Single Electron Transfer to Multielectron Catalysis. Acc. Chem. Res. 2015, 48, 840-850. (13) Sartorel, A.; Miró, P.; Salvadori, E.; Romain, S.; Carraro, M.; Scorrano, G.; Di, V. M.; Llobet, A.; Bo, C.; Bonchio, M. Water Oxidation at a Tetraruthenate Core Stabilized by Polyoxometalate Ligands: Experimental and Computational Evidence to Trace the Competent Intermediates. J. Am. Chem. Soc. 2009, 131, 1605116053. (14) Weinstock, I. A.; Barbuzzi, E. M.G.; Wemple, M. W.; Cowan, J. J.; Reiner, R. S.; Sonnen, D. M.; Heintz, R. A.; Bond, J. S.; Hill, C. L. Equilibrating Metal-oxide Cluster Ensembles for Oxidation Reactions Using Oxygen in Water. Nature 2001,414, 191-195. (15) Lu, Y. H.; Zhou, M.; Zhang, C.; Feng, Y. P. Metal-Embedded Graphene: A Possible Catalyst with High Activity. J. Phys. Chem. C 2009, 113, 20156-20160. (16) Mahmood, J.; Lee, E. K.; Jung, M.; Shin, D.; Jeon, I. P.; Jung, S. M.; Choi, H. J.; Seo, J. M.; Bae, S. Y.; Sohn, S. D.; et al. Nitrogenated Holey Two-Dimensional Structures. Nat. Commun. 2015, 6, 6486. (17) Mahmood, J.; Jung, S. M.; Kim, S. J.; Park, J.; Yoo, J. W.; Baek, J. B. Cobalt Oxide Encapsulated in C2N-h2D Network Polymer as a Catalyst for Hydrogen Evolution. Chem. Mater. 2015, 27, 4860-4864. (18) Qiao, B.; Wang, A.; Yang, X.; Allard, L. F.; Jiang, Z.; Cui, Y.; Liu, J.; Li, J.; Zhang, T. Single-Atom Catalysis of CO Oxidation Using Pt1/FeOx. Nat. Chem. 2011, 3, 634-641. (19) Suntivich, J.; Gasteiger, H. A.; Yabuuchi, N.; Nakanishi, H.; Goodenough, J. B.; Shao-Horn, Y. Design Principles for OxygenReduction Activity on Perovskite Oxide Catalysts for Fuel Cells and Metal-Air Batteries. Nat. Chem. 2011, 3, 546-550. (20) Kresse, G.; Furthmuller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-wave Basis Set. Comput. Mater. Sci. 1996, 6, 15-50. (21) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (22) Zhou, J.; Sun, Q. Magnetism of Phthalocyanine-Based Organometallic Single Porous Sheet. J. Am. Chem. Soc. 2011, 133, 15113-15119. (23) Henkelman, G.; Uberuaga, B. P.; Jonsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901. (24) Zhang, R. Q.; Li, B.; Yang, J. L. Effects of Stacking Order, Layer Number and External Electric Field on Electronic Structures of FewLayer C2N-h2D. Nanoscale 2015, 7, 14062-14070. (25) Qi, W. H.; Wang, M. P. Size Effect on the Cohesive Energy of Nanoparticle. J. Mater. Sci. Lett. 2002, 21, 1743-1745. (26) Gong, K. P.; Du, F.; Xia, Z. H.; Durstock, M.; Dai, L. M. NitrogenDoped Carbon Nanotube Arrays with High Electrocatalytic Activity for Oxygen Reduction. Science 2009, 323, 760-764. (27) Zhang, X.; Lu, Z.; Xu, G.; Wang, T.; Ma, D.; Yang, Z.; Yang, L. Single Pt Atom Stabilized on Nitrogen Doped Graphene: CO Oxidation Readily Occurs via the Tri-molecular Eley-Rideal Mechanism. Phys. Chem. Chem. Phys. 2015, 17, 20006-20013. (28) Du, C.; Lin, H.; Lin, B.; Ma, Z.; Hou, T.; Tang, J.; Li, Y. MoS2 Supported Single Platinum Atoms and Their Superior Catalytic

Activity for CO Oxidation: A Density Functional Theory study. J. Mater. Chem. A 2015, 3, 23113-23119. (29) Liu, K.;, Kattel, S.; Mao, V.; Wang, G. F. Electrochemical and Computational Study of Oxygen Reduction Reaction on Nonprecious Transition Metal/Nitrogen Doped Carbon Nanofibers in Acid Medium. J. Phys. Chem. C 2016, 120, 1586-1596. (30) Duan, Z.; Wang, G. A. First Principles Study of Oxygen Reduction Reaction on a Pt(111) Surface Modified by a Subsurface Transition Metal M (M = Ni, Co, or Fe). Phys. Chem. Chem. Phys. 2011, 13, 20178-20187. (31) Moreno, S. L.; Romero, A. H. Atomic and Molecular Oxygen Adsorbed on (111) Transition Metal Surfaces: Cu and Ni. J. Chem. Phys. 2015, 142, 154702. (32) Bai, S.; Yang, L.; Wang, C. L.; Lin, Y.; Lu, J. L.; Jiang, J.; Xiong. Y. J. Boosting Photocatalytic Water Splitting: Interfacial Charge Polarization in Atomically Controlled Core-Shell Co-Catalyst. Angew. Chem. Int. Ed. 2015, 54, 14810-14814. (33) Bai, Y.; Zhang, W. H.; Zhang, Z. H.; Zhou, J.; Wang, X. J.; Wang, C. M.; Huang, W. X.; Jiang, J.; Xiong, Y. J. Controllably Interfacing with Metal: A Strategy for Enhancing CO Oxidation on Oxide Catalysts by Surface Polarization. J. Am. Chem. Soc. 2014, 136, 14650-14653. (34) Kattel, S.; Wang, G. F. Reaction Pathway for Oxygen Reduction on FeN4 Embedded Graphene. J. Phys. Chem. Lett. 2014, 5, 452-456. (35) Wang, Y.; Yuan, H.; Li, Y.; Chen, Z. Two-Dimensional IronPhthalocyanine (Fe-Pc) Monolayer as a Promising Single-AtomCatalyst for Oxygen Reduction Reaction: A Computational Study. Nanoscale 2015, 7, 11633-11641.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 6 of 9

Page 7 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 1. Top and side view of the structure of the monolayer C2N (h2D form) in a 2×2 supercell (a), and the structures of TM-N-hole in TM@C2N (b) and TM2@C2N (c).

Figure 2. (a) The optimized configurations of O2 adsorbed on TM@C2N (TM= Pt, Co, Ni, Cu); (b) The partial density of states (PDOS) of O2 on TM@C2N; (c)The distribution of exchange/polarized charges between O2 and TM@C2N with the isovalue of 0.003 e Å-3 (computed as the difference charge of systems before and after O2 adsorption). Cyan and yellow bubbles represent positive and negative charges, respectively. 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 9

Figure 3. (a) The optimized configurations of O2 adsorbed on TM2@C2N (TM= Pt, Co, Ni, Cu); (b) The partial density of states (PDOS) of O2 on TM2@C2N; (c) The distribution of exchange/polarized charges between O2 and TM2@C2N with the isovalue of 0.003 e Å-3 (computed as the difference charge of systems before and after O2 adsorption). Cyan and yellow bubbles represent positive and negative charges, respectively.

8 ACS Paragon Plus Environment

Page 9 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 4. (a) The dependence of O2 adsorption energy on TM polarization charge (positive) in TM@C2N and TM2@C2N; (b) The dependence of O2 dissociation barrier on its polarization charge (negative) induced by adsorption on TM@C2N and TM2@C2N. (c) The O2 and OOH dissociation initiated pathways for ORR catalyzed by Co2@C2N. The reaction energy (∆E) and activation energy (Ea) for all steps are given in parentheses with the form of (∆E, Ea).

9 ACS Paragon Plus Environment