pith. sign in

arxiv: 2604.18212 · v2 · pith:IMB5ALD2new · submitted 2026-04-20 · 🪐 quant-ph · cond-mat.mes-hall

Davies-Morris-Shore Framework for Multilevel Quantum Batteries: Dark and Funnel States in Interacting Qutrit Systems

Yi\u{g}it Per\c{c}in , \"Ozg\"ur E. M\"ustecapl{\i}o\u{g}lu This is my paper

Pith reviewed 2026-05-21 00:51 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.mes-hall
keywords quantum batteriesdark statesfunnel statesqutritsopen quantum systemsDavies master equationMorris-Shore decompositiondissipation
0
0 comments X

The pith

Dark and funnel states in two interacting qutrits confine decay to protected manifolds for stable multilevel quantum batteries.

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

The paper develops a systematic way to locate long-lived excited states in open multilevel systems by combining the Davies master equation with a Morris-Shore decomposition of the dissipative couplings. For the minimal case of two qutrits that interact via exchange and share a common bath, the authors analytically build dark, bright, and funnel states whose decay is channeled exclusively into stable manifolds. Robustness of these states is controlled by the ratio of interaction strength to anharmonicity, and numerical checks confirm extended energy storage times under realistic loss. The multilevel ladder plus exchange interaction opens storage possibilities that qubit ensembles cannot reach. High-energy funnel states are singled out as practical targets because their decay routes are structured and point toward protection.

Core claim

In the two-qutrit model the Davies-Morris-Shore decomposition isolates dissipative blocks so that certain excited states decay only into protected manifolds; the lifetime of these states is set by the interaction-to-anharmonicity ratio, and numerical evolution shows that the resulting energetic storage persists far longer than in the absence of the constructed protection.

What carries the argument

Morris-Shore decomposition of dissipative coupling blocks, which partitions the open-system Liouvillian into blocks whose decay pathways lead only into protected subspaces.

If this is right

  • Multilevel ladder structure and exchange interactions produce storage states unavailable in qubit models.
  • High-energy funnel states direct decay toward protected manifolds and therefore serve as concrete design targets.
  • The interaction-to-anharmonicity ratio supplies a quantitative knob for engineering robustness.
  • Knowledge of the structured decay pathways supplies a basis for future protection and control protocols.

Where Pith is reading between the lines

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

  • The same decomposition could be applied to chains or arrays of qutrits to identify scalable protected subspaces.
  • The framework may extend to other anharmonic multilevel systems such as superconducting transmons or trapped ions.
  • Experimental sweeps of the interaction-to-anharmonicity ratio in present-day hardware would directly test the predicted lifetime scaling.

Load-bearing premise

The model assumes exactly two qutrits coupled to a common bath so that the Davies master equation plus Morris-Shore decomposition cleanly separates all dissipative channels without extra decoherence or higher-order terms.

What would settle it

Direct measurement of the funnel-state lifetime in a two-qutrit device showing that lifetime remains short even when interaction strength greatly exceeds anharmonicity would contradict the derived robustness conditions.

Figures

Figures reproduced from arXiv: 2604.18212 by \"Ozg\"ur E. M\"ustecapl{\i}o\u{g}lu, Yi\u{g}it Per\c{c}in.

Figure 1
Figure 1. Figure 1: FIG. 1. Dressed-state decay structure of the two-qutrit sys [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Energy storage dynamics and dark-state overlap fol [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Comparison between qubit and qutrit quantum bat [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Robustness of the dark manifold as a function of [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Decay dynamics of different initial states with a re [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Active stabilization via dynamical decoupling. We [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
read the original abstract

Dark and subradiant states have emerged as a promising resource for stabilizing open quantum batteries against dissipation, but existing studies are largely limited to qubit ensembles and symmetry-based constructions. Here we introduce a systematic, thermodynamically consistent framework for identifying long-lived energy storage states in interacting multilevel quantum batteries, combining the Davies master equation with a Morris-Shore (MS)-type decomposition of dissipative coupling blocks. Focusing on a minimal model of two interacting qutrits coupled to a common bath, we analytically construct dark, bright, and funnel states-excited states that decay exclusively into protected manifolds. We also derive quantitative robustness conditions governed by the ratio of interaction strength to anharmonicity. We show that multilevel ladder structure and exchange interactions enable energetic storage states beyond the qubit case. Numerical simulations confirm that these states exhibit long-lived energy storage under realistic dissipation. Finally, we show that high-energy funnel states provide a natural design target for multilevel quantum batteries, as their decay pathways are highly structured and directed toward protected manifolds. Knowledge of these pathways offers a principled basis for developing future protection and control strategies in superconducting multilevel platforms.

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

1 major / 2 minor

Summary. The paper introduces a Davies-Morris-Shore framework for multilevel quantum batteries in a minimal model of two exchange-coupled qutrits coupled to a common bath. It analytically constructs dark, bright, and funnel states (excited states decaying exclusively into protected manifolds) by combining the Davies master equation with a Morris-Shore-type decomposition of dissipative blocks, derives quantitative robustness conditions controlled by the ratio of interaction strength J to anharmonicity Δ, and reports numerical simulations showing long-lived energy storage under realistic dissipation. The work emphasizes that multilevel ladder structure and exchange interactions enable storage states beyond the qubit case, with high-energy funnel states as natural design targets.

Significance. If the central analytical constructions and robustness conditions hold exactly, the result would be significant for open quantum batteries: it supplies a thermodynamically consistent, systematic method to identify protected multilevel states rather than relying on symmetry arguments alone, and the funnel-state concept offers a concrete design principle for superconducting platforms. The combination of Davies generator with MS decomposition and the explicit derivation of J/Δ-governed conditions are strengths; numerical confirmation of long storage times adds practical support.

major comments (1)
  1. [Morris-Shore decomposition and analytical construction of dark/funnel states] The Morris-Shore decomposition section and the derivation of funnel/dark states: the central claim treats the block diagonalization of the collective dissipator as cleanly separating dark/funnel subspaces from the bright manifold. However, finite anharmonicity Δ introduces off-diagonal matrix elements between these subspaces that are suppressed only by the ratio J/Δ. If the analytic construction sets these elements to exactly zero rather than treating them perturbatively, the derived quantitative robustness conditions become approximate; any residual leakage rate would impose an upper bound on storage time not captured by the reported numerics. Please provide the explicit matrix elements of the dissipator in the MS basis and state whether the separation is exact or holds only to leading order in J/Δ.
minor comments (2)
  1. [Numerical simulations] Numerical simulations section: the abstract states that simulations confirm long-lived storage, but the manuscript lacks explicit error bars, convergence checks with respect to Hilbert-space truncation, and a clear statement of how parameter regimes were chosen to avoid post-hoc selection. Adding these details would strengthen reproducibility.
  2. [Model and Hamiltonian] Notation and definitions: the robustness condition is stated in terms of the ratio of interaction strength to anharmonicity, but the precise definition of the anharmonicity parameter Δ (relative to the qutrit ladder frequencies) should be given explicitly in an equation early in the text.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address the major comment below and describe the revisions we will implement.

read point-by-point responses

  1. Referee: [Morris-Shore decomposition and analytical construction of dark/funnel states] The Morris-Shore decomposition section and the derivation of funnel/dark states: the central claim treats the block diagonalization of the collective dissipator as cleanly separating dark/funnel subspaces from the bright manifold. However, finite anharmonicity Δ introduces off-diagonal matrix elements between these subspaces that are suppressed only by the ratio J/Δ. If the analytic construction sets these elements to exactly zero rather than treating them perturbatively, the derived quantitative robustness conditions become approximate; any residual leakage rate would impose an upper bound on storage time not captured by the reported numerics. Please provide the explicit matrix elements of the dissipator in the MS basis and state whether the separation is exact or holds only to leading order in J/Δ.

    Authors: We thank the referee for this precise observation. In the manuscript the Morris-Shore transformation is applied to the collective dissipator blocks of the Davies generator. The resulting block-diagonal structure is exact only in the limit J/Δ → ∞; for finite anharmonicity the transformation leaves residual off-diagonal couplings whose magnitude is O(Δ/J) (or higher-order in the perturbative expansion). The quantitative robustness conditions we derive are precisely the statements that these couplings remain negligible inside the regime J ≫ Δ, so that leakage out of the dark and funnel manifolds stays perturbatively small. The numerical simulations solve the full master equation and therefore already include any residual leakage; the reported storage times are consistent with the analytic bounds. To make this explicit we will add, in a revised Section III and a new appendix, the full matrix representation of the dissipator in the Morris-Shore basis, together with the scaling of the off-diagonal blocks. We will also insert a short paragraph clarifying that the separation is exact only to leading order in the J/Δ expansion. revision: yes

Circularity Check

0 steps flagged

Derivation chain from Davies master equation and Morris-Shore decomposition remains self-contained with no reduction to inputs by construction.

full rationale

The paper derives dark, bright, and funnel states analytically by applying the Davies master equation to the Hamiltonian of two exchange-coupled qutrits and performing a Morris-Shore-type block decomposition of the collective dissipator. The resulting subspaces and robustness conditions (governed by the J/Δ ratio) follow directly from the structure of the Lindblad operators and the anharmonic ladder without any fitted parameters being relabeled as predictions or any central premise collapsing to a self-citation. Numerical simulations serve as independent verification rather than input. No self-definitional, fitted-input, or load-bearing self-citation steps appear in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based on abstract only; no explicit free parameters, axioms, or invented entities are detailed. The framework relies on standard open-quantum-system assumptions such as the validity of the Davies master equation for the bath coupling.

pith-pipeline@v0.9.0 · 5754 in / 1219 out tokens · 31466 ms · 2026-05-21T00:51:58.322002+00:00 · methodology

Review history (2 revisions) →

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

39 extracted references · 39 canonical work pages · 1 internal anchor

  1. [1]

    Therefore, the spectator state is an exact dark state only in the harmonic limitα= 0

    Restricting the Hamiltonian (3) to theN= 2 manifold yields the Hamiltonian (21), seeking an eigenstate of the symmetric form |ψ⟩ ∝(1, x,1),(34) one finds thatxsatisfies the quadratic equation √ 2J x2 −αx−2 √ 2J= 0.(35) For weak anharmonicityα≪J, the physically relevant solution is x≃ − √ 2 + √ 2 4 α J +O α2 J2 .(36) Equation (36) quantifies the deviation ...

    work page

  2. [2]

    (8) from which one constructs the decay graph between eigenstates

    Identification of funnel states Let{|ϕ k⟩}denote eigenstates of the interacting Hamil- tonianH S, HS |ϕk⟩=E k |ϕk⟩.(37) The dissipative structure is determined by the Davies jump operators, which were given in Eq. (8) from which one constructs the decay graph between eigenstates. A state|ϕ k⟩is classified as •dark ifA(Ω)|ϕ k⟩= 0 for all Ω, •funnel if the ...

    work page

  3. [3]

    Let |ϕF ⟩denote a selected funnel state

    Targeted funnel-state charging Once the funnel states are identified, the charging problem reduces to preparing a chosen funnel state. Let |ϕF ⟩denote a selected funnel state. The objective is to design a control Hamiltonian Hch(t) = X ℓ uℓ(t)Hℓ (38) such that U(T)|ψ 0⟩ ≈ |ϕ F ⟩, U(T) =Texp " −i Z T 0 (HS +H ch(t))dt # ,(39) where|ψ 0⟩is the initial groun...

    work page

  4. [4]

    Importantly, this comparison is per- formedafteridentifying funnel states, rather than search- ing over arbitrary eigenstates

    Selection among multiple funnel states If multiple funnel states exist, one may compare them using their stored energy EF =⟨ϕ F |H S |ϕF ⟩,(42) and an effective decay rate ΓF = X Ω>0 γ(Ω)⟨ϕ F |A †(Ω)A(Ω)|ϕ F ⟩.(43) The optimal target is typically the highest-energy fun- nel state whose decay graph terminates in the protected dark manifold. Importantly, th...

    work page

  5. [5]

    Scalable recipe The proposed scalable protocol for efficient charging of quantum batteries follows these steps: First, the sys- tem Hamiltonian (HS) is diagonalized to generate Davies jump operators decomposed into Bohr frequencies. Next, the SVD-based Morris-Shore transform is applied to the dissipative matrix, which represents the transitions be- tween ...

    work page

  6. [6]

    dark-state preparation

    Two-qutrit example For the two-qutrit model discussed above, the an- tisymmetric two-excitation state|F 2⟩acts as the fun- nel state, while the antisymmetric single-excitation state |D1⟩forms the protected dark manifold. The charging protocol therefore aims to prepare|F 2⟩, after which dis- sipation loads the long-lived storage state|D 1⟩. In our approach...

    work page

  7. [7]

    K. V. Hovhannisyan, M. Perarnau-Llobet, M. Huber, and A. Ac´ ın, Phys. Rev. Lett.111, 240401 (2013)

    work page 2013

  8. [8]

    Alicki and M

    R. Alicki and M. Fannes, Phys. Rev. E87, 042123 (2013)

    work page 2013

  9. [9]

    Campaioli, S

    F. Campaioli, S. Gherardini, J. Q. Quach, M. Polini, and G. M. Andolina, Rev. Mod. Phys.96, 031001 (2024)

    work page 2024

  10. [10]

    F. C. Binder, S. Vinjanampathy, K. Modi, and J. Goold, Phys. Rev. E91, 032119 (2015)

    work page 2015

  11. [11]

    Campaioli, F

    F. Campaioli, F. A. Pollock, F. C. Binder, L. C´ eleri, J. Goold, and K. Modi, Phys. Rev. Lett.118, 150601 (2017)

    work page 2017

  12. [12]

    Tejero, J

    A. Tejero, J. Thingna, and D. Manzano, J. Phys. Chem. C125, 7518 (2021)

    work page 2021

  13. [13]

    J. Q. Quach and W. J. Munro, Phys. Rev. Appl.14, 024092 (2020)

    work page 2020

  14. [14]

    J. Liu, D. Segal, and G. Hanna, J. Phys. Chem. C123, 18303 (2019)

    work page 2019

  15. [15]

    Finkelstein-Shapiro, S

    D. Finkelstein-Shapiro, S. Felicetti, T. Hansen, T. o. Pul- lerits, and A. Keller, Phys. Rev. A99, 053829 (2019)

    work page 2019

  16. [16]

    Dou and F.-M

    F.-Q. Dou and F.-M. Yang, Phys. Rev. A107, 023725 (2023)

    work page 2023

  17. [17]

    Yang and F.-Q

    F.-M. Yang and F.-Q. Dou, Phys. Rev. A109, 062432 (2024)

    work page 2024

  18. [18]

    Gemme, M

    G. Gemme, M. Grossi, S. Vallecorsa, M. Sassetti, and D. Ferraro, Phys. Rev. Res.6, 023091 (2024)

    work page 2024

  19. [19]

    A. A. Rangelov, N. V. Vitanov, and B. W. Shore, Phys. Rev. A74, 053402 (2006)

    work page 2006

  20. [20]

    P.-Y. Sun, H. Zhou, and F.-Q. Dou, (2025), 10.48550/arXiv.2412.01442

    work page doi:10.48550/arxiv.2412.01442 2025

  21. [21]

    J. R. Morris and B. W. Shore, Phys. Rev. A27, 906 (1983)

    work page 1983

  22. [22]

    J. Koch, T. M. Yu, J. M. Gambetta, A. A. Houck, D. I. Schuster, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. A76, 042319 (2007)

    work page 2007

  23. [23]

    M. J. Peterer, S. J. Bader, X. Jin, F. Yan, A. Kamal, T. J. Gudmundsen, P. J. Leek, T. P. Orlando, W. D. Oliver, and S. Gustavsson, Phys. Rev. Lett.114, 010501 (2015)

    work page 2015

  24. [24]

    Zeytino˘ glu, M

    S. Zeytino˘ glu, M. Pechal, S. Berger, A. A. Abdumalikov, A. Wallraff, and S. Filipp, Phys. Rev. A91, 043846 (2015)

    work page 2015

  25. [25]

    E. B. Davies, Com. in Math. Phys.39, 91 (1974)

    work page 1974

  26. [26]

    Breuer and F

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

    work page 2002

  27. [27]

    Saadati-Niari, An

    M. Saadati-Niari, An. of Phys.372, 138 (2016)

    work page 2016

  28. [28]

    H. Kim, Y. Song, H.-g. Lee, and J. Ahn, Phys. Rev. A 91, 053421 (2015)

    work page 2015

  29. [29]

    K. N. Zlatanov, G. S. Vasilev, and N. V. Vitanov, Phys. Rev. A102, 063113 (2020)

    work page 2020

  30. [30]

    C. Brif, R. Chakrabarti, and H. Rabitz, New J. of Phys. 12, 075008 (2010)

    work page 2010

  31. [31]

    Khaneja, T

    N. Khaneja, T. Reiss, C. Luy, and S. J. Glaser, J. of Mag. Res.172, 296 (2005)

    work page 2005

  32. [32]

    Massa, F

    N. Massa, F. Cavaliere, and D. Ferraro, Batteries11, 240 (2025)

    work page 2025

  33. [33]

    Elghaayda, A

    S. Elghaayda, A. Ali, and S. Haddadi, Adv. Quantum Technol. , 2400651 (2024)

    work page 2024

  34. [34]

    R. H. Dicke, Phys. Rev.93, 99 (1954)

    work page 1954

  35. [35]

    Viola and S

    L. Viola and S. Lloyd, Phys. Rev. A58, 2733 (1998)

    work page 1998

  36. [36]

    B.-B. Liu, R. Chen, J.-L. Wu, G. Chen, and S.-L. Su, (2025), 10.48550/arXiv.2603.23009

    work page doi:10.48550/arxiv.2603.23009 2025

  37. [37]

    Dynamically stabilized decoherence-free states in non-Markovian open fermionic systems

    H.-N. Xiong, W.-M. Zhang, M. W.-Y. Tu, and D. Braun, (2012), 10.48550/arXiv.1207.1538

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1207.1538 2012

  38. [38]

    X. Wang, M. Byrd, and K. Jacobs, Phys. Rev. A87, 012338 (2013)

    work page 2013

  39. [39]

    Evangelakos, E

    V. Evangelakos, E. Paspalakis, and D. Stefanatos, Sci. Rep.15, 45626 (2025)

    work page 2025