# Structural determination of bilayer graphene on SiC(0001) using synchrotron radiation photoelectron diffraction

I. Razado-Colambo, J. Avila, D. Vignaud, S. Godey, X. Wallart, D P Woodruff, M.C. Asensio

► To cite this version:

I. Razado-Colambo, J. Avila, D. Vignaud, S. Godey, X. Wallart, et al., Structural determination of bilayer graphene on SiC(0001) using synchrotron radiation photoelectron diffraction. Scientific Reports, 2018, 8 (1), pp.10190. 10.1038/s41598-018-28402-0. hal-04438362

HAL Id: hal-04438362

https://hal.science/hal-04438362v1

Submitted on 5 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Structural determination of bilayer graphene on SiC(0001) using synchrotron radiation photoelectron diffraction

I. Razado-Colambo<fnref n="1" />,², J. Avila<fnref n="1" />, D. Vignaud³, S. Godey³, X. Wallart³, D. P. Woodruff⁴ & M. C. Asensio<fnref n="1" />

In recent years there has been growing interest in the electronic properties of 'few layer' graphene films. Twisted layers, different stacking and register with the substrate result in remarkable unconventional couplings. These distinctive electronic behaviours have been attributed to structural differences, even if only a few structural determinations are available. Here we report the results of a structural study of bilayer graphene on the Si-terminated SiC(0001) surface, investigated using synchrotron radiation-based photoelectron diffraction and complemented by angle-resolved photoemission mapping of the electronic valence bands. Photoelectron diffraction angular distributions of the graphene C 1s component have been measured at different kinetic energies and compared with the results of multiple scattering simulations for model structures. The results confirm that bilayer graphene on SiC(0001) has a layer spacing of 3.48 Å and an AB (Bernal) stacking, with a distance between the C buffer layer and the first graphene layer of 3.24 Å. Our work generalises the use of a versatile and precise diffraction method capable to shed light on the structure of low-dimensional materials.

Graphene has attracted considerable interest in recent years in the scientific, technological and industrial communities due to its novel physical and electronic properties, which offer the possibility of fabricating new devices for future carbon-based nanoelectronics<fnref n="1" />. Graphene comprises a single layer of sp²-bonded carbon atoms that are organised in an open hexagonal network. The electronic band structure of graphene shows linear dispersion around the Fermi energy at the K-point of the surface Brillouin zone, instead of the parabolic dispersion that is observed in graphite. However, its detailed electronic band structure depends on several factors such as the number of layers, the symmetry of the lattice, the stacking ordering and their interlayer distances²⁻⁴. Here we report the results of a determination of the structural parameters, namely the stacking order and interlayer distances, of bilayer graphene grown on the Si-terminated face of SiC(0001). While bulk graphite is known to have the AB stacking structure and an interlayer distance of 3.35 Å, these parameters have been the subject of some debate for bilayer and multilayer graphene.

In recent years, a number of different methods of graphene production have emerged, each method resulting in different structural and electronic properties. Conceptually, the simplest such method is mechanical exfoliation, pioneered by Novoselov and Geim⁵. Prepared in this way the material has the same AB stacking between layers as bulk graphite. Graphene can be also grown by chemical vapour deposition (CVD) onto a range of substrates such as Cu and Ni, the resulting structural properties being strongly dependent on the growth conditions. Epitaxial growth of graphene on SiC is an alternative approach capable of producing large scale, high quality, graphene films that cannot be obtained by mechanical exfoliation. Graphene grown on the two polar faces (Si-face and C-face) of SiC had been thought to have different stacking sequences that would result in different electronic properties. Low energy electron diffraction (LEED) patterns of graphene on the Si-terminated Si-face

Received: 15 March 2018

Accepted: 14 June 2018

Published online: 05 July 2018

<footnote n="1">¹Synchrotron SOLEIL & Université Paris-Saclay, 91192, L'Orme des Merisiers, Saint Aubin-BP 48, Gif sur Yvette Cedex, France. ²Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, Laguna, 4031, Philippines. ³Institut d'Electronique, de Microélectronique et de Nanotechnologie (IEMN UMR 8520), Université Lille, CNRS, Centrale Lille, ISEN, Université Valenciennes, Villeneuve d'Ascq Cedex, 59652, France. ⁴Physics Department, University of Warwick, Coventry, CV4 7AL, United Kingdom. Correspondence and requests for materials should be addressed to I.R.-C. (email: ircolambo@up.edu.ph) or M.C.A. (email: asensio@synchrotron-soleil.fr)</footnote>
appear sharp, but LEED patterns obtained from graphene on the C-face are generally smeared out into a continuous diffraction ring (albeit with some evidence of preferred rotation angles), implying contributions from either twisted adjacent layers or adjacent AB-stacked grains of different azimuthal directions. The electronic valence band structure of bilayer graphene on the Si-face shows splitting into two π band branches, but bilayer graphene on the C-face was commonly found to show linear dispersion near the Dirac point, as in monolayer graphene. However, angle-resolved photoelectron spectroscopy (ARPES) from single µm-sized grains has recently shown that multilayer graphene on the C-face has the same stacking and electronic structure as on the Si-face and also shown how different stacking sequences strongly influence the electronic structure⁶⁻⁸. There have been different theoretical predictions regarding the structural model of stacked graphene namely; AB stacking, AA stacking and twisted graphene layers⁹,¹⁰. Here we have investigated the alternative AB and AA stacking models and show that only AB stacking of bilayer graphene on the Si-face is consistent with our experimental data. In AB stacking, the two adjacent layers are displaced laterally exactly as in bulk graphite, whereas in AA stacking the carbon atoms of two adjacent layers have exactly the same lateral positions

Another important structural parameter in graphene is the spacing between the graphene layers. Density functional theory (DFT) calculations by Guo et al.¹¹ revealed the possibility of tuning the energy gap of bilayer graphene via its interlayer spacing. The calculation showed that a relatively small variation of the interlayer spacing can result in significant changes in the field-induced gap. It is therefore interesting to determine the interlayer distances between graphene layers in order to provide key information required to understand the electronic properties of multilayer graphene as well as of graphene heterostructures with other two-dimensional materials¹². There have been a number of theoretical and experimental studies that have sought to establish the interlayer spacings of multilayer graphene. DFT calculations based on the generalized gradient approximation (GGA) indicated that the interlayer distance of AB stacked bilayer graphene is larger than that of bulk graphite, and was predicted to be 3.58 Å¹³. Another DFT study, which included the effect of van der Waals dispersion forces, obtained an interlayer distance of 3.49 Å for AB-stacked bilayer graphene¹⁴, while the value found in the calculations of Guo et al. was 3.34 Å¹¹. The only experimental investigation sharing the same goals involved surface x-ray reflectivity by Hass et al.¹⁵ was directed particularly to understanding the structure of the SiC/graphene interface; the results included a determination of the interlayer distance between the carbon buffer layer (that lies between the SiC and the graphene) and the first graphene layer of 2.32 Å, but also the interlayer spacings of the first and second graphene layers (3.5 Å) and between successive graphene layers (3.35 Å). However, the signal from the graphene layers was particularly weak compared to the spots originated by the massive SiC substrate, naturally indissociable in experiment. Ferrah et al.¹⁶ reported the results of an X-ray photoelectron diffraction (XPD) of 'few layer graphene' (FLG) on SiC(0001) but their study focused on the epitaxial relationship between FLG and the SiC substrate and not on the interlayer spacings. Another photoelectron diffraction study by Lima et al.¹⁷ was described by the authors as XPD, but was actually a study of the angular distribution of photoelectrons at kinetic energies of only ~160 eV, comparable to the work we report here; the XPD name is normally reserved for experiments with much higher photoelectron energies when the angular distribution is dominated by zero-order forward 'focussing' diffraction, whereas at much lower energies both forward scattering and backscattering can be important. However, this work by Ferrah et al. was focused on the structure of the carbon buffer layer, although they also studied one preparation with monolayer graphene. Their results showed that the buffer layer is formed by two domains rotated by 60° with respect to each other, the distance between the buffer layer and the graphene layer being found to be 3.75 Å.

The objective of the investigation reported here was to determine the stacking sequence of bilayer graphene and the associated interlayer spacings in an effort to reconcile the different values found in existing theoretical and experimental literature. To achieve this goal we used angle-scan photoelectron diffraction (PED) to study single-domain graphene bilayers grown on the Si-terminated SiC(0001) surface. In this technique, the intensity of a core level photoemission peak is measured as a function of the emission angle at fixed values of the photoelectron kinetic energy. Coherent interference of the directly emitted photoelectron wavefield with other components elastically scattered by the surrounding atoms leads to an angular distribution measured outside the surface that is characteristic of the local atomic environment of the emitter atoms. This structure can be determined by comparison of the experimentally-measured angular distributions with the results of multiple scattering simulations for different possible trial structures. Trial and error iteration of the model structure is then undertaken to achieve the optimal agreement between theory and experiment.

## Results and Discussion

Figure 1 shows the LEED pattern from the sample obtained after the graphitisation procedure described in the previous section. The diffraction pattern shows the $(6\sqrt{3} \times 6\sqrt{3})R30°$ reconstruction typically observed for monolayer and bilayer graphene grown on the Si-face of SiC(0001). The large unit cell of this reconstruction is required to achieve commensuration of the lattice parameters of graphene (2.46 Å) and the SiC substrate (3.08 Å), but is already created by the carbon 'buffer layer' that forms prior to the growth of the overlayer graphene. The $(6\sqrt{3} \times 6\sqrt{3})R30°$ mesh implies that the buffer layer and succeeding graphene layers are all rotated by 30° with respect to the SiC substrate. The buffer layer is well-established to serve as a precursor to the subsequent growth of graphene and comprises C atoms arranged in a honeycomb lattice (see Fig. 2). But unlike graphene, it is strongly bonded to the underlying SiC and does not possess the typical electronic properties of graphene, notably the linear band dispersion at the K point. This layer may be reversibly turned into graphene by hydrogen intercalation¹⁸.

The quality and character of the graphene deposition in our experiments, including the number of layers, is provided by an ARPES measurement of the valence band structure. The ARPES 'image' of the valence band structure of bilayer graphene measured perpendicular to the Γ-K direction at a photon energy of 100 eV is shown in Fig. 2. The two bilayer graphene π bands are shown in this figure, together with the fitted tight-binding calculated bands. Splitting of the π band originates from the interlayer interaction between the two graphene layers.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (6√3 × 6√3)R30° LEED pattern of the bilayer graphene on SiC(0001) at 91 eV electron energy. The primitive translation vectors of the reciprocal meshes of the SiC (S) and graphene (Gr) are superimposed while two associated first-order diffracted beams are circled.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** ARPES image of the valence band structure measured at the K point of the graphene surface Brillouin Zone along the direction perpendicular to the Γ-K line using 100 eV photon energy. The Dirac point is located 320 meV below E<sub>F</sub>.</figcaption></figure>

The Dirac point ($E_D$) is located –320 meV relative to the Fermi level ($E_F$) due to charge transfer from the substrate to the graphene layer that causes it to be n-doped. There is a 120 meV gap, as expected in the bilayer graphene on SiC(0001), due to the electric dipole that exists between the graphene and SiC interface; this induces an electrostatic asymmetry between the layers<sup>4,19</sup>. These electronic features displayed in Fig. 2 confirm that the graphene sample is a single domain bilayer.

Figure 3 shows a typical C 1s photoemission spectrum taken at a photon energy of 600 eV from the bilayer graphene. The spectrum is similar to that reported by Riedl *et al.*<sup>18</sup>. There are three components, labelled SiC, Gr, and BL. The SiC component arises from C atoms in the underlying bulk SiC, Gr is attributed to emission from the graphene layer, while BL is attributed to emission from the C in the interfacial boundary layer. The Gr component was fitted with a Doniach-Sunjic function to account for the metallic behaviour and the other components were fitted with a Voigt function. The Lorentzian widths were all set to 230 meV and the asymmetry parameter used for the graphene component was 0.06. The values for the Gaussian width are 267 meV (Gr), 460 meV (SiC), and 1.36 eV (BL). Photoelectron diffraction data were extracted from the intensity modulations of the graphene C 1s component (Gr) measured at different azimuthal and polar emission angles.

The atomic plane arrangement of bilayer graphene with AB stacking is depicted in Fig. 4. The labelled polar angles of the near-neighbour interlayer forward focusing directions correspond to the 3.354 Å interlayer spacing of bulk graphite, with in-plane C-C distances of 1.42 Å. The AB stacking structure consists of planes of C atoms
<figure><img src="image_5.png" /><figcaption>**Figure 3.** C 1s core-level photoemission spectrum obtained at photon energy of 600 eV. The best-fit curve and its chemically resolved components are superimposed.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Atomic plane arrangement of bilayer graphene with AB stacking showing the different near-neighbour interlayer forward scattering directions. The angles shown correspond to the layer spacing of bulk graphite.</figcaption></figure>

each forming a hexagonal net stacked in a manner such that half of the C atoms are located directly above each other in adjacent planes, while the other half are located above the centres of the hexagon in the adjacent plane. The open hexagonal net structure leads to two inequivalent C atoms per unit cell, referred to as the α and β type atoms. The distance between the buffer layer and the first graphene layer is denoted as $d_1$ and the distance between the first and second graphene layers is defined as $d_2$. The AB stacking model is one of the two possible stacking sequences considered in the multiple scattering PED calculations, discussed in the succeeding sections.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Experimental photoelectron diffraction angular distributions from the graphene C 1s component at photoelectron kinetic energies (KE) (a) 326.1 and (b) 638.6 eV. White spots are superimposed at polar emission angles of 12.2°, 20.2°, 26.7°, 32.7° and 37.5°, from the inner to the external part of the pattern, respectively.</figcaption></figure>

The PED measurements of the graphene C 1s component presented here were obtained at two photon energies (hv), 600 and 900 eV, which correspond to 326.1 and 638.6 eV photoelectron kinetic energies, respectively. The photoelectron diffraction technique is typically performed in two different ranges of photoelectron kinetic energy that strongly influences the dominant scattering paths that contribute to the measured angular distribution. At high kinetic energies backscattering is weak and the dominant scattering processes are forward scattering zero-order diffraction with scattering angles close to zero; these 'forward focussing' peaks in the angular distribution correspond to interatomic directions. At low kinetic energies backscattering is much more significant (though forward scattering is still strong), allowing the technique to be used to locate adsorbate atoms relative to the underlying solid on which they are adsorbed. 'High' and 'low' typically means greater than, or less than, approximately 500 eV, although there is no sharp change in behaviour at this energy. The data we present here are thus recorded at intermediate energies in this range, but one may expect the higher-energy data to be more strongly influenced by forward scattering²⁰. Summarizing, as at high kinetic energy, forward scattering is dominant, the evaluation of the photoelectron diffraction patterns of conveniently selected carbon emitters at the graphene bilayer provide precise structural information of the lattice bilayer symmetry. On the contrary, at low kinetic energy, since the backscattering is significant, the photodiffraction modulations are particularly sensible to rather small variations of the interlayer spacing of the graphene bilayer and the distance to the buffer layer.

The PED angular distributions of the graphene component of the C 1s photoemission recorded at the two different photon energies are shown in Fig. 5. These were constructed from data sets recorded in either π or 2π/3 azimuthal ranges, sufficient to confirm the presence of symmetrically equivalent regions. The PED distribution at the lower kinetic energy of 326.1 eV is shown in Fig. 5(a), while Fig. 5(b) shows the distribution at the photoelectron kinetic energy of 638.6 eV. Both clearly show the expected symmetry of a well-ordered graphene bilayer crystal; at first glance this appears to be 6-fold symmetric, but more careful inspection, possible in the extracted azimuthal plots shown below, reveals the true symmetry to be 3-fold with a mirror plane (3 m symmetry). The character of the distributions is rather different, with the lower energy distribution showing bright streaks whereas the higher energy pattern has more clearly-defined narrow peaks. The higher energy distribution does not show very bright sharp peaks common in many XPD experiments and associated with interatomic directions, but this can be partly accounted for by the relatively low energy but also the weak scattering of the low atomic number C atoms. Indeed, even at higher kinetic energy (964 eV) C 1s XPD from graphite does not display very strong forward scattering peaks²¹. The lower energy distribution is expected to be more strongly influenced by backscattering that can prove to be sensitive to the detailed structure.

In order to extract more quantitative structural information, azimuthal angle-dependent plots of the lower kinetic energy graphene C 1s intensity modulations at specific polar emission angles (marked as white spots in Fig. 5 and selected where strong modulations are seen) were extracted from these full angular distributions. These azimuthal distributions are shown in Fig. 6. Careful inspection of these plots reveals the 3 m symmetry remarked upon above. For example, the shape of the peaks at azimuthal angles of ~150° and ~210° in the plot at a polar angle of 37.5° are not identical, but rather are mirror images of each other. Indeed, most of these azimuthal plots reveal a mirror plane at an azimuth of ~180°.

In order to achieve a more reliable structure determination, based on comparisons with simulated azimuthal distributions for different model structures, the data set was enlarged to cover eight polar emission angles of θ = 12.2°, 20.2°, 26.7°, 32.7°, 37.5°, 46°, 52.5° and 56.8°. Multiple scattering simulations of these photoelectron diffraction azimuthal spectra were based on the structural model depicted in Fig. 7. This comprises two graphene layers stacked in AB sequence above the buffer layer that is rotated by 30° with respect to the SiC substrate. The interlayer distance of 2.0 Å between the SiC substrate and the buffer layer employed in the calculation was based on two previous theoretical studies²²,²³ although other calculations have led to values as large as 2.6 Å²⁴,²⁵. The scattering cluster had a radius of 10 Å and contained more than 600 atoms. The two inequivalent α and β type atoms (with and without carbon atom directly below it, respectively) in the hexagonal unit mesh of graphene were
<figure><img src="image_8.png" /><figcaption>Figure 6. Experimental modulations of the graphene C 1s spectral peak as a function of azimuth at different polar emission angles with a photoelectron energy of 326.1 eV.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Figure 7. Side view (left) and tow view (right) of the structural model showing the optimised distances.</figcaption></figure>

used as emitters. Within the graphene bilayer unit mesh there are a total of four C emitters, two each from the first and second graphene layers, that contribute to the graphene C 1s photoemission component. Scattering potentials were calculated within the muffin-tin potential for atoms in a crystal lattice according to the prescription of Mattheis²⁶. Atomic mean square displacements of 0.003 Å² were used to account for thermal vibrations by means of Debye-Waller factors. The inner potential was set to V₀ = -16.5 eV, the value used in the earlier XPD study of this system²¹; because of the relatively high kinetic energies used in our study, small changes in the inner potential have no significant effect on our structural conclusions (see Supplementary Information).

The values of the parameters d₁ (the distance between buffer layer and the first graphene layer) and d₂ (the distance between the first and second graphene layers), as illustrated in Fig. 4, were varied to find the values that gave the best fit of the multiple scattering calculations to the experimental azimuthal plots. The objective test of this best fit was the minimisation of the same reliability factor (R) that has been used extensively in scanned-energy mode photoelectron diffraction²⁷ defined as

$$R = \frac{\sum(\chi_{th} - \chi_{exp})^2}{\sum(\chi_{th}^2 + \chi_{exp}^2)} \qquad (1)$$
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Reliability factor for the AB stacking model when varying the distances between the second and first graphene layers and the buffer layer. (a) Shows a two-dimensional false colour map of the R-factor, while (b) and (c) shows the variation with respect to one layer spacing at a fixed value of the other layer spacing as described in the text.</figcaption></figure>

where $\chi_{th}$ and $\chi_{exp}$ are the theoretical (simulated) and experimental modulation amplitudes and the summations are over all data points in all of the azimuthal plots. Notice that this definition is such that perfect agreement leads to a value of zero while completely uncorrelated values of the simulated and experimental modulations leads to a value of unity.

Simulations were performed over a two-dimensional grid scan of $d_1$ and $d_2$ in the range 3.00 Å to 3.90 Å, in steps of 0.1 Å, in order to identify the best-fit values. Figure 8(a) shows a false-colour contour plot of the R-factor over this grid, and clearly identifies the optimum values of $d_1$ and $d_2$ as approximately 3.3 Å and 3.5 Å, respectively. Figure 8(b) shows a one-dimensional grid scan of the variation of R as a function of the value of $d_1$ with the value of $d_2$ fixed at 3.50 Å, while Fig. 8(c) shows a similar plot of the variation of R as a function of the value of $d_2$ for a fixed value of $d_1$ of 3.30 Å. Based on these plots we conclude that the best-fit values of the $d_1$ and $d_2$ parameters are 3.24 Å and 3.48 Å respectively. In order to estimate the precision of these values we also follow the approach used in scanned-energy mode photoelectron diffraction²⁷ of defining a variance in the minimum value of the R-factor²⁸

$$Var(R_{\min}) = R_{\min} \sqrt{\overline{(2/N)}} \qquad (2)$$

where N is the number of independent pieces of structural information within the data being compared. Notice that N is not simply the number of data points; the azimuthal plots that form the basis of our experiment/theory comparison clearly show that several data points are required to define a modulation, and it is the number of modulations, or more particularly the possible number of modulations (an absence of a modulation is significant) that define N. In scanned-energy mode photoelectron diffraction, the value of N can be estimated, as proposed by Pendry²⁹ in LEED intensity-energy spectra, by dividing the total energy range by a characteristic peak width related to imaginary part of the inner potential that describes the inelastic scattering. In the present case an equivalent procedure would be to divide the total azimuthal angle range of the data (8 × 180°) by the half-angle of the narrowest observed features, i.e. ~5°. This gives a value of N of 288 so as $R_{\min} \approx 0.26$ the variance is $\approx 0.022$ and all values of the structural parameters leading to R value less than 0.302 lie within these nominal precision limits. This leads to final values for the two parameters of $d_1 = 3.24 \pm 0.20$ Å and $d_2 = 3.48 \pm 0.10$ Å. The much lower precision in the determination of $d_1$ is clearly reflected in Fig. 8(a and b). This can be attributed to the fact that the
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Comparison of experimental (data points) and simulated (continuous lines) intensity modulations of the graphene C 1s photoelectron peak as a function of azimuth at different polar emission angles with a photoelectron energy of 326.1 eV for two different structural models: (a) AB stacking and (b) AA stacking.</figcaption></figure>

PED from the Gr C 1s emitters is dependent on the location of the buffer layer only through backscattering, and at an electron energy of 326.1 eV the backscattering from the low atomic number C atoms is weak. By contrast, the data are significantly more sensitive to the graphene bilayer interlayer spacing, because this parameter is influenced by forward scattering from the lower-layer emitters as well as backscattering from the outer layer emitters. The optimised model obtained from the simulation is depicted in Fig. 7 showing all the interlayer distances.

Our conclusion that the spacing of the buffer layer and the first graphene layer, $d_1$, is less than the interlayer spacing in bulk graphite (3.35 Å), and that the interlayer spacing of the graphene bilayer, $d_2$, is larger than the value for bulk graphite, is broadly consistent with the results of previous (mainly theoretical) studies. Specifically, our value of $d_1$ of 3.24 ± 0.20 Å can be compared with the experimental X-ray reflectivity value of Hass et al.¹⁵ of 2.32 Å, and the theoretical values of Mattausch and Pankratov²⁴ (3.3 Å) and of Kim et al.³⁰ (3.35 Å), while our value of $d_2$ of 3.48 ± 0.10 Å may be compared with Hass et al.'s experimental value of 3.5 Å and the theoretically derived values of values of 3.49 Å¹⁴, 3.58 Å¹³ and 3.34 Å¹¹.

Based on previous theoretical and experimental studies, AB stacking, and not AA stacking, is the equilibrium structure for graphene grown on the Si-terminated SiC(0001) surface. To establish whether our experiments are consistent with this view we have also performed multiple scattering calculations on a model having the AA stacking sequence. Figure 9 shows a comparison of the experimental azimuthal dependence of the graphene C 1s emission at a kinetic energy of 326.1 eV with simulations for both stacking sequences. A very good agreement is achieved by using the optimised $d_1$ and $d_2$ values for the AB stacking. By contrast, the simulations for AA stacking clearly show major discrepancies between theory and experiment at polar emission angles of 20.2° and 32.7°; these angles correspond approximately to the first and second nearest neighbour interlayer scattering directions, respectively, for AB stacking (c.f. Fig. 4 in which the smaller value of $d_2$ of bulk graphite puts these angles at 23° and 36°).

This clear difference in visual comparison is also found in the associated R-factor values of 0.26 for AB stacking and 0.45 for AA stacking. Figure 10 shows a similar comparison between experiment and theory for AB and AA stacking models, for a set of graphene C 1s azimuthal scans at the higher photoelectron energy of 638.6 eV. For this data set, too, the AB stacking model yields simulations in very good agreement with experiment, whereas for the simulations based on AA stacking there is particularly poor agreement for the data recorded at 20° and 33° polar emission angles, similarly corresponding approximately to the first and second nearest neighbour forward scattering directions. Clearly these data further reinforce the existing wisdom that bilayer graphene on the Si-face of SiC(0001) has the AB stacking.

## Methods

The photoemission and photoelectron diffraction experiments were carried out at the ANTARES beamline of Synchrotron SOLEIL in Saint-Aubin, France. Photoelectron spectra were obtained using a Scienta R4000 concentric hemispherical analyser equipped with a 2-D detector to allow parallel measurement of a 25° range of collection angles. The sample was placed in the analysis chamber on a sample manipulator with five degrees of freedom, allowing precise orthogonal translation in x, y and z as well as rotation both in polar (θ) and azimuthal
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Comparison of experimental (data points) and simulated (continuous lines) intensity modulations of the graphene C 1s photoelectron peak as a function of azimuth at different polar emission angles with a photoelectron energy of 638.6 eV for two different structural models: (a) AB stacking and (b) AA stacking.</figcaption></figure>

(φ) angles. The geometry of the experiments was as follows. The surface normal, the incident beam direction and the electron emission detection direction define a common plane. The angle between the incident beam and the electron energy analyzer is 45°. Measurements of the complete photoelectron diffraction angular distributions were performed by rotating the azimuth of the sample within the 180° range in steps of 1° at specific polar angles. In order to cover a total polar angle range of around 60° this procedure was repeated at polar angles of 19° and 43°. All measurements were recorded with the sample cooled by liquid nitrogen to 120 K. The photoemission and photoelectron diffraction experiments were obtained using an incident photon beam spot size of 140 µm.

The bilayer graphene sample was prepared by the standard graphitisation procedure described in more detail elsewhere³¹. After an initial degassing at 900 °C and surface flattening at 1050 °C, the graphitisation was achieved by holding the SiC substrate at 1220 °C for six minutes, leading to high quality bilayer graphene on the Si-terminated SiC(0001) surface.

Multiple scattering simulations of the measured photoelectron diffraction data were performed using a computer program, based on the magnetic quantum-number expansion formalism, developed by Fritzsche³². The computational scheme takes into account the finite energy and angular resolution of the electron energy analyser. In the present calculations single and double scattering events were included. This program is embedded in software that searches for the best agreement between experiment and theory in a multi-dimensional parameter space using an adapted Marquardt algorithm.

## Conclusions

In summary, we have exploited the combination of photoelectron diffraction angular distributions and multiple scattering simulations to show that bilayer graphene grown on the Si-terminated SiC(0001) face has an AB stacking and has an interlayer distance between the buffer layer and the first graphene layer of 3.24 ± 0.20 Å while the interlayer spacing of the bilayer graphene is found to be 3.48 ± 0.10 Å. The smaller spacing between the buffer layer and the first graphene layer than in bulk graphite (3.35 Å), and the fact that the interlayer graphene bilayer spacing is larger than that in bulk graphite, are consistent with several previous theoretical studies. Even if the present results are more precise they are still in agreement with the only other experimental investigation, based on X-ray reflectivity.

## References

1. Geim, A. K. & Novoselov, K. S. The rise of graphene. *Nat. Mat.* **6**, 183 (2007).

2. Castro Neto, A. H., Guinea, F., Peres, N. M. R., Novoselov, K. S. & Geim, A. K. The electronic properties of graphene. *Rev. Mod. Phys.* **81**, 109 (2009).

3. Latil, S. & Henrad, L. Charge carriers in few layer graphene. *Phys. Rev. Lett.* **97**, 036803 (2006).

4. Ohta, T., Bostwick, A., Seyller, T., Horn, K. & Rotenberg, E. Controlling the electronic structure of bilayer graphene. *Science* **313**, 951 (2006).

5. Novoselov, K. S. et al. Electric field effect in atomically thin carbon films. *Science* **306**, 666 (2004).

6. Razado-Colambo, I. et al. Probing the electronic properties of graphene on C-face SiC down to single domains by nanoresolved photoelectron spectroscopies. *Phys. Rev. B* **92**, 035105 (2015).
7. Bao, C. et al. Electronic structure of trilayer graphene resolved by nanospot angle-resolved photoelectron spectroscopy. *Nano Letters* **17**, 1564 (2017).

8. Johansson, L. I. et al. Multiple pi-bands and Bernal stacking of multilayer graphene on C-face SiC, revealed by nano-angle resolved photoemission. *Scientific Reports* **4**, 4157 (2014).

9. Shallcross, S., Sharma, S., Kandelaki, E. & Pankratov, O. A. Electronic structure of turbostratic graphene. *Phys. Rev. B* **81**, 165105 (2010).

10. Lobato, I. & Partoens, B. Multiple Dirac particles in AA-stacked graphite and multilayers of graphene. *Phys. Rev. B* **83**, 165429 (2011).

11. Guo, Y., Guo, W. & Chen, C. Tuning field-induced energy gap of bilayer graphene via interlayer spacing. *Appl. Phys. Lett.* **92**, 243101 (2008).

12. Geim, A. K. & Grigorieva, I. V. Van der Waals heterostructure. *Nature* **499**, 419 (2013).

13. Yoshizawa, K., Yumura, T., Yamabe, T. & Bandow, S. The role of orbital interactions in determining the interlayer spacing in graphite slab. *J. Am. Chem. Soc.* **122**, 11871 (2000).

14. Alam, M. S., Lin, J. & Saito, M. First-principles calculation of the interlayer distance of the two-layer graphene. *Jp. J. App. Phys.* **50**, 080213 (2011).

15. Hass, J., Millán-Otoya, J. E., First, P. N. & Conrad, E. H. Interface structure of epitaxial graphene grown on 4H-SiC(0001). *Phys. Rev. B* **78**, 205424 (2008).

16. Ferrah, D., Penuelas, J., Bottela, C., Grenet, G. & Ouerghi X-ray, A. photoelectron spectroscopy (XPS) and diffraction (XPD) study of a few layers of graphene on 6H-SiC(0001). *Surf. Sci.* **615**, 47 (2013).

17. de Lima, L. H. et al. Atomic surface structure of graphene and its buffer layer on SiC(0001): A chemical-specific photoelectron diffraction approach. *Phys. Rev. B* **87**, 081403(R) (2013).

18. Riedl, C., Coletti, C. & Starke, U. Structural and electronic properties of epitaxial graphene on SiC(0001): a review of growth, characterization, transfer doping and hydrogen intercalation. *J. Phys. D: Appl. Phys.* **43**, 374009 (2010).

19. Ohta, T. et al. Interlayer interaction and electronic screening in multilayer graphene investigated with angle-resolved photoemission spectroscopy. *Phys. Rev. Lett.* **98**, 206802 (2007).

20. Woodruff, D. P. Surface structural information from photoelectron diffraction. *J. Electron Spectrosc. Rel. Phenom.* **178–179**, 186 (2010).

21. Küttel, O. M., Agostino, R. G., Fasel, R., Osterwalder, J. & Schlapbach, L. X-ray photoelectron diffraction study of diamond and graphite surfaces. *Surf. Sci.* **312**, 131 (1994).

22. Varchon, F. et al. Electronic structure of epitaxial graphene layers on SiC: Effect of the substrate. *Phys. Rev. Lett.* **99**, 126805 (2007).

23. Qi, Y., Rhim, S. H., Sun, G. F., Weinert, M. & Li, L. Epitaxial graphene on SiC(0001): More than just honeycombs. *Phys. Rev. Lett.* **105**, 085502 (2010).

24. Mattausch, A. & Pankratov, O. *Ab initio* study of graphene on SiC. *Phys. Rev. B* **99**, 076802 (2007).

25. Lampin, E., Prtiester, C., Krzeminski, C. & Magaud, L. Graphene buffer layer on Si-terminated SiC studied with an empirical interatomic potential. *J. App. Phys.* **107**, 103514 (2010).

26. Mattheis, L. F. Energy bands for solid argon. *Phys. Rev. A* **133**, 1399 (1964).

27. Woodruff, D. P. Adsorbate structure determination using photoelectron diffraction: Methods and applications. *Surf. Sci. Rep.* **62**, 1 (2007).

28. Booth, N. A. et al. Structure determination of ammonia on Cu(110)-a low-symmetry adsorption site. *Surf. Sci.* **387**, 152 (1997).

29. Pendry, J. B. Reliability factors for LEED calculations. *J. Phys. C: Solid State Phys.* **13**, 937 (1980).

30. Kim, S., Ihm, J., Choi, H. J. & Son, Y.-W. Origin of anomalous electronic structures of epitaxial graphene on silicon carbide. *Phys. Rev. Lett.* **100**, 176802 (2008).

31. Moreau, E. Elaboration de graphene par epitaxie par jets moleculaires et caracterisation, *Ph.D. Thesis, University de Lille* 1 (2011).

32. Fritzsche, V. A new spherical-wave approximation for photoelectron diffraction, EXAFS and MEED. *J. Phys. Condensed Matter* **2**, 9735 (1990).

## Acknowledgements

The authors acknowledge the support services of Synchrotron SOLEIL, especially the valuable help of Stephane Lorcry. The Synchrotron SOLEIL is supported by the Centre National de la Recherche Scientifique (CNRS) and the Commissariat à l'Énergie Atomique et aux Énergies Alternatives (CEA), France. DV, SG and XW acknowledge the financial support of the Renatech network. This work was also funded by a public grant from the French National Research Agency (ANR) as part of the "Investissements d'Avenir" (Reference No. ANR-17-CE09-0016-05).

## Author Contributions

D.V., S.G. and X.W. carried out the growth of the samples. I.R.C., J.A. and M.C.A., performed the photoemission and the photoelectron diffraction experiments. I.R.C., J.A. and M.C.A. have conducted the research design and data analysis. I.R.C. performed the multiple scattering simulations. The manuscript was written by I.R.C., M.C.A., J.A and D.P.W. All authors discussed the results and commented on the manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-28402-0.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_13.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018
Supplementary Information: Structural determination of bilayer
graphene on SiC(0001) using synchrotron radiation photoelectron
diffraction

I. Razado-Colambo,¹, ²#* J. Avila,¹ D. Vignaud,³ S. Godey,³ X. Wallart,³ D.P. Woodruff,⁴ and M.C. Asensio¹*

¹Synchrotron SOLEIL, L'Orme des Merisiers, Saint Aubin-BP 48, 91192 Gif sur Yvette Cedex, France

²Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, Laguna 4031, Philippines

³ Institut d'Electronique, de Microélectronique et de Nanotechnologie (IEMN UMR 8520), Université Lille, CNRS, Centrale Lille, ISEN, Université Valenciennes, F-59000, France

⁴Physics Department, University of Warwick, Coventry CV4 7AL, United Kingdom

*To whom correspondence should be addressed:

<footnote>* I. Razado-Colambo and M.C. Asensio, corresponding authors, with email addresses ircolambo@up.edu.ph and asensio@synchrotron-soleil.fr</footnote>

<footnote>\# I. Razado-Colambo permanent address: Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, Laguna 4031, Philippines</footnote>
<caption>a)</caption>

<figure><img src="image_14.png" /><figcaption>Stereographic projection</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>b)</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>Fermi Surface of
graphene on SiC(0001)</figcaption></figure>

<figure><img src="image_17.png" /></figure>

**Figure S1.** a) Graphene hexagonal unit cell showing the real-space lattice unit vectors and azimuthal directions, b) graphene surface Brillouin zone showing the reciprocal lattice unit vectors and high-symmetry directions, c) stereographic projection of the graphene C 1s photoelectron diffraction pattern from the bilayer graphene on SiC(0001) sample, and d) Fermi surface mapping from exactly the same sample as (c). These images show the unequivocal correspondence of the high-symmetry directions of the real space photoelectron diffraction pattern and the reciprocal space Fermi surface mapping of the same bilayer graphene sample used in this study.
To evaluate the effect of the inner potential on the azimuthal intensity modulations, an investigation of the variation of the *R*-factor as a function of the inner potential wasperformed. Figure S2 shows the dependence of the *R*-factor on the inner potential in the range from -18.0 eV to -13.0 eV. These calculations were performed using a subset of 5 of the total set of 8 azimuthal plots (a selection that leads to a lower minimum *R*-factor). The lowest value of the *R*-factor corresponds to V₀= -16.5 eV, but within the range -17.5 eV to -16.0 (the range covered in other publications [i, ii]), the *R*-factor varies by less than 1%. The resulting effect on the calculated photoelectron diffraction curves is therefore minimal, and any change in the resulting structural parameter values is much smaller than our estimated precision.

<figure><img src="image_18.png" /><figcaption>**Figure S2.** Variation of the *R*-factor as a function of the inner potential (V₀).</figcaption></figure>

i A. Grüneis, C. Attaccalite, T. Pichler, V. Zabolotnyy, H. Shiozawa, S. L. Molodtsov, D. Inosov, A. Koitzsch, M. Knupfer, J. Schiessling, R. Follath, R. Weber, P. Rudolf, L. Wirtz, and A. Rubio, Electron-Electron Correlation in Graphite: A Combined Angle-Resolved Photoemission and First-Principles Study. *Phys. Rev. Lett.* **100**, 037601 (2008).

ii F. Matsui, H. Nishikawa, H. Daimon, M. Muntwiler, M. Takizawa, H. Namba, and T. Greber, The 4πkz periodicity in photoemission from graphite. *Phys. Rev. B* **97**, 045430 (2018).