sign in sign up
<<
Home , Polariton

arXiv:2004.05823v2 [cond-mat.str-el] 23 Jun 2020 rqec fAMmgo ae ohtesiuiadprobe and stimuli di the more both much makes higher AFM the fields of Furthermore, electromagnetic weaker. the becomes to AFM magnetic of net of a coupling fre- of gigahertz the absence the the moment, in to is due However, material range. FM quency a sublattic in spin magnon to the the due between whereas range frequency interaction terahertz exchange the strong in the is The AFM FM. an a than in dynamics magnon spin faster antif much an has for (AFM) magnet application, spintronics future for promising hs esn,fwwr a eotdo h magnon- the on reported was work few reasons, those fcvt annplrtn[ polariton realizat magnon the cavity on of based a technologies spintronics processing for information proposed been [ have frequency applications resonance Possible cavity the mag- approach external the to via field where tuned netic were resonator cavity resonance microwave magnon (FM) a ferromagnet inside in specifically years recent in ( ef- systems only coupling realized strong was from arising fect splitting mode eminent The ogtm.Hwvr u otedi the to f due However, studied and time. anticipated long been a has material, ordered ically ton, p magnon the dif- , coupling Among light- technology. ferent information and prospe control new up for opening used thus matter, be of a can excitation wave collective electromagnetic polariton, the as between to information matte transfer referred to of being excitation state, hybrid collective The the or di wave levels sta electromagnetic energy hybrid with a produce can or matter of excitation collective and eynro rqec ag nterssrh ad [ bands reststrahl the within in change range reflectance frequency the narrow very from a inferred polari mainly magnon was The ton measurement. optical traditional by seen e hybridization clear elec a and wave, magnet a tromagnetic of spins between condition coupling strong niermgei aeil r xetdt emore be to expected are materials Antiferromagnetic h togculn ewe ih reetoantcwave electromagnetic or light between coupling strong The eeln lr-togmgo-htnculn naplra polar a in coupling magnon-photon ultra-strong Revealing i.e. h yrdlgtadmgo xiaini magnet- in excitation magnon and light hybrid the e.g. eaet esrmn namlieri oa antiferroma polar multiferroic a on is condition measurement coupling terahertz strong format the However, because (AFM) years. ferromagnet recent in pheno cavities The microwave magnon-polariton. inside as named quasiparticle photon w oe tfeunisaoeadblwteelectric-activ the below and above frequencies at modes two h omto fAMmgo-oaio.A lr-togspi ultra-strong An magnon-polariton. AFM of formation the u okrvasta h oeetmgo omto rvnb splitting. driven mode formation polariton magnon detect coherent the that reveals work Our nrylvlsltig oee,teAMmagnon-polariton AFM the However, splitting. level energy ffi h rdtoa Mytimio ant(YIG)) garnet iron yttrium FM traditional the 1 nentoa etrfrQatmMtras colo Phys of School Materials, Quantum for Center International 2 togculn ewe annadeetoantcwv ca wave electromagnetic and magnon between coupling Strong utb ovninleetoismto.For method. electronics conventional by cult eateto hsc n srnm,Uiest fManitob of University Astronomy, and of Department .Y Shi, Y. L. 3 7 olbrtv noainCne fQatmMte,Beijin Matter, Quantum of Center Innovation Collaborative , 12 ]. 1 .Wu, D. ffi ff ut faheigthe achieving of culty rn rmete the either from erent ff c ol ehardly be could ect oantrhrzspectroscopy terahertz domain 1 .X Wang, X. Z. 3 olari- 1 – erro- , 11 ion 1 2 cts es, nd nd or te ]. ]. r. .Lin, T. - - nothdo n erhdo fOwt di with O of tetrahedron and octahedron in iuehge hnteF I nieamcoaecavity microwave a mag- inside of YIG resonance order FM mode one the resonator[ the about than over is higher gap which nitude energy 8%, as split reaches defined the frequency, strength, of coupling ratio spin-photon the The wave. electrom terahertz netic of propagation the in magnon-polariton issihl bv n eo h lcrcatv annfr magnon electric-active below the quency below frequen- and at above modes slightly two cies between phenomenon beating clearly ueeto utfrocplrAMFe mea- AFM terahertz polar time-domain multiferroic present a we on letter, surement this In from dows. responses dynamical the that di is advantage be advantages. particular can certain A has domain measurement trans- frequency domain Fourier in Time By spectrum obtained. the domain. E(t), time sampling the of (EO) in formation electric-optical crystal by nonlinear a detected on be can of E(t) waves duration field time THz electroma with The stable pulses netic are . laser AFM fetosecond of from generated dynamics tool the powerful a detect as to emerged has spectroscopy (THz) terahertz opigi F ytm ofr[ far so systems AFM in coupling iia otepeiu reports[ previous the to similar lcrcplrzto ln h -xs h igecrystal single The c-axis. Fe the along polarization electric hw nFig. in shown en n antzto esrmn.Teltieadmag- Fe and of lattice structure The netic measurement. magnetization di and X-ray terns by characterized were crystals polariton detect to splitting. way mode new a provides technique the domain that time illustrates in- study transform The spectrometer. Fourier (FTIR) using fur- frared fre- by We the measurement in absent domain magnon. is quency magnon-polariton AFM AFM this the that of show ther nature electric-active the to 1 ntFe gnet annfeunybelow frequency magnon e o fmgo-oaio srrl enfra anti- an for seen rarely is magnon-polariton of ion ff ihteavneo lrfs ae ehiu,tm domain time technique, laser ultrafast of advance the With 2 Fe o aiyflle.Hr epeettime-domain present we Here fulfilled. easily not .M Hu, M. C. eo swl tde o ermgei systems ferromagnetic for studied well is menon rn xiain a esprtdi di in separated be may excitations erent Mo sasn ntefeunydmi measurement. domain frequency the in absent is h lrsotTzplepoie e a to way new a provides pulse THz ultrashort the y -htnculn e coupling n-photon c,Pkn nvriy ejn 081 China 100871, Beijing University, Peking ics, 2 Mo 3 ,Wnie,Mntb 3 N,Canada 2N2, R3T Manitoba Winnipeg, a, O 2 edt h omto faculdspin- coupled a of formation the to lead n Mo 3 8 9 O .Teutasrn opige coupling ultra-strong The ]. eegonb hmclvprtasotmethod transport vapor chemical by grown were 8 3 O 2 1 sawl-nw oa F antwt the with magnet AFM polar well-known a is tfroantFe ntiferromagnet T 8 n .L Wang L. N. and efidcerybaigbetween beating clearly find We . N a.Teeaetotpso ein locating ions Fe of types two are There (a). hc eosrtstefraino AFM of formation the demonstrates which , 081 China 100871, g 2 Mo ff c sdrvdbsdo the on based derived is ect 3 O T 8 N hc easg to assign we which , nteAModrdpaeare phase ordered AFM the in 14 ,3, 1, , 13 15 ∼ ∗ ]. .Tersle hexagonal resulted The ]. 2 s h H electric THz The ps. 1 Mo ff ff atos aepat- Laue ractions, 2 3 c a elinked be may ect Mo O ff ff rn iewin- time erent rn magnetic erent 8 3 ytime by O 8 efind We . THz of s ag- g- e- 2

10

4 K

(b)

(a) 4 K 25 K

8

45 K

(a) (b) 60 K c 6

a b 4

濢 (c) 4 K 濠瀂 2

濙濸澳澻濜澼 0 Amplitude (arb. units) (arb. Amplitude

濙濸澳澻濜濜澼 (arb.units) Intensity THz

0 5 10 15

0.5 1.0 1.5 2.0 2.5

Time delay (ps)

Frequency (THz)

FIG. 2. (a) The time traces of the transmitted THz electric field at FIG. 1. (a) Crystal and magnetic structure of Fe2Mo3O8. (b) The temperature dependent magnetization measurement for H//c-axis at several selected temperatures. The lines have been vertically shifted 0.2 T. for clarity. (b) Fourier transformation of the time trace of all the time range at 4 K, and (c) Fourier transformation of the time trace at4K after the main peak (in the time window of dashed rectangular from 5 ps to 17 ps). moments. The octahedrons and tetrahedrons share their cor- ners to form a honeycomb lattice in the ab-plane. Figure 1 (b) shows the magnetic susceptibility as a function of temper- from the reversible energy exchange between THz wave and ature for H//c-axis measured in a Physical Property Measure- magnon mode, therefore it is an indication of the Rabi oscilla- ment System under the magnetic field of 0.2 T. The data are tion from hybridized magnon-polariton state. We can get the similar to the reported results [16], indicating an AFM phase spectrum in frequency domain through Fourier transformation transition at TN = 60 K with ordered moments along c-axis. of E(t) to study the interaction between the magnon mode and Due to the lack of inversion symmetry, an non-zero electric THz electromagnetic wave through the beating frequencies. polarization is present in the compound. A considerable elec- There are two ways to get the spectrum in the frequency do- tromagnetic effect was reported in the multiferroic state below main. One is to perform Fourier transformation of the E(t) in TN [16, 17]. the whole time delays. In this way we get the spectrum of We performed time domain THz measurement on the crys- a complete absorption, including the major absorption in the tals in a home built spectroscopy system with the terahertz main peak and the absorption in the coherent collective exci- beam normal to the ab-plane of the crystals. Details of THz tation process after the main peak. The other one is to do the system setup and measurements are presented in the Supple- Fourier transformation of the time domain spectrum after the mental Materials [18]. Figure 2 (a) shows the waveforms of main peak. Then, the result is dominated by coherent collec- the transmitted THz electric field at several selected tempera- tive modes. We perform Fourier transformation in both ways. tures. Below TN, we find long live oscillations after the main The results for the lowest temperature at 4 K are shown in peak, which corresponds to the coherent magnon mode at fre- Fig. 2 (b) and (c), respectively. In the frequency domain spec- quency 1.25 THz. This collective excitation mode was studied trum transformed from all time delays, we can see only one in previous reports, however the time domain THz waveform dip at frequency of 1.25 THz, related to the absorption of the within a sufficiently long time delay was not presented [19– magnon excitation. In the frequency spectrum transformed 21]. The polarization dependent measurements indicated that from the dashed rectangular in Fig. 2 (a), which is the time it is an electric active magnon of the antiferromagnetic spin window excluding the main peak, we can see two split peaks system[19]. In general, the magnon excitation can be induced at 1.2 and 1.3 THz, respectively. Such beating phenomenon byaneffective interaction through dS i/dt = [S i, He f f ], where was rarely seen before. To the best of our knowledge, there He f f reflects the interaction (e.g. Heisenberg interaction). was only one work reporting similar beating in TmFeO3 [23]. Here, the magnetic excitation possesses in-plane (or perpen- There were two related studies on the hybridized coupling dicular to the c axis) oscillation of electric polarization, which between magnon mode of Fe3+ and the paramag- 3+ was proposed to be induced by the inverse Dzyaloshinskii- netic resonance of Er in ErFeO3 [24] or simultaneously hy- Moriya interaction and/or single-site anisotropy [19]. Be- bridized coupling of magnon--polarition in specially cause the mode is electric active, i.e. driven by the electric designed LiNbO3-ErFeO3 chip devices [25]. field rather than the magnetic field of the THz electromag- For a comparison, we also performed a frequency domain netic wave, the magnon strength is considerably higher than transmission spectroscopy measurement by using a FTIR the usual AFM magnon triggered by the magnetic in- spectrometer Bruker 80V for the sample at the lowest tem- teraction, where coherent oscillations are usually observed af- perature 8 K. The absorption spectrum is shown in Fig. 3. ter subtracting the waveform above the N´eel temperature [22]. We find a single narrow absorption peak at 1.25 THz. Within In our experiment, a clear beating phenomenon from two the signal-to-noise ratio and spectral resolution, no splitting is modes close in energy was observed in the time domain spec- detected. trum below TN as presented in Fig. 2 (a). The signal inten- The different observations in the time domain THz and sity increases as the sample temperature decreases. We shall FTIR measurements are unexpected. It is of great interest to elaborate later that the beating can be explained to originate explore the underlying physics. As we mentioned in the intro- 3

6 usually in THz frequency range. For this reason, the con- ventional cavity method based on the electronic technique 4 in the microwave frequency regime is not applicable. Obvi- ously, new techniques have to be developed to study the phe-

2 ) (arb. units) (arb. )

w nomenon. In the present time domain study, when the THz (

a pulse generated by the ultrafast laser pulse transmits through

0

0.8 1.0 1.2 1.4 1.6 the sample, as illustrated in Fig. 4 (a), it transforms energy Frequency (THz) to the spin system and coherently excites the magnon. So the recorded main peak contains major absorption at the magnon FIG. 3. The high resolution absorption spectrum measured by frequency, the spectra in other frequencies would be transpar- Fourier transform spectrometer. ent for the sample. Because the processions of spins in the sample are triggered by the THz pulse, they have the same phase at the same wavefront of THz pulse, namely, their pro- duction, the light matter interaction can lead to a hybrid state. cession phase is locked to the phase of THz pulse. The os- Nevertheless, detecting the hybridized energy gap would de- cillations after the main peak in Fig. 2 (a) are the THz emis- pend on the coupling strength. Because the coupling between sions from the coherent spin processions, i.e. magnon excita- electromagnetic wave and magnetic excitation through mag- tion in the time domain. The emitted THz wave frequency is netic dipole interaction is usually weak, it is not easy to ob- equal to the magnon frequency, which can be absorbed again serve hybrid excitation by traditional optical techniques. In by the spin system to excite magnon. In this way the energy recent years, much progress has been made for realizing the circling is built up, and a strong coupling condition is estab- magnon polarition in FM systems inside microwave cavities. lished. While the waves will eventually escape, there is no The cavities were used to block the escape of microwave and emitted THz signal at magnon frequency at sufficiently long to limit the energy of microwave around the magnetic matter. time delay, then we can expect that overall absorption in the In such cavity magnet system, the magnon excitation absorbs whole time delay comes from the magnon, as shown in Fig. 2 the energy from the microwave inside the cavity. The excited (b). But in a proper time window of picoseconds time scale magnon will decay by re-emitting [26]. While the a beating emerges as a consequence of the strong coupling, photons are limited in the cavity, they can be absorbed by the reflecting two split modes with close frequencies, as in Fig. magnon again, and this leads to an energy circulation between 2 (c) . This is similar to a case of light interacting with a magnon and microwave. In a real system, the energy may two-level system, resulting in a periodic change between ab- escape from the circulation. The imperfect mirrors will loss sorption and emission of photons, which leads to Rabi split- photons from the cavity, and the magnon may decay by other ting. Such coupled spin-photon system can be described by channels. If the energy exchange between magnon and mi- the two-level Hamiltonian [7, 9, 25]. crowave is more efficient than the escape of the microwave and the energy loss in the matter, the system is in a strong qc 0 0 0 0 Ω/2 coupling regime, resulting in cavity magnon polariton. H = Hp+Hm+Hint = ~ +~ +~ . 0 0! 0 ω0 iδ! Ω/2 0 ! − (1) (a) (b) Here, Hp and Hm are the bare Hamiltonians of photon and

g magnon, Hint denotes their coupling, q the photon wavevector, THz c the speed of light in vacuum, ω0 the magnon frequency, δ the S width of the magnon mode, Ω the splitting energy gap. The g eigenfrequencies are

+ qc ω0 δ 1 2 2 ω = i Ω + ((qc ω0) + iδ) . (2) ± 2 − 2 ± 2 q − As a result, the coupling leads to the splitting of two branches FIG. 4. (a) A schematic diagram for the magnon-light interaction. s between the magnon and the THz wave, i.e. magnon- denotes the spin. A strong coupling is realized under the case when polariton. The splitting frequency at the resonant frequency Ω the energy transformation between magnon and THz photon (g) is is the Rabi frequency, being related to the coupling strength. more efficient than the escape of the THz wave and the magnon decay at a rate of γ. (b) A plot for the of the AFM magnon A plot of the dispersion of the coupled magnon-polaritonfrom polariton from the real part of Eq. (2). The wave vector is defined Eq. (2) with ω0=1.25 THz is shown in Fig. 4 (b). The two as q=1/λ. The red curves indicate the coupled state. The two dashed dashed lines in Fig. 4 (b) indicate the dispersion of photon and lines indicate the dispersion of free photon and magnon. The two magnon without coupling. The light-matter coupling leads modes exhibit avoided crossing when approaching to the resonance to the avoided crossing at resonant frequency, as shown by frequency. the red line. The spin-photon coupling strength can be esti- mated from the ratio of the split energy gap over the mode For an AFM, the magnon appears at much higher energy, resonance frequency, and it reaches 8%. This value is signif- 4 icantly higher than the value for FM YIG inside a microwave splitting does not depend on the pulse power of THz radia- cavity resonator[9], reflecting an ultra-strong coupling effect. tion. We performed power dependent measurement on a sepa- rate sample with different THz radiation powers in a THz time

domain system with amplified laser system and confirmed that 4 (a) (b) splitting does not change (the results are presented in the sup- 3

Sample 2 Sample 2 plementary file [18]). 2

1 Sample 1 To summarize, we observed an ultra strong coupling effect

Sample 1 0 between the THz electromagnetic wave and the electric ac- Amplitude (arb. units) (arb. Amplitude

-1

THz Intensity (arb. units) (arb. Intensity THz tive AFM magnon in Fe2Mo3O8. The level splitting or Rabi 0 5 10 15 0.5 1.0 1.5 2.0 2.5 oscillation is invisible in frequency domain measurement but Time delay (ps) Frequency (THz) can be monitored by time domain THz spectroscopy experi- FIG. 5. (a) THz waveform for two samples with thickness of 1 mm ment in a proper time window. A clear beating between two (sample 1) and 0.4 mm (sample 2) in time domain. (b) The Fourier modes at frequencies slightly above and below the electric- transformed spectra in the frequency domain for the two samples. active magnon frequency is observed below TN. We elab- The beating period and signal strength vary with the thickness of the orate that the coherent formation of the spin procession en- samples. The dash lines indicate the splitting of the modes for the ables the energy circulation in the propagation of the THz 1mm thickness sample. A smaller splitting is seen for the thinner sample. electromagnetic wave in the picoseconds time scale, thus es- tablishing the strong coupling between the AFM magnon and THz light without a cavity. The coupling strength is related On the other hand, the splitting is not observed in FTIR ex- to the assemble of spins, therefore, the thickness of samples. periment. In such experiment the continuous waves from the Our study further illustrates that the electric active magnon far-infrared light source, i.e. mercury lamp, transmit through in multiferroics with strong electromagnetic effect is a par- the sample, and the far-infrared bolometer detector measures ticularly suitable candidate to achieve the strong coupling of the power spectrum within a certain integration time. We be- AFM magnon polarition, and the THz time domain technique lieve that the coupling is still described by the same Hint, how- offers new route to directly monitor the energy gap splitting. ever, the light waves emitted by a mercury light source are incoherent in nature. Even though a similar beating should This work was supported by National Natural Sci- be formed in a time window of picoseconds for a specific ar- ence Foundation of China (No. 11888101), the National ray of wave passing through the sample, the average of differ- Key Research and Development Program of China (No. ent arrays of waves emitted by a mercury light source would 2017YFA0302904, 2016YFA0300902). give a zero amplitude in the oscillations. Additionally, the re- sponse and the integration time of bolometer detector is much longer than the time scale of tens of picoseconds, therefore, the Fourier transformation of the recorded far-infrared signal ∗ [email protected] would give only a single absorption from the AFM magnon. [1]D. L. Mills and E. Burstein, This explains why the strong coupling polariton effect cannot Reports on Progress in Physics 37, 817 (1974). be observed in FTIR measurement. [2] M. R. Jensen, S. A. Feiven, T. J. Parker, and R. E. Camley, In fact, the coupling strength is enhanced for a assemble Journal of Physics Condensed Matter 9, 7233 (1997). [3] O. O. Soykal and M. E. Flatt´e, of N spins system. The coupling strength is known to be en- Phys. Rev. Lett. 104, 077202 (2010). hanced by a factor of √N [4, 9, 24, 27]. In the present study, [4] H. Huebl, C. W. Zollitsch, J. Lotze, F. Hocke, M. Greifen- the number of spins interacting with the THz electromagnetic stein, A. Marx, R. Gross, and S. T. B. Goennenwein, wave transmitting through the sample should depend on the Phys. Rev. Lett. 111, 127003 (2013). thickness l of the sample. Thus, we expect that the splitting of [5] Y. Tabuchi, S. Ishino, T. Ishikawa, R. Yamazaki, K. Usami, and the modes Ω √l. To testify the relation between the number Y. Nakamura, Phys. Rev. Lett. 113, 083603 (2014). of spins and the∼ coupling strength, we performed the experi- [6] X. Zhang, C.-L. Zou, L. Jiang, and H. X. Tang, Phys. Rev. Lett. 113, 156401 (2014). ff ments on samples with di erent thicknesses. Figure 5 shows [7] L. Bai, M. Harder, Y. P. Chen, X. Fan, J. Q. Xiao, and C.-M. the THz spectra for two samples with thickness of 1 mm and Hu, Phys. Rev. Lett. 114, 227201 (2015). 0.4 mm in both time and frequency domains at 4 K. Indeed, [8] P. Hyde, L. Bai, M. Harder, C. Dyck, and C.-M. Hu, we observe a clear increase of beating period in the time do- Phys. Rev. B 95, 094416 (2017). main and a reduction of the splitting in frequency domain for [9] C. Match, M. Harder, L. Bai, P. Hyde, and C.-M. Hu, the thinner sample. The ratio of the splitting energy Ω of the Phys. Rev. B 99, 134445 (2019). two samples is about 0.6, being consistent with the ratio of the [10] D. Zhang, X. Q. Luo, Y. P. Wang, T. F. Li, and J. Q. You, Nature Communications 8, 1 (2017). square root of thickness 0.63. The results confirm the pos- ∼ [11] M. Harder and C. M. Hu, itive correlation between spin number and coupling strength, State Physics - Advances in Research and Applications, one of the representative characteristics of magnon polariton 1st ed., Vol. 69 (Elsevier Inc., 2018) pp. 47–121, coupling physics. We would also like to emphasize that the arXiv:arXiv:1601.06049v1. 5

[12] H. Huebl and S. Goennenwein, Physics 8, 51 (2015). Phys. Rev. Lett. 119, 077206 (2017). [13] M. Białek, A. Magrez, and J. P. Ansermet, [21] S. Yu, B. Gao, J. W. Kim, S.-W. Cheong, M. K. L. Physical Review B 101, 024405 (2020). Man, J. Mad´eo, K. M. Dani, and D. Talbayev, [14] P. Strobel, Y. L. Page, and S. McAlister, Phys. Rev. Lett. 120, 037601 (2018). Journal of Solid State Chemistry 42, 242 (1982). [22] L. Y. Shi, Y. Q. Liu, T. Lin, M. Y. Zhang, S. J. [15] P. Strobel and Y. L. Page, Zhang, L. Wang, Y. G. Shi, T. Dong, and N. L. Wang, Journal of Crystal Growth 61, 329 (1983). Phys. Rev. B 98, 094414 (2018). [16] Y. Wang, G. L. Pascut, B. Gao, T. A. Tyson, [23] K. Grishunin, T. Huisman, G. Li, E. Mishina, T. Rasing, A. V. K. Haule, V. Kiryukhin, and S.-W. Cheong, Kimel, K. Zhang, Z. Jin, S. Cao, W. Ren, G. H. Ma, and R. V. Scientific Reports 5, 12268 (2015). Mikhaylovskiy, ACS 5, 1375 (2018). [17] T. Kurumaji, S. Ishiwata, and Y. Tokura, [24] X. Li, M. Bamba, N. Yuan, Q. Zhang, Y. Zhao, M. Xiang, Phys. Rev. X 5, 031034 (2015). K. Xu, Z. Jin, W. Ren, G. Ma, S. Cao, D. Turchinovich, and [18] See Supplemental Material for the detailed information of the J. Kono, Science 361, 794 (2018). THz system setup and additional time domain measurements [25] P. Sivarajah, A. Steinbacher, B. Dastrup, J. Lu, M. Xi- with different excitation powers. ang, W. Ren, S. Kamba, S. Cao, and K. A. Nelson, [19] T. Kurumaji, Y. Takahashi, J. Fujioka, R. Ma- Journal of Applied Physics 125, 213103 (2019). suda, H. Shishikura, S. Ishiwata, and Y. Tokura, [26] K. Yamaguchi, T. Kurihara, Y. Minami, M. Nakajima, and Phys. Rev. B 95, 020405 (2017). T. Suemoto, Phys. Rev. Lett. 110, 137204 (2013). [20] T. Kurumaji, Y. Takahashi, J. Fujioka, R. Ma- [27] H. Y. Yuan and X. R. Wang, suda, H. Shishikura, S. Ishiwata, and Y. Tokura, Appl. Phys. Lett. 110, 082403 (2017).