pith. machine review for the scientific record. sign in

arxiv: 2602.01688 · v2 · submitted 2026-02-02 · 🪐 quant-ph

Recognition: 2 theorem links

· Lean Theorem

Optimal Control to Minimize Dissipation and Fluctuations in Open Quantum Systems Beyond Slow and Rapid Regimes

Yuki Kurokawa , Yoshihiko Hasegawa
Authors on Pith no claims yet

Pith reviewed 2026-05-16 08:42 UTC · model grok-4.3

classification 🪐 quant-ph
keywords optimal controlopen quantum systemsdissipationfluctuationsquantum thermodynamicsspin-boson modelquantum dotnumerical optimization
0
0 comments X

The pith

Optimal driving protocols for minimizing dissipation and fluctuations in open quantum systems switch discontinuously at intermediate timescales.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper uses numerical optimization to find driving protocols that reduce both dissipated work and its fluctuations in open quantum systems, focusing on speeds between the slow and rapid limits where standard approximations fail. In the coherent spin-boson model, the best protocol jumps abruptly between different shapes depending on how much weight is given to dissipation versus fluctuations. In a single-level quantum dot with a fermionic bath, the optimized protocol instead forms a series of distinct steps. These structures allow tighter control over energy loss and noise in quantum devices operating at realistic intermediate speeds.

Core claim

By numerically optimizing the driving protocols, the authors demonstrate that open quantum systems exhibit distinct optimal structures not captured by the conventional limits. Specifically, in the coherent spin-boson model, the optimal protocol switches discontinuously between distinct locally optimal solutions as the relative weight between dissipation and fluctuations is varied. Furthermore, for a single-level quantum dot coupled to a fermionic reservoir, the optimized protocol develops a characteristic multi-step structure.

What carries the argument

Numerical optimization of time-dependent driving protocols that minimize a weighted sum of dissipated work and work fluctuations (via the two-point measurement scheme) in the coherent spin-boson model and the single-level quantum dot model.

If this is right

  • The optimal protocol can be tuned by changing the relative weight to select among qualitatively different driving shapes.
  • Multi-step protocols emerge as the natural optimum in systems with fermionic reservoirs at intermediate driving times.
  • Control strategies must account for discontinuous changes rather than assuming smooth interpolation between slow and fast limits.
  • Intermediate-speed operation can achieve lower combined dissipation and fluctuations than either extreme regime allows.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • These discontinuous and stepped protocols could be tested in superconducting qubit or trapped-ion platforms where spin-boson dynamics are realizable.
  • The same numerical approach might reveal similar structural features when applied to quantum heat engines or refrigerators running at practical cycle times.
  • Universal features of optimal protocols at intermediate timescales may appear across a wider class of open quantum systems beyond the two models studied.

Load-bearing premise

Numerical optimization reliably locates the globally optimal protocols rather than local minima, and the chosen models capture the essential physics of the target open quantum systems.

What would settle it

An experiment realizing the coherent spin-boson model and varying the relative weight between dissipation and fluctuations would show an abrupt change in the measured optimal protocol shape at the predicted switching points.

Figures

Figures reproduced from arXiv: 2602.01688 by Yoshihiko Hasegawa, Yuki Kurokawa.

Figure 1
Figure 1. Figure 1: FIG. 1. Optimized control field [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Optimized control field [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Optimized control field u(t) obtained by the GRAPE [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Pareto front for the quantum-dot model at [PITH_FULL_IMAGE:figures/full_fig_p004_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Optimized control protocols [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
read the original abstract

Optimal control is a central problem in quantum thermodynamics. When minimizing dissipated work and work fluctuations defined via the two-point measurement scheme in open quantum systems, existing approaches largely focus on the rapid- and slow-driving limits, leaving the behavior at intermediate timescales elusive. In this work, by numerically optimizing the driving protocols, we demonstrate that the open quantum systems exhibit distinct optimal structures not captured by the conventional limits. Specifically, in the coherent spin-boson model, we find that the optimal protocol switches discontinuously between distinct locally optimal solutions as the relative weight between dissipation and fluctuations is varied. Furthermore, for a single-level quantum dot coupled to a fermionic reservoir, the optimized protocol develops a characteristic multi-step structure.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript numerically optimizes driving protocols to minimize a weighted combination of dissipated work and work fluctuations in open quantum systems at intermediate timescales. In the coherent spin-boson model, the authors report that optimal protocols switch discontinuously between distinct locally optimal solutions as the relative weight between dissipation and fluctuations is varied. In a single-level quantum dot coupled to a fermionic reservoir, the optimized protocols exhibit a characteristic multi-step structure.

Significance. If the reported protocols are globally optimal, the work would usefully demonstrate that intermediate-regime optimal control exhibits structures absent from the slow- and fast-driving analytic limits, potentially guiding future theory. The numerical approach is a reasonable way to explore the non-convex landscape, but the absence of convergence diagnostics and global-optimality checks reduces the strength of the central claims.

major comments (2)
  1. [Abstract and Numerical Optimization section] The central claims rest on numerical minimization of the weighted cost functional, yet no description is given of the optimization algorithm, control discretization, convergence tolerances, or multi-start/basin-hopping diagnostics that would establish the reported discontinuous switches as global rather than local minima (see Abstract and the spin-boson results paragraph).
  2. [Results for the coherent spin-boson model] For the spin-boson model, the discontinuous switch between locally optimal solutions is presented as a key finding, but without comparison to known analytic limits at the slow- and rapid-driving boundaries or exhaustive sampling statistics, it is unclear whether the discontinuity survives a more thorough search of the control landscape.
minor comments (1)
  1. [Abstract] The abstract should explicitly name the weighting parameter and the two models studied to improve immediate readability.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We have revised the paper to address the concerns about missing numerical details and validation against analytic limits, which we believe strengthens the central claims regarding the distinct optimal structures at intermediate timescales.

read point-by-point responses

  1. Referee: [Abstract and Numerical Optimization section] The central claims rest on numerical minimization of the weighted cost functional, yet no description is given of the optimization algorithm, control discretization, convergence tolerances, or multi-start/basin-hopping diagnostics that would establish the reported discontinuous switches as global rather than local minima (see Abstract and the spin-boson results paragraph).

    Authors: We agree that the original manuscript omitted key details on the numerical procedure. In the revised version, we have added a dedicated subsection 'Numerical Methods' that specifies the optimization algorithm (a hybrid of gradient descent using the Adam optimizer with learning rate 0.01 and basin-hopping for escape from local minima), control discretization (piecewise-constant protocols with 200 time steps), convergence tolerances (relative change in the cost functional below 10^{-8} for 50 consecutive iterations), and multi-start diagnostics (200 independent random initializations per parameter value, with the lowest-cost solution retained and statistics on cost variance reported). These additions confirm that the reported discontinuous switches are robust across restarts. revision: yes

  2. Referee: [Results for the coherent spin-boson model] For the spin-boson model, the discontinuous switch between locally optimal solutions is presented as a key finding, but without comparison to known analytic limits at the slow- and rapid-driving boundaries or exhaustive sampling statistics, it is unclear whether the discontinuity survives a more thorough search of the control landscape.

    Authors: We have added explicit comparisons to the analytic limits in the revised Figure 2: as the driving time approaches the slow-driving regime, the optimal protocols converge to the linear-response prediction, while in the rapid-driving limit they recover the sudden-quench form. These limits are now overlaid on the numerical results. We have also included exhaustive sampling statistics in a new appendix (500 random starts plus 1000 basin-hopping iterations per weight value), which show that the discontinuity in the optimal protocol persists and is not an artifact of incomplete search. We acknowledge that, for non-convex landscapes, numerical methods cannot provide an absolute mathematical guarantee of global optimality. revision: partial

standing simulated objections not resolved
  • Absolute mathematical proof of global optimality for the numerically obtained protocols in a non-convex control landscape

Circularity Check

0 steps flagged

No significant circularity; results from direct numerical optimization

full rationale

The paper derives its claims about optimal protocols in the spin-boson and quantum-dot models exclusively via numerical minimization of a weighted cost functional over driving protocols. No algebraic derivation chain is presented that reduces a claimed prediction back to a fitted input, self-defined quantity, or self-citation load-bearing premise. The enumerated circularity patterns (self-definitional, fitted-input-called-prediction, uniqueness-imported-from-authors, etc.) are absent; the central results are computational outputs rather than tautological re-expressions of inputs. This is the normal, self-contained case for a numerical optimal-control study.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The work rests on the assumption that the two model Hamiltonians adequately represent generic open quantum systems and that the numerical optimizer finds protocols that are at least locally optimal for the chosen cost functions.

free parameters (1)
  • relative weight between dissipation and fluctuations
    Varied continuously to observe the discontinuous switch in optimal protocol.
axioms (1)
  • domain assumption The coherent spin-boson and single-level quantum-dot models capture the essential open-system dynamics relevant to the optimization problem.

pith-pipeline@v0.9.0 · 5417 in / 1094 out tokens · 21757 ms · 2026-05-16T08:42:24.357251+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages · 2 internal anchors

  1. [1]

    (A7)–(A9)

    For a given trial controlv(t), propagate the state variablesu(t),y(t), andx(t) forward in time using Eqs. (A7)–(A9)

    work page

  2. [2]

    (A10)-(A12) att=T

    Set the terminal values of the adjoint variables accord- ing to Eqs. (A10)-(A12) att=T

    work page

  3. [3]

    (A13) and (A14)-(A16)

    Propagate the adjoint variables backward in time us- ing the adjoint equations and the Pontryagin Hamil- tonianH pmp(t) in Eqs. (A13) and (A14)-(A16)

    work page

  4. [4]

    Computeδv(t) =∂H pmp/∂v(t) and update the control according tov(t)→v(t)−η δv(t), with a suitable step sizeη >0

    work page

  5. [5]

    Repeat steps 1–4 until convergence of the costJ. In step 1, the dynamics of the state variables is then expressed as ˙u(t) =v(t),(A7) ˙y(t) =A(u(t))y(t) +v(t)s x(t), u(t) ,(A8) ˙x(t) =A(u(t))x(t).(A9) Next, in step 2, we associate adjoint variablesp(t), Λ(t), and Π(t) withu(t),y(t), andx(t), respectively. The terminal values of the adjoint variables are d...

    work page

  6. [6]

    For a given trial controlu(t), propagate the state vari- ableq(t) forward in time by solving ˙q(t) =f q(t), u(t) , q(0) =q 0.(D8)

    work page

  7. [7]

    Set the terminal value of the adjoint variable att=T as p(T) = ∂d(T, q) ∂q ⊤ q=q(T) .(D9)

    work page

  8. [8]

    Propagate the adjoint variablep(t) backward in time using the adjoint equation ˙p(t) =− ∂H ∂q q(t), p(t), u(t) ⊤ .(D10)

    work page

  9. [9]

    Compute the gradient g(t) = ∂H ∂u q(t), p(t), u(t) , and update the control according to u(t)→u(t)−η g(t), with a suitable step sizeη >0

    work page

  10. [10]

    Appendix E: calculation of∂ uA(u) We analytically calculate∂ uA(u) beforehand for the GRAPE algorithm

    Repeat steps 1–4 until convergence of the cost. Appendix E: calculation of∂ uA(u) We analytically calculate∂ uA(u) beforehand for the GRAPE algorithm. ∂uA(u) =∂ uC(u) +∂ u R(−θ)L(E)R(θ) =∂ uC(u) + ∂uR(−θ) L(E)R(θ) +R(−θ) ∂uL(E) R(θ) +R(−θ)L(E) ∂uR(θ) . (E1) As defined in Eqs.(14), dE du = u E , d du cosθ= ∆2 E3 , d du sinθ=− u∆ E3 . (E2) The coherent cont...

    work page

  11. [11]

    Stochastic thermodynamics, fluctuation the- orems and molecular machines,

    U. Seifert, “Stochastic thermodynamics, fluctuation the- orems and molecular machines,” Rep. Prog. Phys.75, 126001 (2012)

    work page 2012

  12. [12]

    Optimal finite-time pro- cesses in stochastic thermodynamics,

    T. Schmiedl and U. Seifert, “Optimal finite-time pro- cesses in stochastic thermodynamics,” Phys. Rev. Lett. 98, 108301 (2007)

    work page 2007

  13. [13]

    Computing the optimal protocol for finite-time processes in stochastic thermodynamics,

    H. Then and A. Engel, “Computing the optimal protocol for finite-time processes in stochastic thermodynamics,” Phys. Rev. E77, 041105 (2008)

    work page 2008

  14. [14]

    Steps min- imize dissipation in rapidly driven stochastic systems,

    S. Blaber, M. D. Louwerse, and D. A. Sivak, “Steps min- imize dissipation in rapidly driven stochastic systems,” Phys. Rev. E104, L022101 (2021)

    work page 2021

  15. [15]

    Thermodynamic metrics and optimal paths,

    D. A. Sivak and G. E. Crooks, “Thermodynamic metrics and optimal paths,” Phys. Rev. Lett.108, 190602 (2012)

    work page 2012

  16. [16]

    Phase transition in pro- tocols minimizing work fluctuations,

    A. P. Solon and J. M. Horowitz, “Phase transition in pro- tocols minimizing work fluctuations,” Phys. Rev. Lett. 120, 180605 (2018)

    work page 2018

  17. [17]

    Binder, L

    F. Binder, L. A. Correa, C. Gogolin, J. Anders, and G. Adesso (eds.),Thermodynamics in the Quantum Regime: Fundamental Aspects and New Directions, Springer, Cham (2018)

    work page 2018

  18. [18]

    A Quantum Fluctuation Theorem

    J. Kurchan, “A quantum fluctuation theorem,” arXiv:cond-mat/0007360 (2000)

    work page internal anchor Pith review Pith/arXiv arXiv 2000

  19. [19]

    Jarzynski Relations for Quantum Systems and Some Applications

    H. Tasaki, “Jarzynski relations for quantum systems and some applications,” arXiv:cond-mat/0009244 (2000)

    work page internal anchor Pith review Pith/arXiv arXiv 2000

  20. [20]

    Fluctuation the- orems: Work is not an observable,

    P. Talkner, E. Lutz, and P. H¨ anggi, “Fluctuation the- orems: Work is not an observable,” Phys. Rev. E75, 050102(R) (2007). 11

    work page 2007

  21. [21]

    Colloquium: Quantum fluctuation relations: Foundations and appli- cations,

    M. Campisi, P. H¨ anggi, and P. Talkner, “Colloquium: Quantum fluctuation relations: Foundations and appli- cations,” Rev. Mod. Phys.83, 771 (2011)

    work page 2011

  22. [22]

    Nonequi- librium fluctuations, fluctuation theorems, and count- ing statistics in quantum systems,

    M. Esposito, U. Harbola, and S. Mukamel, “Nonequi- librium fluctuations, fluctuation theorems, and count- ing statistics in quantum systems,” Rev. Mod. Phys.81, 1665 (2009)

    work page 2009

  23. [23]

    No-Go Theorem for the Char- acterization of Work Fluctuations in Coherent Quantum Systems,

    M. Perarnau-Llobet, E. B¨ aumer, K. V. Hovhannisyan, M. Huber, and A. Ac´ ın, “No-Go Theorem for the Char- acterization of Work Fluctuations in Coherent Quantum Systems,” Phys. Rev. Lett. 118, 070601 (2017)

    work page 2017

  24. [24]

    Quantum extension of the Jarzynski re- lation: Analogy with stochastic dephasing,

    S. Mukamel, “Quantum extension of the Jarzynski re- lation: Analogy with stochastic dephasing,” Phys. Rev. Lett.90, 170604 (2003)

    work page 2003

  25. [25]

    Unified treatment of the quantum fluctua- tion theorem and Jarzynski equality in terms of micro- scopic reversibility,

    T. Monnai, “Unified treatment of the quantum fluctua- tion theorem and Jarzynski equality in terms of micro- scopic reversibility,” Phys. Rev. E72, 027102 (2005)

    work page 2005

  26. [26]

    Moments of work in the two-point mea- surement protocol for a driven open quantum system,

    S. Suomela, P. Solinas, J. P. Pekola, J. Ankerhold, and T. Ala-Nissila, “Moments of work in the two-point mea- surement protocol for a driven open quantum system,” Phys. Rev. B90, 094304 (2014)

    work page 2014

  27. [27]

    Work fluctuations in slow processes: Quantum signatures and optimal control,

    H. J. D. Miller, M. Scandi, J. Anders, and M. Perarnau- Llobet, “Work fluctuations in slow processes: Quantum signatures and optimal control,” Phys. Rev. Lett.123, 230603 (2019)

    work page 2019

  28. [28]

    Optimal control of dissipation and work fluctuations for rapidly driven systems,

    A. Rolandi, M. Perarnau-Llobet, and H. J. D. Miller, “Optimal control of dissipation and work fluctuations for rapidly driven systems,” New J. Phys.25, 073005 (2023)

    work page 2023

  29. [29]

    Breuer and F

    H.-P. Breuer and F. Petruccione,The Theory of Open Quantum Systems, Oxford University Press, Oxford (2002)

    work page 2002

  30. [30]

    Dynam- ics of the dissipative two-state system,

    A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, “Dynam- ics of the dissipative two-state system,” Rev. Mod. Phys. 59, 1 (1987)

    work page 1987

  31. [31]

    Completely positive dynamical semigroups ofN-level systems,

    V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, “Completely positive dynamical semigroups ofN-level systems,”J. Math. Phys.17, 821–825 (1976)

    work page 1976

  32. [32]

    On the generators of quantum dynamical semigroups,

    G. Lindblad, “On the generators of quantum dynamical semigroups,”Commun. Math. Phys.48, 119–130 (1976)

    work page 1976

  33. [33]

    A short introduction to the Lindblad mas- ter equation,

    D. Manzano, “A short introduction to the Lindblad mas- ter equation,” AIP Adv.10, 025106 (2020)

    work page 2020

  34. [34]

    Rivas and S

    A. Rivas and S. F. Huelga,Open Quantum Systems: An Introduction, Springer, Berlin (2012)

    work page 2012

  35. [35]

    Thermodynamic length in open quantum systems,

    M. Scandi and M. Perarnau-Llobet, “Thermodynamic length in open quantum systems,” Quantum3, 197 (2019)

    work page 2019

  36. [36]

    Op- timal thermodynamic control in open quantum systems,

    V. Cavina, A. Mari, A. Carlini, and V. Giovannetti, “Op- timal thermodynamic control in open quantum systems,” Phys. Rev. A98, 012139 (2018)

    work page 2018

  37. [37]

    Introduction to the Pontryagin maximum principle for quantum optimal control,

    U. Boscain, M. Sigalotti, and D. Sugny, “Introduction to the Pontryagin maximum principle for quantum optimal control,” PRX Quantum2, 030203 (2021)

    work page 2021

  38. [38]

    Optimal control of coupled spin dy- namics: Design of NMR pulse sequences by gradient as- cent algorithms,

    N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbr¨ uggen, and S. J. Glaser, “Optimal control of coupled spin dy- namics: Design of NMR pulse sequences by gradient as- cent algorithms,” J. Magn. Reson.172, 296–305 (2005)

    work page 2005

  39. [39]

    Quantum master equation for electron transport through quantum dots and single molecules,

    U. Harbola, M. Esposito, and S. Mukamel, “Quantum master equation for electron transport through quantum dots and single molecules,” Phys. Rev. B74, 235309 (2006)

    work page 2006

  40. [40]

    Fluc- tuations of work from quantum subensembles: The case against quantum work-fluctuation theorems,

    A. E. Allahverdyan and T. M. Nieuwenhuizen, “Fluc- tuations of work from quantum subensembles: The case against quantum work-fluctuation theorems,” Phys. Rev. E71, 066102 (2005)

    work page 2005

  41. [41]

    Quan- tum adiabatic Markovian master equations,

    T. Albash, S. Boixo, D. A. Lidar, and P. Zanardi, “Quan- tum adiabatic Markovian master equations,” New J. Phys.14, 123016 (2012)

    work page 2012

  42. [42]

    Pareto-optimal cycles for power, efficiency and fluctuations of quantum heat engines using reinforcement learning,

    P. A. Erdman, A. Rolandi, P. Abiuso, M. Perarnau- Llobet, and F. No´ e, “Pareto-optimal cycles for power, efficiency and fluctuations of quantum heat engines using reinforcement learning,” Phys. Rev. Research5, L022017 (2023)

    work page 2023