光束正入射至界面时的自旋-轨道相互作用及其增强*

罗慧玲 凌晓辉 † 周新星 罗海陆

1) (衡阳师范学院物理与电子工程学院, 智能信息处理与应用湖南省重点实验室, 衡阳 421002)

2) (湖南师范大学物理与电子科学学院, 长沙 410081)

3) (湖南大学物理与微电子科学学院, 长沙 410082)

光束正入射至均匀突变界面时的自旋−轨道相互作用表现为拓扑荷数为±2的、自旋可控的涡旋相位.然而, 该涡旋相位的物理来源以及界面的性质在自旋−轨道相互作用过程中起到何种作用, 这些问题还有待解决.首先建立一个简洁的菲涅耳琼斯矩阵来描述这种自旋−轨道相互作用, 并揭示其中的涡旋相位其实是一种贝里(Berry)几何相位, 它来源于光束本身的拓扑结构, 而界面的性质影响自旋−轨道相互作用的转换效率.一般情况下, 转换效率极低, 限制了其应用.因此, 基于上述理论, 提出采用光轴平行于界面法线方向的单轴薄层材料, 来极大地增强这种自旋−轨道相互作用.

1 引 言

光既可以具有自旋角动量, 又可携带轨道角动量.自旋角动量与光的偏振有关, 如左右旋圆偏振光子分别携带的自旋角动量.轨道角动量有两类[1,2], 一类为内禀(intrinsic)的轨道角动量, 与涡旋光场有关, 每个光子携带的轨道角动量, 其中l为涡旋相位的拓扑荷数; 另一类为外禀(extrinsic)的轨道角动量, 和光束传播的轨迹有关, 定义为坐标原点到光束中心的距离与线动量的叉乘, 与经典粒子的机械角动量类似.光的自旋角动量和轨道角动量之间的相互转换和耦合被称为自旋−轨道相互作用 (spin−orbit interaction, SOI)或耦合[1,2].它是光学中的一种基本效应, 广泛存在于界面的反射和折射、非均匀各向异性介质、强聚焦、粒子散射、表面波和消逝波等体系中, 在光学、纳米光子学和等离子光学等领域扮演越来越重要的角色, 并在精密测量与探测、信息存储与处理、微粒操纵以及各种功能光子器件设计等方面显示出巨大的应用潜力[1−9].在旋转对称的系统中, 光的SOI表现为自旋可控的涡旋相位的产生(内禀轨道角动量); 在旋转对称性破缺的系统中, 它表现为自旋霍尔效应(外禀轨道角动量)[1,2].

光的自旋霍尔效应存在于很多体系中, 如光束在突变界面的斜入射[1,10−17]、一维的潘查拉特南−贝里 (Pancharatnam−Berry, PB)相位元件[18−21]、各向同性的非均匀材料[21−23]等; 自旋可控的涡旋相位也在方位变化的PB相位元件[24−27]、强聚焦[28,29]、单轴晶体中的传输[30,31]等体系中出现.然而, 有趣的是, 当光束正入射至均匀的、各向同性的突变界面时, 也能产生自旋相关的涡旋相位[32−34].光束正入射时, 极小的一部分入射光束发生自旋反转(左旋变为右旋或者右旋变为左旋), 并获得拓扑荷数为±2的涡旋相位(图1(a)).其内在机制被认为是SOI, 但这种相位的物理来源、为什么拓扑荷数为±2以及界面在其中究竟扮演何种角色等一系列的问题, 目前并不清楚.另外, 该SOI与光束通过方位变化的各向异性PB相位元件[24−27]时产生涡旋相位的过程极为相似.光束入射到方位PB相位元件时, 一部分入射光束发生自旋反转并获得2倍于元件拓扑荷数(q)的涡旋相位因子2qφ, 其中φ是PB相位元件的局部的光轴方向, 是坐标位置的函数.也就是说这种相位因子来源于PB相位元件的非均匀的各向异性.而前文所提到的界面是各向同性且均匀的, 这与PB相位元件的情况又有何联系和区别?

首先建立一个菲涅耳琼斯矩阵来描述光束正入射至突变界面的透(反)射光束, 并发现其产生的涡旋相位因子来源于光束本身的拓扑结构, 具有几何性, 是一种自旋重构的贝里(spin−redirection Berry)相位.界面的性质影响光束中各平面波分量的菲涅耳系数, 并决定发生了SOI的那部分光束的转换效率.而对于PB相位元件, 光束也是部分发生自旋反转, 经历SOI, 并获得PB涡旋相位, 但这种相位来源于外部材料的各向异性.进一步研究发现, 光束正入射至突变界面的SOI的转换效率取决于光束中各斜入射的平面波的TM和TE分量的菲涅耳系数之差.对于传统材料来说, 这种效应非常弱, 转换效率极低, 这限制了它的应用, 目前也没有这方面的实验见诸报道.Ciattoni等[33]从理论上提出, 采用各向同性的、介电常数近零的薄层来增强这种效应, 但因各向同性材料可调的自由度非常有限, 其转换效率最高也只能达到20%左右.由于各向异性材料具有更多调控的自由度,因此提出用光轴方向平行于界面法线方向的单轴晶体薄层, 来极大地增强正入射时的SOI, 使转换效率在某些条件下可达100%.

图1 光束正入射至各向同性的突变界面时SOI的示意图 (a) 左旋圆偏振光束正入射至界面后, 部分光束发生自旋反转变成右旋光, 并获得拓扑荷数为2的涡旋相位(两个小图分别表示一种典型的涡旋光束的强度和相位分布); 注意, 未发生SOI的那部分光束并没有在图中画出; 和 分别表示左、右旋圆偏振; (b) 光束中各平面波分量的自旋与局部坐标的旋转耦合的示意图, 其中圆锥代表光束的角谱, 绿色的箭头线代表任意的两支平面波的波矢, 橙色带箭头的小圆圈表示各平面的偏振矢量在实验室坐标上的投影(均为圆偏振), Ω ξ 为坐标旋转的空间旋转Fig.1.Schematic illustration of the SOI for a light beam normally impinging onto a sharp isotropic interface.(a) When a left−circu−larly polarized beam normally passes through the interface, part of the incident beam converts into a right−circularly polarized beam, and carries a vortex phase with a topological charge of 2.Note that the spin−maintained portion is not shown in the picture.and denotes the left− and right−handed polarization, respectively.(b) Schematic illustration of rotational coupling between the local coordinates and the spin of the plane wave components within the beam spectra.The cone represents the angular spec−trum of the beam.The two green arrows represent the wave vectors of arbitrary two plane waves.The orange circles with arrows indicate the projection of polarization vectors of each plane wave on the laboratory coordinates (all circularly polarized).Ω ξ is the spatial coordinate rotation.

2 涡旋相位的物理来源：理论与模型

2.1 建立透射与入射光束之间的菲涅耳琼斯矩阵

考虑一个单色的有限宽傍轴光束, 正入射至一个由各向同性的、均匀的、无吸收材料构成的突变界面, 如图1(a)所示.由于反射和透射光束具有相似的行为, 因此本文以透射为例来分析其中的SOI过程.由角谱理论可知, 有限宽的光束可以看成由许多具有略微不同传播方向的平面波相干叠加而成.众所周知, 平面波在界面的透射可以由菲涅耳公式来描述, 因此, 光束在界面的透射场由其所有平面波分量的透射场相加(积分)而成.根据角谱理论, 可以把入射(i)和透射(t)电场统一地写成如下傅里叶积分形式(a = i, t)[35]：

为了使表达式看起来简洁, 以下将 tTM,TE(ϑi) 简写为 tTM,TE.由于透射光和入射光的观察面上的场均是指垂直于中心波矢的横向场分布, 因此将非中心平面波的偏振矢量投影到中心平面波的偏振矢量上, 并忽略纵向z分量得

将上式写成矩阵形式, 并计算得

联立(2)和(4)式, 可得

通过矩阵相乘得：

其中 ζ =cosϑt/cosϑi.至此, 在圆偏振基下, 用菲涅耳琼斯矩阵来建立之间的关系.很显然, (6)式的矩阵中, 反对角元是自旋反转的项, 而对角元是自旋不变的项, 可分别称之为反常(abnormal)模式和寻常(normal)模式.这种矩阵不仅比已有的方法[32,33]简洁, 而且清晰地展现了光束本身拓扑结构(矩阵)与界面性质(矩阵)的不同的物理贡献.也就是说, 光束本身的拓扑结构贡献了涡旋相位因子 e±i2φ; 而与界面性质有关的菲涅耳系数贡献了反常模式和寻常模式的振幅 ( ζtTM± tTE)/2 , 并决定了透射光束中反常模式和寻常模式的“比重”(即SOI中的转换效率).在左旋圆偏振光束(其圆偏基下的角谱分布为入射下, 根据(6)式可得

将(7)式代入(1)式即可得到透射光束的电场分布.

2.2 Berry相位: 涡旋相位的物理来源

下面分析涡旋相位因子的物理来源.很显然,(7)式中的反常模式携带一个拓扑荷数为2的涡旋相位.由(4)式可知, 它来源于非中心平面波与中心平面波的投影操作.任意非中心平面波的偏振矢量投影到中心平面波的偏振矢量上后, 产生了一个自旋相关的涡旋相位因子 e±iφ.这种相位因子与各平面波的入射面的方位角有关, 而中心平面波的入射面无法确定, 即TM和TE分量的偏振矢量也无法确定, 因此是涡旋相位的奇点.对于旋转不变的光束(如高斯光束和贝塞尔光束), 该相位因子都存在.从本质上看, 该相位因子来源于光束本身的拓扑结构, 是几何性的, 本身是一个不可观测的量,只有透射和入射光束的相位差才是可观测量.对于反常模式, 由于发生了自旋反转, 相位差为 ± 2φ ;对于寻常模式, 相互抵消, 相位差为0.

本文用自旋角动量与坐标旋转之间的耦合[1]或者光学科里奥利(Coriolis)效应[36]来解释这种相位的物理来源.由于要求入射光束在横截面上是均匀的圆偏振(实验中容易产生), 即每个非中心平面波的偏振矢量在投影到中心平面波的横向面之后都是圆偏振的(图1(b)), 因此光束在z方向上的平均光子自旋角动量为 J =σ(其中σ = +1和—1分别表示左旋和右旋圆偏振).以每个非中心平面波各自的TM和TE偏振矢量构成的局部坐标框架, 投影到实验室坐标上后, 相对于实验室坐标的空间旋转率为 Ωξ=dφ/dξ , 其中 ξ 是坐标框架的旋转路径, 旋转轴为z轴.因此, 推导出一个几何相位[1,36]：

它体现为自旋( σa)与坐标旋转( φ )之间的耦合,只与坐标旋转的路径有关, 因此是几何性的.而φ又是螺旋相位因子, 与内禀轨道角动量有关, 因此又体现为自旋与内禀轨道角动量之间的耦合,即SOI.这里的坐标旋转, 是指各非中心平面波所在的局部坐标系在实验室坐标上的投影, 相对于实验室坐标系的旋转.这种几何相位与光束中各平面波的传播方向的SO(3)旋转有关, 因此它是自旋重构的贝里相位[1,36−38].

实际上, 最终透射光束的寻常和反常模式的相位是透射光束的几何相位与入射光束的几何相位之差由于反常模式是自旋反转的结果,即 σt=-σi, 因此反常模式的几何相位为而寻常模式的几何相位为0.这个结果与(7)式的计算结果是一致的, 也与文献[32−34]中的结果相同.这种相位因子还与PB相位元件中产生的涡旋相位在形式上极为相似.在PB相位元件中, PB相位来源于材料外部的各向异性; 而这里的几何相位来源于光束本身的拓扑结构以及光束中各平面波分量由于斜入射造成的 tTM与 tTE之间的不同(也可看作是一种“各向异性”).这两种情况在原理上是不同的, 但在形式上又是一致的, 可以用上述的自旋角动量与坐标旋转的耦合模型来统一地理解[1,36].

还可以从角动量守恒的角度来考虑.由于本文中研究的界面是关于z轴旋转对称的、无吸收的体系, 因此参与SOI的那部分光束在z方向上的总角动量必须守恒(诺特定理).入射光束中各平面波分量投影到z方向上的自旋角动量(光子的平均自旋角动量)为 σ ℏ , 且不携带轨道角动量.透射后, 部分光束发生自旋反转, 其自旋角动量变为 - σℏ , 且同时获得与入射自旋相关的、 2 σℏ 的额外轨道角动量.此时这部分光束的总的角动量还是 σ ℏ , 并无增减.因此, z方向的总角动量是守恒的.

3 增强光束正入射时的SOI效率的方法

3.1 现有方法的SOI效率

上文建立了一个由各向同性的、均匀的、无吸收材料构成的界面的菲涅耳琼斯矩阵, 分析了光束本身的拓扑结构和构成界面的材料性质的各自贡献, 即涡旋相位来源光束本身的拓扑结构, 而界面性质影响SOI的转换效率.现以左旋圆偏振贝塞尔光束的正入射为例来具体讨论.零阶贝塞尔光束的横向的电场可写为

其中A0为任意振幅, Δ k=2π/w0为横向的谱半宽度, w0为光束束腰半宽度,表示第一类n阶贝塞尔函数.其角谱分布为一个冲激函数的形式：

这意味着, 贝塞尔光束的角谱实际上呈旋转不变的、空心圆锥状分布(图1(b)), 中心轴垂直于界面,且所有平面波分量均为斜入射, 虽方位角不同, 但入射角均为 ϑi=sin-1(Δk/ki) , 即 ϑi由 w0决定.由此, 将代入(7)式并联立(1)式, 可得透射光束的电场为

式中透射光束分为两部分, 一部分与入射光束相同, 是寻常模式; 另一部分表现出自旋反转现象并携带拓扑荷数为2的涡旋相位, 是反常模式.图2给出了寻常模式和反常模式的光强和相位分布.反常模式(图2(a))光斑中心是光强为0的空心区域,相位在方位方向变化4π, 即拓扑荷数为2的涡旋相位; 寻常模式(图2(b)) 光斑中心是实心区域, 不携带方位方向的涡旋相位.当寻常模式强度为0,即只有反常模式时, 该SOI过程中的转换效率为100%.然而一般情况下, 对于传统材料构成的界面, | tTM-tTE| 实际上是一个非常小的值, 尤其在入射角较小时.这种情况意味着SOI的转换效率是极低的, 常规的实验精度难以被观察到.这也是这种效应迄今为止没有在实验上被观察到的原因之一.

SOI中的转换效率可定义为透射光束的反常模式的功率与入射光束的功率之比.当考虑贝塞尔光束正入射时, SOI的转换效率为[33]

因此, 转换效率取决于贝塞尔光束中各平面波的TM和TE分量菲涅耳系数之差.对傍轴光束( Δ k≪k )来讲, ζ =cosϑt/cosϑi是一个接近于1的值.若考虑入射介质和出射介质折射率相等, 如一个放置于自由空间中的薄层, ζ ≡1 .此时,SOI的转换效率 η =|tTM-tTE|2/4 .由于贝塞尔光束中的平面波分量的入射角 ϑi取决于光束束腰半宽度w0的大小, 且具有一个确定的值, 因此, 只要找到合适的材料, 使 | tTM-tTE| 在感兴趣的 ϑi范围内具有尽可能大的值, 就能获得尽可能增强的SOI.实际理想情况下, 如果能够使 tTM和 tTE的值一个为1, 另一个为—1, 则能获得100%的效率.

图2 左旋圆偏振贝塞尔光束正入射至一个界面时, 透射光束的反常模式(a)和寻常模式(b)的归一化光强分布, 其中两个小图分别表示为对应的相位分布, 在计算中, 取入射光束的波长 λ =1 且w0=20λFig.2.Normalized intensity distribution of the abnormal mode (a) and normal mode (b) of transmitted light beam under the nor−mal incidence of a left−handed circularly polarized Bessel beam at a sharp interface.The insets represent the phase distribution of corresponding modes.Here, we take the working wavelength as λ =1 and w0=20λ .

图3 三种放置于自由空间的单层薄膜材料的透射系数, 以及SOI的转换效率 η =|tTM-tTE|2/4 (a) ε =2.25 , (b) ε =0.01 ,(c) εx= εy=1 且 εz=0.01 ; 计算中, 取入射光束的波长 λ =1 , 三种材料厚度 h =2λ ; (d)和(e)分别是 ϑi 和 εz , ϑi 和h同时变化时的转换效率, 在(d)中, 取 εx= εy=1 , h =1λ ; 在(e)中, 取 εx= εy=1 ,εz=0.01Fig.3.Transmission coefficients and conversion efficiency ( η =|tTM-tTE|2/4 ) of three optically thin films placed in free space：(a) ε =2.25 , (b) ε =0.01 , (c) εx= εy=1 and εz=0.01 , where we take λ =1 and h =2λ ; (d) conversion efficiencies versus ϑi and εz of a uniaxial layer with εx= εy=1 and h =1λ ; (e) conversion efficiencies versus ϑi and h of a uniaxial layer with εx= εy=1and εz=0.01 .

然而, 传统的材料(如空气、玻璃等)构成的界面, SOI的转换效率极低.考虑一个放置于自由空间、介电常数为ε = 2.25、厚度为h (波长量级)的非磁性(磁导率为1)各向同性薄膜, 其TM和TE平面波分量的透射系数分别为[39]

鉴于传统材料的SOI极弱, Ciattoni等[33]在2017年从理论上提出采用介电常数近零的各向同性薄层来增强这种效应(图3(b)).对于介电常数近零材料薄层, TM波在入射角很小时就可以满足法布里−珀罗共振, 使 | tTM| 达到1, 并同时使|tTM-tTE|达到较大的值.然而, 其转换效率最高也只可达20%左右.其原因是, 对于各向同性材料,tTM和 tTE同时受到介电常数的影响, 无法独立地调控, 很难使 | tTM-tTE| 接近于2, 因此转换效率难以达到100%.

3.2 通过单轴薄层增强SOI的效率至100%

各向异性材料比各向同性材料具有更多的自由度, 有望获得接近100%的效率.考虑厚度为h的非磁性(磁导率为1)单轴晶体薄层, 其介电常数张量为

即此单轴层的光轴与z轴平行, 也就是光束的传输方向.这种情况下, 体系依然具有旋转不变性, 上文的所有理论仍然适用.由于该材料具有 εx,y和εz两个可以独立调控的介电常数, 所以相比于各向同性材料多了一个调控的自由度.写出其透射系数为[40]

注意, 其中TM和TE波的交叉偏振透射系数

当 εx= εy= εz时, (13)式回到各向同性材料的情况, 即(11)式.

要想获得100%的转换效率, 必定要使 tTM和tTE的值一个为1, 另一个为—1, 当然也就没有反射.首先, 令 εx= εy=1 , 即平行于界面方向的介电常数与自由空间相同, 这保证了TE波的透射系数为tTE=exp(iqoh), 其模值 | tTE|≡1 (图3(c)), 即任意入射角下, TE波均全部透射.然后, 令 εz→ 0 (远小于自由空间的介电常数), 使TM波的全反射临界角(即满足 εz=sin2ϑi条件时)变得很小.大于临界角时, tTM=0 , | tTM-tTE|=1 , 转换效率恒为25%; 小于临界角时, 由于TM波可能满足法布里−珀罗共振条件 2 qeh=2mπ (m为整数), 而出现全透射的情况( tTM=±1 ), 使得 | tTM-tTE| 的值可能为2, 也可能为0, 以及它们之间任意的中间值(两个共振峰之间).也即, 在TM波满足法布里−珀罗共振时, 有可能实现100%的效率.

下 面 具 体 计 算 εx= εy=1 且 εz=0.01 时 , 厚度h = 2λ的单轴薄层的透射系数和转换效率.如图3(c)所示, | tTE| 在任意入射角度时恒为1; 而TM波在 ϑi＞ 5.8°时发生全反射( | tTM|=0 ), 当ϑi＜5.8°时, 在某些角度发生法布里−珀罗共振, 出现tTM=±1的情况, 使SOI的转换效率η=|tTM-tTE|2/4在 ϑi=3.8°和 5 .6°附 近 达 到 100%.当 εz和 厚 度h发生改变时, 影响法布里−珀罗共振出现的角度,因此出现100%效率的角度 ϑi也随之改变, 如图3(d)和图3(e)所示.

虽然具有这种极端的介电常数的单轴材料很难在自然界中找到, 但近年来随着微纳光学, 特别是超构材料和超构表面领域的飞速发展, 具有上述等效介电常数的人工合成材料有望通过合适地设计超构材料或超构表面来实现, 比如双曲超构材料(hyperbolic metamaterials)[41,42].最后还需要指出的是, 本文所建立的全波理论适合于旋转不变的体系, 比如由各向同性材料构成的界面、各向同性材料与光轴平行于界面法线方向的单轴晶体构成的界面.当系统不具有旋转不变性时, 比如光轴方向平行于界面时的单轴晶体, 需对本文理论作较大修正才能适用.

4 结 论

首先建立了菲涅耳琼斯矩阵来描述光束正入射至突变界面时的SOI, 分析和澄清了光束本身的拓扑结构和界面性质在SOI中所扮演的不同角色,并揭示了其中所产生的涡旋相位的物理根源.简而言之, 这种涡旋来源于光束本身的几何拓扑结构,在本质上是一种自旋重构的贝里相位; 而界面的性质影响SOI的转换效率.在形式上, 这种SOI与PB相位元件中的SOI极为相似, 但PB相位来源于材料外部的各向异性.由于该效应在一般情况下极弱, 限制了其应用.因此提出用光轴平行于界面法线方向的单轴晶体薄层来有效地增强它, 使之在一定条件下达到100%的转换效率.本文的研究不但为这种SOI建立了简洁明晰的理论框架, 而且揭示了现象背后的物理机理, 并进一步给出了增强这种效应的可行方案, 为未来的潜在应用指明了方向.