Highly selective oxygen reduction to hydrogen peroxide on transition metal single atom coordination

Shifting electrochemical oxygen reduction towards 2e– pathway to hydrogen peroxide (H2O2), instead of the traditional 4e– to water, becomes increasingly important as a green method for H2O2 generation. Here, through a flexible control of oxygen reduction pathways on different transition metal single atom coordination in carbon nanotube, we discovered Fe-C-O as an efficient H2O2 catalyst, with an unprecedented onset of 0.822 V versus reversible hydrogen electrode in 0.1 M KOH to deliver 0.1 mA cm−2 H2O2 current, and a high H2O2 selectivity of above 95% in both alkaline and neutral pH. A wide range tuning of 2e–/4e– ORR pathways was achieved via different metal centers or neighboring metalloid coordination. Density functional theory calculations indicate that the Fe-C-O motifs, in a sharp contrast to the well-known Fe-C-N for 4e–, are responsible for the H2O2 pathway. This iron single atom catalyst demonstrated an effective water disinfection as a representative application.


M
olecular oxygen (O 2 ) can be electrochemically reduced to water (H 2 O) via a 4e -pathway, or hydrogen peroxide (H 2 O 2 ) with 2e -transferred in aqueous solutions.The former pathway is preferred in fuel cell applications to maximize the energy conversion efficiencies [1][2][3] , and the latter one represents a green synthetic method for H 2 O 2 [4][5][6][7] .Compared to the traditional energy, capital, and waste intensive anthraquinone process 6,8,9 , electrochemical synthesis of H 2 O 2 becomes a promising alternate with significant advantages including: (1) mild reaction conditions under room temperature and ambient pressure; (2) renewable electricity as the energy source without fossil fuel consumptions; and (3) green precursors starting with water and air.Although a variety of highly active catalysts driving the 4e -oxygen reduction reaction (ORR) have been developed to improve the performance of fuel cells [10][11][12][13][14] , there are less known catalysts that can selectively reduce O 2 to H 2 O 2 including noble metals 5,[15][16][17][18][19] and carbon materials [20][21][22][23][24][25][26][27][28][29] .A molecular level understanding of the elementary reaction steps could provide important guidance to the design of ORR catalysts for different pathways.The critical knob for searching efficient ORR catalysts relies on the proper binding strength between the reaction site and O species [27][28][29] .Taking the selective ORR to H 2 O 2 , where the O-O bond in O 2 needs to be preserved, as an example here, a toostrong interaction between the reaction site and O species could easily dissociate the O 2 molecule and direct the selectivity towards H 2 O, while a too-weak one may create a high reaction barrier to overcome 4,27 .Therefore, a materials platform with flexible tunability in electronic structures is highly desired for systematic control of ORR pathways as well as improvements in catalytic activities.
Transition-metal (TM) single atom motifs coordinated in welldefined carbon matrix, with a variety of tuning knobs such as the different metal atom centers and the varied adjacent coordinative dopants, have been attracting considerable interests in heterogeneous catalysis field [30][31][32][33] .This is mainly due to their unique electronic properties compared to their bulk metal counterparts for extraordinary activities, as well as their capability in tuning the binding strength with reaction intermediates for boosting desired catalytic pathways [34][35][36] .One representative example is the Ni single atom catalyst in our previous CO 2 reduction studies where the binding with CO was significantly weakened on isolated Ni atomic sites compared to that of bulk Ni surface to facilitate CO evolution [37][38][39] .With potentially a wide range of tunability in binding with O species via different TM single atom coordination, this class of materials as a powerful platform thus provides us with great opportunities in exploring highly selective and active catalysts for H 2 O 2 generation, as well as a flexible control of ORR pathways.
Here we report the TM single atom coordination motifs for a full range control of ORR pathways from the 2e -reduction selectively to H 2 O 2 towards the 4e -to H 2 O.A series of TM single atoms including Fe, Pd, Co, and Mn are anchored into carbon nanotube (TM-CNT) vacancies (Fe-CNT, Pd-CNT, Co-CNT, and Mn-CNT, respectively) with neighboring C, O, or N coordination for pathway tuning.Among the different TMs, Fe-CNT presents the state-of-the-art performance towards H 2 O 2 generation in terms of activity and selectivity.An unprecedented onset potential to reach 0.1 mA cm −2 H 2 O 2 generation current is achieved at only 0.822 V versus reversible hydrogen electrode (vs.RHE) in 0.1 M KOH on rotating ring-disc electrode (RRDE), while a maximum H 2 O 2 selectivity of more than 95% is delivered in both alkaline and neutral pH.With the O 2 mass transport facilitated by a gas diffusion layer (GDL) electrode, the H 2 O 2 generation rate by Fe-CNT can reach to 43 mA cm −2 with a 95.4% selectivity under only 0.76 V.By switching the neighboring O with N coordination, the 2e -ORR pathway is successfully shifted towards 4e -of H 2 O, demonstrating a wide range of reaction tunability in this materials platform.Density functional theory (DFT) calculations suggest that the catalytically active C and Fe sites in Fe-C-O and Fe-C-N motifs are responsible for the H 2 O 2 and H 2 O pathways, respectively.In a variety of Fe-C-O motifs calculated, the incorporation of Fe atoms significantly improves their catalytic activities for H 2 O 2 generation compared to those with only O dopants.As a prototype demonstration of potential applications, this high-performance H 2 O 2 generation catalyst enables an effective water disinfection of >99.9999% bacteria removal at a treating rate of 125 L h −1 m −2 electrode .

Results
Synthesis and characterizations of single atom catalysts.A small quantity of TM cations (~0.1 at%) were first dispersed onto commercial surface-functionalized CNTs as the carbon matrix suspended in water (Supplementary Fig. 1), followed by freezedry and thermo annealing under Ar flow at 600 °C (Methods, Supplementary Note 1 and Supplementary Figs.2-3) 38 .All four types of TM-CNT samples, including Fe, Pd, Co, and Mn, have shown similar structures in Fig. 1 by transmission electron microscopy (TEM) and aberration-corrected high-angle annular dark-field scanning TEM (HAADF-STEM).No nanoparticles or clusters were observed in the bright field TEM images by different scales (Fig. 1a-d; Supplementary Figs.4-7), suggesting a well dispersion of TM atoms.Isolated TM atoms can be resolved by HAADF-STEM due to their high Z-contrast compared to those neighboring light elements such as C or O 32 .While all four isolated metal atoms were observed as the white dots in Fig. 1e-h, Pd-CNT presents the most distinguishable single atoms due to its heaviest atomic mass compared to the other three metal elements.
Although no obvious TM signals were detected by X-ray photoelectron spectroscopy (XPS, Supplementary Fig. 8) due to the low mass loading, synchrotron-based X-ray absorption spectroscopy (XAS) on the K-edge of Fe 40 , Pd 41 , Co 42 , and Mn 43 can instead provide direct evidence of the valence states and coordination environments 38 .The X-ray absorption near-edge structure (XANES) spectra in Supplementary Fig. 9 indicate the partial-oxidation states of these metal dopants where the pre-edge peaks are located between the corresponding metal foils and metal oxides.This is mainly due to the strong M-C and M-O bonds via which the electrons from the metal atomic centers are partially depleted to the neighboring C and O sites, in a good agreement with the simulated charge density distributions and Balder charge analysis (Supplementary Table 1) 44 .In addition, the oxidation state of coordinated Fe is lower than simply adsorbed Fe on CNT, suggesting the different chemical environment between the adsorption case and coordination case (Supplementary Note 2, Supplementary Fig. 10).The corresponding Fourier transforms of extended X-ray absorption fine structure (EXAFS) for each metal are plotted in Fig. 1i-l.It is noted that the signals from bulk metal foils are dominated by M-M bonds, i.e., 2.17 Å for Fe-Fe, 2.55 Å for Pd-Pd, 2.18 Å for Co-Co and 2.31 Å for Mn-Mn, respectively [40][41][42][43] .In contrast, M-C or M-O bonds at ~1.4-2.0Å are the major features in M-CNTs, reinforcing the isolated distribution of these TM single atoms in the carbon matrix.A shoulder peak at 2.05 Å is observed for Mn-CNT, which is corresponding to the Mn-O coordination 43 .
Due to the capabilities in both three dimensional (3D) imaging and chemical composition measurements at the atomic scale, atom probe tomography (APT) becomes a complementary characterization to STEM and XAS to reveal additional information about the atomic structure in catalysts 37,[45][46][47][48] .Fe-CNT was analyzed as a representative of other M-CNTs.We used a "sandwich" approach for the sample tip preparation (Fig. 2).We drop-casted dispersed Fe-CNT onto a nickel-coated substrate, sandwiched in between two 25-nm gold layers deposited by sputter-coating; the gold serves as a marker for the original surface of the CNT powder (Fig. 2a, Supplementary Fig. 11).The Au-CNT-Au sandwich is itself enclosed in layers of sputterdeposited Ni, to give sufficient material supports for APT tip preparation.Focused ion beam (FIB) procedures are applied to cut the sample into pieces, with a corresponding cross-section structure shown in the SEM image of Fig. 2b.Preparation of an APT-compatible needle-shaped specimen is carried out by FIB as shown in Fig. 2c.APT was then performed on a Cameca LEAP 4000 HR as detailed in the Methods, with a total of 9 million atoms collected.The reconstructed tomography of the sample tip is presented in Supplementary Movie 1, with the projected 2D image shown in Fig. 2d.Each dot represents the position of the ion detected by mass spectrometry, with a maximum detection efficiency of 38% (reflectron atom probe).Au layer as a marker in yellow is observed in the upper left corner, indicating the catalyst region for analysis.A few hundreds of isolated 56 Fe atoms (identified only from the 56 Da peak due to molecular fragment overlaps at the other Fe isotopes and charge states) were well dispersed across the sample tip and highlighted in green.A smaller region attached to the Au marker was selected for detailed analysis as shown by the 2D contour map of carbon atom distribution (Fig. 2e), where the shape of a CNT was resolved.Fe atoms were distributed uniformly along the CNT as shown in the corresponding 2D projection (Fig. 2f), with no neighboring Fe atoms observed with a distance less than 3 Å.This is in consistent with STEM and XAS characterizations and further confirms the atomic dispersion of Fe in CNT.While a closer observation around Fe atoms is not able to provide the exact atomic structures due to the loss of ion signals during the detection process as well as the APT resolution (0.1-0.3 nm), it could still give us important hints of the possible coordination environment: the side view suggests a few atomic layers where Fe single atoms were anchored; the top view of one atomic layer, although with vacant space since a portion of atoms were not detected, shows that the isolated Fe atoms may have neighboring coordination with both C and O, suggesting a possible coordination of Fe-C-O motifs which will be further discussed in our theoretical calculations.
Electrocatalytic ORR characterizations.The ORR performances were evaluated in 0.1 M KOH by casting a thin catalyst layer onto rotation ring disk electrode (RRDE), with the collection efficiency     16) 5,18,[22][23][24][25] , such as Pd-Hg (0.72 V) 15 , Au-Pd (0.54 V) 17 , Pt single atoms (0.71 V) 16 , and highly oxidized CNTs (0.75 V) 27 , representing a facile ORR kinetics with negligible overpotential for O 2 -to-H 2 O 2 conversion 5,19,[23][24][25][26] .By switching the metal dopants from Fe to Pd, Co, and Mn, the H 2 O 2 selectivity was changed to 90.  19).The catalytic activity of Fe-CNT in both RRDE test and bulk electrolysis presents significant improvements compared to previous studies of oxidized carbon catalysts (Supplementary Fig. 20).The performance stability of Fe-CNT single atom catalyst was also demonstrated on RRDE in Fig. 3c, with a stable H 2 O 2 selectivity of above 90% over the 8-h continuous operation (Methods).Post-catalysis XAS analysis of Fe K-edge XANES overlaps well with that of pristine Fe-CNT, suggesting that the electronic structure and coordination of Fe single atoms remains unchanged (Supplementary Fig. 21).The corresponding Fourier transformed EXAFS spectrum of post-catalysis Fe-CNT reveals that Fe atoms still maintain an atomic dispersion.The reaction pathway can also be tuned by maintaining the metal center while switching its neighboring metalloid coordination, which combined with the corresponding changes in catalytic performances could further reveal the possible active coordination motifs for H 2 O 2 generation.The H 2 O 2 selectivity of Fe-CNT was significantly decreased to a maximum of 60% when the catalyst was annealed in forming gas with Fe-C-O coordination reduced (Red.Fe-CNT); the 4e -ORR pathway was preferred when O was replaced with N to form Fe-C-N coordination, with the electron transfer number boosted from 2.09 of Fe-CNT to 3.71 of Fe-N-CNT 13,[49][50][51][52] , and even to 3.90 when the mass loading was increased to a typical fuel cell test condition (Supplementary Fig. 22) 13,53 .
Theoretical calculations.The above experiments suggested that the product selectivity of ORR could be related with the coordination motifs of Fe-C-O or Fe-C-N in those single atom catalysts.While based on current experimental techniques it is still challenging to directly resolve all of the exact atomic structures, density functional theory (DFT) calculations of different configurations (a total of 18 configurations as shown in Fig. 4a) become of great importance to understand the origin of selectivity change, and for rational design of new catalysts in the future.We used a two-dimensional graphene to model the CNT support used here, since the difference in binding free energies of reaction intermediates at CNT with a diameter larger than 20 nm and at the two-dimensional graphene is <0.1 eV 54 .We considered carbon vacancies up to six, tested all the possible sites of metal atom binding, and selected the most stable configurations.Representative formation energies (Supplementary Note 4) of single atom configurations and surface adsorption of metal atoms suggest that the metal single atoms prefer to chemically bind with C and O in CNT vacancies, forming stable M-C-O configurations other than being simply adsorbed on CNT (Supplementary Table 2, Supplementary Fig. 23).In addition, the effects of O and N were carefully investigated.As a prototypical system, we investigated Fe doped system, which was found to be the most active and selective for H 2 O 2 production in the experiments.It has been shown in the past that the oxygen reduction catalytic activity and selectivity can be related to the binding free energies of reaction intermediates 4,5,[27][28][29]   ) j (mA cm -2 ) j HO À 2 in pH > 11.6), the catalytic activity can be determined by the binding strength of the key reaction intermediate, *OOH.For the catalysts with weak tendency toward *OOH adsorption, O 2 activation (Eq. 1) is the potential limiting step (PDS), while desorption of *OOH (Eqs.2.1 and 2.2) is the PDS for strongly *OOH binding catalysts.
We therefore use ΔG *OOH as a descriptor to estimate the limiting potential (U L ) of O 2 reduction to H 2 O 2 .The maximum U L , i.e., zero overpotential at the top of the volcano, is 0.68 V vs. RHE originated from the experimental formation free energies of H 2 O 2 .There is a small difference in the formation free energy of its deprotonated anion of HO À 2 , shifting the equilibrium potential positively by 70 mV to 0.75 V vs. RHE in alkaline solutions with pH higher than 11.6.We note that an onset potential around 0.8 V, which is higher than the thermodynamic potential of 0.75 V, was observed in our experimental results.This phenomenon has been reported in previous studies 5,[26][27][28]55 . Themain driving force for such difference has been suggested to be the Nernstrelated potential shift due to the low concentration of H 2 O 2 in the electrolyte.We present the volcano plots in the form of H 2 O 2 with the 0.68 V equilibrium potential throughout this work unless stated otherwise to represent a wider pH range for general mechanism understanding.It should be noted that catalysts positioned at the left leg of the volcano plot are expected to be less selective toward H 2 O 2 .This is due to their strong *OOH adsorption energy, which favors dissociating the O-O bond and driving the 4e -pathway with H 2 O as the major product 5,27-29 .Our results show that in the examined Fe-C systems without O, Fe, and the C atoms around Fe present U L much smaller than 0.68 V and therefore do not contribute to the ORR activity (Supplementary Fig. 24).Once oxygen is included into the systems, most of the C atoms in the vicinity of the Fe-C-O motif (Fig. 4b) show very high activity for H 2 O 2 production and are positioned near the top of the volcano.The fact that most of the C sites are mainly located at the right leg of the volcano, implies that the C atoms of Fe-C-O motifs could also be selective for 2e - product of H 2 O 2 over 4e -to H 2 O.This trend remains valid in strong alkaline solutions with the formation of HO À 2 as shown in Supplementary Fig. 25, as in the theoretical model the difference in the equilibrium potential only results in the change of the peak position by 0.07 eV.In addition, we found out that the incorporation of Fe atoms in Fe-C-O motifs, compared to those with only O dopants, can generally strengthen the *OOH binding on C sites (Supplementary Fig. 26).This as a result significantly improves the H 2 O 2 activity towards the top of the volcano, suggesting the critical roles of both Fe and O atoms in the coordination.We also calculated the ORR activity of Fe sites in different Fe-C-O model systems, and realized that those Fe sites generally bind the *OOH too strongly, resulting in O-O bond dissociation (Fig. 4b).The calculated U L s of the Fe in variety of Fe-C-O are <0.3V, indicating a very low ORR activity.Therefore, we suspect that the C site in Fe-C-O motif is where O 2 gets selectively reduced to H 2 O 2 .This is in consistent with our in-Operando XANES spectra where the oxidation state of Fe remained unshifted, suggesting no additional Fe-O bonding formed during the reaction conditions (Supplementary Fig. 27).
To investigate the effect of local coordination environment on selectivity, we also modeled variety of Fe-C-N configurations, where O atoms are replaced by N (Supplementary Fig. 24c).Our results show that the Fe sites in all the Fe-C-N motifs are rather active for ORR to H 2 O with *OOH adsorption energies positioned in the range of interests for 4e -pathway 4 .This is in a very good agreement with our experimental results as well as previous literature, that Fe-C-N system is highly active for ORR in fuel cell applications 13,56 .Contrary to the Fe-C-O coordination, C sites of Fe-C-N system bind *OOH very weakly (Supplementary Fig. 24c), indicating that C sites are not active for 2e -nor 4e -ORR processes.The Fe active site in Fe-N-CNT was again confirmed by our in-Operando XANES experiment that, the oxidation state of Fe was slightly increased due to new Fe-O bonds formed under reaction conditions (Supplementary Fig. 27) 50,[57][58][59] .Additional calculations of different M-C-O motifs in M-CNT catalysts as well as CNT support were presented by volcano plots in Supplementary Fig. 28, where the general trend of H 2 O 2 selectivity matches with our experimental results.
Water disinfection by Fe-CNT catalyst.A wide range of practical applications could be realized with this low-cost Fe-CNT catalyst for highly efficient H 2 O 2 generation.Since H 2 O 2 has been widely used in killing bacteria 60 , one promising field is the delocalized or green-route water disinfection with accessible inputs including sunlight for electricity, air for O 2 , and water as shown in the schematic in Fig. 5a.Here we performed a prototype experiment to test the catalyst's disinfection effectiveness.Neutral pH needs to be used instead of alkaline solutions to mimic the practical applications, therefore the ORR selectivity of Fe-CNT was first evaluated in 0.1 M PBS electrolyte using RRDE as shown in Figs.5b and 5c.H 2 O 2 generation starts at ~0.53 V and maintains a high selectivity above 90% from 0.5 to 0.3 V.The practical electrolysis was performed in an H-cell where Fe-CNT catalyst was casted onto a GDL electrode (0.5 mg cm −2 catalyst loading), with the catalytic performance plotted in Fig. 5d.The potential to deliver a 20 mA cm −2 constant current for H 2 O 2 generation remained unchanged over the course of electrolysis (Fig. 5e).Around 1613 ppm H 2 O 2 was generated within 210 min electrolysis as determined by the colorimetric quantification method, representing an average Faradaic Efficiency of 90.8%.With those performance metrics obtained, electrolyte with Escherichia coli (E.coli) was then used as a model system at a bacteria concentration of ∼10 7 colony forming units (c.f.u.) mL −1 .The disinfection process was monitored by picking up several droplets during the 20 mA cm −2 chronopotentiometric measurement, followed by serially dilution and spread plating onto LB agar for overnight culture 61,62 .The dark-field CCD images of agar plates with cultured bacteria colony are shown in Fig. 5f, with the calculated killing rate plotted in Fig. 5g.Fe-CNT demonstrates a rapid disinfection efficiency for E. coli, delivering a 43% bacteria inactivation in 5 min and more than 99.9999% in 120 min (equals to a 125 L h −1 m 2 electrode processing rate) with no recovery observed (Supplementary Fig. 32).We further suspended E coli bacteria into 2-h electrolyzed 0.1 M PBS solution (ca.1000 ppm of generated H 2 O 2 concentration) with the image of overnight cultured plate shown in Supplementary Fig. 29, which clearly excludes any effects of disinfection from applied potential.This device can be further upgraded into a flow cell with tunable water flow rate and current density for a wellbalanced water disinfection rate and efficiency for specific applications in the future.

Discussion
Our experimental and theoretical results highlight that the TM single atom coordination motifs can effectively tune the ORR pathways and product selectivity.Among different catalysts examined, Fe-C-O coordination was identified as highly active and selective motif for O 2 reduction to H 2 O 2 .Given the many choices of TM atoms, supports, as well as non-metal dopants, this wide range of coordination combinations could enrich the pool of catalytic active sites for improving, controlling, and understanding catalytic reactions beyond ORR.Future efforts also could be focused on electrolysis cell design for pure H 2 O 2 solution generation without solute mixtures.In addition, by integrating renewable electricity from solar or wind with water and air, our earth-abundant Fe-CNT H 2 O 2 catalyst could therefore pave a green way, different from the traditional anthraquinone process, for H 2 O 2 generation, water treatment, chemical synthesis, and many other important applications in the future.

Methods
Material synthesis and characterization.M-CNT catalysts were prepared by the impregnation and reduction method.In a typical synthesis of Fe-CNT, a 7.5-mM iron nitrate stock solution was first prepared by dissolving Fe(NO 3 ) 3 •9H 2 O (ACS Grade, Alfa Aesar) into Millipore water (18.2MΩ•cm).A carbon suspension was prepared by mixing 50 mg multi-walled carbon nanotubes (Carbon Nanotubes Plus GCM389, used as received) with 20 mL of Millipore water, and tip sonicated (Branson Digital Sonifier) for 30 min till a homogeneous dispersion.Then 200 µL of Fe 2+ solution, given a raw atomic ratio of Fe:C to be ~0.1 at.%, was dropwise added into CNT solution under vigorous stirring, followed by a quickly frozen in liquid nitrogen.The as-prepared Fe(NO 3 ) 3 /CNT powder was heated up in a tube furnace to 600 °C at a pressure of 1 Tor and a gas flow of 100 sccm Ar (UHP, Airgas) within 20 min, and kept at same temperature for another 40 min before cooling down to room temperature.Other Pd-, Co-, and Mn-CNTs were prepared in a similar way to Fe-CNT except for various metal salt precursors, i.e., Pd (NO 3 ) 2 •2H 2 O, Co(NO 3 ) 2 •6H 2 O, and Mn(NO 3 ) 2 •4H 2 O (Puriss or ACS Grade, Sigma-Aldrich), respectively.N doped Fe-N-CNT was prepared by heating up the above-mentioned Fe(NO 3 ) 3 /CNT powder under a same temperature program with Fe-CNT but within a mixed gas flow of 50 sccm NH 3 (anhydrous, Airgas) + 100 sccm Ar.
Probe-corrected FEI Titan Themis 300 S/TEM with ChemiSTEM technology was used for S/TEM characterization.Drift correction was applied during acquisition.SEM was performed on a Zeiss Supra55VP field emission scanning electron microscope with in-lens detector.X-ray photoelectron spectroscopy was obtained with a Thermo Scientific K-Alpha ESCA spectrometer, using a monochromatic Al Kα radiation (1486.6 eV) and a low energy flood gun as neutralizer.Thermo Avantage V5 program were employed for surface componential content analysis.
Ex situ XAS spectra on Co, Mn K-edge were acquired at the SXRMB beamline of Canadian Light Source using a 4-element Si drift fluorescence detector, and Pd K-edge spectra were acquired in fluorescence mode with a 32-element Ge detector at the HXMA beamline.The powder sample was spread onto doublesided, conducting carbon tape.In situ electrochemical XAS measurement on Fe K-edge was carried at Beamline 8-ID, National Synchrotron Light Source II, Brookhaven National Laboratory, using a Si(111) monochromator and a Lytle detector.For in situ spectroelectrochemical tests, continuous O 2 flow was delivered into a homemade Teflon H-cell filled with O 2 -saturated 0.1 M KOH, a Kapton film covered gas diffusion layer (Fe-CNT/GDL) working electrode was served as the Xray window for synchrotron radiation.Analyses of both the near edge (in energy scale) and extended range (in R space) XAS spectra were performed using Athena software.
Nanofabrication and APT characterization.Specimens for the 3D Atom Probe Tomography were prepared using an FEI Helios 660 Nanolab Dual-Beam FIB/ SEM.Liftout was performed following standard APT sample prep procedures; sharpening of the tips was performed first with an ion beam accelerating voltage of 30 kV, and then (for final shaping and cleanup) at 2 kV, to produce tips with apices of ~30 nm diameter.3D APT was performed using a Cameca LEAP 4000 HR instrument in laser pulsing mode, with the following operating conditions: laser pulse energy of 100 pJ and repetition rate of 100 kHz; base temperature of 40 K; detection rate of 0.5%.The dataset shown in the main text consisted of slightly more than 9M ions, of which the majority were Ni from the capping layer above the Au-CNT-Au sandwich.No specimens survived to show the lower Au layer, likely due to the high porosity of the CNT powder material.Analysis of that data was performed using Cameca IVAS 3.16.4software, with reconstructions based on SEM images of the tip profile.Mass spectrum peaks were ranged using the fullwidth-at-one-tenth-maximum method.Only one peak was ranged as Fe, that being the peak at 56 Da; no other Fe isotopes could clearly be identified both due to mass overlaps with other species and due to the low number of detected atoms in this peak (sufficiently low that the other Fe isotopes would be expected to be undetectable below the background noise even without mass overlaps).This implies that the stated conclusions about Fe in this system, as determined by APT, represent an upper limit on the possible number and distribution of Fe atoms around the CNT.
Electrochemical measurements.A BioLogic VMP3 work station was employed to record the electrochemical response.Certain amounts of KOH (Reagent Grade, Sigma-Aldrich) or K 2 HPO 4 /KH 2 PO 4 (ACS Grade, Sigma-Aldrich) was dissolved in Millipore water to prepare the 0.1 M electrolyte.The rotation ring disk electrode (RRDE) measurements were run at 25 °C in a typical three-electrode cell.A platinum foil (99.99%,Beantown Chemical) and a saturated calomel electrode (SCE, CH Instruments) were used as the counter and reference electrode, respectively.A RRDE assembly (AFE6R1PTPK, Pine Instruments) consisting of a glassy carbon rotation disk electrode (Φ = 5.0 mm) and a Pt ring (Φ = 15.0 mm) was used, with a theoretical collection efficiency of 25%.Experimentally, the apparent collection efficiency (N) was determined to be 24.1% in the ferrocyanide/ferricyanide half reaction system at a rotation rate between 400 and 2025 rpm (Supplementary Fig. 9) 63 .To prepare M-CNT cast working electrode, typically, 3.3 mg of as-prepared M-CNT catalyst was mixed with 1 mL of ethanol and 10 µL of Nafion 117 solution (5 %, Sigma-Aldrich), and sonicated for 20 min to get a homogeneous catalyst ink.6 µL of the ink was pipetted onto glassy carbon disk (0.196 cm 2 area, 0.1 mg/cm 2 mass loading), got vacuum dried prior to usage.As the catalyst can be dispersed very well in ethanol solutions, uniform catalyst coating can be made on the disc electrode without obvious pin holes or uncovered edge.All potentials measured against SCE was converted to the reversible hydrogen electrode (RHE) scale in this work using E (vs.RHE) = E (vs.SCE) + 0.244 V + 0.059*pH, where pH values of electrolytes were determined by Orion 320 PerpHecT LogR Meter (Thermo Scientific, i.e., 13.0 for 0.1 M KOH and 13.9 for 1 M KOH).This SCE reference electrode was further calibrated to freshly prepared RHE prior to usage (Supplementary Fig. 13), in good agreement with the calculated values.The diffusion limited current of our catalyst is at around 3 mA cm −2 and is higher than that of bare glassy carbon, which may due to the low density of active sites on mirror polished glassy carbon electrode.A gradual degradation of ring current was observed during the continuous RRDE stability test, which was mainly due to the surface oxidation of Pt ring electrode constantly operated at high potential (1.2 V, Fig. 3c) and can be readily recovered by rapid cyclic voltammetry at low potentials to reduce PtO x .H 2 O 2 selectivity was calculated using the following equation: H 2 O 2 (%)=200 I Ring =N I Disk þI Ring =N , and the electron transfer number (n) at the disk electrode during ORR was calculated using n ¼ 4 I Disk j j I Disk þI Ring =N , where I Ring is the ring current, I Disk is the disk current and N is the collection efficiency.
Bulk H 2 O 2 production in 1 M KOH was carried out in a customized H-cell electrolyzer, with 0.5 mg cm −2 Fe-CNT air-brushed onto a 1×2.5 cm 2 Freudenberg GDL electrode (Fuel Cell Store) as ORR cathode, with the anode of 0.2 mg cm −2 IrO 2 /GDL for water oxidation.A Fumasep FAA-3-PK-130 anion exchange membrane (Fuel Cell Store) was employed to separate the chambers.H 2 O 2 concentration was quantified by cerium sulfate titration based colorimetric method (2Ce 4+ + H 2 O 2 → 2Ce 3+ + 2 H + + O 2 ) 27 .The H 2 O 2 concentration-absorbance curve was calibrated by mixing known amount of commercial H 2 O 2 solution with 1 mM Ce(SO 4 ) 2 .The absorption at 320 nm wavelength was measured on a Cary 5000 UV-Vis-NIR spectrometer (Agilent) and used to determine the Ce 4+ /H 2 O 2 concentration (Supplementary Fig. 18).To fit the linear range of calibration curve, the electrolyte collected was further diluted by 10 to 100 times in 0.5 M H 2 SO 4 .
Computational details.Density functional theory calculations were performed using the Vienna Ab Initio Simulation Package (VASP) 64,65 .We used BEEF-vdW exchange-correlation functional, which has been shown to accurately describe chemical and physical interactions between adsorbates and graphene 66 .We set an energy cutoff, a convergence criteria for self-consistent iteration and ionic relaxation to be 500 eV, 10 −4 eV and 0.05 eV/A ̊, respectively.Bulk graphene unit cell was optimized using (12 × 12 × 1) k-points mesh, resulting in C-C distance to be 1.424 Å. (7 × 7) supercell consisting of 98 carbon atoms with 15 Å of a vacuum perpendicular to the graphene plane was employed to model carbon materials, and we considered up to six carbon vacancies.Various sites for a single metal atom adsorption were assessed, and the effects of oxygen or nitrogen near the single metal atom were taken into account as well.For the adsorption calculations, (2 × 2 × 1) k-points mesh was utilized.
We analyzed *OOH adsorption on the catalysts, where the possible catalytic active site is either the single metal or nearby carbon sites.For single atom catalysts, the single metal atom is significantly under-coordinated compared to their bulk counterparts resulting in markedly strong binding of adsorbates, and it is highly likely that the site is pre-occupied by other adsorbates, such as *O and *OH, under the ORR conditions.In this sense, for metal site adsorptions, we first determined the most relevant coverage of the metal atom at 0.7 V RHE by constructing a surface Pourbaix diagram (Supplementary Figs.29-31).For carbon site adsorptions, we considered all carbon sites near the single metal atom and reported the most stable binding free energies.
To construct a free energy diagram of ORR to H 2 O 2 , we corrected the calculated electronic energies by adding zero-point energy, enthalpy and entropy of adsorbate at 300 K obtained from a harmonic oscillator approximation using Atomic Simulation Environment 67 .Since O 2 molecule is poorly described by standard DFT calculations, we used the calculated energies of gas-phase H 2 O and H 2 molecules, and experimental formation free energies of H 2 O and H 2 O 2 as references.To take into account the effect of the electrode potential, computational hydrogen electrode (CHE) method was employed 68 .In this method the chemical potential of proton and electron pair μ(H + +e -) is equal to one half of that of gas-phase hydrogen molecule ð1=2μðH 2 ÞÞ at standard conditions, and the effect of the electrode potential is included by shifting the electron free energy by −eU elec , where e and U elec are elementary charge and electrode potential, respectively.Water disinfection.A standard lab bacterial Escherichia coli strain was kindly provided by Howard Berg's laboratory.E. coli was cultured to stationary phase in LB broth for 14 h at 37 °C, harvested by centrifugation at 800 × g, washed three times with 0.1 M PBS solution and suspended in 0.1 M PBS to ~7.1 × 10 6 c.f.u.mL −1 .The electrochemical disinfection measurements were run at 25 °C in a H-type glass cell separated by a Fumasep FBM bipolar membrane (Fuel Cell Store).25 mL of the prepared E. Coli in 0.1 M PBS was injected into the cathodic chamber, with a 1 cm 2 Fe-CNT/GDL (0.5 mg cm −2 catalyst loading) serving as the working electrode.A chrono-potentiometric measurement at a fixed current density of 20 mA cm −2 was performed to ascertain water disinfection.Another control experiment was run in 25 mL of bacteria-free electrolyte to calculate the overall H 2 O 2 concentration and generation rate.The pH value was noted to maintain ~7.2 before and after the continuous electrolysis in PBS buffer.Bacterial concentrations and killing rates were measured at different time points during electrolysis using standard spread plating techniques.Each sample was serially diluted and each dilution was plated in triplicate onto LB agar plates, and incubated at 37 °C for 12 h.The images of overnight cultured plates were taken with a custom-built dark-field imager equipped with a CCD camera (Point Grey Chameleon).

Fig. 1 Fig. 2
Fig. 1 TEM and XAS characterizations of M-CNT catalysts.High resolution TEM and aberration-corrected HAADF-STEM images of a, e Fe-, b, f Pd-, c, g Co-, d, h Mn-CNT are shown in upper and middle panels, respectively, in which multiple rolled graphene layers can be clearly distinguished with no metal clusters/particles observed in each sample.The bright dots in HAADF-STEM (marked by red arrows) represent some typical metal single atoms.Scale bars, 5 nm.i-l Corresponding Fourier transformed EXAFS spectra of M-CNTs are plotted in R space at the bottom panel in comparison to their bulk metal foil and metal oxide control samples

Fig. 4 2 6
Fig. 4 DFT calculations of the ORR activity and selectivity on different motifs.a All examined configurations for single Fe atom coordinated in twodimensional carbon material with and without O species.We used the most stable surface coverage as a reference (Supplementary Figs.29-31).Green, red, and gray colors denote Fe, O (or N), and C atoms, respectively.b The calculated ORR activity volcano plot for 2e -pathway to H 2 O 2 .Red and blue symbols indicate *OOH adsorption at C and Fe, respectively.The equilibrium potential of O 2 /H 2 O 2 is shown as a black dashed line.Some of the points are not shown in the volcano plot since their binding energies are out of the range

Fig. 5
Fig. 5 Disinfection performance of Fe-CNT in neutral pH. a Schematic of electrochemical synthesis of H 2 O 2 for water disinfection, with green inputs such as sunlight, air, and water.b, c LSV of Fe-CNT catalyst on RRDE at 1600 rpm in 0.1 M PBS (pH 7.2) with the corresponding H 2 O 2 selectivity under different potentials.d LSV of Fe-CNT catalyst on GDL electrode (0.5 mg cm −2 mass loading) in an H-cell electrolyzer.e Bulk electrolysis at a constant current density of 20 mA cm −2 in 0.1 M PBS containing ~10 7 c.f.u.mL −1 E coli bacteria.The H 2 O 2 concentration of 1613 ppm at 210 min was determined from a bacteria-free control experiment.f CCD photos of overnight cultured plates with spread droplets taken from different time slots during the electrolysis.Dilution factor is labeled at the right bottom corner of each image.No bacteria colonies were observed after 2-h treatment.g The disinfection efficiency as a function of treatment time.The error bar represents two identical cultured plates RRDE, with a maximal H 2 O 2 selectivity of more than 95%, and a high potential of 0.822 V vs. RHE to deliver a 0.1 mA cm −2 H 2 O 2 onset current.This early onset is superior to the so-far reported H 2 O 2 catalysts as listed in the performance chart (Supplementary Fig. Fig. 3 ORR performance of M-CNT catalysts cast RRDE in 0.1 M KOH. a Linear sweep voltammetry (LSV) of CNT background and Fe-, Pd-, Co-, Mn-CNT catalysts recorded at 1600 rpm and a scan rate of 5 mV s −1 , together with the detected H 2 O 2 currents on the ring electrode (upper panel) at a fixed potential of 1.2 V vs. RHE.b Calculated H 2 O 2 selectivity and electron transfer number during potential sweep.c Stability measurement of Fe-CNT at a fixed disk potential of ~0.71 V.The ring electrode was refreshed several times by rapid scan at negative potentials to remove the accumulated PtO x in continuous operation.d, e LSV and corresponding H 2 O 2 selectivity comparison on Fe-CNT, Fe-N-CNT, and forming gas reduced Fe-CNT (Red.Fe-CNT) catalysts performance evaluated by O 2 activity and selectivity.Compared to bare CNT, the performance was gradually increased with the increase of Fe atom loading, but dramatically dropped once Fe clusters was formed (Supplementary Note 1 and Supplementary Fig.2), demonstrating the critical role of atomically dispersed Fe.Fe-CNT maintains its high H 2 O 2 selectivity and activity when applied onto a GDL electrode with facilitated O 2 mass diffusion for large current densities in electrolyzer, where the colorimetric quantification of H 2 O 2 was employed instead (Supplementary Fig.18, Methods).In 1 M KOH, the catalyst delivered a steady-state H 2 O 2 partial current of 43 mA cm −2 at 0.76 V with a Faradaic efficiency of 95.4%, corresponding to a H 2 O 2 production rate of ~1.6 mol g −1 h −1 or 8 mol m −2 h −1 (Supplementary Fig.
. In the 2e -reduction of O 2 to H 2 O 2 (H 2 O 2 in pH < 11.6 or its deprotonated anion but