A. T. Costa,Mikhail I. Vasilevskiy,J. Fernández-Rossier,Nuno M. R. Peres
Strongly Coupled Magnon–Plasmon Polaritons in Graphene-Two-Dimensional Ferromagnet Heterostructures
magnons,plasmons,polaritons,FMR,2D magnets
Magnons and plasmons are different collective modes, involving the spin and charge degrees of freedom, respectively. Formation of hybrid plasmon–magnon polaritons in heterostructures of plasmonic and magnetic systems faces two challenges, the small interaction of the electromagnetic field of the plasmon with the spins, and the energy mismatch, as in most systems plasmons have energies orders of magnitude larger than those of magnons. We show that graphene plasmons form polaritons with the magnons of two-dimensional ferromagnetic insulators, placed up to to half a micrometer apart, with Rabi splittings in the range of 100 GHz (dramatically larger than cavity magnonics). This is facilitated both by the small energy of graphene plasmons and the cooperative super-radiant nature of the plasmon–magnon coupling afforded by phase matching. We show that the coupling can be modulated both electrically and mechanically, and we propose a ferromagnetic resonance experiment implemented with a two-dimensional ferromagnet driven by graphene plasmons.
Strongly Coupled Magnon–Plasmon Polaritons in Graphene-Two-Dimensional Ferromagnet Heterostructures

Magnons and plasmons are different collective modes, involving the spin and charge degrees of freedom, respectively. Formation of hybrid plasmon–magnon polaritons in heterostructures of plasmonic and magnetic systems faces two challenges, the small interaction of the electromagnetic field of the plasmon with the spins, and the energy mismatch, as in most systems plasmons have energies orders of magnitude larger than those of magnons. We show that graphene plasmons form polaritons with the magnons of two-dimensional ferromagnetic insulators, placed up to to half a micrometer apart, with Rabi splittings in the range of 100 GHz (dramatically larger than cavity magnonics). This is facilitated both by the small energy of graphene plasmons and the cooperative super-radiant nature of the plasmon–magnon coupling afforded by phase matching. We show that the coupling can be modulated both electrically and mechanically, and we propose a ferromagnetic resonance experiment implemented with a two-dimensional ferromagnet driven by graphene plasmons.

Magnons are the elementary excitations of every magnetically ordered system, governing their low-energy properties. Magnons have attracted renewed interest for several reasons. They can transport spin currents for applications in nondissipative spintronics, host topological order with chiral edge states, form exotic collective states such as Bose condensates and spin superfluids; most important for the scope of this work, they can couple to photons. Magnons play a particularly important role in two-dimensional (2D) magnets, as their uncontrolled thermal proliferation prevents long-range order. Thus, most prominent examples of 2D ferromagnets (2DFM), such as VI3, CrI3, and Fe3GeTe2, have a sizable gap in the magnon energy spectrum. Experimental techniques that are very successful in producing and probing magnons in bulk ferromagnets are not easily adaptable to 2D systems due to the intrinsically small sample volume. For instance, the sensitivity of the ferromagnetic resonance is limited by the ratio between sample and detector sizes. Recent proposals, such as ferromagnetic resonance force spectroscopy, address the challenge of probing submicrometer-size samples but are a long way from monolayer van der Waals magnets. Cavity magnonics has also emerged as a way of enhancing the coupling between exciting/probing fields and the magnetic sample. Rabi splittings on the order of 100 MHz have been obtained for micrometer-sized spheres on resonant microwave cavities. Further enhancement in coupling strength, leading to Rabi splittings of a few GHz, has been achieved for macroscopic-sized ferromagnets in optical and superconducting cavities. In this context, exciting and probing magnons efficiently in 2D ferromagnets remain a challenge. There are three main bottlenecks for the existing techniques. One is having a driving field of the right frequency: magnons in 2DFM have frequencies in the range ∼0.25–1 THz, whereas the highest frequencies achieved in ferromagnetic resonance (FMR) experiments are ∼700 GHz. This stems from a combination of the scarcity of microwave sources of higher frequencies and the need to match the resonance frequency of a cavity. This brings forward the second challenge, the strength of the photon–magnon coupling. The interaction of the magnetic field of light with matter is notoriously much weaker than that of the electric field. Placing the ferromagnetic sample in a resonant cavity enhances the coupling between the magnon and the cavity modes. The frequencies of those modes, however, decrease as the cavity volume increases, whereas the enhancement factor goes in the opposite direction. There is, thus, a compromise between enhancement factor and resonance frequency that limits the sensitivity of setups of this kind. This links to the third challenge, which is detector sensitivity. Again, this is limited by the smallness of light’s coupling to magnetic dipoles and puts a constraint on the minimum enhancement factor needed. With regard to the frequency of the driving field, graphene plasmons come to mind as prime candidates. Their frequencies can be tuned essentially continuously, by gating graphene away from charge neutrality. Current experimental limits on such control set the spectral range of graphene plasmons to a few THz within the wavelength range of interest to us. Graphene plasmons have been shown to form various kinds of polaritons in van der Waals heterostructures. Coupling to graphene plasmons has been proposed recently as a way to probe collective excitations in superconductor surfaces, 2D superconductors, and excitons in insulators. The common theme of those works is the coupling between the strongly confined electric field associated with the graphene plasmon and the charges of the electrons in the nearby system. The coupling to spins is more subtle, since it relies on the much smaller magnetic-dipolar nature. It is known, however, that momentum and frequency matching can enhance dramatically the coupling between light and an ensemble of quantum objects. With this in mind, we have studied the coupling between graphene plasmons and 2D magnons in a van der Waals heterostructure. Long-range ferromagnetic order in 2D is only possible in the presence of magnetic anisotropy, on account of the Mermin–Wagner theorem. Spin–orbit coupling breaks spin rotation symmetry, stabilizes long-range magnetic order, and opens up a gap in the magnon spectrum at zero wave vector, q = 0. In many cases of interest, the magnon gap in 2D ferromagnets is much larger than typical values in 3D. For instance, the magnon gap of CrI3 monolayers, one of the most prominent 2D magnents, has been reported to be in the 0.3–1.0 meV range. For some materials this value can exceed 5 meV, putting the lowest energy magnon in the terahertz region. On the other hand, the energy and wave vector of graphene plasmons may be tuned to match those of magnons in a 2DFM by adjusting the charge density of the graphene sheet. Thus, van der Waals heterostructures composed of 2DFM and plasmonic materials, such as graphene, may provide a platform to bridge the terahertz gap in optoelectronics. Previous attempts in this direction have been aimed at the coupling between light and the orbital magnetic moments of electrons in conducting materials, but here we focus on the spin magnetic moment, which is associated with quantum magnetism. We consider a van der Waals heterostructure, depicted schematically in Figure 1, composed of a 2D ferromagnet with an off-plane easy axis and a graphene sheet, separated by a dielectric, such as hexagonal boron nitride, of thickness z and relative dielectric constant ε. The 2D ferromagnet is described with a spin Hamiltonian in the linear spin wave approximation. The in-plane magnetic field of the plasmon is coupled to the local spins of the ferromagnet via a Zeeman interactionwhere R⃗l is the 2D vector marking the position of unit cell l in the 2D ferromagnet, z is the vertical distance between the graphene sheet and the 2D ferromagnet, and σ⃗ are the dimensionless Pauli spin matrices that relate to the spin angular momentum through . Using the expression for the quantized field of the graphene plasmon given in ref (31), the Zeeman interaction with the transverse-magnetic (TM) plasmon magnetic field readswhere aq⃗† is the creation operator for a plasmon with wave vector q⃗ parallel to the graphene sheet. We note that the plasmon magnetic field lies in-plane, so that it generates a torque on the static magnetization. At the microscopic level, this entails the creation of magnons. The coupling strength F(q,z) is given bywhere ℏωpl(q) is the energy of a plasmon with wave vector q⃗, , Λ(q) is the mode length of the plasmon (see the Supporting Information), and A is the area of the graphene sheet. The plasmon field decays exponentially, but for the range of wave vectors relevant to this work the decay length is of the order of several micrometers, thus presenting no practical concern. We note that the coupling strength of a plasmon mode with wave vector q with an atomic spin S⃗l is vanishingly small, as it scales with the inverse of √A. In contrast, for collective excitations such as magnons, it makes sense to transform the spins to a plane-wave basiswhere k⃗ is a wave vector in the Brillouin zone of the 2DFM and N is the number of unit cells. After applying this transformation to eq 2 we obtainCompared to the case of atomic spins, the magnon–plasmon coupling is enhanced by a factor √N, where N is the number of spins, resulting in a Rabi coupling that no longer depends on system size, as N ∝ A. Thus, magnon–plasmon coupling is enhanced due to the phase-matching of the plasmon field to a macroscopic number of phase-locked precessing spins. The quantized Hamiltonian for plasmons in graphene readswhere their energy dispersion curve is given byHere, EF is graphene’s Fermi energy, ε is the average dielectric constant of the two media surrounding the graphene sheet, α is the fine structure constant, c is the speed of light, and q is the plasmon’s (in-plane) wave vector. To study the effect of plasmon–magnon coupling, we adopt a description of magnons in terms of linearized Holstein–Primakoff (HP) bosonswhere σ̂+,– are the ladder operators acting on the spin located at site l; their magnitude S is assumed to be the same throughout the whole material. The operators bl† and bl respectively create and annihilate a localized spin flip excitation at site l. Assuming translation symmetry in the 2D ferromagnet, we can rewrite the HP bosons in reciprocal space. Then, the Hamiltonian for bare magnons has the formThe wave vectors k⃗ span the Brillouin zone of the 2D ferromagnet. The function ℏωmag(k⃗) is the dispersion relation for the bare magnons. For small momenta, we have ℏω ≃ ℏω0 + ρk2, where the first term is the magnon gap and the second provides the dispersion due to the exchange-driven spin stiffness ρ. For plasmons with energies ℏωpl ≈ 1 meV (thus close to that of uniform magnons in typical 2DFM) and typical graphene doping levels (EF ≈ 100 meV), q ≲ 0.1 μm–1. This is tiny compared to the linear dimensions of the magnon Brillouin zone (∼104 μm–1), so that the dispersion of the magnon states is negligible in that wave-vector window. After transforming the Zeeman Hamiltonian to the HP representation in reciprocal space, it readsThe coupling strength is given bywhere z is the distance between the graphene sheet and the 2DFM, N is the number of spins in the 2DFM, and k(+) ≡ kx + iky. The function F(k,z) has been defined in eq 3. Notice that the plasmon–magnon coupling is diagonal in wave vector, meaning that each bare plasmon of wave vector k⃗ couples only to magnons with the same wave vector. We notice in passing that the quantized treatment of both plasmons and magnons that we adopt here can be used to model intrinsically quantum phenomena, such as spontaneous emission of a plasmons due to a spin flip, or the dynamics of non-Glauber states, which are relevant for cavity quantum magnonics. These are out of reach of the classical description adopted in previous works. If the terms proportional to b–k⃗ak⃗ and bk⃗†ak⃗† in eq 10 are neglected, the remaining Hamiltonian can be mapped onto a single-particle problem, leading to approximate analytic forms for the dispersion relations of the two hybrid plasmon-magnon modeswhere . This equation predicts a gap opening of magnitude ℏΩk⃗(z) at the crossing frequency where the plasmon–magnon detuning ω– vanishes. In the following we treat the complete magnon–plasmon Hamiltonian, including the nonconserving terms ba and a†b†, by analyzing the plasmon Green function (see the Supporting Information for details), which can be probed in near-field optical experiments. Since the plasmon–magnon coupling is linear, the equations of motion for all Green functions can be solved analytically. Their explicit expressions are given in the Supporting Information. Here we will highlight the most relevant features by plotting the plasmon spectral density, , that is of course affected by the coupling to magnons. In Figure 2a is shown the spectral density for the case where the magnon gap has ωmag(0) = 3 meV, corresponding to a 2DFM such as Fe3GeTe2. In that figure we also show the dispersion curves for the bare plasmon and magnon. The formation of a plasmon–magnon polariton with a Rabi splitting larger than 100 GHz, dramatically larger than the values reported in cavity magnonics, is apparent. Interestingly, the magnitude of the Rabi coupling Ω can be tuned mechanically, by controlling the graphene–ferromagnet distance z, as we show in Figure 2b. In this energy range the plasmon decay rate in the direction perpendicular to the graphene layer is small, which means that the plasmon–magnon coupling is sizable even for graphene–2DFM distances on the order of 1 μm, where interlayer exchange is completely negligible. The Rabi coupling can be further tuned electrically, controlling the graphene Fermi energy EF with a back gate, as we show in Figure 2b,c for three different 2D ferromagnets. For a given magnon energy, there is an optimal value of EF that maximizes the Rabi coupling strength, as we show Figure 2d. We thus see that in a wide range of experimentally relevant parameters, the intrinsic magnon–plasmon Rabi coupling due to Zeeman coupling can be the larger than 50 GHz. The estimated Rabi coupling is a lower bound, coming from the intrinsic Zeeman interaction, and additional contributions to the Rabi coupling can occur when the magnon anisotropy gap is sensitive to the plasmon electric field. We now propose to take advantage of the magnon–plasmon coupling to carry out ferromagnetic resonance of monolayers using an attenuated total reflection setup (see Figure 3a). Exciting plasmons directly with optical beams is impossible due to the kinematic mismatch between plasmons and propagating light. By placing a prism of a high dielectric constant material on top of the hBN layer, it is possible to generate evanescent waves within the hBN that will excite the surface polaritons of the heterostructure. Whenever the in-plane component of the wave vector of light matches that of a polariton with the same frequency, there is a dip in the reflected intensity. The in-plane wave vector can be controlled via the incidence angle. With this setup, it is possible to excite polaritons whose wave vectors and frequencies lie between the light cones inside hBN and the dielectric of which the prism is made. Germanium, for instance, would be a convenient material to use for the prism. It is transparent to electromagnetic radiations of frequencies below 1 THz and its relative dielectric constant within the same frequency range is εGe ≈ 16. In Figure 3a we plot the spectral density for a magnon gap of ℏω0 = 3 meV and graphene gating voltage corresponding to EF = 500 meV. The dispersions for light inside hBN and germanium are plotted as dashed lines, to mark the spectral region probed by the experiment. The black dot in Figure 3b marks the point at which the plasmon–magnon spectral density is probed. By applying an external magnetic field perpendicular to the structure, we shift the magnon energy, thereby changing the spectral density and, consequently, the device’s reflection coefficient. In Figure 3c we plot the spectral density, for a fixed wave vector and frequency, as a function of the external magnetic field. The sharp peak heralds the magnetic nature of the polariton being probed in this setup. Finite plasmon and magnon lifetimes will broaden the resonance peak and reduce the contrast in the reflection coefficient of the device we propose. Recent experimental results indicate that graphene plasmon-polaritons can have lifetimes in the vicinity of 40 ps at a temperature of 60 K. Magnons in bulk CrI3, on the other hand, can have lifetimes of up to 2 ns for temperatures below the ferromagnetic transition (∼60 K). To assess the impact of finite excitation lifetime on our results, we introduced a phenomenological broadening to the bare plasmon and magnon Green functions, to mimic the experimentally estimated lifetimes. The result is shown in Figure 3d. As expected, the resonance peak is considerably broader than before, but the contrast is still large enough to allow detection of the resonance, provided the range of magnetic fields is wide enough. Thus, our proposal would permit tackling the three issues that make FMR in 2D magnets challenging and could open a new venue to explore collective spin excitations in 2D systems.