---
title: "复算一份高斯相干模型，再读前沿实验"
description: "固定 2026 年相干过滤预印本的高斯密度模型，用独立数值积分核对纯度公式，并区分模型复算、实验重建与开放问题。"
canonical_url: "https://sandboxphysics.com/zh/stories/time-frequency-quantum-optics/paper-benchmarks/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/time-frequency-quantum-optics/paper-benchmarks/"
last_modified: "2026-10-28"
---

[时频模式量子光学](https://sandboxphysics.com/zh/stories/time-frequency-quantum-optics/)

L4 · 第 8 章

# 复算一份高斯相干模型，再读前沿实验

Sandbox Physics · 2026-10-28

要读懂相干过滤实验，先把一份能复算的模型算透。这里固定一篇预印本的高斯谱密度矩阵，用直接积分计算纯度；再回到实验图，看哪些结论还需要装置、采集和重建来支撑。

## 开始前需要什么

第三章的纯度，第六章的重建与留出，第七章的转换核。高斯积分在本章逐步说明；附脚本可复算数值。

[回看上一章](https://sandboxphysics.com/zh/stories/time-frequency-quantum-optics/conversion-kernels/)

- 从谱宽与相干宽推导高斯模型纯度。
- 解释独立积分与网格收敛检查。
- 辨认论文理论图、测量图和重建图。

## 先固定论文版本和要复算的那一项

Bhattacharjee 等人在 2026 年八月提交的相干时间模式过滤预印本，比较闲频端相干预告与强度过滤。阅读对象固定为 arXiv 的第一版。主文图三比较高斯滤波条件下的理论分布、测得分布与重建纯度；补充材料第四节给出一份可以单独检验的部分相干高斯模型。

下面复算的是补充式十二至十三的密度矩阵积分关系。没有从主文图三抄取实验点，也没有把无量纲宽度比替换成实际器件的太赫兹带宽。这样能先检验公式与数值方法，再决定需要哪些实验资料来做更完整的比较。

对应阅读：A. Bhattacharjee et al. · Sandbox Physics FR10 independent numerical benchmarks

## 沿对角线看谱宽，离开对角线看相干

取无量纲频率 $x=(\omega-\omega_0)/\sigma_s$，并定义宽度比 $r=\sigma_s/\sigma_c$。$\sigma_s$ 是高斯振幅包络的宽度参数，$\sigma_c$ 控制不同频率之间相干随间隔衰减的尺度。它们不是任意互换的强度半高全宽。

沿矩阵对角线，两频率相同，相干衰减因子等于一，所以归一化只由谱包络确定。离开对角线后，频率间隔越大、相干宽越窄，非对角元素就衰减得越快。纯度是整个矩阵模平方的积分，因而会随非对角相干减少而下降。

$$
\Gamma(x,x')=\frac{1}{\sqrt\pi}\exp\!\left[-\frac{x^2+x'^2}{2}\right]\exp\!\left[-\frac{r^2(x-x')^2}{2}\right]
$$

在这份无量纲定义下 $\int\Gamma(x,x)\,dx=1$。这是论文补充式十二经变量缩放和显式归一化后的同一模型。

## 把双重积分拆成两个普通高斯积分

在纯度积分中令 $u=(x+x')/\sqrt2$、$v=(x-x')/\sqrt2$。这是旋转坐标，面积元不变。平方后的指数成为 $-u^2-(1+2r^2)v^2$，两个方向可以分别积分。利用高斯积分随宽度缩放的规律，得到下式。

独立脚本直接在原来的两个频率坐标上对密度矩阵模平方作梯形积分，再除以数值迹的平方。它没有调用实验台函数。先从负六到正六分成二百四十个间隔，再从负八到正八分成四百八十个间隔，分别计算，检查截断和离散误差。五个宽度比的最大绝对误差约为十的负十三次方。

$$
P=\iint|\Gamma(x,x')|^2\,dx\,dx'=\frac{1}{\sqrt{1+2r^2}}
$$

这项结果对应论文补充式十三，另由直接积分计算。相干宽无限大时为纯态；谱宽相对相干宽增大时纯度下降。

$$
B=\Delta\Omega\,\Delta T=\sqrt{1+r^2},\qquad P=\frac{1}{\sqrt{2B^2-1}}
$$

这两式按同一高斯模型的补充式十五至十六计算，是解析关系的代入展示；本次没有另做 chronocyclic Q 函数边缘宽度的独立数值积分。使用时需要先确认高斯形状与相干假设。

纵轴：纯度 横轴：谱宽与相干宽之比

曲线为论文高斯模型的解析式，点为五个无量纲宽度比上的原创直接积分。两者都是模型计算，没有实验测点。

原创确定性复算；每个数都由保存的模型与脚本生成。

| 谱宽与相干宽之比 | 纯度 | 边缘宽度积 |
| --- | --- | --- |
| 0 | 1.000000 | 1.000000 |
| 0.5 | 0.816497 | 1.118034 |
| 1 | 0.577350 | 1.414214 |
| 2 | 0.333333 | 2.236068 |

对应阅读：A. Bhattacharjee et al. · Sandbox Physics FR10 independent numerical benchmarks

## 回到图三，先问每一幅图从哪里来

主文图三的理论分布是前向计算；测得分布来自模式探针的符合扫描；谱密度矩阵和纯度来自带正定与归一化条件的重建。三者各回答一层问题。高斯宽度积是有用的直观量，但论文的纯度估计依靠密度矩阵重建，不能把两条过程混写。

补充材料还显示，有限扫描窗口可能让不同重建态在已测区域给出很相近的分布。拟合图像很像，仍可能留下状态上的不唯一。第六章的留出投影让这个问题在两模式里变得可操作；真实连续态还需要足够宽的测量窗口和经过标定的探测效应。

对应阅读：A. Bhattacharjee et al.

## 接下来要研究的尺度

可编程多输出分选把问题从“选一个模式”推进到“可靠地区分一组模式”。2024 年预印本所述工作后来发表于 Optics Express；它的器件测量效应、输出分辨率与输入编码是继续阅读的入口。下一步应把已标定效应代入条件态预测，再用未参与标定的模式检验。

更高增益还会引入多对发射；漂移会把长时间采集中的不同状态混在一起；实际损失可能依赖模式。这些都需要增加状态或测量模型。两模式算例已经把相位、条件选择和计数分开，扩展时可以逐项追踪，比较哪里首先偏离理想预测。

对应阅读：L. Serino et al.

## 动手检验预测

### [用留出模式检查一次条件态预测](https://sandboxphysics.com/zh/labs/quantum-pulse-gate/)

1. 选一份源和一个相干模式，先写下纯度与成功率预测。
2. 采集合成计数，用三个基重建，再比较留出投影。
3. 选择第二转换模式预设并重采，解释变化属于理想模型、器件选择还是有限次数。

## 练习与解释性答案

1. 若高斯模型的谱宽等于相干宽，纯度与边缘宽度积各是多少？

   展开解释

   宽度比为一，纯度为 $1/\sqrt3$，边缘宽度积为 $\sqrt2$。代入宽度积公式可再次得到同一纯度。这个一致性检查仍在同一高斯假设内。

2. 数值积分与解析式相符到十三位，是否证明实验装置的纯度也是这个数？

   展开解释

   它证明所定义高斯模型的这项计算相符。实验态是否符合该模型、仪器是否校准、扫描窗口是否充分，都要由测量与误差分析回答。把这两层证据区分开，复算才对读论文有帮助。

[上一章：真实脉冲门怎样接近单模式投影](https://sandboxphysics.com/zh/stories/time-frequency-quantum-optics/conversion-kernels/)

## 来源与计算范围

1. [L. Serino et al., Programmable time-frequency mode-sorting of single photons with a multi-output quantum pulse gate, arXiv:2410.03606v1; DOI 10.1364/OE.544206](https://arxiv.org/html/2410.03606v1)

   第二节和图二说明转换核及理想投影、非理想测量效应的区别；第三至四节讨论探测器层析与输出分辨率。该工作已有 Optics Express 期刊 DOI；这里固定核读 2024 年预印本版本。

2. [A. Bhattacharjee et al., Coherent temporal filtering of multimode parametric down-conversion using a quantum pulse gate, arXiv:2608.12544v1 (12 August 2026 preprint)](https://arxiv.org/html/2608.12544v1)

   核读正文式一至五、图一至四和补充材料第一至七节，以确定条件过滤与重建。有限复算选补充式十二至十三；式十五至十六只作解析代入。未复现图三实验点、论文层析算法或完整器件。

3. [Sandbox Physics FR10 independent numerical benchmarks](https://github.com/yangzichao/sandboxphysics/blob/main/docs/frontier/fr10/benchmarks.json)

   原创归一化高斯谱密度直接积分与独立相位矩阵、特征值计算。reproduce-benchmarks.mjs 保存窗口、分辨率和误差；数据不含历史测量。