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)λ_0vtv是活塞的运动速度。对于固定宽度λ系统的本征态和本征能量为 φ_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 膨胀过程的热力学特征在膨胀过程中(λ_01.0→λ_τ3.0T5.0)我们观察到以下现象低速区域(|v|≲1)驱动接近绝热系统保持在瞬时本征态。功分布P(W)集中在离散的负功值附近反映了能级随势阱扩大而降低。高速区域非绝热效应显著出现能级间跃迁。功分布展宽出现正功成分。对于全同玻色子跃迁振幅涉及单粒子传播子子矩阵的积和式增强了粒子间的关联。平均功⟨W⟩随速度单调增加始终大于平衡自由能差ΔF_th符合热力学第二定律。耗散功W_diss⟨W⟩-ΔF_th随速度增加。4.2 压缩过程的热力学特征压缩过程(λ_03.0→λ_τ1.0T0.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指数增长。因此光子处理器为研究多体量子热力学提供了独特优势。潜在应用包括量子热机设计与优化量子制冷协议评估量子电池性能表征量子信息处理中的热管理未来的工作将探索更复杂的驱动协议、更高粒子数情况以及与其他量子平台的结合。
