[论文解读] Asymptotic evaluation of the matrix permanents giving the bosonic probability amplitudes in linear unitary networks

V. S. Shchesnovich|arXiv (Cornell University)|Apr 24, 2013

Optical Network Technologies被引用 1

一句话总结

该论文提出了一种渐近解析方法，用于在玻色子数 $N$ 远大于模式数 $M$ 的情况下，估算有限模式线性酉网络中的玻色子振幅。通过在多维积分表示上应用鞍点法，该方法通过求解酉矩阵的标度问题，将误差控制在 $\mathcal{O}(1/N)$ 以内，从而得到二模和三模网络中简单鞍点的显式公式。

ABSTRACT

An asymptotic analytical approach is proposed for bosonic probability amplitudes in unitary linear networks, such as the optical multiport devices for photons. The asymptotic approach applies for large number of bosons $N\gg M$ in the $M$-mode network, where $M$ is finite. The probability amplitudes of $N$ bosons unitarily transformed from the input modes to the output modes of a unitary network are expressed through a multidimensional integral with the integrand containing a large parameter (N) in the exponent. The integral representation allows an asymptotic estimate of bosonic probability amplitudes up to a multiplicative error of order 1/N by the saddle point method. The estimate depends on solution of the scaling problem for the $M imes M$-dimensional unitary network matrix: to find the left and right diagonal matrices which scale the unitary matrix to a matrix which has specified row and column sums (equal, respectively, to the distributions of bosons in the input and output modes). The scaled matrices give the saddle points of the integral. For simple saddle points, an explicit formula giving the asymptotic estimate of bosonic probability amplitudes is derived. Performance of the approximation and the scaling of the relative error with N are studied for two-mode network (the beam-splitter), where the saddle-points are roots of a quadratic and an exact analytical formula for the probability amplitudes is available, and for three-mode network (the tritter).

研究动机与目标

• 开发一种在 $N \gg M$ 条件下估算线性酉网络中玻色子概率振幅的渐近方法。
• 解决在包含大量玻色子的量子光学网络中计算高维类似永久行列式振幅的挑战。
• 为这些振幅提供一种具有可控误差量级的系统性近似方法。
• 推导出如分束器和三工器等简单网络中渐近估计的显式公式。

提出的方法

• 将玻色子概率振幅表示为指数中含有大参数 $N$ 的多维积分。
• 应用鞍点法近似该积分，识别出由鞍点主导的贡献。
• 提出一个标度问题：寻找左右对角矩阵，使酉矩阵经标度后其行和与列和分别匹配输入和输出玻色子分布。
• 求解该标度问题，以确定复平面上鞍点的位置。
• 当鞍点为简单情况时，推导出振幅的显式渐近公式。
• 通过与精确解析结果对比，在二模（分束器）和三模（三工器）网络中对近似结果进行数值验证。

实验结果

研究问题

• RQ1 当玻色子数 $N$ 远大于模式数 $M$ 时，如何近似有限模式线性酉网络中的玻色子概率振幅？
• RQ2 在大 $N$ 极限下，振幅积分的主导贡献结构是什么？
• RQ3 酉矩阵的标度问题如何与渐近近似中鞍点的位置相关联？
• RQ4 该渐近近似的精度如何？相对误差随 $N$ 如何变化？
• RQ5 该方法能否为如分束器和三工器等简单网络生成振幅的显式公式？

主要发现

• 通过鞍点法，玻色子振幅的渐近近似实现了 $\mathcal{O}(1/N)$ 量级的乘法误差。
• 鞍点由一个矩阵标度问题确定，该问题通过平衡酉矩阵的行和与列和，使其匹配输入和输出玻色子分布。
• 在二模分束器中，鞍点对应于一个二次方程的根，且渐近公式与精确振幅一致，误差可控。
• 在三模三工器中，该方法成功实现了振幅的估计，且误差量级可定量预测。
• 随着 $N$ 增大，近似结果的相对误差减小，证实了该方法的渐近有效性。
• 当鞍点非退化且为简单情况时，推导出了渐近振幅的显式解析公式。