稀土正铁氧体中THz自旋波的相干调控与强耦合研究进展*

金钻明 阮舜逸 李炬赓 林贤 任伟 曹世勋‡ 马国宏† 姚建铨

1)(上海大学理学院物理系,上海 200444)

2)(天津大学精密仪器与光电子工程学院,天津 300110)

1 引 言

对于反铁磁材料而言,相邻原子的交替磁矩导致了净磁矩为零,使得反铁磁体对外界磁场不敏感[1,2].同时,它同样导致了反铁磁材料中没有漏磁场.相比于铁磁体,这些性质能够极大地减小储存1比特信息所需的最小体积.另一方面,反铁磁体对外加磁场的不敏感性使得有效操控和探测反铁磁体的自旋态极具挑战性.克服这个挑战的关键是如何有效地激发反铁磁体中的自旋共振模式以及在皮秒量级的时间尺度上实现磁矢量的超快翻转[3-5].基于以上特点,反铁磁材料在自旋电子学中引起了广泛的兴趣.

基于圆偏振的超短激光脉冲所携带的角动量与反铁磁NiO晶体间相互作用的逆法拉第效应(inverse Faraday effect,IFE),Satoh等[6]观测到频率在1.07 THz和0.14 THz的THz辐射.他们指认高频模式为面外磁模式,低频模式为面内磁模式[7].2012年,Nishitani等[8]利用延迟时间可调谐的飞秒激光双脉冲,实现了面外磁模式的光学相干控制.IFE是一种非线性光学效应,属于冲击受激Raman散射效应.2015年,Mikhaylovskiy等[9]研究了通过交换相互作用实现光脉冲电场与自旋的耦合.通过对铁氧体(具有倾角自旋结构)的自旋共振超快光激发所产生的太赫兹(THz)发射,实验上证实了各向同性的光磁效应—逆磁折射效应.根据THz辐射强度可以估算,能量密度约为1mJ/cm2的激光脉冲可以实现交换相互作用的亚皮秒调制,相当于一个0.01 T的脉冲有效磁场.

超短激光脉冲直接照射到材料上,不可避免地存在“热”效应,光诱导的电子动力学过程在一定程度上会调制激光所诱导的自旋波.理想情况下,自旋动力学的超快激发和调制是“非热”的,不希望有任何形式的“热”加载到材料或器件上[10].除了利用飞秒激光脉冲诱导反铁磁晶体中的自旋波,近年来THz电磁脉冲也被广泛地应用于反铁磁材料的自旋动力学,尤其是自旋模式的共振激发[11].2011年,Kampfrath等[12]利用强THz脉冲的磁场分量直接激发反铁磁NiO的自旋极化波(面外磁模式),并通过同步飞秒脉冲的磁光Kerr效应探测该磁模式的进动与弛豫动力学.此外,利用时间延迟的双THz脉冲,实现面外磁模式的相干控制.

畸变钙钛矿结构的稀土正铁氧体RFeO3(R为Y和稀土元素)属于正交晶系,空间群为Pbnm,每个单胞中有四个等价的Fe3+离子.Fe3+与O2－形成FeO6八面体,稀土离子填充在八面体形成的空隙中.由于晶格畸变,Fe—O—Fe之间的键角小于180°,RFeO3属于倾角反铁磁体.其独特的磁性能表现为铁的磁矩来自3d电子,稀土元素的磁矩来自4f电子.稀土正铁氧体中发现的激光诱导超快自旋重取向和IFE产生的有效磁场实现自旋波调控,使RFeO3成为凝聚态和材料物理中研究的热点[2].2010年Suemoto课题组[13]研究了THz辐射实现YFeO3单晶准铁磁和准反铁磁模式的选择性激发和相干控制.2013年,Yamaguchi等[14]利用THz时域光谱研究了 E rFeO3中的自旋重取向相变过程.近年来关于RFeO3的超快自旋重取向,以及自旋波的相干控制研究激发了研究兴趣[15-21].自旋重取向以及超快磁相变不仅具有重要的学术价值,同时必将加速超快光磁效应在自旋电子器件中的应用.

结合本课题组最近几年的研究工作,本文主要介绍国内外在利用THz时域光谱研究反铁磁晶体中的自旋动力学的进展.首先总结了RFeO3晶体THz自旋波的共振激发以及THz光谱应用于自旋重取向温区的标定.简述了利用THz脉冲实现自旋波的相干控制.介绍了上海大学和Rice大学合作在实验上证实了固体中的Dicke协作耦合[22],并定量给出了 F e3+-Er3+对称和反对称的交换耦合常数.这些常数不仅有助于理解3d-4f的磁耦合,也是许多化合物中诸如新的磁相变、磁电效应、电磁振子、非线性自旋激发等奇异现象产生的原因.最后,介绍了上海大学与荷兰Radboud大学合作,在TmFeO3单晶中观测到THz光子与反铁磁共振模式的免交叉现象,证明了RFeO3单晶的法布里-珀罗腔可以实现THz光子与反铁磁自旋波模式的强耦合效应,产生光子与磁振子的杂化元激发-磁极化子.

2 基于瞬态Zeeman转矩的磁共振激发

如图1(a)所示,RFeO3为正交的钙钛矿结构,用空间群表示[23].最强的磁相互作用是Fe3+-Fe3+交换相互作用,导致在约650 K以下Fe3+的自旋为反铁磁有序.反对称的Dzyaloshinskir-Moriya(DM)相互作用和磁各向异性产生的综合效应使 F e3+亚晶格产生了很小的倾斜,偏离完美的反平行构置,从而导致一个宏观的磁化矢量M,RFeO3表现出弱铁磁性.

图1 (a)RFeO3反铁磁晶体晶体结构与自旋结构,邻近的 F e3+ 离子沿着a轴反平行有序排列;(b)THz激发的准铁磁模式(qFM)和准反铁磁(qAFM)模式Fig.1.(a)Lattice and spin structure of RFeO3,adjacent Fe3+ions are antiparallel and ordered along the a axis;(b)THz excitation of qFM mode and qAFM mode.

对RFeO3磁结构的描述,通常采用方向相反的两套铁离子的亚晶格M1和M2,且|M1|=|M2|=M0.分别定义“铁磁矢量”(FM)M=M1+M2和“反铁磁矢量”(AFM)L=M1-M2.M和L需同时满足[24]:

任何电磁波辐射不仅有电场分量也携带有磁场分量,当电磁波与磁有序物质相互作用时会发生强烈的耦合,这在人工电磁材料[25]和置于微波腔中的铁磁性YIG[26]中表现得十分明显.对于频率范围为0.1—10 THz的THz电磁辐射,相应的能量跨度在0.4—40 meV范围.当THz脉冲入射到晶体中,THz脉冲的磁场分量类似于一个亚皮秒时间尺度的磁脉冲.宏观上,THz瞬态磁场(HTHz)的方向正交于材料的宏观磁化矢量M,必将产生一个作用于M的瞬态Zeeman转矩,T∝M×HTHz,这一瞬态转矩T使得反铁磁介质中的磁化矢量M偏离平衡位置.当瞬态磁场HTHz消失后,M将围绕材料的有效磁场以Larmor频率进动,并最终弛豫到平衡位置.在磁化矢量M的进动过程中,旋转的磁偶极子会辐射圆偏振的电磁波,ETHz∝∂M2/∂t2,此电磁波的频率等于磁化矢量的进动频率,这一过程称为自由感应衰减.对于反铁磁材料来说,磁共振模式的频率可以简单表示为其中各向异性场HA< 1 T,而交换场Hex≈100 T.相比于铁磁体中的ωF≈γ ·HA,两者大小相差了两个数量级.

通过THz时域光谱,利用线偏振的THz脉冲可以选择性地激发 F e3+系统中的两类磁子共振模式.当HTHz垂直于M,THz波激发RFeO3的准铁磁模式(quasi FM-mode,qFM).当HTHz平行于M,THz波激发RFeO3的准反铁磁模式(quasi AFM-mode,qAFM).图1(b)给出了THz波磁场分量对RFeO3中的qFM和qAFM模式的激发示意图.对于 F e3+的 qFM 模式,M1和M2同相位振动,然而它们之间的夹角保持不变.因此简化为M绕着a轴进动.然而,qAFM模式则表现为M沿着a轴以 “应切”方式辐射线偏振的自由感应衰减信号.

以浮区法生长的b切 D yFeO3单晶为例.样品双面抛光,厚度为 0.72 mm.在 Néel温度(TN=645 K)以下[27],D yFeO3晶体在c轴方向表现出宏观的弱铁磁性.图2(a)和图2(b)分别是当入射的THz脉冲的电场矢量垂直于a轴(ETHz⊥a,HTHz⊥c)和电场矢量平行于a轴(ETHz‖a,HTHz‖c)情况下,透过 D yFeO3单晶后的THz脉冲时域波形.实验测得THz波形包括两个部分,一个是在零延迟时间处的宽带THz主脉冲部分,二是THz主脉冲后指数衰减的长寿命振荡分量.这部分振荡信号突出显示于图2(a)和图2(b)的插图中.qFM和qAFM模式的振荡信号的周期分别为2.63 ps和1.92 ps.图2(c)和图2(d)为振荡分量(时域范围从22—36 ps)的傅里叶变换光谱.两个共振峰分别位于0.38 THz和0.52 THz,分别对应于DyFeO3单晶的qFM(Ma,Mb和Lc的振荡频率)和qAFM(La,Lb和Mc的振荡频率)自旋共振模式.通过改变HTHz与M间的夹角,可有效调节两种磁共振模式的幅度[28].

图2 (a),(b)分别为激发qFM模式(红色)与qAFM模式(蓝色)时的THz时域透射谱,插图为振荡部分的放大图;(c),(d)分别为振荡部分的傅里叶变换光谱Fig.2.(a),(b)THz time-domain transmission spectrum of qFM mode(red curve)and qAFM mode(blue curve),the insets are enlarged versions of the oscillating sections;(c),(d)Fourier transform spectra of the oscillating parts.

表1 RFeO3的qFM和qAFM模式的自旋共振频率Table 1. The qFM and qAFM mode resonance frequencies of rare earth orthoferrite.

表2 SmxDy1—xFeO3单晶在40K的qFM模式的自旋共振频率与Sm浓度的关系[35]Table 2. Resonance frequencies of qFM mode for the single crystals versus Sm concentration at 40 K[34].

对于不同的稀土元素,R-Fe之间的相互作用能也就不同,RFeO3的相变过程及相变温区也不尽相同,表3总结了不同稀土正铁氧体的自旋重取向 温 区,如 Y bFeO3[37],N dFeO3[16],T bFeO3[38],HoFeO3[33],ErFeO3[14],T mFeO3[3],D yFeO3,S mFeO3[34]等.此外,如表4所列,通过改变稀土元素共掺杂的比例,也可以有效地调控重取向温区.在温度诱导的自旋重取向过程中,有的样品只发生二级相变过程(Γ4↔Γ2),宏观磁矩只在晶体的ac面内发生旋 转 ,如 N dFeO3[16],ErFeO3[14],S mFeO3[34],TmFeO3[39],Y bFeO3[40]等.值得注意的是,DyFeO3样品发生一级相变(Γ4↔Γ1)[41].由于RFeO3共振吸收对THz的磁场分量的偏振方向敏感,如果R为非磁性金属,如Y或 P r3+,不含4f电子,则RFe之间的相互作用可以忽略,几乎不会产生温度诱导的Fe离子的自旋重取向过程.然而,利用其他外部条件,如外加磁场可以实现 Y FeO3单晶的自旋重取向过程[42].THz时域光谱为研究磁场调控的 N dFeO3自旋重取向磁相变动力学提供了研究方法[43].

图3 (a),(b)40 K和300 K时 ErFeO3 的THz时域透射谱,插图为振荡部分的放大图;(c),(d)分别为振荡部分的傅立叶变换光谱,插图为 Γ 2 与 Γ 4 的示意图Fig.3.(a),(b)THz time-domain transmission spectra of ErFeO3 at 40 K and 300 K;(c),(d)Fourier transform spectra of the oscillating signals.Insets: schematic diagram of Γ 2 and Γ4 .

表3 稀土正铁氧体的自旋重取向温区总结Table 3. The spin reorientation temperature region of rare earth ferrite.

表4 SmxDy1—xFeO3的自旋重取向温度与Sm浓度的关系Table 4. SRT temperatures for the SDFO single crystal family versus Sm concentration.

3 自旋波的相干控制

相干控制主要有两种途径,一是采用多束不同频率且具有一定相位关系的激光同时作用在物质上,形成特定的干涉场,以控制系统演化的路径及速率.另一种是利用系统在退相时间内能态间的相位关系,采用两个或两个以上相同频率且具有一定时间间隔的脉冲,控制系统的演化.采用第二种方法,通过一对具有时间延迟的飞秒脉冲激发ZnTe晶体,产生两个THz 脉冲序列.当第二个THz脉冲与第一个THz脉冲相位相同,此时所激发自旋波的相位也相同,将相干地加强qFM或qAFM模式; 当第二个THz脉冲与第一个THz脉冲产生的自旋波的相位相反,此时第二个THz控制脉冲会相干抵消第一个脉冲产生的效应,实现对自旋波的相干相消.全THz相干操控技术可以在反铁磁晶体中选择性地激发或湮灭任意模式的自旋波,从而实现自旋的相干操控.

图4(a)—(c)所示为双THz脉冲激发下得到的THz透过 D yFeO3单晶后的时域波形图,两个THz脉冲的间隔时间为 Δt.图4(a)和图4(b)分别 对应 qFM和qAFM模式.当ω·Δt/π=3(奇数),两个自旋波将相干相消,如图4(a)和图4(b)中的黑色曲线所示.当THz脉冲的电场分量与样品的a轴呈45°时,可以同时观察到0.38 THz和0.52 THz两个频率,表明单个THz脉冲同时激发qFM和qAFM模式.采用双THz脉冲技术,改变Δt为 qFM模式振荡周期的1.5倍(Δt=1.5TqFM)时,如图4(c)中的蓝色曲线,此时第二个THz脉冲将抑制qFM模式,同时加强qAFM模式,如图4(d)中蓝色曲线所示; 当 Δt=1.5TqAFM时,qAFM模式被抑制,qFM模式得到加强,如图4(d)中红色曲线所示.

此外,利用 Y FeO3晶体本征的介电各向异性,通过选择THz脉冲和晶体光轴间的夹角θ,可实现单THz脉冲在晶体中的相干控制.图5(a)为在不同方位角θ的情况下,THz波透过 Y FeO3晶体后时域谱上的电场分量.θ定义为THz的偏振方向与晶体b轴间的夹角.θ=0°和θ=90°分别对应于THz的电场分量沿着b轴和a轴.Y FeO3晶体的双折射会导致沿着晶体的b和a轴的THz脉冲间存在时间延迟.当θ=45°时,THz脉冲将历经最大的各向异性.从实验数据得到 Y FeO3晶体的双折射 Δn=na-nb≈0.2 .图5(b)为 Y FeO3晶轴和入射THz脉冲偏振方向间不同夹角时的振荡信号,实线为通过指数衰减的正弦函数的拟合结果.实验结果表明,准铁磁共振信号(0.3 THz)的振幅和相位都强烈依赖于方位角θ,FID信号的振幅改变了近一半,相位移动了近π/2 .

图4 (a)THz脉冲激发qFM模式;(b)THz脉冲激发qAFM模式;(c)THz脉冲同时激发qFM和qAFM模式,当Δt为qFM(qAFM)振荡周期的1.5倍时,该自旋进动被有效地抑制;(d)图(c)中振荡部分的傅里叶变换光谱Fig.4.(a)THz pulses excite qFM mode;(b)THz pulse excited the qAFM mode;(c)THz pulse excites both qFM mode and qAFM mode,as the interval time is 1.5 times of the qFM(qAFM)oscillation period,the spin precession is suppressed;(d)Fourier transform spectra of the oscillating parts in(c).

图5 当 θ=0°,45°,90°时,透过样品后THz波的电场强度,θ 的定义如插图所示;(b)振荡部分的放大,其时域区间为10-25 ps范围内的THz电场,实线是单指数衰减拟合;(c),(d)qFM模式自旋振荡部分的傅立叶变换得到振幅和相位随方位角 θ 的变化[17]Fig.5.(a)THz electric fields transmitted through the YFeO3,as θ=0°,45° and 90°;(b)the oscillating parts of the(a)from 10 ps to 25 ps,the solid line is a single exponential decay fitting;(c)amplitude and(d)phase varies with the angle θ ,by using the Fourier transform of the spin oscillating of qFM mode [17].Reproduced with permission from Ref.[17].

4 稀土元素的THz激发跃迁

除了对 Fe3+的自旋波激发与相干控制,Er3+离子(4 f11)占据了晶体中低对称性的位置.晶体场将导致 Er3+离子的能级劈裂,形成了Kramers双重态,每一个双重态由一对时间反演态组成.当在低温下,施加一个外加的静态磁场HDC,最低的两个时间反演态之间的跃迁可以解释为离子的电子顺磁共振(EPR).在经典物理图像中,EPR代表自旋绕HDC的拉莫进动.

图6(a)描述的是Er3+离子在多重效应下的能级图.自由空间中的Er3+离子能级标记为2S+1LJ,其中S和L分别是总的自旋和轨道量子数,J=L+S为总的角动量量子数.由于晶体结构具有低的对称性,除了Er3+离子的Kramers简并,晶体场效应使有关磁量子数m的能级退简并.图中用|i=1〉,|i=2〉,|i=3〉代表每一个晶体场能级,它们由一对时间反演态组成.即使在零磁场下,由于破坏时间反演对称性的 F e3+-Er3+磁交换作用,他们也将发生交换劈裂成为非简并的.进一步沿着不同的晶轴施加外磁场HDC,出现不同的Zeeman劈裂图案.

图6(b)为不同磁场下c切 ErFeO3的吸收光谱,实验条件为HTHz‖b,HDC‖c,T=30 K.Er3+的吸收谱线可以分成两类.图6(b)中标记为“A”的谱线在低的磁场下处于THz系统的可探测带宽之下.随HDC的增加而发生蓝移,在高的外磁场下变得十分清晰.谱线“A”是 Er3+离子的 EPR,对应于|1-〉→|1+〉能级跃迁.在吸收谱上并没有观察到|2-〉→|2+〉的跃迁(图6(c)中的“B”谱线),这是由于当前的实验温度30 K,|i=2〉能级上还没有足够的布居数.此外,图6(b)中标记为1,2,3,4的一组谱线描述的是,零磁场下位于1.4 THz处的强吸收带在磁场HDC作用下分裂成多个峰,如图6(c)所示.这一组谱线代表的是两个Kramers双重态之间的晶体场跃迁,即|1-〉→|2-〉跃迁(线 1),|1+〉→|2-〉跃迁(线 2),|1-〉→|2+〉跃迁(线 3),和|1+〉→|2+〉跃迁(线4).此外,在HDC=0时,中心频率为1.9 THz的晶体场跃迁谱线属于|i=2〉→|i=3〉跃迁.图中,0.45 THz的吸收谱线(以白色虚线标记)是 F e3+的准FM磁子.由于极强的磁各向异性,它在外加磁场的作用下并不发生移动.除了 Er3+离子,在 T mFeO3[39]和PrFeO3中都观察到了稀土离子由于THz诱导的晶体场跃迁.

图6 (a)ErFeO3中的Er3+离子由于多重效应的能级分裂示意图;(b)从0 T到10 T不同磁场下的吸收系数谱,白色虚线为Fe3+的qFM磁振模式;(c)计算得到不同磁场下晶体场的双重态|i=1〉和|i=2〉 [22]Fig.6.(a)Energy level splitting scheme of Er3+ ions due to multiple effects;(b)absorption coefficient spectra at various magnetic fields from 0 T to 10 T,the white dashed line is the F e3+ qFM magnon mode;(c)calculated energy levels for the|i=1〉 and|i=2〉crystal-field doublets as a function of magnetic field[22].

5 磁相互作用中Dicke协作耦合

N个二能级原子与单个光场的相互作用是量子光学中研究最广泛的多体问题,Dicke首次分析了其在超辐射背景下的相互作用.这类系统的一个共同特征是其耦合强度加强了倍.美国Rice大学Kono教授和上海大学曹世勋教授团队合作,发现在 Y3+掺杂的 ErFeO3单晶样品(ErxY1-xFeO3)中存在类似于典型的光与物质相互作用的效应.以物质-物质相互作用形式出现协作耦合加强,实验结果表明找到了自旋与真空磁子的关联交换耦合的实验证据.这与在一个光学腔中,二能级原子系统与真空光场耦合类似.

如图7所示,实验系统研究了Y3+掺杂比例、温度、外加磁场依赖的THz吸收光谱.图7(a)—(k)所示为不同温度 ErxY1-xFeO3晶体的吸收光谱,实验构置为HTHz‖b,HDC‖c.可以明显观察到在吸收光谱上,磁子与EPR的免交叉行为.当未耦合的磁子和EPR模式的频率相等时,定义此时外加磁场的大小为零调谐磁场H0,此时的频率为零调谐跃迁频率ω0.当磁场为H0时,光谱中杂化分支间的频率分裂称为拉比劈裂Ω(H0),大小等于磁子-EPR耦合强度Λ的两倍.Λ随温度的降低而增加,表明Fe3+-Er3+耦合变得更强.

为了系统研究Λ如何随晶体中 Er3+自旋密度的变化,实验中采用的是 Y3+掺杂的ErxY1-xFeO3体系.先前的文献已经报道,通过改变比例系数x,非磁性的 Y3+离子能简单地减少 Er3+的自旋密度,且并不改变样品的晶体结构和磁结构.如图7(f)—(h),(i)—(k)显示的是x=0.75(x=0.5)样品温度依赖的免交叉光谱线,同样发现 Λ 随温度的降低而增加的规律.此外,比较图7(c),7(f),7(i)可见,当样品温度为20 K,Er3+离子比例高的样品具有更大的Λ.当样品温度为10 K(图7(d),7(g)和7(j))和5 K(图7(e),7(h)和7(k)),都能观察到掺杂浓度的依赖关系具有相同的变化趋势.值得注意的是,x=0.5的样品体系在温度为5 K时具有最大的Λ/ω0=0.18,此时的材料系统进入超强耦合范畴[44,45].实验发现,Λ正比于两个量乘积的平方根,即式中的ηspin是对EPR产生贡献的 Er3+净自旋密度(与T,H0和x有关),ωFM是Fe3+的qFM磁子频率.

如图7(l)插图所示,ηspin由两种机制所决定.首先,低温下的 Er3+离子可以近似认为是二能级系统,遵循玻尔兹曼分布.温度为T时,自旋向上(|1+〉)和自旋向下(|1-〉)态的热力学分布遵循居里定律,其中〈µEr〉是每个Er3+自旋的平均磁矩,µsat是每个Er3+自旋的饱和磁矩; 相应于T=0时的情形,此时所有的Er3+离子占据|1-〉态;ΔE(T,H0)是|1+〉和|1-〉态之间温度依赖和H0依赖的能级间隔;kB为玻尔兹曼常数.因此,〈µEr〉随温度的下降而单调地增加.此外,通过稀释效应(调节比例因子x),非磁性的Y3+掺杂将降低Er3+的自旋密度.结合这两种机制,所得的净Er3+自旋密度随温度、磁场、掺杂浓度的关系写为

对于图7(a)—(k)所示的光谱曲线,提取出所有免交叉谱线中可分辨的Λ值.在图7(l)给出了Λ与计算得到的的一一对应关系.如图7(l)所示,所有的数据点落在一条穿过原点的直线上.这一比例行为证实了 Er3+自旋关联地与 F e3+的真空磁子场耦合.

为了从实验数据中获得更多的微观信息,Li等[46]发展了一种微观理论模型,对于 F e3+子系统,遵循Herrmann的处理方法.哈密顿量HˆFe考虑两个子晶格自旋S2l-1和S2l的对称交换相互作用、反对称交换相互作用和磁各向异性,其中晶胞指数l=1,2…Nuc(Nuc是晶胞的总数).对于Er3+子系统,可以以第l个晶胞内 Er3+自旋Rl的拉莫进动来模拟跃迁.EPR的哈密顿量其中磁矩表示为各向异性朗德因子g和R1的点积,是真空磁导率.F e3+-Er3+的耦合哈密顿量写成

其中J1和J2是对称交换常数,D1和D2为反对称交换常数.由总的哈密顿量得到

由于对称参数J1=J2=J,γ是旋磁比,JFe是Fe3+的自旋之间各向同性的交换常数,ℏ为普朗克常数除以2π.通过(6)式拟合 Λ 的斜率得到交换常数J=2.95 meV.

6 THz磁子-极化子

最近,莫斯科科技大学的Grishunin、上海大学马国宏以及内梅亨大学的Kimel等人的合作研究表明,稀土正铁氧体中的时间分辨THz磁共振光谱必须考虑电磁波在材料中的传播效应[47].以THz波在TmFeO3中的传播为例,实验发现当THz波的频率与qAFM模式的自旋共振频率接近时,磁子-极化子(magnon-polaritons)扮演了一个重要的角色.

典型数据如图8(a)和(b)所示,TmFeO3的温度为40 K,图中的振荡曲线为THz透射脉冲及其回波之间的THz自由感应衰减信号的时域和频域信号.红色圆圈为实验结果,黑色实线为理论计算结果.图8(c)和(d)显示的是TmFeO3样品薄片在激光辐照下发射的THz脉冲电场的时域谱和频谱.THz透射和发射实验结果都表明,样品在40 K的自旋动力学表现为略高于和略低于qAFM共振频率两个模式的拍频信号.通过建立模型描述TmFeO3的自旋与THz脉冲的相互作用,理论结果表明这两个频率为磁极化子的上下边带,如图8(e)所示.

图8 (a)透过TmFeO3晶体的THz时域波形及其(b)傅里叶变换谱;(c)TmFeO3薄片的THz产生波形及其相应的(d)傅里叶变换谱;(e)自旋共振附近的磁子-极化子色散关系[47]Fig.8.(a)THz waveforms transmitted through the TmFeO3 sample and(b)its Fourier transforms of the time traces;(c)THz generation in a TmFeO3 slab by a laser pulse and(d)its Fourier spectra shown in(c);(e)the magnon-polariton dispersion in the vicinity of the spin resonance[47] .Reproduced with permission from Ref.[47].

当THz光斑尺寸大于THz波长,利用平面波近似,设电磁波的波矢k垂直于界面.THz脉冲透过一个厚度为d的样品薄层,其相对介电常数和磁导率分别是ε和µ(ω).在频谱上的耦合波方程写为

式中Ex为电场时域信号的傅里叶变换,波矢最终求得透过样品的光谱和入射样品的光谱之比为

式中k0是空气中的波矢,k1为材料中的波矢,d是样品的厚度,E(ω)为入射THz脉冲的归一化光谱振幅.对于材料的磁导率,用标准Drude-Lorentz形式表示为

其中ω0为共振频率,Δω为吸收线宽,Δµ为有效振荡强度.当ω0=0.86THz,d=1.5mm,其他参数为Δω=8.6GHz,Δµ=10-3,ε=22-0.26i,并将入射的THz脉冲的光谱取逆傅里叶变换代入(9)式中,可以得到透过样品的THz脉冲的时域和频域波形,实验结果表明理论模型与实验结果符合得很好.

由于磁极化子在空间的传播效应,有必要研究拍频信号和晶体厚度的依赖关系.实验结果表明,当晶体的厚度越厚,拍频信号变得更加显著.对于THz透射光谱,当晶体的厚度小于1 mm,极化子的拍频和THz脉冲的Fabry-Perot回声场发生干涉,这将阻碍观察到明显的拍频信号.

7 总 结

随着THz光谱技术的不断发展,基于RFeO3晶体的交换相互作用的超快光学调控,对磁学的基础研究具有重要的意义.本文主要介绍了THz脉冲的磁场分量实现反铁磁晶体中两类自旋模式的“非热”激发、相干控制、实现自旋重取向的可视化的研究进展.特别介绍了最近发现固体中稀土铒离子自旋系综与铁离子磁振子间的协同耦合效应表现为真空拉比劈裂特性.Fe3+和Er3+耦合强度与稀土离子Er3+的浓度依赖关系满足迪克模型.由此得出的Fe3+和Er3+交换耦合常数为人们更加深刻地理解3d-4f电子的磁耦合导致的新奇效应提供了重要的实验和理论证据.此外,介绍了磁子-极化子(magnon-polaritons)对于 THz波在 TmFeO3中的传播具有重要的影响.这两项工作为利用量子光学中的概念去理解和控制凝聚态物质中的新物相提供了新的途径.