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

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

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

\[\varphi = 2\sigma\theta\]

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

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

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

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

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

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\)（归一化）。叠加后的强度：

当 \(\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 相位，就是把"旋转半波片"做成"旋转纳米柱阵列"：

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

图 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

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

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