Abstract
Recent demonstrations of moiré magnetism, featuring exotic phases with noncollinear spin order in the twisted van der Waals (vdW) magnet chromium triiodide CrI3, have highlighted the potential of twist engineering of magnetic (vdW) materials. However, the local magnetic interactions, spin dynamics, and magnetic phase transitions within and across individual moiré supercells remain elusive. Taking advantage of a scanning single-spin magnetometry platform, here we report observation of two distinct magnetic phase transitions with separate critical temperatures within a moiré supercell of small-angle twisted double trilayer CrI3. By measuring temperature-dependent spin fluctuations at the coexisting ferromagnetic and antiferromagnetic regions in twisted CrI3, we explicitly show that the Curie temperature of the ferromagnetic state is higher than the Néel temperature of the antiferromagnetic one by ~10âK. Our mean-field calculations attribute such a spatial and thermodynamic phase separation to the stacking order modulated interlayer exchange coupling at the twisted interface of moiré superlattices.
Similar content being viewed by others
Introduction
Recently, twisted two-dimensional (2D) van der Waals (vdW) magnets have emerged as a new member in the suite of moiré quantum materials, where the spin degree of freedom is controlled through the local stacking order of moiré superlattices1,2,3,4,5,6,7,8,9,10,11. As opposed to moiré quantum electronic matter such as twisted graphene12,13,14,15 and transition metal dichalcogenides14,15,16 that feature enhanced electronic interactions from moiré flat bands, moiré magnets rely on twist-engineered magnetic interaction competitions to realize unconventional magnetic orders and spin excitations over moiré wavelengths, for example, in twisted vdW magnet chromium triiodide CrI31,2,3,4,5,7,8,9,11.
Few-layer CrI3 shows stacking-dependent interlayer magnetic exchange coupling that is ferromagnetic (FM) for the rhombohedral stacking and antiferromagnetic (AFM) for the monoclinic stacking geometry8,17,18,19, whereas the intralayer exchange coupling is FM with a strong out-of-plane easy-axis anisotropy8,17,18. In twisted CrI3 with spatially distributed rhombohedral and monoclinic stacking geometries, the competition between uniform FM intralayer exchange coupling and FM-AFM modulated moiré interlayer exchange interaction drives the formation of a range of novel magnetic orders, such as co-existing FM and AFM states within moiré supercells in twisted bilayer and double trilayer CrI31,2, and emergent magnetization and noncollinear spins in twisted double bilayer CrI33,4,5. So far, the ongoing research on twisted CrI3 has mainly focused on investigating spatially modulated magnetic orders of the ground states1,2,3,4,5,11, which represents only a subset of information about moiré magnetism. The local magnetic interactions, spin dynamics, and magnetic phase transitions within and across moiré supercells await exploration and are necessary for developing a comprehensive picture of moiré magnetism.
Here, we report scanning single-spin quantum sensing20,21,22 of both static magnetization and dynamic spin fluctuations of moiré magnetism hosted by twisted double trilayer (tDT) CrI3 across the second-order magnetic phase transition points (Tc). We show that the FM region within individual moiré supercells formed in small-twist-angle tDT CrI3 exhibit a higher Tc up to ~58âK in comparison with that of ~48âK for their AFM counterparts resulting in a nanoscale co-existing paramagnetic(PM)-ferromagnetic (FM) phase in an intermediate temperature regime (48âK \( < \) T \( < \) 58âK), while such a phenomenon is absent in the large-twist-angle regime. Our experimental results are well explained by a proposed mean-field theoretical model of layer-resolved magnetic phases of tDT CrI3 taking account of stacking engineered exchange interactions at the twisted interface. The current work highlights twist engineering as a promising tuning knob to realize local control of magnetic responses at individual stacking sites, which could contribute to a broad range of emerging 2D electronic applications23,24. The new insights on moiré magnetism presented in this study further highlight the potential of quantum metrology tools20 in exploring unconventional spin-related phenomena in correlated magnetic quantum states of matter.
Results
We first briefly review the pertinent material properties of tDT CrI3, which provides a reliably high-quality moiré magnet platform for the current study2. Figure 1a shows a moiré superlattice structure formed by stacking two CrI3 trilayers with a small twist angle. The local atomic registry exhibits a periodic modulation in real space, leading to spatially alternating stacking geometries on a length scale of moiré wavelengths1,2,3. At the monoclinic (ABâ) stacking site, the two CrI3 trilayers are coupled by a positive exchange interaction JM, leading to local AFM order at the twisted interface with a fully compensated net magnetic moment in the ground state2,8. In contrast, FM order with a net magnetic moment is established at the rhombohedral (AB) stacking site driven by a negative interlayer exchange interaction JR2,8. Using first-principles calculations, a theoretical study has reported that the magnitude of JR could be one order of magnitude larger than that of JM due to the orbital-dependent exchange coupling8. Fundamentally, the sign and magnitude of magnetic interactions in twisted CrI3 not only determine its local magnetic phases but also affect their dynamic responses to external perturbations. The former has been predicted and experimentally demonstrated recently2,8, and the latter is the focus of the current work.
We fabricated tDT CrI3 devices by the standard âtear-and-stackâ technique and encapsulated them with hexagonal boron nitride (hBN) nanoflakes1,2,3,12,25. The samples for transmission electron microscopy and quantum sensing measurements were fabricated separately because they required different substrates for individual measurement purposes (See Methods Section for details). Selected area electron diffraction (SAED) patterns (Fig. 1b) of a surveyed sample area of a tDT CrI3 device shows sets of the fifth-order Bragg peaks with sixfold rotation symmetry. The local mean twist angle is measured to be 0.8° \(\pm\) 0.1° by fitting 2D Gaussians to the diffraction peaks, giving a moiré period of ~50ânm. Such an intermediate twist angle ensures a decent, regular moiré lattice structure formed in tDT CrI3 that is visible in transmission electron microscopy measurements. Figure 1c presents a bright-field transmission electron microscopy (BF-TEM) image, which shows the characteristic hexagonal superlattice structures with a periodicity commensurate with the moiré wavelength. The distortions from the expected moiré lattice patterns could be induced by lattice strain, relaxation, and local structural inhomogeneities (see Supplementary Information Note 1 for details)26.
Next, we utilize scanning nitrogen-vacancy (NV) microscopy2,18,27,28,29 to spatially resolve the moiré magnetism hosted by tDT CrI3 as illustrated in Fig. 2a. Scanning NV magnetometry exploits the Zeeman effect to quantitatively detect local magnetic stray fields longitudinal to the NV spin axis21. The magnitude of the magnetic field is directly related to the splitting of NV spin energies, which can be readout by optically detected magnetic resonance measurements2,18 (see Supplementary Information Note 2 for details). The spatial resolution of scanning NV magnetometry is primarily determined by the NV-to-sample distance30, which is ~70ânm in our measurements (see Supplementary Information Note 3 for details). In the current study, we report scanning NV quantum sensing measurements of a total of three twisted vdW magnet samples: 0.15° tDT CrI3, 0.25° tDT CrI3, and 15° tDT CrI3. For the brevity of our narrative, 0.15° and 0.25° refer to the small-twist angle, and 15° refers to the large-twist angle in our description.
Figure 2b shows an optical microscope image of a prepared tDT CrI3 sample with a calibrated local twist angle \(\alpha\)â=â0.15° and an expected moiré period of ~250ânm. Note that in order to visualize nanoscale magnetic patterns within individual moiré supercells, here \(\alpha\) is chosen to be smaller than ~0.5°, ensuring that the resulting moiré period is larger than our NV spatial sensitivity (~70ânm). Figure 2c presents a stray field BF map measured on a selected sample area of the 0.15° tDT CrI3 device at 2âK. An external magnetic field of ~2000âG is applied along the NV spin axis in this measurement. The tDT CrI3 sample shows clear multidomain features with stray fields of opposite polarity emanating from individual domains. The reconstructed out-of-plane magnetization (mz)2,18 map (Fig. 2d) manifests alternating FM and AFM patches on a length scale comparable with the moiré period (see Supplementary Information Note 4 for details). The local net magnetization is measured to be 0 and ~30âμB/nm2 for the AFM and FM domains, respectively, exhibiting the key feature of stacking-induced co-existing magnetic phases of moiré magnetism1,2 (see Supplementary Information Note 5 for details). The measured local magnetization of FM order in 0.15° tDT CrI3 also agrees with the theoretical value 29.4âμB/nm2,18. By performing an autocorrelation operation2 on the magnetization map, the periodic hexagonal shaped magnetic patterns of tDT CrI3 is revealed in Fig. 2e, from which the local mean twist angle (\(\alpha\)â=â0.15°) and moiré period (~250ânm) are confirmed. The notable alternating FM-AFM state is also observed in another 0.25° tDT CrI3 device (Fig. 2fâh). It is worth mentioning that a larger twist angle naturally results in spatially more compact magnetic moiré patterns (Fig. 2g) with a reduced moiré period of ~150ânm (Fig. 2h).
We now present systematic scanning NV magnetometry measurements to show distinct magnetic phase transition temperatures for the observed FM and AFM states within individual moiré supercells. Figure 3a presents a zoomed-in magnetization map of a selected sample area of the 0.15° tDT CrI3 device, showing two neighboring nanoscale FM and AFM domains. Due to the fully compensated net magnetic moment, temperature driven second-order magnetic phase transition of the AFM domain in tDT CrI3 is challenging to access by measuring its emanating magnetic flux. To circumvent this issue, here we employ NV relaxometry method27,31,32,33,34,35 to probe the intrinsic spin fluctuations at the local AFM (FM) regions. Figure 3b illustrates the mechanism of scanning NV relaxometry measurements, which takes advantage of the dipole-dipole interaction between local spin fluctuations of AFM (FM) spin density and a proximal NV center contained in a diamond cantilever. Spin fluctuations in a magnetically correlated system are driven by its time-dependent spin density distribution, for example, due to the dynamic imbalance in the thermal occupation of magnon bands with opposite chiralities31,34,36. For both AFM and FM systems with (un)compensated net static magnetic moment, spin-spin correlation-induced time-dependent fluctuations of the average spin density do not vanish and are expected to reach a maximum intensity around the magnetic phase transition points33,34. The emanating fluctuating magnetic fields at the NV electron spin resonance (ESR) frequencies will induce NV spin transitions from the msâ=â0 to msâ=â\(\pm\)1 state, resulting in enhancement of the corresponding NV spin relaxation rates27,32,33,36. By measuring the spin-dependent NV photoluminescence, the occupation probabilities of NV spin states can be quantitatively obtained, allowing for extraction of NV spin relaxation rate that is proportional to the magnitude of the local fluctuating magnetic fields transverse to the NV axis (see Supplementary Information Note 6 for details)27,32,33,34. Figure 3c shows the measured temperature-dependent NV spin relaxation rate \(\Gamma\) when the NV center is positioned right above the FM (AFM) domain formed in the 0.15° tDT CrI3 sample (Fig. 3a). Due to the divergent magnetic susceptibility, the measured NV spin relaxation rate shows a clear enhancement across the second-order magnetic phase transition points of tDT CrI3, from which the Curie and Néel temperatures of the co-existing FM-AFM states are measured to be ~58âK and ~48âK, respectively. We would like to highlight that similar experimental signatures are also observed in the 0.25° tDT CrI3 sample as detailed in Supplementary Information Note 7. It is worth noting that twist-induced lattice reconstructions for small-twist angles typically happen at the boundary of different stacking areas within individual moiré supercells in tDT CrI337, so the potential strain effect would play a marginal role in our measurements here.
At the monoclinic (ABâ) stacking sites, tDT CrI3 features the same type AFM interlayer interaction across the six CrI3 layers, thus it is not surprising that the formed local AFM moment manifests a similar Tc with that of the atomically thin pristine CrI3 crystals studied in previous work (see Supplementary Information Note 8 for details)8,17,38. In contrast, the FM order established in tDT CrI3 is driven by rhombohedral (AB) stacking featuring an enhanced, negative interlayer exchange coupling at the twisted interface8,39. Intuitively, such a stronger magnetic interaction will render the magnetic ordering more robust against thermal perturbations, resulting in an increased Tc of the local FM state formed in tDT CrI3. It is instructive to note that the Tc (or Tc-equivalent) defined in the current manuscript for the FM order in tDT CrI3, an âinhomogeneousâ magnetic system along the thickness direction, has considered the overall magnetic contributions from all the six CrI3 monolayers. Our one-dimensional (1D) scanning NV relaxometry measurements across the FM-AFM domains (Fig. 3d) further confirm this point. When Tâ=â48âK, one can see that the measured 1D NV spin relaxation spectrum shows a peak value at the corresponding AFM domain site. As the temperature increases to 58âK, the observed peak of NV relaxation rate shifts to the FM domain side, demonstrating that the co-existing FM-AFM states show distinct magnetic phase transition temperatures. Note that such an experimental feature is absent in NV relaxation measurements performed at 38âK and 65âK, when the temperature is away from the Curie (Néel) points of FM (AFM) domains in 0.15° tDT CrI3.
After showing the nanoscale stacking engineered Tc of moiré magnetism, next, we present temperature-dependent magnetization maps to provide an alternative perspective to examine the second-order magnetic phase transitions in tDT CrI3. Figure 4a, b presents reconstructed out-of-plane magnetization (mz) maps of the 0.15° and 0.25° tDT CrI3 devices measured at 38âK. It is evident that rhombohedral (AB) stacking induced FM state in tDT CrI3 features uncompensated net magnetization. As the temperature increases, the FM moment in tDT CrI3 gradually decreases and remains robust at 46âK (Fig. 4d, e). Figure 4g plots the temperature dependence of the average out-of-plane magnetization \({\bar{m}}_{{{{{{\rm{z}}}}}}}\) of FM domains formed in the selected sample areas of 0.15° and 0.25° tDT CrI3 devices, from which the corresponding Curie point is obtained to be ~58âK, in agreement with our NV relaxometry results. Note that in the small-twist-angle regime, local moiré lattices relax to the monoclinic (ABâ) and rhombohedral (AB) stacking sites whose magnetic ground states follow their naturally preferred state. Under this condition, the interlayer exchange energy in tDT CrI3 is dominated by the local stacking order, and the (small) twist angle plays a secondary role. Thus, the interlayer exchange energy of 0.15° and 0.25° tDT CrI3 are basically the same, leading to the (almost) identical temperature-dependent magnetic phase transition behaviors. To further investigate stacking order-dependent magnetic response in tDT CrI3, Fig. 4c presents a magnetization map of a 15° tDT CrI3 device measured at 38âK. Notably, the co-existing FM-AFM phase disappears while a pure collinear FM ground state emerges in the large-twist-angle regime1. It is worth mentioning that the two CrI3 trilayers are weakly ferromagnetically coupled at the twisted interface of 15° tDT CrI3, resulting in a reduced Curie temperature (~48âK) in comparison with that of the small-twist-angle tDT CrI3. One can see that the 15° tDT CrI3 sample enters the paramagnetic phase showing zero net FM moment at 48âK (Fig. 4f). To better illustrate this point, Fig. 4h presents the histograms of Curie temperatures of individual FM domains formed in 0.15°, 0.25°, and 15° tDT CrI3 devices. Statistically, it is evident that FM domains in small-twist-angle tDT CrI3 show clearly higher Curie temperatures than their counterparts in large-twist-angle tDT CrI3 (see Supplementary Information Note 9 for details).
We now present a mean-field theoretical model of atomically layer-resolved magnetic order to explain the nanoscale stacking-dependent magnetic phase transitions observed in small-twist-angle tDT CrI3. Our model consists of layer-uniform Ising spins representing out-of-plane magnetic moments of Cr atoms in tDT CrI3 as shown in Fig. 5a. The intralayer magnetic exchange interaction Ji, stacking-dependent monoclinic (ABâ) type interlayer exchange interaction JM, and rhombohedral (AB) type interlayer exchange interaction JR in tDT CrI3 are obtained to be \(-\)1.32âmeV, 0.086âmeV, and \(-\)0.99âmeV by comparing mean-field theory results with existing experimental transition temperatures and/or theoretically predicted values (See Method Section and Supplementary Information Note 10 for details)40. We compare the mean-field phases of tDT CrI3 with uniform monoclinic stacking (bottom panel of Fig. 5a) and rhombohedral stacking at the twisted interface (top panel of Fig. 5a), approximating the situation for regions deep inside different stacking domains. The case of monoclinic stacking features the same exchange coupling between all neighboring layers and should, therefore, show a single Néel transition at the temperature scale set by Ji and JM. However, for the rhombohedral stacking case, the much stronger JR (in comparison with JM) between the two middle CrI3 layers (layer 3 and layer 4 shown in Fig. 5a) is expected to drive the system through an FM transition at a higher temperature. To corroborate the above physical picture, we have solved the mean-field equations of our model to get the temperature dependence of the normalized out-of-plane magnetization of individual CrI3 layers as shown in Fig. 5b, c (See Supplementary Information Note 10 for details). In the rhombohedral (AB) stacking region (Fig. 5b), one can see that the middle CrI3 layer (layer 3) indeed exhibits a higher Tcâ~â58âK in comparison with that of ~48âK for the monoclinic (ABâ) stacking case (Fig. 5c). Moreover, in the AB stacking case the mean-field magnetic moments of the two top CrI3 monolayers (layer 1 and layer 2) undergo a smooth crossover as temperature decreases below Tc before reaching the saturated value but not another phase transition as dictated by the Lee-Yang theorem41. By considering the overall contributions from all the six CrI3 layers, the effective Curie temperature (Tc-equivalent) of stacking induced FM order in small-twist-angle tDT CrI3 is calculated to be 58âK, in agreement with our NV measurement results. In contrast, in the monoclinic (ABâ) stacking case (Fig. 5c), the top and middle CrI3 monolayers (layer 1 and layer 3 in Fig. 5a) exhibit the same temperature dependence with a Tc of ~48âK.
To further highlight the role of stacking engineering in affecting the local magnetic order and transition temperatures of small-twist-angle tDT CrI3, Fig. 5d plots a mean-field phase diagram of the normalized out-of-plane magnetization of CrI3 layer 4 (mz4) as a function of temperature T and Jtwist/JM using our model. Here, Jtwist is a variable interlayer exchange interaction at the twisted interface of tDT CrI3, and the sign of mz4 is defined relative to the direction of the out-of-plane magnetization mz3 of the CrI3 layer 3. When Jtwist/JM \( > \) 0, the system orders antiferromagnetically layer-wise below the Tc and mz4 is oppositely aligned with mz3 as shown in Fig. 5a. When Jtwist/JM \( < \) 0, the two twist-interfaced CrI3 monolayers are ferromagnetically coupled and uncompensated net magnetic moment is formed in tDT CrI3 below the magnetic critical temperature. Notably, as the magnitude of the FM-like coupling Jtwist increases, the calculated Tc enhances from 48âK to 58âK when Jtwist/JM reaches \(-\)11.5, corresponding to the rhombohedral stacking case (âJtwistâ=âJR). In the high-temperature regime, small-twist-angle tDT CrI3 enters the paramagnetic phase where the long-range FM (AFM) order vanishes.
Discussion
In summary, we have utilized scanning NV magnetometry techniques to investigate spatial and thermodynamic phase separation of moiré magnetism hosted by twisted CrI3. By using NV spin relaxometry methods to probe local spin fluctuations, we explicitly show that the co-existing FM-AFM phases within individual moiré supercells manifest distinct second-order magnetic phase transition points. The measured Curie temperature of the rhombohedral (AB) stacking-driven FM state is significantly higher than the Néel temperature of monoclinic (ABâ) stacking-driven AFM due to the spatially modulated exchange interaction at the twisted interface. In addition, we have directly visualized the stray field and magnetization maps across the magnetic phase transition points showing that twist engineering can effectively control the Curie temperature of local FM order in tDT CrI3. Our results are well rationalized by a proposed mean-field theoretical model, which captures the layer-resolved Tc in small-twist-angle tDT CrI3. When extending to other combinations of atomically thin twisted odd number of CrI3 layers, we expect that the observed stacking order-dependent magnetic phase transition with separate critical temperatures remains observable by the scanning NV quantum microscopy techniques. The presented results further highlight the opportunities provided by quantum spin sensors for investigating the local spin-related phenomena in moiré quantum magnets. The stacking and temperature-driven âintracellâ magnetic phase separations observed in twisted 2D magnets may also find relevant applications in realizing local control of functional material properties through carefully engineered proximity effect, advancing the current state of the art of vdW spintronic devices24,42.
Methods
Materials and device fabrications
CrI3 crystals used in this study were grown by the chemical vapor transport method, as reported in previous literature3. The tDT CrI3 devices are fabricated by the âtear-and-stackâ method and encapsulated by hBN nanoflakes. The samples for SAED/TEM and scanning NV measurements were fabricated separately as described below, and different samples were used for SAED/TEM and NV studies. We first exfoliated the bulk crystals onto SiO2/Si substrates to obtain trilayer CrI3 and few-layer hBN. Next, the top hBN and one part of trilayer CrI3 were picked up by a poly(bisphenol A carbonate) stamp by sequence. The other part of trilayer CrI3 remained on the Si/SiO2 substrate and was rotated by a well-controlled angle and then picked up. The two CrI3 trilayer flakes were stacked with each other to form a twisted device and finally encapsulated by the bottom hBN flake. The layer number of atomically thin CrI3 flakes was determined by thickness-dependent optical contrast and confirmed by magnetic circular dichroism measurements5. The entire 2D device fabrication processes were performed inside a nitrogen-filled glovebox with oxygen level below 0.1âppm and water level below 0.5âppm.
For scanning NV measurements, the final hBN/tDT CrI3/hBN vdW stack was released onto quartz (SiO2) substrates with prepatterned Au striplines. The quartz substrate was chosen due to its ideal insulating nature and dielectric properties for microwave performance. The Au stripline was utilized to deliver on-chip microwave currents to realize local control of the NV spin sensor. We have prepared multiple tDT CrI3 samples to ensure the consistency of the presented NV results. In this paper, we reported scanning NV magnetometry studies of four CrI3 devices including 0.15° tDT CrI3, 0.25° tDT CrI3, 15° tDT CrI3, and a pristine trilayer CrI3 control sample (presented in Supplementary Information). Note that samples with other layer thicknesses such as small-angle twisted double bilayer CrI3 and twisted bilayer CrI3, are not studied here due to the lack of net ferromagnetism or potential device quality issues. For SAED/TEM measurements, prepared hBN/tDT CrI3/hBN stacks were transferred onto TEM grids with a 10-nm-thick SiN membrane. TEM images of tDT CrI3 with different local twist angles are presented in Supplementary Information Note 2.
Scanning NV magnetometry measurements
The presented quantum sensing measurements were performed using a scanning NV magnetometry system consisting of a home-built confocal and a custom-designed atomic force microscope (AFM) operating in a cryo-free cryostat (attocube Inc.). A commercially available diamond cantilever (Q-Zabre LLC) containing individually addressable NV centres was glued to a quartz tuning fork for force-feedback AFM operations, and a window on top of the cryostat provided optical access for NV measurements. A sample holder with coplanar waveguides was fixed onto a stack of piezo-based positioners and scanners to engage with the diamond cantilever and perform 2D scanning measurements. We applied continuous green laser and microwave signals to carry out NV optically detected magnetic resonance measurements36. NV spin states were addressed by measuring NV photoluminescence using an avalanche photodiode. Microwave signals delivered to the on-chip Au stripline were supplied by a Stanford Research Systems SG386 signal generator, and the external magnetic field applied in NV measurements was generated by a three-axis superconducting vector magnet.
NV spin relaxometry measurements
Pulsed NV spin relaxometry measurements were performed using the scanning NV magnetometry system presented above. An external magnetic field of ~1200âG was applied along the NV spin axis in these measurements, and the corresponding NV ESR frequency for the msâ=â0 to msâ=â\(-\)1 spin transition is ~0.5âGHz. Green-laser pulses used for NV initialization and readout were generated by an electrically driven 515-nm laser. The laser power entering the objective was ~0.5âmW. The trigger pulses to the optical modulator and photon counter were generated by a programmable pulse generator, and microwave signals were modulated by a switch (Minicircuits ZASWA-2-50DR+). The top panel of Supplementary Fig. 7b shows the details of the pulsed optical and microwave sequences for our NV spin relaxometry measurements. A 1.5-μs-long green laser pulse was first applied to initialize the NV spin to the msâ=â0 state. After a delay time t, we measured the occupation probabilities of the NV spin at the msâ=â0 and msâ=â\(-\)1 state by applying a microwave \(\pi\) pulse on the corresponding NV ESR frequencies and measuring the spin-dependent NV photoluminescence during the first ~600âns of the green-laser readout pulse. By measuring the integrated photoluminescence intensity as a function of the delay time t and fitting the data with a three-level model (See Supplementary Information Note 6 for details)36, NV spin relaxation rates can be quantitatively measured.
SAED and TEM measurements
We utilized Thermo Fisher Talos, operated at 200âkV and equipped with Gatan OneView camera, to perform the SAED and TEM measurements to characterize microscopic lattice structures of tDT CrI3 devices encapsulated by hBN nanoflakes. The local mean twist angle of the sample was obtained by fitting 2D Gaussians to the Bragg peaks. Real-space structural image (Fig. 1c) was acquired by averaging DF-TEM images from three fifth-order Bragg peaks 120° apart. For tDT CrI3, the moiré domains appear â3 smaller in the TEM image(s) in comparison to the expected domain size for the corresponding twist angle \(\alpha\).
Theoretical modeling
We used mean-field theory and minimal Ising models to investigate layer-resolved magnetic phases of small-twist-angle tDT CrI3 with different stacking orders. Each Ising spin, representing a Cr ion, has three intralayer nearest neighbors with exchange coupling \({J}_{{{{{{\rm{i}}}}}}}\) and one nearest neighbor in an adjacent layer with interlayer exchange coupling \({J}_{{{{{{\rm{o}}}}}}}\). The Hamiltonian of the system is given by:
where \(j,l\) label sites within each layer and \(a,b\) label layers, \({\sigma }_{{ja}}=\pm\)1, \({J}_{{{{{{\rm{i}}}}}}} \, < \, 0\) (FM), and \({J}_{{{{{{\rm{o}}}}}},{ab}}\) can be negative (AB or rhombohedral stacking, denoted by \({J}_{{{{{{\rm{R}}}}}}}\) below) or positive (ABâ or monoclinic stacking, denoted by \({J}_{{{{{{\rm{M}}}}}}}\) below) depending on the stacking order. Since we do not consider magnetic ordering beyond the FM order within each layer, \({\sigma }_{{ja}}\equiv {m}_{{za}}\) is independent of \(j\). We fix \({J}_{{{{{{\rm{i}}}}}}}\), \({J}_{{{{{{\rm{R}}}}}}}\), and \({J}_{{{{{{\rm{M}}}}}}}\) by requiring the mean-field \({T}_{{{{{{\rm{c}}}}}}}\) of the Ising model to be consistent with the experimental (or theoretically predicted) values, from which the following exchange couplings: \({J}_{{{{{{\rm{i}}}}}}}\)â=â\(-\)1.32âmeV, \({J}_{{{{{{\rm{R}}}}}}}\)â=â\(-\)0.99âmeV, and \({J}_{{{{{{\rm{M}}}}}}}\)â=â0.086âmeV can be obtained. To mimic the local stacking order of small-twist-angle tDT CrI3, we have considered different 6-layer models with ABâ stacking between adjacent layers in the top 3 (labeled by layers 1, 2, 3) and bottom 3 (labeled by layers 4, 5, 6), but either AB or ABâ stacking between layers 3 and 4. Temperature-dependent mean-field values of the layer-resolved Ising spins are solved numerically from the following nonlinear mean-field equation:
Data availability
All data supporting the findings of this study are available from the corresponding author upon reasonable request.
References
Xu, Y. et al. Coexisting ferromagneticâantiferromagnetic state in twisted bilayer CrI3. Nat. Nanotechnol. 17, 143â147 (2022).
Song, T. et al. Direct visualization of magnetic domains and moiré magnetism in twisted 2D magnets. Science 374, 1140â1144 (2021).
Xie, H. et al. Twist engineering of the two-dimensional magnetism in double bilayer chromium triiodide homostructures. Nat. Phys. 18, 30â36 (2022).
Cheng, G. et al. Electrically tunable moiré magnetism in twisted double bilayers of chromium triiodide. Nat. Electron. 6, 434â442 (2023).
Xie, H. et al. Evidence of non-collinear spin texture in magnetic moiré superlattices. Nat. Phys. 19, 1150â1155 (2023).
Tong, Q., Liu, F., Xiao, J. & Yao, W. Skyrmions in the Moiré of van der Waals 2D magnets. Nano. Lett. 18, 7194â7199 (2018).
Akram, M. et al. Moiré skyrmions and chiral magnetic phases in twisted CrX3 (X = I, Br, and Cl) bilayers. Nano. Lett. 21, 6633â6639 (2021).
Sivadas, N. et al. Stacking-dependent magnetism in bilayer CrI3. Nano. Lett. 18, 7658â7664 (2018).
Wang, C., Gao, Y., Lv, H., Xu, X. & Xiao, D. Stacking domain wall magnons in twisted van der Waals magnets. Phys. Rev. Lett. 125, 247201 (2020).
Hejazi, K., Luo, Z.-X. & Balents, L. Noncollinear phases in moiré magnets. Proc. Natl Acad. Sci. 117, 10721â10726 (2020).
Huang, M. et al. Revealing intrinsic domains and fluctuations of moiré magnetism by a wide- field quantum microscope. Nat. Commun. 14, 5259 (2023).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43â50 (2018).
Balents, L., Dean, C. R., Efetov, D. K. & Young, A. F. Superconductivity and strong correlations in moiré flat bands. Nat. Phys. 16, 725â733 (2020).
Andrei, E. Y. et al. The marvels of moiré materials. Nat. Rev. Mater. 6, 201â206 (2021).
Kennes, D. M. et al. Moiré heterostructures as a condensed-matter quantum simulator. Nat. Phys. 17, 155â163 (2021).
Wu, F., Lovorn, T., Tutuc, E. & MacDonald, A. H. Hubbard model physics in transition metal dichalcogenide moiré bands. Phys. Rev. Lett. 121, 026402 (2018).
Huang, B. et al. Layer-dependent ferromagnetism in a van der Waals crystal down to the monolayer limit. Nature 546, 270â273 (2017).
Thiel, L. et al. Probing magnetism in 2D materials at the nanoscale with single-spin microscopy. Science 364, 973â976 (2019).
Niu, B. et al. Coexistence of magnetic orders in two-dimensional magnet CrI3. Nano. Lett. 20, 553â558 (2020).
Degen, C. L., Reinhard, F. & Cappellaro, P. Quantum sensing. Rev. Mod. Phys. 89, 035002 (2017).
Rondin, L. et al. Magnetometry with nitrogen-vacancy defects in diamond. Rep. Prog. Phys. 77, 056503 (2014).
Doherty, M. W. et al. The nitrogen-vacancy colour centre in diamond. Phys. Rep. 528, 1â45 (2013).
Geim, A. K. & Grigorieva, I. V. Van der Waals heterostructures. Nature 499, 419â425 (2013).
Burch, K. S., Mandrus, D. & Park, J.-G. Magnetism in two-dimensional van der Waals materials. Nature 563, 47â52 (2018).
Kim, K. et al. van der Waals heterostructures with high accuracy rotational alignment. Nano. Lett. 16, 1989â1995 (2016).
Sung, S. H. et al. Torsional periodic lattice distortions and diffraction of twisted 2D materials. Nat. Commun. 13, 7826 (2022).
Finco, A. et al. Imaging non-collinear antiferromagnetic textures via single spin relaxometry. Nat. Commun. 12, 767 (2021).
Ku, M. J. H. et al. Imaging viscous flow of the Dirac fluid in graphene. Nature 583, 537â541 (2020).
Pelliccione, M. et al. Scanned probe imaging of nanoscale magnetism at cryogenic temperatures with a single-spin quantum sensor. Nat. Nanotechnol. 11, 700â705 (2016).
Casola, F., van der Sar, T. & Yacoby, A. Probing condensed matter physics with magnetometry based on nitrogen-vacancy centres in diamond. Nat. Rev. Mater. 3, 17088 (2018).
McCullian, B. A. et al. Broadband multi-magnon relaxometry using a quantum spin sensor for high frequency ferromagnetic dynamics sensing. Nat. Commun. 11, 5229 (2020).
Andersen, T. I. et al. Electron-phonon instability in graphene revealed by global and local noise probes. Science 364, 154â157 (2019).
Flebus, B. & Tserkovnyak, Y. Quantum-impurity relaxometry of magnetization dynamics. Phys. Rev. Lett. 121, 187204 (2018).
McLaughlin, N. J. et al. Quantum imaging of magnetic phase transitions and spin fluctuations in intrinsic magnetic topological nanoflakes. Nano. Lett. 22, 5810â5817 (2022).
Ariyaratne, A., Bluvstein, D., Myers, B. A. & Jayich, A. C. B. Nanoscale electrical conductivity imaging using a nitrogen-vacancy center in diamond. Nat. Commun. 9, 2406 (2018).
Du, C. et al. Control and local measurement of the spin chemical potential in a magnetic insulator. Science 357, 195â198 (2017).
Yoo, H. et al. Atomic and electronic reconstruction at the van der Waals interface in twisted bilayer graphene. Nat. Mater. 18, 448â453 (2019).
Jin, W. et al. Raman fingerprint of two terahertz spin wave branches in a two-dimensional honeycomb Ising ferromagnet. Nat. Commun. 9, 5122 (2018).
Chen, L. et al. Topological spin excitations in honeycomb ferromagnet CrI3. Phys. Rev. X 8, 041028 (2018).
Stanley, H. E. Introduction to phase transitions and critical phenomena (Oxford University Press, 1971).
Lee, T. D. & Yang, C. N. Statistical theory of equations of state and phase transitions. II. lattice gas and Ising model. Phys. Rev. 87, 410â419 (1952).
ŽutiÄ, I., Fabian, J. & Das Sarma, S. Spintronics: fundamentals and applications. Rev. Mod. Phys. 76, 323â410 (2004).
Acknowledgements
The authors are grateful to Rainer Stöhr and Sreehari Jayaram for valuable discussions on cryogenic NV measurement techniques. This work was primarily supported by the U.S. Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under award No. DE-SC0024870. Development of the cryogenic scanning NV microscopy was supported by the Air Force Office of Scientific Research (AFOSR) under Grant no. FA9550-20-1-0319 and its Young Investigator Program under Grant no. FA9550-21-1-0125. C.R.D. also acknowledges the support from the Office of Naval Research (ONR) under Grant no. N00014-23-1-2146. L.Z. acknowledges the support from the U.S. Department of Energy (DOE), Office of Science, Basic Energy Science (BES), under award No. DE-SC0024145 (for the moiré magnets sample fabrication) and AFOSR YIP Grant no. FA9550-21-1-0065 and Alfred. P. Sloan Foundation (for the instrument development of 2D vdW structure fabrication platform). A.S. and H.C. acknowledge the support from the National Science Foundation CAREER Grant no. DMR-1945023. R.H. acknowledges support from the Department of Energy, Basic Energy Sciences (DE-SC0024147). H.C.L. was supported by the National Key R&D Program of China (Grant no. 2018YFE0202600 and 2022YFA1403800), the Beijing Natural Science Foundation (Grant no. Z200005), the National Natural Science Foundation of China (Grants no. 12174443), and the Beijing National Laboratory for Condensed Matter Physics.
Author information
Authors and Affiliations
Contributions
S.L. performed the NV measurements and analyzed the data with N.M., M.H., H.L., and H.W.Z.S., and L.Z. prepared the CrI3 devices. A.S. and H.C. performed the mean-field theoretical calculations. N.A., S.H.S., and R.H. performed the transmission electron microscopy characterizations. S.Y. and H.C.L. provided bulk CrI3 crystals. C.R.D. supervised this project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Li, S., Sun, Z., McLaughlin, N.J. et al. Observation of stacking engineered magnetic phase transitions within moiré supercells of twisted van der Waals magnets. Nat Commun 15, 5712 (2024). https://doi.org/10.1038/s41467-024-49942-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-49942-2
This article is cited by
-
Nanoscale diamond quantum sensors for many-body physics
Nature Reviews Physics (2024)