We describe measurements of the uniaxial magnetic anisotropy and spin-flop rotation of the Néel vector in antiferromagnetic CuMnAs thin films using neutron diffraction. The suppression of the magnetic (100) peak under magnetic fields is observed for films as thin as 20 nm, indicating that they undergo a spin-flop transition. Good agreement is found between neutron diffraction and electrical transport measurements of the spin-flop rotation in the same layer with similar shape and hysteresis of the obtained curves, while the neutron measurements provide a quantitative determination of the spin flop extending throughout the antiferromagnet layer.
There is growing interest in antiferromagnetic (AF) materials for applications in spintronics. Predictions and demonstrations of current-induced switching of the AF Néel vector have stimulated the development of all-electrical AF memory devices with potential for high speed operation, robustness against external fields, and multi-level neuromorphic outputs.1–5 Electrical readout of the AF state typically utilizes anisotropic magnetoresistance (AMR) or spin Hall magnetoresistance, but such effects are small and can easily be obscured by structural changes induced by Joule heating.6 Direct determination of the Néel vector is, therefore, important. This has been achieved by imaging local modifications of AF domains and domain walls affected by current pulses or magnetic fields using x-ray photoemission electron microscopy7–11 and more recently by optical birefringence,12 spin-Seebeck,13 or magneto-Seebeck14 microscopies.
Neutron diffraction is a well-established technique to directly measure the atomic arrangement of magnetic moments in a bulk AF crystal. It provides a quantitative determination of the Néel vector and, unlike electron microscopy techniques, is routinely performed in high magnetic fields. It is, therefore, well-suited to study spin-flop transitions, which typically occur in magnetic fields of several Tesla. In AF thin films and devices, spin-flop rotations of the Néel vector are invaluable for characterizing extraordinary magnetotransport11 or magnon transport phenomena15,16 and exploring multi-stability of magnetic memory states.17,18 However, there are few direct studies of spin-flop rotation in AF thin films using neutron diffraction due to the very small volume of the sample film compared to a typical neutron experiment, which is a flux limited technique.
Here, we demonstrate that neutron diffraction can be used in an unambiguous total volume determination of the spin-flop rotation in CuMnAs films with thickness down to 20 nm. We also show good agreement between neutron diffraction and electrical transport measurements performed on the same sample. CuMnAs is a collinear antiferromagnet with a Néel temperature of around 480 K.19 It can exhibit pronounced modifications of its AF domain structure in response to applied electrical or optical pulses, which may be accompanied by significant electrical readout signals.5 Understanding of its magnetic and transport properties is, therefore, important for the design of spintronic memory devices.
The CuMnAs 20 and 45 nm films were grown on GaP(001) substrates using molecular beam epitaxy.20 X-ray diffraction measurements confirmed the tetragonal crystal structure of the films, illustrated in Fig. 1(a), with the CuMnAs[001] axis perpendicular to the plane and the CuMnAs principal axes at to the substrate principal axes. Previous neutron diffraction and x-ray magnetic linear dichroism measurements on similar films have shown that the magnetic moments align in the ab plane with a competition between a uniaxial CuMnAs[010] anisotropy and a biaxial CuMnAs[110]/[10] anisotropy.19 The uniaxial in-plane anisotropy arises because of the substrate-film interface, where the symmetry is broken by the P-terminated unit cell of the GaP substrate.
(a) CuMnAs unit cell with arrows representing the expected uniaxial anisotropy along the [010] crystal direction. (b) Schematic diagram of the neutron diffraction experiment showing the horizontal scattering plane and the vertical magnetic field. The incident beam is represented by the black arrow, and the detector arrays covering in-plane scattering angles from 10° to 170° are depicted as gray semi-circles. (c) and (d) Sketches showing the magnetic moment orientations for an uniaxial antiferromagnet with respect to an applied magnetic field below (c) and above (d) the spin-flop field Hsf. (e) Plots of the (100) peak for the 45 nm film with three different field strengths applied along the easy axis. The integration regions used in background subtraction are shown.
(a) CuMnAs unit cell with arrows representing the expected uniaxial anisotropy along the [010] crystal direction. (b) Schematic diagram of the neutron diffraction experiment showing the horizontal scattering plane and the vertical magnetic field. The incident beam is represented by the black arrow, and the detector arrays covering in-plane scattering angles from 10° to 170° are depicted as gray semi-circles. (c) and (d) Sketches showing the magnetic moment orientations for an uniaxial antiferromagnet with respect to an applied magnetic field below (c) and above (d) the spin-flop field Hsf. (e) Plots of the (100) peak for the 45 nm film with three different field strengths applied along the easy axis. The integration regions used in background subtraction are shown.
The neutron diffraction measurements utilized the time of flight cold neutron beam line WISH at the ISIS facility.21 Figure 1(b) represents the geometry of the experiment with a horizontal scattering plane and a vertical magnetic field generated by a superconducting split-pair solenoid magnet. For both films, the samples were mounted in two orientations: in one case, the uniaxial easy axis ([010]) was parallel to the applied magnetic field, and in the other, the field was parallel to the [100] axis. They corresponded to the (h0l) and (0kl) scattering planes. The position of the samples in respect of the incident beam was chosen to optimize the neutron flux for the (100) and (010) magnetic reflections as well as to access the (001) reflection at the same time. It was achieved when the surface plane of the samples was at to the beam direction. The measurements were conducted at a sample temperature of 100 K, and the sample position was the same for both films investigated.
For a magnetic field applied parallel to the easy axis of a uniaxial antiferromagnet [Fig. 1(c)], there is a threshold field strength known as the spin flop field, , which causes a 90° reorientation of the antiferromagnetic Néel vector. Above this field, the moments are canted around the external field direction [Fig. 1(d)]. Due to higher-order magnetic anisotropies as well as domain walls and other inhomogeneities, the sharp transition at is typically broadened into an s-shaped curve.
Nuclear contribution to the (100) and (010) reflections is forbidden in tetragonal CuMnAs (P4/nmm space group), and therefore, they contain only magnetic scattering, which is sensitive to the component of magnetic moments perpendicular to the scattering vector. Consider the case where the [010] direction is a uniaxial easy magnetic axis, and the field is to be applied along this direction. In the zero field state where , the magnetic moments are parallel to the easy axis, and the (100) peak is observable. Under a field , the moments are perpendicular to the easy axis and aligned along the [100] direction with small-angle canting of the moments in the direction of the field, and therefore, the (100) peak is extinguished. Hence, the magnitude of the peak provides a direct measurement of the spin flop transition in the sample. If only a partial rotation of the magnetic moments occurs, the magnitude of the peak will reduce in proportion, so the ratio of the peak amplitude at field compared to the original state will give a direct measurement of the extent of spin flop in the sample.
The (100) peak for a 45 nm CuMnAs film is shown in Fig. 1(e) for three different magnetic field strengths. The peak is significantly reduced at 1 T applied field compared to its zero-field amplitude and has almost vanished at 8 T. To quantify the spin-flop transition, the measured intensity is integrated across the peak and background regions indicated in Fig. 1(e). The combined magnetic and structural (001) peak is also measured to allow calibration of the alignment and detection efficiency. Because the moments are expected to stay within the (001) plane, the (001) peak is unaffected by spin flop and is used as a control measurement. Consistent with this, the peak amplitude is found to be independent of the magnetic field within the experiment uncertainty (see the supplementary material).
Plots of the ratios of integrated peak amplitudes vs magnetic field for 45 and 20 nm thick CuMnAs films are shown in Fig. 2. For the 45 nm sample [Fig. 2(a)], the (100) peak is larger than the (010) peak by a factor of 5.4 ± 0.9. This indicates that the magnetic moments in the sample are predominantly aligned with the [010] axis. A sharp decrease in intensity is observed between zero and 1 T for the (010) peak, and between zero and 2 T for the (100) peak, due to the spin flop rotation of the Néel vector into an axis which is perpendicular to the field. In high fields, the (100) peak is decreased to 2.2% ± 0.7% of its zero-field amplitude, showing conclusively that the vast majority of the CuMnAs layer undergoes spin flop in this field range. After the high field measurements, the sample recovers its original state. This observation in combination with the well-defined critical fields strongly suggests that the reorientation of the magnetic moments is a flop-transition rather than due to continuous moving domain walls.
Ratio of in plane (IP) with the sample surface to out-of-plane (OOP) integrated peak amplitudes vs magnetic field for CuMnAs films of thicknesses (a) 45 and (b) 20 nm. The red circles and blue diamonds correspond to the measurements of the (100) and (010) in-plane peaks with the field along the [010] and [100] axes, respectively. Open and filled symbols represent decreasing and increasing field sweeps.
Ratio of in plane (IP) with the sample surface to out-of-plane (OOP) integrated peak amplitudes vs magnetic field for CuMnAs films of thicknesses (a) 45 and (b) 20 nm. The red circles and blue diamonds correspond to the measurements of the (100) and (010) in-plane peaks with the field along the [010] and [100] axes, respectively. Open and filled symbols represent decreasing and increasing field sweeps.
The 20 nm CuMnAs layer [Fig. 2(b)] exhibited a (100) to (010) peak amplitude ratio of 10.5 ± 3.4 and, therefore, demonstrates even stronger anisotropy than the thicker layer. As with the thicker layer, the (100) peak intensity drops sharply between 1 and 2 T applied fields. However, there is significant intensity remaining between 2 and 6 T, which may be attributed to the presence of magnetic domain walls and defect-driven domain pinning. The intensity then drops to zero between 6 and 8 T. Due to time constraints and the weak signal from this thin film, the (010) peak was not measured over the whole magnetic field range.
Complementary electrical transport measurements were conducted on the 45 nm sample, which was patterned into an L-shaped Hall bar structure as illustrated in Fig. 3(a). This allows currents to be applied simultaneously along two in-plane crystal axes, at and to an applied magnetic field, respectively, while recording the longitudinal and transverse voltages and . Each arm of the structure is of width m. The magnetoresistance measurements were conducted with the sample submerged in liquid helium at 4 K inside a superconducting magnet. Magnetic fields were applied either parallel to [010] or [100] with applied current Iprobe = 2 mA.
(a) Schematic of an L-bar device. The orange represents the 100 μm wide CuMnAs bar, and the gold squares represent contact pads to which wires are bonded. The arrows represent the current flow in the device, and the voltages drawn correspond to the voltages described in the text where 1 and 2 correspond to the first and second arm the current flows through, respectively. (b) Transverse AMR and (c) longitudinal AMR, as defined in Eqs. (1) and (2), respectively, vs the magnetic field for an L-bar device fabricated from the 45 nm thick CuMnAs film. Red and blue lines are for fields along the [010] and [100] axes, respectively. The plots are offset for clarity to show zero AMR at zero field.
(a) Schematic of an L-bar device. The orange represents the 100 μm wide CuMnAs bar, and the gold squares represent contact pads to which wires are bonded. The arrows represent the current flow in the device, and the voltages drawn correspond to the voltages described in the text where 1 and 2 correspond to the first and second arm the current flows through, respectively. (b) Transverse AMR and (c) longitudinal AMR, as defined in Eqs. (1) and (2), respectively, vs the magnetic field for an L-bar device fabricated from the 45 nm thick CuMnAs film. Red and blue lines are for fields along the [010] and [100] axes, respectively. The plots are offset for clarity to show zero AMR at zero field.
Due to anisotropic magnetoresistance (AMR), the longitudinal and transverse voltages depend on the angle between the applied current and the orientation of the Néel vector. Neglecting crystalline terms, this dependence is given by and . Hence, the rotation of the Néel vector from to with respect to the current direction should generate no change in the longitudinal voltages and a maximum change in the transverse voltages.
Figures 3(b) and 3(c) show the measured normalized transverse and longitudinal anisotropic magnetoresistances, defined as
where the factor of 2 for accounts for the number of squares between the longitudinal voltage contacts. For both field directions, a broad s-shaped change in is observed below 2 T. A hysteresis is observed between the increasing and decreasing field sweeps, but the values return to the same point at zero field. The variations in are much smaller and are ascribed to small variations in the sample temperature during the field sweep.
Figure 4 compares the neutron diffraction and measurements over the same magnetic field range. The spin flop extent is defined as , where P(H) is the integrated peak amplitude at field H and is the initial value at zero field. It can be seen that the electrical transport measurement gives a qualitatively accurate measure of the spin flop extent with the shape, hysteresis, and spin flop field all matching well. Note that the neutron diffraction measurements were performed at 100 K while the AMR was measured at 4 K. Previous studies have shown that this difference in temperature creates a small change in amplitude of the AMR signal due to changes in the base resistance but does not qualitatively change the signal.11
(a) Spin flop extent extracted from neutron diffraction and b) transverse AMR for the 45 nm thick CuMnAs film with an external magnetic field applied along the [010] axis. Arrows indicate the direction of the magnetic field sweep.
(a) Spin flop extent extracted from neutron diffraction and b) transverse AMR for the 45 nm thick CuMnAs film with an external magnetic field applied along the [010] axis. Arrows indicate the direction of the magnetic field sweep.
We have demonstrated the use of magnetic neutron scattering as a tool for probing the antiferromagnetic spin-flop regime in CuMnAs. From our results, we were able to show that 98% of the magnetic order is reoriented under the action of an applied magnetic field along the magnetic easy axis in a strongly uniaxial sample. Any remnant signal could be attributed to domain walls and domain pinning due to defects. Complementary electrical transport measurements were congruent with the neutron scattering signal in both shape and hysteretic behavior. Overall, we propose that this technique should be applicable for measuring the magnetic Néel vector magnitude and direction of any such compatible magnetic crystals.
See the supplementary material for plots of the peak amplitudes that were used in the normalization process to make Fig. 2. It also includes waterfall plots of the diffraction peaks used to find the amplitudes in the aforementioned plots.
This work was supported by the EU FET Open RIA via Grant No. 766566 and the UK Engineering and Physical Sciences Research Council (Grant No. EP/P019749/1). We wish to acknowledge ISIS Neutron and Muon Source for the beamtime at WISH.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Stuart Forbes Poole: Formal analysis (lead); Investigation (equal); Methodology (supporting); Visualization (lead); Writing – original draft (lead); Writing – review and editing (supporting). Peter Wadley: Conceptualization (lead); Formal analysis (supporting); Funding acquisition (equal); Methodology (supporting); Supervision (lead); Validation (supporting); Writing – review and editing (supporting). Luke Barton: Investigation (equal); Methodology (supporting); Visualization (supporting). Mu Wang: Formal analysis (supporting); Investigation (equal); Methodology (supporting); Visualization (supporting); Writing – review and editing (supporting). Pascal Manuel: Data curation (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Writing – review and editing (supporting). Dmitry Khalyavin: Data curation (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Writing – review and editing (supporting). Sean Langridge: Conceptualization (supporting); Data curation (equal); Investigation (equal); Methodology (equal); Project administration (equal); Supervision (supporting); Writing – review and editing (supporting). Kevin William Edmonds: Conceptualization (supporting); Funding acquisition (equal); Project administration (supporting); Supervision (supporting); Writing – original draft (supporting); Writing – review and editing (lead). Richard Paul Campion: Resources (equal). Vit Novak: Resources (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request. The data that support the findings of this study are openly available in https://doi.org/10.5286/ISIS.E.RB1920313 at Science and Technology Facilities Council, Ref. 22.