1. 量子热力学与Jarzynski等式基础
量子热力学是传统热力学在微观尺度的延伸，它研究量子系统在非平衡过程中的能量转换与熵产生。与经典系统不同，量子系统表现出独特的相干性和纠缠性，这使得量子热力学过程展现出丰富的物理现象。在量子领域，功和热的概念需要重新定义，因为量子测量本身会干扰系统状态。
Jarzynski等式由Christopher Jarzynski于1997年提出，它建立了非平衡过程中功的统计分布与平衡态自由能差之间的联系。这个等式可以表示为：
⟨e^(-βW)⟩ = e^(-βΔF)
其中β=1/(k_B T)是逆温度，W是过程中系统所做的功，ΔF是初末态的自由能差。这个等式的惊人之处在于，它适用于任意速度的驱动过程，无论是准静态还是高度非平衡的情况。
在量子领域，Jarzynski等式需要特别考虑量子测量带来的影响。量子功的定义通常采用"两次投影测量"方案：首先在过程开始时测量系统的能量，然后让系统经历非平衡演化，最后在过程结束时再次测量能量。两次测量结果的差值即为量子功。
2. 光子处理器实验设计
2.1 量子活塞模型
量子活塞是一个理想化的模型，用于研究受限量子粒子的非平衡热力学。在我们的实验中，量子活塞由一个随时间变化的无限深方势阱模拟，其哈密顿量为：
Ĥ(t) = -ℏ²/2m ∂²/∂x², 0 < x < λ(t)
其中λ(t)是随时间变化的势阱宽度。我们考虑线性驱动协议λ(t)=λ_0+vt，v是"活塞"的运动速度。
对于固定宽度λ，系统的本征态和本征能量为： φ_n(x;λ) = √(2/λ) sin(nπx/λ) E_n(λ) = ℏ²π²n²/(2mλ²)
2.2 光子处理器实现
我们使用可编程的集成光子处理器来模拟量子活塞动力学。处理器基于硅光子学平台，包含12×12的Clements型干涉仪网络。每个基本单元是一个马赫-曾德尔干涉仪(MZI)，由两个50:50多模干涉(MMI)耦合器和两个热光相位调制器组成。
实验的核心挑战是将非幺正的量子活塞演化嵌入到幺正的光学干涉仪中。我们采用准幺正嵌入技术：首先对截断的4×4活塞矩阵进行奇异值分解，然后将其扩展为5×5幺正矩阵，其中新增的模式代表泄漏到高能级的部分。
3. 实验方法与测量技术
3.1 双光子态制备
实验使用II型自发参量下转换(SPDC)过程产生纠缠光子对。泵浦激光(785nm)入射到周期极化KTP晶体，产生信号光和闲频光光子对(约1570nm)。通过精心调节相位匹配条件，我们获得了高纯度、高亮度的光子对源。
光子对的二阶关联函数g^(2)(0)测量值为0.014，表明良好的单光子特性。通过温度控制晶体，我们实现了约20%的预示效率，两光子符合计数率约20kHz(2ns符合窗)。
3.2 状态重构与测量
光子处理器输出端的10个探测通道通过超导纳米线单光子探测器(SNSPD)阵列进行测量。由于探测器数量限制，我们采用三组不同的配置覆盖所有必要的符合测量组合。
对于每个活塞参数设置，我们记录所有两光子符合事件，并通过效率校准和归一化处理，重构出完整的双光子Fock态输出分布。实验与理论的吻合度通过Bhattacharyya系数量化：
B = Σ_n √[P_exp(n)P_th(n)]
其中P_exp(n)和P_th(n)分别是实验测量和理论预测的输出概率。
4. 实验结果与分析
4.1 膨胀过程的热力学特征
在膨胀过程中(λ_0=1.0→λ_τ=3.0，T=5.0)，我们观察到以下现象：
低速区域(|v|≲1)：驱动接近绝热，系统保持在瞬时本征态。功分布P(W)集中在离散的负功值附近，反映了能级随势阱扩大而降低。
高速区域：非绝热效应显著，出现能级间跃迁。功分布展宽，出现正功成分。对于全同玻色子，跃迁振幅涉及单粒子传播子子矩阵的积和式，增强了粒子间的关联。
平均功⟨W⟩随速度单调增加，始终大于平衡自由能差ΔF_th，符合热力学第二定律。耗散功W_diss=⟨W⟩-ΔF_th随速度增加。
4.2 压缩过程的热力学特征
压缩过程(λ_0=3.0→λ_τ=1.0，T=0.5)展现出不同的行为：
即使准绝热情况下，功值也较大且为正，反映了压缩过程固有的不可逆性。
高速驱动时，功分布呈现明显的高功尾巴，这是能级快速靠近导致非绝热耦合增强的结果。
对于全同玻色子，积和式增强的跃迁振幅进一步放大了这种重分布。
4.3 Jarzynski等式的验证
通过Jarzynski自由能估计器：
ΔF_exp = -T ln⟨e^(-W/T)⟩
我们在膨胀和压缩过程中都观察到了ΔF_exp与理论值ΔF_th的出色吻合。在膨胀过程中，偏差不超过0.07T；在压缩过程中，即使⟨W⟩达到2757T的极端值，ΔF_exp仍能准确恢复ΔF_th。
这一结果验证了Jarzynski等式在强驱动量子系统中的有效性，表明即使在高不可逆条件下，通过指数重加权，罕见的低功轨迹仍能精确补偿平均超额功。
5. 量子干涉对功统计的影响
全同玻色子的量子统计对功分布有显著影响。与可区分粒子相比，玻色子的跃迁振幅涉及单粒子传播子子矩阵的积和式(permanent)，这导致：
- 增强的能级间跃迁概率
- 功分布的更宽展宽
- 正功成分的增强
这种量子干涉效应可以通过比较相同条件下可区分粒子和全同玻色子的功分布来量化。实验结果显示，玻色子统计导致功分布的方差增加约30%，证实了量子干涉对热力学过程的显著影响。
6. 循环协议与不可逆性
我们研究了由膨胀和压缩组成的完整热力学循环。定义循环耗散功：
W_diss,cyc = (⟨W⟩_exp - ΔF_exp) + (⟨W⟩_comp - ΔF_comp)
即使在最低驱动速度下，W_diss,cyc也不为零，这是量子涨落的直接结果。随着速度增加，非绝热激发进一步增加了不可逆性。
通过Bhattacharyya重叠B_cycle量化初末态的一致性，我们发现：
- 低速(|v|≲0.5)：B_cycle > 0.99，接近可逆过程
- 中速(|v|≈2.5)：B_cycle < 0.95
- 高速(|v|≈4.1)：B_cycle < 0.90
这种退化反映了膨胀过程中产生的高激发态在压缩过程中无法完全退激发。
7. 技术细节与误差分析
7.1 准幺正嵌入误差
由于将非幺正演化嵌入幺正光学网络的限制，我们引入了准幺正误差：
ε_unitary(%) = 100|1-F|, F = (1/d)Tr[√(UU†)]
实验要求ε_unitary < 5%。通过分析，我们确定：
- 95%状态保真度要求ε_unitary ≲14.8%
- 90%保真度容忍ε_unitary ≈31.6%
7.2 系统校准与稳定性
实验系统采用多项稳定措施：
- 光子芯片安装在帕尔贴冷却平台上，温度稳定性±0.01°C
- 每个MZI相位调制器通过精密电流源控制，分辨率740nA
- 定期用经典激光校准干涉仪，确定电流-相位映射关系
- SNSPD工作在2.6K，探测效率>70%，暗计数<10Hz
这些措施确保了实验数据的高重复性和可靠性。
8. 扩展与应用前景
本研究建立的框架可扩展到更多玻色子情况。N玻色子的功分布涉及N×N矩阵的积和式，其经典计算复杂度随N指数增长。因此，光子处理器为研究多体量子热力学提供了独特优势。
潜在应用包括：
- 量子热机设计与优化
- 量子制冷协议评估
- 量子电池性能表征
- 量子信息处理中的热管理
未来的工作将探索更复杂的驱动协议、更高粒子数情况，以及与其他量子平台的结合。
