Abstract
The recent observation of the anomalous Hall effect (AHE) without notable magnetization in antiferromagnets has suggested that ferromagnetic ordering is not a necessary condition. Thus, recent theoretical studies have proposed that higher-rank magnetic multipoles formed by clusters of spins (cluster multipoles) can generate the AHE without magnetization. Despite such an intriguing proposal, controlling the unconventional AHE by inducing these cluster multipoles has not been investigated. Here, we demonstrate that strain can manipulate the hidden Berry curvature effect by inducing the higher-rank cluster multipoles in spin-orbit–coupled antiferromagnets. Observing the large AHE on fully strained antiferromagnetic Nd2Ir2O7 thin films, we prove that strain-induced cluster T1-octupoles are the only source of observed AHE. Our results provide a previously unidentified pathway for generating the unconventional AHE via strain-induced magnetic structures and establish a platform for exploring undiscovered topological phenomena via strain in correlated materials.
INTRODUCTION
The anomalous Hall effect (AHE) is a fundamental transport phenomenon that has been universally observed in time-reversal symmetry broken systems. AHE can arise from two different forms of mechanism (1): extrinsic mechanism, such as skew scattering or side jump due to magnetic impurities, and intrinsic mechanism originating from Berry curvature in momentum space. Since the fundamental topological properties of electronic wave functions are encoded in the Berry curvature, AHE is considered as a powerful tool for probing the topological properties of materials (2, 3). In addition to its fundamental interest, AHE can be applied for memory devices (4).
Conventionally, AHE has been observed mostly in itinerant ferromagnets. Its magnitude is known to be proportional to the magnetization (5), which is a measure of broken time-reversal symmetry. Recently, a large AHE has been unexpectedly found in noncollinear antiferromagnets, such as Mn3X (X = Sn, Ge) (6–8) and GdPtBi (9), which do not exhibit spontaneous magnetization. This unconventional response indicates that ferromagnetism is not a necessary condition for AHE and suggests a possible alternative origin of AHE. A recent theory proposed that higher-rank magnetic multipole (cluster multipole) moments formed from spin clusters in antiferromagnet can induce a nonzero AHE, beyond the conventional dipoles of ferromagnets (10). Subsequently, the anomalous Nernst (11) and magneto-optical Kerr effects (12) in Mn3Sn have also been attributed to its cluster octupoles. However, since antiferromagnets are not easily coupled to both magnetic and electric fields (13), it is very difficult to manipulate the higher-rank cluster multipoles. This imposes substantial limitations on controlled experiments on the cluster multipoles and associated AHE.
Here, we demonstrate that the strain can generate the AHE by inducing the higher-rank cluster multipoles, by using antiferromagnetic Nd2Ir2O7 (NIO) thin film. Further investigation reveals that biaxial strain on the pyrochlore lattice can modulate the spin structure and induce certain magnetic octupoles. The induced cluster octupoles can generate the net Berry curvature effect hidden in the bulk, leading to a finite AHE. We expect that our method could be widely applied to other spin-orbit–coupled topological magnets (10) and antiferromagnetic spintronics (4, 14).
RESULTS
Strain-induced cluster multipoles in a pyrochlore lattice
The NIO belongs to the pyrochlore iridates family, R2Ir2O7 (R, rare-earth ions). The members of the family are geometrically frustrated magnets with complex lattice structures. As shown in Fig. 1A, R2Ir2O7 is composed of linked tetrahedrons with R and Ir at each vertex. In R2Ir2O7, strong electron correlations and large spin-orbit coupling of Ir d electrons result in unique antiferromagnetic spin structures, called all-in-all-out (AIAO) ordering (15, 16). As shown in the circle in Fig. 1B, the spins in one tetrahedron point inward and those in the neighboring tetrahedron point outward. The Néel temperatures of the Ir and Nd sublattices of bulk NIO are
(A) Pyrochlore lattice structure. Yellow (red) circles depict Nd (Ir) ions. Note that oxygen ions are not shown. (B) Schematic diagram of epitaxial NIO thin film on the yttria-stabilized zirconia (YSZ) substrate with biaxial strain along [111]. The deformed pyrochlore lattice is schematically displayed in the circle. The blue arrow at each site denotes the spin direction in the AIAO antiferromagnetic configuration. (C) The spin arrangement of a tetrahedron in the undistorted (i.e., bulk) AIAO configuration. According to the cluster magnetic multipole theory, the AIAO magnetic ordering can be represented as an A2-octupole. (D) Schematic diagram of the strained magnetic ground state. Biaxial strain along the [111] direction will distort the AIAO configuration, which can be represented as the superposition of a cluster dipole (M), an A2-octupole, and a T1-octupole (ω). M represents a ferromagnetic ordering, while ω represents an antiferromagnetic ordering other than AIAO. On the basis of symmetry analyses (see section S1 and table S1), we demonstrated that only the T1-octupole can induce the AHE without magnetization.
However, since AIAO ordering preserves the cubic crystalline symmetry, the net Berry curvature effect vanishes when we integrate over the Brillouin zone (BZ). Unless the cubic crystalline symmetry is broken, AHE cannot be observed in this system. To break the cubic symmetry, a magnetic field was applied to pressured NIO single crystals (19) and Pr-doped bulk samples (20, 21). However, the spin structures modulated by the magnetic field are fragile and easily return once the magnetic field is turned off. Thus, a stable method to break the cubic symmetry is highly desirable; here, we choose a strain engineering approach and investigate the associated AHE.
As shown in Fig. 1B, the biaxial strain elongates the unit tetrahedra along the [111] direction. This will naturally break the cubic symmetry of the system. Since the deformation modulates magnetic anisotropy (22), the Ir spin directions should be changed. To systematically describe the change of spin direction, we adopted the cluster multipole theory. Since the conduction electrons come from Ir d orbitals, we considered Ir sublattice only (16). In the cubic pyrochlore lattice, all spin ordering patterns can be classified into five different irreducible representations, carrying 12 distinct cluster multipoles (18). Among them, certain cluster multipoles that break the cubic symmetry are responsible for the AHE (see section S1).
In a bulk NIO, the AIAO ordering corresponds to a higher-rank magnetic multipole called the A2-octupole (Fig. 1C). Since the A2-octupole preserves the cubic symmetry, it cannot generate AHE. However, in a strained NIO (s-NIO), the AIAO spin structure becomes modulated under the strain. The resulting spin configuration is denoted by strained AIAO (s-AIAO), composed of a superposition of three kinds of cluster multipoles, namely, a dipole, an A2-octupole, and a T1-octupole (Fig. 1D). Note that the dipole is just the ferromagnetic ordering, while the T1-octupole is an antiferromagnetic ordering other than AIAO. Only the T1-octupole can induce the AHE without magnetization since it breaks the cubic symmetry as the dipole does.
Characterizations of relaxed and fully s-NIO thin films
To investigate the strain-induced magnetic multipole and associated AHE, we prepared two kinds of NIO thin films on the yttria-stabilized zirconia (YSZ) substrates: relaxed and fully strained films. The biaxial strain can arise from the lattice mismatch between the R2Ir2O7 film and the YSZ substrate (see Fig. 1B) (23, 24). Since the lattice parameter of YSZ is smaller than those of NIO bulk, the NIO film should be compressively strained. We estimated the strain
Despite the substantial past efforts (25–27), the in situ growth of high-quality R2Ir2O7 thin film is notoriously difficult. Under the proper crystalline growth conditions for pyrochlore oxides (28), the corresponding R2Ir2O7 phase becomes extremely unstable because of the formation of a gaseous IrO3 phase (29). To avoid this instability, many studies have used the “solid-phase epitaxy (SPE)” (25, 27) method, which involves the initial growth of amorphous R2Ir2O7 films at a lower temperature (T) followed by ex situ thermal annealing in a sealed tube. Although SPE can provide a method for the growth of single-phase R2Ir2O7 films, it usually produces relaxed films (25, 26). Therefore, we developed a previously unknown in situ film growth method, i.e., repeated rapid high-temperature synthesis epitaxy (RRHSE; see section S2 and Materials and Methods) (30).
The RRHSE method made us successfully grow the fully s-NIO films on YSZ (111) substrates. Figure 2A shows an x-ray diffraction θ-2θ scan. The NIO (lll) and YSZ (lll) (l: integer) peaks can be seen, indicating epitaxial growth of NIO single phase. Particularly, the satellite peaks near the NIO (222) peaks are observed, which is commonly called “thickness fringes.” These interference peaks indicate the high quality of a sharp interface between film and substrate. Figure 2B shows x-ray reciprocal space mapping around the NIO (662) and YSZ (331) Bragg peaks of a 9-nm-thick NIO film. The lattice parameter of the (662)-plane, d(662), of bulk NIO is 1.19 Å, and the d(331) of YSZ is 1.18 Å. Note that the NIO (662) Bragg peak has the same Qx value as the YSZ (331) peak, demonstrating that our film becomes fully strained (~1% compressive strain) by the substrate.
(A) X-ray diffraction θ-2θ scans of an epitaxial NIO film grown on a YSZ (111) substrate. The scans reveal that the NIO film was grown coherently with the YSZ substrate. (B) Reciprocal space mapping around the YSZ (331) and NIO (662) diffraction peaks. The Qx values of both peaks are the same, indicating that the film is under ~1% strain. (C) Scanning transmission electron microscopy image with the zone axis parallel to [1-10]; a clear interface between the film and the substrate can be seen. The distinctly colored circles indicate the pyrochlore structural ordering of the Nd and Ir ions. Images of selected areas in (C) were fast Fourier–transformed for (D) the NIO thin film, and (E) the YSZ substrate. The three vertical dotted lines between (D) and (E) are plotted without changing the scale. These lines indicate that the in-plane lattice constants of NIO and YSZ are the same, providing further direct evidence for the fully s-NIO thin film.
Figure 2C shows a scanning transmission electron microscopy image that indicates the high quality of our film. The NIO pyrochlore phase is formed with few structural defects or disordered structures. Figure 2 (D and E) shows fast Fourier transform patterns from the selected areas in the film and substrate, respectively, marked in Fig. 2C. As demonstrated by the red dotted lines, the as-grown NIO film has the same inverse lattice constant as the YSZ substrate, which also confirms that our film is fully strained.
Electronic structures of relaxed and fully s-NIO thin films
We compared these fully s-NIO films grown by RRHSE with the relaxed NIO (r-NIO) films grown by the SPE (see section S3). The resistivity ρ (T) curve of a 9-nm-thick s-NIO film exhibits a semimetallic behavior at most T. As shown in Fig. 3A, the s-NIO film has ρ (T) an order of magnitude smaller than that of the r-NIO film. The ρ (T) curve of an 80-nm-thick r-NIO film exhibits a metal-insulator transition around ~30 K (black dashed line in Fig. 3A), in agreement with its bulk counterpart (17, 31). The strong upturn of the r-NIO film is due to its insulating nature below
(A) Orange solid (black dashed) line indicates the resistivity, ρxx, of the fully s(r)-NIO film, prepared by the RRHSE (SPE) method. ρxx reveals that the electronic structure of NIO could be changed under strain. The calculated band structure of the (B) bulk and (C) 1% biaxial s-NIO is shown. The insulator-to-semimetal transition can explain the large change in ρxx in (A). (D) Measured AHC
To understand the corresponding electronic structure changes, we performed mean-field calculations using the Hubbard model (see section S4). The previous study shows that the most valence and conduction bands near the Fermi energy come mainly from Ir 5d electrons (16). Our calculated electronic structure of the bulk (the r-NIO film in our case) explains its insulating nature. The energy gap opens with a value of about 13 meV (Fig. 3B), which agrees well with the bulk value (32). Under 1% compressive strain, the valence and conduction bands move, which slightly increases the direct gap at most k regions. However, some valence and conduction bands become crossed with Fermi level; thus, small electron and hole pockets develop near L1,2,3,4 (Fig. 3C), creating a semimetallic state. These model calculations can explain why the s-NIO film has a much smaller ρ (T) than the r-NIO film.
Large AHE in fully s-NIO thin films
Besides, the s-NIO film shows a much larger anomalous Hall conductivity (AHC) compared to the r-NIO film. Figure 3D shows the magnetic field (H)–dependent AHC
To cross-check, we compared our magnetotransport property values with those of ferromagnets. For example, (Ga, Mn)As (34) and CuCr2Se4–xBrx (35) typically exhibit SHC with
Accordingly, we calculated the Berry curvature effect on AHC from the band structure obtained from the mean-field calculations mentioned above (see section S4). AHC can be obtained by integrating the Berry curvature
Antihysteresis of AHC
Another notable feature of s-NIO film is that its
To understand our antihysteresis-like
(A) AHC
Figure 4C summarizes the results of the AHC fitting at H = −9 T with the
AHE from strain-induced T1-octupoles
To reveal the relation of AHE and cluster multipoles under strain, we should compare M and
(A) The M/H curve (purple squares) of s-NIO film at 3 K, overlaid with the experimental AHC (orange circles). Note that M = 0 but
To elucidate how T1-octupole emerges under the strain, we calculated the spin structure from the spin model. Since both Nd and Ir spins play important roles, we included the Heisenberg, Dzyaloshinskii-Moriya, anisotropic spin-exchange interactions between Ir spins (36), the f-d exchange interaction between the Ir and Nd spins (17), and the Zeeman energy for both the Ir and Nd spins (for details, see section S7). On the basis of the calculated spin structure, we obtained the cluster multipoles (table S1 in section S1). Figure 5C shows the calculated dipole (M, green circles) and T1-octupole (ω, blue circles) as a function of the effective Zeeman energy h in the r-NIO. According to our calculation, r-NIO does not have a finite M or ω value for the Ir sublattice at h = 0. The zero values of M and ω can explain the negligible SHC of the r-NIO film (see Fig. 3D). Figure 5D shows the calculated M and ω of s-NIO, which are finite even for h = 0. Particularly, the hysteresis curve of ω looks similar to the
DISCUSSION
Our work demonstrates that the strain-engineering of an antiferromagnet can generate the net Berry curvature effects by modulating its cluster multipoles. In particular, our findings highlight that the strain-induced T1-octupole is closely connected with the topological properties of NIO. We can further extend this strain-engineering approach to search for the other novel topological phenomena in correlated magnets. For example, the strain engineering approach on numerous series of R2Ir2O7 can realize novel correlated topological phases, such as Weyl semimetal, axion insulator (16), a strong topological insulator (18), and line-node semimetal (19, 21) by properly modifying their magnetic structure. In this perspective, we believe that our strain study on NIO could provide a cornerstone to discover strain-engineered novel topological phenomena in oxides and to understand their fundamental mechanisms.
MATERIALS AND METHODS
Film growth and structural characterization
Fully s-NIO films were in situ grown on insulating YSZ substrates using the RRHSE method. This film growth method is a modified form of pulsed laser deposition, based on repeated short-term thermal annealing processes using an infrared laser. RRHSE consists of two key steps in one thermal cycle. During the first step, amorphous stoichiometric NIO and IrO2 layers were deposited by a KrF excimer laser (λ = 248 nm, 5 Hz) at T ~ 600°C with PO2 ~ 50 mtorr. The additional IrO2 layer was deposited to compensate for the Ir loss that would unavoidably occur later during the synthesis process. During the second step, the pyrochlore phase is formed by rapidly raising T to 800°C (up to ~400°C min−1). We must expose the sample to the high T for a period that is sufficiently long to synthesize the pyrochlore phase but short enough to minimize the formation of IrO3. Last, we repeated these deposition and thermal synthesis processes until the desired film thickness was obtained. During the growth, the reflection high-energy electron diffraction pattern was monitored and the intensity of the oscillation was recorded. After growth, NIO films were characterized by an x-ray diffractometer (Bruker Corp.) and an atomic-resolution high-angle annular dark-field scanning transmission electron microscope (JEM-ARM200F; JEOL) equipped with an energy-dispersive x-ray spectrometer.
Transport and magnetic properties
Magnetotransport properties were measured via a standard four-point probe method using a commercial physical property measurement system (PPMS, Quantum Design), which has a base T of 2 K and a maximum magnetic field of 9 T. During the measurements, the current was applied along the [1-10] direction, and H was applied along the [111] direction. Magnetization data were obtained using a commercial superconducting quantum interference device magnetometer (MPMS, Quantum Design) with the magnetic field applied normal to the film.
The AHC value
Self-consistent mean-field Hubbard model
We developed the Hubbard model for the s-NIO thin film under the magnetic field and acquired the ground state and electronic structure by the self-consistent mean-field method. We adopted 24 × 24 × 24 and 32 × 32 × 32 k-point mesh and found that the results are consistent. We calculated the AHC by integrating the Berry curvature, adopting a 48 × 48 × 48 k-point mesh of the entire BZ. Details of the calculation are provided in section S4.
Spin model calculation
We developed the spin model including Heisenberg exchange, Dzyaloshinskii-Moriya interaction, anisotropic exchange, the f-d exchange between Nd and Ir electrons, and the Zeeman effect by applying second-order perturbation theory to the Hubbard model and referring to previous works. We calculated the ground state by the iterative minimization method, which repeatedly aligns spins to the effective field direction until each spin is fixed. Details of the calculation are provided in section S7.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/29/eabb1539/DC1
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 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).