超表面 PB 相位原理：庞加莱球上的几何相位

引言：三个环环相扣的问题

超表面里有一类奇妙的相位：不改尺寸、不改材料，只把纳米柱原地旋转一个角度 θ，出射的圆偏振光就会获得一个相位，而且这个相位恰好是旋转角的两倍：

$$\varphi = 2\sigma\theta$$

其中 $\sigma = \pm 1$ 对应两种圆偏振。这就是 PB 相位（Pancharatnam–Berry 相位）。

这个公式看似简单，背后却藏着三个环环相扣的问题：

1. 旋转角凭什么变成相位？ 旋转是一个几何动作，相位是波的振荡参量，两者之间没有天然的等价关系；
2. 为什么恰好是两倍（2θ），而不是 θ 或者别的倍数？
3. 这个相位为什么与波长无关？

本文用两条路线回答前两个问题：先走琼斯矩阵的代数路线（最直接、可逐行验算），再走庞加莱球的几何路线（给出图像直觉）；第三个问题（无色散）在文末回答。两条路线不是平行的两个答案，而是互相印证：代数告诉我们”发生了什么”，几何告诉我们”为什么必然如此”。

一、琼斯矩阵路线：旋转波片如何”变出”相位

1.1 一根纳米柱 = 一个微型波片

一根各向异性的纳米柱有两个正交主轴（快轴与慢轴），光沿两个主轴传播时折射率不同，因此两个分量的相位延迟不同。它本质上就是一个微型波片：相位延迟为 $\delta$。

设快轴沿 $x$ 方向，波片对电场 $(E_x, E_y)$ 的作用可用琼斯矩阵表示：

$$J_0 = \begin{pmatrix} e^{i\delta/2} & 0 \cr 0 & e^{-i\delta/2} \end{pmatrix}$$

含义很直白：$x$ 分量（沿快轴）获得相位 $+\delta/2$，$y$ 分量（沿慢轴）获得相位 $-\delta/2$，两者相位差正好是 $\delta$。当 $\delta = \pi$ 时就是半波片：两个正交分量相差半个波长，于是线偏振方向被旋转、圆偏振旋向被翻转。

先记住一个关键点：如果波片一直不转（$\theta=0$），这个矩阵不会产生任何与旋转有关的相位——它只是分别给两个分量加相位。旋转才是把”角度”变成”相位”的开关。接下来就看这个开关如何起作用。

1.2 旋转 θ 角：J(θ) = R(θ)J_0R(−θ)

现在把波片绕光轴旋转 θ，使快轴指向 $\theta$ 方向。计算分两步：先把电场旋转到波片自己的坐标系（乘旋转矩阵 $R(-\theta)$，相当于问”这个电场在快轴、慢轴上各有多少分量”），在波片坐标系里应用 $J_0$，再旋转回实验室坐标系（乘 $R(\theta)$）。矩阵乘积从右往左读，先作用的是 $R(-\theta)$：

$$J(\theta) = R(\theta) J_0 R(-\theta), \qquad R(\theta) = \begin{pmatrix} \cos\theta & -\sin\theta \cr \sin\theta & \cos\theta \end{pmatrix}$$

把矩阵乘开（这一步值得自己验算一遍），利用二倍角公式 $\cos^2\theta = \frac{1+\cos2\theta}{2}$、$\sin^2\theta = \frac{1-\cos2\theta}{2}$、$\sin\theta\cos\theta = \frac{\sin2\theta}{2}$，得到：

$$J(\theta) = \cos\frac{\delta}{2} \begin{pmatrix} 1 & 0 \cr 0 & 1 \end{pmatrix} + i\sin\frac{\delta}{2} \begin{pmatrix} \cos2\theta & \sin2\theta \cr \sin2\theta & -\cos2\theta \end{pmatrix}$$

注意这里发生了什么：旋转角 θ 只以 $2\theta$ 的形式出现。这是”两倍”的第一个来源——旋转矩阵的元素本身含 $\cos\theta\sin\theta$、$\cos^2\theta$ 等，二倍角公式把它们组合成了 $\cos2\theta$、$\sin2\theta$。物理上更本质的原因是：波片有”二重对称性”，快轴转 $180°$ 后元件与原来不可区分，所以一切响应都按 $2\theta$ 周期变化。

1.3 换到圆偏振基：相位因子 e^{±2iθ} 现身

上面的矩阵是在线偏振基 $(E_x, E_y)$ 下写的。PB 相位关心的是圆偏振，所以换成圆偏振基。右旋圆偏振 $\sigma=+1$ 与左旋圆偏振 $\sigma=-1$ 的基矢（在 $x,y$ 线偏振基下的坐标）为：

$$|\sigma_+\rangle = \frac{1}{\sqrt2}\begin{pmatrix}1 \cr i\end{pmatrix}, \qquad |\sigma_-\rangle = \frac{1}{\sqrt2}\begin{pmatrix}1 \cr -i\end{pmatrix}$$

这里需要先立下一个旋向约定，因为”哪个是右旋”在不同文献里可能互换：本文采用物理光学常用的约定——时间因子取 $e^{-i\omega t}$，从接收端（面向光源）看电场矢量顺时针旋转为右旋。在此约定下，右旋圆偏振的琼斯矢量正是 $(1,i)/\sqrt2$（即 $E_y = iE_x$，两个分量等幅、相位差 $90°$），左旋为 $(1,-i)/\sqrt2$。若改用 $e^{+i\omega t}$ 或从发射端看的约定，$(1,i)$ 与 $(1,-i)$ 的命名会互换，但下面推导的物理结论不变。

这一步的目的是让旋转角从”角度”变成”相位”：在 1.2 的 $J(\theta)$ 里，$\theta$ 还躲在 $\cos2\theta$、$\sin2\theta$ 这样的三角函数里，它是个角度，不是相位；而 PB 相位关心的是圆偏振，圆偏振天然用指数 $e^{i\varphi}$ 描述。把 $J(\theta)$ 换到圆偏振基后，非对角元变成 $i\sin\frac{\delta}{2}e^{\pm 2i\theta}$——$\theta$ 从三角函数里”钻”进了指数，直接以相位因子 $e^{\pm 2i\theta}$ 的形式出现：

$$J_{\text{circ}}(\theta) = \begin{pmatrix} \cos\frac{\delta}{2} & i\sin\frac{\delta}{2} e^{-2i\theta} \cr i\sin\frac{\delta}{2} e^{2i\theta} & \cos\frac{\delta}{2} \end{pmatrix}$$

非对角元描述的是右旋与左旋之间的相互转换，所以”旋转角变成相位”这件事，恰好发生在圆偏振互相转换的那个分量上——这正是 PB 相位的代数来源。不过这一节只是”看出”了相位因子的存在，下一节回到线偏振基，用右旋圆偏振的琼斯矢量 $(1,i)/\sqrt2$ 代入算一遍，验证它确实落到了输出光上。

1.4 半波片极限：读出 φ = 2σθ

当纳米柱满足半波片条件 $\delta = \pi$ 时，$\cos\frac{\pi}{2}=0$、$\sin\frac{\pi}{2}=1$，1.2 的展开式退化为：

$$J(\theta) = i \begin{pmatrix} \cos2\theta & \sin2\theta \cr \sin2\theta & -\cos2\theta \end{pmatrix}$$

现在把右旋圆偏振的琼斯矢量 $(1,i)/\sqrt2$ 直接代入：

$$\frac{i}{\sqrt2} \begin{pmatrix} \cos2\theta & \sin2\theta \cr \sin2\theta & -\cos2\theta \end{pmatrix} \begin{pmatrix} 1 \cr i \end{pmatrix} = \frac{i}{\sqrt2} \begin{pmatrix} \cos2\theta + i\sin2\theta \cr \sin2\theta - i\cos2\theta \end{pmatrix} = i e^{2i\theta} \cdot \frac{1}{\sqrt2} \begin{pmatrix} 1 \cr -i \end{pmatrix}$$

输出是 $(1,-i)/\sqrt2$，即左旋圆偏振（自旋翻转），并且带相位因子 $e^{2i\theta}$（常数 $i$ 是全局相位，所有输出都有，可吸收掉）。同理，把左旋琼斯矢量 $(1,-i)/\sqrt2$ 代入，会得到右旋 $(1,i)/\sqrt2$ 并带 $e^{-2i\theta}$。合起来就是：

$$\varphi = 2\sigma\theta$$

到此，第一个问题的代数答案齐了：旋转角通过波片琼斯矩阵的二倍角结构，在圆偏振基下变成 $e^{\pm 2i\theta}$ 的相位因子，半波片条件下全部转换为自旋翻转分量，于是得到 $2\sigma\theta$。

但琼斯矩阵只是”算出”了 $2\theta$，并没有解释”为什么是 $2\theta$”——它忠实描述了结果，却没有告诉我们旋转角为什么以这种几何的方式进入相位。要看清原因，需要换一个视角：把偏振态本身当作球面上的几何对象。这就是下一节的庞加莱球。

二、为什么相位”不可传递”：Pancharatnam 连接

PB 相位还有一层更深的身份：它是几何相位。要理解这一点，必须回答一个看起来更基础的问题——两束不同偏振的光，凭什么比较相位？

2.1 什么算”同相”：内积的辐角

设两束光偏振态分别为 $|A\rangle$、$|B\rangle$（归一化）。叠加后的强度：

$$I = \left| |A\rangle + |B\rangle \right|^2 = 2 + 2\left|\langle A|B\rangle\right|\cos\bigl(\arg\langle A|B\rangle\bigr)$$

当 $\cos(\arg\langle A|B\rangle)=1$ 时干涉最强，我们说这两束光”同相”。于是自然地定义：$A$、$B$ 之间的相位差就是内积 $\langle A|B\rangle$ 的辐角。这就是 Pancharatnam 在 1956 年给出的连接（connection）。

这个定义有一个惊人的后果：“同相”关系不可传递。如果 $B$ 与 $A$ 同相、$C$ 与 $B$ 同相，那么 $C$ 与 $A$ 却不一定同相。

这乍听起来很反直觉——“同相”听上去像一种等价关系，应该可以传递。关键在于：这里的”同相”不是两个态之间绝对的相等，而是沿着某条路径建立的连接；路径不同，结果就可能不同。下面的球面图像会把这个”反直觉”变成”理所当然”。

2.2 几何原因：球面上的平行输运

为什么？先做一个地球仪实验，但这次不要用指南针——指南针始终指向地球上的磁北极，绕任何一圈回来它还是指北，不会偏转。我们要追踪的是另一种东西：你的朝向。

从北极出发，面朝正南，沿一条经线（大圆弧）走到赤道；然后沿赤道向东走 $90°$ 经度；再沿另一条经线面朝正北走回北极。三段都是大圆弧，恰好围成一个闭合的测地三角形。回到北极后你发现：你的朝向与出发时相差了 $90°$——出发时你面对的是经线 A 的方向，回来时面对的是经线 B 的方向，两条经线在北极的夹角正好是 $90°$。

这个 $90°$ 不是走路走错的误差，它等于三角形在球面上围住的立体角 Ω：赤道上相距 $90°$ 的两点与北极围成的三角形，恰好占球面的八分之一，立体角正是 $4\pi/8 = \pi/2$（$90°$）。球面是弯的，所以”保持朝向不变地沿大圆走”这件事本身就会累积一个转角——这正是平行输运的几何后果。

偏振态的情形完全一样：把三个偏振态画在庞加莱球上，它们围成一个测地三角形，沿球面平行输运（每一步都保持”同相”地前进）走完一圈，回到起点时，偏振态相对初态多了一个相位。你偏转的是朝向，偏振态多出来的是相位——两者都来自球面的弯曲。

图 1：庞加莱球上的测地三角形与平行输运。三个偏振态 A、B、C 围成测地三角形；沿闭合路径平行输运一周，方向旋转量等于立体角 Ω，旋量对应相位 −Ω/2。

偏振态不是普通向量，而是旋量（spinor）：庞加莱球上的一个点对应一个旋量态。旋量有个特性——方向旋转 $360°$ 后它变成自己的相反数，旋转 $720°$ 才复原；所以方向旋转 $\Omega$ 对应相位因子 $e^{\pm i\Omega/2}$，而不是 $e^{\pm i\Omega}$。于是沿闭合路径走一圈，旋量相位变化：

$$\Delta\varphi = -\frac{\Omega}{2}$$

把三条边逐段”同相”地连起来：$A$ 到 $B$ 同相、$B$ 到 $C$ 同相，沿 $C$ 回到 $A$ 时，闭合路径的总相位差不是 $0$，而是 $-\Omega/2$。这正是”不可传递”的精确表述：

$$\langle A|B\rangle\langle B|C\rangle\langle C|A\rangle = e^{-i\Omega/2}\cdot\left|\langle A|B\rangle\langle B|C\rangle\langle C|A\rangle\right|$$

（这个恒等式可以逐点验算，例如取北极与赤道上经度相差 $60°$ 的两点，$\Omega = \pi/3$，乘积相位的大小正好是 $\pi/6 = \Omega/2$；符号随旋量内积约定取正或负，Berry 原文取 $-\Omega/2$。可验算的是”乘积相位由立体角决定、大小恰为 $\Omega/2$”，这是不依赖约定的物理内核。）

2.3 一般化：闭合路径的立体角定理

把上面的三角形推广到任意闭合路径 $C$（Berry 1987 年用旋量代数严格证明）：圆偏振态沿庞加莱球上闭合路径演化一周回到起点后，相位因子为

$$\langle A|A’\rangle = \exp\left(-\frac{i}{2} \Omega(C)\right)$$

其中 $\Omega(C)$ 是路径围成的立体角。这个相位只依赖路径围出的面积，与演化快慢、经过哪些中间点都无关——只要闭合路径围出的立体角相同，相位就相同。它因此叫几何相位：几何决定一切，动力学细节无关紧要。注意这里的负号与旋量内积的约定有关，不同文献可能写成 $+\Omega/2$，两种写法只是圆偏振正负号的约定不同，物理内容完全一致。

三、几何路线：在球面上”算”出 2σθ

代数路线已经走完，但还没回答”为什么是几何”。现在用几何语言重走一遍半波片的推导，和琼斯矩阵的结果对账——两条路线指向同一个 $2\theta$，本身就是一次交叉验证。

3.1 波片 = 球面上的旋转

庞加莱球上，一个相位延迟 $\delta$、快轴方向 $\alpha$ 的波片，等价于绕赤道上经度 $2\alpha$ 处的轴旋转 $\delta$ 角。对半波片（$\delta=\pi$），这是绕赤道某直径旋转 $180°$：

• 入射圆偏振（北极或南极）被翻转到对极点——旋向反转（自旋翻转）；
• 快轴方向 $\alpha$ 决定旋转轴在赤道上的位置，也就是决定从北极到南极走哪条经线：$\alpha$ 变化 $\theta$，球面经度变化 $2\theta$（线偏振方向 $\theta$ 与 $\theta+\pi$ 是同一个，映射到球面必须用 $2\theta$）。

3.2 两条经线夹出”月牙”：Ω = 4Δθ

把旋转角 $\theta_1$、$\theta_2$ 的两种纳米柱路径画在一起（下图）。北极是入射圆偏振，南极是出射的自旋翻转圆偏振；半波片决定了从北极到南极”走哪条经线”，而旋转角决定经线的位置。两条经线从北极出发、在南极汇合，中间夹出一个”月牙”——这个月牙的面积，就是几何相位的来源。

图 2：庞加莱球上的 PB 相位路径。快轴旋转角 $\theta$ 对应球面经度 $2\theta$；两条经线夹出的月牙立体角 $\Omega = 2\times 2\Delta\theta = 4\Delta\theta$。

两条经线在赤道上的经度差为 $2(\theta_2-\theta_1)=2\Delta\theta$。球面上两条夹角为 $\delta$ 的经线围成的月牙（球面二角形），立体角是 $2\delta$，所以：

$$\Omega = 2 \times 2\Delta\theta = 4\Delta\theta$$

代入立体角定理：

$$\varphi = \pm\frac{\Omega}{2} = \pm 2\Delta\theta$$

这里的 $\pm$ 来自闭合路径的环绕方向与旋量内积约定：走的方向相反，立体角变号，相位也跟着变号。与琼斯矩阵对照：右旋入射（$\sigma=+1$）对应 $+2\theta$，左旋（$\sigma=-1$）对应 $-2\theta$，几何路线的 $\pm$ 正好由 $\sigma$ 吸收。旋转角差 $\Delta\theta$ 先被球面几何放大成 $4\Delta\theta$ 的立体角，再被 $1/2$ 因子折回 $2\Delta\theta$ 的相位差——两倍系数在这里第二次出现，与琼斯矩阵的结果完全一致。把符号约定写成 $\sigma$（两种圆偏振互为相反数），就是 $\varphi = 2\sigma\theta$。

到这里，两条路线汇合了：琼斯矩阵从”代数上”算出 $2\theta$，庞加莱球从”几何上”解释 $2\theta$。前者回答”它确实是 $2\theta$”，后者回答”它为什么必然是 $2\theta$”。

需要说明的是，$\sigma$ 的正负约定在不同文献里可能相反（取决于从哪一端看光的旋向）。本文约定 $\sigma=+1$ 对应右旋圆偏振、$\sigma=-1$ 对应左旋圆偏振（从接收端看）。无论哪种约定，物理内核不变：两种圆偏振获得方向相反、大小相等的相位。

3.3 两倍系数从哪来：一张图收束

把琼斯矩阵和球面几何两条路线的关键对应关系放在一起：

表 1：旋转角 θ 如何变成相位 2σθ——两条路线的对照。

环节	琼斯矩阵路线	庞加莱球路线
旋转角出现的位置	旋转矩阵的二倍角：$\cos2\theta$、$\sin2\theta$	快轴方向 $\theta$ 映射到经度 $2\theta$
相位因子	圆偏振基非对角元 $e^{\pm 2i\theta}$	月牙立体角 $\Omega=4\Delta\theta$
最后一步	半波片 $\delta=\pi$ 全转给自旋翻转分量	立体角定理 $\varphi=-\Omega/2$
结果	$\varphi = 2\sigma\theta$	$\varphi = 2\sigma\theta$

四、进阶：Berry 联络、Berry 曲率与”自旋×旋转角”

前面用”指南针绕三角形走一圈偏转 Ω”给出了几何直觉。如果想再往前走一步，把”相位差密度”本身写成严格的数学对象，就得到微分几何的语言。设偏振态 $|\psi(\theta,\phi)\rangle$ 是庞加莱球上的旋量，定义 Berry 联络

$$\mathcal{A} = i \langle\psi| \mathrm{d}|\psi\rangle$$

它给出球面上相邻偏振态之间的”相位差密度”。对旋量计算可得 $\mathcal{A}$ 正比于 $-\frac{1}{2}(1-\cos\theta) \mathrm{d}\phi$，其旋度（Berry 曲率）$\mathcal{F}=\nabla\times\mathcal{A}$ 正好是一个强度为 $-1/2$ 的磁单极场——电荷就是偏振态，磁通就是立体角。闭合路径的相位：

$$\varphi_{\text{PB}} = \oint_C \mathcal{A} = -\frac{1}{2}\oint_C \mathrm{d}\Omega = -\frac{\Omega(C)}{2}$$

这正是上一节的结果，只不过现在”立体角”变成了磁单极通量，几何相位的拓扑性一目了然。如果觉得这一节太快，可以直接跳过——它不影响前面所有结论，只是给”为什么是 $-1/2$”一个更精确的答案。

这里值得区分另一类容易混淆的相位——自旋重定向相位（spin-redirection phase）。它发生在光传播方向（波矢 $\boldsymbol{k}$）在空间变化时：偏振矢量必须始终垂直于 $\boldsymbol{k}$，沿弯曲路径平行输运时偏振面会相对球坐标旋转，从而产生几何相位。两者的共同本质是”自旋 × 旋转角”：PB 相位来自各向异性结构旋转坐标系，自旋重定向相位来自波矢方向在 $\boldsymbol{k}$ 球面上的演化。前者是超表面可以”定制”的，后者是各向同性介质中固有的自旋-轨道相互作用。这也解释了为什么 PB 相位特别适合超表面：我们很难让光在器件里”拐弯”来获得自旋重定向相位，却可以自由旋转纳米柱来获得 PB 相位。

五、为什么无色散，以及如何在超表面上实现

5.1 无色散：几何量 vs 材料量

PB 相位只由旋转角决定，不依赖材料的折射率色散曲线。对比一下两种相位机制：

表 2：传播相位与 PB 相位的对比。

维度	传播相位	PB 相位
调控方式	调整纳米柱尺寸	调整纳米柱旋转角
相位来源	光程 / 等效折射率	偏振态的球面几何演化
波长依赖	强（色散明显）	无（纯几何量）
偏振要求	任意偏振	圆偏振入射，取自旋翻转分量
典型器件	常规超透镜	宽带超透镜、涡旋光、偏振复用

只要旋转角 $\theta(x,y)$ 的空间分布给定，PB 相位提供的波前在所有波长上都一样。这就像一张与颜色无关的”相位地形图”：无论红光还是蓝光经过，等高线（等相位线）的位置都相同。宽带消色差超透镜正是利用这一点提供基础相位，再用其他机制补偿色差，从而在可见光全波段聚焦成像。

5.2 实现：旋转纳米柱阵列

在超表面上实现 PB 相位，就是把”旋转半波片”做成”旋转纳米柱阵列”：

1. 单元设计：每根纳米柱是椭圆或矩形截面的双折射单元，调节长轴/短轴尺寸使两个主轴间相位差为 $\pi$（半波片条件），此时自旋翻转效率最高；
2. 空间编码：在每个位置让纳米柱旋转 $\theta(x,y)$，出射的自旋翻转分量在该点的相位就是 $2\sigma\theta(x,y)$；
3. 波前成型：给定想要的相位面，反解出每个位置的旋转角，就能设计超透镜、偏转器、涡旋光发生器、偏振复用全息等器件。

图 3：旋转纳米柱阵列与 PB 相位面。（a）纳米柱的旋转角 $\theta(x,y)$ 随位置变化（俯视图，颜色表示不同旋转角）；（b）对应的相位面示例，这里画出的是产生涡旋光（轨道角动量）的螺旋相位分布。

一个关键限制也由此而来：PB 相位只作用于自旋翻转分量。入射圆偏振经过纳米柱后，一部分保持原自旋（不获得 PB 相位），另一部分翻转自旋（获得 $2\sigma\theta$）。理想半波片条件下翻转分量占主导、效率接近 $100%$，但实际效率会随波长与材料损耗变化，这也是宽带器件需要额外补偿的原因。

PB 相位超表面已经催生了一批代表性应用：宽带消色差超透镜、携带轨道角动量的涡旋光束、左旋/右旋圆偏振各看一套相位的偏振复用全息与多通道成像，以及基于量化 PB 相位（$0$、$\pi/2$、$\pi$、$3\pi/2$）的偏振相关衍射元件。回看全文：琼斯矩阵告诉我们旋转角如何变成 $2\sigma\theta$，庞加莱球告诉我们为什么必然是 $2\sigma\theta$，而”只与旋转角有关”这一几何本质，正是它无色散、适合宽带器件的根源。一把旋转角，就是一把通往任意波前的钥匙——而钥匙的几何原理，就刻在庞加莱球上。

参考资料

• Pancharatnam, S. “Generalized theory of interference and its applications.” Proc. Indian Acad. Sci. A 44: 247 (1956). DOI
• Berry, M. V. “Quantal phase factors accompanying adiabatic changes.” Proc. R. Soc. Lond. A 392: 45 (1984). DOI
• Berry, M. V. “The adiabatic phase and Pancharatnam’s phase for polarized light.” J. Mod. Opt. 34: 1401 (1987). DOI
• Bomzon, Z., Kleiner, V., Hasman, E. “Pancharatnam–Berry phase in space-variant polarization-state manipulations with subwavelength gratings.” Opt. Lett. 26: 1424 (2001). DOI
• Hasman, E. et al. “Polarization dependent focusing lens by use of quantized Pancharatnam–Berry phase diffractive optics.” Appl. Phys. Lett. 82: 328 (2003). DOI
• Chen, W. T. et al. “A broadband achromatic metalens for focusing and imaging in the visible.” Nat. Nanotechnol. 13: 220 (2018). DOI
• Zhan, Q. et al. “Geometric phase encoded metasurfaces.” Nanophotonics（综述）. DOI
• Bliokh, K. Y. et al. “Spin-orbit interactions of light.” Nat. Photon. 9: 796 (2015). DOI