Abstract
The combination of two-dimensional (2D) materials into heterostructures enables the formation of atomically thin devices with designed properties. To achieve a high-density, bottom-up integration, the growth of these 2D heterostructures via van der Waals epitaxy (vdWE) is an attractive alternative to the currently mostly employed mechanical transfer, which is problematic in terms of scaling and reproducibility. Controlling the location of the nuclei formation remains a key challenge in vdWE. Here, a focused He ion beam is used to deterministically place defects in graphene substrates, which serve as preferential nucleation sites for the growth of insulating, 2D hexagonal boron nitride (h-BN). Therewith a mask-free, selective-area vdWE (SAvdWE) is demonstrated, in which nucleation yield and crystal quality of h-BN are controlled by the ion beam parameters used for defect formation. Moreover, h-BN grown via SAvdWE is shown to exhibit electron tunneling characteristics comparable to those of mechanically transferred layers, thereby lying the foundation for a reliable, high-density array fabrication of 2D heterostructures for device integration via defect engineering in 2D substrates.
Similar content being viewed by others
Introduction
Through the combination of two-dimensional (2D) materials, novel van der Waals heterostructures with tailored properties can be designed for applications in atomically thin and transparent optoelectronic, electronic, and sensing devices1,2,3,4. However, the currently most employed integration via mechanical transfer remains challenging with regard to reliability as well as scalability for a future ultrasmall-size device fabrication. Therefore, the growth of 2D materials on top of each other via van der Waals epitaxy (vdWE) is an attractive alternative to form atomically thin heterostructures1,5,6.
Within 2D heterostructures, hexagonal boron nitride (h-BN) plays a central role in two ways: Electronically, its large band-gap enables efficient electronic confinement (as required in ultra-thin layers)3, while its chemical inertness (i.e. the lack of dangling bonds), together with its insulating nature and high mechanical strength, make it an ideal substrate or encapsulation layer for other 2D materials7. A detailed overview of various techniques used for the synthesis of h-BN can be found elsewhere8. Scalable fabrication techniques such as chemical vapor deposition or molecular beam epitaxy (MBE) were used for the growth of h-BN/graphene heterostructures9,10,11,12. The weak out-of-plane interactions between 2D materials lead to preferential nucleation at defects and morphological features in the 2D substrates in vdWE, giving rise to non-uniform growth12,13,14,15. Up to now there is no feasible way to control the exact location and orientation for the nucleation of 2D materials, making it one of the key challenges to achieve high-quality material for electronic applications16. Previous studies on the selective heteroepitaxy on vdW substrates concentrated on the growth of multiple 2D layers on mica or 3D nanostructures on graphene using photo- or electron beam lithography and plasma etching for processing a mask17,18,19,20. However, these processes have a limited resolution (particularly in the case of photolithography) and the masking process threatens to contaminate or damage the 2D substrates.
Instead of a mask, we used a focused ion beam (FIB) within a He ion microscope to deliberately create atomic-scale defects in epitaxial graphene (EG) on SiC with perfect position control. The use of EG as a substrate is advantageous for two reasons: (i) the sublimation of Si from SiC leads to a reproducible formation of epitaxially aligned graphene with low intrinsic defect density and (ii) it can be directly synthesized on semi-insulating SiC with no need for a transfer process. According to molecular dynamics simulations, irradiation of 2D materials with ions of noble gases with energies in the keV range results in the formation of defects with control down to single vacancies21,22,23. In this regard, a He FIB was recently used to write optical active defects in MoS224. In our approach, the irradiation-induced defects (IIDs) in EG exhibit dangling bonds, which act as preferential nucleation sites for atomically thin h-BN grown from elemental B and N sources via plasma-assisted MBE (PAMBE) (Fig. 1aâc). Our mask-free SAvdWE of h-BN is superior to photo- or electron beam lithography, as no residues from the mask threatens to contaminate the surface of the 2D substrate. Ideally, only a single C atom needs to be removed from pristine graphene, as graphene itself serves as a mask through its chemically inert nature. Unlike heavier Ga ions, which have been used for the formation of nanoholes in graphite flakes and the selective growth of few-layer h-BN within, the lighter He ions are better suitable to form atomic-scale defects in 2D substrates and therefore for the selective growth of atomically thin heterostructures21,25.
Results and discussion
Deterministic placement of defects in graphene via He FIB
The processing of EG on semi-insulating SiC substrates results in step edges on the surface with up to 20ânm height (visible in the dark-filed micrograph in Fig. 1d) and several μm wide terraces in between (see Supplementary Fig. 1), which are covered with single- (SLG) and bi-layer graphene (BLG). Using a He FIB with an energy of 30âkeV and focal spot sizes below 1ânm, the IIDs were written in a pattern, which was aligned to marker structures on the SiC substrate (see method section as well as Supplementary Fig. 1). The He FIB patterns are located in between three cross-shaped markers in the upper left and lower right corners, encompassing a structured area of 40âÃâ40âµm2, as visible in the darkfield micrograph in Fig. 1d. The 16 fields of 10âÃâ10âµm2 are separated by additional, single cross-shaped markers written by He FIB. Each of the upper 12 fields contains 361 spot irradiations in a square grid of 19âÃâ19. The spot irradiations define the positions of the intended defect locations and are arranged with a grid spacing of either 250 or 500ânm with two repetitions to cover different areas on the EG (see annotations in Fig. 1d).
Based on molecular dynamics simulations the sputtering rate of C atoms from graphene upon He ion irradiation with an energy of 30âkeV is expected to be very low23. However, the defect formation in 2D materials is not only determined by the ion mass and its energy, but is also affected by the interactions with the substrate, e.g. by back-scattered ions or atoms sputtered from the substrate. Hence, despite the low sputtering rate, by using enough He ions the FIB could remove more than only one C atom per intended location of the IID. Therefore, the same square grids were written not only with varying distances, but also with varying ion numbers (10,000, 20,000, or 40,000) per intended location of the IID, to study the influence of the amount of He ions on the defect formation and the consequent growth of h-BN. The two 10âÃâ10âµm2 fields in the lower right of the pattern remained blank without IIDs, serving as a reference for an unstructured sample area within the pattern. In the two squares to the lower left more complex shapes for the IIDs were drawn with the He FIB (not shown here).
Figure 1e compares the intensities of the D (~1375âcmâ1) and G (~1600âcmâ1) Raman peaks of graphene in the different areas, with the mean spectra for the denser pattern with varying number of ions used per IID shown in Fig. 1f (the mean spectra of the pattern with 500ânm spacing is shown in Supplementary Fig. 2). For higher ion numbers a small increase of the D/G peak intensity ratio can be observed, indicating a slightly higher defect density26. Since the number of IID locations is constant in each pattern, this indicates the formation of more than one defect per intended defect location for higher amounts of He ions. For the G and 2D peak a shift to higher wavenumbers can be observed for the structured samples, which we attribute to charge transfer from adsorbates at the IIDs to the graphene after exposure to ambient conditions27. Mapping the D peak intensity, we can distinguish the structured area of 19âÃâ19 defects in a dense pattern (with 40,000âions per IID) from the surrounding pristine EG (inset of Fig. 1f).
Selective-area van der Waals epitaxy of h-BN on graphene
Following the IID preparation, the substrates were introduced to the MBE system for the SAvdWE of h-BN islands (see Method section). Figure 2a presents an atomic force microscopy (AFM) image of an area of 19âÃâ19 IIDs, each written with 20,000âions, in a distance of 250ânm (areas with 500ânm distance can be seen in Supplementary Fig. 2). Areas with a rougher surface can be distinguished from smoother areas close to the step edges, which we attribute to overgrown SLG and BLG, respectively, as shown in previous publications12,14. In EG, SLG is known to have a different chemical reactivity than BLG, due to an atomic corrugation of the C-rich and partially sp3-bonded interface layer to SiC, which is screened by an additional graphene layer in BLG, leading to the different growth behavior28,29. As the focus of the current study lies on the SAvdWE of h-BN islands at IIDs, we will concentrate in the discussion on the growth of h-BN at IIDs on BLG, where the influence of the interface layer, which is specific for EG on SiC, is negligible. The corresponding conductive AFM (cAFM) image of the same area (Fig. 2b) shows a clear contrast between graphene and h-BN islands, which appear as insulating areas around the intended locations for the IIDs within the markers. Besides the nucleation at the IIDs, insulating h-BN can be also seen nucleating at the step edge and scarcely at natural defects, similar to what was observed in the previous studies14. In the detailed AFM image shown in Fig. 2c, a mean diameter of 150ânm of the h-BN islands can be observed on the BLG areas. Their average height of 0.34ânm (see also profile along the blue dotted line in Fig. 2e) corresponds well to previously reported values for h-BN islands grown on graphene and graphite11,14. The distance of ~250ânm between the islands and their square arrangement precisely follow the FIB pattern. For single-crystalline h-BN islands either a triangular or hexagonal shape would have been expected30, but the islands here show a non-regular shape. Additional h-BN layers can be observed close to the center of the h-BN islands, together with ~2ânm high nanoparticles (NPs), similar to the formation reported from naturally occurring defects in EG14. Fig. 2d shows the corresponding cAFM image of the same area, where a clear contrast is observed not only between graphene and h-BN islands, but also within the islands, where the central area with the additional layers and NPs shows no tunneling current at a bias of 100âmV.
The Raman spectra of these samples shown in Fig. 2f reveal a substantially increased intensity at ~1360âcmâ1, where the E2g peak of h-BN is expected31. Similar to previous studies a shift of the G and 2D peak can also be observed, which can be attributed to a change in strain and doping of graphene12. In aligned h-BN/graphene heterostructures fabricated via mechanical exfoliation a strained-induced broadening of the 2D peak (as we observe it here) was also reported32. Comparing the mean spectra of the patterned areas with non-patterned areas, a slightly higher intensity of the h-BN related Raman peak can be observed, which we attribute to the denser nucleation of h-BN.
Controlling selective nucleation and orientation
In Fig. 3a, b we compare the cAFM images of the areas structured with 40,000 and 10,000âions per IID (corresponding AFM images are shown in Supplementary Fig. 3). For the higher amount of He ions in every location of an IID on the terraces a h-BN island can be found, while for fewer ions only in 86% of the IIDs a h-BN island nucleated (with most islands missing on BLG, see also Supplementary Fig. 3). The h-BN islands nucleating at IIDs written with 40,000âions per IID have a diameter of ~165ânm and a non-hexagonal shape (see Fig. 3c). Furthermore, the islands contain additional layers and a cluster of NPs in the center, indicating a similar formation as in the area studied before (see Fig. 2).
The morphology of the h-BN forming at IIDs written with 10,000âions is shown in Fig. 3d. Several mono-layer thick h-BN islands on BLG present a hexagonal shape, with well-defined edges and without wrinkles. Their mean diameter amounts to ~135ânm and they have a common alignment among each other (indicated by green dashed arrows). Furthermore, considering the well-known epitaxial relationship between EG and SiC (EG[1010] // SiC[1120]) and the orientation of the surface step edges (see Fig. 1d), the results suggest that the edges of the h-BN islands parallel to the green arrows follow the [1120] crystalline direction in EG. Unlike the non-hexagonal islands (see cAFM image in Fig. 3e), these oriented islands exhibit a Moiré pattern in their conductivity (Fig. 3f). The detailed cAFM image of an exemplary island (Fig. 3g) and the profile along the orange dotted line (Fig. 3h) reveal a periodicity of the Moiré pattern of ~14ânm. This value is expected for a relative rotation angle of 0° between graphene and h-BN, and indicates epitaxially aligned h-BN on graphene, with EG[1010]//h-BN[1010], which is considered thermally stable33,34. A similar periodicity can be observed in all hexagonal-shaped islands. Hence, single-crystalline and epitaxially aligned h-BN islands can be grown via SAvdWE by defect engineering in the EG substrate using a He FIB. It is assumed, that this approach is not only limited to h-BN on EG, but should be also applicable to other 2D materials, such as transition-metal dichalcogenides, which have been shown to preferentially nucleate at defects in 2D substrates as well15,35.
Furthermore, unlike in the non-hexagonal-shaped islands, the centers of the epitaxially aligned islands do not always coincide with the location of the IIDs. This could indicate that the defects written with 10,000 He ions per IID may heal out or move during heating of the substrate to the growth temperature36. In contrast, we expect that using 20,000 or 40,000 He ions per defect location resulted in the removal of multiple C atoms per individual IID and therefore in the formation of laterally extended defects (e.g. nanoholes), which are more stable upon thermal treatment and allow for differently oriented h-BN domains to nucleate at the same IID, resulting in the non-hexagonal-shaped h-BN islands. Further studies will elucidate the nature of the IIDs and develop protocols for their stabilization.
The NPs on the well-aligned islands shown in Fig. 3d also got partially removed during the cAFM measurements (see Supplementary Fig. 4), indicating their relatively weak bonding to the h-BN as compared to the clustered NPs observed in the islands nucleating at IIDs written with 20,000 and 40,000 He ions. This demonstrates that the deterministic placement of defects in the 2D substrate does not only affect the growth of the first 2D layer, but also the formation of additional layers.
Lateral growth of h-BN layer from IIDs
The vertical and lateral formation of the h-BN islands from the IIDs via SAvdWE was studied further by using an increased deposition rate of B and longer growth time to form a coalesced h-BN layer between the IIDs. The AFM image in Fig. 4a depicts an overgrown surface area in BLG, which was patterned with 20,000âions per defect location, in a distance of 250ânm to each other. Inside the patterned area a completely coalesced h-BN mono-layer formed due to the higher nucleation yield at the IIDs, while outside of the pattern BLG remains still largely uncovered (see Supplementary Fig. 5). The chemical bonds in this sample were further studied via X-ray photoemission spectroscopy (see Supplementary Fig. 6) with values for the B1s and N1s peaks in agreement with values for h-BN on EG14.
Apart from ~8ânm high NPs some additional h-BN layers with diameters of up to 100ânm formed close to some IIDs, which have a relatively lower tunneling current, as shown in the cAFM image in Fig. 4b (some NPs got removed during the cAFM scan). Furthermore, between some differently oriented domains we observe lower tunneling currents at the grain boundaries, where the morphology shows a network of wrinkles. This suggests the formation of overlapping h-BN layers upon the coalescence of the different domains in these areas, similar to observations in h-BN grown via catalytic chemical vapor deposition on metal foils37. In other areas we also observe the formation of atomically smooth grain boundaries (see Fig. 4c), where differently oriented grains merged seamlessly. No increased tunneling current could be detected from either of these grain boundaries, which is a key aspect if such h-BN layers are to be employed as tunneling barriers.
The h-BN grain on the right-hand side of the cAFM image in Fig. 4c presents a Moiré pattern with a periodicity of 14ânm, indicating the thermally stable relative rotation angle of 0° between graphene and h-BN (as also observed in Fig. 3g). A Moiré periodicity of 3ânm, as we observed it in the grain to the left-hand side of Fig. 4c, indicates a relative rotation angle of ~5°, which is thermally unstable34. We attribute the appearance of such unstable orientations in our sample to thermal expansion during cool down in combination with fixation of the h-BN domains at the IIDs.
In comparison to standard growth, another advantage of SAvdWE is the potential of controlling the grain boundary location by IID pattern design. Furthermore, neither the overlapping nor the seamless grain boundaries induced preferential formation of additional h-BN layers. Instead, additional layers can be found in the vicinity of the IIDs. This demonstrates that the artificially created defects can act as preferential nucleation points for the formation of additional h-BN layers on top of the first coalesced one. Thereby it appears plausible, that also other 2D materials could be connected to the IIDs upon completion of the first 2D layer, allowing more complex 2D heterostructures to form via SAvdWE. For the future device integration of such heterostructures, the active region could then be designed such, that it only incorporates the area in between the IIDs, where the heterostructures grow smoothly.
Figure 4d shows representative IâV curves from different locations on the samples with varying numbers of h-BN layers (the measurement configuration is depicted in the inset). The tunnel conductance shown in Fig. 4e was extracted from the linear regime of these IâV curves and the observed exponential decay with increasing layer number is expected for direct tunneling of electrons. A similar tunneling behavior was observed for h-BN flakes, exfoliated from high-quality h-BN bulk crystals, showing that the electronic quality of the h-BN grown via SAvdWE can compete with heterostructures fabricated via mechanical transfer38,39.
In conclusion, we demonstrated a mask-free selective-area epitaxy of two-dimensional h-BN on graphene using a deterministic placement of defects via He FIB. The SAvdWE allowed an unprecedented control over the formation of epitaxially aligned and single-crystalline 2D heterostructures via a bottom-up approach. Additionally to the nucleation location, the lateral formation of h-BN on EG was controlled via the growth parameter during the growth process, resulting in a 2D h-BN layer with tunneling characteristics identical to exfoliated material. Therefore, the SAvdWE contributes to precise and scalable fabrication of 2D heterostructures by establishing a reliable processing route for the site-selective, epitaxial bottom-up manufacturing with the potential for a reproducible fabrication of nano-devices based on 2D materials.
Methods
Substrate preparation
The markers for the alignment of the FIB pattern were prepared on the 1âÃâ1âcm² large 6H-SiC substrates with a miscut of ~0.2° using photolithography, followed by reactive ion etching using SF6. After chemically cleaning, the SiC(0001) was etched for 15âmin at 1450â°C in forming gas (5âat % H2 in Ar) using a flux of 500 sccm at a pressure of 900âmbar. In a second step EG formed on the etched SiC via sublimation of Si for 30âmin at 1650â°C in Ar at a pressure of 900âmbar. Further information on the etching and sublimation process can be found elsewhere40,41. For the thermal coupling during the growth of h-BN via MBE the backside of the SiC substrate was coated with 1âμm of Ti via electron beam evaporation.
He FIB patterning
The samples were loaded to the FIB system one day in advance to the actual patterning to allow for desorption of adhering species during pumping over night. The background pressure during patterning was around 1âÃâ10â7 mbar. After adjusting column and beam all stage moves and patterning tasks were carried out without imaging. Hence, unintended ion beam irradiation of the patterned regions could be avoided. The pattern was designed as described in the discussion of Fig. 1, with varying dwell times and distance between IIDs. Further information on the patterning and beam control can be found elsewhere42. The IIDs were prepared using a 30âkeV He ion beam at an ion current of 1.8 pA. With the settings for the FIB used in this study, the dwell times were varied between 4, 2, and 1âms corresponding to approximately 40,000, 20,000 and 10,000âions per intended defect location, respectively. The appropriate ion number for creating defects which initiate nucleation was determined on test samples, where the ion dose was varied over several orders of magnitude covering the range from one ion per IID up to 130,000âions. A typical test pattern for this purpose contained 6 individual IIDs as well as IIDs written along 6 lines, each 1âµm long, to ensure sufficient statistics (see Supplementary Fig. 7).
PAMBE of h-BN
After loading the patterned EG samples into the MBE system, they were thermally cleaned in a preparation chamber for 30âmin in UHV (1âÃâ10â9 mbar) at 400â°C. Consecutively, the h-BN islands (coalesced h-BN layers) were grown for 300âmin (480âmin) at a substrate temperature of 850â°C using a high temperature effusion cell operated at ~1850â°C (1900â°C) as a B source, whereas N was supplied using an Addon RFX450 RF plasma source operated at 350âW and a N2 flux of 0.2 sccm. The temperature was monitored using an optical pyrometer, and the heating and cooling rates were ~20â°C/min. Further information on the growth behavior of h-BN on EG at similar conditions can be found elsewhere12,14.
Micro-Raman spectroscopy
Raman measurements were conducted in back-scattering configuration using an excitation wavelength of 473ânm, with the laser focused on the sample using a Ã100 objective with a NA of 0.9, resulting in a spot size of ~1âμm. The reflected light was analyzed using a grating with 600âgr.âmmâ1 and a liquid N2 cooled Si CCD camera. Raman spectra were recorded for each pattern in a 5âÃâ5âμm2 map containing 36 measurement spots (each with 2âÃâ60âs integration time) on the patterned sample with and without h-BN.
AFM and cAFM
The morphology of the samples was monitored using AFM in intermittent contact mode (tip diameter, 2ânm; spring constant, 0.4âN/m) with the cantilever deflection used as feedback43. AFM scans were aligned with regard to the marker in the FIB pattern (e.g., see Fig. 1e and Supplementary Fig. 1). The cAFM measurements were conducted in contact mode at ambient conditions using Pt/Ir coated tips (tip diameter, 25ânm; spring constant, 0.4âN/m). A bias between 100 and 200âmV was applied to the EG on the semi-insulating SiC as described in the text, and the tunneling current between graphene and AFM tip was measured at a constant force of ~10ânN. IâV curves were recorded in N2 atmosphere at 144 locations in each of the areas shown in Fig. 4a, b with five repetitions per location.
X-ray photoemission spectroscopy
XPS measurements presented in Supplementary Fig. 6 were carried out with a Scienta Omicron Argus CU electron analyzer using a monochromated Al Ka X-ray source (Al Kα: 1486.6âeV), with a spot size of a few millimeters in diameter.
Data availability
The data supporting the findings of this work are available from the corresponding authors upon reasonable request.
References
Novoselov, K. S., Mishchenko, A., Carvalho, A. & Castro Neto, A. H. 2D materials and van der Waals heterostructures. Science 353, aac9439âaac9439 (2016).
Liu, C. et al. Two-dimensional materials for next-generation computing technologies. Nat. Nanotechnol. 15, 545â557 (2020).
Yankowitz, M., Ma, Q., Jarillo-Herrero, P. & LeRoy, B. J. van der Waals heterostructures combining graphene and hexagonal boron nitride. Nat. Rev. Phys. 1, 112â125 (2019).
Lyu, J., Pei, J., Guo, Y., Gong, J. & Li, H. A new opportunity for 2D van der Waals heterostructures: making steep-slope transistors. Adv. Mater. 32, 1906000 (2020).
Koma, A. Fabrication of ultrathin heterostructures with van der Waals epitaxy. J. Vac. Sci. Technol. B: Microelectron. Nanom. Struct. 3, 724 (1985).
Meng, J., Wang, D., Cheng, L., Gao, M. & Zhang, X. Recent progress in synthesis, properties, and applications of hexagonal boron nitride-based heterostructures. Nanotechnology 30, 074003 (2019).
Mayorov, A. S. et al. Micrometer-scale ballistic transport in encapsulated graphene at room temperature. Nano Lett. 11, 2396â2399 (2011).
Lopes, J. M. J. Synthesis of hexagonal boron nitride: from bulk crystals to atomically thin films. Prog. Cryst. Growth Charact. Mater. 67, 100522 (2021).
Liu, Z. et al. Direct growth of graphene/hexagonal boron nitride stacked layers. Nano Lett. 11, 2032â2037 (2011).
Wang, M. et al. A platform for large-scale graphene electronicsâCVD growth of single-layer graphene on CVD-grown hexagonal boron nitride. Adv. Mater. 25, 2746â2752 (2013).
Cho, Y. et al. Hexagonal boron nitride tunnel barriers grown on graphite by high temperature molecular beam epitaxy. Sci. Rep. 6, 34474 (2016).
Heilmann, M. et al. Influence of proximity to supporting substrate on van der waals epitaxy of atomically thin graphene/hexagonal boron nitride heterostructures. ACS Appl. Mater. Interfaces 12, 8897â8907 (2020).
Lin, Y.-C. et al. Direct synthesis of van der Waals solids. ACS Nano 8, 3715â3723 (2014).
Heilmann, M., Bashouti, M., Riechert, H. & Lopes, J. M. J. Defect mediated van der Waals epitaxy of hexagonal boron nitride on graphene. 2D Mater. 5, 025004 (2018).
Zhang, X. et al. Defect-controlled nucleation and orientation of WSe2 on hBN: a route to single-crystal epitaxial monolayers. ACS Nano 13, 3341â3352 (2019).
Briggs, N. et al. A roadmap for electronic grade 2D materials. 2D Mater. 6, 022001 (2019).
Li, H. et al. Controlled synthesis of topological insulator nanoplate arrays on mica. J. Am. Chem. Soc. 134, 6132â6135 (2012).
Lin, M. et al. Controlled growth of atomically thin In 2 Se3 flakes by van der Waals epitaxy. J. Am. Chem. Soc. 135, 13274â13277 (2013).
Zhou, Y. et al. Epitaxy and photoresponse of two-dimensional gase crystals on flexible transparent mica sheets. ACS Nano 8, 1485â1490 (2014).
Munshi, A. M. et al. Selective area growth of AlGaN nanopyramid arrays on graphene by metal-organic vapor phase epitaxy. Appl. Phys. Lett. 113, 263102 (2018).
Lehtinen, O. et al. Effects of ion bombardment on a two-dimensional target: Atomistic simulations of graphene irradiation. Phys. Rev. B: Condens. Matter Mater. Phys. 81, 153401 (2010).
Ghorbani-Asl, M., Kretschmer, S., Spearot, D. E. & Krasheninnikov, A. V. Two-dimensional MoS2 under ion irradiation: From controlled defect production to electronic structure engineering. 2D Mater. 4, 025078 (2017).
Kretschmer, S. et al. Supported two-dimensional materials under ion irradiation: the substrate governs defect production. ACS Appl. Mater. Interfaces 10, 30827â30836 (2018).
Klein, J. et al. Site-selectively generated photon emitters in monolayer MoS2 via local helium ion irradiation. Nat. Commun. 10, 2755 (2019).
Yun, J. et al. Selective-area heteroepitaxial growth of h-BN micropatterns on graphene layers. 2D Mater. 5, 015021 (2018).
Cancado, L. G. et al. Quantifying defects in graphene via Raman spectroscopy at different excitation energies. Nano Lett. 11, 3190â3196 (2011).
Schedin, F. et al. Detection of individual gas molecules adsorbed on graphene. Nat. Mater. 6, 652â655 (2007).
Sclauzero, G. & Pasquarello, A. Carbon rehybridization at the graphene/SiC(0001) interface: effect on stability and atomic-scale corrugation. Phys. Rev. B: Condens. Matter Mater. Phys. 85, 161405(R) (2012).
Bellucci, L., Cavallucci, T. & Tozzini, V. From the buffer layer to graphene on silicon carbide: exploring morphologies by computer modeling. Front. Mater. 6, 1â6 (2019).
Stehle, Y. et al. Synthesis of hexagonal boron nitride monolayer: control of nucleation and crystal morphology. Chem. Mater. 27, 8041â8047 (2015).
Reich, S. et al. Resonant Raman scattering in cubic and hexagonal boron nitride. Phys. Rev. B 71, 205201 (2005).
Wang, Z. et al. Composite super-moiré lattices in double-aligned graphene heterostructures. Sci. Adv. 5, eaay8897 (2019).
Yankowitz, M. et al. Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382â386 (2012).
Wang, D. et al. Thermally induced graphene rotation on hexagonal boron nitride. Phys. Rev. Lett. 116, 126101 (2016).
Zhang, F. et al. Full orientation control of epitaxial MoS2 on hBN assisted by substrate defects. Phys. Rev. B 99, 155430 (2019).
Li, D. W. et al. Controlled defect creation and removal in graphene and MoS 2 monolayers. Nanoscale 9, 8997â9008 (2017).
Bayer, B. C. et al. Introducing overlapping grain boundaries in chemical vapor deposited hexagonal boron nitride monolayer films. ACS Nano 11, 4521â4527 (2017).
Lee, G. H. et al. Electron tunneling through atomically flat and ultrathin hexagonal boron nitride. Appl. Phys. Lett. 99, 243114 (2011).
Britnell, L. et al. Electron tunneling through ultrathin boron nitride crystalline barriers. Nano Lett. 12, 1707â1710 (2012).
Oliveira, M. H., Schumann, T., Ramsteiner, M., Lopes, J. M. J. & Riechert, H. Influence of the silicon carbide surface morphology on the epitaxial graphene formation. Appl. Phys. Lett. 99, 111901 (2011).
Galves, L. A. et al. The effect of the SiC (0001) surface morphology on the growth of epitaxial mono-layer graphene nanoribbons. Carbon 115, 162â168 (2017).
Deinhart, V. et al. The patterning toolbox FIB-o-mat: exploiting the full potential of focused helium ions for nanofabrication. Beilstein J. Nanotechnol. 12, 304â318 (2021).
Kaemmar, S. B. Introduction to Brukerâs ScanAsyst and PeakForce Tapping AFM Technologyâ, Application Note 133, Rev. A0 (2011).
Acknowledgements
The authors would like to acknowledge fruitful discussions with Dr. Arkady Krasheninnikov on the irradiation-induced defect formation in 2D materials. M.H, A.T., and J.M.J.L. would also like to thank Werner Seidel for his help in the preparation of the patterned SiC substrates. Furthermore, they are grateful to Carsten Stemmler, Hans-Peter Schönherr, and Claudia Herrmann, for their dedicated maintenance of the PAMBE system and appreciate the critical reading of the manuscript by Dr. Alberto Hernández-MÃnguez. V.D. and K.H. acknowledge financial support from DFG under grant no. HO 5461/3-1 and thank the European COST action CA19140 âFIT4NANOâ (www.fit4nano.eu). The He ion beam patterning was performed in the Corelab Correlative Microscopy and Spectroscopy at Helmholtz-Zentrum Berlin.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
M.H., K.H., and J.M.J.L. conceived and designed the experiments. M.H. and A.T. designed the mask for the photolithography used for marking the SiC substrates. A.T. supervised the photolithography processes. M.H. fabricated the EG samples on the marked SiC substrates. M.H., V.D., and K.H designed the mask for the FIB structuring. V.D. and K.H. performed the FIB structuring of the defects and characterized the samples using HIM and SEM. M.H. performed the MBE growth of h-BN and characterized the samples using OM, AFM, TUNA, and Raman measurements. J.M.J.L. performed XPS measurements, supervised the project and advised on experiments and data analysis. M.H. wrote the manuscript. All authors discussed the results and commented on the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
41699_2021_250_MOESM1_ESM.pdf
Supplementary Information: Spatially controlled epitaxial growth of 2D heterostructures via defect-engineering using a focused He ion beam
Rights and permissions
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/.
About this article
Cite this article
Heilmann, M., Deinhart, V., Tahraoui, A. et al. Spatially controlled epitaxial growth of 2D heterostructures via defect engineering using a focused He ion beam. npj 2D Mater Appl 5, 70 (2021). https://doi.org/10.1038/s41699-021-00250-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41699-021-00250-z
This article is cited by
-
Understanding the impact of heavy ions and tailoring the optical properties of large-area monolayer WS2 using focused ion beam
npj 2D Materials and Applications (2023)
-
Pinhole-seeded lateral epitaxy and exfoliation of GaSb films on graphene-terminated surfaces
Nature Communications (2022)