# Fileset

[eaay8897.full.pdf](https://mdr.nims.go.jp/filesets/835c6a2e-a60f-48e4-9cff-38298cd4b749/download)

## Creator

Zihao Wang, Yi Bo Wang, J. Yin, E. Tóvári, Y. Yang, L. Lin, M. Holwill, J. Birkbeck, D. J. Perello, Shuigang Xu, J. Zultak, R. V. Gorbachev, A. V. Kretinin, [T. Taniguchi](https://orcid.org/0000-0002-1467-3105), [K. Watanabe](https://orcid.org/0000-0003-3701-8119), S. V. Morozov, M. Anđelković, S. P. Milovanović, L. Covaci, F. M. Peeters, A. Mishchenko, A. K. Geim, K. S. Novoselov, Vladimir I. Fal’ko, Angelika Knothe, C. R. Woods

## Rights

[Creative Commons BY Attribution 4.0 International](https://creativecommons.org/licenses/by/4.0/)

## Other metadata

[Composite super-moiré lattices in double-aligned graphene heterostructures](https://mdr.nims.go.jp/datasets/c1af8a5a-7de0-4219-b06d-aef3c89660a6)

## Fulltext

Science Journals — AAASWang et al., Sci. Adv. 2019; 5 : eaay8897     20 December 2019S C I E N C E  A D V A N C E S  |  R E S E A R C H  A R T I C L E1 of 7M A T E R I A L S  S C I E N C EComposite super-moiré lattices in double-aligned graphene heterostructuresZihao Wang1*, Yi Bo Wang1*, J. Yin1,2, E. Tóvári3, Y. Yang1,3, L. Lin3, M. Holwill3, J. Birkbeck3, D. J. Perello1,3, Shuigang Xu1,3, J. Zultak3, R. V. Gorbachev1,3,4, A. V. Kretinin3,5, T. Taniguchi6, K. Watanabe6, S. V. Morozov7, M. Anđelković8, S. P. Milovanović8, L. Covaci8, F. M. Peeters8, A. Mishchenko1,3, A. K. Geim1,3, K. S. Novoselov1,3,9,10, Vladimir I. Fal’ko1,3,4,  Angelika Knothe3†, C. R. Woods1,3†When two-dimensional (2D) atomic crystals are brought into close proximity to form a van der Waals hetero-structure, neighbouring crystals may influence each other’s properties. Of particular interest is when the two crystals closely match and a moiré pattern forms, resulting in modified electronic and excitonic spectra, crystal reconstruction, and more. Thus, moiré patterns are a viable tool for controlling the properties of 2D materials. However, the difference in periodicity of the two crystals limits the reconstruction and, thus, is a barrier to the low-energy regime. Here, we present a route to spectrum reconstruction at all energies. By using graphene which is aligned to two hexagonal boron nitride layers, one can make electrons scatter in the differential moiré pattern which results in spectral changes at arbitrarily low energies. Further, we demonstrate that the strength of this potential relies crucially on the atomic reconstruction of graphene within the differential moiré super cell.INTRODUCTIONVan der Waals heterostructures allow combining different two-dimensional (2D) materials into functional stacks (1, 2), which has already produced a range of interesting electronic (3, 4) and opto-electronic (5–8) devices and resulted in observation of exciting phys-ical phenomena. The large variety of the heterostructures is mainly due to the large selection of 2D materials. However, the assembly of van der Waals heterostructures allows one extra degree of freedom: Apart from the selection of the sequence of the 2D crystals, the indi-vidual crystals can be differently oriented with respect to each other. Previously, such control over the rotational alignment between crys-tals resulted in the observation of the resonant tunneling (9–11), renormalization of exciton binding energy (12), and insulating (13) and superconducting (4) states.Probably one of the most spectacular results of the rotational alignment between different 2D crystals is the observation of the band reconstruction due to electron scattering on the moiré pattern in graphene aligned with hexagonal boron nitride (hBN). Because the lattice constants of graphene and hBN are relatively close to each other, the alignment between the two crystals leads to the formation of a moiré pattern (14, 15) with a relatively small wave vector, which results in the appearance of the secondary Dirac points (16–18) in the electronic spectrum. Furthermore, the strong van der Waals interaction also leads to the atomic reconstruction of the graphene lattice (19–22). Unfortunately, the characteristic energies at which the electronic spectrum can be reconstructed are given by the difference between the lattice constants of graphene and hBN, which does not allow changes to be made to the low-energy part of the spectrum.Here, we demonstrate how we can gain further control over the band reconstruction of graphene using the differential between two moiré patterns (super-moiré) created by top and bottom hBN in hBN/graphene/hBN heterostructures. Such super-moiré patterns are not related to the difference in the lattice constants between the two crystals and thus can be of any arbitrary wave number, which makes it possible to arrange the spectrum reconstruction at arbitrary low energies.RESULTSTo this end, we created encapsulated graphene devices where the graphene layer was aligned to both bottom and top hBN layers (align-ment angles  and ). The fabrication and transfer procedures have been previously described in (23) with the exception that not only the bottom but also the top hBN is now crystallographically aligned to the graphene. Briefly, we started by identifying the top hBN layer on SiO2. We then use a thin film of polypropylene carbonate (PPC) on Polydimethylsiloxane (PDMS) to lift the hBN from its substrate. This film then facilitates bringing the top hBN into contact with a graphene crystal (Fig. 1, A and B). We use very long and straight edges of the crystals to identify crystallographic axes and align them using a commercially available transfer rig (24). The graphene can then be lifted away from its substrate (Fig. 1B). At this point, we per-form atomic force microscopy (AFM) (14, 15, 19) and Raman spec-troscopy (25) experiments on the hBN/graphene bilayer to confirm the alignment. One such AFM image is presented in Fig. 1E, showing the characteristic hexagonal pattern (the Fourier transformation 1Department of Physics and Astronomy, University of Manchester, Oxford Road, Manchester M13 9PL, UK. 2Institute of Nano Science, Nanjing University of Aero-nautics and Astronautics, Nanjing 210016, China. 3National Graphene Institute, University of Manchester, Oxford Road, Manchester M13 9PL, UK. 4Henry Royce Institute for Advanced Materials, Oxford Road, Manchester M13 9PL, UK. 5Department of Materials, University of Manchester, Oxford Road, Manchester M13 9PL, UK. 6National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan. 7Institute of Microelectronics Technology RAS, Chernogolovka 142432, Russia. 8Department of Physics, University of Antwerp, Groenenborgerlaan 171, Antwerp, Belgium. 9Centre for Advanced 2D Materials, National University of Singapore, Singapore 117546, Singapore. 10Chongqing 2D Materials Institute, Liangjiang New Area, Chongqing 400714, China.*These authors contributed equally to this work.†Corresponding author. Email: colin.woods@manchester.ac.uk (C.R.W.); angelika.knothe@manchester.ac.uk (A.K.)Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY). on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/Wang et al., Sci. Adv. 2019; 5 : eaay8897     20 December 2019S C I E N C E  A D V A N C E S  |  R E S E A R C H  A R T I C L E2 of 7is shown in Fig.  1H). The crystals are then aligned and brought into contact with a second thin hBN layer (typically less than 1.5 nm or five atomic layers thick; Fig. 1B). This layer is also lifted away from its substrate, leaving a triple layer on the thin polymer film (Fig. 1C). We then perform AFM and Raman characterization again. Figure 1F and Fig. 1I are an example of one of our double-aligned moiré AFM images and its Fourier transformation for the case of the second hBN layer being one atomic layer thick, which allows one to see both moiré patterns simultaneously. Although not immedi-ately clear in the real-space image, the Fourier transformation shows two sets of peaks corresponding to two hexagonal patterns (red and green dashed hexagons), as also schematically shown in Fig. 1D. The triple layer is then misaligned (~15°) and placed on top of a thick substrate hBN (Fig. 1C). Last, we use standard lithographic techniques to create the Hall bar geometry.The longitudinal resistance (Rxx) as a function of carrier concen-tration is shown in Fig. 2A. Here, apart from the resistance peak associated with the main Dirac point (26), several additional peaks can be seen. Most of such peaks correspond to the change of sign of the transversal (Hall) resistance (Rxy) measured in nonquantized magnetic field (Fig. 2B). Typically, if graphene is aligned with only one hBN, then a single moiré pattern is produced, and only one secondary Dirac point for electrons and one for holes can be seen at concentrations that correspond to the wave vector determined by the periodicity of the moiré pattern (14–18). Aligning graphene to both the top and the bottom hBN will produce two moiré patterns (if  and  are not equivalent), which should result in two second-ary Dirac points for electrons and two for holes. However, if elec-trons can feel potential from both moiré patterns simultaneously, then second-order processes can be allowed, which would result in the reconstruction of the electronic spectrum at many other wave vectors.In quantized magnetic fields, Landau fans can be seen to origi-nate from these peaks (Fig. 2C). The Landau fans for most peaks Fig. 1. Device fabrication and characterization. (A) Step 1: A thick hBN layer is aligned and used to pick up graphene (Gr). (B) Step 2: A thin hBN layer is aligned and picked up forming a triple layer. (C) Step 3: The heterostructure is placed on top of a thick hBN layer at 15° rotation angle (the substrate). (D) Illustration of the two indi-vidual moiré patterns and super-moiré pattern for three overlapping hexagonal lattices. (E) AFM image of the moiré pattern after the graphene is picked up. This shows only one moiré periodicity. Scale bar, 100 nm [(F) shares this scale]. (F) AFM image of the moiré patterns after the second thin hBN layer is picked up. Here, periodicities due to both hBN layers are visible. (G) Reciprocal space image of graphene’s first Brillouin zone. G1 (blue), ​​g​1​ ​​ (gray), and ​​g​1​  ​​ (purple) are the reciprocal lattice vectors for graphene and  and  hBN layers, respectively.  and  are at angles  and  relative to graphene. ​​b​1​  ​​  (green) is the moiré between graphene and the  hBN layer. ​​b​1​ ​​ (red) is the moiré between graphene and the  hBN layer. (H) Fourier transformation of the image in (E), displaying only one hexagonal periodic pattern (red dashed hexagon). Scale bar, 0.2 nm−1 [(I) shares this scale]. (I) Fourier transformation of the image in (F), showing two sets of distinct hexagonal patterns (red and green dashed hexagons). on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/Wang et al., Sci. Adv. 2019; 5 : eaay8897     20 December 2019S C I E N C E  A D V A N C E S  |  R E S E A R C H  A R T I C L E3 of 7exhibit both positive and negative indices (positive and negative slopes in Fig. 2C), which suggests that those are originating from the additional Dirac points and not from the higher-order zone edges (27).To interpret these additional peaks, we recall that the perfectly aligned graphene on hBN should produce moiré pattern with ap-proximately 14-nm periodicity, which corresponds to carrier concen-tration n ≈ 2.3 × 1012 cm−2. Thus, for our double-aligned graphene, we interpret the most pronounced peaks at ±1.98 × 1012 and ±2.34 × 1012 cm−2 as coming from the moiré patterns from the top and bot-tom hBN layers.The periodicity of a moiré pattern, L, can be related to the carrier concentration required to reach the edge of its first Brillouin zone by ​​n​ SDP​​  =  ​  8 _ ​√ _ 3 ​​L​​ 2​​​. Using this, we get periodicities of 15.3 and 14.0 nm, respectively, for the two most prominent features. The moiré peri-odicities are dependent on both the lattice constant mismatch and the alignment angle, as given by ​​b​​ ,​ = ∣​b​n​ ,​∣ = ​ 4 _ ​√ _ 3 ​a​ ​√ _ ​​​ 2​ + ​​​ ,​​​ 2​​ ​​, where ​​b​n​ ,​  =  ​G​ n​​ − ​g​n​ ,​​ (n = 1,…, 6) are the moiré reciprocal lattice vector between the  or  hBN layer (with reciprocal lattice vectors ​​g​n​ ,​​) and graphene (with reciprocal lattice vectors Gn),  is the graphene-hBN lattice constant mismatch, ,  is the misalignment angle for  or , and a is graphene’s lattice constant. One of the observed periods is larger than that could be expected for graphene aligned with hBN [~15.3 nm; see (14–19)]. We attribute this slightly larger moiré period to stretching of graphene as it interacts more strongly with the two aligned hBN layers. Since the angle is zero, or sufficiently close to zero ( ≫ ), we may calculate a new . In this case, the lattice mismatch to achieve the periodicity of 15.3 nm should be ~1.64%. This cor-responds to ~0.16% strain in the graphene crystal.Then, we would like to notice the small peaks at ±0.35 × 1012 cm−2, which corresponds to the largest differential moiré pattern. From the carrier concentration, we can infer a periodicity of 35 nm. Further still, there is a pronounced peak at approximately ±0.90 × 1012 cm−2, which would yield a period of 22.7 nm. These features represent pre-viously impossible periodicities for the graphene/hBN moiré pattern.In Fig. 3A, we schematically describe the geometric origin of the super-moiré features. ​​b​m​  ​​ and ​​b​k​ ​​ (red and green vectors; m = 1, 2,…, 6, k = 1, 2,…, 6) are the  and  moiré patterns. Their combination produces six new super-moiré patterns by the combinations of the vectors. In Fig. 3A, we highlight the ​​b​1​ ​ − ​b​k​ ​​ vectors (blue). In Fig. 3B, we present the position of the moiré and super-moiré zone edges in carrier concentration as a function of the angle between the second hBN layer () and graphene for the case when the first hBN layer is held at zero angle mismatch ( = 0) and = 1.64%, as calculated. For  = 0.4°, the features correspond exactly to the observed peaks in Rxx (as shown by the dashed lines connecting Fig. 3B and Fig. 3C) and sign reversal of Rxy, thus revealing the presence of new second-ary Dirac points in the low-energy electronic spectrum. Such low-energy peaks originate from the differential super-moirés, ​​b​1​ ​ − ​b​1​ ​​, ​​b​1​ ​ − ​b​6​ ​​, and ​​b​1​ ​ − ​b​2​ ​​. Further, in Fig. 3D, we show the position of the Rxx peak in carrier concentration for the ​​b​1​ ​ − ​b​2​ ​​ super-moiré, against the calculated angle between the two hBNs (∣ − ∣). This peak is unique because its period is independent of the graphene sheet as it is geometrically identical to moiré pattern between the two hBN layers. As expected, Fig. 3D shows that the peak position for four of our samples (blue circles) exactly follows expectation (red line).To check that all these peaks originate from the spectrum recon-struction because of scattering on the additional periodic potential, we measured the Brown-Zak oscillations at elevated temperatures (T) where cyclotron oscillations are suppressed (Fig. 4). At T > 70 K, oscillations independent of the carrier concentration can be clearly seen. At low fields B < 2.5 T (Fig. 4B), the oscillations are periodic in 1/B with the fundamental field Bf = 9.3 T. Assuming a hexagonal unit cell, such fundamental field can be calculated to correspond to a moiré periodicity of 22.7 nm, which corresponds to the peak in Rxx at n = ±0.90 × 1012 cm−2 (Fig. 2A).The behavior at high fields is more complex, as several Brown-Zak oscillations that originate from different periodicities overlap. However, by taking the periodicities, which correspond to the most prominent Fig. 2. Transport properties of double-aligned hBN/graphene/hBN device. (A) Rxx as a function of n for one of our devices with ​​b​1​ α​​ ≈ 15.3 nm (​​φ​​ α​​ = 0°), ​​b​1​ β ​​≈ 14.0 nm (​​φ​​ β​​ = 0.4°). Lattice mismatch, , is taken as 1.64%. The moiré and super-moiré peaks are marked by arrows and also labeled with their periods on the hole side (in nanometers). The position of the peaks is symmetric with respect to holes and elec-trons. The top axis mark is labeled in the size of the moiré pattern that corresponds to the particular carrier concentration. Dashed lines correspond to the arrows and can be traced to those in (B). (B) Rxy for the same device measured at B = 0.2 T, symmetrized to avoid contributions from Rxx. The dashed lines can be traced to the arrows in (A) and correspond to the particular peaks in Rxx. (C) Fan diagram xx(n, B) for the same device (scale black to white, 0.5 to 70 e2/hour). All measurements are performed at T = 1.7 K. on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/Wang et al., Sci. Adv. 2019; 5 : eaay8897     20 December 2019S C I E N C E  A D V A N C E S  |  R E S E A R C H  A R T I C L E4 of 7peaks in Rxx at B = 0 (15.3, 14, and 11.2 nm; see Fig. 2A), we could identify most oscillations in terms of fractions of the flux quantum per the corresponding plaquette, labeled in Fig. 4C. Thus, the graphene-hBN moiré periods (14.0 and 15.3 nm) and super-moiré periods (11.2, 22.7, and 36.3 nm) each give features attributable to secondary Dirac points at well-understood values of carrier concen-tration. In addition, they produce clear Brown-Zak oscillations for unitary flux through moiré unit cells for the 11.2-, 14.0-, 15.3-, and 22.7-nm periods.In addition to our previous observations, we would like to note that there are several unexplained features in Rxx (Fig. 2A) and Rxy (Fig. 2B), most pronounced at ne ≈±3.2 × 1012 and ±4.1 × 1012 cm−2. One possible explanation for these features is higher-order moiré periodicities, that is, moiré patterns between super-moiré periods. However, the probability of such multiple scattering events dimin-ishes strongly. Likewise, there could exist features due to super-moiré patterns between further zone edges (second, third,… Brillouin zone edges of the graphene/hBN moirés) or a more exotic superlattice (SL) phenomenon.DISCUSSIONTheoretically (28), moiré effects on graphene can be described in terms of a periodic SL potential applied to Dirac electrons produced by incommensurable lattices of two (top and bottom) hBN flakes​​ ​​​ ​​   H ​  =  vp • 𝛔 + ​Σ​ j=±​​ ​Σ​ n=0…5​​​[​​ ​U​0​ j ​ + ​​(​​ − 1​)​​​​ n​​(​​ ​iU​3​ j ​ ​σ​ 3​​ + ​U ​1​ j ​ ​​a​ n​​ • σ _ a  ​​)​​​]​​ ​   × e​​ i​b​n​ j ​•​(​​r+j​R _ 2 ​​)​​​ ​e​​ i​G​ n ​​• u​(​​r, R​)​​​​​  (1)Here, 3 and  = (1, 2) are Pauli matrices acting in the sublat-tice space of graphene’s Bloch states; j = ± identifies layers  (+) and  (−); ​​U​0​ j ​​, ​​U​1​ j ​​, and ​​U​3​ j ​​ parameterize a smoothly varying moiré poten-tial, the asymmetric sublattice on-site energies, and hopping between A and B sublattices, respectively [based on the earlier studies (29–31), ​​U ​0​ j ​  ≈  8.5 meV​, ​​U ​1​ j ​  ≈  − 17 meV​, and ​​U ​3​ j ​  ≈  − 15 meV​ for  j ≪ ]. Vec-tor R describes the phase shift between moiré produced by hBN flakes  and .From this, the super-moiré periods for the individual moiré SLs (, ) in graphene originate in two ways. One is due to the quantum mechanical interference, which appears in the second-order per-turbation theory. In this case, Eq. 1 allows for the electron scatter-ing from the combined j = ±SLs with the Bragg vectors ​​b​m​  ​ − ​b​k​ ​​. These are composed of different moiré SL reciprocal vectors, leading effectively to the SLs with Fourier components described in Figs. 1G and 3A. The second possibility is the reconstruction of graphene, which leads to a displacement field, u(r, R), generating mixing of the moiré SL’s reciprocal vectors. When this is considered, the longest-period super-moiré SL (m = 1, k = 1) that determines the Fig. 3. Super-moiré geometry. (A) Reciprocal-space image of the area around graphene’s K point. ​​b​m​  ​​ (green) and ​​b​m​  ​​ (red) for m = 1, 2,…, 6 are the graphene-hBN moiré reciprocal lattice vectors. ​​b​1​ ​ − ​b​m​  ​​ (blue) are the six super-moiré reciprocal lattice vectors. The blue hexagonal area indicates the ​​b​1​ ​ − ​b​1​  ​​ first Brillouin zone. (B) Carrier concentration of the first Brillouin zone edge for the two moiré and four lowest-energy super-moiré features as a function of ( = 1.64%,  = 0°). (C) Rxx peak positions in carrier concentration. Dashed lines connect values of carrier concentration for  = 0.4 in (B) to the position in (C). Each line matches a peak. (D) Carrier concentration of the ​​b​1​ ​ − ​b​1​  ​​ super-moiré feature versus ∣ − ∣  for four of our samples (blue circles) and by calculation (red line).Fig. 4. Brown-Zak oscillations in one of our double-aligned hBN/graphene/hBN devices. (A) Map of xx(n, B); scale of blue to red is 0.5 to 70 e2/hour. (B) Zoom-in view of the low-field part of the map, marked by the yellow dashed rectangle in (A); scale of black to white is 7 to 37 e2/hour. The Brown-Zak oscillations correspond to a moiré structure with a periodicity of 22.7 nm, and the fundamental field is 9.3 T. (C) Zoom-in view of the high-field part of the map, marked by the white dashed rectangle in (A); scale of blue to red is 7 to 37 e2/hour. The Brown-Zak oscillations that correspond to moiré structures of different periodicities are marked by dotted lines of different colors. Black, 15.3 nm (BF = 20.5 T); green, 14.0 nm (BF = 24.2 T); brown, 11.2 nm (BF = 38 T). All measurements are performed at T = 70 K. on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/Wang et al., Sci. Adv. 2019; 5 : eaay8897     20 December 2019S C I E N C E  A D V A N C E S  |  R E S E A R C H  A R T I C L E5 of 7low-energy part of graphene spectrum can be described by the fol-lowing SL potential ​​​​    H ​​1,1​ P  ​  ≈  − 12 ​U​ 3​​ ​ω​ as​​ ​σ​ 3​​ − ​Σ​ m​​​[​​2 ​U​ 3​​ ​ω​ as​​ ​σ​ 3​​ + ​ 4 ​U​ 3​​ ​U​ 1​​ ─ vb ​  + i ​ 2 ​U​1​ 2​ ─ vb ​ ​  ​b​ n​​ • 𝛔 ─ b  ​​]​​ × ​e​​ i​b​m​ α ​R​ ​e​​ i​(​​​b​m​ α ​−​b​m​ β ​​) ​​• r​​​​​​​   H ​​1,1​ AP ​  ≈ ​ ∑ m​ ​​​[​​i ​​(​​ − 1​)​​​​ m​ 2 ​U​ 0​​ ​ω​ as​​ − ​4 ​U​ 3​​ ​U​ 1​​ _ vb ​  + i ​2 ​U​1​ 2​ _ vb ​ ​​b​ n​​ • 𝛔 _ b  ​​]​​ × ​e​​ i​b​m​ α ​R​ ​e​​ i​(​​​b​m​ α ​−​b​m​ β ​​) ​​• r​​​     (2)This expression was derived for  ≪  for both parallel (P) and antiparallel (AP) mutual orientations of the two hBN crystals (hence, approximately ​ ​b​m​  ​ − ​b​m​  ​  ⊥ ​ b​m​ ,​​), and, here, as parameterizes the am-plitude of inversion asymmetric component of strain.Vital to this description is the understanding that the displace-ment field u(r, R) develops because of the competition between stacking-dependent van der Waals adhesion, of graphene and hBN, and elasticity of graphene. Previous work (19, 20) has identified that the crystals form a 2D fixed-density commensurate state when graphene and hBN are close to perfect alignment. The effect relies upon strain and, thus, u(r, R), modulating with a period matching that of the moiré pattern to minimize adhesive and elastic energy. The commensurate state is characterized by hexagonal domains with increasingly sharp domain walls near  = 0, observed in PeakForce AFM (32), and broadening of the 2D peak in the Raman spectrum (19, 25).To evaluate the degree of strain within our super-moiré samples, we have used Raman spectroscopy. In Fig. 5A, we show the 2D peak and its full width at half maximum (FWHM) for a typical unaligned sample (black circles), a sample aligned to one hBN (blue triangles), and one of our double-aligned samples (red squares). In the case that the two hBN layers may be treated entirely independently, the signature in the Raman spectrum would remain unchanged from the single-aligned case [FWHM (2D), ~36 cm−1]. However, it is clear that when a second aligned hBN layer is added, the width of the 2D peak increases by a factor of ~2 (Fig. 5A). We attribute this to re-structuring of strain within the super-moiré unit cells. This observa-tion supports the proposed model of the two moiré patterns mixing through strain fields. Therefore, u(r, R) should have periodicities de-scribed by ​​b​m​  ​​ (-moiré), ​​b​k​ ​​ (-moiré), and ​​b​m​  ​ − ​b​k​ ​​ (super-moirés). This is also supported by molecular dynamics simulations of the re-laxation in the double-aligned systems, as shown in Fig. 5B, where the Raman spectra of simulated relaxed configurations is presented (see the Supplementary Materials).CONCLUSIONTo conclude, graphene’s electronic spectrum is substantially altered by scattering from super-moiré structures described by the preexisting moiré between graphene and its substrate and encapsulating hBN layers. These alterations may be considered in two ways: as double-scattering events from both graphene-hBN moiré patterns or as single-scattering events from a reconstructed graphene layer. Such super-moiré poten-tial can be of an arbitrarily small wave vector (unlike moiré potential from single hBN aligned with graphene), which allows modification of the graphene band structure at arbitrarily low energies.MATERIALS AND METHODSFabricationIn addition to the sample preparation described previously, the gen-eral process is as follows. We used PPC, or sometimes polymethyl methacrylate, spun on a thick PDMS membrane to facilitate moving and orientating the crystals. This membrane was used to pick up the first hBN layer. The crystal was then positioned and aligned to graphene before the two were brought into contact. We removed the membrane quickly which lifted the graphene off from its substrate. At this point, we inverted the membrane and perform various characterization techniques on the half-assembled hetero-structure. We then picked up a further thin hBN crystal, before repeating our characterization. These previous steps could have been repeated multiple times to produce increasingly complex hetero-structures with a variety of crystals (although, here, we limited ourselves to graphene and hBN). Last, the stack of crystals was positioned and brought into contact with a final “substrate” hBN. The membrane was removed slowly so that all the crystals were left on a SiO2 wafer.Atomic force microscopySL-resolution AFM images of the moiré patterns were taken in PeakForce Quantitative Nanomechanical Mapping (PF-QNM) mode on a Bruker Icon AFM. PF-QNM allows the capture of individual tip-sample force curves for each pixel in the image. These force curves were used to extract additional elastic information about the Fig. 5. Strain distribution in the aligned graphene-hBN heterostructures. Raman spectra (2D peak region) for an unaligned sample (gray), single-aligned sample (blue), and double-aligned sample (red). (A) Experimental results. (B) Molecular dynamics relaxation simulations. a.u., arbitrary units. on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/Wang et al., Sci. Adv. 2019; 5 : eaay8897     20 December 2019S C I E N C E  A D V A N C E S  |  R E S E A R C H  A R T I C L E6 of 7sample as well as topographic information. We used “ScanAsyst Fluid+” probes whose resolution is regularly better than 2 nm.Raman spectroscopyRaman spectroscopy measurements were performed on a HORIBA XploRA PLUS Raman spectrometer. The laser wavelength was 532 nm with a power of 0.5 mW through a ×100 objective.Moiré Fourier analysisThe analysis of moiré patterns in double-aligned heterostructures was performed following an approach described in (33). The method-ology was based on the Fourier analysis of moiré patterns appearing when two hexagonal lattices are combined.Molecular dynamics simulationsMolecular dynamics simulations were performed for the single-aligned hBN/graphene and the double-aligned hBN/graphene/hBN by allowing the relaxations of both hBN and graphene layers. We used the bond-order Brenner potentials for the graphene layer, Tersoff potentials for the B-N interaction in the hBN layers, and the Morse potential developed in (34) for the interlayer interactions. The simulations are performed within the “large-scale atomic/molecular massively parallel simulator” (35, 36) by considering a disk of radius 120 nm. We fixed the atoms in a boundary region of 2 nm but allowed the relaxation of all other atoms. The total energy was minimized until the forces are below 10−6 eV/Å.Raman spectrum simulationUsing the relaxed structures, we calculated the shift of the 2D peak and the entire Raman spectrum of the graphene layer. Calculations were performed using the prescriptions given in (37–39).SUPPLEMENTARY MATERIALSSupplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/12/eaay8897/DC1More examples of double alignmentDifferent fundamental frequencies of Brown-Zak oscillationsGap opening at the main Dirac pointAFM of other double-aligned samplesUniformity in heterostructuresAnalysis of super-moiré peaksTight-binding modelMolecular dynamics simulations and Raman shift calculationsTable S1.  and  for each device.Fig. S1. Transport properties of double-aligned encapsulated graphene devices.Fig. S2. Brown-Zak oscillations in sample 1.Fig. S3. Brown-Zak oscillations for sample 4.Fig. S4. Electron-hole symmetry in super-moiré features.Fig. S5. Gap opening in one of our double-aligned samples.Fig. S6. Examples of double-aligned samples.Fig. S7. Uniformity in double-aligned heterostructures.Fig. S8. Frequency analysis of different harmonics of hBN-graphene-hBN structure.Fig. S9. Super-moiré periods corresponding to different harmonics.Fig. S10. Molecular dynamics simulations of bond lengths in graphene-hBN superlattices.Reference (40)REFERENCES AND NOTES  1.  A. K. Geim, I. V. Grigorieva, Van der Waals heterostructures. Nature 499, 419–425 (2013).  2.  K. S. Novoselov, A. Mishchenko, A. Carvalho, A. H. Castro Neto, 2D materials and van der Waals heterostructures. Science 353, aac9439 (2016).  3.  L. Britnell, R. V. Gorbachev, R. Jalil, B. D. Belle, F. Schedin, A. Mishchenko, T. Georgiou, M. I. Katsnelson, L. Eaves, S. V. Morozov, N. M. R. Peres, J. Leist, A. K. Geim, K. S. Novoselov, L. A. Ponomarenko, Field-effect tunneling transistor based on vertical graphene heterostructures. Science 335, 947–950 (2012).  4.  Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, P. Jarillo-Herrero, Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).  5.  L. Britnell, R. M. Ribeiro, A. Eckmann, R. Jalil, B. D. Belle, A. Mishchenko, Y. J. Kim, R. V. Gorbachev, T. Georgiou, S. V. Morozov, A. N. Grigorenko, A. K. Geim, C. Casiraghi, A. H. C. Neto, K. S. Novoselov, Strong light-matter interactions in heterostructures of atomically thin films. Science 340, 1311–1314 (2013).  6.  F. Withers, O. Del Pozo-Zamudio, A. Mishchenko, A. P. Rooney, A. Gholinia, K. Watanabe, T. Taniguchi, S. J. Haigh, A. K. Geim, A. I. Tartakovskii, K. S. Novoselov, Light-emitting diodes by band-structure engineering in van der Waals heterostructures. Nat. Mater. 14, 301–306 (2015).  7.  F. Withers, O. Del Pozo-Zamudio, S. Schwarz, S. Dufferwiel, P. M. Walker, T. Godde, A. P. Rooney, A. Gholinia, C. R. Woods, P. Blake, S. J. Haigh, K. Watanabe, T. Taniguchi, I. L. Aleiner, A. K. Geim, V. I. Fal’ko, A. I. Tartakovskii, K. S. Novoselov, WSe2 light-emitting tunneling transistors with enhanced brightness at room temperature. Nano Lett. 15, 8223–8228 (2015).  8.  T. Low, T. Low, A. Chaves, J. D. Caldwell, A. Kumar, N. X. Fang, P. Avouris, T. F. Heinz, F. Guinea, L. Martin-Moreno, F. Koppens, Polaritons in layered two-dimensional materials. Nat. Mater. 16, 182–194 (2016).  9.  T. L. M. Lane, J. R. Wallbank, V. I. Fal’ko, Twist-controlled resonant tunnelling between monolayer and bilayer graphene. Appl. Phys. Lett. 107, 203506 (2015).  10.  J. R. Wallbank, D. Ghazaryan, A. Misra, Y. Cao, J. S. Tu, B. A. Piot, M. Potemski, S. Pezzini, S. Wiedmann, U. Zeitler, T. L. M. Lane, S. V. Morozov, M. T. Greenaway, L. Eaves, A. K. Geim, V. I. Fal’ko, K. S. Novoselov, A. Mishchenko, Tuning the valley and chiral quantum state of Dirac electrons in van der Waals heterostructures. Science 353, 575–579 (2016).  11.  A. Mishchenko, J. S. Tu, Y. Cao, R. V. Gorbachev, J. R. Wallbank, M. T. Greenaway, V. E. Morozov, S. V. Morozov, M. J. Zhu, S. L. Wong, F. Withers, C. R. Woods, Y.-J. Kim, K. Watanabe, T. Taniguchi, E. E. Vdovin, O. Makarovsky, T. M. Fromhold, V. I. Fal’ko, A. K. Geim, L. Eaves, K. S. Novoselov, Twist-controlled resonant tunnelling in graphene/boron nitride/graphene heterostructures. Nat. Nanotechnol. 9, 808–813 (2014).  12.  E. M. Alexeev, D. A. Ruiz-Tijerina, M. Danovich, M. J. Hamer, D. J. Terry, P. K. Nayak, S. Ahn, S. Pak, J. Lee, J. I. Sohn, M. R. Molas, M. Koperski, K. Watanabe, T. Taniguchi, K. S. Novoselov, R. V. Gorbachev, H. S. Shin, V. I. Fal’ko, A. I. Tartakovskii, Resonantly hybridized excitons in moiré superlattices in van der Waals heterostructures. Nature 567, 81–86 (2019).  13.  Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, R. C. Ashoori, P. Jarillo-Herrero, Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80–84 (2018).  14.  J. Xue, J. Sanchez-Yamagishi, D. Bulmash, P. Jacquod, A. Deshpande, K. Watanabe, T. Taniguchi, P. Jarillo-Herrero, B. J. LeRoy, Scanning tunnelling microscopy and spectroscopy of ultra-flat graphene on hexagonal boron nitride. Nat. Mater. 10, 282–285 (2011).  15.  M. Yankowitz, J. Xue, D. Cormode, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, P. Jarillo-Herrero, P. Jacquod, B. J. LeRoy, Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382–386 (2012).  16.  L. A. Ponomarenko, R. V. Gorbachev, G. L. Yu, D. C. Elias, R. Jalil, A. A. Patel, A. Mishchenko, A. S. Mayorov, C. R. Woods, J. R. Wallbank, M. Mucha-Kruczynski, B. A. Piot, M. Potemski, I. V. Grigorieva, K. S. Novoselov, F. Guinea, V. I. Fal’ko, A. K. Geim, Cloning of Dirac fermions in graphene superlattices. Nature 497, 594–597 (2013).  17.  C. R. Dean, L. Wang, P. Maher, C. Forsythe, F. Ghahari, Y. Gao, J. Katoch, M. Ishigami, P. Moon, M. Koshino, T. Taniguchi, K. Watanabe, K. L. Shepard, J. Hone, P. Kim, Hofstadter’s butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598–602 (2013).  18.  B. Hunt, J. D. Sanchez-Yamagishi, A. F. Young, M. Yankowitz, B. J. LeRoy, K. Watanabe, T. Taniguchi, P. Moon, M. Koshino, P. Jarillo-Herrero, R. C. Ashoori, Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427–1430 (2013).  19.  C. R. Woods, L. Britnell, A. Eckmann, R. S. Ma, J. C. Lu, H. M. Guo, X. Lin, G. L. Yu, Y. Cao, R. V. Gorbachev, A. V. Kretinin, J. Park, L. A. Ponomarenko, M. I. Katsnelson, Y. N. Gornostyrev, K. Watanabe, T. Taniguchi, C. Casiraghi, H.-J. Gao, A. K. Geim, K. S. Novoselov, Commensurate–incommensurate transition in graphene on hexagonal boron nitride. Nat. Phys. 10, 451–456 (2014).  20.  C. R. Woods, F. Withers, M. J. Zhu, Y. Cao, G. Yu, A. Kozikov, M. Ben Shalom, S. V. Morozov, M. M. van Wijk, A. Fasolino, M. I. Katsnelson, K. Watanabe, T. Taniguchi, A. K. Geim, A. Mishchenko, K. S. Novoselov, Macroscopic self-reorientation of interacting two-dimensional crystals. Nat. Commun. 7, 10800 (2016). 21.  M. Yankowitz, K. Watanabe, T. Taniguchi, P. San-Jose, B. J. LeRoy, Pressure-induced commensurate stacking of graphene on boron nitride. Nat. Commun. 7, 13168 (2016). on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/cgi/content/full/5/12/eaay8897/DC1http://advances.sciencemag.org/cgi/content/full/5/12/eaay8897/DC1http://advances.sciencemag.org/Wang et al., Sci. Adv. 2019; 5 : eaay8897     20 December 2019S C I E N C E  A D V A N C E S  |  R E S E A R C H  A R T I C L E7 of 7  22.  M. Yankowitz, J. Jung, E. Laksono, N. Leconte, B. L. Chittari, K. Watanabe, T. Taniguchi, S. Adam, D. Graf, C. R. Dean, Dynamic band-structure tuning of graphene moiré superlattices with pressure. Nature 557, 404–408 (2018).  23.  T. Uwanno, Y. Hattori, T. Taniguchi, K. Watanabe, K. Nagashio, Fully dry PMMA transfer of graphene on h-BN using a heating/cooling system. 2D Mater. 2, 041002 (2015).  24.  Y. Cao, A. Mishchenko, G. L. Yu, E. Khestanova, A. P. Rooney, E. Prestat, A. V. Kretinin, P. Blake, M. B. Shalom, C. Woods, J. Chapman, G. Balakrishnan, I. V. Grigorieva, K. S. Novoselov, B. A. Piot, M. Potemski, K. Watanabe, T. Taniguchi, S. J. Haigh, A. K. Geim, R. V. Gorbachev, Quality heterostructures from two-dimensional crystals unstable in air by their assembly in inert atmosphere. Nano Lett. 15, 4914–4921 (2015).  25.  A. Eckmann, J. Park, H. Yang, D. Elias, A. S. Mayorov, G. Yu, R. Jalil, K. S. Novoselov, R. V. Gorbachev, M. Lazzeri, A. K. Geim, C. Casiraghi, Raman fingerprint of aligned graphene/h-BN superlattices. Nano Lett. 13, 5242–5246 (2013).  26.  K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, A. A. Firsov, Electric field effect in atomically thin carbon films. Science 306, 666–669 (2004).  27.  R. K. Kumar, X. Chen, G. H. Auton, A. Mishchenko, D. A. Bandurin, S. V. Morozov, Y. Cao, E. Khestanova, M. Ben Shalom, A. V. Kretinin, K. S. Novoselov, L. Eaves, I. V. Grigorieva, L. A. Ponomarenko, V. I. Fal’ko, A. K. Geim, High-temperature quantum oscillations caused by recurring Bloch states in graphene superlattices. Science 357, 181–184 (2017).  28.  J. R. Wallbank, A. A. Patel, M. Mucha-Kruczynski, A. K. Geim, V. I. Falko, Generic miniband structure of graphene on a hexagonal substrate. Phys. Rev. B 87, 245408 (2013).  29.  M. Lee, J. R. Wallbank, P. Gallagher, K. Watanabe, T. Taniguchi, V. I. Falko, D. Goldhaber-Gordon, Ballistic miniband conduction in a graphene superlattice. Science 353, 1526–1529 (2016).  30.  R. K. Kumar, A. Mishchenko, X. Chen, S. Pezzini, G. H. Auton, L. A. Ponomarenko, U. Zeitler, L. Eaves, V. I. Fal’ko, A. K. Geim, High-order fractal states in graphene superlattices.  Proc. Natl. Acad. Sci. U.S.A. 115, 5135–5139 (2018).  31.  J. R. Wallbank, R. Krishna Kumar, M. Holwill, Z. Wang, G. H. Auton, J. Birkbeck, A. Mishchenko, L. A. Ponomarenko, K. Watanabe, T. Taniguchi, K. S. Novoselov, I. L. Aleiner, A. K. Geim, V. I. Fal’ko, Excess resistivity in graphene superlattices caused by umklapp electron-electron scattering. Nat. Phys. 15, 32–36 (2019).  32.  B. B. Pittenger, N. Erina, C. Su, in Application Note #128 (2012).  33.  P. Zeller, S. Günther, What are the possible moiré patterns of graphene on hexagonally packed surfaces? Universal solution for hexagonal coincidence lattices, derived by a geometric construction. New J. Phys. 16, 083028 (2014).  34.  G. Argentero, A. Mittelberger, M. Reza Ahmadpour Monazam, Y. Cao, T. J. Pennycook, C. Mangler, C. Kramberger, J. Kotakoski, A. K. Geim, J. C. Meyer, Unraveling the 3D atomic structure of a suspended graphene/hBN van der Waals heterostructure. Nano Lett. 17, 1409–1416 (2017).  35.  S. Plimpton, Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).  36.  S. J. Plimpton, A. P. Thompson, Computational aspects of many-body potentials.  MRS Bull. 37, 513–521 (2012).  37.  S. Reich, J. Maultzsch, C. Thomsen, P. Ordejón, Tight-binding description of graphene. Phys. Rev. B 66, 035412 (2002).  38.  T. Vincent, V. Panchal, T. Booth, S. R. Power, A.-P. Jauho, V. Antonov, O. Kazakova, Probing the nanoscale origin of strain and doping in graphene-hBN heterostructures. 2D Mater. 6, 015022 (2019).  39.  T. M. G. Mohiuddin, A. Lombardo, R. R. Nair, A. Bonetti, G. Savini, R. Jalil, N. Bonini, D. M. Basko, C. Galiotis, N. Marzari, K. S. Novoselov, A. K. Geim, A. C. Ferrari, Uniaxial strain in graphene by Raman spectroscopy: G peak splitting, Grüneisen parameters, and sample orientation. Phys. Rev. B 79, 205433 (2009).  40.  D. A. Cosma, J. R. Wallbank, V. Cheianov, V. I. Fal’ko, Moiré pattern as a magnifying glass for strain and dislocations in van der Waals heterostructures. Faraday Discuss. 173, 137–143 (2014).Acknowledgments Funding: We acknowledge support from the EU Graphene Flagship Project (contract CNECTICT-604391); the European Research Council Synergy Grant Hetero2D; the Royal Society; EPSRC grants EP/S019367/1, EP/S030719/1, and EP/N010345/1; the U.S. Army Research Office (W911NF-16-1-0279); and the FLAG-ERA project TRANS2DTMD. S.V.M. was supported by RFBR (17-02-01129a). J.Y. and A.M. acknowledge the support of EPSRC Early Career Fellowship EP/N007131/1. S.P.M. was supported by the Flemish Science Foundation (FWO). E.T. acknowledges the support of the Marie Skłodowska–Curie grant agreement  no. 751883. Author contributions: K.S.N. and C.R.W. developed the concept of the experiment. T.T. and K.W. created hBN. Y.B.W., Y.Y., L.L., M.H., J.B., D.J.P., S.X., J.Z., and C.R.W. performed device fabrication. Y.B.W., Y.Y., L.L., M.H., and C.R.W. carried out Raman spectroscopy and AFM measurements. Z.W., J.Y., E.T., and S.V.M. performed transport measurements. Z.W., K.S.N., and C.R.W. performed experimental analysis. Z.W., K.S.N., V.I.F., A.K., and C.R.W. interpreted results with help from R.V.G., A.V.K., A.M., and A.K.G. M.A., S.P.M., L.C., F.M.P., V.I.F., and A.K. performed theoretical calculations. M.A., S.P.M., L.C., and F.M.P. performed molecular dynamics simulations. Z.W., K.S.N., V.I.F., A.K., and C.R.W. wrote the manuscript and the Supplementary Materials with additional contributions from all authors. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.Submitted 26 July 2019Accepted 22 October 2019Published 20 December 201910.1126/sciadv.aay8897Citation: Z. Wang, Y. B. Wang, J. Yin, E. Tóvári, Y. Yang, L. Lin, M. Holwill, J. Birkbeck, D. J. Perello, S. Xu, J. Zultak, R. V. Gorbachev, A. V. Kretinin, T. Taniguchi, K. Watanabe, S. V. Morozov, M. Anđelković, S. P. Milovanović, L. Covaci, F. M. Peeters, A. Mishchenko, A. K. Geim, K. S. Novoselov, V. I. Fal'ko, A. Knothe, C. R. Woods, Composite super-moiré lattices in double-aligned graphene heterostructures. Sci. Adv. 5, eaay8897 (2019). on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/Composite super-moiré lattices in double-aligned graphene heterostructuresPeeters, A. Mishchenko, A. K. Geim, K. S. Novoselov, Vladimir I. Fal'ko, Angelika Knothe and C. R. WoodsGorbachev, A. V. Kretinin, T. Taniguchi, K. Watanabe, S. V. Morozov, M. Andelkovic, S. P. Milovanovic, L. Covaci, F. M. Zihao Wang, Yi Bo Wang, J. Yin, E. Tóvári, Y. Yang, L. Lin, M. Holwill, J. Birkbeck, D. J. Perello, Shuigang Xu, J. Zultak, R. V.DOI: 10.1126/sciadv.aay8897 (12), eaay8897.5Sci Adv ARTICLE TOOLS http://advances.sciencemag.org/content/5/12/eaay8897MATERIALSSUPPLEMENTARY http://advances.sciencemag.org/content/suppl/2019/12/16/5.12.eaay8897.DC1REFERENCEShttp://advances.sciencemag.org/content/5/12/eaay8897#BIBLThis article cites 39 articles, 9 of which you can access for freePERMISSIONS http://www.sciencemag.org/help/reprints-and-permissionsTerms of ServiceUse of this article is subject to the  is a registered trademark of AAAS.Science AdvancesYork Avenue NW, Washington, DC 20005. The title (ISSN 2375-2548) is published by the American Association for the Advancement of Science, 1200 NewScience Advances BY).Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of on December 26, 2019http://advances.sciencemag.org/Downloaded from http://advances.sciencemag.org/content/5/12/eaay8897http://advances.sciencemag.org/content/suppl/2019/12/16/5.12.eaay8897.DC1http://advances.sciencemag.org/content/5/12/eaay8897#BIBLhttp://www.sciencemag.org/help/reprints-and-permissionshttp://www.sciencemag.org/about/terms-servicehttp://advances.sciencemag.org/