Abstract
Scalable fabrication of magnetic 2D materials and heterostructures constitutes a crucial step for scaling down current spintronic devices and the development of novel spintronic applications. Here, we report on van der Waals (vdW) epitaxy of the layered magnetic metal Fe3GeTe2 (FGT)—a 2D crystal with highly tunable properties and a high prospect for room temperature ferromagnetism (FM)—directly on graphene by employing molecular beam epitaxy. Morphological and structural characterization confirmed the realization of large-area, continuous FGT/graphene heterostructure films with stable interfaces and good crystalline quality. Furthermore, magneto-transport and x-ray magnetic circular dichroism investigations confirmed a robust out-of-plane FM in the layers, comparable to state-of-the-art exfoliated flakes from bulk crystals. These results are highly relevant for further research on wafer-scale growth of vdW heterostructures combining FGT with other layered crystals such as transition metal dichalcogenides for the realization of multifunctional, atomically thin devices.
Export citation and abstract BibTeX RIS
![](https://content.cld.iop.org/journals/2053-1583/8/4/041001/revision4/tdmac171dlicense.gif)
Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 license. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
Layered magnetic materials such as Cr2Ge2Te6 [1], CrI3 [2], CrTe2 [3], and Fe3GeTe2 (FGT) [4] are considered to be promising building blocks for the development of ultra-compact spintronic devices with faster response and low-power dissipation. Moreover, their study is expected to open new perspectives on a more versatile modulation of magnetic properties, beyond what can be achieved in traditional 3D magnetic thin films [5–7]. As an example, combining these layered magnets with other 2D crystals such as graphene, transition metal dichalcogenides (TMDCs), or metal monochalcogenides to form van der Waals (vdW) heterostructures offers great potential to tailor magnetism via proximity-induced phenomena [6, 8–10].
2D magnets and related vdW heterostructures promise a large impact in future applications if they can sustain magnetic order up to room temperature [11]. Within the current catalogue of available materials, one of the most promising candidates to fulfill this requirement is FGT [4, 12]. FGT is a ferromagnetic conductor with a robust out-of-plane anisotropy [13] and, most remarkably, highly tunable properties. It has been shown that the Curie temperature (TC) in few-layer thick FGT can be boosted (from around 100 K) to room temperature, or even beyond that (up to 400 K), via electrostatic doping [4], or when grown onto a topological insulator [12]. Nearly room temperature ferromagnetism (FM) could also be achieved by modulating the composition during growth to realize the Fe-rich phases Fe4GeTe2 and Fe5−x GeTe2 [14, 15]. Besides these possibilities, strain engineering and substitutional doping have also been identified as effective routes to exert control over the magnetic properties of FGT [16–18]. Finally, although a native surface oxide forms upon air exposure [19], degradation in FGT is not as severe as in chromium halides [20], which should facilitate its realistic exploration in devices.
As for other magnetic 2D materials, experimental research on FGT has mainly concentrated on micrometer-sized flakes exfoliated from bulk crystals. This type of sample allows for exploring fundamental properties as well as device concepts based on vdW heterostacks [21]. However, it is not compatible with standard device fabrication, which usually requires a wafer-scale, uniform form of material. Recent progress has been made in this direction, with a few studies reporting on molecular beam epitaxy (MBE) of FGT films on Al2O3(0001) [12, 22], GaAs(111) [22], and Ge(111) [23]. MBE growth of FGT/MnTe [22] and FGT/Bi2Te3 [12] heterostructures, both on Al2O3(0001) wafers, has also been demonstrated. In general, the collection of data contained in these studies suggests—as usual in conventional epitaxy—that lattice-matching between FGT and its growth template is crucial to obtain large-area films with good crystalline quality. However, such requirement also reduces the number of choices in terms of suitable materials that can be employed as substrate for FGT, or combined with it to realize magnetic superlattices. Interestingly, this seems to hold true even for Bi2Te3-buffered Al2O3(0001) templates [12], on which FGT growth should in principle proceed via vdW epitaxy [24]. The detection of strain in the FGT film has been associated with the lattice mismatch between both materials (the (a = b) in-plane lattice constant of Bi2Te3 and FGT is 4.38 and 3.99 Å, respectively) [25, 26], indicating that interlayer interactions beyond vdW-bonding take place at the FGT/Bi2Te3 interface. In spite of this observation, the weak nature of vdW interactions is generally believed to alleviate the need for lattice matching conditions [27], so that vdW epitaxy of FGT on top of 2D crystals such as graphene, hexagonal BN (h-BN), and TMDCs should enable the transfer-free, scalable realization of 2D heterostructures exhibiting atomically smooth vdW interfaces. Nevertheless, the intrinsically weak bonding between 2D materials also creates difficulties namely uncontrolled nucleation on morphological defects (e.g. wrinkles) [28, 29], which can in turn result in non-uniform growth of low-quality material. Hence, exploiting epitaxial growth of FGT on other 2D crystals is greatly demanded for the development of atomically thin magnetic vdW materials to be implemented into future applications.
In this work, we assess the feasibility of growing FGT thin films directly on graphene surfaces by using MBE. We employed different characterization tools to investigate the structure and magnetic properties of the FGT/graphene heterostructures. Large-scale, homogeneous FGT films exhibiting a smooth surface morphology could be realized on graphene with a preferential epitaxial orientation and low in-plane mosaicity. Additionally, no structural damage or modification took place in graphene due to FGT growth, indicating the formation of a stable vdW interface between the materials. Magnetism in the samples were probed macro- and microscopically, revealing a ferromagnetic behavior for the FGT films that is similar to state-of-the art exfoliated flakes as well as thin films grown on 3D substrates. Overall, our results demonstrate that vdW epitaxy of FGT on graphene is a powerful approach for the wafer-scale fabrication of magnetic vdW heterostructures combining dissimilar 2D crystals for spintronics.
2. Experimental section
FGT films were synthesized by MBE (base pressure around 5 × 10−11mbar) using elemental Fe, Ge, and Te evaporated from Knudsen cells. For all growth experiments the flux ratio of Fe, Ge, and Te was kept at ∼2:1:20. The flux for each element was obtained by measuring the respective beam equivalent pressure employing a pressure gauge. Note that the flux ratio was chosen as it resulted in the epitaxial 2D growth of FGT films exhibiting structural and magnetic properties similar to those reported for stoichiometric FGT bulk crystals and thin films, as it will be shown in the next section. Differently from bulk crystal growth, for which a stoichiometric mixture of Fe, Ge, and Te powders is usually employed as starting material [13], the difference in sticking coefficient for each element on the graphene surface at the growth temperature is the main reason for the utilization of a ratio differing from the nominal composition of 3:1:2. A similar growth strategy was adopted for MBE growth of FGT on the 3D substrates Al2O3 and Ge(111) [12, 23]. Continuous films with thicknesses around 10 and 20 nm (i.e. from 12 to 25 quintuple layers (QLs); 1 QL ∼ 0.8 nm formed by sequential Te/Fe/FeGe/Fe/Te slabs) were obtained using growth times of 1 and 2 h, respectively. As substrates, 1 cm × 1 cm large epitaxial graphene on semi-insulating, on-axis 4H-SiC(0001) were utilized. They were fabricated via SiC surface graphitization at high temperatures in an Ar atmosphere, following well-established synthesis protocols [30, 31]. FGT films were also synthesized on 1 cm × 1 cm Al2O3(0001) pieces using the same conditions for comparison. The graphene/SiC(0001) substrates were outgassed at 450 °C (550 °C for sapphire) for at least 1 h and then cooled to 300 °C for FGT growth. All substrates were coated with ∼1 µm of Ti on the backside via electron beam evaporation to allow non-contact heating by radiation. In-situ growth monitoring was performed by reflection high energy electron diffraction (RHEED). Finally, in some cases the FGT films were capped in-situ with a ∼5 nm thick amorphous Te layer deposited after sample cooling to room temperature. This procedure was adopted in order to reduce FGT surface oxidation upon prolonged air exposure [23], which could hinder the successful use of the surface-sensitive techniques x-ray absorption spectroscopy (XAS), x-ray magnetic circular dichroism (XMCD), and grazing incidence x-ray diffraction (GID) (see below).
Atomic force microscopy (AFM) in tapping mode was employed to probe the surface morphology of the FGT films. Raman spectroscopy in back scattering configuration was performed with the aim of obtaining information about the structure of FGT and graphene. Raman measurements were also carried out on commercially available FGT bulk crystals (HQ graphene) for comparison. The Raman spectra were excited at a wavelength of 473 nm with the laser beam focused onto the sample surface to a 1 µm diameter spot. In order to probe the azimuthal dependence of various in-plane lattice parameters, we have used synchrotron-based GID. This experiment was performed at the SpLine beamline BM25 at the European Synchrotron Radiation Facility (ESRF) in Grenoble. An incidence angle of the illuminating x-rays of 0.2° sufficiently suppresses strong scattering by the substrate, and thus makes this method highly surface sensitive. A primary x-ray photon energy of 17 keV (λ = 0.7293 Å) enables the inspection of a comparatively large area in reciprocal space, which is important to probe also multiple reflections of the same lattice plane family.
Magneto-transport characterization was performed within the 4 K–300 K temperature range and with magnetic field up to 8 T, using a large-area van der Pauw (vdP) geometry with Al contact wires bonded at the edges of the 1 cm × 1 cm FGT film. XAS and XMCD measurements were employed to investigate the magnetic properties microscopically [32]. XAS spectra at the Fe L3 and L2 edges were recorded in the beam line BL-29 (BOREAS) at the ALBA synchrotron (Spain), which provides an ultra-high-vacuum sample environment with a base temperature of ∼3 K and magnetic field up to 6 T. Measurements used total-electron-yield detection, where the drain current was measured from the sample to the ground. The magnetic field was applied along the x-ray beam at normal incidence relative to the sample plane. XMCD was obtained by subtracting XAS spectra with the photon helicity vector antiparallel and parallel to the magnetic field (see supplementary materials (available online at stacks.iop.org/2DM/8/041001/mmedia)).
3. Results and discussion
AFM analyses revealed the formation of continuous FGT films on graphene/SiC(0001) without strong thickness inhomogeneity due to island formation or localized out-of-plane-growth. Figure 1(a) shows a typical AFM height image of a ∼10 nm thick FGT film which demonstrates its coverage over the micrometer-large graphene/SiC surface terraces. Very similar results were obtained for 20 nm thick FGT. Discontinuities (holes) in the film are observed mostly close to step edges (indicated by red arrows), resulting in relatively small areas with exposed graphene (see figure 1(b)). The absence of FGT in these regions is confirmed by phase-contrast imaging (not shown), as well as by the depth profile shown in figure 1(b) (inset).
Figure 1. (a) AFM height image of a ∼10 nm thick FGT grown on graphene/SiC(0001). The red arrows indicate the position of surface steps in graphene/SiC(0001). (b) AFM image of a surface area close to a step edge in graphene/SiC(0001) where a discontinuity in FGT is evident. The inset shows the profile corresponding to the red line in the main figure. (c) In-situ RHEED patterns taken perpendicular to the SiC<100> direction, before (upper panel) and after FGT growth (lower panel). (d) Raman spectra showing the graphene and SiC components for samples with (blue spectrum) and without (green spectrum) a ∼20 nm thick FGT film on top. The inset depicts the A1g and A1g + E2g Raman components associated with the FGT film (blue spectrum). For comparison, the Raman spectrum collected from a bulk FGT crystal is also plotted (green spectrum).
Download figure:
Standard image High-resolution imageInterestingly, a small step of ∼0.4 nm is visible in the valley region of the line profile, which is known to form at the border between n and n + 1 graphene layers as a result of the SiC decomposition process during graphene formation [33]. The appearance of such morphological irregularities in the FGT film is probably related to an inability of FGT to grow—under the employed synthesis conditions—over the bi- to few-layer thick graphene ribbons and patches located on step edge regions [30, 31, 33, 34]. They exhibit a much lower chemical reactivity in comparison to that of monolayer graphene covering surface terraces. This variation is known to affect the overgrowth of 2D materials [35]. Further optimization of FGT growth parameters and graphene substrate preparation are anticipated to mitigate this problem. The root-mean-square (RMS) roughness of the FGT film was measured to be ∼1 nm in 1 μm × 1 μm on a surface terrace, which is higher than that measured for an area of same size in graphene (RMS ∼ 0.2 nm). Note that the AFM images depicted here were acquired in non-capped FGT films exposed to air for several days. Thus, natural oxidation of the topmost few layers of the FGT film [19], as well as water and other surface contaminants, are expected to affect the surface roughness.
Figure 1(c) exhibits a typical RHEED pattern obtained during FGT growth (bottom panel), where the narrow streaks indicate epitaxial growth of a 2D film with smooth surface and good structural quality. In addition, the shift in their location with respect to the pattern obtained for the graphene/SiC template (upper panel) reveals the evolution to a larger in-plane lattice constant, as expected for FGT growth (see GID later).
Raman spectroscopy (see figure 1(d)) confirmed that the graphene layer did not suffer structural changes due to FGT synthesis. Besides the SiC background signal, the graphene-related 2D and G peaks remain mostly unaltered. The small variation in the position and width of the 2D peak can be associated with excitation of surface areas containing a different number of graphene layers, e.g. close to step edge regions [30, 31, 35]. Additionally, other effects that might take place at the FGT/graphene interface such as charge transfer, hybridization, as well as local strain fluctuations, are also anticipated to promote variations in the graphene Raman modes [35, 36]. In order to disentangle the potential contribution of each of them, future investigations are required. Importantly, the intensity of the defect-related D peak does not increase after growth. This confirms the structural integrity of graphene and indicates that no rehybridization takes place in graphene, e.g. due to the formation of covalent bonds with FGT. Finally, the Raman measurements also served to confirm the formation of FGT. The inset in figure 1(d) shows the spectral region that is relevant for this material. The two components appearing at ∼120 and ∼155 cm−1 are associated with the A1g and the sum of the A1g + E2g phonon modes in FGT [37]. The A1g mode for the FGT is slightly blue-shifted with respect to that of bulk FGT. This might be associated with the strain in the FGT film (see GID results below), or due to different degrees of degradation of the materials due to surface oxidation and/or laser exposure [37]. Investigating these issues in detail requires dedicated experiments. The strong peak located around 205 cm−1 originates from the SiC substrate [38].
Further information on the structure of FGT/graphene heterostacks was obtained by GID. Figure 2(a) shows an in-plane reciprocal space map of a Te-capped, ∼10 nm thick FGT film, in which the color-coded scattered intensity in reciprocal space is plotted as a function of reciprocal lattice units (rlu), referring to the hexagonal lattice of the SiC(0001) substrate. The corresponding absolute length (in units of Å−1) of the scattering vector is additionally provided along the Hhex axis. There are highly localized substrate reflections, namely the SiC(2.0) reflection and its higher orders. Over imposed to the map itself there is a line scan crossing the SiC(2
.0) reflection. Note that we use the four-component vector notation for hexagonal symmetry, i.e. (hklm) = (hk.m) = (hk−(h + k)m). As a guide for the eye we have drawn different radii corresponding to net planes distances from graphene, G(10.0) and G(11.0), and FGT, i.e. the FGT(n0.0) with n = [2..4] and FGT(nn.0) with [n = 1..3]. From the radial scan one can extract quantitative information on the various in-plane lattice spacing values. Based on the three most intense contributions from FGT [i.e. (11.0), (30.0) and (22.0)] one can deduce an average in-plane lattice parameter for FGT of 4.011 Å, which indicates a slight tensile strain of about 0.5% (taking as reference value dFGT = 3.991 Å for bulk crystals [26]). This tensile strain is much smaller than what has been reported for FGT grown on Bi2Te3 (1.7% tensile strain) [12], in agreement with the anticipated weak interaction between FGT and graphene. We speculate that the observed strain might be related to structural and morphological deformations in the FGT layers, such as bending taking place around surface irregularities (graphene wrinkles) and the SiC surface steps.
Figure 2. (a) GID in-plane reciprocal space map of a Te-capped 10 nm thick FGT on graphene/SiC(0001), reflecting the azimuthal dependence of various in-plane lattice parameters. The white line shows a linear scan along the SiC[2.0] direction. (b) Two azimuthal intensity profiles intersecting the maxima caused by the FGT lattice.
Download figure:
Standard image High-resolution imageFollowing the intensity distribution along the FGT arcs [e.g. the FGT(22.0)], one recognizes a pronounced preferential azimuthal orientation of the FGT layer in a way that the FGT[11.0] net planes are parallel to SiC[11.0]. Two of those angular profiles crossing the FGT(11.0) and (22.0) reflection are plotted in figure 2(b). They prove a full width at half maximum of about 4.7°, indicating a FGT lattice well-aligned with respect to the substrate. Additionally, there are weaker maxima 30° off revealing a smaller fraction of the FGT layer which appears twinned to the major variant. Complementary to that, arcs corresponding to FGT(n0.0) (see figure 2(a)) have their maxima at K = 0, which fully confirms the previous statement. Moreover, there are additional low-intensity maxima symmetric to the FGT(22.0) at 0° and ±60°, which are about 10° off the (2.0) direction. These highly symmetric sub-features are most probably due to a coincidence lattice between the FGT layer and the underlying graphene and/or SiC substrate [34]. The very same coincidence lattice might be responsible for the additional arcs next to the SiC(2
.0) and (4
.0), which contain pronounced maxima azimuthally 20° off the [2
.0] direction, or in other words: they are 10° off the [10.0] direction.
The magnetic and electrical characteristics of FGT/graphene heterostacks were investigated by transverse resistance (RXY) measurements. In ferromagnetic materials, RXY under the application of an out-of-plane magnetic field H is given by the superposition of the ordinary (OHE) and anomalous (AHE) Hall effects and can be expressed as:
![Equation (1)](https://content.cld.iop.org/journals/2053-1583/8/4/041001/revision4/tdmac171deqn1.gif)
where µ0 is the vacuum permeability and Mz the out-of-plane component of the magnetization. RO and RA are the ordinary and anomalous Hall coefficients, respectively [39]. Whereas the OHE (ROHE) contains information about the electrical characteristics, the AHE (RAHE) can be utilized to investigate the magnetic order and the magnetization reversal in ferromagnetic films. Figure 3(a) displays RXY for a 20 nm thick FGT on graphene measured at different temperatures during subsequent downward and upward sweeps of an external magnetic field. At low temperature (5 K), the AHE is clearly detected as a square shape hysteresis loop superimposed on an OHE contribution with its linear dependence on the external field (see equation 1) [23, 40]. The occurrence of the AHE reflects the ferromagnetic order in the FGT film, with the observed remanence providing evidence for a strong perpendicular magnetic anisotropy. In agreement with [22] and [23], both the coercive fields HC as well as the saturation resistance ΔRAHE decrease with increasing temperature. Thereby, ΔRAHE is defined as the value of RAHE (see equation (1)) at magnetization saturation. The temperature dependencies of HC and ΔRAHE are shown in figure 3(b). Note that one has to be cautious in regarding the absolute value of HC as a physical property of FGT due to its possible dependence on domain wall pinning effects. Nevertheless, the extrapolation of both quantities to zero indicates in average a Curie temperature of about 220 K that agrees well with what has been reported for FGT nanoflakes [40].
Figure 3. (a) Transverse resistance RXY of a 1 cm × 1 cm FGT/graphene vdP structure during downward and upward magnetic field sweeps at different temperatures. The RXY curves for temperatures above 5 K are successively upshifted by 2.5 Ω for clarity. (b) Coercive field HC and saturation resistance ΔRAHE extracted from the AHE contribution as a function of temperature. The solid lines are guides to the eye. (c) Transverse resistance RXY of the same FGT/graphene vdP structure as a function of magnetic field in the range from −4 to −8 T. (d) Landau-level index of the quantum Hall effect (QHE) plateaus in (c) as a function of the inverse magnetic field. The solid line represents the result of a linear fitting.
Download figure:
Standard image High-resolution imageFrom the slope of the OHE contribution to the Hall curves shown in figure 3(a) a low-temperature carrier density of 6.8 × 1013cm−2 is extracted which, in principle, can be ascribed to highly n-type doped graphene on SiC(0001). In contrast, FGT films are expected to exhibit a hole-like transport character with an opposite ROHE slope and have a much larger carrier density (on the order of 2 × 1015cm−2), as determined for films on sapphire substrates (see supplementary materials (available online at stacks.iop.org/2DM/8/041001/mmedia)), which is in agreement with [22]. The sheet resistivity of the FGT/graphene heterostructure (150 Ω sq−1) is found to be almost one order of magnitude smaller than the one commonly observed for our epitaxial graphene films (∼2 kΩ sq−1). These findings are explained consistently by considering parallel conduction through two transport channels for the analysis of the OHE [41, 42]. In this case, the sign as well as the absolute value of the ROHE slope depend on the relative conductivities of the FGT and graphene films. Therefore in general they do not reflect the carrier type and density of only one of the two transport channels [42].
Figure 3(c) reveals that the transverse resistance RXY actually exhibits plateaus which can be explained by the QHE in graphene or precursors of the QHE occurring at magnetic fields with magnitudes above 4 T [43]. Please note that the RXY data shown in figure 3(c) contain a certain contribution of the magnetic-field dependent longitudinal resistance due the utilized vdP geometry. The positions of the plateaus can be associated with Landau-level indices and exhibit the expected linear dependence on the inverse field (µ0 H)−1 [43, 44], as can be seen in figure 3(d). The occurrence of the QHE plateaus underlines the high quality and carrier mobility of the graphene film underneath FGT over the large-area vdP structure.
XMCD/XAS results obtained for a Te-capped, 20 nm thick FGT film on graphene/SiC(0001) are shown in figure 4. The magnetic hysteresis loop shown in figure 4(a) was recorded at base temperature of 3 K with the beam energy fixed at the maximum of the Fe L3 XMCD signal (706.4 eV) (see inset figure 4(b)). The signal is normalized to an energy point at the pre-edge region (703 eV) [45], measured on the fly with a magnetic field ramp of 1.5 T min−1 at high B (from ±6 to ±2 T) and of 1.0 T min−1 at low B (from ±2 to ∓2 T). In accordance with the magnetotransport results, the square shape of the loop reveals the strong out-of-plane ferromagnetic character of the Fe atoms in the FGT film.
Figure 4. (a) Hysteresis loop recorded at the Fe L3 edge for a Te-capped, ∼20 nm thick FGT film on graphene/SiC(0001) at 3 K. The XMCD signal is normalized to its value in saturation. (b) Temperature dependence of the Fe L3 XMCD maximum in saturation (6 T, blue circles) and remanence (red triangles). Inset: XMCD signal, expressed in percentage of the average XAS, recorded at 3 K under an out-of-plane magnetic field of 6 T (blue) and at remanence (red).
Download figure:
Standard image High-resolution imageThe temperature dependence of the maximum of the Fe L3 XMCD signal in saturation (6 T) and remanence (0 T) is plotted in figure 4(b). The remanent magnetization extrapolates to zero at about ∼220 K, which agrees closely with the TC found by transport measurements. The 6 T XMCD signal persists at room temperature, demonstrating a strong paramagnetic response even though long range ferromagnetic order has already vanished. Application of the magneto-optical sum rules for XAS and XMCD [46–48] allows us to extract the orbital, spin and total moments (morb, mspin and mTotal) of Fe 3d electrons in the ground state of FGT. The number of 3d holes was set to 4, as reported for the Fe 3d6 configuration in FGT [49], and the value of the magnetic dipole operator term was neglected. From the experimental Fe L2,3 edges in the spectra at 3 K and 6 T we obtain morb = 0.047 ± 0.004 μB/Fe; mspin = 1.13 ± 0.09 μB/Fe and mTotal = 1.18 ± 0.10 μB/Fe. The orbital-to-spin moment ratio is 0.042 ± 0.008, which is comparable to that of elemental Fe(0.043) [48], somewhat larger than what is reported for FGT(0.03) [49], but lower than that of Fe5GeTe2(0.056) [50].
In summary, vdW epitaxy of FGT magnetic films on graphene was investigated by employing MBE. The realization of large-area FGT/graphene heterostructure films with high structural and interface quality was demonstrated via a thorough structural characterization, as well as magneto-transport and XMCD investigations, which revealed a robust perpendicular magnetic anisotropy in FGT and the observation of plateaus associated with the QHE effect in epitaxial graphene. These results are highly relevant for further research on wafer-scale growth of atomically thin, all-epitaxial vdW heterostructures with multifunctional properties. By combining FGT with other selected 2D crystals including TMDCs, metal monochalcogenides, and h-BN, which can be grown by precise methods such as MBE [51–55], we envision the realization of spintronic devices with dimensions and performance not achievable with conventional bulk materials.
Acknowledgments
The authors would like to thank H-P Schönherr, C Hermann, and C Stemmler for their dedicated maintenance of the MBE system, as well as A Riedel, W Seidel, S, Rauwerdink, and A Tahraoui for technical support. Furthermore, the authors appreciate the critical reading of the manuscript by Yukihiko Takagaki. They also acknowledge the provision of beamtime under the project HC-4068 at the European Synchrotron Radiation Facility (ESRF), located in Grenoble (France). ICN2 researchers acknowledge support from the European Union Horizon 2020 research and innovation programme under Grant Agreement No. 881603 (Graphene Flagship).
Data availability statement
The data that support the findings of this study are available upon reasonable request from the authors.