开放量子系统中的最优耗散路径：最小化环境热耗散的变分框架

《Frontiers in Physics》：Optimal dissipative pathways in open quantum systems: a variational framework for environmental heat dissipation minimization

　　

本文研究了开放量子系统中的最优控制问题，重点关注一个两级逻辑量子比特系统，其时间演化由Lindblad (GKSL)形式体系描述的马尔可夫动力学所支配。研究目标是在状态转移过程中最小化传递给环境的净热量，同时明确考虑外部控制场所做的功。在此背景下，研究人员提出了一个变分最优控制框架，用于在保真度约束下最小化耗散热，该框架在有界控制和下半连续性的标准假设下，确立了至少存在一个有界控制，能够在保持终态保真度高于指定阈值的同时，最小化环境热耗散泛函。该分析结合了变分法的直接方法与庞特里亚金最大值原理。变分法的直接方法（在L^∞空间中）通过证明容许控制的紧致性和成本泛函的下半连续性来保证解的存在性。其次，开放量子系统的庞特里亚金最大值原理确立了基本的必要最优性条件，导出了耦合的状态-伴随方程和定义最优控制的切换函数。研究人员通过在Lindblad框架内定义内能、热流和功，为该优化问题奠定了物理基础，从而得出了耗散能的热力学相干表达式。该公式既允许进行解析理解，又允许基于MDR-GRAPE算法进行实际的数值近似。为了证明数值一致性，研究人员提出了MDR-GRAPE算法，这是一种改进的基于梯度的算法，将GRAPE扩展到耗散环境，集成了保持密度矩阵正性的技术。数值实验证明了在研究中考虑的五个独立初始化中具有可重复的收敛性，为持续收敛于可行的低耗散热解提供了经验证据。

物理系统从激发态向平衡态或基态的转变是量子力学、热力学、凝聚态物理以及生物系统等领域的基础过程。在此类转变中，能量通过辐射、声子发射等耗散通道传递给环境。尽管最终的平衡态由系统哈密顿量明确界定，但连接初态和终态的动力学路径并未由最小环境热耗散原理唯一规定，而是受限于系统的固有弛豫率及其与环境耦合的支配。因此，自然发生的弛豫过程在环境热耗散目标上通常不是最优的。这种热力学低效在纳米级量子器件、超导量子比特处理器等对能量高度敏感的量子技术中尤为关键，因为这些设备中每次操作微小的耗散减少都会累积成显著的性能提升。现有理论框架无法完全解释这一问题：最小作用量原理不显式关注状态转变中环境热耗散的减少；非平衡态热力学虽量化熵产生却不提供最小化耗散的构造性原理；传统量子最优控制理论（如GRAPE算法）旨在最大化态转移保真度，将热耗散视为副作用。针对这一理论空白，研究人员在《Frontiers in Physics》上发表了相关研究，提出并建立了一个针对耗散量子系统的变分最优控制框架。该研究证明了在满足保真度约束的条件下，存在至少一个物理上容许的控制路径，能够以严格低于自然弛豫的热耗散将系统驱动至基态。这一研究不仅具有基础理论意义，对于量子计算和量子传感等领域的能耗优化也具有重要的工程指导价值。

研究人员主要采用了以下关键技术方法开展研究：首先，基于GKSL（Gorini-Kossakowski-Sudarshan-Lindblad，描述马尔可夫开放量子系统动力学标准形式）形式体系构建了两级量子比特系统的马尔可夫开放系统动力学模型，包含幅度阻尼和纯退相干通道；其次，结合变分法的直接方法和庞特里亚金最大值原理，在L^∞空间中建立了环境热耗散泛函的最小化优化模型；最后，开发了MDR-GRAPE（Minimum Dissipation Relaxation-Gradient Ascent Pulse Engineering，最小耗散弛豫-梯度上升脉冲工程）算法，这是一种扩展自标准GRAPE并适用于耗散环境的梯度优化算法，通过引入保正性积分技术，确保密度矩阵在数值演化过程中的物理有效性。

研究结果部分总结如下：

3.1 Sensitivity of the fidelity penalty（保真度惩罚的敏感性）
通过分析软惩罚权重在优化过程中的作用，研究表明增加惩罚权重能持续促进算法收敛至满足保真度阈值的解，而自动缩放的惩罚机制在保持保真度和最小化能量之间提供了平衡的折中。实验得出自然弛豫和优化后的环境热耗散值，证明优化后的轨迹在保持0.998640高保真度的同时，热耗散相对自然弛豫减少了87.79%，验证了软惩罚策略的鲁棒性。

3.2 Parameter sweep（参数扫描）
通过对不同幅度阻尼率γ₁（范围从0.05到0.30）进行参数扫描并重新计算自然和最优耗散，研究得出在所考虑的所有数值实验参数组合中，MDR优化目标均得到满足。这表明优化框架在测试的阻尼参数范围内保持了有效性，尽管结果不应被解释为与阻尼率独立相关。

3.3 Entropy evolution（熵演化）
通过追踪非受控弛豫路径和MDR-GRAPE优化轨迹的冯·诺依曼熵，研究发现自然路径中熵从0升至最大值0.693（对应最大混合态）再随系统接近基态而下降；而在最优轨迹中，控制场重塑了熵的分布，最终态保留了相对更高的纯度。该两阶段熵行为符合幅度阻尼开放量子系统的预期物理行为，未违反热力学原理。

讨论部分总结指出，研究人员的数值实验结果与提出的变分框架一致，证明自然GKSL弛豫轨迹通常并非最小环境热耗散的路径。MDR-GRAPE最优控制引导量子比特偏离自然阻尼路径，穿越局部抑制瞬时耗散率的状态空间区域。尽管庞特里亚金分析提示了最优控制结构，但实际数值计算受软惩罚和二次正则化影响，呈现平滑近似而非严格的继电控制解。熵曲线的先升后降反映了优化器先建立相干性以获得控制优势，再利用相干性高效驱动系统。该MDR框架与传统最大化保真度的方法在根本上有别，其核心目标是热力学层面的。研究存在一定局限性，如局限于马尔可夫近似下的两能级系统，且使用一阶梯度方法在非凸景观上不保证全局最优。

研究结论部分翻译如下：本研究构建了最小耗散弛豫（MDR）优化框架，该理论结果表明，在有界控制和GKSL动力学的标准假设下，所提出的优化问题至少存在一个满足规定保真度要求的最小化解。该命题通过两种互补的方法进行阐述：变分法的直接方法利用紧致性和成本泛函的下半连续性确立了解的存在性，庞特里亚金最大值原理导出了最优性的必要条件并刻画了最优控制的结构。为证明可构造性，研究人员开发了MDR-GRAPE算法，通过保正性积分确保所有计算状态的物理有效性。数值实验表明，该优化框架在测试的阻尼参数范围内依然有效。由于当前研究并非旨在独立改变阻尼率，未来的工作将进行系统的固定时间研究以分离阻尼率的影响。由于外部控制场提供的功未纳入优化目标，报告的结果不应被解释为最小化受控过程的完整热力学成本。未来向多级系统、非马尔可夫环境以及多量子比特架构的扩展，以及超导量子比特和囚禁离子平台的实验一致性验证，是后续研究的重要方向。