Indirect observation of phase conjugate magnons from non-degenerate four-wave mixing

A phase conjugate mirror utilising four-wave mixing in a magnetic system is experimentally realised for the first time. Indirect evidence of continuous-wave phase conjugation has been observed experimentally and is supported by simulations. The experiment utilises a pump-probe method to excite a four-wave mixing process. Two antennae are used to pump a region of a thin-film yttrium iron garnet waveguide with magnons of frequency f1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{1}$$\end{document} to create a spatio-temporally periodic potential. As the probe magnons of fp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{\mathrm {p}}$$\end{document} impinge on the pumped region, a signal with frequency fc=2f1-fp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{\mathrm {c}} = 2f_{1}-f_{\mathrm {p}}$$\end{document} is observed. The amplitude of the nonlinear signal was highly dependent on the applied magnetic field H. Width modes of the probe magnons and standing wave modes of the pump magnons were shown to affect the amplitude of the signal at fc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{{\rm c}}$$\end{document}. Experimental data are compared with simulations and theory to suggest that fc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{{\rm c}}$$\end{document} is a phase conjugate of fp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{{\rm p}}$$\end{document}.


Introduction
The insatiable appetite for smaller, more powerful computing devices is leading to the inevitable breakdown of Moore's Law [29]. A considerable deviation from Dennard's scaling is already underway [18] opening the door for alternative computational paradigms. In recent years, there has been increasing interest in magnon-based computing, or magnonics, [5,7,8,11,12,24,26] as a solution to particular problems [34] facing the future of conventional complementary metal-oxide-semiconductor (CMOS) computing.
An advantage of wave computing not restricted to magnonics is the ability to encode information in two variables: amplitude and phase. The ability to perform operations on these parameters is therefore of fundamental importance to the future of wave computing. The study of nonlinear systems [1] opens the door for the development of these novel operations. One useful operation, in particular, is phase conjugation; a process that exactly reverses the propagation direction and phase factor for every plane wave in an arbitrary wavefront [10]. The result is the creation of a phase conjugate mirror (PCM) that reflects any beam along the same path by which it arrived at the mirror, irrespective of incident angle. This remarkable property can lead to aberration correction of waves after passing through a nonuniform distorting medium; a process with many useful applications including image processing, encryption and spectroscopy [4,9,27].
First experimentally realised in the 1970s, phase conjugate mirrors are a well-established phenomenon in the optical community. A preferred method for creating a PCM is by way of the third-order nonlinear process: four-wave mixing (FWM) [2]. Experiments of this nature have a general form wherein the confluence of two 'pump' beams and a 'probe' beam in a region of nonlinear medium causes the appearance of a fourth beam which is the phase conjugate of the probe signal [3]. In addition to creating a PCM, this method of optical FWM has lead to a number of novel physical processes, including squeezing [17], vortices [32], and the emergence of upto fifth-order nonlinear effects [30].
Until now, phase conjugation of spin waves has only been achieved using methods of parametric pumping and second-order nonlinear processes [20,21,25]. In this work, we report on the creation of a PCM in a magnon waveguide. Our experiments differ from previous work for three reasons: (1) We achieve a phase conjugate with the process of FWM, rather than parametric pumping. (2) Our experiments observe a phase conjugate using the continuous-wave (CW) regime, in contrast to the pulsed experiments that utilise spin wave bullets. (3) In contrast to similar optical experiments [33], we exploit the large magnon beam widths, the strength of the nonlinear interaction and the isotropic magnon dispersion relation to implement a geometry that has orthogonal pump and probe beams.
In general, for FWM to occur certain physical conditions must be met. Energy conservation dictates that where 1,2 are the angular frequencies of the two pump waves, p is the probe wave angular frequency, and c is the conjugate wave angular frequency. Similarly, the condition arises from momentum conservation, where is a magnon wave vector and the subscripts follow the same convention as the angular frequencies. These constraints are met by wave vectors depicted by Fig. 1b. It is clear that if 1 and 2 serve as counter-propagating pumps of same frequency in an isotropic medium and p serves as the probe then by equation (2), the condition p = − c must be satisfied.
We now consider FWM in a nonlinear magnetic system. By examining the Landau-Lifshitz equation with a perturbation expansion of the magnetisation, it has been shown that there exists a driving term for a third-order spin wave that depends on the product of m 3 [16,19], where m is a component of the transverse magnetisation.
To understand the origin of the expected phase conjugate magnon, consider the m 3 term in the location of the pumped region with the probe magnons also present. This term may be expanded into its constituent parts. Following a mathematics analogous to the derivation of optical phase conjugation [31], we expand m 3 and express it as a real quantity: Here m 1,2 and m p represent the amplitude of the transverse magnetisation of the pumps and probe, respectively. Since m i is a complex amplitude, it also contains the phase information. Upon expansion of equation (3), we obtain 56 cross-terms. Terms with phase factors that have combinations of and k that are forbidden by the dispersion relation may be neglected. Of the remaining terms, there is one of particular significance: There are a number of points to note about this term. Firstly, it is proportional to m * p and is therefore the phase conjugate of m p . Secondly, this is the only possible term possessing a wavevector antiparallel to the original probe beam [22,31]. These two properties define the resulting spin wave as a phase conjugate reflection. Finally, we note that this term is proportional to m 2 1 assuming the pumps are of equal magnitude, that is m 1 = m 2 .

Experiment
For the magnon waveguide, we utilised a yttrium iron garnet (YIG) film of thickness 7.8 μm on a gallium gadolinium garnet (GGG) substrate, 2.1 mm wide and 18 mm long, with corners cut at 45 • to minimise reflections. A schematic of the set-up is shown in Fig. 1a. The waveguide was mounted on a printed circuit board (PCB) with three antennae. The antennae were in the formation of a meander structure to suppress the coupling to the ferromagnetic resonance (FMR) mode. The antennae were 4 mm long comprising three legs 40 μm wide, spaced 50 μm apart. The two pump antennae were parallel to the long edge are the pumps, p is the probe, and c is the phase conjugate of the waveguide. The probe antenna was placed 3 mm from the pumps, across the width as shown in Fig. 1a. The probe antenna is used for both transmission and detection of magnons interacting with the pumped region, while the pump antennae excite counter-propagating spin waves that exploit the strong intrinsic nonlinearity of the magnon system to generate a periodic mesoscopic texture. In our experiment, spin waves travelling in both the x and y directions were utilised. For this reason, an isotropic magnon dispersion was required. Such a condition is offered by exciting forward volume magnetostatic spin waves (FVMSWs) which have a uniquely isotropic dispersion, compared to the highly anisoptropic dispersion relations for other magnetostatic spin wave modes [23]. An electromagnet was used to apply the external magnetic field perpendicular to the plane of the film, a field geometry that is necessary for the excitation of FVMSWs.
The pump antennae were excited at f 1 = f 2 = 3.915 GHz by a Hewlett-Packard HP8672A microwave source. A second source (HP8671A) was used to excite the probe antenna at f p = 3.91825 GHz . This detuning Δf = 3.25 MHz was required to use the spectrum analyser to discriminate between input and output signals. Introducing Δf and ensuring Equations (1) & (2) are satisfied leads to a small phase mismatch Δk which manifests itself experimentally as a reduced efficiency [10,22].
The power of the probe magnons immediately after insertion was approximately 8 μW , which attenuated to approximately 2 μW after propagation towards the pumped region, which was populated by spin waves of power approximately 40 μW . The combined spin wave power is above the reported nonlinear threshold for systems satisfying the Lighthill criteria [13,15] as ours does.
The circulator shown in Fig. 2 allows the probe to act as a transmitter-receiver antenna. Reflections from the pumped region will propagate back towards the probe antenna and be detected by the spectrum analyser (ZHL Rhode & Schwarz). The magnetic field was measured with a Hall probe which was connected to a data acquisition computer. Figure 3a, b shows the measured spectra for two different field configurations. When there is no external field applied as in (a), the only signals measured are the input signals. The large signal at f p is due to the impedance mismatch between the microwave transmission line and the probe antenna causing electrical reflections to be detected by the spectrum analyser, while the peak at f 1 is due to direct coupling from the pump antennae to the probe antenna. As the magnetic field was increased to 3077 Oe,  Spectrum measured with applied magnetic field of 3077 Oe. c Spectrum from simulated pump-probe system at 3077 Oe. The large width of the peak at f p is due to limited computational resolution and conversion to a logarithmic scale. Both measured and calculated spectra show peaks at expected phase conjugate frequency f c the spectrum 3b was observed. Of note is the signal at the expected phase conjugate frequency f c = 2f 1 − f p . This is not inconsistent with the notion that at the correct field strength, phase conjugate reflections are occurring from the pumped region. Also present in the spectrum is a term at frequency f = 2f p − f 1 which is due to a third-order process resulting from reflections of probe magnons from the waveguide edge farthest from the probe antenna.

Results and discussion
Investigating further, a sweep of the applied magnetic field was performed, the results of which are shown in Fig. 4. The dashed yellow line marks the field at which the spectrum in Fig. 3b was measured. The strong lines at f p and f 1 appear to have minimal field dependence since the field dependent contribution to the signal is small compared to the electrical response described above.
More interesting, however, is the white line at f c = 3.91175 GHz showing an obvious field dependence. The nature of the dependence of the boxed region can be seen in more detail in Fig. 5. The amplitude reaches the noise floor at a field of approximately 3093 Oe. This is due to the excited FVMSWs approaching the FMR above which no more spin waves are excited. Below this field however, the signal oscillates as a function of field, with a period of approximately 5 Oe. This oscillation corresponds to different standing wave modes across the width of the pumped region. That is, when there are an integer number of halfwavelengths across the width of the waveguide, a standing wave is present and the pump amplitude is larger. This in turn, amplifies the phase conjugate signal, since it scales with m 1 m 2 as shown in equation (4).
In addition to this fast oscillation, Fig. 5 also shows a field dependence on a larger scale. This may be explained by considering width modes [14] of the probe magnons. At approximately 3040 Oe, the width modes excited by the probe antenna destructively interfere, inhibiting the propagation of the probe spin wave, thus diminishing the interaction with the pumped region.

Simulations
To aid with the interpretation of our results, simulations of the experiment were carried out using MuMax3 [28], a micromagnetic simulation software package. Physical parameters were set such that saturation magnetisation, M sat = 197 kA/m, exchange stiffness, A ex = 3.5 × 10 −12 J/m, and Gilbert damping, = 5 × 10 −5 [6]. As in the experiment, the probe antenna was 3 mm from the pumps, while a separate detection region was defined 1 mm from the pumped region. The pumps were driven at 3.915 GHz and the probe driven at 3.91825 GHz with each field configuration simulated for 2 μs before performing a Fourier transform to investigate the frequency response of each system. The secondary effects caused by reflections of spin waves from the short ends of the waveguide were controlled for by increasing the value for in these regions by factor of 300. The regions of waveguide edge directly in contact with the pump antennae were also assigned the increased damping value. By suppressing these edgereflections, any output signals may be attributed purely to the nonlinear interaction of the excited pump and probe magnons. Figure 6 shows amplitudes of pump and probe magnons for different magnetic fields. The white regions represent the simulated waveguides, while the red and blue represent the amplitude of the x component of the transverse magnetisation, m. In the simulations, the probe magnons are launched from an antenna placed at the leftmost edge of the waveguide. At every point along the waveguide, the magnons excited at H 1 = 3051 Oe are weaker than those excited at H 2 = 3082 Oe. It is clear from the many more nodes across the width of the waveguide that (a) shows a higher width mode than in (c), and that its propagation efficiency is reduced. This phenomenon manifests itself  experimentally as highlighted in Fig. 5, where a reduced phase conjugate signal is evident at H 1 compared to H 2 .
The simulation also illuminated the pump behaviour. Figure 6b and d both shows a standing spin wave across the width of the waveguide. At lower field, there are 18 nodes across the width, compared to the 12 nodes present at H 2 . As the field increases from H 1 to H 2 , the intensity of the standing wave will oscillate with every node that is removed. Given Eq. (4) this phenomenon explains the 5 Oe oscillation in Fig. 5, with the simulation matching the experiment well. Indeed, increasing from H 1 to H 2 the calculation shows a difference of 6 nodes, while for the same measured fields, the intensity goes through 6 oscillations.
A time-domain Fourier transform was performed on the simulated data. A typical example is shown in Fig. 3c where the applied field is 3077 Oe as it was for the data measured in (b). The large peak at f p is due to the detection region being placed between the probe and pump antennae, therefore picking up the original probe signal. The width of this peak is an artefact of both the limited computational resolution and the conversion to a logarithmic scale. Also of note is the comparatively small power of the pump frequency at f 1 which is due to leakage of pump magnons, which are generally well confined between the antennae as shown in Fig. 6b and d. We also observe a small bump at 3.925 GHz which arises from a higher order mixing term.
The significance of Figs. 3b and c is a prominent peak at f c in both simulation and measurement. Because the simulated system eliminates reflections from the ends of the waveguide, this signal must necessarily be reflecting from the pumped region of periodic potential. Since this signal was generated by FWM and has angular frequency c = 1 + 2 − p , it must also have wavevector c in order to satisfy equation (4) confirming that it is indeed a phase conjugate of the probe signal. Furthermore, the absence of the peak at 2f p − f 1 in the simulated spectrum supports the notion that it was due to a third-order nonlinear effect caused by reflections from the waveguide edge.

Conclusion
In summary, this work demonstrates through experiments and simulations the generation of a phase conjugate magnon from a non-degenerate four-wave mixing process. The phase conjugate signal is enhanced when the applied magnetic field strength is such that the pumps form a standing wave across the width of the waveguide. This standing wave causes a large pump amplitude which significantly increases the nonlinearity of region. The geometry of the experiment and simulation ensure that any return signal at f c must be a phase conjugate signal.
This new phase conjugate differs from previous observations in magnonic systems in three ways: (1) It utilises CW signals rather than spin wave bullets or pulsed signals (2) It uses a third-order FWM process in contrast to threewave parametric pumping. (3) We utilise the 2D nature of the waveguide with perpendicularly travelling spin waves. Future work would involve the creation of degenerate FWM, though discriminating between monochromatic signals poses different challenges. Our work opens the door for this type of phase conjugation as yet another process that is exploitable in novel magnon-based computational paradigms.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide