Introduction

Large and bistable magnetoresistance signals have been observed in tunnelling devices with an antiferromagnetic (AFM) IrMn layer on one side and a non-magnetic metal on the other side of the tunnel barrier1,2. The work has experimentally demonstrated the feasibility of a spintronic concept3,4,5 in which the electronic device characteristics are governed by the staggered magnetization axis in an AFM. It has been shown that the AFM moments can be manipulated via an exchange-coupled ferromagnet (FM)1,2,6 and that the AFM magnetoresistance signals can persist to room temperature6. No stray fields and the relative insensitivity to external magnetic fields are among the features that make AFMs attractive complements to the conventionally utilized FMs in the design of spintronic devices, as highlighted in a recent study of AFM linear chains of a few Fe atoms7 and in the demonstration of the concept of an AFM memory8. The possibility to design spintronic elements based on AFMs becomes even more attractive in the context of magnetic counterparts of conventional compound semiconductors, which may enable new devices combining spintronic and nanoelectronic functionalities. In this paper, we report the discovery of a new member of this family: tetragonal epitaxial CuMnAs, a room-temperature AFM whose in-plane lattice constant matches GaP or Si.

Over the past two decades, the introduction of magnetism into common semiconductor hosts has driven a number of new research areas in spintronics. Inspired by spintronics research and applications based on transition metal ferromagnets (FMs), the focus has been on magnetic counterparts of semiconductors with the FM order. (Ga,Mn)As and related (III,Mn)V compounds9,10,11 have become archetypes among these materials, resulting in discoveries of new spin-related physical phenomena and device functionalities. For example, the understanding of spin-orbit-coupling phenomena has advanced because of experiments in (Ga,Mn)As from the ohmic transport regime to new realizations in tunnelling and Coulomb blockade devices12,13. Low Curie temperatures have prevented the direct integration of (III,Mn)V FM semiconductors into spintronic technologies. Nevertheless, new phenomena discovered in (Ga,Mn)As have been subsequently reported in the room-temperature metal FMs, relevant not only for basic science but also for practical spintronic applications14,15,16.

Although a ferromagnetic ground state is rare within magnetic materials derived from semiconductor compounds and their Curie temperatures are below room temperature, high Néel temperature (TN) AFMs can be found, for example, among the magnetic counterparts of the I–II–V semiconductors (see also Supplementary Table S1 and Supplementary Note 1)4,5. The structures of the non-magnetic semiconductors such as LiZnAs, NaZnAs, CuZnAs and AgZnAs are closely related to the III–V zincblende structure, as shown in Fig. 1a,b. By splitting the group III atom into two elements from groups I and II, their combined valence equals that of the group III atom. One of the two elements resides on the zincblende site, whereas the other occupies one of the tetrahedral interstitial sites that is empty in the zincblende structure17. An extensively studied magnetic compound with the same filled zincblende (also called half-Heusler) crystal structure is CuMnSb (Fig. 1c)18,19. The AFM coupling of FM (111) planes of this cubic crystal has a frustration in the magnetic interactions, with half of the Mn nearest-neigbours coupled antiferromagnetically and half ferromagnetically20. This frustration leads to the relatively low TN=50 K of CuMnSb and shows that keeping the same cubic structure as in the non-magnetic parent I–II–V compound is unfavourable for robust antiferromagnetism. Moreover, lowering the symmetry from cubic to, for example, tetragonal enhances the magnetocrystalline anisotropy phenomena and is therefore favourable for the concept of AFM spintronics.

Figure 1: Unit cell structures of non-magnetic semiconductors and their magnetic counterparts.
figure 1

(a) III–V zincblende structure. (b) I–II–V half-Heusler structure. (c) Half-Heusler CuMnSb. (d) Tetragonal LiMnAs. (e) Orthorhombic CuMnAs. (f) Tetragonal CuMnAs. The bonds in (c), highlighted in green and red, show the AFM and FM nearest-neighbour-exchange coupling of Mn present in the half-Heusler CuMnSb. The bonds in (f), highlighted in green, show that all Mn nearest neighbours are coupled antiferromagnetically in the tetragonal CuMnAs, which is favorable for high TN.

LiMnAs and NaMnAs are among the examples of AFM I–II–V semiconductors4,21,22,23,24, whose equilibrium crystal structure changes from the cubic half-Heusler lattice of their non-magnetic counterparts to a layered tetragonal structure (Fig. 1d). This more-anisotropic-crystal arrangement removes the frustration in the nearest-neighbour magnetic coupling, resulting in TN’s far above the room temperature. Single-crystal thin films of the LiMnAs AFM semiconductor were recently prepared by molecular beam epitaxy (MBE) on lattice-matched InAs substrate4,24,25. The inclusion of alkali metal elements represents, however, a challenge both in terms of the growth and the stability of devices. Group Ib transition metal elements may represent the solution to this problem: NaZnAs and AgZnAs are known to be twin compounds with identical crystal structure and lattice constants, and CuZnAs is also very similar17. This has motivated our interest in CuMnAs as a suitable I–Mn–V AFM compound for spintronics20. The bulk equilibrium phase of CuMnAs displays room-temperature AFM ordering26, as confirmed in our recent study of chemically synthesized bulk samples20. However, the orthorhombic crystal structure (Fig. 1e) is not compatible with conventional compound semiconductor substrates.

In the present work, we demonstrate that CuMnAs films can be prepared by MBE in the tetragonal form reminiscent of NaMnAs (Fig. 1f). In analogy to the Ia–Mn–V compounds, all nearest neighbours in the tetragonal CuMnAs films couple antiferromagnetically, resulting in a high TN. We demonstrate relaxed epilayers of the tetragonal phase of CuMnAs grown on GaAs substrate and fully strained CuMnAs films grown on the closely lattice-matched GaP. The functionality of this new AFM I–II–V compound is demonstrated in all-epitaxial CuMnAs/Fe AFM/FM bilayers, in which we observe large room-temperature-exchange bias effects.

Results

Synthesis and structure of CuMnAs epilayers

To stabilize the tetragonal phase of CuMnAs, layers were grown on the zincblende III–V semiconductor substrates GaAs(001) and GaP(001). In situ reflection high-energy electron diffraction (RHEED) and ex situ X-ray diffraction (XRD) measurements, shown in Fig. 2, revealed a structure incompatible with the chemically synthesized bulk material studied by Mundelein and Schuster26. The space group of this structure was determined using a combined experimental and computational method that enables the precise determination of the atomic positions in a thin film using conventional Cu Kα X-ray sources only27. The space group was found to be that of the tetragonal Cu2Sb family, P4/nmm. The lattice parameters of epitaxial CuMnAs were found to be a=b=3.820 Å and c=6.318 Å for the relaxed material grown on GaAs.

Figure 2: Growth and structural characterization of CuMnAs.
figure 2

(a,b) RHEED images of the surface of a CuMnAs layer during growth on a GaP substrate taken along the substrate [110] and [10], respectively. (c) Diagram depicting the CuMnAs unit cell structure and spin arrangement, with Mn (purple), Cu(blue) and As(yellow). (d) A birds-eye view of how the layer grows on GaAs and GaP substrates, the green balls represent the As or P atoms. (e) X-ray 2Θ/ω scans of the CuMnAs (003) diffraction peak for layers grown on GaP and GaAs. The experimental data are shown as grey dots, whereas the red line is a mosaic block model fitting for the layer grown on GaAs. This model is not applicable for the layers on GaP. The corresponding (003) reciprocal space maps are shown in the insets. (f,g) Z-contrast STEM image of CuMnAs along the [100] direction. The As positions are clearly visible as the brightest atoms, and intensity analysis along the rows reveal a highly ordered compound in complete agreement with the overlaid model structure in (g). The CuMnAs unit cell, as depicted in panel (c), is highlighted by a white dashed box in (g).

Figure 2c shows the unit cell of CuMnAs and Fig. 2d a birds-eye view of its registry with the substrate. The closest match to the two substrates used, GaAs and GaP, is with the CuMnAs principal in-plane axis matching along the half diagonal of the GaAs/GaP unit cell, and this is indeed the orientation it robustly adopts for a range of temperatures and stoichiometries. In the case of the GaAs substrate, the match is only preserved for the first few monolayers whereupon the layer relaxes, as evidenced by a dimming in the RHEED pattern. However, it continues to grow with the same relative orientation, and the RHEED pattern reappears several nanometres later. For the GaP substrate, a mismatch of ~0.3% is observed along the same geometry and the layer remains strained up to the thickest layers grown so far (130 nm).

Figure 2e shows X-ray 2Θ/ω linescans and reciprocal space maps around the CuMnAs(003) diffraction peak for CuMnAs grown on GaAs(001) and GaP(001). The broadening of the reciprocal space map, as evidenced in Fig. 2e (left inset) of the CuMnAs(003), is attributed to mosaicity introduced by the poor lattice matching of CuMnAs to the GaAs substrate. The high strain and subsequent relaxation introduce defects, leading to mosaic blocks. This is supported by the fitting of a mosaic block model to the CuMnAs/GaAs that shows good agreement with a model assuming ellipsoid blocks with lateral dimensions of 50 nm, vertical dimensions of 20 nm and RMS misalignment of 1.5° (red line in Fig. 2e)28,29. The much sharper spot produced by the CuMnAs on GaP (right inset) indicates a more coherent layer that no longer fits a mosaic block model well.

The image shown in Fig. 2f is a Z-contrast scanning transmission electron microscopy (STEM) slice along the CuMnAs [100] direction. In the enlarged region shown in Fig. 2g, the brighter As atoms are clearly visible in a triangular arrangement that matches exactly to the X-ray-determined structural model along the same [100] direction, overlaid on the top right side of the image. The Cu and Mn, which are much closer in size, are not clearly distinguishable by eye; however, when the intensity is profiled along the atomic planes, the Cu and Mn can be distinguished and they show a high degree of chemical order.

Magnetic characterization of CuMnAs

Magnetic properties of the CuMnAs layers were examined using a combination of neutron diffraction and superconducting quantum interference device magnetometry. Neutron diffraction measurements were taken using the WISH instrument at the ISIS neutron source. Two scattering geometries were used to access the (h00) and (00l) sets of reflections and revealed the presence of the structurally forbidden (100) magnetic reflection. The absence of this peak structurally was confirmed using in-plane XRD, which as expected reveals the presence of the (200) and (400) reflections but not the forbidden (100) reflection. Fig. 3a shows the (100) magnetic peak, which confirms the presence of long-range AFM order with the same dimensions as the structural unit cell. This peak is completely absent at 500 K, whereas the nuclear (200) peak (not shown) and the structural contribution of the (001) peak are still present (Fig. 3b), confirming that the layer has been neither destroyed nor undergone a phase change, but has lost long-range AFM order.

Figure 3: Magnetic characterization.
figure 3

(a) Neutron-diffraction data showing the structurally forbidden CuMnAs (100) diffraction peak at 300 K (filled symbols) and 500 K (open symbols). (b) Neutron-diffraction from the CuMnAs (001) peak measured at the same temperatures. At 300 K, this reflection combines both nuclear and magnetic contributions. (c) Hysteresis loops performed on a 3-nm Fe film on CuMnAs at 300 K after cooling in a 1-T field from 430 K. A stable exchange bias is observed symmetrically for positive (blue) and negative (orange) field cool directions.

An analysis of the possible spin configurations, classified according to the irreducible representations associated with the -point (k=0) of the P4/nmm space group, reveals that only models involving basis functions of the bidimensional representation are consistent with the neutron-diffraction results. These models imply AFM coupling of the two symmetry-related Mn sites, with spins confined to the (ab) plane (Fig. 2c). Precise determination of the direction of the spin axes within this plane is complicated by the expected magnetic domains. The room-temperature data normalized to a vanadium spectrum and scaled to the nuclear (200) peak allowed us to estimate the magnetic moments of Mn to be 3.6(2) μB, assuming that the possible magnetic domains are equally populated.

Superconducting quantum interference device magnetometry performed on stoichiometric CuMnAs layers revealed no ferro- or paramagnetic response, consistent with a pure AFM. In order to provide further support for the magnetic properties of the samples, CuMnAs layers were grown with 3 nm of epitaxial Fe on top and capped to prevent oxidation. By field-cooling these layers through 430 K, it is possible to induce an exchange bias shift at and above room temperature. This exchange bias is evident in the hysteresis loops shown in Fig. 3c, and is robust and reproducible across many samples (Supplementary Figs S1 and S2, and Supplementary Note 2). The effect remains present after training, and there is also a large change in the coercive field of the Fe layer compared with Fe on GaAs alone30. Heating above 440 K was found to damage the Fe/CuMnAs interface in a similar manner to that observed for Fe/GaAs heterostructures30,31,32, and after this permanent damage it is no longer possible to induce exchange bias (Supplementary Fig. S2a). The permanent changes to the interface were confirmed using temperature-dependent X-ray reflectivity and XRD, where roughening and density changes are observed in the interfacial region but the bulk crystal remains intact.

Transport measurements on CuMnAs

Transport measurements on stoichiometric CuMnAs samples are shown in Fig. 4. The resistivity at 4 K is about 90 μΩ cm rising to about 160 μΩ cm at 300 K. These values are compatible with the semimetallic-like band structure suggested by the GGA+U calculations (detailed in the next section), with strongly suppressed density of states around the Fermi energy. The measured positive Hall coefficient of 6 × 10−8 Ω cm T−1 interpreted in a single-carrier model would give a carrier density of 1.1 × 1022 cm−3. However, as discussed in Maca et al.20 for the orthorhombic bulk CuMnAs and confirmed by our calculations for the tetragonal CuMnAs presented below, the compound’s electronic structure is at the transition from a semiconductor to a semimetal. The nature of electronic states near the Fermi energy may therefore be sensitive to relatively small changes of intrinsic and extrinsic parameters of the CuMnAs films. In this situation, where the band structure is complex, the single-carrier interpretation is inappropriate and can greatly overestimate the number of carriers if (as the calculations indicate) the numbers of hole-like and electron-like carriers are comparable.

Figure 4: Transport in CuMnAs.
figure 4

(a) shows the temperature dependence of resistivity and (b) the Hall resistivity as a function of applied magnetic field.

Calculations

From GGA+U band-structure calculations, shown in Fig. 5, we confirm that the AFM phase in CuMnAs has lower energy than the FM phase, and the difference is 109 meV f.u.−1. For comparison, calculations for the cubic half-Heusler CuMnSb give the difference of only 4 meV f.u.−1. This suggests that the Néel temperature of the tetragonal CuMnAs is much higher than the TN=50 K of CuMnSb (ref. 33), as indeed is observed in our experiments. The robust AFM order is retained when the CuMnAs crystal transforms from the equilibrium bulk orthorhombic structure26 to the tetragonal structure stabilized in our thin film epilayers.

Figure 5: Density of states calculations.
figure 5

GGA+U-calculated density of states of tetragonal CuMnAs for (a) antiferromagnetic and (b) ferromagnetic configurations. Upper and lower sections are for up and down spins, respectively. Blue and red lines in (a) are the local density of states for Mn in each sublattice. Blue lines in (b) are the local density of states of Mn. The black lines in both panels are the total density of states.

The GGA+U calculations show that tetragonal CuMnAs in the AFM ground state has a strongly suppressed density of states around the Fermi energy (see Fig. 5a), indicating an electronic structure of a semimetal or a narrow bandgap semiconductor. This is in a contrast to the FM phase, which is strongly metallic (see Fig. 5b). Our results in CuMnAs illustrate the generic trend that AFMs are more compatible with semiconductor/semimetal electronic structure than FMs. We also performed relativistic GGA+U calculations of magnetic anisotropy of the tetragonal CuMnAs in the AFM ground state. Consistent with the experiment, we find that the Mn moments, ordering ferromagnetically within the plane and antiferromagnetically between the planes, are oriented in the (ab) plane.

Finally, we note that the Ia–Mn–V compounds, NaMnAs and KMnAs, crystallize in a tetragonal structure very similarly to the structure of CuMnAs, with the only difference being that the group Ia element occupies the Mn position in the crystal structure of CuMnAs and Mn occupies the Cu position. For this reason, we compared total energies of CuMnAs when Cu and Mn switched places in the lattice. The lower energy state is observed for the configuration corresponding to the one inferred from our structural measurements of CuMnAs and shown in Fig. 2c, with the difference of 0.56 eV f.u.−1.

Discussion

The results presented above demonstrate CuMnAs as a tetragonal room-temperature AFM with high crystalline quality and chemical fidelity. The exchange bias and broadening of the Fe hysteresis loop induced by CuMnAs provides not only additional evidence for the high-temperature AFM ordering in the tetragonal CuMnAs epilayers, but also clearly demonstrates the potential of this new material for AFM-based spintronic functionalization. The same type of AFM/FM exchange-coupling phenomena was utilized in the demonstration of the first AFM-tunnelling anisotropic magnetoresistance devices1,2,6.

Tunnelling anisotropic magnetoresistance is one of a family of spin-orbit-coupling phenomena in which electronic properties and spin can influence each other. This is manifested in the realization of large magnetoresistances because of anisotropies in the density of states and chemical potential and in the realization of carrier-induced spin reorientation phenomena. Many of these relativistic spintronic effects were first observed in the ferromagnetic Mn-doped zincblende GaAs12,13,34,35. CuMnAs shares the spin orbit, broken space-inversion symmetry character of the electronic structure with these magnetic zincblende compounds, and the tetragonal distortion of its lattice further enhances the magnetic anisotropy. The in-plane lattice parameter closely matches to GaP and Si. CuMnAs layers can be grown in a conventional III–V MBE system and therefore can be readily combined with established semiconductor heterostructure technologies. All these characteristics, complemented by the room-temperature collinear magnetic order, make CuMnAs a favourable material candidate for exploring the phenomena and functionalities within the emerging concept of spintronics based on AFMs.

Methods

Growth of CuMnAs films

The CuMnAs films were grown in a Veeco GEN-III MBE system. Standard effusion cells were used as sources for elemental Mn and Cu, and a valved cracker cell (in the As4 regime) was used as the source of As. The epitaxial growths were carried out at a substrate temperature of 300 °C, with the material fluxes set to cover a range of ratios around the stoichiometric point of Cu:Mn:As=1:1:1. For growth on GaAs(001), after opening the cell shutters the (2 × 4) RHEED image of the clean GaAs surface became diffuse within several seconds. Typically, a new streaky (2 × 2) pattern emerged during the next 20 min (corresponding to ~10 nm of the epitaxial CuMnAs), which then remained stable till the end of the growth. In contrast, for growth on GaP(001), the transition period of the diffuse RHEED image was completely missing; instead, a (1 × 2) pattern arose immediately after opening the shutters, transforming smoothly into the streaky (2 × 2) within the first 5 min of growth (Fig. 2a,b).

Structural characterization

High-resolution XRD measurements to determine the lattice parameters were performed on a Panalytical X'pert Materials Research Diffractometer (MRD) equipped with a standard Cu X-ray tube, a parabolic X-ray mirror, a 4 × 220 Ge Bartels monochromator and a 3 × 220 Ge channel-cut crystal analyser with point detector. For the space group refinement, a combined experimental and computational procedure, which utilized a Bruker area detector in combination with the aforementioned MRD, was used. The details of this technique are published elsewhere27.

High-resolution aberration-corrected STEM was performed at Oak Ridge National Laboratory in a Nion UltraSTEM column, operated at 100 kV and equipped with a fifth-order Nion aberration corrector. Samples were prepared by conventional thinning, grinding, dimpling and Ar ion milling.

Calculations

We used the full-potential linearized-augmented-plane-wave method (WIEN2k package36) to calculate the magnetic structure of the tetragonal CuMnAs. In all calculations we used the PBE-GGA exchange-correlation potential36. Plane wave cutoff was given by RMT × kmax=7.5, where RMT is the smallest muffin-tin radius. Brillouin zone was sampled by a mesh composed of up to 20,000 k-points. To describe the correlations of Mn 3d electrons, we employed the GGA+U approximation with U=4.6 eV and J=0.5 eV. The lattice constants and z positions of As and Mn were taken from the experiment.

Additional information

How to cite this article: Wadley, P. et al. Tetragonal phase of epitaxial room-temperature antiferromagnet CuMnAs. Nat. Commun. 4:2322 doi: 10.1038/ncomms3322 (2013).