[论文解读] Multiphoton Emission

Guillermo Díaz-Camacho, Eduardo Zubizarreta Casalengua+5 more|arXiv (Cornell University)|Sep 24, 2021

Mechanical and Optical Resonators被引用 7

一句话总结

本文提出了一套基本的、可解析求解的多光子发射与探测框架，严格推导出频率滤波后 Fock 态自发发射的完整时间结构。研究证明，空腔中两能级系统发出的多光子遵循检测概率的二项分布，给出了光子计数统计和联合时间分布的精确闭式表达式，揭示了滤波如何塑造量子统计特性，并实现超越标准自发发射的鲁棒发射。

ABSTRACT

We describe the emission, detection and structure of multiphoton states of light. We include the effect of frequency filtering, which describes, at a fundamental level, physical detection of a quantum emitter. The case of the spontaneous emission of Fock states is treated fully and analytically. We stress this picture by contrasting it to the numerical simulation of two-photon bundles emitted from a two-level system in a cavity. We show that dynamical factors exist that allow for a more robust multiphoton emission than spontaneous emission. We also describe how this relates to thermal light.

研究动机与目标

• 提供多光子发射与探测的严格、基础性描述，包括频率滤波在物理测量中的作用。
• 推导频率滤波后 Fock 态自发发射中光子计数概率的精确闭式表达式。
• 将自发发射与连续波及热发射进行对比，突出滤波在塑造量子统计特性中的作用。
• 表征多光子脉冲的完整时间结构，包括联合与边缘发射时间分布。
• 探讨滤波如何将非热场（如单光子发射）转化为有效热类似态，并识别由有限带宽探测引入的非热光谱与关联特征。

提出的方法

• 使用 Mandel 公式推导完整的光子计数概率 p(n, T; N)，引入时间积分强度算符 Ω，结合探测效率 ξ 和频率滤波。
• 应用 Lindblad 主方程建模两能级系统发射 N 个 Fock 态光子的自发发射过程，解析计算多时间关联函数。
• 引入滤波场算符 ς(t) 和时间平均强度算符 ΩΓ(T)，以模拟具有有限带宽的物理探测，从而得到广义探测效率 TΓ(T)。
• 推导 N 个光子检测时间的联合概率密度函数 φ(N)Γ(t1,…,tN)，表明其为受滤波与衰减动力学调制的指数衰减项的乘积。
• 计算单个光子发射时间的边缘分布 φ(N)Γ,k(tk)，从而可计算平均值与方差。
• 分析未滤波（Γ→∞）与零带宽（Γ→0）的极限情况，揭示其与广义调和数及非洛伦兹谱形的关联。

实验结果

研究问题

• RQ1 频率滤波后，N 光子 Fock 态自发发射的光子计数概率的精确解析形式是什么？
• RQ2 频率滤波如何影响多光子发射脉冲的时间结构与检测统计特性？
• RQ3 多光子脉冲中单个光子的统计特性（如平均发射时间、方差与关联）是什么？
• RQ4 滤波如何将非热场（如单光子发射）转化为有效热类似态？其与真实热行为的偏离特征为何？
• RQ5 有限探测带宽对观测统计特性的影响如何，特别是在滤波带宽极低或极高时？

主要发现

• 从自发 N 光子发射中，在时间窗口 T 内检测到 n 个光子的概率服从二项分布：p(n, T; N) = (N choose n) × T(T)^n × (1 - T(T))^(N-n)，其中 T(T) = ξ(1 - e^(-γa T)) 为有效探测效率。
• 对于频率滤波发射，探测效率变为 TΓ(T) = Γ/Γ+ - Γ² e^(-γa T) + Γγa e^(-ΓT) / (Γ² - Γ+)，该式推广了未滤波情况，且考虑了有限探测器带宽的影响。
• N 个光子检测时间的联合概率密度函数为 φ(N)Γ(t1,…,tN) = N! γa^N × (Γ/Γ-)^2N × ∏_{i=1}^N (e^(-Γ ti/2) - e^(-γa ti/2))^2 × 1_{[t_{i-1}, t_{i+1}]}(ti)，完整描述了光脉冲的时间结构。
• 第 k 个光子发射时间的边缘分布为 φ(N)Γ,k(tk) = - (Γ/Γ-)^2N × γa^k × (N choose k) × g(tk)^(N-k) × (g(0) - g(tk))^(k-1) × g’(tk)，其中 g(t) 定义为 e^(-γa t)/γa + e^(-Γ t)/Γ - 4 e^(-Γ+ t/2)/Γ+。
• 在未滤波极限（Γ→∞）下，脉冲总持续时间 τN 的平均值与方差为 ⟨τN⟩γa = H_{N-1} 与 σ∞,τN = √(H_{N-1,2}) / γa，其中 H_{N-1} 与 H_{N-1,2} 为广义调和数。
• 用窄带洛伦兹滤波器处理热场并不会得到热场；相反，它产生非洛伦兹谱 Sth,Γ(ω) 与二阶相关函数 g(2)th,Γ(τ)，即使 g(2)(0) = 2，其形式也偏离热场行为。