Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

Phase-Space Framework for Noisy Intermediate-Scale Quantum Optical Neural Networks

T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Pith's one-line read The paper claims that bosonic quantum reservoir networks, simulated with the positive-P method, perform best at five to seven nodes and degrade at larger sizes because input photons are diluted across modes.

desk verdict A solid extension of positive-P to cascaded bosonic reservoirs with a plausible scaling result, but the missing trajectory counts and many-mode convergence checks leave a genuine artifact risk that referees should pin down. read the letter →

arxiv 2507.07684 v1 pith:I2HHG2UM submitted 2025-07-10 quant-ph cond-mat.dis-nncond-mat.quant-gas

classification quant-phcond-mat.dis-nncond-mat.quant-gas
keywords positive-PmethodquantumreservoircomputingbosonicneuralnetworksstateclassificationfeaturepredictionKerrnonlinearitynoisyintermediate-scaledriven-dissipativeBose-Hubbardmodel
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper establishes that the positive-P phase-space method can simulate quantum optical reservoir networks with dozens of bosonic modes, a regime where the conventional density-matrix description would require an impractical amount of memory. Using this framework, the authors show that reservoir performance on quantum state classification and squeezing-parameter prediction does not improve monotonically with the number of modes. For weak Kerr nonlinearity, accuracy peaks at five to seven reservoir nodes, then falls at fifteen and fifty nodes. The paper attributes the decline to the dilution of effective nonlinearity: the injected photons are spread across more modes, so the average occupation per node drops and the Kerr interaction acts more weakly. If this is right, it provides a concrete design rule for photonic quantum neuromorphic hardware: optimize the reservoir size against the input photon budget and nonlinearity strength rather than simply scaling up.

What carries the argument

The engine of the paper is the positive-P phase-space representation applied to a driven-dissipative Bose-Hubbard lattice with on-site Kerr nonlinearity $(U/2)\hat{b}^{\dagger 2}\hat{b}^2$ and nearest-neighbour hopping. In this representation the full density matrix is replaced by a positive probability distribution over doubled phase-space variables and the Lindblad master equation by stochastic differential equations whose trajectories can be evolved independently; observables such as mode occupations are recovered as averages over trajectories. The input quantum state is injected unidirectionally through a cascade-coupling term, and input states (coherent, squeezed vacuum, optical Schrödinger cat) are prepared by sampling their canonical positive-P distributions. The physical mechanism carrying the argument is that the effective nonlinearity at a node is set by the product of the Kerr strength $U$ and the local photon occupation, so spreading a fixed input photon budget over more modes lowers the effective nonlinearity and makes the reservoir behave more linearly. This is what lets the authors reach reservoirs of fifty modes, a scale whose density-matrix representation would be impossible to store.

What would settle it

Repeat the $N=15$ and $N=50$ classification runs with trajectory counts $S$ spanning at least an order of magnitude ($10^3$, $10^4$, $10^5$) and check that test accuracy and mean-squared error do not drift; the paper does not report $S$, so this check settles whether the large-reservoir decline is physical or a sampling artifact.

Watch

Extended reading notes

Core claim

The central discovery is that a bosonic quantum reservoir's usefulness is not monotonic in its size. In quantum state classification (Schrödinger-cat, squeezed, and coherent states) and in predicting the complex squeezing parameter of a squeezed vacuum, reservoirs with five to seven nodes give the best test accuracy and lowest mean-squared error at Kerr nonlinearity $U/\gamma = 0.02$, while $N=15$ and $N=50$ perform markedly worse. The mechanism the authors identify is dilution of effective nonlinearity: the input quantum state's photons are distributed across all reservoir modes, so the average occupation per mode falls as $N$ grows, and the on-site Kerr interaction, which is what makes the reservoir response nonlinear enough to encode quantum correlations, becomes effectively weaker. Higher Kerr nonlinearity compensates for this dilution and shifts the optimum to larger reservoirs. In the authors' formulation, reservoir performance is governed by the balance between reservoir size and input field density, not by size alone.

Load-bearing premise

The central claim rests on the assumption that the stochastic averages used for the reservoir occupations are converged and unbiased; if rare diverging trajectories or finite-sample noise biased the large-reservoir data, the performance decline could be a numerical artifact rather than a physical dilution of nonlinearity.

Editorial extensions

If this is right

  • Reservoir design should not blindly scale up: at weak Kerr nonlinearity the useful range sits near five to seven modes, and the optimum shifts upward only when $U/\gamma$ is increased.
  • The positive-P framework can serve as a benchmarking tool for intermediate-scale bosonic quantum reservoirs, including $N=50$ configurations whose density matrix would require an impractical amount of memory, before experimental implementation.
  • Even at zero Kerr nonlinearity, mode-occupation dynamics carry enough information to separate coherent states from squeezed and cat states, so linear reservoirs solve simple tasks while nonlinearity is needed for harder ones.
  • The simulation cost depends strongly on the input state: experimentally relevant cat states with mean occupation between 1 and 3 can be sampled accurately, while larger non-Gaussian states require $10^5$ to $10^6$ trajectories.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: if the dilution mechanism is the whole story, then for a fixed input photon budget the optimal reservoir size should grow roughly in proportion to the Kerr nonlinearity or the input intensity; this scaling law is not derived in the paper but is a direct testable consequence.
  • Editorial inference: the same dilution argument predicts that tasks requiring less nonlinearity, such as separating coherent states only, would tolerate larger reservoirs; the $U=0$ classification curve reported here is consistent with that prediction.
  • Editorial inference: a photonic-lattice experiment with $N=5$, $7$, and $15$ coupled modes, fixed injected photon number, and fixed Kerr strength could verify the non-monotonic accuracy curve, since the paper's numbers give a concrete prediction to aim for.
  • Editorial inference: extending the framework to time-multiplexed readout of a small reservoir is an open direction; the paper deliberately uses simultaneous readout of all modes, so whether the dilution penalty appears under time multiplexing is not yet known.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper develops a positive-P phase-space framework for simulating driven-dissipative bosonic quantum reservoir neural networks, and applies it to two tasks: classification of coherent, squeezed, and Schrödinger-cat states, and prediction of the complex squeezing parameter of squeezed states. The central numerical claim is that reservoir performance is non-monotonic in the number of modes, with an optimum at N=5-7 for low Kerr nonlinearity U/γ=0.02 and degradation at N=15 and N=50, attributed to dilution of input photons across modes and a consequent reduction of effective nonlinearity. The manuscript includes a derivation of the cascade-coupling terms in positive-P form, a semi-implicit numerical integration scheme, a stability analysis for a single mode, and an open data repository.

Significance. If the non-monotonic scaling result is correct, it is a useful and non-obvious design principle for quantum optical reservoir computers, and the positive-P framework would open a route to simulating intermediate-scale bosonic reservoirs that are inaccessible to exact density-matrix methods. The paper provides a self-contained derivation of the stochastic equations and makes its data publicly available. However, the central numerical observation is only as strong as the convergence and unbiasedness of the positive-P stochastic averages, and the manuscript does not currently establish these for the actual reservoir sizes studied.

major comments (3)
  1. [§III and Eq. (12)] The number of stochastic trajectories S used in the QSC and QFP reservoir simulations is never stated, despite Eq. (12) defining observables as averages over S realizations. For a fixed S, the stochastic sampling error does not decrease with N, while the injected photon number per reservoir node decreases roughly as 1/N; undersampling would therefore degrade large-N points selectively and could produce exactly the reported non-monotonic accuracy curve. The authors must report S, and provide a convergence study for the actual reservoir sizes (at least N=5, 7, 15, 50 at U/γ=0.02) showing that accuracy and MSE are stable when S is increased, with error bars computed over independent trajectory batches.
  2. [Appendix C2, Fig. 7] The stability analysis is performed for a single driven-dissipative mode without the cascade source term, reservoir hopping, or the many-mode lattice, and the sentence 'consistent with having no diverging trajectories (also in the many-mode runs)' is asserted without quantitative evidence. Since rare divergent or boundary-violating trajectories can bias the stochastic averages in Eq. (12), the authors should report, for each reservoir size and nonlinearity used in the main text, the fraction of trajectories that diverge or exceed a numerical threshold, and describe how such trajectories are treated in the reported averages.
  3. [§III, generally] No comparison is made against an exact or independent numerical method for the observables that feed the machine-learning tasks. A benchmark for small reservoirs (e.g., N=2 or 3) against a master-equation solution of Eq. (2) or Eq. (10) for the time-integrated occupations ⟨n_i⟩, and ideally for the resulting classifier accuracy, would establish that the positive-P averages are unbiased in this parameter regime. Without such a check, the central non-monotonic performance curve is not fully separated from possible method bias at large N.
minor comments (5)
  1. [Abstract and Conclusions] The phrase 'does not improve monotonously' should be 'does not improve monotonically'.
  2. [Eq. (10) and Appendix C] The normalization condition for η is given as η = Σ_j (W_in_j)^2 in the main text, but in Appendix C near Eq. (C6) it appears as Σ_i (W_in_i)^*; the latter is a typo and should be the squared modulus.
  3. [Appendix B, Eqs. (B3) and (B5)] The addition of +iU/2 as an 'autocorrelation correction' is introduced without explanation; a sentence clarifying its origin and its relation to the Ito representation in Eq. (B1) would substantially improve reproducibility.
  4. [Appendix C2, Fig. 7] The stability diagram would be easier to interpret if the color scale and the definition of 'fraction of convergent trajectories' (e.g., the threshold criterion for divergence) were specified in the caption or text.
  5. [§IV C] The statement that a reservoir density matrix for the sizes studied here would require 'over 10^37 TB' appears inconsistent with the stated d^{2N} scaling for the reported values of d and N; the numerical estimate should be rechecked (for d=10 and N=50 the density matrix dimension is 10^100, which is many orders of magnitude larger than 10^37 TB).

Circularity Check

0 steps flagged · score 2.0 of 10

No equation-level circularity: the non-monotonic reservoir-size result is generated from the model rather than fitted to itself; self-citations are tool citations, not load-bearing reductions.

full rationale

The central simulation pipeline is self-contained: input states are sampled from their positive-P distributions, evolved through the stochastic equations (6) and (11), reduced to observables via Eq. (12), and then read out by a linear classifier. The reported optimum at N=5-7 emerges from these simulations and is not used to construct them. The paper does not fit any parameter to the classification accuracy or MSE and then rename that fit as a prediction. The 'effective nonlinearity' explanation is an interpretation of the simulated trend, not an input assumption that forces the trend. The authors cite their own prior phase-space work (e.g., refs. [4], [40]) as methodological tools, including the single-mode stability simplification in Appendix C2, but these citations are used to justify the numerical framework, not to import the central claim. The unspecified trajectory number S and the asserted absence of many-mode diverging trajectories are robustness or correctness concerns rather than circularity: they question whether the simulation is converged, not whether the derivation reduces to its inputs. No equation or fitted quantity is equivalent by construction to the claimed result, so the paper is not significantly circular.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central scaling result is conditional on hand-chosen nonlinearity values, an input-occupation window around 1-3 photons, a specific disorder model, and the unverified assumption that finite-trajectory positive-P sampling is converged in the many-mode regime. No new physical entities are postulated.

free parameters (4)
  • Kerr nonlinearity U/γ = 0.0, 0.02, 0.05, 0.1
    Hand-chosen values scanned in Figs. 3-4; the optimal reservoir size and the onset of performance degradation depend on U, so the central scaling claim is conditional on this choice.
  • Input state parameter ranges = |β| in (1.03,1.38), r in (0.9,1.1), giving ⟨n⟩ ~ 1-3
    Chosen to create overlapping average occupations across the three classes (Sec. IIIA); the abstract lists 'average occupation of the input mode' as a driver of the scaling, so conclusions apply in this window.
  • Reservoir disorder distributions = Δ_j uniform in (0,0.1γ); J_ij uniform in (-1,1), normalized by spectral radius
    Random connectivity and detunings are drawn from chosen distributions; the performance curves are averages over these realizations, so the conclusions are specific to this disorder model.
  • Integration time step and window = τ = 0.05, T = 25
    Set by trial and error (App. B); the feature vector n is a time integral over this window, so results depend on these numerical choices.
assumptions (4)
  • standard math The canonical positive-P distribution exists and can be sampled for coherent, squeezed, cat, and thermal states.
    Used in Sec. IID and App. C1; based on Drummond-Gardiner representation [28] and canonical construction [41,42]; exact in the S→∞ limit.
  • domain assumption The open system dynamics are Markovian and described by the Lindblad master equation with local dissipation γ and no thermal noise.
    Eqs. (2) and (10); standard for photonic reservoirs, excludes non-Markovian and finite-temperature effects.
  • domain assumption The source-reservoir coupling is unidirectional (cascade) with no back-action.
    Sec. IID, following [27]; required for the input-injection protocol used in all numerical experiments.
  • ad hoc to paper Single-mode stability analysis (Fig. 7) transfers to the many-mode lattice.
    App. C2 asserts no diverging trajectories in many-mode runs based on single-mode checks; this transfer is assumed, not proven, and it is load-bearing for the reliability of large-N results.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Phase-Space Framework for Noisy Intermediate-Scale Quantum Optical Neural Networks." pith.science (2026). https://pith.science/paper/I2HHG2UM

@misc{pith2026250707684,
  author       = {Pith},
  title        = {Pith review of: Phase-Space Framework for Noisy Intermediate-Scale Quantum Optical Neural Networks},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/I2HHG2UM}},
  note         = {Machine review of arXiv:2507.07684}
}
read the original abstract

Quantum optical neural networks (QONNs) enable information processing beyond classical limits by exploiting the advantages of classical and quantum optics. However, simulation of large-scale bosonic lattices remains a significant challenge due to the exponential growth of the Hilbert space required to describe a quantum network accurately. Consequently, previous theoretical studies have been limited to small-scale systems, leaving the behaviour of multimode QONNs largely unexplored. This work presents an efficient computational framework based on the phase-space positive-P method for simulating bosonic neuromorphic systems. This approach provides a view to previously inaccessible regimes, allowing the validation of large-scale bosonic networks in various quantum machine learning tasks such as quantum state classification and quantum state feature prediction. Our results show that the performance of a large quantum reservoir does not improve monotonously with the number of bosonic modes, instead following a complex dependence driven by the interplay of nonlinearity, reservoir size, and the average occupation of the input mode. These findings are essential for designing and optimising optical bosonic reservoirs for future quantum neuromorphic computing devices.

Figures

Figures reproduced from arXiv: 2507.07684 by the authors.

Figure 1
Figure 1. Schematic illustration of a quantum reservoir applied for quantum state recognition. (a), Structure of the QRNN, depicting input, reservoir, and output layers. The reservoir contains nine nonlinear bosonic nodes that transform input quantum states into real-valued observables given by average occupations ⟨nˆi⟩, where i indexes the nodes. (b), samples from the training dataset consisting of various quantum states, re… view at source ↗
Figure 2
Figure 2. Comparison of Wigner W(q, p) (top row) and positive-P (α, α˜ ∗ ) distributions (bottom row) for selected quantum states. To visualise the positive-P distributions, the complex parameters α and α˜ are separated into their real and imaginary parts, forming two-dimensional projections. Blue and purple points represent positions of α˜ and α, respectively (in reduced phase space). The number of stochastic trajectories us… view at source ↗
Figure 3
Figure 3. Quantum state recognition with quantum reservoirs. (a) Distribution of average occupation ⟨nˆ⟩ of source quantum states used during the training phase, illus￾trating the sampling over the quantum state parameter space. (b) Average (over three runs) test accuracy as a function of reservoir size for four different nonlinearity values, depicted on the plot. Reservoirs examined in this study were composed of N ∈ {2, 3, … view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: d illustrates the MSE of the predicted squeez￾ing parameter as a function of reservoir size for various Kerr nonlinearities. For all reservoir sizes, the error de￾creases with increasing nonlinearity, consistent with the expected behaviour of reservoir networks. Predic…
Figure 5
Figure 5. Figure 5: Complexity of QRC with fermionic and bosonic reservoirs. This figure illustrates the computa￾tional complexity of QRC, showing how the Hilbert space dimension of the reservoir scales with the number of nodes for fermionic (black dashed line) and bosonic (violet solid l…
Figure 6
Figure 6. Figure 6: Training process for a reservoir of size N = 7 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]
Figure 7
Figure 7. Figure 7: The stability of the positive-P method for simu￾lating the dynamics of a single reservoir mode for different parameter values. Panel a shows the stability as a function of driving field amplitude F and Kerr nonlinearity U. The increase in nonlinearity significantly dec…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Optimizing Quantum Photonic Integrated Circuits using Differentiable Tensor Networks

    quant-ph 2025-09 unverdicted novelty 6.0 of 10

    Gradient-based optimization of quantum photonic circuits is achieved via differentiable tensor networks that model nonlinear unitary gates and stochastic losses at low photon numbers.

Reference graph

Works this paper leans on

61 extracted references · 45 canonical work pages · cited by 1 Pith paper

  1. [1]

    (C9) We have used this to generate cat states

    Details on quantum state sampling The canonical distribution corresponding to any density matrix ˆρ is given constructively by [28] Pcan(α, ˜α∗) ≡ 1 4π2 e− 1 4 |α− ˜α|2 D (α+ ˜α) 2 ˆρ (α+ ˜α) 2 E . (C9) We have used this to generate cat states. For other quan- tum states portrayed in Figure (2), the canonical distri- bution functions, together with a conv...

  2. [2]

    A significant contribution of the noise term and low dissi- pation can contribute to large fluctuations propagating

    Stability analysis Section B of the Appendix shows that the amplitude of Gaussian noise depends on the Kerr nonlinearity and the amplitude of the corresponding field (e.g.α, ˜α). A significant contribution of the noise term and low dissi- pation can contribute to large fluctuations propagating. Therefore, corresponding trajectories may diverge. It is impo...

  3. [3]

    Labay-Mora, J

    A. Labay-Mora, J. García-Beni, G. L. Giorgi, M. C. So- riano, and R. Zambrini, Neural networks with quantum states of light, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sci- ences 382, 20230346 (2024)

  4. [4]

    S.-Y. Ma, T. Wang, J. Laydevant, L. G. Wright, and P. L. McMahon, Quantum-limited stochastic optical neu- ral networks operating at a few quanta per activation, Nature Communications16, 359 (2025)

  5. [5]

    Matuszewski, A

    M. Matuszewski, A. Prystupiuk, and A. Opala, Role of all-optical neural networks, Phys. Rev. Appl.21, 014028 (2024)

  6. [6]

    Ghosh, K

    S. Ghosh, K. Nakajima, T. Krisnanda, K. Fujii, and T. C. H. Liew, Quantum neuromorphic computing with reservoir computing networks, Advanced Quantum Tech- nologies 4, 2100053 (2021)

  7. [7]

    Ghosh, A

    S. Ghosh, A. Opala, M. Matuszewski, T. Paterek, and T. C. H. Liew, Quantum reservoir processing, npj Quan- tum Information5, 35 (2019)

  8. [8]

    Marković and J

    D. Marković and J. Grollier, Quantum neuromorphic computing, Applied Physics Letters117, 150501 (2020)

Show all 61 references
  1. [9]

    Krisnanda, S

    T. Krisnanda, S. Ghosh, T. Paterek, W. Laskowski, and T. C. Liew, Phase measurement beyond the standard quantum limit using a quantum neuromorphic platform, Phys. Rev. Appl.18, 034011 (2022)

  2. [10]

    Mujal, R

    P. Mujal, R. Martínez-Peña, J. Nokkala, J. García-Beni, G. L. Giorgi, M. C. Soriano, and R. Zambrini, Oppor- tunities in quantum reservoir computing and extreme learning machines, Advanced Quantum Technologies4, 2100027 (2021)

  3. [11]

    Yao, C.-J

    Y.-C. Yao, C.-J. Lee, Y.-J. Chen, J.-Z. Feng, H. Oh, C.-S. Lue, J.-K. Sheu, and Y.-J. Lee, All-inorganic per- ovskite quantum-dot optical neuromorphic synapses for near-sensor colored image recognition, Advanced Science 12, 2409933 (2025)

  4. [12]

    Krisnanda, S

    T. Krisnanda, S. Ghosh, T. Paterek, and T. C. Liew, Cre- ating and concentrating quantum resource states in noisy environments using a quantum neural network, Neural Networks136, 141 (2021)

  5. [13]

    D. Ko, S. Świerczewski, A. Opala, M. Matuszewski, and A. Rahmani, Estimation of the second-order coherence function using quantum reservoir and ensemble methods (2025), arXiv:2504.18205 [quant-ph]

  6. [14]

    Ghosh, T

    S. Ghosh, T. Paterek, and T. C. H. Liew, Quantum neu- romorphic platform for quantum state preparation, Phys. Rev. Lett.123, 260404 (2019)

  7. [15]

    Krisnanda, T

    T. Krisnanda, T. Paterek, M. Paternostro, and T. C. H. Liew, Quantum neuromorphic approach to efficient sens- ing of gravity-induced entanglement, Phys. Rev. D107, 086014 (2023)

  8. [16]

    H. Xu, T. Krisnanda, R. Bao, and T. C. H. Liew, The roles of Kerr nonlinearity in a bosonic quantum neural network, New Journal of Physics25, 023028 (2023)

  9. [17]

    T.Krisnanda, H.Xu, S.Ghosh,andT.C.H.Liew,Tomo- graphic completeness and robustness of quantum reser- voir networks, Phys. Rev. A107, 042402 (2023)

  10. [18]

    Kirchmair, B

    G. Kirchmair, B. Vlastakis, Z. Leghtas, S. E. Nigg, H. Paik, E. Ginossar, M. Mirrahimi, L. Frunzio, S. M. Girvin, and R. J. Schoelkopf, Observation of quantum state collapse and revival due to the single-photon Kerr effect, Nature495, 205 (2013)

  11. [19]

    Krisnanda, K

    T. Krisnanda, K. Dini, H. Xu, W. Verstraelen, and T. C. Liew, Wisdom of crowds in quantum machine learning, Phys. Rev. Appl.19, 034010 (2023)

  12. [20]

    Ghosh, T

    S. Ghosh, T. Krisnanda, T. Paterek, and T. C. H. Liew, Realising and compressing quantum circuits with quan- tum reservoir computing, Communications Physics 4, 105 (2021)

  13. [21]

    Preskill, Quantum Computing in the NISQ era and beyond, Quantum2, 79 (2018)

    J. Preskill, Quantum Computing in the NISQ era and beyond, Quantum2, 79 (2018). 15

  14. [22]

    Grimm, N

    A. Grimm, N. E. Frattini, S. Puri, S. O. Mundhada, S. Touzard, M. Mirrahimi, S. M. Girvin, S. Shankar, and M. H. Devoret, Stabilization and operation of a Kerr-cat qubit, Nature584, 205 (2020)

  15. [23]

    Dawid, J

    A. Dawid, J. Arnold, B. Requena, A. Gresch, M. Płodzień, K. Donatella, K. A. Nicoli, P. Stornati, R. Koch, M. Büttner, R. Okuła, G. Muñoz-Gil, R. A. Vargas-Hernández, A. Cervera-Lierta, J. Carrasquilla, V. Dunjko, M. Gabrié, P. Huembeli, E. van Nieuwenburg, F. Vicentini, L. Wa...

  16. [24]

    Kounalakis, C

    M. Kounalakis, C. Dickel, A. Bruno, N. K. Langford, and G. A. Steele, Tuneable hopping and nonlinear cross-Kerr interactions in a high-coherence superconducting circuit, npj Quantum Information4, 38 (2018)

  17. [25]

    Cheng, X.-H

    B. Cheng, X.-H. Deng, X. Gu, Y. He, G. Hu, P. Huang, J. Li, B.-C. Lin, D. Lu, Y. Lu, C. Qiu, H. Wang, T. Xin, S.Yu, M.-H.Yung, J.Zeng, S.Zhang, Y.Zhong, X.Peng, F. Nori, and D. Yu, Noisy intermediate-scale quantum computers, Frontiers of Physics18, 21308 (2023)

  18. [26]

    Wigner, On the quantum correction for thermody- namic equilibrium, Phys

    E. Wigner, On the quantum correction for thermody- namic equilibrium, Phys. Rev.40, 749 (1932)

  19. [27]

    J. C. L. Carreño and F. P. Laussy, Excitation with quan- tum light. i. exciting a harmonic oscillator, Phys. Rev. A 94, 063825 (2016)

  20. [28]

    Cramer, A

    M. Cramer, A. Bernard, N. Fabbri, L. Fallani, C. Fort, S. Rosi, F. Caruso, M. Inguscio, and M. B. Plenio, Spa- tial entanglement of bosons in optical lattices, Nature Communications 4, 2161 (2013)

  21. [29]

    Tanaka, T

    G. Tanaka, T. Yamane, J. B. Héroux, R. Nakane, N. Kanazawa, S. Takeda, H. Numata, D. Nakano, and A. Hirose, Recent advances in physical reservoir comput- ing: A review, Neural Networks115, 100 (2019)

  22. [30]

    Deuar and P

    P. Deuar and P. D. Drummond, First-principles quan- tum dynamics in interacting Bose gases: I. the positive p representation, Journal of Physics A: Mathematical and General 39, 1163 (2006)

  23. [31]

    P. D. Drummond and C. W. Gardiner, Generalised P- representations in quantum optics, Journal of Physics A: Mathematical and General13, 2353 (1980)

  24. [32]

    Deuar and P

    P. Deuar and P. D. Drummond, Gaugep representations for quantum-dynamical problems: Removal of boundary terms, Phys. Rev. A66, 033812 (2002)

  25. [33]

    P. D. Drummond and J. F. Corney, Quantum dynamics of evaporatively cooled Bose-Einstein condensates, Phys. Rev. A60, R2661 (1999)

  26. [34]

    C. W. Gardiner and P. Zoller, Quantum Noise: A handbook of Markovian and non-Markovian Quantum stochastic methods with applications to Quantum Optics (Springer, 2010)

  27. [35]

    S. J. Carter, P. D. Drummond, M. D. Reid, and R. M. Shelby, Squeezing of quantum solitons, Phys. Rev. Lett. 58, 1841 (1987)

  28. [36]

    Perrin, C

    A. Perrin, C. M. Savage, D. Boiron, V. Krachmalnicoff, C. I. Westbrook, and K. V. Kheruntsyan, Atomic four- wave mixing via condensate collisions, New Journal of Physics 10, 045021 (2008)

  29. [37]

    Carusotto, Nonlinear atomic Fabry-Perot interferom- eter: From the mean-field theory to the atom blockade effect, Phys

    I. Carusotto, Nonlinear atomic Fabry-Perot interferom- eter: From the mean-field theory to the atom blockade effect, Phys. Rev. A63, 023610 (2001)

  30. [38]

    Deuar and P

    P. Deuar and P. D. Drummond, Correlations in a BEC collision: First-principles quantum dynamics with 150000 atoms, Phys. Rev. Lett.98, 120402 (2007)

  31. [39]

    Ng and E

    R. Ng and E. S. Sørensen, Exact real-time dynamics of quantum spin systems using the positive-P representa- tion,JournalofPhysicsA:MathematicalandTheoretical 44, 065305 (2011)

  32. [40]

    are simple, with a zero diffusive term. Since stabil- ity was shown to depend almost exclusively on the rela- tionshipbetweenoccupation, decayandnonlinearity[40], one can expect similar stability regions to the dissipative Bose-Hubbard model, which turns out to be the case – s...

  33. [41]

    K. V. Kheruntsyan, J.-C. Jaskula, P. Deuar, M. Bon- neau, G. B. Partridge, J. Ruaudel, R. Lopes, D. Boiron, and C. I. Westbrook, Violation of the Cauchy-Schwarz in- equality with matter waves, Phys. Rev. Lett.108, 260401 (2012)

  34. [42]

    R. J. Lewis-Swan and K. V. Kheruntsyan, Proposal for demonstrating the Hong–Ou–Mandel effect with matter waves, Nature Commun.5, 3752 (2014)

  35. [43]

    Deuar, A

    P. Deuar, A. Ferrier, M. Matuszewski, G. Orso, and M. H. Szymańska, Fully quantum scalable description of driven-dissipative lattice models, PRX Quantum 2, 010319 (2021)

  36. [44]

    S. L. Braunstein, C. M. Caves, and G. J. Milburn, Inter- pretation for a positive p representation, Phys. Rev. A 43, 1153 (1991)

  37. [45]

    Olsen and A

    M. Olsen and A. Bradley, Numerical representation of quantum states in the positive-P and Wigner represen- tations, Optics Communications282, 3924 (2009)

  38. [46]

    P. D. Drummond and B. Opanchuk, Initial states for quantum field simulations in phase space, Phys. Rev. Re- search 2, 033304 (2020)

  39. [47]

    X. Pan, Z. Lu, W. Wang, Z. Hua, Y. Xu, W. Li, W. Cai, X. Li, H. Wang, Y.-P. Song, C.-L. Zou, D.-L. Deng, and L. Sun, Deep quantum neural networks on a super- conducting processor, Nature Communications14, 4006 (2023)

  40. [48]

    Gilchrist, C

    A. Gilchrist, C. W. Gardiner, and P. D. Drummond, Pos- itive P representation: Application and validity, Phys. Rev. A55, 3014 (1997)

  41. [49]

    S. A. Khan, F. Hu, G. Angelatos, and H. E. Türeci, Phys- ical reservoir computing using finitely-sampled quantum systems (2021), arXiv:2110.13849 [quant-ph]

  42. [50]

    Ghosh, A

    S. Ghosh, A. Opala, M. Matuszewski, T. Paterek, and T.C.H.Liew,Reconstructingquantumstateswithquan- tum reservoir networks, IEEE Transactions on Neural Networks and Learning Systems32, 3148 (2021)

  43. [51]

    Suzuki, Q

    Y. Suzuki, Q. Gao, K. C. Pradel, K. Yasuoka, and N. Ya- mamoto, Natural quantum reservoir computing for tem- poral information processing, Scientific Reports12, 1353 (2022)

  44. [52]

    M. J. Steel, M. K. Olsen, L. I. Plimak, P. D. Drummond, S. M. Tan, M. J. Collett, D. F. Walls, and R. Graham, Dynamical quantum noise in trapped Bose-Einstein con- densates, Phys. Rev. A58, 4824 (1998)

  45. [53]

    Sinatra, C

    A. Sinatra, C. Lobo, and Y. Castin, The truncated Wigner method for Bose-condensed gases: limits of va- lidity and applications1, Journal of Physics B: Atomic, Molecular and Optical Physics35, 3599 (2002)

  46. [54]

    G. Orso, J. Zakrzewski, and P. Deuar, Self-organized cav- ity bosons beyond the adiabatic elimination approxima- tion, Phys. Rev. Lett.134, 183405 (2025)

  47. [55]

    Verstraelen, R

    W. Verstraelen, R. Rota, V. Savona, and M. Wouters, Gaussian trajectory approach to dissipative phase transi- tions: The case of quadratically driven photonic lattices, Phys. Rev. Res.2, 022037 (2020)

  48. [56]

    A. J. D. and, Quantum trajectories and open many-body quantum systems, Advances in Physics 63, 77 (2014), https://doi.org/10.1080/00018732.2014.933502

  49. [57]

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

  50. [58]

    Pan and Z

    G. Pan and Z. Y. Meng, The sign problem in quan- tum monte carlo simulations, in Encyclopedia of Con- densed Matter Physics (Second Edition) , edited by T. Chakraborty (Academic Press, Oxford, 2024) second edition ed., pp. 879–893

  51. [59]

    Świerczewski, V

    S. Świerczewski, V. W., P. Deuar, B. Piętka, T. C. H. Liew, M. Matuszewski, and A. Opala, Phase-space framework for noisy intermediate-scale quantum optical neural networks, https://github.com/StanislawSw/ Phase-Space-Framework-for-Noisy-Intermediate-Scale-QONN (2025). 16

  52. [60]

    Drummond, Central partial difference propagation al- gorithms, Computer Physics Communications 29, 211 (1983)

    P. Drummond, Central partial difference propagation al- gorithms, Computer Physics Communications 29, 211 (1983)

  53. [61]

    Deuar, Multi-time correlations in the positive-P, Q, and doubled phase-space representations, Quantum 5, 455 (2021)

    P. Deuar, Multi-time correlations in the positive-P, Q, and doubled phase-space representations, Quantum 5, 455 (2021)

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.