pith. sign in

arxiv: 2605.11103 · v2 · pith:MVL6QFEOnew · submitted 2026-05-11 · 🪐 quant-ph · cond-mat.mes-hall· cond-mat.quant-gas

Permutation-symmetric quantum trajectories

Elliot W. Lloyd , Aleksandra A. Ziolkowska , Jonathan Keeling This is my paper

Pith reviewed 2026-05-20 21:57 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.mes-hallcond-mat.quant-gas
keywords stochastic unravelingpermutation symmetryquantum trajectoriesopen quantum systemscavity QEDemitter ensemblescollective dynamicscomputational scaling
0
0 comments X

The pith

A stochastic unraveling that respects weak permutation symmetry reduces the simulation cost for N identical emitters coupled to a cavity from O(N^5) to O(N^2) for two-level systems.

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

The paper shows a method for stochastic unraveling of open quantum systems with many identical emitters that preserves a weak permutation symmetry among them. This symmetry allows equivalent emitter configurations to be grouped, shrinking the space in which trajectories must be computed. For two-level emitters the cost drops from fifth power to second power in the number of emitters, and further refinements reach linear scaling. The averaged results still match the exact master equation dynamics. The same approach works for d-level emitters, with cost scaling as the power d(d-1)/2, and it already permits large-N runs when d equals three.

Core claim

We show how one may perform a stochastic unraveling which respects weak permutation symmetry for models of N emitters coupled to a common system (e.g. a cavity mode). For problems involving 2-level emitters, such an unravelling reduces the computational cost from O(N^5) to O(N^2), and with additional refinements, allows reduction to O(N). This significantly increases the range of system sizes for which one can model exact quantum dynamics of such systems. We further show how the method can also be applied to d-level systems, with computational effort scaling as O(N^{d(d-1)/2}), and we show it allows large-N simulations for d=3.

What carries the argument

The weak-permutation-symmetric stochastic unraveling, which groups symmetric emitter states so that trajectories evolve in a reduced effective space while reproducing the exact open-system statistics upon averaging.

If this is right

  • Exact quantum trajectory statistics are recovered for any model that admits the required symmetry.
  • Simulations of collective effects in emitter-cavity systems become feasible at emitter numbers previously out of reach.
  • The scaling improvement extends to three-level emitters, opening large-N studies for more complex atomic structures.
  • Refinements that reach linear scaling further widen the accessible system sizes without sacrificing exactness.

Where Pith is reading between the lines

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

  • The same symmetry grouping could be combined with trajectory reweighting or importance sampling to handle even larger or less symmetric cases.
  • Applications to superradiance, synchronization, or many-body quantum optics problems would benefit directly from the reduced cost.
  • The approach suggests analogous symmetry reductions might be developed for other identical-particle models outside cavity QED.

Load-bearing premise

The dynamics or initial state must permit a weak permutation symmetry that the unraveling can respect while still reproducing the exact open-system evolution.

What would settle it

Direct numerical comparison, for N equal to four or five two-level emitters, between the ensemble-averaged density matrix obtained from the symmetric trajectories and the solution of the full Lindblad master equation, checking agreement within sampling error.

Figures

Figures reproduced from arXiv: 2605.11103 by Aleksandra A. Ziolkowska, Elliot W. Lloyd, Jonathan Keeling.

Figure 2
Figure 2. Figure 2: FIG. 2. Superradiant emission from the Tavis-Cummings [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Inversionless lasing of a three-level system [ [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
read the original abstract

We show how one may perform a stochastic unraveling which respects weak permutation symmetry for models of $N$ emitters coupled to a common system (e.g. a cavity mode). For problems involving 2-level emitters, such an unravelling reduces the computational cost from $\mathcal{O}(N^5)$ to $\mathcal{O}(N^2)$, and with additional refinements, allows reduction to $\mathcal{O}(N)$. This significantly increases the range of system sizes for which one can model exact quantum dynamics of such systems. We further show how the method can also be applied to d-level systems, with computational effort scaling as $\mathcal{O}(N^{d(d-1)/2})$, and we show it allows large-N simulations for d=3.

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 / 1 minor

Summary. The manuscript presents a stochastic unraveling technique for open quantum systems of N emitters coupled to a common mode (e.g., cavity) that respects weak permutation symmetry. For two-level emitters it claims an exact reduction in cost from O(N^5) to O(N^2), with further refinements reaching O(N); the method is generalized to d-level emitters with claimed scaling O(N^{d(d-1)/2}), enabling large-N simulations at least for d=3.

Significance. If the scaling and exactness claims are substantiated, the approach would meaningfully enlarge the range of N for which exact quantum trajectory simulations of collective open-system dynamics are feasible, which is valuable for quantum-optics models exhibiting superradiance or cavity-mediated interactions.

major comments (1)
  1. [Abstract and d-level scaling discussion] Abstract and d-level generalization: the claimed scaling O(N^{d(d-1)/2}) yields O(N^3) for d=3, whereas the dimension of the symmetric subspace is binom(N+d-1,d-1) ~ O(N^{d-1}) = O(N^2). The manuscript must clarify whether the weak-symmetry unraveling employs a larger effective space, dense operations, or an alternative parameterization that produces this cost while still reproducing the exact master-equation evolution.
minor comments (1)
  1. The phrase 'with additional refinements' for the O(N) scaling in the two-level case would benefit from a short pointer to the relevant subsection or equation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for identifying the need to clarify the scaling claims for the d-level generalization. We address this point directly below and have revised the manuscript to improve the explanation of our parameterization and its relation to the symmetric subspace dimension.

read point-by-point responses

  1. Referee: [Abstract and d-level scaling discussion] Abstract and d-level generalization: the claimed scaling O(N^{d(d-1)/2}) yields O(N^3) for d=3, whereas the dimension of the symmetric subspace is binom(N+d-1,d-1) ~ O(N^{d-1}) = O(N^2). The manuscript must clarify whether the weak-symmetry unraveling employs a larger effective space, dense operations, or an alternative parameterization that produces this cost while still reproducing the exact master-equation evolution.

    Authors: We thank the referee for this observation, which highlights an important point of clarification. Our weak-permutation-symmetry unraveling does not operate in a larger effective space; it exactly preserves and evolves within the symmetric subspace. The claimed scaling O(N^{d(d-1)/2}) arises from the specific alternative parameterization we employ for the state and the associated operators in the stochastic unraveling. Rather than using the standard occupation-number basis (of dimension O(N^{d-1})), we represent the permutation-symmetric states via a structured multi-index tensor parameterization whose update rules and jump operators admit efficient structured multiplications. This yields operation costs scaling as O(N^{d(d-1)/2}) per time step while remaining exactly equivalent to the master-equation dynamics restricted to the symmetric subspace. For d=3 this produces O(N^3) cost, which is lower than a naive dense implementation on the O(N^2)-dimensional space. We have revised the abstract and added a dedicated explanatory paragraph in the d-level section to describe this parameterization and its relation to the symmetric subspace dimension. revision: yes

Circularity Check

0 steps flagged

No circularity: new algorithmic construction for symmetry-respecting unraveling

full rationale

The paper presents a stochastic unraveling method that respects weak permutation symmetry for N emitters coupled to a common system, claiming exact reproduction of open-system dynamics with reduced cost O(N^2) for two-level systems (or O(N) with refinements) and O(N^{d(d-1)/2}) for d-level systems. This is framed as an independent algorithmic development rather than a derivation from equations or parameters that reduce to the inputs by construction. No self-definitional steps, fitted inputs renamed as predictions, load-bearing self-citations, or ansatzes smuggled via prior work appear in the provided abstract and context. The derivation chain is self-contained as a novel technique under stated symmetry assumptions, with the skeptic's scaling concern relating to optimality rather than circular reduction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method extends standard quantum trajectory techniques without introducing new fitted parameters or postulated entities.

axioms (1)
  • standard math Standard Lindblad master equation and stochastic unraveling framework for open quantum systems.
    The symmetry-respecting variant is built on established quantum trajectory methods.

pith-pipeline@v0.9.0 · 5658 in / 1152 out tokens · 76555 ms · 2026-05-20T21:57:57.104115+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

42 extracted references · 42 canonical work pages

  1. [1]

    Haken, The semiclassical and quantum theory of the laser, inQuantum Opt., edited by Kay and Maitland (1970) p. 201

    work page 1970

  2. [2]

    Scully and M

    M. Scully and M. Zubairy,Quantum Optics, Quantum Optics (Cambridge University Press, 1997)

    work page 1997

  3. [3]

    Bonifacio and G

    R. Bonifacio and G. Preparata, Coherent spontaneous emission, Phys. Rev. A2, 336 (1970)

    work page 1970

  4. [4]

    Hepp and E

    K. Hepp and E. Lieb, Equilibrium statistical mechanics of matter interacting with the quantized radiation field, Phys. Rev. A8, 2517 (1973)

    work page 1973

  5. [5]

    Y. K. Wang and F. T. Hioe, Phase Transition in the Dicke Model of Superradiance, Phys. Rev. A7, 831 (1973)

    work page 1973

  6. [6]

    Abouelela, M

    A. Abouelela, M. Turaev, R. Kramer, M. Janning, M. Kajan, S. Ray, and J. Kroha, Stabilizing open pho- ton condensates by ghost-attractor dynamics, Phys. Rev. Lett.135, 053402 (2025)

    work page 2025

  7. [7]

    J. A. ´Cwik, S. Reja, P. B. Littlewood, and J. Keeling, Polariton condensation with saturable molecules dressed by vibrational modes, Eur. Lett.105, 47009 (2014)

    work page 2014

  8. [8]

    Galego, F

    J. Galego, F. J. Garcia-Vidal, and J. Feist, Cavity- induced modifications of molecular structure in the strong-coupling regime, Phys. Rev. X5, 041022 (2015)

    work page 2015

  9. [9]

    Herrera and F

    F. Herrera and F. C. Spano, Cavity-Controlled Chem- istry in Molecular Ensembles, Phys. Rev. Lett.116, 238301 (2016)

    work page 2016

  10. [10]

    Buˇ ca and T

    B. Buˇ ca and T. Prosen, A note on symmetry reductions of the Lindblad equation: transport in constrained open spin chains, New J. Phys.14, 73007 (2012)

    work page 2012

  11. [11]

    V. V. Albert and L. Jiang, Symmetries and conserved quantities in lindblad master equations, Phys. Rev. A 89, 022118 (2014)

    work page 2014

  12. [12]

    B. A. Chase and J. M. Geremia, Collective processes of an ensemble of spin-1/2 particles, Phys. Rev. A78, 052101 (2008)

    work page 2008

  13. [13]

    M. Xu, D. A. Tieri, and M. J. Holland, Simulating open quantum systems by applying su(4) to quantum master equations, Phys. Rev. A87, 062101 (2013)

    work page 2013

  14. [14]

    Richter, M

    M. Richter, M. Gegg, T. S. Theuerholz, and A. Knorr, Numerically exact solution of the many emitter–cavity laser problem: Application to the fully quantized spaser emission, Phys. Rev. B91, 035306 (2015)

    work page 2015

  15. [15]

    Kirton and J

    P. Kirton and J. Keeling, Superradiant and lasing states in driven-dissipative Dicke models, New J. Phys.20, 15009 (2018)

    work page 2018

  16. [16]

    Shammah, S

    N. Shammah, S. Ahmed, N. Lambert, S. De Liberato, and F. Nori, Open quantum systems with local and col- lective incoherent processes: Efficient numerical simula- tions using permutational invariance, Phys. Rev. A98, 063815 (2018)

    work page 2018

  17. [17]

    Freter, P

    L. Freter, P. Fowler-Wright, J. Cuerda, B. W. Lovett, J. Keeling, and P. T¨ orm¨ a, Theory of dynamical super- radiance in organic materials, Nanophotonics14, 5323 (2025)

    work page 2025

  18. [18]

    X. Cao, A. Tilloy, and A. D. Luca, Entanglement in a fermion chain under continuous monitoring, SciPost Phys.7, 24 (2019)

    work page 2019

  19. [19]

    Y. Li, X. Chen, and M. P. A. Fisher, Quantum Zeno effect and the many-body entanglement transition, Phys. Rev. B98, 205136 (2018)

    work page 2018

  20. [20]

    Skinner, J

    B. Skinner, J. Ruhman, and A. Nahum, Measurement- Induced Phase Transitions in the Dynamics of Entangle- ment, Phys. Rev. X9, 31009 (2019)

    work page 2019

  21. [21]

    S. Choi, Y. Bao, X.-L. Qi, and E. Altman, Quan- tum Error Correction in Scrambling Dynamics and Measurement-Induced Phase Transition, Phys. Rev. Lett.125, 30505 (2020)

    work page 2020

  22. [22]

    M. J. Gullans and D. A. Huse, Dynamical Purification Phase Transition Induced by Quantum Measurements, Phys. Rev. X10, 41020 (2020)

    work page 2020

  23. [23]

    Gisin and I

    N. Gisin and I. C. Percival, The quantum-state diffu- sion model applied to open systems, J. Phys. A25, 5677 (1992)

    work page 1992

  24. [24]

    Dalibard, Y

    J. Dalibard, Y. Castin, and K. Mølmer, Wave-function approach to dissipative processes in quantum optics, Phys. Rev. Lett.68, 580 (1992)

    work page 1992

  25. [25]

    R. Dum, P. Zoller, and H. Ritsch, Monte Carlo simulation of the atomic master equation for spontaneous emission, Phys. Rev. A45, 4879 (1992)

    work page 1992

  26. [26]

    H. J. Carmichael, Quantum trajectory theory for cas- caded open systems, Phys. Rev. Lett.70, 2273 (1993)

    work page 1993

  27. [27]

    M. B. Plenio and P. L. Knight, The quantum-jump ap- proach to dissipative dynamics in quantum optics, Rev. Mod. Phys.70, 101 (1998)

    work page 1998

  28. [28]

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

    work page 2014

  29. [29]

    Vovk and H

    T. Vovk and H. Pichler, Entanglement-optimal trajec- tories of many-body quantum markov processes, Phys. Rev. Lett.128, 243601 (2022)

    work page 2022

  30. [30]

    Barberena, Generalized holstein-primakoff mapping and 1/nexpansion of collective spin systems undergo- ing single particle dissipation (2025), arXiv:2508.05751 [quant-ph]

    D. Barberena, Generalized holstein-primakoff mapping and 1/nexpansion of collective spin systems undergo- ing single particle dissipation (2025), arXiv:2508.05751 [quant-ph]

    work page arXiv 2025

  31. [31]

    Supplementary Material for Permutation symmetric quantum trajectories

    work page

  32. [32]

    Kirton and J

    P. Kirton and J. Keeling, Suppressing and restoring the dicke superradiance transition by dephasing and decay, Phys. Rev. Lett.118, 123602 (2017)

    work page 2017

  33. [33]

    A. J. Daley, J. M. Taylor, S. Diehl, M. Baranov, and P. Zoller, Atomic three-body loss as a dynamical three- body interaction, Phys. Rev. Lett.102, 040402 (2009)

    work page 2009

  34. [34]

    (1) these are ˆx i,1 =σ x i /2, ˆXc,1 = (g/ √ N)(ˆa+ˆa†) and ˆxi,2 =σ y i /2, ˆXc,2 =i(g/ √ N)(ˆa−ˆa†)

    In the notation of Eq. (1) these are ˆx i,1 =σ x i /2, ˆXc,1 = (g/ √ N)(ˆa+ˆa†) and ˆxi,2 =σ y i /2, ˆXc,2 =i(g/ √ N)(ˆa−ˆa†)

    work page

  35. [35]

    Schack, T

    R. Schack, T. A. Brun, and I. C. Percival, Quantum state diffusion, localization and computation, Journal of Physics A: Mathematical and General28, 5401 (1995)

    work page 1995

  36. [36]

    M¨ uller and W

    K. M¨ uller and W. T. Strunz, Quantum trajectory method for highly excited environments in non-markovian open quantum dynamics, Phys. Rev. A112, 033719 (2025)

    work page 2025

  37. [37]

    Bastin and J

    T. Bastin and J. Martin, Permutationally invariant pro- cesses in open multiqudit systems, J. Phys. A:Math. Gen. 58, 275301 (2025)

    work page 2025

  38. [38]

    Ceccherini-Silberstein, F

    T. Ceccherini-Silberstein, F. Scarabotti, and F. Tolli, Representation Theory of the Symmetric Groups: The Okounkov-Vershik Approach, Character Formulas, and Partition Algebras, Cambridge Studies in Advanced Mathematics (Cambridge University Press, 2010)

    work page 2010

  39. [39]

    I. M. Gelfand and M. L. Tsetlin, Finite-dimensional rep- resentations of the group of unimodular matrices, Dokl. Akad. Nauk SSSR71, 825 (1950)

    work page 1950

  40. [40]

    Werren, E

    N. Werren, E. M. Gauger, and P. Kirton, A quantum model of lasing without inversion, New J. Phys.24, 93027 6 (2022)

    work page 2022

  41. [41]

    Suess, A

    D. Suess, A. Eisfeld, and W. T. Strunz, Hierarchy of stochastic pure states for open quantum system dynam- ics, Phys. Rev. Lett.113, 150403 (2014)

    work page 2014

  42. [42]

    PERMUTATION-SYMMETRIC QUANTUM TRAJECTORIES

    K. M¨ uller, K. Luoma, and C. Sch¨ afer, A hierarchical ap- proach to quantum many-body systems in structured en- vironments (2025), arXiv:2405.05093 [quant-ph]. 1 SUPPLEMENTAL MATERIAL FOR: “PERMUTATION-SYMMETRIC QUANTUM TRAJECTORIES” EXPRESSIONS FOR COEFFICIENTSf ˆX In the main text, we noted that the results of Refs. [12, 30] can be written in terms of...

    work page arXiv 2025