SMR and XMLD measurements. We first show in Fig. 2a SMR signals of a control sample Fe2O3(12)/Pt(4) (units in nanometers), where the magnetic field (H) and current (I) is along x-axis and the spin polarization generated by the spin Hall effect of Pt is along y-axis. As expected, comparatively low resistance states at high magnetic fields reflect that the Néel vector (n) of Fe2O3 is perpendicular to H (I) due to the spin-flop at high fields and deviates towards H (I) at low fields, which is quite characteristic for negative SMR of AFM17–20. The resistance peaks at approximately µ0H = ±0.35 T owing to the deviation of n from the spin-flop state appears at a negative field as sweeping the field from positive to negative (black line), indicating that the Néel vector almost keeps the spin-flop state at zero-field16. Note that Fe2O3 with the thickness below tens of nanometers maintains easy-plane anisotropy without Morin transition16,18,20. Similar SMR signals are obtained in another control sample Cr2O3(4.4)/Fe2O3(4)/Pt(4) (Fig. 2b), where Fe2O3 was grown on a Cr2O3 buffer to ensure a closer scenario as the top Fe2O3 in the Fe2O3/Cr2O3/Fe2O3 junction. The SMR signals of the control samples are simulated and shown in Supplementary Note 2, where the hysteresis is due to the existence of Dzyaloshinskii-Moriya interaction (DMI) in Fe2O331. The antiferromagnetic Cr2O3 buffer possesses a high spin-flop field higher than 6 T9, which does not contribute to the observed SMR signals.
Figure 2c displays a representative HAADF-STEM (high-angle annular dark-field scanning transmission electron microscopy) image of the Fe2O3(12)/Cr2O3(4.4)/Fe2O3(4) cross-section, reflecting the epitaxial growth of the junction (Supplementary Note 3). Figure 2d presents the SMR curves of the Fe2O3/Cr2O3/Fe2O3 junction, which was covered by 4 nm-thick Pt. Four typical H [(i) → (iv)] are denoted in the inset. The most eminent feature is that two resistance peaks exist when sweeping H from positive to negative (black line) or reverse (red line), which is different from the SMR signals of a single Fe2O3 in Fig. 2a,b. The resistance peak appears before H = 0, which violates the principle of thermodynamics, indicating the existences of coupling effect. A low resistance is obtained for n ⊥ I (n is parallel to spin polarization) at the spin-flop state. As H is swept downward, the first resistance peak (high resistance state) at a positive H [µ0H = +0.3 T, (i)] reveals that n deviates from the spin-flop state and is unexpectedly aligned along n // H (I). This observation indicates that another effect suppresses the magnetic field effect. We attribute the overwhelming effect to the interlayer coupling between two Fe2O3 layers through the Cr2O3 spacer. The AFM coupling generates an orthogonal (90º) arrangement of n in two Fe2O3 layers. The coupling between net moment in Cr2O3 at high temperature and n in Fe2O3 is excluded (Supplementary Note 4). Based on both magnetic field (Fig. 2a,b) and angle dependent SMR measurements (Supplementary Note 5) in Fe2O3/Pt and Cr2O3/Fe2O3/Pt control samples, we find the top thinner Fe2O3 possesses a lower spin-flop field than its bottom thicker counterpart, in analogy to a soft ferromagnet with a small coercivity. Because of the relatively lower spin-flop field and smaller Zeeman energy of top thinner Fe2O3, the n in the top Fe2O3 has the priority to deviate from the spin-flop state as a result of the interlayer coupling, resulting in the resistance peak before zero-field. This is bolstered by the simulation based on calculating the energy profile of different magnetic configurations in Fig. 2e (Supplementary Note 2).
As H sweeps to the negative side, the SMR signal decreases and a resistance valley appears at negative H (ii), which is almost the same as the location of resistance peak in Fig. 2a. The direction of n in the bottom Fe2O3 is n // H, and the interlayer coupling drives the Néel vector in the top Fe2O3 to n ⊥ H (I), again giving rise to the orthogonal configuration [(ii)’ in Fig. 2e]. In this case, the spin current is reflected at the interface between Pt/top Fe2O3, leading to the relatively low resistance. The SMR valley in Fe2O3/Cr2O3/Fe2O3 occurs at the magnetic field that is much larger than that of net moment reversal (Supplementary Note 6), excluding possible coupling between net moment in Cr2O3 at high temperature and n in Fe2O3 as well as artifacts due to positive SMR from weak ferromagnetism (caused by defects or uncompensated interface). Then n in the bottom Fe2O3 rotates towards the spin-flop state (n ⊥ H) due to the increasing negative H, and n in the top Fe2O3 deviates towards n // H (I) [(iii)’ in Fig. 2e], resulting in the absorption of spin current and the second resistance peak (iii). It should be clarified that the second peak can exist when the coupling energy is large enough to overcome the Zeeman energy of the top Fe2O3 at the valley (ii), otherwise the n (top Fe2O3) will maintain spin-flop state rather than deviating towards n // H. The magnitude of the second peak is smaller than the first one can be ascribed to less component of n along x-axis. In contrast, the SMR in inverted sandwich, Fe2O3(4)/Cr2O3(4.4)/Fe2O3(12) (Supplementary Note 7), does not present resistance peak signal before H = 0, demonstrating that the n in the 12 nm-thick Fe2O3 maintains spin-flop state rather than deviating towards H at low magnetic field because of the large Zeeman energy, indicating the existence of the orthogonal coupling.
Besides magneto-transport measurements, we further confirm the interlayer coupling by direct Néel vector characterizations. Fe L-edge XMLD spectra were used to detect the n of the top Fe2O3 (several nanometers-thick sensitivity) in the Fe2O3/Cr2O3/Fe2O3 junction, where 2 nm-thick Pt was deposited on top. The XMLD spectra were recorded at zero-field after applying a high magnetic field along the x-axis for the sake of the non-volatile feature of n in easy-plane Fe2O316. X-ray was vertically incident to the film and the polarized direction was parallel to the film plane. XMLD signals are obtained as XMLD = XAS⊥− XAS//, where XAS// and XAS⊥ denote the x-ray absorption spectroscopy (XAS) recorded with polarization parallel with x-axis (//) and y-axis (⊥), respectively. Corresponding data are presented in Fig. 2f, where L2-edge is highlighted because it is generally used for the Fe-based XMLD spectra analysis32,33. Remarkably, Fe L2-edge XMLD spectrum exhibits a zero–positive–negative–zero feature, which is quite a characteristic for the n along the parallel direction (n // x-axis)32,33, rather than the spin-flop direction (y-axis). The n (top Fe2O3) aligned along H confirms the interlayer coupling, which is also corroborated by a series of XMLD measurements with sample rotation (Supplementary Note 8). The scenarios differ dramatically for the control samples Fe2O3 and Cr2O3/Fe2O3, in which identical experiments were carried out, but an opposite polarity at L2-edge (Fig. 2g,h, respectively), namely zero–negative–positive–zero, was observed, suggesting that the n in Fe2O3 is mainly aligned along the spin-flop direction (n // y-axis) without interlayer coupling.
Temperature dependence of interlayer coupling. We now turn towards the temperature dependence of SMR measurements in Fe2O3/Cr2O3/Fe2O3/Pt samples. Figure 3a shows the SMR results at various temperatures (Supplementary Note 9). At a relatively high temperature (T = 270 K), there exists two resistance peaks as we discussed above for T = 300 K, but the intensity of the second peak is lower than that at T = 300 K. Such a tendency continues with further decreasing temperature to 250 K, producing a tiny peak (or just a protruding), accompanied by the absence the second peak at 200 K. Also visible is that the location of the first resistance peak shifts towards zero-field with decreasing T but maintains at positive H, reflecting n // H in the top Fe2O3 before zero-field. This behavior discloses that although the interlayer coupling persists at low temperatures, the coupling energy decreases, resulting in the dominant spin-flop state and the disappearance of the second resistance peak. This phenomenon is similar to the temperature dependence of the spin fluctuation around equilibrium position in Cr2O3 spacer13. And we also demonstrate that the interlayer coupling does not depend on the magnetic field direction and the magnitude of applied reading current (Supplementary Note 10 and 11).
In addition to the H-dependent SMR we have explored the interlayer coupling between antiferromagnets by in-plane angle(α)-dependent SMR. Corresponding data of the Fe2O3/Cr2O3/Fe2O3/Pt sample at T = 300 K measured at two typical field of 0.5 T and 1 T are shown in Fig. 3b where α = 0º means H // I. For µ0H = 1 T, the SMR signals exhibit a negative polarity with the valley at α = 0º, which is a typical feature for the antiferromagnetic SMR at spin-flop state18,19. The situation differs dramatically for µ0H = 0.5 T. The SMR curve exhibits a positive polarity, indicating that the Néel vector of the top Fe2O3 maintains n // H due to the dominant interlayer coupling. This finding coincides with the results of field dependent SMR. In contrast, the polarity of SMR keeps negative in the control sample Cr2O3/Fe2O3/Pt (Supplementary Note 4), reflecting the antiferromagnetic feature of Fe2O3 and the absence of the interlayer coupling. Identical angle-dependent measurements were carried out in the Fe2O3/Cr2O3/Fe2O3/Pt sample with µ0H = 0.5 T at various temperatures. The polarity of SMR is positive at high temperatures (T = 350 and 300 K). With decreasing temperature to 250 K, the SMR signals become quite weak or even noisy, because of a competition between the interlayer coupling (n // H) and the H-induced spin-flop (n ⊥ H). This is accompanied by the typically negative SMR induced by the spin-flop with further decreasing temperature to 230 K and 200 K. The polarity of the control sample Cr2O3/Fe2O3/Pt is always negative at different temperatures (Supplementary Note 4), reflecting the absence of the interlayer coupling and further eliminating the existences of coupling between net moment in Cr2O3 and Néel vectors in top Fe2O3.
Interlayer coupling strength. It is significant to characterize the interlayer coupling strength. Considering that the existence of the first peak is the compromise between the interlayer coupling and spin-flop state, its location (µ0HCoupling) as a function of temperature for different Cr2O3 thicknesses (t) is summarized in Fig. 4a (Supplementary Note 12) to reflect the coupling strength. The first peak persists at a positive field for all of the measured SMR curves, suggesting the orthogonal antiferromagnetic interlayer coupling when the t ranges 3–4.4 nm. The maximum coupling strength increases with decreasing Cr2O3 thickness. For thin Cr2O3 (t = 3.0 and 3.5 nm), µ0HCoupling emerges from about 10 K, increases with increasing temperature and gets saturated at around 150 K. A plateau of the coupling strength exists from 150 K to 300 K. Then the coupling strength drops just above room temperature, which coincides with the spin fluctuation in Cr2O3 (bulk Néel temperature ~307 K)13. While for thick Cr2O3 (t = 4.1 and 4.4 nm), µ0HCoupling has an onset temperature of about 100 K, increases with increasing temperature and reaches the maximum just above room temperature, then drops, without showing a plateau. For thicker Cr2O3 (t = 6 and 12 nm), Fe2O3/Cr2O3/Fe2O3/Pt samples show almost the same SMR signals at 300 K as samples with only one Fe2O3 layer (Supplementary Note 13), indicating the absence of the interlayer coupling when the Cr2O3 layer is too thick. In addition, the interlayer coupling is observed in the Fe2O3/NiO/Fe2O3 junction with antiferromagnetic NiO spacer (Supplementary Note 14), indicating that the coupling effect is not restricted to a certain spaced material.
The distinct behavior for the samples with thin and thick Cr2O3 shows a significant role of spacer thickness. Since the out-of-plane anisotropy decrease rapidly with smaller sample size34, the most possible origin of the two distinct types of temperature dependence is that at low temperature, the moments in the thin Cr2O3 already have large in-plane component34–37, but the thick Cr2O3 have stable out-of-plane Néel order. The coupling arises from the fluctuating magnetic moments, hence thin Cr2O3 mediates coupling in low temperature, and thick Cr2O3 can only mediate coupling above 100 K.
As the coupling field shifts towards H = 0 with decreasing temperature, the magnitude of second resistance peak, which is also related to the deviation of Néel vector in top Fe2O3, also changes (Fig. 3a). We then summarize in Fig. 4b the temperature dependent magnetoresistance (MR) for the Fe2O3/Cr2O3/Fe2O3/Pt samples with various t. MR is related to the magnitude of the second resistance peak and is defined as MR = [R(second peak)–R(lowest)]/R(lowest), where the lowest is the minimum of the SMR curves. It can be seen that the MR exhibits a similar temperature dependence as the coupling field, strongly suggesting that the MR is also relevant to the interlayer coupling. As we discussed above (Fig. 2d), the second resistance peak is the result of competition between interlayer coupling and spin-flop state in the top Fe2O3 at the valley. When the temperature is low, the coupling energy is relatively small as compared with the Zeeman energy. Therefore, the Néel vector in the top Fe2O3 maintains spin-flop state (n ⊥ H), causing the vanishment of the second peak. The vanishing temperature (MR = 0) is obviously higher than its counterpart for no coupling effect (µ0HCoupling = 0), such as T = 225 K and T = 125 K for t = 4.4 nm, respectively. With further increase of temperature, the coupling energy is enhanced and exceeds the Zeeman energy, resulting in the deviation of more n towards H and the resultant rapid rise of the second resistance peak (MR). As the spin correlation is partially destroyed at high temperatures (T > 340 K for t = 4.4 nm), the coupling is reduced, accompanied by the decreasing of MR. An analogical situation occurs for the other samples (t = 3.0, 3.5, and 4.1 nm), while the samples with thin (t = 3.0 and 3.5 nm) and thick (t = 4.1 and 4.4 nm) Cr2O3 are divided into two groups according to the temperature dependence of the coupling strength.
Analysis on the magnetic ordering. In the following we discuss the role of magnetic ordering in the Cr2O3 spacer in the interlayer coupling of Fe2O3. The energies related to the Cr2O3 magnetic ordering are the interfacial coupling F(Nt,Fe, Mt,Fe, Nt,Cr, Mt,Cr), F(Nb,Fe, Mb,Fe, Nb,Cr, Mb,Cr), and the magnetic energy U in the Cr2O3 spacer, where N is the Néel vector and M is the net magnetization, t, b label the top and bottom surfaces, respectively. For thinner sample, the exchange energy makes it harder to let (Nt,Cr, Mt,Cr) ≠ (Nb,Cr, Mb,Cr). Thus, the observed interlayer coupling which increases with decreasing Cr2O3 thickness is hard to be explained by the non-uniform distribution of the magnetic order in Cr2O3 in the thickness direction. Also, it is known that Cr2O3 is lack of inter-unit cell DMI31, which favor out-of-plane spiral spin structure and may cause (Nt,Cr, Mt,Cr) ≠ (Nb,Cr, Mb,Cr). If (Nt,Cr, Mt,Cr) = (Nb,Cr, Mb,Cr), assuming F(Nt,Fe, Mt,Fe, Nt,Cr, Mt,Cr) ≤ F(Nb,Fe, Mb,Fe, Nb,Cr, Mb,Cr), one can lower the total energy by rotating (Nb,Fe, Mb,Fe) towards (Nt,Fe, Mt,Fe). The process above is solid even when the interfacial coupling at different interfaces has a different magnitude, as long as the interfacial coupling has the same form. Thus, by considering the Cr2O3 magnetic ordering which is uniform in the film plane, the lowest energy state always has (Nb,Fe, Mb,Fe) = (Nt,Fe, Mt,Fe), i.e., no orthogonal interlayer coupling can be generated.
Quasi-long range order mediated interlayer coupling. Having excluded the magnetic ordering which is uniform in the film plane as the origin of the interlayer coupling, we thus consider magnetic ordering which is non-uniform in the film plane, i.e., quasi-long range order (QLRO) enabled by fluctuations. It is known that in extremely thin magnets with in-plane moments, long-range magnetic order could be absent easily38, while spin correlation still exists39, which distinguishes from disorder paramagnetic phase [Fig. 4(e)]. Indeed, such QLRO emerges when the film gets thinner and the system drops from three to two dimension, consistent with the larger coupling strength in thinner films. Meanwhile, the temperature-dependent behavior of µ0HCoupling in Fig. 4a also strongly suggests that the interlayer coupling is correlated to the spin fluctuation in the Cr2O3 spacer. For QLRO state in the Cr2O3 sandwiched by two Fe2O3 layers, the parallel and perpendicular arrangements of the n in the two Fe2O3 layers means different boundary conditions. The perpendicular arrangement provides potential with \({D}_{2d}\) symmetry for QLRO rather than the \({D}_{2h}\) symmetry in parallel arrangement condition. Since the in-plane component of the moments point to all in-plane directions in QLRO state, the perpendicular arrangement provides more symmetric environment for such excitations than the parallel arrangement. Hence, we assume a smaller energy in the perpendicular arrangement, with the same entropy in the two arrangements since the entropy only depends on the QLRO itself40. This smaller energy makes the perpendicular arrangement favorable in the free energy, leading to the orthogonal interlayer coupling in the QLRO state.
Based on the above considerations, the schematic free energy diagrams at zero field for the thin and thick Cr2O3 cases are shown in Fig. 4(c) and (d), respectively. The ordered ground state has a free energy of F0, and all other states are represented by their free energy difference with the ground state F – F0. In thin Cr2O3 cases, no coupling exists in the ground state. At T1, QLRO emerges under the perpendicular condition, while the Cr2O3 under the parallel condition remains in the ground state. Therefore, the coupling emerges with a coupling strength represented by F0 – F([email protected]°) and increasing from T1 to T2. At T2, the Cr2O3 under the parallel condition switches from the ground state to the QLRO state, and the coupling strength saturates to F([email protected]°) – F([email protected]°). Then the Cr2O3 under the parallel condition switches from the QLRO state to paramagnetic state at T3, and the coupling strength F(paramagnet) – F([email protected]°) starts to decrease. The coupling strength finally vanishes at T4, where the Cr2O3 under the perpendicular condition switches to the disorder paramagnetic phase. Similar process happens to the thick Cr2O3 cases, with a larger T1 because the perpendicular anisotropy stabilizes the Néel order in the ground state. The coupling strength increases up to T3*, where the Cr2O3 under the parallel condition switches from the ground state directly to the disorder state, without entering the QLRO state. Then the coupling strength decreases and finally vanishes at T4.
The phase transition from the ordered ground state to the QLRO state can be phenomenally described by the following free energy

where T is the temperature, b is a parameter to stabilize the Néel Order, and λ is a parameter related to the entropy difference between the ground state and the QLRO. Δ reflects the energy difference between the ground state and the QLRO, with values Δ⊥ and Δ// for the perpendicular and the parallel conditions, respectively. Δ⊥ is smaller than Δ// and Δ// – Δ⊥ increases with decreasing Cr2O3 thickness. The Néel Order n and the volume of the QLRO state q satisfies

The first term in (1) describes the blue and red lines in Fig.4(c) and (d), and a combination of (1) and (2) gives the n-related part of the free energy ~ (Δ/λ–T) n2 + bn4, yielding a Néel temperature Δ/λ and n ~
near the Néel temperature40.The coupling energy Ec, the free energy difference under the two conditions (Supplementary Section 15),is displayed in Fig.4(f), obtaining a correspondence of 0.01 meV(per unit cell) to the coupling field 0.86 T, which semi-quantitatively agrees with the experimental curves (≤ 200 K).