# Fileset

[sciadv.adk1874.pdf](https://mdr.nims.go.jp/filesets/70fd193e-ef04-47dc-9a45-6ab4a98d58aa/download)

## Creator

Yichao Zhang, Ji-Hwan Baek, Chia-Hao Lee, Yeonjoon Jung, Seong Chul Hong, Gillian Nolan, [Kenji Watanabe](https://orcid.org/0000-0003-3701-8119), [Takashi Taniguchi](https://orcid.org/0000-0002-1467-3105), Gwan-Hyoung Lee, Pinshane Y. Huang

## Rights

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

## Other metadata

[Atom-by-atom imaging of moiré transformations in 2D transition metal dichalcogenides](https://mdr.nims.go.jp/datasets/4ba05f93-e32e-4847-8b05-d77906c3a479)

## Fulltext

Atom-by-atom imaging of moiré transformations in 2D transition metal dichalcogenidesZhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e1 of 8C O N D E N S E D  M AT T E R  P H Y S I C SAtom-by-atom imaging of moiré transformations in 2D transition metal dichalcogenidesYichao Zhang1†, Ji-Hwan Baek2†, Chia-Hao Lee1, Yeonjoon Jung2, Seong Chul Hong2,  Gillian Nolan1, Kenji Watanabe3, Takashi Taniguchi4, Gwan-Hyoung Lee2, Pinshane Y. Huang1,5*Understanding the atomic-scale mechanisms that govern the structure of interfaces is critical across materials systems but particularly so for two-dimensional (2D) moiré materials. Here, we image, atom-by-atom, the ther-mally induced structural evolution of twisted bilayer transition metal dichalcogenides using in situ transmission electron microscopy. We observe low-temperature, local conversion of moiré superlattice into nanoscale aligned domains. Unexpectedly, this process occurs by nucleating a new grain within one monolayer, whose crystal orien-tation is templated by the other. The aligned domains grow through collective rotation of moiré supercells and hopping of 5|7 defect pairs at moiré boundaries. This provides mechanistic insight into the atomic-scale interac-tions controlling moiré structures and illustrates the potential to pattern interfacial structure and properties of 2D materials at the nanoscale.INTRODUCTIONA rich array of emergent properties has been found in twisted two-dimensional (2D) moiré systems, including superconductivity (1), strongly correlated states (2–6), moiré excitons (7–10), magnetism (11, 12), and interfacial ferroelectricity (13, 14). The interfacial align-ment, such as the relative rotation, translation, and lattice mismatch between adjacent layers, is a critical knob to modulate the properties of moiré materials. While the average, global alignment of 2D moirés can be tuned through various fabrication and processing techniques (15–18), understanding and controlling the local moiré structures remains a challenge. The structure and properties of 2D multilayers are often highly heterogeneous (19), varying widely both between and within individual samples in response to defects (20, 21), strain (22–24), and pinning sites (25). Moreover, rather than being static, 2D moirés are known to reconfigure readily under external stimuli, such as heating (26–29), bending (30), and shear (31). As a result, understanding the atomic-scale mechanisms that govern the struc-tural evolution of 2D moirés is key to designing the properties of 2D multilayers.Here, we use in situ transmission electron microscopy (TEM) and scanning TEM (STEM) to directly image the thermally driven struc-tural rearrangements of twisted bilayer transition metal dichalcogen-ide (TB-TMD) homo- and heterostructures down to the atomic scale. We prepare single-crystal monolayers of WSe2 and MoSe2 us-ing chemical vapor deposition (CVD) and then fabricate graphene-encapsulated TB-TMDs using the pick-up technique (32, 33). We then transfer the samples onto microelectromechanical system (MEMS)–based E-chips (Fig.  1A) for in  situ heating experiments (33). The MEMS-based heating chips heat and cool rapidly (up to 1000°C/ms) (34), making it possible to access and image metastable states. The graphene encapsulation produces a solid-state reaction chamber that mitigates electron-beam damage and prevents sample sublimation at high temperatures, preserving the atomic-resolution capability of STEM (35). This approach enables acquisitions of mul-tiple (S)TEM images to comprehensively survey the structural trans-formation during heating.RESULTSFigure 1 (A and B) shows our experimental design. We characterize the samples using optical microscopy and STEM imaging before heat-ing to verify the sample quality and stacking angle (fig. S1). We study a total of four samples, two WSe2/WSe2 homo-bilayers and two MoSe2/WSe2 hetero-bilayers with twist angles between 3.3° and 45.0° [see (33) and fig. S2]. Here, 0° and 60° are defined as parallel (AB′ or 3R) and antiparallel (AA′ or 2H) stacking (36), respectively. We apply a sequence of short, 0.5-s heat pulses from 100° to 1000°C to the TB-TMDs (see heating profile in Fig.  1B) and acquire (S)TEM images between pulses to reduce thermal drift.Nucleation of nanoscale aligned domainsFigure 1B shows the interlayer alignment of a MoSe2/WSe2 sample before and after heating. The sample consists of two bilayer regions, one with a 7.6° interlayer twist angle (Fig. 1C) and the other with an 18.2° twist as measured from selected-area electron diffraction (see fig. S3 for TEM images and diffraction patterns). These two different twist regions occur because of a partial tear in the MoSe2 layer during transfer. After the sequence of heat pulses, we observe the formation of nanoscale aligned areas via atomic-resolution annular dark-field (ADF)–STEM (Fig.  1D and fig.  S2, C and D) and dark-field TEM (DFTEM) (fig. S3). Notably, only 0.04 μm2 of the 16.76 μm2 WSe2 lat-tice becomes aligned after the heating (fig. S3), indicating local recon-struction instead of the global, whole-flake rotations that have been previously observed in other 2D systems (3, 26–28). The fractional conversion is likely due to the short overall heating time (<10 s in to-tal, see Fig. 1B). Longer thermal annealing at high temperature (e.g., 800°C for 1 hour) does lead to full conversion of twisted to aligned structure (37). We do not observe any tears between the moiré do-mains and the newly formed aligned domains. Instead, the aligned 1Department of Materials Science and Engineering, University of Illinois Urbana-Champaign, Urbana, IL 61801, USA. 2Department of Materials Science and Engi-neering, Seoul National University, Seoul 08826, Korea. 3Research Center for Electronic and Optical Materials, National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan. 4Research Center for Materials Nanoarchitectonics, National Insti-tute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan. 5Materials Re-search Laboratory, University of Illinois Urbana-Champaign, Urbana, IL 61801, USA.*Corresponding author. Email: pyhuang@​illinois.​edu†These authors contributed equally to this work.Copyright © 2024 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 NonCommercial License 4.0 (CC BY-NC). Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024mailto:pyhuang@​illinois.​eduhttp://crossmark.crossref.org/dialog/?doi=10.1126%2Fsciadv.adk1874&domain=pdf&date_stamp=2024-03-27Zhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e2 of 8domains are produced through the nucleation of a new crystal grain within one layer that is rotationally aligned to the other. For example, in the 7.6° twisted region, the WSe2 produces a new grain aligned with the MoSe2 layer, while the reverse is true in the 18.2° twisted region [see (33) and fig. S4].After the sequence of heat pulses, the bilayer structures are highly heterogeneous. Figure 1E shows an atomic-resolution ADF-STEM image of a WSe2/WSe2 homo-bilayer post-heating, with a global twist angle of 3.3°. In this ~300 nm2 field of view, a range of inhomogeneous stacking arrangements is visible, including twist moirés (3.3°), an AB′-stacked region (0°), and stacking faults. Un-expectedly, we also observe an ~5-nm AA′ stacked (60°) area (yel-low triangle in Fig. 1D), although the initial twist angle is closer to 0° rather than 60°.By comparing the structural transformation of all four samples studied, we can extract general trends about the formation of aligned regions [see (33) and fig. S5 for statistics of aligned regions in all sam-ples]. In each sample, several distinct regions, from several nanome-ters to micrometers across, convert into aligned bilayers after heating [see (33) and figs. S3 and S5]. We observe a larger aligned domain in the 18.2°-twisted region than the 7.6°. Higher twist angles correspond to larger interfacial stacking energy and correspondingly a larger en-ergy difference between twisted and aligned states (38). Therefore, we would expect larger aligned domains to form in the bilayer with large twist angle, although other effects such as pinning, strain, and sample edges may also play substantial roles. Typically, the final stacking ar-rangement of the aligned structure depends on the initial twist angle such that a minimum angular change occurs (i.e., <30° converts to 0° Fig. 1. In situ heating of twisted bilayer WSe2/WSe2 and MoSe2/WSe2 under pulsed heating. (A) Schematic of graphene encapsulated twisted bilayer transition metal dichalcogenides (TMDs) on a MEMS-based heating chip. (B) Heating profile and crystallographic alignment of MoSe2/WSe2 sample before and after heating. The sample un-dergoes a series of 0.5-s heat pulses between 100° and 1000°C. The cartoons and pie charts depict the specimen geometry and relative crystal orientation, before (left) and after (right) heating. Before heating, the sample consists of two bilayer domains with 7.6° and 18.2° relative twist; after heating, a portion of the sample converts to aligned bilayer. DFTEM images associated with the cartoons are provided in the Supplementary Materials. (C and D) Fourier-filtered ADF-STEM images of 7.6°-twisted MoSe2/WSe2 hetero-bilayer before and after in situ heating. (E) Post-heating Fourier-filtered ADF-STEM image of an antiparallel (AA′) domain in 3.3°-twisted WSe2/WSe2 homo-bilayer. Dashed yellow lines outline position of the AA′ domain. Additional stacking orders including moiré, AB′, and shear-misaligned are labeled. Scale bars, 5 nm [(C) to (E)].Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024Zhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e3 of 8and >30° converts to 60°; see fig. S5). Although in our experiments most of the twist angle changes result in aligned domains, we also oc-casionally observe conversion from one misaligned twist to another; for example, a 56.4°-twisted domain formed next to a 45°-twist do-main (fig. S6).In situ growth of aligned domainsTo visualize the nucleation and growth of the aligned domains, we analyze DFTEM images on a MoSe2/WSe2 twisted bilayer sample. Both layers are initially single crystals, as verified by electron diffrac-tion (Fig. 2A). After each heat pulse, a pair of complimentary DFTEM images is acquired using the Bragg spots of each monolayer (see movies S1 and S2). Formation of an aligned domain appears as a bright region in one of the DFTEM images (Fig.  2, B to E, and movie S1) along with a corresponding dark region at the same posi-tion in the complimentary DFTEM image (movie S2). An aligned bilayer domain is visible in the DFTEM images starting at 200°C (see movie S1). Notably, this temperature is in the 90° to 350°C Fig. 2. In situ growth and boundary structure of aligned domain in a MoSe2/WSe2 hetero-bilayer. (A) Selected-area electron diffraction pattern of 7.6°-twisted MoSe2-WSe2 hetero-bilayer. Colored circles mark the positions of objective aperture used to generate DFTEM images of MoSe2 (red) in (B) to (E) and WSe2 (yellow) layers in the Supplementary Materials. Scale bar, 5 nm−1. Inset cartoon illustrates geometry of the hetero-bilayer. (B to E) DFTEM image series under heating, obtained with objective aperture position 2. Bright region represents the formation of a new grain in the WSe2 layer, aligned with the MoSe2 lattice. Images are smoothed using Gaussian blurring with a radius of 2 pixels. (F) Line traces of aligned domain boundaries measured from DFTEM images. The color indicates the temperature of the heat pulse after which the image is acquired. The brown arrow illustrates the growth direction of the aligned domain. Scale bars, 100 nm (B to F). (G) Fourier-filtered ADF-STEM image of interface between the aligned domain and the moiré lattice. A 7.6° rotational grain boundary is located in the WSe2 layer. Scale bar, 2 nm. (H and I) Structural cartoons of MoSe2 and WSe2 layers of the region shown in (G). W atoms in three 5|7 defect pairs are labeled in blue. (J) ADF-STEM image overlaid with grain boundary structure.Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024Zhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e4 of 8range commonly used to anneal 2D heterostructures (1, 39), indicat-ing that nanoscale aligned regions may be introduced during typical 2D device processing.Figure 2F shows how the positions of the domain boundaries of the aligned region evolve during heating, extracted via edge detection on DFTEM images (33). After the full heating sequence, the aligned domain is ~0.03 μm2, with an irregular and anisotropic shape. The aligned domain growth approximately follows an Arrhenius relation with respect to temperature from 200° to 1000°C [see (33) and fig. S7], consistent with the temperature dependence of grain boundary mo-bility in bicrystals (40). From the Arrhenius plot, we extract an activa-tion energy of this reconstruction to be 3.54 ± 0.03 eV/nm2. This is equivalent to four bonds broken per moiré supercell, considering the 4.3-eV bond dissociation energy in WSe2 (41).Figure  2G shows a post-heating, atomic-resolution ADF-STEM image of the domain boundary in a MoSe2/WSe2 heterostructure. We solve the structure by separating the two monolayers in the image using Fourier filtering (Fig. 2, H to J; see Materials and Methods). As illustrated in Fig. 2 (H to J), the domain boundary is formed by a grain boundary in the WSe2 layer consisting of 5|7 defect pairs. We ob-serve a similar domain boundary structure in the WSe2/WSe2 bilayer (fig. S9). The grain boundary structures are consistent with previous observations of the atomic structure of grain boundaries in 2D TMDs (42, 43).The formation and growth of an aligned domain is a competition between the interfacial stacking energy (which decreases with the area of the aligned domain) and the grain boundary energy (which increases with the length and twist angle of the grain boundary). Ac-cording to first-principles studies (38), WSe2 twisted bilayers are 0.2 to 0.3 eV/nm2 higher energy than aligned bilayers, with similar values expected for other TMD bilayers (44). Meanwhile, the energetic cost associated with a grain boundary is 1.4 to 3 eV/nm (45), for the grain boundary angles in our study. By balancing these two energetic terms, we can estimate a critical size above which growth of the aligned do-mains is favored. For twisted bilayer WSe2, we find that the critical domain radius, which is a function of the twist angle, ranges between 1 and 20 nm [see (33) and fig. S8].Next, we study the atomic-scale growth mechanism for the aligned domains using in situ aberration-corrected STEM. Figure 3 (A to D) shows an atomic-resolution ADF-STEM image series detailing the grain boundary migration (see movie S3 for complete image series) in a 45.0°-twisted WSe2/WSe2 homo-bilayer. After six 0.5-s heat pulses between 400° and 900°C, the grain boundary, indicated by colored lines in Fig. 3, displaces up to 5 nm. We also observe the formation of vacancies, multi-atom voids (Fig.  3A), and small nearly aligned (57.2°-twisted) domains (fig. S10) during heating; such defects as well as strain may alter the activation energies and equilibrium struc-tures formed.As shown in Fig. 3 (B to D), the 5|7 dislocation cores are preferen-tially located at the edges of the moiré supercells in each frame; these positions are represented as the unshaded regions, where the inter-layer lattice misalignment is greatest. For example, Fig. 3 (B and C) shows that, between the 500° and 600°C frames, two 5|7 pairs hop to approximately equivalent sites at the adjacent moiré supercell, thus converting an area of the size of one moiré unit cell into aligned bi-layer (gray diamond in Fig. 3C). Correspondingly, the displacement of the domain boundary progresses in discrete units, approximately equal to the radius of a moiré supercell (1.1 nm at θ  =  45.0°; see fig. S11 for additional regions). This mechanism is notably different from grain boundary migration in monolayer TMDs, which proceeds by continuous gliding of defect pairs (32, 46, 47).Figure  4 analyzes the atom-by-atom transformation during the aligned domain growth. Focusing on a small subregion from Fig. 3A, we obtain atomic coordinates from each STEM image in the tempera-ture series and track the atoms between consecutive image frames us-ing a global best-fit algorithm (33, 48). Figure 4A shows the atomic displacements of the reconstructing WSe2 monolayer after the 700°C heat pulse, overlaid on a STEM image acquired just before the heat pulse. We observe three main types of atomic motions, depicted in Fig. 4 (B to D), which correspond to the three main regions in the sample. The most notable rearrangements occur in region 1 (Fig. 4A), which converts from twisted to aligned bilayer after the pulse. Here, we observe that the atoms in one WSe2 monolayer rearrange, produc-ing rotational displacements with the periodicity of the moiré lattice (Fig. 4B). In contrast, in region 2 (Fig. 4A), which contains aligned bilayer both before and after the heat pulse, the atoms remain station-ary (Fig. 4C) to within our measurement error. In region 3, which represents twisted bilayer in both time steps, we observe large, 80- to 154-pm translations (Fig. 4D), corresponding to sliding of the WSe2 layer. The magnitude of such translation is largest (~150 pm) at the 700°C grain boundary and decreases linearly away from the grain boundary (Fig. 4E), indicating that it accommodates changes in stress resulting from the dislocation motion.We calculate the local change in rotation and shear between frames using the spatial derivatives of the displacement fields (Fig. 4, E and F) (33, 48). The rotation maps in Fig. 4F show a uniform counterclock-wise rotation of +15° within each moiré supercell, consistent with the angle difference between the initial 45° twist angle and final AA′ aligned stacking (60°). Figure 4G shows that shear concentrates at the supercell edges, where displacement vectors across the boundary point in opposite directions. Note that values in Fig. 4 (F and G) are defined as changes in strain between frames, not an absolute strain relative to the equilibrium atomic structure.Together, the data in Figs. 3 and 4 indicate that, to allow each moi-ré supercell to rotate while maintaining a continuous 2D layer, bonds break and reform at the edges of the moiré supercells, while atoms stay bonded within the supercells. This mechanism minimizes the number of covalent bonds broken per unit area during the growth of the aligned bilayer and is consistent with the small activation energy we measured using DFTEM image series. Because the interlayer alignment progresses via local rotation of nanometer-scale regions, this mechanism should be favored over whole-flake rotation when the activation energy for the creation and migration of dislocations is smaller than the sliding energy barrier integrated over the system size. We expect this mechanism to be applicable to twist angles within the rigid moiré and part of the small-angle reconstruction regimes (in-cluding 7.6°, rigid moiré and the 3°, partially reconstructed samples measured here). In the case of extremely small twist angles (0° to 2°), a large fraction of the moiré superlattice is already aligned. In the small twist angle limit, where the system can be described as an aligned lattice containing a dilute soliton network, thermally activated soliton motion has previously been shown to dominate lattice recon-structions (20).DISCUSSIONIn summary, using in situ (S)TEM in a solid-state graphene reac-tion chamber, we report a mechanism for the nucleation and Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024Zhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e5 of 8growth of nanometer-scale aligned domains in TB-TMDs. Here, bilayers align through the formation and migration of grain boundary dislocations in one of the 2D layers, combined with the collective motion of atoms within a moiré supercell. This process enables the nucleation of aligned domains in both homo- and hetero-bilayers with a wide range of twist angles at temperatures as low as 200°C, within the range of typical processing tempera-tures for 2D heterostructures. Our results can explain some of the inconsistent and heterogeneous properties measured in moiré-based devices and may lead to new, bottom-up routes to fabricate and design structures to control and pattern the twist angle on the nanoscale.MATERIALS AND METHODSCVD growth of MoSe2 and WSe2 monolayersTo synthesize monolayers of WSe2 and MoSe2 by CVD, we used SiO2 (285 nm)/Si substrates in 2-inch quartz tube under ambient pressure. For monolayer WSe2, we used a quartz boat with WO2.9 powder (99.99%, Alfa Aesar) and KI powder (99.99%, Alfa Aesar) placed at the center of the furnace and the SiO2/Si substrate suspended face-down on top of the boat. Another quartz boat with selenium powder (99.99%, Sigma-Aldrich) was located upstream of the quartz tube at a distance of 19.5 cm from the center of the furnace. The furnace tem-perature was raised to 800°C at a rate of 50°C/min, maintained for 15 min, and cooled to room temperature. We flew 300 standard cubic centimeters per minute (sccm) of Ar during the whole growth process and 10 sccm of H2 during the cooling process. For monolayer MoSe2, we used MoO3 powder (99.97%, Sigma-Aldrich). The furnace tem-perature was raised to 750°C at a rate of 50°C/min, maintained for 20 min, and cooled to room temperature. All other parameters are the same as used for WSe2.Fabrication of graphene encapsulated twisted bilayer TMD homo/heterostructuresWe exfoliated flakes of hexagonal boron nitride (hBN) and graphene from bulk crystals onto SiO2 (285 nm)/Si substrates. Monolayer WSe2 and MoSe2 were either mechanically exfoliated from bulk crystals Fig. 3. Atomic-resolution imaging of in situ growth of aligned domain in a 45°-twisted WSe2/WSe2 homo-bilayer. (A) ADF-STEM image of 45.0° moiré-AA′ interface at 900°C. Colored lines mark the positions of the interface from 400° to 900°C. Nanometer-scale voids and trilayer patches are also visible as darker and brighter regions. Scale bar, 5 nm. (B to D) Magnified ADF-STEM images of the interface from the region marked with the white box in (A) after heat pulses at 500°, 600°, and 700°C. Dashed colored lines mark the positions of the interface. Colored pentagons and heptagons mark the positions of the 5|7 defect pairs at the interface. Colored blocks mark moiré supercells. Gray pentagons and heptagons in (C) mark positions of 5|7 defect pairs after 500°C heat pulse. Gray diamond in (C) marks a single moiré cell transformed to aligned region after 600°C heat pulse. Scale bar, 2 nm. All images are aligned by cross-correlating the ADF-STEM images in the large (~100 nm) field of view.Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024Zhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e6 of 8(from HQ Graphene) or synthesized from CVD. We identified the thickness of monolayers with optical contrast, photoluminescence measurements, and atomic force microscopy. Next, we assembled the hBN/Gr/WSe2/MoSe2/Gr heterostructures with polycarbonate (PC)/polydimethylsiloxane (PDMS) stamp using dry-transfer method (49). The pick-up temperature for hBN, TMDs, and graphene was around 115°, 75°, and 75°C. After stacking, we transferred the heterostruc-tures onto a clean SiO2 (285 nm)/Si substrate by releasing the PC film from the PDMS lens at 180°C. Transferred samples were placed in chloroform for 3 hours to remove PC. We used few-layer graphene to prevent degradation of heterostructures during annealing and hBN as top layer to detach graphene from the substrate (50). We then used XeF2 gas to selectively etch hBN (51) to obtain graphene encapsulated TB-TMD samples.Sample preparation for in situ (S)TEM experimentsWe prepared TEM samples using a poly(methyl methacrylate) (PMMA)–based wet transfer method. We transferred graphene en-capsulated TB-TMD/PMMA on MEMS-based heating chips (Pro-tochips, E-FHDN-VO-10) and removed PMMA by acetone bath for 24 hours.Electron microscopy imagingThe samples were imaged in a Thermo Fisher Scientific Themis Z aberration–corrected (S)TEM. The microscope was operated at 80 kV to minimize knock-on damage of the 2D materials. The aberration-corrected ADF-STEM images and DFTEM images were obtained us-ing the STEM and TEM capabilities of the instrument, respectively. We used a convergence angle of 25.2 mrad.In situ heatingWe used a heating holder (Fusion Select, Protochips) to apply heat pulses of 0.5 s from 100° to 1000°C. DFTEM and/or ADF-STEM im-ages were acquired after each 0.5-s heat pulse was applied, and the samples were quenched to room temperature. The temperatures for DFTEM image series are 100° to 400°C at 100°C temperature steps and 400° to 1000°C at 50°C temperature steps. The temperatures for ADF-STEM imaging are 400° to 900°C at 100°C temperature steps.Supplementary MaterialsThis PDF file includes:Materials and MethodsFig. 4. Atomic displacements in a moiré-to-antiparallel reconstruction. (A) Atomic displacement field showing the lattice transformation before and after 700°C pulse. A Fourier-filtered ADF-STEM image, acquired before the 700°C pulse, is overlaid with the vectors (green arrows) showing atomic displacements after the 700°C pulse. Vectors are 1:1 scale. Colored dashed lines label positions of the grain boundary at 600° and 700°C. Pentagons and heptagons mark 5|7 defect pairs. Scale bar, 1 nm. (B to D) Magnified regions of interest (ROIs) displaying rotational, close to stationary, and translational atomic displacements, respec-tively. (E) Line profile of average magnitude of the displacement vectors in ROI 3 as a function of position. Zero position is defined at the 700°C grain boundary. Inset is cropped from ROI 3 in (A). Yellow arrow marks the position and direction of line profile. (F and G) Strain tensor components derived from the displace-ment field in (A) showing rotation (ω) and shear [εxy + 12(εxx − εyy)]. In (F), positive rotations are defined as counterclockwise. Color wheel in the phase map in (G) represents the shear direction (color) and magnitude (intensity). The strain and rotation fields are relative to the previous frame rather than equilibrium atomic structures.Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024Zhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e7 of 8Figs. S1 to S12Tables S1 and S2Legends for movies S1 to S3ReferencesOther Supplementary Material for this manuscript includes the following:Movies S1 to S3REFERENCES AND NOTES  1.  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).  2.  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).  3.  F. Wu, T. Lovorn, E. Tutuc, I. Martin, A. H. Macdonald, Topological insulators in twisted transition metal dichalcogenide homobilayers. Phys. Rev. Lett. 122, 086402 (2019).  4. L . Wang, E. M. Shih, A. Ghiotto, L. Xian, D. A. Rhodes, C. Tan, M. Claassen, D. M. Kennes,  Y. Bai, B. Kim, K. Watanabe, T. Taniguchi, X. Zhu, J. Hone, A. Rubio, A. N. Pasupathy,  C. R. Dean, Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 19, 861–866 (2020).  5. E . C. Regan, D. Wang, C. Jin, M. I. Bakti Utama, B. Gao, X. Wei, S. Zhao, W. Zhao, Z. Zhang,  K. Yumigeta, M. Blei, J. D. Carlström, K. Watanabe, T. Taniguchi, S. Tongay, M. Crommie,  A. Zettl, F. Wang, Mott and generalized Wigner crystal states in WSe2/WS2 moiré superlattices. Nature 579, 359–363 (2020).  6.  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).  7.  K. L. Seyler, P. Rivera, H. Yu, N. P. Wilson, E. L. Ray, D. G. Mandrus, J. Yan, W. Yao, X. Xu, Signatures of moiré-trapped valley excitons in MoSe2/WSe2 heterobilayers. Nature 567, 66–70 (2019).  8.  Z. Wang, D. A. Rhodes, K. Watanabe, T. Taniguchi, J. C. Hone, J. Shan, K. F. Mak, Evidence of high-temperature exciton condensation in two-dimensional atomic double layers. Nature 574, 76–80 (2019).  9.  K. Tran, G. Moody, F. Wu, X. Lu, J. Choi, K. Kim, A. Rai, D. A. Sanchez, J. Quan, A. Singh,  J. Embley, A. Zepeda, M. Campbell, T. Autry, T. Taniguchi, K. Watanabe, N. Lu,  S. K. Banerjee, K. L. Silverman, S. Kim, E. Tutuc, L. Yang, A. H. MacDonald, X. Li, Evidence for moiré excitons in van der Waals heterostructures. Nature 567, 71–75 (2019).  10.  X. Zhao, J. Qiao, X. Zhou, H. Chen, J. Y. Tan, H. Yu, S. M. Chan, J. Li, H. Zhang, J. Zhou, J. Dan, Z. Liu, W. Zhou, Z. Liu, B. Peng, L. Deng, S. J. Pennycook, S. Y. Quek, K. P. Loh, Strong moiré excitons in high-angle twisted transition metal dichalcogenide homobilayers with robust commensuration. Nano Lett. 22, 203–210 (2022).  11.  A. L. Sharpe, E. J. Fox, A. W. Barnard, J. Finney, K. Watanabe, T. Taniguchi, M. A. Kastner,  D. Goldhaber-Gordon, Emergent ferromagnetism near three-quarters filling in twisted bilayer graphene. Science 365, 605–608 (2019).  12.  M. Serlin, C. L. Tschirhart, H. Polshyn, Y. Zhang, J. Zhu, K. Watanabe, T. Taniguchi, L. Balents, A. F. Young, Intrinsic quantized anomalous Hall effect in a moiré heterostructure. Science 367, 900–903 (2020).  13.  A. Weston, E. G. Castanon, V. Enaldiev, F. Ferreira, S. Bhattacharjee, S. Xu, H. Corte-León,  Z. Wu, N. Clark, A. Summerfield, T. Hashimoto, Y. Gao, W. Wang, M. Hamer, H. Read,  L. Fumagalli, A. V. Kretinin, S. J. Haigh, O. Kazakova, A. K. Geim, V. I. Fal’ko, R. Gorbachev, Interfacial ferroelectricity in marginally twisted 2D semiconductors. Nat. Nanotechnol. 17, 390–395 (2022).  14.  M. Lv, X. Sun, Y. Chen, T. Taniguchi, K. Watanabe, M. Wu, J. Wang, J. Xue, Spatially resolved polarization manipulation of ferroelectricity in twisted hBN. Adv. Mater. 34, 2203990 (2022).  15.  Y. Cao, J. Y. Luo, V. Fatemi, S. Fang, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi,  E. Kaxiras, P. Jarillo-Herrero, Superlattice-induced insulating states and valley-protected orbits in twisted bilayer graphene. Phys. Rev. Lett. 117, 116804 (2016).  16.  X.-D. Chen, W. Xin, W.-S. Jiang, Z.-B. Liu, Y. S. Chen, J.-G. Tian, High-precision twist-controlled bilayer and trilayer graphene. Adv. Mater. 28, 2563–2570 (2016).  17.  A. J. Mannix, A. Ye, S. H. Sung, A. Ray, F. Mujid, C. Park, M. Lee, J. H. Kang, R. Shreiner,  A. A. High, D. A. Muller, R. Hovden, J. Park, Robotic four-dimensional pixel assembly of van der Waals solids. Nat. Nanotechnol. 17, 361–366 (2022).  18. H . Zheng, H. Guo, S. Chen, B. Wu, S. Li, J. He, Z. Liu, G. Lu, X. Duan, A. Pan, Y. Liu, Strong interlayer coupling in twisted transition metal dichalcogenide moiré superlattices. Adv. Mater. 35, 2210909 (2023).  19. T . Benschop, T. A. de Jong, P. Stepanov, X. Lu, V. Stalman, S. J. van der Molen, D. K. Efetov, M. P. Allan, Measuring local moiré lattice heterogeneity of twisted bilayer graphene. Phys. Rev. Res. 3, 013153 (2021).  20.  J. S. Alden, A. W. Tsen, P. Y. Huang, R. Hovden, L. Brown, J. Park, D. A. Muller, P. L. McEuen, Strain solitons and topological defects in bilayer graphene. Proc. Natl. Acad. Sci. U.S.A. 110, 11256–11260 (2013).  21. D . Edelberg, H. Kumar, V. Shenoy, H. Ochoa, A. N. Pasupathy, Tunable strain soliton networks confine electrons in van der Waals materials. Nat. Phys. 16, 1097–1102 (2020).  22. L . Huder, A. Artaud, T. Le Quang, G. T. de Laissardière, A. G. M. Jansen, G. Lapertot,  C. Chapelier, V. T. Renard, Electronic spectrum of twisted graphene layers under heterostrain. Phys. Rev. Lett. 120, 156405 (2018).  23.  Z. Bi, N. F. Q. Yuan, L. Fu, Designing flat bands by strain. Phys. Rev. B. 100, 035448 (2019).  24.  S. Zhu, P. Pochet, H. T. Johnson, Controlling rotation of two-dimensional material flakes. ACS Nano 13, 6925–6931 (2019).  25. L . Brown, R. Hovden, P. Huang, M. Wojcik, D. A. Muller, J. Park, Twinning and twisting of tri- and bilayer graphene. Nano Lett. 12, 1609–1615 (2012).  26. 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).  27. D . Wang, G. Chen, C. Li, M. Cheng, W. Yang, S. Wu, G. Xie, J. Zhang, J. Zhao, X. Lu, P. Chen,  G. Wang, J. Meng, J. Tang, R. Yang, C. He, D. Liu, D. Shi, K. Watanabe, T. Taniguchi, J. Feng,  Y. Zhang, G. Zhang, Thermally induced graphene rotation on hexagonal boron nitride. Phys. Rev. Lett. 116, 126101 (2016).  28.  F. Peymanirad, S. Kumar Singh, H. Ghorbanfekr-Kalashami, K. S. Novoselov, F. M. Peeters, M. Neek-Amal, Thermal activated rotation of graphene flake on graphene. 2D Mater. 4, 025015 (2017).  29.  X. Wang, Y. Zhao, X. Kong, Q. Zhang, H. K. Ng, S. X. Lim, Y. Zheng, X. Wu, K. Watanabe,  Q.-H. Xu, T. Taniguchi, G. Eda, K. E. J. Goh, S. Jin, K. P. Loh, F. Ding, W. Sun, C. H. Sow, Dynamic tuning of moiré superlattice morphology by laser modification. ACS Nano 16, 8172–8180 (2022).  30.  J. Yu, E. Han, M. A. Hossain, K. Watanabe, T. Taniguchi, E. Ertekin, A. M. van der Zande,  P. Y. Huang, Designing the bending stiffness of 2D material heterostructures. Adv. Mater. 33, 2007269 (2021).  31.  Y. Liu, J. N. B. Rodrigues, Y. Z. Luo, L. Li, A. Carvalho, M. Yang, E. Laksono, J. Lu, Y. Bao, H. Xu, S. J. R. Tan, Z. Qiu, C. H. Sow, Y. P. Feng, A. H. C. Neto, S. Adam, J. Lu, K. P. Loh, Tailoring sample-wide pseudo-magnetic fields on a graphene–Black phosphorus heterostructure. Nat. Nanotechnol. 13, 828–834 (2018).  32.  X. Zou, M. Liu, Z. Shi, B. I. Yakobson, Environment-controlled dislocation migration and superplasticity in monolayer MoS2. Nano Lett. 15, 3495–3500 (2015).  33.  Materials and methods are available as supplementary materials.  34.  M. Briceno, K. Hattar, J. Damiano, D. Nackashi, I. M. Robertson, In-situ TEM observations on the sintering process of colloidal gold using an ultra-fast heating stage. Microsc. Microanal. 14, 1336–1337 (2008).  35. C .-H. Lee, H. Ryu, G. Nolan, Y. Zhang, Y. Lee, S. Oh, H. Cheong, K. Watanabe, T. Taniguchi,  K. Kim, G.-H. Lee, P. Y. Huang, In situ imaging of an anisotropic layer-by-layer phase transition in few-layer MoTe2. Nano Lett. 23, 677–684 (2023).  36.  J. He, K. Hummer, C. Franchini, Stacking effects on the electronic and optical properties of bilayer transition metal dichalcogenides MoS2, MoSe2, WS2, and WSe2. Phys. Rev. B 89, 075409 (2014).  37.  J.-H. Baek, H. G. Kim, S. Y. Lim, S. C. Hong, Y. Chang, H. Ryu, Y. Jung, H. Jang, J. Kim,  Y. Zhang, K. Watanabe, T. Taniguchi, P. Y. Huang, H. Cheong, M. Kim, G.-H. Lee, Thermally induced atomic reconstruction into fully commensurate structures of transition metal dichalcogenide layers. Nat. Mater. 22, 1463–1469 (2023).  38.  Z. Zhang, J. Li, G. Yang, G. Ouyang, Interface engineering of band evolution and transport properties of bilayer WSe2 different electric fields. J. Phys. Chem. C 123, 19812–19819 (2019).  39. L . Wang, I. Meric, P. Y. Huang, Q. Gao, Y. Gao, H. Tran, T. Taniguchi, K. Watanabe,  L. M. Campos, D. A. Muller, J. Guo, P. Kim, J. Hone, K. L. Shepard, C. R. Dean, One-dimensional electrical contact to a two-dimensional material. Science 342, 614–617 (2013).  40.  K. Chen, J. Han, D. J. Srolovitz, On the temperature dependence of grain boundary mobility. Acta Mater. 194, 412–421 (2020).  41.  A. Sevy, R. F. Huffaker, M. D. Morse, Bond dissociation energies of tungsten molecules: WC, WSi, WS, WSe, and WCl. J. Phys. Chem. A 121, 9446–9457 (2017).  42.  S. Najmaei, Z. Liu, W. Zhou, X. Zou, G. Shi, S. Lei, B. I. Yakobson, J.-C. Idrobo, P. M. Ajayan,  J. Lou, Vapour phase growth and grain boundary structure of molybdenum disulphide atomic layers. Nat. Mater. 12, 754–759 (2013).  43.  W. Zhou, X. Zou, S. Najmaei, Z. Liu, Y. Shi, J. Kong, J. Lou, P. M. Ajayan, B. I. Yakobson,  J.-C. Idrobo, Intrinsic structural defects in monolayer molybdenum disulfide. Nano Lett. 13, 2615–2622 (2013).  44.  A. A. Puretzky, L. Liang, X. Li, K. Xiao, B. G. Sumpter, V. Meunier, D. B. Geohegan, Twisted MoSe2 bilayers with variable local stacking and interlayer coupling revealed by low-frequency Raman spectroscopy. ACS Nano 10, 2736–2744 (2016).  45. L . Yang, J. Liu, Y. Lin, K. Xu, X. Cao, Z. Zhang, J. Wu, Strengthening and weakening by dislocations in monolayer MoS2. Chem. Mater. 33, 8758–8767 (2021).  46.  A. Azizi, X. Zou, P. Ercius, Z. Zhang, A. L. Elías, N. Perea-López, G. Stone, M. Terrones,  B. I. Yakobson, N. Alem, Dislocation motion and grain boundary migration in two-dimensional tungsten disulphide. Nat. Commun. 5, 4867 (2014).Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024Zhang et al., Sci. Adv. 10, eadk1874 (2024)     27 March 2024S c i e n c e  A d v a n c e s  |  R e s e ar  c h  A r t i c l e8 of 8  47.  X. Zhao, Y. Ji, J. Chen, W. Fu, J. Dan, Y. Liu, S. J. Pennycook, W. Zhou, K. P. Loh, Healing of planar defects in 2D materials via grain boundary sliding. Adv. Mater. 31, 1900237 (2019).  48.  P. Y. Huang, S. Kurasch, J. S. Alden, A. Shekhawat, A. A. Alemi, P. L. McEuen, J. P. Sethna,  U. Kaiser, D. A. Muller, Imaging atomic rearrangements in two-dimensional silica glass: Watching silica’s dance. Science 342, 224–227 (2013).  49.  F. Pizzocchero, L. Gammelgaard, B. S. Jessen, J. M. Caridad, L. Wang, J. Hone, P. Bøggild,  T. J. Booth, The hot pick-up technique for batch assembly of van der Waals heterostructures. Nat. Commun. 7, 11894 (2016).  50. H . Rokni, W. Lu, Direct measurements of interfacial adhesion in 2D materials and van der Waals heterostructures in ambient air. Nat. Commun. 11, 5607 (2020).  51.  J. Son, J. Kwon, S. P. Kim, Y. Lv, J. Yu, J. Y. Lee, H. Ryu, K. Watanabe, T. Taniguchi,  R. Garrido-Menacho, N. Mason, E. Ertekin, P. Y. Huang, G.-H. Lee, A. M. van der Zande, Atomically precise graphene etch stops for three dimensional integrated systems from two dimensional material heterostructures. Nat. Commun. 9, 3988 (2018).  52.  J. Schindelin, I. Arganda-Carreras, E. Frise, V. Kaynig, M. Longair, T. Pietzsch, S. Preibisch,  C. Rueden, S. Saalfeld, B. Schmid, J.-Y. Tinevez, D. J. White, V. Hartenstein, K. Eliceiri,  P. Tomancak, A. Cardona, Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).  53.  A. Nalin Mehta, J. Mo, G. Pourtois, A. Dabral, B. Groven, H. Bender, P. Favia, M. Caymax,  W. Vandervorst, Grain-boundary-induced strain and distortion in epitaxial milayer MoS2 lattice. J. Phys. Chem. C 124, 6472–6478 (2020).  54.  A. MacEri, Theory of Elasticity (Springer, 2010).  55.  A. De Backer, K. H. W. van den Bos, W. Van den Broek, J. Sijbers, S. Van Aert, StatSTEM: An efficient approach for accurate and precise model-based quantification of atomic resolution electron microscopy images. Ultramicroscopy 171, 104–116 (2016).  56. D . Blair, E. Dufresne, The MATLAB Particle Tracking Code Repository (2008), p. 2019.  57.  J. Crocker, D. Grier, Methods of digital video microscopy for colloidal studies. J. Colloid Interface Sci. 179, 298–310 (1996).Acknowledgments: Part of this work was carried out in the Materials Research Laboratory at the University of Illinois Urbana-Champaign. K.W. and T.T. acknowledge support from World Premier International Research Center Initiative (WPI), MEXT, Japan. G.-H.L. acknowledges support from the Creative-Pioneering Researchers Program, the Research Institute of Advanced Materials (RIAM), Institute of Engineering Research (IER), Institute of Applied Physics (IAP), and Inter-University Semiconductor Research Center (ISRC) at the Seoul National University. Funding: This work was supported by Department of Energy grant DE-SC0020190 (P.Y.H.); JSPS KAKENHI grants 20H00354, 21H05233, and 23H02052 (K.W. and T.T.); and Korean National Research Foundation grants 2021R1A2C3014316, 2021M3F3A2A01037858, and 2017R1A5A1014862 (G.-H.L.). Author contributions: Conceptualization: P.Y.H. and G.-H.L. Methodology: Y.Z., J.-H.B., C.-H.L., Y.J., S.C.H., G.N., K.W., and T.T. Investigation: Y.Z., C.-H.L., and G.N. Visualization: Y.Z. Funding acquisition: P.Y.H. and G.-H.L. Project administration: P.Y.H. and G.-H.L. Supervision: P.Y.H. and G.-H.L. Writing—original draft: Y.Z. and P.Y.H. Writing—review and editing: Y.Z., J.-H.B., C.-H.L., Y.J., S.C.H., G.N., K.W., T.T., G.-H.L., and P.Y.H. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions of the paper are present in the paper, Supplementary Materials, and Zenodo: https://doi.org/10.5281/zenodo.10058870.Submitted 7 August 2023 Accepted 21 February 2024 Published 27 March 2024 10.1126/sciadv.adk1874Downloaded from https://www.science.org at National Institute for Materials Science on March 29, 2024https://doi.org/10.5281/zenodo.10058870 Atom-by-atom imaging of moiré transformations in 2D transition metal dichalcogenides INTRODUCTION RESULTS Nucleation of nanoscale aligned domains In situ growth of aligned domains DISCUSSION MATERIALS AND METHODS CVD growth of MoSe2 and WSe2 monolayers Fabrication of graphene encapsulated twisted bilayer TMD homo/heterostructures Sample preparation for in situ (S)TEM experiments Electron microscopy imaging In situ heating Supplementary Materials This PDF file includes: Other Supplementary Material for this manuscript includes the following: REFERENCES AND NOTES Acknowledgments