Thickness-dependent magnetic order in CrI3 single crystals.

Two-dimensional (2D) materials with intrinsic ferromagnetism provide unique opportunity to engineer new functionalities in nano-spintronics. One such material is CrI3, showing long-range magnetic order in monolayer with the Curie temperature (Tc) of 45 K. Here we study detailed evolution of magnetic transition and magnetic critical properties in response to systematic reduction in crystal thickness down to 50 nm. Bulk Tc of 61 K is gradually suppressed to 57 K, however, the satellite transition at T * = 45 K is observed layer-independent at fixed magnetic field of 1 kOe. The origin of T * is proposed to be a crossover from pinning to depinning of magnetic domain walls. The reduction of thickness facilitates a field-driven metamagnetic transition around 20 kOe with out-of-plane field, in contrast to the continuous changes with in-plane field. The critical analysis around Tc elucidates the mean-field type interactions in microscale-thick CrI3.

Introduction

Layered materials, when thinned to atomic limits in 2D exhibit novel properties, different from the bulk counterpart. Recent discoveries of intrinsic 2D ferromagnetism in atomically thin CrI3 and Cr2Ge2Te6 open up new opportunities for studying fundamental 2D magnetism and show great potential in spintronic applications[1-3]. CrI3 and Cr2Ge2Te6 are ferromagnetic (FM) below T ~ 61 K in bulk[4-6]. First-principle calculations predict a robust 2D ferromagnetism with T ~ 57.2 or 106 K in monolayer Cr2Ge2Te6[7,8], However, the scanning magneto-optic Kerr microscopy experiment shows that the T monotonically decreases with decreasing thickness, showing 30 K in bilayer and the absence of T in monolayer[3].

In bulk CrI3, the Cr atoms in each layer form a honeycomb structure, and each Cr atom is surrounded by six I atoms in an octahedral coordination [Fig. 1(a)]. In contrast to Cr2Ge2Te6, ferromagnetism in CrI3 persists in monolayer with T of 45 K[2]. Intriguingly, the magnetism in CrI3 is layer-dependent, from FM in monolayer to antiferromagnetic (AFM) in bilayer, and back to FM in trilayer, providing great opportunities for designing magneto-optoelectronic devices[2]. The monolayer CrI3 can be described by the Ising model[2,9]. Remarkably, the electrostatic doping can modify the saturation magnetization, the coercive force, and the T in monolayer as well as change the interlayer magnetic order in bilayer CrI3[10,11]. In addition, giant tunneling magnetoresistance is observed in few-layer CrI3, exhibiting multiple states as a function of magnetic field[12-15]. This highlights CrI3 as a potential magnetic tunnel barrier for van der Waals heterostructure spintronic devices and also shows demand for magnetic order investigation at all length scales from bulk to monolayer.

Figure 1

(a)Low temperature structure of CrI3 in the space group. (b) Optical microscopy of bulk S0 sample (thickness ~100 μm) and cross-section of S2 sample (thickness ~2 μm). (c) Transmission electron microscopy (TEM) image of S3 sample in Cu grid along with the low-loss and core-loss electron energy loss spectra (EELS) taken from top right area marked by red circle. The thickness of S3 was calculated from the low-loss EELS based on log-ratio method in Digital Micrograph program to be 44 nm. The average thickness for many areas measured is 49 ± 8 nm.

Here we report a detailed study of thickness-dependent magnetic transition and critical behavior in CrI3. An additional satellite transition at T* was observed just below T. The value of T* is found thickness-independent at fixed magnetic field, in contrast to gradually suppressed T with reduction in thickness. Thin CrI3 crystals show an increase in coercive field and distinct field-driven metamagnetic transition around 20 kOe with out-of-plane field but not in-plane field. The critical behavior suggests the long-range mean-field type interactions in mesoscale-thick CrI3 crystals.

Methods

Bulk CrI3 single crystals were fabricated by the chemical vapor transport method and characterized as described previously[16]. A series of thicknesses were obtained step-by-step mechanical exfoliating bulk crystal down to nanometer scale [Fig. 1(b,c)]. Considering the reactivity of thin CrI3 flakes, all the samples were prepared in an argon-filled glove-box and protected using scotch tape on both sides when transferring for magnetization measurement. The dc magnetization was measured in Quantum Design MPMS-XL5 system. An optical microscopy equipped with 100× objective lens in Witec alpha 300 confocal Raman microscope was used for imaging the cross-section of S2 sample to determine its thickness. Transmission electron microscopy (TEM) sample of S3 was prepared by putting the exfoliated S3 sample in Cu grid. The thickness (t) of S3 sample was measured by the low loss EELS spectrum in JEM-ARM200F microscope [Fig. 1(c)], t = ln(I/I0) × λ, where I0 and I are integrated EELS signal intensity under the zero-loss peak and the total spectrum, respectively; λ is a total mean free path for all inelastic scattering which can be calculated based on the composition of the sample[17].

Results and Discussion

Figure 2(a–d) present the temperature-dependent magnetization M(T) with zero-field cooling (ZFC) and field-cooling (FC) modes measured in out-of-plane field H = 1 kOe for the samples (S0–S3) with typical thickness from 100 μm to 50 nm. There is an apparent increase in M(T), in line with the paramagnetic (PM) to FM transition[1]. The T can be determined from the minima of the first derivative of the M(T) curves [insets in Fig. 2(a–d)]. Given the large magnetocrystalline anisotropy in CrI3[18], the divergence of ZFC and FC curves is observed below T in H = 1 kOe. With subsequent reduction of crystal thickness this discrepancy tends to be larger, indicating an increasing anisotropy in thin crystals. Intriguingly, an additional satellite transition T* just below T is observed. A similar phenomenon was also observed in Fe3GeTe2, arising from the emergence of antiparallel spin arrangement along the c axis between different Fe3Ge layers[19]. This satellite transition T* is defined by the maxima of the first derivative of the M(T) curves [insets in Fig. 2(b–d)], which is more apparent in thinner crystal. The thickness dependence of T and T* is summarized in Fig. 2(e). The T shows linear decrease with logarithm of thickness. In contrast, the T* seems thickness-independent and shows a value of 45 K identical to the T in monolayer measured by magneto-optic Kerr microscopy in H = 1.5 kOe[2]. The isothermal magnetization at T = 2 K for the corresponding samples is shown in Fig. 2(f–i). With H//c, the magnetization of bulk CrI3 crystal saturates at a relatively low magnetic field H ~ 2 kOe and in line with the previous reports[1,16]. With reduction of thickness, there is an increase of coercive field H as shown in insets of Fig. 2(f–i), from about 80 Oe for S0 sample to 920 Oe for S3 sample, still identified as a relatively soft ferromagnet. It is interesting to note that an apparent magnetization jump around H ~ 20 kOe is observed in thin crystals. The calculated value δ = [M(25kOe) − M(20kOe)]/M(20kOe) increases from 5% for S1 to 43% for S3. This field-driven metamagnetic transition is first identified by the step-like change of magnetization in the hysteresis curves, which is similar with the giant tunneling magnetoresistance step around the same field observed in few-layer devices[15]. The tunneling magnetoresistance step in few-layer devices is attributed to spin filtering effect arising from collapse of interlayer AFM in thin CrI3[12-15]. Above this critical field, a fully spin polarized FM state is reached in S0–S2 samples, however, a large hysteresis is observed in S3 sample.

Figure 2

(a–d) Temperature dependence of magnetization for CrI3 at various crystal thicknesses measured with H//c at 1 kOe. Insets show the first derivative of the M(T) curves. (e) Thickness-dependent T* and T. (f–i) Isothermal magnetization at T = 2 K for corresponding S0-S3 samples. Insets show the enlarged low field part and the derived coercive field H. (j) Thickness dependence of H.

To well understand the nature of magnetism in CrI3 and its properties from bulk to monolayer, we then focus on the S2 sample with mesoscale thickness of 2 μm. Figure 3 shows the temperature- and field-dependent magnetization M(T,H) with out-of-plane and in-plane fields, respectively. The M(T) is nearly isotropic in H = 50 kOe, whilst anisotropic magnetic response is observed in low fields. In H = 10 kOe, a monotonic increase of M(T) with decreasing temperature is observed for H//c [Fig. 3(a)], however, the M(T) decreases with decreasing temperature below 50 K for H//ab [Fig. 3(d)]. Similar feature was also observed in bulk CrI3 crystal, associated with its temperature dependent magnetocrystalline anisotropy[18,20]. The satellite transition T* reveals itself by the emergence of a weak kink in the dM/dT curve with H = 10 kOe. With decreasing field, the magnetic anisotropy results in larger divergence between ZFC and FC magnetization and a more abrupt satellite transition, as shown in Fig. 3(b,e). Both T and T* are field-dependent but with an opposite tendency, i.e., with decreasing field the T decreases whereas the T* gradually increases. The isothermal magnetization at 2 K increases discretely for out-of-plane field [Fig. 3(c)], but monotonically for in-plane field [Fig. 3(f)]. Although an A-type AFM ground state was proposed in few-layer devices[15], a FM ground state is still expected in S2 sample with thickness ~2 μm. Therefore, the field-driven metamagnetic transition is most likely caused by a movement or depinning of magnetic domains, rather than flipping the spins from an antiparallel to a parallel configuration. A crossover from pinning to depinning of magnetic domain walls is also proposed as the reason for the satellite T* kink. When the sample is cooled in zero field, the magnetic domains start to be pinned below the crossover temperature T*. Small applied field at base temperature is insufficient to move the pinned domains. With increasing temperature, thermal fluctuations gradually weaken the pinning force and finally completely depin the domains above T*. This also explains why the T* kink is more apparent in the ZFC curves. In the FC process the domains are always pinned with the effective FM moment aligned along the cooling field.

Figure 3

(a,b)Temperature-dependent magnetization measured at various magnetic fields and (c) field-dependent magnetization measured at T = 2 K for H//c. (d,e) Temperature-dependent magnetization measured at various magnetic fields and (f) field-dependent magnetization measured at T = 2 K for H//ab.

Near a second-order phase transition, the free energy of a ferromagnet . By ∂F(M)/∂M = 0, the equation of state is obtained as H/M = a + bM2 for the behavior of magnetization near T. According to the scaling hypothesis, the spontaneous magnetization M below T, the inverse initial susceptibility H/M above T, and the M(H) at T can be characterized by a series of critical exponents[21]:where M0, h0/m0, D and ε = (T − T)/T are the critical amplitudes and the reduced temperature, respectively[22].

The isothermal magnetization curves with out-of-plane field are depicted in Fig. 4(a,b). The magnetization jump gradually moves to lower field and finally disappears at 50 K for S2. The Arrott plot of M2 vs H/M [Fig. 4(c)], with β = 0.5 and γ = 1.0[23,24], gives a set of quasi-straight lines and the isotherm at T pass through the origin, suggesting the mean-field type magnetic interactions in S2 sample. In addition, the positive slope of the straight lines indicates that it is a second-order transition based on the Banerjee’s criterion[25]. The Arrot-Noaks equation of state gives a general modified Arrott plot (H/M)1/ = aε + bM1/, where a and b are constants[26]. Taking into consideration the three-dimensional (3D) critical behavior in bulk CrI3 (S0 sample)[16], and 2D Ising-like in monolayer[2], a set of possible exponents belonging to different models are used to build the modified Arrott plots[27]. Comparison the normalized slope (NS = S(T)/S(T)) with ideal value “1” can determine the most suitable model, where S(T) = dM1//d(H/M)1/. Figure 4(d) presents the NS vs T curves for indicated models, confirming that the mean-field model is the best and the 2D Ising model shows the largest deviation. That is to say, the critical behavior of mesoscale S2 sample is quite different when compared to bulk and monolayer[2,16].

Figure 4

(a,b) Initial isothermal magnetization at various temperatures with out-of-plane field. (c) The Arrott Plot of S2 sample with β = 0.5 and γ = 1.0. (d) Temperature-dependent NS = S(T)/S(T) for indicated models.

The linearly extrapolated M and H/M are plotted from in Fig. 5(a). The solid lines are fitted lines according to Eqs (1) and (2). The precise critical exponents β = 0.45(2), with T = 60.9(2) K, and γ = 0.90(8), with T = 60.6(2) K, are obtained, close to the parameters of mean-field model. According to Eq. (3), the M(H) at T should be a straight line in log-log scale with the slope of 1/δ [inset in Fig. 5(a)]. Such fitting yields δ = 3.12(3) for T = 60 K and 2.88(4) for T = 61 K, respectively. The Kouvel-Fisher (KF) method can be used to double check the obtained critical exponents,where M/(dM/dT) and (H/M)/[d(H/M)/dT] are linear in T and the slopes are 1/β and 1/γ, respectively[28]. Such fitting gives β = 0.41(2), with T = 60.5(3) K, and γ = 0.87(8), with T = 60.7(3) K, respectively, in line with the values obtained by the modified Arrott plot.

Figure 5

(a)Temperature-dependent M (left) and H/M (right) with solid fitting curves. Inset shows the isotherms M(H) in log-log scale near T with linear fitting curves. (b) Temperature-dependent M(dM/dT)−1 (left) and (H/M)[d(H/M)/dT)]−1 (right) with linear fitting curves. (c) Scaling plot of m vs h with m ≡ ε−M(H, ε) and h ≡ ε−(H, respectively. (d) MH−1/ vs εH−1/(.

The reliability of obtained critical exponents β, γ, δ, and T can be further checked by a scaling analysis. Near the critical region, the magnetic equation of state can be expressed aswhere f± is the regular function with f− for T < T, and f+ for T > T respectively. This equation can be further expressed as m = f± (h), where m ≡ ε−β M (H, ε) and h ≡ ε−( H are the rescaled magnetization and magnetic field, respectively. Figure 5(c) presents the scaled m vs h in log-log scale, falling on two separate branches above and below T, confirming the true scaling treatment and the intrinsic critical values of β, γ, and δ. The scaling equation of state can also be expressed aswhere the MH−1/ vs εH−1/( should collapse into a single curve with the T locating at horizontal zero point [Fig. 5(d)].

The critical exponents of CrI3 (S0 and S2) with different thickness are listed in Table 1 for comparison with the theoretical models[2,16,21,27,29]. For a 2D magnet, the value of critical exponent β should be within a window ~0.1 ≤ β ≤ 0.25[30]. As we can see, the critical exponents of bulk CrI3 crystal (S0) are between the values of theoretical tricritical mean-field and 3D Ising model, exhibiting a clear 3D critical phenomenon, in contrast to Cr2(Si,Ge)2Te6 showing 2D Ising-like behavior[31,32]. The monolayer CrI3 is claimed as a 2D Ising-like ferromagnetism as a result of large magnetic anisotropy removing the Mermin-Wagner restriction[2,33]. However, the critical exponents of β, γ, and δ for mesoscale CrI3 (S2) are well fit with the mean-field model (Table 1), which is quite different from the behavior in bulk and monolayer. It should be noted that the precise value of T for S2 almost equals to T = 61 K for S0 based on the critical analysis, in contrast to the strong thickness dependence of β, γ, and δ. We assume that the T of 45 K in monolayer is most likely stabilized by the satellite transition (T* ~ 45 K) in bulk CrI3, calling for further in-depth theoretical study.

Table 1

The critical exponents of CrI3 (S0 and S2) with different thickness obtained by various methods such as the modified Arrott plot (MAP), the Kouvel-Fisher plot (KFP), the critical isotherm (CI) and the magneto-optical Kerr (MOK) effect, in comparison with different theoretical models.

Thickness	Theoretical model	Reference	Technique	Tc (K)	β	γ	δ
S0 ~100 μm		[16]	MAP	60.45 (3)/60.38 (6)	0.284 (3)	1.146 (11)	5.04 (1)
		[16]	KFP	60.05 (13)/60.43 (4)	0.260 (4)	1.136 (6)	5.37 (4)
		[16]	CI	60			5.32 (2)
	3D Heisenberg	[29]	Theory		0.365	1.386	4.8
	3D XY	[29]	Theory		0.345	1.316	4.81
	3D Ising	[29]	Theory		0.325	1.24	4.82
	Tricritical mean field	[27]	Theory		0.25	1.0	5
S2 ~2 μm		This work	MAP	60.9 (2)/60.0 (2)	0.45 (2)	0.90 (8)	3.00 (9)
		This work	KFP	60.5 (3)/60.7 (3)	0.41 (2)	0.87 (8)	3.12 (9)
		This work	CI	60			2.88 (4)
		This work	CI	61			3.12 (3)
monolayer	Mean field	[21]	Theory		0.5	1.0	3.0
		[2]	MOK	45
	2D Ising	[27]	Theory		0.125	1.75	15

Conclusion

In summary, we have investigated the thickness-dependent magnetism in CrI3 from bulk to 50 nm thickness. In addition to suppression of bulk T, we observe an additional satellite transition T* at 45 K which is independent of thickness and corresponds to the T observed in monolayer. Analysis of critical behavior shows mean-field interactions in mesoscale-thick CrI3, in contrast to that in bulk and monolayer, which is helpful to understand the nature magnetism in in CrI3 and is of great importance for future magneto-optoelectronic devices.