🔬 optics

Phase and coherence retrieval from near- and far-field intensities

本文介绍了一种通过两种受 Gerchberg-Saxton 算法启发的算法——一种高精度的 4D 张量方法和一种计算高效的蒙特卡洛变体——从近场和远场强度测量中检索部分相干光空间相干性的新范式，这些算法已通过在大规模激光阵列上的模拟和实验得到了验证。

原作者： Eran Bernstein, Amit Pando, Asher A. Friesem, Nir Davidson

发布于 2026-06-11

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Eran Bernstein, Amit Pando, Asher A. Friesem, Nir Davidson

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，你正试图弄清楚一群复杂舞团的形态与动作，但你只能看到两样东西：舞者在舞台上的模糊快照（“近场”）以及远处墙上投影出的影子轮廓（“远场”）。你看不见舞者的脸庞，也看不见他们的具体舞步，只能看到他们反射出的光线总亮度。

这就是光学科学家在试图理解部分相干光时面临的问题。不像完美的激光束那样每个光子都步调一致，许多光源（比如大型激光阵列或穿过雾气的光）都带有一点混沌性。有些部分是同步的，而有些则不是。为了完整理解这种光，科学家不仅需要知道光在哪里，还需要知道光的不同部分是如何相互连接的。这种连接被称为互强度或相干性。

几十年来，只有当光是完全同步时，科学家才能重建出完整的图像。如果光是混乱的（部分相干），数学模型就会失效，他们无法仅通过观察两个位置的亮度来推导出隐藏的细节。

新的“魔术技巧”

本文的作者开发了一种解决这个难题的新方法。他们称之为**“Gerchberg–Saxton (GS) 框架”**，这本质上是一个不断重复进行的智能猜谜游戏，直到答案自然浮现。

把它想象成尝试重建一面破碎的镜子。你有一张镜子正面的模糊照片和一张背面的模糊照片。你先对镜子的样子做一个随机猜测。然后，你运行你的猜测通过一个“模拟器”，来看看这些照片“应该”长什么样。

如果模拟的正向照片与真实的正面照片不符，你就微调你的猜测。
如果模拟的背面照片与真实的背面照片不符，你再次进行微调。
你重复这个循环数千次。最终，你的猜测会变得如此精确，以至于能完美匹配这两张照片，从而揭示出镜子的隐藏结构（光的相位和相干性）。

两种游戏玩法

团队创建了两个版本的猜谜游戏，以应对不同的情况：

“超级计算机”版 (Tensor GS)：
想象一下尝试解决一个 4D 谜题，其中每一块碎片都与其它所有碎片相连。这种方法正是如此。它将光视为一个巨大的四维物体，并同时计算每一个点之间的连接。
- 优点： 它极其精确，能给出完美的图像。
- 缺点： 它需要大量的计算能力，就像试图解开一个城市规模的魔方一样。
“众包”版 (Monte Carlo GS)：
与其让一个巨大的大脑去解决整个谜题，不如想象雇佣 1000 个人，每人负责解决一小块随机的谜题。然后，你将他们的答案取平均值。
- 优点： 它在计算机上运行得更快、更便宜。
- 缺点： 它是一种近似值。你雇佣的人越多（即进行的“模拟”越多），它就越精确，但它永远不会像“超级计算机”版本那样完美。

成功了吗？

团队通过两种方式测试了他们的方法：

在计算机中（模拟）： 他们创建了包含多达 600 束光束的虚拟激光阵列。他们将光设置为以特定的“混乱”方式运行（某些部分相连，某些则不相连）。
- 结果： 当光是混乱的时候，旧方法惨败。他们的两种新方法都表现出色，即使在光高度无序的情况下也能完美运作。他们能够准确测量光的“连接”程度，即使是在复杂的环形结构中。
在实验室中（真实实验）： 他们构建了一个由 130 个耦合激光器组成的真实三角形阵列。这些激光器被设计用来产生一种非常特定且棘手的、具有“交错”相位（类似于螺旋形）的光模式。
- 结果： 他们使用“超级计算机”方法，仅通过强度测量就重建了光的隐藏属性。结果与理论几乎完全吻合。其相位（光波的“时序”）误差极小，仅相当于 2π/250——这只是一个完整圆周的极小部分。

核心结论

这篇论文引入了一个全新的工具箱，让科学家能够仅通过两个简单的亮度测量，就能“看见”混乱光源中隐藏的连接。它既适用于简单也适用于复杂的的光模式，并提供了在最大精度（如果你有强大的计算机）或速度（如果你需要快速得到答案）之间进行选择。这是理解复杂激光系统和量子光的重要一步，而且无需在开始前预知答案。

技术摘要：近场与远场强度的相位与相干性检索

问题陈述
量化复杂光场的相干特性对于从高功率激光阵列到量子相干系统等应用领域至关重要。虽然迭代相位检索算法（如 Gerchberg–Saxton (GS) 框架）是根据强度测量重建复数场的标准方法，但它们通常假设场是完全相干的。许多重要的物理系统（包括大型激光阵列、X射线/电子束以及结构化量子态）仅表现出部分空间相干性。在这些情况下，场无法用简单的复振幅来描述，而需要一个四维互强度函数 $J(r, r') = \langle E^*(r)E(r') \rangle$ ，它代表了一阶空间相干性。现有的重建方法通常依赖于限制性的假设（例如，均匀相干性、已知的模态基或低秩密度矩阵），或者无法仅通过两个平面（近场和远场）的强度测量来恢复完整的互强度。

方法论
作者提出了一种新的范式，用于从两个傅里叶共轭平面的强度测量中重建部分相干场的完整互强度。该方法改进了 GS 算法，使其能够在近场和远场域之间迭代传播并投影互强度。开发了两种互补的实现方式：

张量 GS 算法（Tensor GS Algorithm）： 该方法通过对四维互强度张量应用交替投影，直接重建互强度。该过程包括：
- 使用四维傅里叶变换在平面之间传播互强度。
- 通过将互强度的对角元素投影到测得的近场和远场强度上，来强制执行目标强度约束。
- 更新非对角元素，以在强制执行埃尔米特对称性和物理一致性（确保 $|\tilde{g}_1(k, k')| \leq 1$ ）的同时，保持底层的相干结构。
- 迭代直至收敛。
蒙特卡洛 (MC) GS 变体： 为了降低与四维傅里叶变换相关的计算成本，该方法将部分相干场表示为 $P$ 个独立的随机复数实现的集合。
- 在平面之间使用标准的二维傅里叶变换传播场。
- 在每一步中，场振幅都会乘以一个共同的函数，以确保系综平均强度与目标强度相匹配。
- 通过对 $P$ 个实现进行系综平均来估计互强度。
- 精度随 $1/\sqrt{P}$ 缩放，允许受控的近似。其计算复杂度为 $PM \log(M)$ （其中 $M$ 是空间模态数），而张量 GS 方法的复杂度为 $M^2 \log(M)$ 。

主要贡献

通用重建框架： 本文提出了一种通用的方法，仅通过两个平面中的强度测量即可重建完整的空间互强度，无需预定义的模态基、低秩参数化或多个检测平面。
双重算法路径： 引入了高精度的张量 GS 方法和计算高效的 MC GS 变体，为不同的系统规模和精度需求提供了灵活性。
处理非均匀性： 这些方法专门设计用于处理空间非均匀的部分相干场，在这一领域，传统的逆傅里叶变换 (IFT) 方法会失效。

结果
作者通过数值模拟和实验演示验证了两种方法：

数值模拟：
- 线性阵列： 对最多包含 $N=600$ 个具有预设高斯衰减相干性的高斯束的线性阵列进行的模拟表明，虽然 IFT 方法仅在短相干长度（ $\xi \ll Nd$ ）下能准确恢复相干性，但 Tensor GS 和 MC GS 方法即使在相干长度超过阵列尺寸（ $\xi > Nd$ ）且互强度呈现强非均匀性时，仍能保持准确性。
- 环形阵列： 对环形阵列（ $N=50$ ）的模拟显示，由于平移不变性假设的失效，IFT 方法出现了定性失败。相比之下，所提出的两种方法都能准确恢复互强度的正确结构，其中 Tensor GS 提供高保真重建，而 MC GS 则提供计算高效的近似。
实验演示：
- 将 Tensor GS 方法应用于由 130 个具有受控耦合和交错手性的耦合激光器组成的三角形阵列。
- 该系统产生了具有非单调互相干衰减和非平凡相位的状态（交错反涡旋和涡旋/反涡旋等量叠加态）。
- 重建结果与理论高度一致，成功复现了预期的相位缠绕（例如 $2\pi/3$ ）和相干强度的高斯衰减。
- 测得交错反涡旋态的均方根 (RMS) 相位误差为 $2\pi/250$ ，等量叠加态的相位误差为 $2\pi/65$ 。

意义与主张
论文声称这些结果表明该方法对于表征大规模部分相干光学系统是实用的。这些方法通过消除对相干结构先验知识或多个检测平面的需求，补充了现有的断层扫描方法。作者指出，尽管从两个傅里叶共轭平面重建互强度的理论唯一性仍是一个开放性问题，但其数值和实验结果展示了向具有物理意义的解稳定收敛的能力，尤其是在面对强非均匀相干结构时，并表现出对实验噪声的鲁棒性。

这项工作暗示了未来的潜在扩展方向，例如结合自由空间传播从非傅里叶共轭平面重建场，以及将该框架应用于量子断层扫描，因为互强度重建与密度矩阵重建之间存在形式上的类比。然而，本文保持了审慎的态度，重点关注在特定且具有挑战性的非均匀系统中重建相位和相干性的能力，而非过度扩张到任意光学场的范畴。

新的“魔术技巧”

两种游戏玩法

成功了吗？

核心结论

类似论文