[论文解读] Out-of-equilibrium operation of a quantum heat engine: The cost of thermal coupling control

Michael H. Wiedmann, Jürgen T. Stockburger+1 more|arXiv (Cornell University)|Mar 27, 2019

Advanced Thermodynamics and Statistical Mechanics参考文献 58被引用 13

一句话总结

本文提出了一套非微扰、微观的非平衡量子热机理论框架，明确考虑了热库耦合与解耦过程中相关的功。研究揭示，这种耦合功 $W_I$ 是能量平衡中的关键且常具破坏性的组成部分——但可通过量子关联实现抑制，从而借助最优控制协议提升效率。

ABSTRACT

Real quantum heat engines lack the separation of time and length scales that is characteristic for classical engines. They must be understood as open quantum systems in non-equilibrium with time-controlled coupling to thermal reservoirs as integral part. Here, we present a systematic approach to describe a broad class of engines and protocols beyond conventional weak coupling treatments starting from a microscopic modeling. For the four stroke Otto engine the full dynamical range down to low temperatures is explored and the crucial role of the work associated with the coupling/decoupling to/from reservoirs in the energy balance is revealed. Quantum correlations turn out to be instrumental to enhance the efficiency which opens new ways for optimal control techniques.

研究动机与目标

• 为解决微观量子热机中因缺乏时间尺度分离而导致经典热力学假设失效的问题。
• 发展一种系统性、非微扰的量子热机处理方法，实现时间控制的热耦合，超越弱耦合近似。
• 量化耦合功 $W_I$ 在能量平衡中的作用及其对效率和功率输出的影响。
• 探索如何利用量子关联减少 $W_I$ 并提升热机性能。
• 为强耦合与有限时间动力学下的量子奥托循环提供理论建模与优化框架。

提出的方法

• 构建了一个包含工作介质、热库及时间依赖耦合算符 $H_{I,c/h}(t)$ 的微观哈密顿量模型。
• 将费曼-弗尔顿影响泛函方法映射为随机刘维尔-冯诺依曼方程，以描述非马尔可夫性与强耦合动力学。
• 该方法引入时间依赖的耦合强度 $\lambda_{h/c}(t)$，实现对热接触时序与持续时间的完全控制。
• 能量平衡被分解为驱动功 $W_d$、耦合功 $W_I$ 和热量交换 $Q$，其中 $W_I$ 通过 $\dot{\lambda}^2(t)$ 显式计算。
• 该方法应用于谐振子与非谐振子模型，模拟具有有限时间冲程的四冲程奥托循环。
• 通过 $qp$-关联量化量子关联，并分析其相位对 $W_I$ 的依赖关系，以识别控制策略。

实验结果

研究问题

• RQ1 耦合功 $W_I$ 如何影响量子热机的能量平衡与效率？
• RQ2 在强耦合、非平衡量子热机中，能否利用量子关联减轻 $W_I$ 的负面影响？
• RQ3 耦合时间 $\tau_I$ 与耦合强度 $\gamma$ 在决定系统作为热机或耗散器运行方面起何作用？
• RQ4 量子奥托热机的效率如何依赖于非谐性与热耦合参数？
• RQ5 通过操控 $qp$-关联，能否通过耦合协议的最优控制提升效率？

主要发现

• 在低温与强耦合条件下，耦合功 $W_I$ 是量子热机能量平衡中的主导且不可忽略的贡献，此时 $|W_I| \sim |W_d|$。
• 系统根据耦合时间 $\tau_I$ 与耦合强度 $\gamma$ 的不同，从量子热机（$\eta > 0$）转变为耗散器（$\eta \leq 0$），相图可明确标识这些区域。
• 为实现热机运行，需满足临界阈值 $\omega_0\tau_I > \frac{2\gamma}{\Delta\omega}\frac{1+R}{1-R}$，其中 $R$ 取决于循环终点的位置涨落。
• 忽略 $W_I$ 会导致对功率输出的误导性预测，包括错误的峰值高度与最优循环时间，如图 5(a) 所示。
• 对于非谐势，刚度增加会抑制位置涨落，从而降低效率，其效应类似于更强的热耦合。
• 通过调节耦合/解耦过程中 $qp$-关联的相位，如图 6 所示，可调控 $W_{I,qm}$：通过引入单位演化段可反转 $W_{I,qm}$ 的符号，展示出显著的控制潜力。