Lukas
Madau߇
a,
Ioannis
Zegkinoglou‡
b,
Henrique
Vázquez Muiños
c,
Yong-Wook
Choi
b,
Sebastian
Kunze
b,
Meng-Qiang
Zhao
d,
Carl H.
Naylor
d,
Philipp
Ernst
a,
Erik
Pollmann
a,
Oliver
Ochedowski
a,
Henning
Lebius
e,
Abdenacer
Benyagoub
e,
Brigitte
Ban-d'Etat
e,
A. T. Charlie
Johnson
d,
Flyura
Djurabekova
c,
Beatriz
Roldan Cuenya
*bf and
Marika
Schleberger
*a
aFaculty of Physics and CENIDE, Universität Duisburg-Essen, 47057 Duisburg, Germany. E-mail: marika.schleberger@uni-due.de
bDepartment of Physics, Ruhr-Universität Bochum, 44780 Bochum, Germany
cDepartment of Physics, University of Helsinki, P.O. Box 43, 00014, Helsinki, Finland
dDepartment of Physics and Astronomy, University of Pennsylvania, Philadelphia, USA
eCIMAP, (CEA-CNRS-ENSICAEN-UNICAEN), blvd Henri Becquerel, Caen, France
fDepartment of Interface Science, Fritz-Haber Institute of the Max Planck Society, 14195 Berlin, Germany. E-mail: roldan@fhi-berlin.mpg.de
First published on 5th November 2018
Two-dimensional molybdenum-disulfide (MoS2) catalysts can achieve high catalytic activity for the hydrogen evolution reaction upon appropriate modification of their surface. The intrinsic inertness of the compound's basal planes can be overcome by either increasing the number of catalytically active edge sites or by enhancing the activity of the basal planes via a controlled creation of sulfur vacancies. Here, we report a novel method of activating the MoS2 surface using swift heavy ion irradiation. The creation of nanometer-scale structures by an ion beam, in combination with the partial sulfur depletion of the basal planes, leads to a large increase of the number of low-coordinated Mo atoms, which can form bonds with adsorbing species. This results in a decreased onset potential for hydrogen evolution, as well as in a significant enhancement of the electrochemical current density by over 160% as compared to an identical but non-irradiated MoS2 surface.
So far, two general approaches have been followed to further optimize MoS2 for the HER: (i) improving the intrinsic reactivity of MoS2 by lowering the hydrogen binding energy20,24 and (ii) increasing the total number of catalytically active sites.10–13 The latter are low-coordinated Mo atoms either located at the edges of MoS2 crystallites11 or in the vicinity of sulfur (S) vacancies on the basal planes of MoS2.12 Extensive synthesis and surface engineering efforts are usually required to achieve these goals. Some of the methods that have been suggested include the synthesis of three-dimensional (3D) shaped building blocks,10 the growth of multi-layered MoS2 sandwich structures,22 and the doping of MoS2 with metal nanoparticles.20
Because of the lower displacement threshold energy of S compared to Mo (roughly by a factor of six23), S may be selectively removed from MoS2 by either electron or ion bombardment.11,12,25 For ions with kinetic energies of up to a few 10 keV, the main interaction mechanism involves direct collisions between particles which result mainly in single or double S vacancies. Note that for 2D materials (in contrast to 3D materials) this limitation cannot be overcome by increasing the kinetic energy because collision cascades are absent in 2D materials. The efficiency is thus low, and as a consequence, a rather high fluence of roughly one incoming ion per removed S atom must be applied.
The alternative approach we present here is engineering defects via electronic energy deposition. This can be achieved by either highly charged or by very fast (swift) heavy ions. For both it has been shown that they can be used for defect engineering of 2D materials such as carbon nano-membranes,26,27 graphene,28–33 hexagonal boron nitride,34 and MoS2.35–37
Swift heavy ions (SHI) excite target atoms along their trajectory and the corresponding energy deposited per track length into the target material is usually given in terms of electronic stopping power Se = dE/dx. They offer an additional advantage that they are not stopped within a few nm as the ions in the keV range. Instead, their projected range extends up to several μm into the material. By using grazing incidence conditions, it is possible to bring a large part of this interaction volume to the surface.38,39 In this way, every single ion modifies extended surface areas, which in the case of ultrathin MoS2 leads to the creation of incisions, i.e. elongated interstices with parallel edges, which are hundreds of nm long and tens of nm wide.
This unique type of modification is due to a combination of direct energy deposition into MoS2 and the interaction of the SHI with the substrate. As was shown earlier,37,39 the dimensions of the nano-incisions agree well with the assumption that the energy transfer of the incident ion to the material can be described as a two-step process: first, the SHI excites the electronic system of the substrate material. In the second step, the energy is transferred via electron–phonon-coupling to the lattice, resulting in a so-called thermal spike in the substrate. This in turn leads to the decomposition of MoS2 and in particular, the vaporization of S. The minimum energy deposited per track length for this process was found to be Se = 2 keV nm−1.37
This work demonstrates by combining experiment and theory that SHI irradiation can be used to introduce artificial defects which enhance the catalytic activity. The SHI-induced defects include both nanometer-scale morphological features (such as incisions and Origami-like structures) and atomic scale lattice modifications (S vacancies). Both types of defects contribute to an increase of the concentration of low-coordinated Mo sites. Through the investigation of the role of defects in MoS2 HER catalysts, a fundamental understanding of structure–reactivity correlations will be obtained.
Fig. 1(b) schematically illustrates the mechanism of defect formation in MoS2 by SHI irradiation. The removal of material around the ion track leaves edges behind consisting of under-coordinated Mo atoms. Furthermore, due to the thermal spike induced by the ion, also adjacent S atoms are likely to evaporate during the process. Both the low-coordinated Mo atoms at the edges of the nano-incisions and the additional S vacancies are expected to contribute to the catalytic activity. Atomic force microscopy (AFM) images of SHI irradiated MoS2 samples, acquired after exposure of the samples to air for extended periods of time, show that the edges of the incisions are heavily decorated with adsorbates, see the inset of Fig. 1(c). This is an indication of an enhanced chemical reactivity and is also found for bilayer MoS2 (BLM) samples.
To test our hypothesis and to investigate the catalytic activity of our defect-engineered samples in detail we irradiated SLM with 91 MeV Xe ions under grazing incidence at an angle of Θ = 1° with respect to the surface and with a fluence (number of ions per area) of 15 ions per μm2. The corresponding energy deposited per track length was Se = 19 keV nm−1 into MoS2, Se = 12 keV nm−1 into SiO2, and Se = 9.3 keV nm−1 into glassy carbon.40 AFM images were obtained on irradiated samples directly after irradiation, see Fig. 1(c). The AFM images reveal elongated irradiation induced incisions in the MoS2 with an average length of 400 nm ± 200 nm on SiO2 (at 1° ion incidence with respect to the surface). Assuming that the edges are completely desulfurized we would arrive at an upper limit for the efficiency of ∼4800 S atoms per incident ion, which corresponds to a sulfur vacancy density of ca. 7.2 × 104 μm−2 for these given irradiation parameters. It should be noted that the length and the density of the incisions can be easily controlled during irradiation by varying the angle with respect to the surface and choosing a higher or lower ion fluence, respectively. At too high fluences however, tracks will start to overlap so that the achievable sulfur vacancy density will saturate at some stage.
The effect of ion irradiation on the morphology of the samples is also evident in scanning electron microscopy (SEM) images. Representative examples of both irradiated and non-irradiated MoS2 surfaces on GC are shown in Fig. 2. In addition to the incisions that are formed along the direction of the ion beam, three-dimensional structures with typical sizes of a few dozens of nm are also formed due to the bombardment of the surface with SHI. The latter consist of folded material and are a result of two incisions running very close to each other.37 These “Origami-like” structures exhibit a large surface area and a large number of edge sites and are thus expected to be beneficial for the catalytic activity of the material. From a comparison of SHI induced incisions in MoS2 deposited on SiO2 (Fig. 1) and on glassy carbon, (Fig. 2), it becomes evident that both exhibit similar morphological properties. We can thus infer that the creation mechanism is the same for the two substrates. No SHI-induced three-dimensional structures were observed in areas of the glassy carbon substrate surface which were not covered by a single layer of MoS2 (Fig. S3†). Thus, these morphological features are characteristic of the MoS2 system and not of the underlying support.
Fig. 2 SEM images of non-irradiated (a), (b) and ion-irradiated (c), (d) MoS2 monolayers deposited on GC. The irradiated surface exhibits nanometer-scale morphological features, such as incisions (∼80 nm wide and several μm long) and three-dimensional “Origami-like” structures (see the example in the circle and ESI Fig. S2 and S3†), which are characterized by a large number of under-coordinated edge sites. The non-irradiated surfaces are flat and featureless except for the intrinsic roughness of the glassy carbon substrate. The MoS2 areas in (a) and (c) are highlighted in blue. The marked sample regions are used in the calculation of the geometric surface area. |
To obtain a microscopic image of the ion induced modifications, we performed molecular dynamics (MD) simulations. Irradiation was simulated with the inelastic thermal spike model41,42 on unsupported MoS2 under grazing incidence (1° with respect to the surface) for 91 MeV Xe ions. Note that the dimensions of the sheet were 200 × 40 nm2 (for details see the ESI†), considerably shorter than the length of the incisions as measured in the experiments.
The results are shown in Fig. 3. One can easily recognize the nano-incision in the MoS2 layer, which has two amorphous edges running parallel to each other, similar to the features seen in the AFM and SEM images in Fig. 1 and 2. The amorphization along the trajectory in Fig. 3(a) and (b) impedes the quantification of the relationship between the catalytic activity and the structural changes related to low-coordinated Mo atoms and sulfur vacancies. We therefore plotted the result in a different way by removing all sputtered Mo and S atoms from the simulation box, together with the amorphous chains on the edges, see Fig. 3(c) and (d). Here, the color of the Mo atoms corresponds to their coordination number (blue = coordination number of 6, green = coordination number of 5, red = coordination number of 4). The obtained MD data suggest that on each side of the incision roughly 29% of adjacent S atoms are missing and that 28% of the edge Mo atoms are under-coordinated (higher if the amorphous chains are also considered). This means we have to refine our initial estimate of the irradiation efficiency to approximately 1400 S atoms per incident ion, corresponding to a sulfur vacancy density of ca. 2.1 × 104 μm−2 at an ion fluence of 15 μm−2.
The width of the incision in the MD simulations is smaller (∼40 nm) in comparison with the experiment (∼80 nm) probably due to neglecting the role of the substrate, which can create defects in the MoS2 in two ways. First, particles may be sputtered from the substrate resulting in further damage. This effect is assumed to be small in our case, as the cross section for nuclear stopping is negligible for SHI. However, the substrate will be heated by the thermal spike and this may also give rise to further damage. This effect cannot be neglected, in particular as the direct interaction of the beam with the 2D material is limited to a few 10 nm depending on the angle of incidence but the incisions can be much longer, i.e. a few 100 nm.
We therefore investigated the effect of the thermal spike in the substrate by applying the Two-Temperature Model (TTM).37–39 This model assumes that the energy deposited into the electronic system is transferred into the phonon system. The flow of energy from the electronic to the phonon system depends on the temperature difference between the two reservoirs and is governed by the electron–phonon-coupling constant of the material. Solving the two coupled heat diffusion equations yields spatial temperature profiles from which one can estimate whether a specific phase transition of the material occurs or not (see the ESI† for more details).
The calculated temperature profile for the SiO2 substrate irradiated with 91 MeV Xe (Se = 12.3 keV nm−1) is shown in Fig. 4. The white dotted lines depict the isotherm at 1458 K which is the decomposition temperature of MoS2.43,44 The required temperature for evaporation of S is even lower, e.g. S2 evaporates at 1140 K and in the presence of water already at 565 K by vaporizing as SO2. This indicates that indeed the substrate plays an important role for the formation of incisions. Note that we did not simulate the different substrates used for this study as we expect no significant differences in the mechanism. Quantitative differences such as deposited energy and length of incisions do exist between the different substrates. The defect creation mechanism is expected to be the same though.
As the simulations support our interpretation that catalytically active sites are indeed efficiently created by SHI irradiation we investigated our samples using X-ray photoelectron spectroscopy (XPS) to obtain information on the elemental composition and chemical state of the catalyst surface. The data for non-irradiated and irradiated MoS2 samples (the same irradiation conditions as before) are shown in Fig. 5. In both cases we measured a signal originating from the Mo 3d5/2 and Mo 3d3/2 peaks of MoS2 (blue) and MoO3 (green).45 In the case of non-irradiated MoS2 (Fig. 5(a)) the MoO3 contributes ∼20% to the overall Mo signal. This can be attributed to the already oxidized edges of the grown MoS2 flakes as for the XPS measurements the samples had to be exposed to ambient conditions. After irradiation, however, the MoO3 constitutes ∼49% of the overall Mo signal. We can easily explain this by assuming that the edges of the incisions created by ions in MoS2 provide additional catalytically active low-coordinated Mo atoms which are likely to oxidize when exposed to air.46
Apart from these low-coordinated Mo atoms, we argue that SHI irradiation also creates S vacancies next to the incisions due to the thermal spike after ion impact because the lower sublimation temperature of S compared to Mo results in a preferential evaporation of S atoms in adjacent rows next to the ion trajectory. When analyzing the stoichiometry of Mo and S in the XPS data prior to and after irradiation using the XPS peak areas and the corresponding relative sensitivity factors (RSFs),47 we indeed see a clear increase of the total Mo:S ratio. Whereas the Mo(MoS2):S ratio does not change upon irradiation maintaining its initial value of 0.45 (corresponding to an atomic composition of Mo0.91S2), the total Mo(MoS2 + MoO3):S ratio strongly increases from 0.56 to 0.88. In agreement with our hypothesis it is the significant increase of the MoO3 concentration which leads to the change in stoichiometry. Since SHI irradiation creates S vacancies, which are occupied by O atoms upon exposure of the sample to air, the increase of the amount of oxides in the XPS spectra after irradiation reflects the higher concentration of S vacancies (compared to the non-irradiated sample) prior to exposure to air. At least part of these O-filled vacancies becomes unoccupied again under the reducing conditions of the HER, thus contributing to the catalytic activity. Further XPS data and analysis can be found in the ESI.†
To test the extent to which the ion irradiation influences the catalytic activity of MoS2, we performed linear sweep voltammetry (LSV) measurements on non-irradiated and irradiated MoS2 in an aqueous 0.5 M H2SO4 solution (Fig. 6(a)). For these experiments glassy carbon was chosen as a substrate due to its high conductivity and inert behaviour during hydrogen evolution reaction (HER) measurements. The ion fluence for glassy carbon/MoS2 samples used for electrochemical characterization was chosen slightly higher, i.e. 40 ions per μm2, while the angle of incidence was kept at 1°. From the AFM analysis these parameters were estimated to yield a sufficiently high sulfur vacancy density without compromising sample integrity due to excessive overlapping of rifts. At a potential of −0.6 V vs. RHE, non-irradiated MoS2 (red curve) displays a current density of −13.3 mA cm−2, whereas irradiated MoS2 (blue curve) exhibits an almost three times higher current density of −35.3 mA cm−2. Large current values at low applied voltages are desirable for the HER as this qualifies the material under investigation as a suitable catalyst.
Another important characteristic of an efficient catalyst is the onset potential, which determines the necessary energy input to initiate the hydrogen evolution. It corresponds to the Gibbs free energy ΔGH of hydrogen bonding to the catalytic sites of MoS2.14,48 A large deviation from this value will lead to a too strong/too weak hydrogen adsorption which negatively affects the hydrogen desorption and electron–proton transfer, respectively, and ultimately the catalytic activity. It is therefore desirable to have materials with a ΔGH close to zero. In our experiment, we measured for pristine MoS2 an onset potential of −220 mV, which is consistent with previous studies for nanostructured MoS2.10,15,16,48,49 For ion irradiated MoS2 there is a strong indication that the onset potential decreases, e.g. in order to obtain a current density of 200 μA cm−2, the applied potential for non-irradiated MoS2 was −300 mV, whereas the potential for irradiated MoS2 was just −200 mV. This observation is in accordance with the work performed by Li et al.50 who theoretically increased the number of S vacancies, hence obtaining a minimum in ΔGH of approx. 0.02 eV at a S-vacancy density of 12.5%. In our experiment, the introduction of S vacancies and S-depleted regions is due to a thermal effect resulting from a heated substrate track along the ion trajectory as shown in TTM simulations and XPS data. Note, that determining an exact value of the onset potential for irradiated MoS2 is rather difficult due to possible non-faradaic behaviour at low currents.
We further analyzed the Tafel-slope and the exchange current density, see Fig. 6(b). The Tafel-slope can be viewed as a measure of the reaction kinetics. It is desirable that the value of the Tafel-slope is as low as possible.51,52
ln(−j) = ln(j0) − η/b | (1) |
Here, j represents the current density, j0 is the exchange current density, b is the Tafel slope, and η is the overpotential. In the case of non-irradiated MoS2, we measured a Tafel-slope of 160 mV dec−1 which is slightly larger compared to values reported in the literature. Irradiated MoS2 however displays two independent regions with different Tafel-slopes, both differing from the slope of the non-irradiated MoS2 sample. The region with a slope of 316 mV dec−1 is likely affected by capacitance effects and shall not be considered further. For applied potentials more negative than −0.4 V vs. RHE we measured a Tafel-slope of 104 mV dec−1, which is in good agreement with the values reported in the literature for MoS2.9,10,15,16,22,24,53 The smaller Tafel slope as compared to the non-irradiated sample is also an indication of improved catalytic properties upon SHI treatment.
The exchange current density should in principle be determined by linear extrapolation of the Tafel plot in the low potential region, i.e. where the Tafel slope for the irradiated MoS2 is equal to 316 mV dec−1 (Fig. 6). However, since this region is likely affected by capacitance effects, the exchange current density determined in this way would be overestimated.
Because our catalytic material is purely two-dimensional every part of the MoS2 is in direct contact with the GC electrode which positively contributes to an efficient electron transfer (see Fig. S5 in the ESI†). Apart from the morphological changes induced by the ion irradiation, a reaction promoting effect is also expected as a result of desulfurization of the basal planes, in line with our MD and XPS data. This effect has been previously seen also for keV Ar+ beam irradiation experiments and corroborated by ab initio molecular dynamics simulations.54,55 The important role of S vacancies in the HER catalytic activity of TMDCs has also been recently demonstrated for bulk pentlandite (Fe4.5Ni4.5S8) materials.56,57 Given that the concentration of S vacancies changes during the reaction and that the interplay between desulfurization and subsequent protonation is dependent on the applied potential, a systematic study of the S vacancy effect on the catalytic activity would require operando spectroscopic studies, which are beyond the scope of this work.
We thank U. Hagemann from the Interdisciplinary Center for Analytics on the Nanoscale (ICAN, core facility funded by the German Research Foundation, DFG) for support with the XPS measurements.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8nr04696d |
‡ Both authors have contributed equally. |
This journal is © The Royal Society of Chemistry 2018 |