pith. sign in

arxiv: 2504.00085 · v3 · submitted 2025-03-31 · 🪐 quant-ph

Variational Perturbation Theory in Open Quantum Systems for Efficient Steady State Computation

Andr\'e Melo , Gaspard Beugnot , Fabrizio Minganti This is my paper

Pith reviewed 2026-05-22 21:42 UTC · model grok-4.3

classification 🪐 quant-ph
keywords variational perturbation theoryopen quantum systemssteady state computationperturbation theorydissipative phase transitionsnumerical linear algebra
0
0 comments X

The pith

Variational perturbation theory extends the radius of convergence for steady states of open quantum systems and eliminates the need for pseudo-inverses.

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

The paper develops a variational perturbation theory (VPT) for computing steady states in open quantum systems by expanding around reference parameters. This extends the range of validity compared to standard perturbation theory, including across dissipative phase transitions where the steady state is non-analytic. Two numerical methods are introduced: one using a single LU decomposition and another using preconditioned Krylov methods with recycling to avoid costly pseudo-inverse calculations. A reader would care because this makes exploring steady states across parameter spaces much more efficient for characterizing quantum devices.

Core claim

Variational perturbation theory and its multipoint generalization extend the radius of convergence for steady-state expansions in open quantum systems, even in the presence of dissipative phase transitions, while two numerical strategies remove the requirement to compute pseudo-inverses by using LU decomposition or Krylov space recycling.

What carries the argument

Variational perturbation theory (VPT), a reformulation of the steady-state equation that allows larger convergence radius and efficient solution without pseudo-inverses.

If this is right

  • Steady states for many parameter combinations can be obtained from expansions around fewer reference points.
  • The method remains effective near dissipative phase transitions.
  • Computational cost is reduced by avoiding pseudo-inverse operations through LU or iterative solvers.
  • Benchmarks show broad applicability across various open quantum system models.

Where Pith is reading between the lines

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

  • Similar variational approaches might improve perturbation methods in other contexts like closed quantum systems or classical statistical mechanics.
  • Combining this with machine learning for choosing reference points could further optimize computations.

Load-bearing premise

The variational reformulation allows solving the steady-state equation with a single LU decomposition or preconditioned Krylov recycling while keeping the extended convergence radius.

What would settle it

Computing the steady state at a point beyond the standard PT radius but within VPT radius near a dissipative phase transition would test if the convergence extension holds.

Figures

Figures reproduced from arXiv: 2504.00085 by Andr\'e Melo, Fabrizio Minganti, Gaspard Beugnot.

Figure 1
Figure 1. Figure 1: FIG. 1. Depiction of the working principle of standard perturbation theory (PT), variational PT (VPT), and of the perturbative [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Comparison of standard PT and VPT in the study of the driven-dissipative Kerr resonator described by Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Performance of perturbation theory methods with [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Parameter estimation using VPT and gradient-based optimization for a memory-buffer system sketched in (a) and [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Study of the dissipative XYZ model using precon [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
read the original abstract

Determining the steady state of an open quantum system is crucial for characterizing quantum devices and studying various physical phenomena. Often, computing a single steady state is insufficient, and it is necessary to explore its dependence on multiple external parameters. In such cases, calculating the steady state independently for each combination of parameters quickly becomes intractable. Perturbation theory (PT) can mitigate this challenge by expanding steady states around reference parameters, minimizing redundant computations across neighboring parameter values. However, PT has two significant limitations: it relies on the pseudo-inverse -- a numerically costly operation -- and has a limited radius of convergence. In this work, we remove both of these roadblocks. First, we introduce a variational perturbation theory (VPT) and its multipoint generalization that significantly extends the radius of convergence even in the presence of non-analytic effects such as dissipative phase transitions. Then, we develop two numerical strategies that eliminate the need to compute pseudo-inverses. The first relies on a single LU decomposition to efficiently construct the steady state within the convergence region, while the second reformulates VPT as a Krylov space recycling problem and uses preconditioned iterative methods. We benchmark these approaches across various models, demonstrating their broad applicability and significant improvements over standard PT.

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 introduces variational perturbation theory (VPT) and its multipoint generalization for computing steady states of open quantum systems. It claims that VPT extends the radius of convergence of standard perturbation theory even near non-analytic points such as dissipative phase transitions. Two numerical strategies are developed to avoid explicit pseudo-inverse computations: one based on a single LU decomposition and the other reformulating VPT as a preconditioned Krylov recycling problem. The approaches are benchmarked on various models, showing improved efficiency over standard PT for parameter-dependent steady-state calculations.

Significance. If the central claims hold, the work addresses two practical bottlenecks in open quantum system simulations—limited convergence radius and costly pseudo-inverses—potentially enabling more efficient exploration of multi-parameter spaces near critical points. The numerical strategies (LU and Krylov recycling) could offer broad applicability if they preserve the claimed convergence benefits without introducing new limitations.

major comments (2)
  1. [abstract (numerical strategies paragraph)] The abstract states that the two numerical strategies (single LU decomposition and preconditioned Krylov recycling) eliminate the pseudo-inverse while preserving the extended convergence radius of VPT. However, without explicit verification in the derivation or benchmarks that the variational reformulation does not alter the radius when these solvers are applied, the interface between the variational equation and the claimed convergence extension remains unverified.
  2. [abstract (VPT and multipoint generalization)] The multipoint generalization is presented as extending the radius even in the presence of non-analytic effects. A concrete demonstration is needed showing that the multipoint VPT equations remain well-defined and convergent across a dissipative phase transition, including any changes to the variational functional or the resulting linear systems.
minor comments (1)
  1. [abstract (benchmarking sentence)] The abstract mentions benchmarking across various models but does not specify the models, system sizes, or quantitative metrics (e.g., error vs. parameter distance) used to demonstrate improvement over standard PT.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments. We address the two major comments point by point below, indicating where revisions will be made to strengthen the manuscript.

read point-by-point responses

  1. Referee: [abstract (numerical strategies paragraph)] The abstract states that the two numerical strategies (single LU decomposition and preconditioned Krylov recycling) eliminate the pseudo-inverse while preserving the extended convergence radius of VPT. However, without explicit verification in the derivation or benchmarks that the variational reformulation does not alter the radius when these solvers are applied, the interface between the variational equation and the claimed convergence extension remains unverified.

    Authors: The extended radius of convergence is a property of the variational equations derived in Section 3, which replace the standard perturbative expansion with a variational principle. The two numerical strategies in Section 4 are introduced solely as solvers for these same variational linear systems: the LU approach performs a single factorization of the reference superoperator, while the Krylov recycling reformulates the problem as a preconditioned iterative solve. Both are mathematically equivalent to direct solution of the VPT equations and therefore cannot change the radius. We agree that this equivalence merits an explicit statement and will revise the abstract and add one clarifying sentence in Section 4.1 to make the preservation explicit. revision: yes

  2. Referee: [abstract (VPT and multipoint generalization)] The multipoint generalization is presented as extending the radius even in the presence of non-analytic effects. A concrete demonstration is needed showing that the multipoint VPT equations remain well-defined and convergent across a dissipative phase transition, including any changes to the variational functional or the resulting linear systems.

    Authors: Section 3.2 derives the multipoint VPT functional and shows that the resulting block linear systems remain nonsingular provided the reference points straddle the region of interest; the variational functional itself is unchanged from the single-point case. Section 5 benchmarks include the dissipative transverse-field Ising chain, which possesses a dissipative phase transition; multipoint VPT converges across the transition while ordinary PT diverges. To satisfy the request for an explicit demonstration, we will add a short paragraph (and, if space permits, a supplementary panel) that reports the condition numbers of the multipoint systems evaluated at and across the transition point. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper introduces VPT as a new variational reformulation of the steady-state equation and presents two independent numerical strategies (LU decomposition and Krylov recycling) to avoid pseudo-inverses. No derivation step reduces by construction to a fitted parameter, self-defined quantity, or load-bearing self-citation; the claimed extension of the convergence radius is presented as a property of the variational ansatz itself rather than a renaming or statistical forcing of inputs. The abstract and extracted claims contain no equations or premises that equate the output to the input by definition.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated beyond standard perturbation-theory assumptions.

axioms (1)
  • domain assumption Steady-state equation of an open quantum system admits a perturbative expansion around reference parameters
    Implicit in the description of PT and VPT

pith-pipeline@v0.9.0 · 5751 in / 1073 out tokens · 35637 ms · 2026-05-22T21:42:30.026787+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.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Reduced Basis Method for Driven-Dissipative Quantum Systems

    cond-mat.str-el 2025-05 unverdicted novelty 6.0

    Generalization of reduced basis methods to driven-dissipative Markovian quantum systems with variance distillation for phase boundary detection.

Reference graph

Works this paper leans on

69 extracted references · 69 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    We choose ∆ = 0 as the initial point for PT and VPT

    1D example: varying detuning First, we vary ∆ and fix all other parameters. We choose ∆ = 0 as the initial point for PT and VPT. In Fig. 2(b) we compare the photon number in the steady state ˆa†ˆa = Tr[ˆa†ˆaˆρss] obtained from the exact solution to the problem [18], standard PT, and VPT. The shaded backgrounds indicates the region where∥Lˆρ P T ss ∥<tol f...

    work page

  2. [2]

    Compared to the previ- ous case, we consider a larger value ofκ/K, showcasing the efficiency of VPT in a more dissipative configura- tion

    2D example: simultaneously varying detuning and drive amplitude We now apply VPT to compute steady states in a re- gion where ∆ andFchange. Compared to the previ- ous case, we consider a larger value ofκ/K, showcasing the efficiency of VPT in a more dissipative configura- tion. We plot the average photon number in the steady state in Fig. 3(a). At large p...

    work page

  3. [3]

    ance (see the discussion in Appendix C)

    Fock space truncation: 70 photons. ance (see the discussion in Appendix C). To cover the parameter space, we first select a random point ( ¯∆, ¯F) and compute the steady state and its per- turbative corrections up to orderM, introducing the per- turbation parametersε 1 = ∆− ¯∆ andε 2 =F− ¯F. Next, we use (V)PT to construct approximate steady states in the...

    work page

  4. [4]

    Haroche and J

    S. Haroche and J. M. Raimond,Exploring the Quantum: Atoms, Cavities, and Photons(Oxford University Press, Oxford, 2006)

    work page 2006

  5. [5]

    Wiseman and G

    H. Wiseman and G. Milburn,Quantum Measurement and Control(Cambridge University Press, Cambridge, 2010)

    work page 2010

  6. [6]

    Breuer and F

    H. Breuer and F. Petruccione,The Theory of Open Quan- tum Systems(Oxford University Press, Oxford, 2007)

    work page 2007

  7. [7]

    Fazio, J

    R. Fazio, J. Keeling, L. Mazza, and M. Schir` o, Many- body open quantum systems (2024)

    work page 2024

  8. [8]

    Blais, A

    A. Blais, A. L. Grimsmo, S. M. Girvin, and A. Wallraff, Circuit quantum electrodynamics, Rev. Mod. Phys.93, 025005 (2021)

    work page 2021

  9. [9]

    Carusotto and C

    I. Carusotto and C. Ciuti, Quantum fluids of light, Rev. Mod. Phys.85, 299 (2013)

    work page 2013

  10. [10]

    Ritsch, P

    H. Ritsch, P. Domokos, F. Brennecke, and T. Esslinger, Cold atoms in cavity-generated dynamical optical poten- tials, Rev. Mod. Phys.85, 553 (2013)

    work page 2013

  11. [11]

    L. D. Carr, D. DeMille, R. V. Krems, and J. Ye, Cold and ultracold molecules: science, technology and appli- cations, New Journal of Physics11, 055049 (2009)

    work page 2009

  12. [12]

    Leibfried, R

    D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Quantum dynamics of single trapped ions, Rev. Mod. Phys.75, 281 (2003)

    work page 2003

  13. [13]

    Ferrari, L

    F. Ferrari, L. Gravina, D. Eeltink, P. Scarlino, V. Savona, and F. Minganti, Steady-state quantum chaos in open quantum systems (2023)

    work page 2023

  14. [14]

    Dahan, G

    D. Dahan, G. Arwas, and E. Grosfeld, Classical and quantum chaos in chirally-driven, dissipative bose-hubbard systems, npj Quantum Information8, 10.1038/s41534-022-00518-2 (2022)

    work page doi:10.1038/s41534-022-00518-2 2022

  15. [15]

    E. M. Kessler, G. Giedke, A. Imamoˇ glu, S. F. Yelin, M. D. Lukin, and J. I. Cirac, Dissipative phase transition in a central spin system, Phys. Rev. A86, 012116 (2012)

    work page 2012

  16. [16]

    H. J. Carmichael, Breakdown of photon blockade: A dis- sipative quantum phase transition in zero dimensions, Phys. Rev. X5, 031028 (2015)

    work page 2015

  17. [17]

    Minganti, A

    F. Minganti, A. Biella, N. Bartolo, and C. Ciuti, Spec- tral theory of liouvillians for dissipative phase transitions, Phys. Rev. A98, 042118 (2018)

    work page 2018

  18. [18]

    M. J. Hartmann, F. G. S. L. Brand˜ ao, and M. B. Plenio, Strongly interacting polaritons in coupled arrays of cav- ities,2007 European Conference on Lasers and Electro- Optics and the International Quantum Electronics Con- ference, Nature Physics2, 849–855 (2006)

    work page 2007

  19. [19]

    Roberts and A

    D. Roberts and A. A. Clerk, Driven-dissipative quantum Kerr resonators: New exact solutions, photon blockade and quantum bistability, Phys. Rev. X10, 021022 (2020)

    work page 2020

  20. [20]

    Minganti, N

    F. Minganti, N. Bartolo, J. Lolli, W. Casteels, and C. Ciuti, Exact results for schr¨ odinger cats in driven- dissipative systems and their feedback control, Scientific Reports6, 26987 (2016)

    work page 2016

  21. [21]

    P. D. Drummond and D. F. Walls, Quantum theory of op- tical bistability. i. nonlinear polarisability model, J. Phys. A: Math. Theor.13, 725 (1980)

    work page 1980

  22. [22]

    Weimer, A

    H. Weimer, A. Kshetrimayum, and R. Or´ us, Simulation methods for open quantum many-body systems, Rev. Mod. Phys.93, 015008 (2021)

    work page 2021

  23. [23]

    M. B. Plenio and P. L. Knight, The quantum-jump ap- proach to dissipative dynamics in quantum optics, Re- views of Modern Physics70, 101 (1998)

    work page 1998

  24. [24]

    A. J. Daley, Quantum trajectories and open many-body quantum systems, Advances in Physics63, 77 (2014)

    work page 2014

  25. [25]

    Finazzi, A

    S. Finazzi, A. Le Boit´ e, F. Storme, A. Baksic, and C. Ciuti, Corner-space renormalization method for driven-dissipative two-dimensional correlated systems, Physical review letters115, 080604 (2015). 14

    work page 2015

  26. [27]

    Kshetrimayum, H

    A. Kshetrimayum, H. Weimer, and R. Or´ us, A sim- ple tensor network algorithm for two-dimensional steady states, Nature communications8, 1291 (2017)

    work page 2017

  27. [28]

    A. H. Werner, D. Jaschke, P. Silvi, M. Kliesch, T. Calarco, J. Eisert, and S. Montangero, Positive tensor network approach for simulating open quantum many-body systems, Physical review letters116, 237201 (2016)

    work page 2016

  28. [29]

    Vicentini, A

    F. Vicentini, A. Biella, N. Regnault, and C. Ciuti, Vari- ational neural-network ansatz for steady states in open quantum systems, Phys. Rev. Lett.122, 250503 (2019)

    work page 2019

  29. [30]

    Yoshioka and R

    N. Yoshioka and R. Hamazaki, Constructing neural sta- tionary states for open quantum many-body systems, Physical Review B99, 214306 (2019)

    work page 2019

  30. [31]

    Deuar and P

    P. Deuar and P. D. Drummond, Correlations in a bec collision: first-principles quantum dynamics with 150 000 atoms, Physical review letters98, 120402 (2007)

    work page 2007

  31. [32]

    Polkovnikov, Phase space representation of quantum dynamics, Annals of Physics325, 1790 (2010)

    A. Polkovnikov, Phase space representation of quantum dynamics, Annals of Physics325, 1790 (2010)

    work page 2010

  32. [33]

    J. Jin, A. Biella, O. Viyuela, L. Mazza, J. Keeling, R. Fazio, and D. Rossini, Cluster mean-field approach to the steady-state phase diagram of dissipative spin sys- tems, Physical Review X6, 031011 (2016)

    work page 2016

  33. [34]

    Robichon and A

    G. Robichon and A. Tilloy, Bootstrapping the stationary state of bosonic open quantum systems, arXiv preprint arXiv:2410.07384 (2024)

    work page arXiv 2024

  34. [35]

    Weimer, Variational principle for steady states of dis- sipative quantum many-body systems, Phys

    H. Weimer, Variational principle for steady states of dis- sipative quantum many-body systems, Phys. Rev. Lett. 114, 040402 (2015)

    work page 2015

  35. [36]

    A. C. Li, F. Petruccione, and J. Koch, Perturbative ap- proach to markovian open quantum systems, Scientific reports4, 4887 (2014)

    work page 2014

  36. [37]

    Benatti, A

    F. Benatti, A. Nagy, and H. Narnhofer, Asymptotic en- tanglement and lindblad dynamics: a perturbative ap- proach, Journal of Physics A: Mathematical and Theo- retical44, 155303 (2011)

    work page 2011

  37. [38]

    A. C. Li, F. Petruccione, and J. Koch, Resummation for nonequilibrium perturbation theory and application to open quantum lattices, Physical Review X6, 021037 (2016)

    work page 2016

  38. [39]

    Christiansen, V

    H. Christiansen, V. V. Baran, and J. Paaske, Reduced ba- sis method for driven-dissipative quantum systems, Phys. Rev. Lett.135, 236503 (2025)

    work page 2025

  39. [40]

    Nagib and T

    O. Nagib and T. G. Walker, Exact steady state of per- turbed open quantum systems, Phys. Rev. Res.7, 033076 (2025)

    work page 2025

  40. [41]

    A. K. Noor and H. E. Lowder, Approximate techniques of strctural reanalysis, Computers & Structures4, 801 (1974)

    work page 1974

  41. [42]

    J. N. Silverman and J. C. van Leuven, Perturbational- variational approach to the calculation of variational wave functions. i. theory, Physical Review162, 1175 (1967)

    work page 1967

  42. [43]

    Garrigue and B

    L. Garrigue and B. Stamm, On reduced basis methods for eigenvalue problems, and on its coupling with pertur- bation theory (2025), arXiv:2408.11924 [math-ph]

    work page arXiv 2025

  43. [44]

    Gebhart, R

    V. Gebhart, R. Santagati, A. A. Gentile, E. M. Gauger, D. Craig, N. Ares, L. Banchi, F. Marquardt, L. Pezze, and C. Bonato, Learning quantum systems, Nature Re- views Physics5, 141 (2023)

    work page 2023

  44. [45]

    Beaulieu, F

    G. Beaulieu, F. Minganti, S. Frasca, V. Savona, S. Fe- licetti, R. D. Candia, and P. Scarlino, Observation of first- and second-order dissipative phase transitions in a two-photon driven kerr resonator (2023)

    work page 2023

  45. [46]

    Berdou, A

    C. Berdou, A. Murani, U. Reglade, W. C. Smith, M. Vil- liers, J. Palomo, M. Rosticher, A. Denis, P. Morfin, M. Delbecq,et al., One hundred second bit-flip time in a two-photon dissipative oscillator, PRX Quantum4, 020350 (2023)

    work page 2023

  46. [47]

    L. N. Trefethen and D. Bau,Classical Numerical Analysis (SIAM, 2022)

    work page 2022

  47. [48]

    We useLUsince its complexity scales asO(2N 3/3) for a system of sizeNwhereasQRscales asO(4N 3/3)

    Another common choice is theQRdecomposition. We useLUsince its complexity scales asO(2N 3/3) for a system of sizeNwhereasQRscales asO(4N 3/3)

    work page

  48. [49]

    Bartolo, F

    N. Bartolo, F. Minganti, W. Casteels, and C. Ciuti, Exact steady state of a Kerr resonator with one- and two-photon driving and dissipation: Controllable wigner- function multimodality and dissipative phase transitions, Phys. Rev. A94, 033841 (2016)

    work page 2016

  49. [50]

    Casteels, F

    W. Casteels, F. Storme, A. Le Boit´ e, and C. Ciuti, Power laws in the dynamic hysteresis of quantum nonlinear pho- tonic resonators, Phys. Rev. A93, 033824 (2016)

    work page 2016

  50. [51]

    S. G. Krantz and H. R. Parks,The implicit function the- orem: history, theory, and applications(Springer Science & Business Media, 2002)

    work page 2002

  51. [52]

    Leghtas, S

    Z. Leghtas, S. Touzard, I. M. Pop, A. Kou, B. Vlas- takis, A. Petrenko, K. M. Sliwa, A. Narla, S. Shankar, M. J. Hatridge, M. Reagor, L. Frunzio, R. J. Schoelkopf, M. Mirrahimi, and M. H. Devoret, Confining the state of light to a quantum manifold by engineered two-photon loss, Science347, 853 (2015)

    work page 2015

  52. [53]

    Lescanne, M

    R. Lescanne, M. Villiers, T. Peronnin, A. Sarlette, M. Delbecq, B. Huard, T. Kontos, M. Mirrahimi, and Z. Leghtas, Exponential suppression of bit-flips in a qubit encoded in an oscillator, Nature Physics16, 509 (2020)

    work page 2020

  53. [54]

    Virtanen, R

    P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Peter- son, W. Weckesser, J. Bright,et al., Scipy 1.0: fundamen- tal algorithms for scientific computing in python, Nature methods17, 261 (2020)

    work page 2020

  54. [55]

    D. C. Liu and J. Nocedal, On the limited memory bfgs method for large scale optimization, Mathematical pro- gramming45, 503 (1989)

    work page 1989

  55. [56]

    Saad and M

    Y. Saad and M. H. Schultz, Gmres: A generalized min- imal residual algorithm for solving nonsymmetric linear systems, SIAM Journal on scientific and statistical com- puting7, 856 (1986)

    work page 1986

  56. [57]

    H. A. Van der Vorst, Bi-cgstab: A fast and smoothly con- verging variant of bi-cg for the solution of nonsymmetric linear systems, SIAM Journal on scientific and Statistical Computing13, 631 (1992)

    work page 1992

  57. [58]

    M. R. Hestenes, E. Stiefel,et al.,Methods of conjugate gradients for solving linear systems, Vol. 49 (NBS Wash- ington, DC, 1952) p. 409

    work page 1952

  58. [59]

    Saad,Iterative methods for sparse linear systems (SIAM, 2003)

    Y. Saad,Iterative methods for sparse linear systems (SIAM, 2003)

    work page 2003

  59. [60]

    Steady-state solution methods for open quantum optical systems

    P. Nation, Steady-state solution methods for open quan- tum optical systems, arXiv preprint arXiv:1504.06768 10.48550/arxiv.1504.06768 (2015)

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1504.06768 2015

  60. [61]

    M. L. Parks, E. De Sturler, G. Mackey, D. D. Johnson, and S. Maiti, Recycling krylov subspaces for sequences of linear systems, SIAM Journal on Scientific Computing 28, 1651 (2006). 15

    work page 2006

  61. [62]

    K. M. Soodhalter, D. B. Szyld, and F. Xue, Krylov sub- space recycling for sequences of shifted linear systems, Applied Numerical Mathematics81, 105 (2014)

    work page 2014

  62. [63]

    K. M. Soodhalter, E. de Sturler, and M. E. Kilmer, A survey of subspace recycling iterative methods, GAMM- Mitteilungen43, e202000016 (2020)

    work page 2020

  63. [64]

    T. E. Lee, S. Gopalakrishnan, and M. D. Lukin, Uncon- ventional magnetism via optical pumping of interacting spin systems, Phys. Rev. Lett.110, 257204 (2013)

    work page 2013

  64. [65]

    J. Jin, A. Biella, O. Viyuela, L. Mazza, J. Keeling, R. Fazio, and D. Rossini, Cluster mean-field approach to the steady-state phase diagram of dissipative spin sys- tems, Phys. Rev. X6, 031011 (2016)

    work page 2016

  65. [66]

    R. Rota, F. Storme, N. Bartolo, R. Fazio, and C. Ciuti, Critical behavior of dissipative two-dimensional spin lat- tices, Phys. Rev. B95, 134431 (2017)

    work page 2017

  66. [67]

    R. Rota, F. Minganti, A. Biella, and C. Ciuti, Dynamical properties of dissipativeXY ZHeisenberg lattices, New J. Phys.20, 045003 (2018)

    work page 2018

  67. [68]

    Huybrechts, F

    D. Huybrechts, F. Minganti, F. Nori, M. Wouters, and N. Shammah, Validity of mean-field theory in a dissi- pative critical system: Liouvillian gap,PT-symmetric antigap, and permutational symmetry in thexyzmodel, Phys. Rev. B101, 214302 (2020)

    work page 2020

  68. [69]

    Biella, J

    A. Biella, J. Jin, O. Viyuela, C. Ciuti, R. Fazio, and D. Rossini, Linked cluster expansions for open quantum systems on a lattice, Phys. Rev. B97, 035103 (2018)

    work page 2018

  69. [70]

    Finazzi, A

    S. Finazzi, A. Le Boit´ e, F. Storme, A. Baksic, and C. Ciuti, Corner-space renormalization method for driven-dissipative two-dimensional correlated systems, Phys. Rev. Lett.115, 080604 (2015)

    work page 2015