Magnetic effects in sulfur-decorated graphene.

The interaction between two different materials can present novel phenomena that are quite different from the physical properties observed when each material stands alone. Strong electronic correlations, such as magnetism and superconductivity, can be produced as the result of enhanced Coulomb interactions between electrons. Two-dimensional materials are powerful candidates to search for the novel phenomena because of the easiness of arranging them and modifying their properties accordingly. In this work, we report magnetic effects in graphene, a prototypical non-magnetic two-dimensional semi-metal, in the proximity with sulfur, a diamagnetic insulator. In contrast to the well-defined metallic behaviour of clean graphene, an energy gap develops at the Fermi energy for the graphene/sulfur compound with decreasing temperature. This is accompanied by a steep increase of the resistance, a sign change of the slope in the magneto-resistance between high and low fields, and magnetic hysteresis. A possible origin of the observed electronic and magnetic responses is discussed in terms of the onset of low-temperature magnetic ordering. These results provide intriguing insights on the search for novel quantum phases in graphene-based compounds.

Due to the two-dimensional nature and Dirac fermionic behaviour of charge carriers, the electronic and magnetic properties of graphene can be easily modified, making it one of the most appealing materials for a variety of disparate applications1234. Such easiness to access novel regimes has shifted the focus in the graphene research from graphene itself to the modification of graphene, providing an exciting and versatile platform for realization of novel phenomena and device functionality5. One of such efforts is to induce magnetic effects in graphene. Indeed a growing number of studies ranges from theoretical predictions on the intrinsic ferromagnetism or spin ordering6789101112 to experimental probes of defect/impurity-induced local magnetic moments1314151617181920. However, these magnetic moments are induced by the change of local crystal structure, e. g., carbon vacancies and deformations from to -type crystal structure, and hence, they have been identified as spin-1/2 paramagnets21. On the other hand, magnetic ordering can be realized in graphene when decorated with sulfur. Upon doping with sulfur, stacked graphene layers, i.e., graphite, exhibit ferromagnetism, which has been claimed to coexist with superconductivity2223.

In this work, we have combined two different but complementary probes such as angle-resolved photoemission spectroscopy (ARPES) to study the evolution of the graphene band structure upon sulfur introduction and magneto-transport to explore the electro-magnetic properties of the graphene/sulfur (G/S) system. Figure 1A shows the procedure adopted to prepare the G/S samples. First, graphene samples are grown epitaxially on n-doped 6 H-SiC(0001) and undoped 4 H-SiC surfaces by silicon sublimation method, as detailed elsewhere242526. The graphene sample and a piece of sulfur are then sealed in a glass ampule, with a vacuum of 10−6 Torr. The ampule is annealed at 230 °C in a furnace for 60 hours, while the pressure inside the ampule increased by vapourized sulfur was ~360 Torr. The presence of sulfur in the samples is confirmed by the observation of: a) sulfur 2 p core electrons in the photoemission spectra and b) Auger electrons corresponding to the sulfur LMM transition in the Auger electron spectroscopy (AES) spectra27. The concentration of sulfur is determined by AES and the value of sulfur/carbon ratio is 1/9.

Figure 1

Sample preparation process and energy spectra near EF.

(A) An SiC substrate is annealed at 1400 °C in ultra-high vacuum (UHV) to grow graphene, followed by another annealing process with a piece of sulfur in a sealed glass ampule at 230 °C. (B,C) ARPES intensity maps of the as-grown sample (G: panel (B)) and the graphene/sulfur compound (G/S: panel (C)), both of them measured at 10 K. The inset in panel (B) is the Fermi surface of the G sample where the black line is the direction that the energy-momentum dispersion of both samples has been taken. (D) Angle-integrated intensity of the energy spectra for the G (blue curve) and G/S (red curve) samples shown in panels (B,C), respectively. The inset shows full width at half maximum (FWHM) of the momentum distribution curves (MDCs) near the Dirac energy for G (blue circles) and G/S (red circles) samples. Both spectra show that the minimum corresponding to shifts towards upon sulfur introduction.

Changes in the electronic structure are monitored by ARPES experiments. Figure 1B,C show ARPES intensity maps as a function of energy and momentum for as-grown graphene (G) and G/S samples. Data are taken at 10 K, near the Brillouin zone corner K, along the black line shown in the inset. In line with previous reports282930, the Dirac energy, , of the G sample lies ~0.35 eV below the Fermi energy, , due to the formation of a Schottky barrier31. Two most obvious effects resulting from the introduction of sulfur are: a) broadening of the ARPES spectra; and b) small shift of towards . The former suggests that sulfur introduces disorder and the latter indicates charge transfer between sulfur and the graphene layer. More specifically, is estimated by the minimum of the integrated ARPES intensity (Fig. 1D), a quantity proportional to the one-dimensional density of states, or the full width at half maximum (FWHM) of the momentum distribution curves (MDCs, momentum spectra at constant energy shown in the inset of Fig. 1D) that shifts towards by with introduction of sulfur. This signifies a decrease of the carrier concentration by 17% from to .

A closer look at the energy distribution curves (EDCs, energy spectra at constant momentum indicated by the vertical blue lines in Fig. 1B,C) reveals a surprising depletion of states at (Fig. 2A). This is a typical signature of energy gap opening at . The temperature dependence of the EDCs for the G/S sample is shown in Fig. 2B. As the temperature T is lowered, a clear depletion of states at or a shift of the leading edge toward higher binding energy is observed. This is clearly in striking contrast to the G sample (Fig. 2C), where, as previously reported32, no gap is observed when measured below 45 K. The magnitude of the energy gap is determined by the leading edge of the EDCs with respect to through the standard procedure33 of fitting the first derivative of each energy spectrum with a Gaussian function. The temperature dependence of the gap is summarized in Fig. 2D. The gap develops below a transition temperature with a maximum shift of at 10 K and appears to be independent from the substrate. The different T-dependence of G/S from G excludes the possibility of thermal smearing as the origin of the depletion of states at . Several mechanism can be accounted for the gap opening from superstructures343536 and charge- or spin-density wave37 to magnetic ordering38 and superconductivity39. The lack of replica (shadow) bands corresponding to the additional periodicities, however, makes the formation of the former two less likely scenarios. Indeed, considering the sulfur/carbon ratio present in our G/S sample, one would expect the gap opening at a momentum of 0.3 ~ 0.4 Å−1, which is far away from 0.05 Å−1 where the energy gap is observed (Fig. 1C). Magnetism and superconductivity are certainly appealing scenarios also in view of the report of coexisting ferromagnetism and superconductivity in sulfur-doped graphite22.

Figure 2

Energy gap at EF of G/S.

(A) Energy distribution curves (EDCs) taken at 10 K for G (blue curve) and G/S (red curve) samples. (B) EDCs taken at several temperatures for G/S. The leading edge shifts away from EF with decreasing temperature. The leading edge gap Δ (roughly half the energy gap) is determined by the position of the leading edge with respect to EF. (C) EDCs taken at several temperatures for G. The leading edge stays at the same energy within the fitting error (~±0.6 meV). (D) Temperature dependence of Δ for three different samples: S1 and S2 (G/S on an SiC(0001) substrate), and C1 (G/S on an SiC substrate).

To gain better insight on the origin of the gap opening, we have performed magneto-transport measurements (Figs 3 and 4). In each case, the introduction of sulfur strongly modifies the response of the G sample. Especially, the resistance R of G/S gradually deviates and sharply increases with respect to the one of the G sample with lowering T (Fig. 3A). Figure 3B displays a ln(R) versus curve for the same data40. The linear T-dependence observed at low temperatures is typical of disordered systems in the variable range hopping regime (VRH)404142, consistent with an increase in the disorder of the G/S sample as discussed in Fig. 1C. This is qualitatively different from that of G (inset of Fig. 3B), where the linearity is not well defined due to the curvature in the whole T range.

Figure 3

Transport properties of G/S.

(A) The R versus T curves of G (blue curve) and G/S (red curve) samples. (B) The ln(R) versus curve of the G/S sample. The black-dashed line is a fit, where T0 = 2.49 K and C is an arbitrary constant. The inset shows the ln(R) versus curve of the G sample, for comparison.

Figure 4

Magnetic hysteresis of G/S.

(A) Magneto-resistance at 2 K of the G sample, with increasing (blue curve) and decreasing (dark blue curve) magnetic field, H. R0 is the maximum resistance (91 Ω) at H = 10 T. (B) Magneto-resistance at 2 K for the G/S sample, with increasing (red curve) and decreasing (dark red curve) H. R0 is the maximum resistance (6.05 kΩ) near H = 0 T. (C) The magnetoresistance of G near H = 0 T denoted by gray-shaded area in panel (A). The hysteretic effect with a coercive field of 0.003 T (blue-dashed lines) is observed. (D) The magnetoresistance of G/S near H = 0 T denoted by gray-shaded area in panel (B), showing magnetic hysteresis with two coercive fields of 0.12 T (green-dashed lines) and 0.02 T (blue-dashed lines).

In Fig. 4, we show the low-temperature magneto-resistance (R versus external magnetic field H) for both samples. The most striking difference is the behaviour of the magneto-resistance at low fields, where its complete reversal is observed upon sulfur introduction. This reversal can be explained with a transition from a weak antilocalization regime to a weak localization regime as previously reported for pure graphene4344. For this case, intervalley scattering is turned on, thus causing an increase in R near zero field. While the magneto-resistance of the G sample (Fig. 4A) is positive (increasing R upon increasing H) over the whole field range, the magneto-resistance of the G/S sample (Fig. 4B) shows a crossover at ~5 T from positive at high fields to negative (decreasing R upon increasing H) at low fields, similar to the case of weak localizations. We note, however, that the crossover field is three orders of magnitude higher than the previous results for weak localizations, where the strength of localization is suppressed very quickly in the presence of magnetic field43. Instead, the negative magneto-resistance at low fields is reminiscent of fluorinated graphene45, identified as a spin-1/2 paramagnet21.

The magneto-resistance near H = 0 T reveals another surprising phenomena of G/S, magnetic hysteresis with coercive fields of 0.12 T and 0.02 T (defined to be half the distance between the two peaks denoted by green and blue dashed lines in Fig. 4D). On the other hand, the G sample measured under the same condition (Fig. 4C) shows one order of magnitude smaller hysteretic effect with a coercive field of 0.003 T, likely induced by a remnant field of the superconducting magnet used for the measurements. Hysteresis in the magneto-resistance normally indicates the formation of magnetic domains or the existence of a magnetic granular system. The domains or grains change polarity above a given coercive field, creating a remanence in the signal. In the case of the G/S sample, the observed hysteresis persists in high fields up to 5 T, suggesting very large magnetic anisotropy and/or saturation fields. The observed magnetic hysteresis in the disordered system (broad ARPES spectrum in Fig. 1C along with the behaviour in Fig. 3B) requires that magnetic moments are correlated with electronic states of graphene in the presence of sulfur.

One of the plausible origins of the observed magnetic moments is magnetized sulfur atoms. Our first principles calculations for G/S (Fig. 5A) suggest that a non-magnetic sulfur atom, after being trapped in between graphene and the buffer layer (a carbidic layer in between graphene and the SiC substrate whose crystal structure is the same as graphene, but π bands are absent due to the interactions with the substrate) exhibits spin-polarized states via the interaction with the buffer layer (see the spin density map in Fig. 5B, where the red and blue iso-surfaces correspond to spin up and down, respectively) exhibiting a magnetic moment of 0.63 μB. This can result in indirect exchange interactions analogous to the Ruderman-Kittel-Kasuya-Yosida interactions in metals driving the spin-dependent VRH4647. Within this picture, spin-dependent hopping between neighbouring sublattices costs an energy J leading to the opening of a gap in the electronic spectra below a certain temperature3847 similar to the double-exchange mechanism in manganites, i. e., magnetism induced by adatoms with a concentration of can open a charge gap that is equal to 2 × (leading edge gap), where J is the strength of exchange interactions48. Indeed, we obtain (with and the nearest neighbour hoping parameter 49), when (at 10 K from the results in Fig. 2 with consistent with C:S = 9:1 that we have determined by Auger electron spectra). This result is similar to what the double-exchange mechanism expects50. Such spin-dependent VRH and energy gap opening are exactly demonstrated in the temperature-dependence of magneto transport signals (Figs 3 and 4) and energy spectra (Fig. 2B,D).

Figure 5

Magnetization of an intercalated sulfur atom.

(A) The top view of the crystal structure of G/S (the SiC substrate is not shown for simplicity). The dotted lines denote the unit cell of G/S compared to the unit cell of graphene. (B) The spin density of G/S. The red and blue iso-surfaces are for spin up and down, respectively.

It is interesting to note that similar behaviour has been observed at the interface between LaAlO3 and SrTiO35152. This interface is a quasi-two-dimensional metallic layer that presents ferromagnetic order and large negative magneto-resistance51. In this case, electronic charge separation has been suggested as the main ingredient for the appearance of the novel behaviour52. At the G/S sample, it is clear that sulfur removes electrons from the graphene layer (Fig. 1D). This similarity suggests that depending on the nature of the charge transfer between sulfur and carbon, nanoscale inhomogeneities are created, leading to nanoscopic droplets with different electronic and magnetic properties, similar to the LaAlO3/SrTiO3 interfaces525354. However, we should stress the stark contrast between these layered oxides and G/S: G/S is metallic not insulating, and hence electronic screening should be more efficient55, Coulomb interactions weaker, and therefore the observed behaviour is even more surprising.

In conclusion, we have revealed an unusual response to temperature and magnetic field from a material that consists of non-magnetic and light elements such as carbon and sulfur. The results here presented place graphene/sulfur as another appealing system where magnetic ordering can be realized in graphene and carbon-based materials in general, and suggest the appealing possibility to realize novel phenomena such as ferromagnetic quantum Hall effect56 for graphene-based spintronic devices operated not only by temperature and doping, but also by magnetic field.

Methods

Experiments. High-resolution ARPES experiments were performed at beamlines 10.0.1 and 12.0.1 of the Advanced Light Source at Lawrence Berkeley National Laboratory using 50 eV photons, with energy and momentum resolutions of 9 meV and 0.01 Å−1, respectively. The magneto-transport properties were measured in a Van-der-Pauw geometry for the G/S and G samples using a custom-built cryogenic vacuum probe equipped with a 10 T superconducting magnet.

Calculations. Ab initio total energy calculations were performed with a plane-wave basis set57 using the Vienna Ab-initio Simulation Package (VASP)585960. The exchange-correlation of electrons was treated within the generalized gradient approximation (GGA) as implemented by Perdew-Berke-Enzelhof 61. The crystal structure is simulated by using the supercell slab approach with vacuum region which separates the two surfaces (top and bottom surfaces). The bottom one is passivated by hydrogens to remove the lone pair of the SiC substrate. The k point sampling is (9 × 9 × 1). The structural optimisation and spin-polarized calculation were performed.

Additional Information

How to cite this article: Hwang, C. et al. Magnetic effects in sulfur-decorated graphene. Sci. Rep. 6, 21460; doi: 10.1038/srep21460 (2016).