Pith. sign in

REVIEW 4 major objections 5 minor 2 cited by

Variational Multi-Gaussian Phase-Space Dynamics via Automatic Differentiation

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

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper claims that a variational ansatz of complex-centered Gaussians for the Wigner function, evolved by the Dirac-Frenkel principle and evaluated by automatic differentiation, matches exact diagonalization and scales to…

desk verdict A genuinely useful variational method for open bosonic systems whose critical-dynamics section overclaims: the 'extracted' dynamical exponents are imported from the literature and the gap extraction lacks documented convergence and fitting details. read the letter →

arxiv 2507.14076 v2 pith:X2G7BSI6 submitted 2025-07-18 quant-ph cond-mat.other

Jacopo Tosca , Francesco Carnazza , Luca Giacomelli , Cristiano Ciuti This is my paper · ORCID
classification quant-phcond-mat.other
keywords openquantumsystemsWignerfunctionvariationalmethodDirac-FrenkelprincipleautomaticdifferentiationBose-HubbardmodelcriticalslowingdownLiouvillianspectralgap
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 introduces a variational method for open quantum bosonic systems that approximates the Wigner function as a sum of complex-centered Gaussian components. It derives equations of motion from the Dirac-Frenkel principle and evaluates all phase-space integrals analytically via automatic differentiation, so accuracy improves systematically as the number of Gaussians is increased. The central claim is that this method matches exact diagonalization on a single-mode benchmark, including negative Wigner fringes, and scales to a 12x12 driven-dissipative Bose-Hubbard lattice. From the relaxation dynamics the authors extract the Liouvillian spectral gap and report a finite-size collapse with the critical exponents of the 2D quantum Ising universality class, implying critical slowing down in the thermodynamic limit. If correct, this provides a scalable route to strongly correlated open bosonic dynamics beyond existing tensor-network and phase-space methods.

What carries the argument

The central object is the Variational Multi-Gaussian ansatz $W_\theta(\xi)=\sum_i \mathrm{Re}[G(\xi;\theta_i)]$ for the Wigner function, where each $G$ is a normalized Gaussian with complex center $\mu = \alpha + i\beta$. Complex centers produce oscillatory modulations of the Gaussian envelope and capture negative interference fringes; the even-parity ansatz adds mirrored Gaussians. The dynamics is governed by the Dirac-Frenkel equations $T\,d\theta/dt = V$, where $T$ is the Wigner-space quantum geometric tensor (overlap of parameter derivatives) and $V$ the Liouvillian gradient; both reduce to generalized Gaussian moments. The key technical move is to evaluate these moments as derivatives of a closed-form generating function $Z[J,\tilde J]$, using Taylor-mode automatic differentiation to handle derivatives up to order six, which keeps the method linear in the number of modes and Gaussians.

What would settle it

Compute the Liouvillian gap for a small 2D lattice, such as 4x4, both with the VMG method and with exact diagonalization in a truncated Fock basis; if the VMG gap does not converge to the exact value as the number of Gaussians increases, the finite-time relaxation estimate is not trustworthy. Additionally, refit the 12x12 gap data with twice as many Gaussian components; if the extracted $z$ shifts outside the reported uncertainty, the critical-exponent claim is not converged.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that a Variational Multi-Gaussian ansatz for the Wigner function, evolved with the Dirac-Frenkel principle and evaluated with Taylor-mode automatic differentiation, gives controlled and scalable open quantum bosonic dynamics. The ansatz represents negative Wigner regions through complex-centered Gaussians, and the paper shows exponential reduction of observable error with the number of Gaussian components in the single-mode Kerr parametric oscillator. Applied to a two-dimensional Bose-Hubbard lattice with two-boson driving and losses, the method reproduces steady-state parity from corner-space renormalization and extends to lattice sizes up to 12x12, i.e. 144 modes. The paper's strongest claim is the finite-size scaling collapse of the Liouvillian gap using exponents $\beta = 0.32641871$, $\nu = 0.62997097$, and $z = 2.0235$, identifying the transition as 2D quantum Ising and demonstrating critical slowing down.

Load-bearing premise

The method's central assumption is that the variational equations of motion stay accurate over very long times on the large lattice, so the decay rate read off from finite trajectories equals the true asymptotic Liouvillian gap; convergence is verified explicitly only on the single-mode test case, not on the 12x12 lattice used for the critical exponents.

Editorial extensions

If this is right

  • The method provides a systematically convergent variational route for open bosonic systems deep in the quantum regime, including transient Wigner negativities.
  • For the studied driven-dissipative Bose-Hubbard model, the Liouvillian gap vanishes in the thermodynamic limit with 2D quantum Ising exponents, so the phase transition shows critical slowing down.
  • The approach extends steady-state benchmarks to lattices of 144 modes, beyond previously accessible corner-space or tensor-network limits.
  • In the single-mode benchmark, increasing the number of Gaussians by a factor of four reduces the relative observable error by roughly four orders of magnitude, indicating controlled convergence with $N_G$.
  • Because the Liouvillian acts as a polynomial differential operator in phase space, the same automatic-differentiation machinery can be applied to other analytical ansatze and phase-space representations.

Reading between the lines

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

  • Beyond the paper's benchmarks, the same machinery should extend to spin or fermionic phase-space representations, since the Liouvillian differential structure is generic; the authors themselves point toward such extensions.
  • The reported $z = 2.0235$ is taken from a five-loop epsilon expansion of the 2D Ising model, not derived by the variational method; a direct extraction of $z$ by fitting the VMG gap data without fixing the Ising value would be a stronger independent test.
  • One could test the method on quench dynamics that create cat states or on regimes with multiple steady states; the paper does not explore these, and they would probe the ansatz's expressivity beyond the studied parameter window.
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

4 major / 5 minor

Summary. The paper introduces a variational method for open quantum bosonic systems based on a Wigner phase-space representation. The ansatz, Eq. (5), is a sum of complex Gaussian components, and the time evolution of its parameters is obtained from the Dirac-Frenkel principle, Eq. (8), with the Wigner quantum geometric tensor and Liouvillian gradient evaluated analytically through Gaussian moments and automatic differentiation. The method is benchmarked on a single-mode driven-dissipative Kerr parametric oscillator against exact diagonalization, including Wigner-function snapshots and exponential convergence with the number of Gaussians (Section V). It is then applied to two-dimensional Bose-Hubbard lattices with two-boson driving and losses (Section VI): the steady-state parity is compared with corner-space renormalization for a 6x6 lattice and computed for lattices up to 12x12, and the Liouvillian gap is extracted from relaxation dynamics (Section VII). The central physics claim is that the finite-size scaling of the gap reveals critical slowing down with dynamical exponents of the 2D quantum Ising universality class, using beta=0.32641871, nu=0.62997097, and z=2.0235.

Significance. The methodological core is valuable and largely well executed. The single-mode benchmarks in Section V, especially the exponential control of the error with the number of Gaussians shown in Fig. 3(c), are convincing and demonstrate that the ansatz can capture strong Wigner negativity. The steady-state comparison against the corner-space renormalization method in Fig. 5(a) provides an independent check for a genuinely many-body lattice. The analytical Gaussian-moment framework combined with Taylor-mode automatic differentiation is a genuine technical contribution, and the reported scalability to a 12x12 lattice with 6920 variational parameters is impressive. If the lattice dynamics and critical-scaling claims are substantiated, this would be a significant advance for open bosonic many-body simulation. The main weakness is that the central physics claim rests on a gap-extraction and scaling-collapse analysis that currently lacks convergence checks, error bars, and a clear fitting protocol; the collapse is a consistency check under externally fixed exponents rather than an independent extraction.

major comments (4)
  1. [Section VII, Fig. 7] The extraction of the Liouvillian gap lambda is not described: no fitting function, fit window, or uncertainty estimates are given, and Fig. 7(a) shows that at G=Gc the parity has not reached its steady state within the simulation time for the larger lattices. The reported lambda is therefore inferred from a transient, and the identification of this rate with the asymptotic Liouvillian gap is unsupported. Please specify the fitting procedure, show that a single-exponential decay describes the data over a stable window, and demonstrate convergence of lambda with respect to the final integration time.
  2. [Section VI and Appendix D] No convergence in the number of Gaussians NG is shown for the lattice dynamics. The exponential convergence of Fig. 3(c) is demonstrated only for the single-mode Kerr oscillator, whereas the 12x12 results use NG=16 with the restricted parametrization of Appendix D, whose covariance matrices contain only single-mode squeezing and no inter-mode blocks. Because this ansatz manifold is far more constrained in the lattice case, the authors should show, at least for a 6x6 or 8x8 lattice, that the steady-state parity and the extracted relaxation rate are stable under increasing NG, for example NG=4, 8, 16, 32.
  3. [Section VII, Fig. 7(c)] The finite-size collapse is a consistency check under assumed exponents rather than an independent extraction: beta, nu, and z are imported from the 2D quantum Ising literature, and Gc about 1.0 gamma is estimated from the same finite-size data that is then rescaled. The manuscript should quantify the sensitivity of the collapse to the choice of Gc and to the assumed exponents, and should report error bars on lambda and Gc; without this, the claim that the method extracts dynamical exponents of the 2D quantum Ising universality class is overstated.
  4. [Section VI and Appendix D] The main text states that the reader is referred to Appendix D for a detailed description of the initial conditions used for the calculated dynamics, but Appendix D contains only the variational parametrization and does not specify the initial conditions. This information is necessary to reproduce the relaxation dynamics and to assess whether the extracted long-time rate depends on the initial preparation.
minor comments (5)
  1. [Section I and Appendix B] There are several typos that should be corrected: 'systems systems' in the introduction, 'powerfull' and 'architechtures' in Section I, and 'Lindbald' in Appendix B.
  2. [Appendix B1, Eq. (B8)] The onsite sum in Eq. (B8) runs to NG rather than to M, which is inconsistent with the main-text Hamiltonian in Eq. (17); this appears to be a typo.
  3. [Section III, Eq. (8)] The Dirac-Frenkel equations require solving a linear system with the quantum geometric tensor T, but the manuscript does not discuss the conditioning or possible rank-deficiency of T, nor any regularization used in the numerical solution; a brief statement on how Eq. (8) is inverted in practice would aid reproducibility.
  4. [Section V, Fig. 3 and Section VII] No data or code availability statement is provided; given the complexity of the implementation, releasing the code or at least the data underlying Figs. 5 and 7 would strengthen the paper.
  5. [Section VII] The phrase 'first calculation of the dynamical critical exponent' should be tempered: since the collapse uses literature exponents and an estimated Gc, the analysis is better described as a consistency check, not a standalone determination of z.

Circularity Check

1 steps flagged · score 4.0 of 10

The VMG method itself is benchmarked externally, but the claimed 'extraction' of 2D quantum Ising dynamical exponents is an import of literature values presented as a numerical prediction.

  1. fitted input called prediction [Section VII (Fig. 7c) and Section VIII (Conclusions)]
    "Remarkably, in Fig. 7(c), we demonstrate that the driven-dissipative Bose-Hubbard model belongs to the Ising universality class by rescaling with the corresponding dynamical critical exponent z = 2.0235 [71]. ... we revealed a vanishing asymptotic decay rate and extracted critical exponents belonging to the 2D quantum Ising universality class."

    The claimed output, the dynamical critical exponent and the Ising universality class, is not computed from the VMG data. The values used for the collapse, beta = 0.32641871, nu = 0.62997097 from refs. [61-63] and z = 2.0235 from ref. [71], are inputs taken from the literature. A data collapse with these fixed exponents is a consistency test under an assumed universality class, not an independent extraction. The conclusion that the paper 'extracted critical exponents belonging to the 2D quantum Ising universality class' reverses the logical direction: the exponents appear as inputs in Eq. (7c) rescaling and are then reported as outputs.

full rationale

The variational machinery is not circular. The Dirac-Frenkel equations of motion (Eq. 8, T dtheta/dt = V) are derived from the ansatz in Eq. (5), with T and V evaluated from the explicit generating-function integrals in Appendix A. The method is benchmarked against external exact diagonalization for the single-mode Kerr parametric oscillator (Figs. 2 and 3), where convergence in the number of Gaussians is demonstrated, and against steady-state Corner Space Renormalization data for the Bose-Hubbard lattice (Fig. 5a). These benchmarks give the VMG approach independent support and do not reduce to the paper's own claims. The self-citations to the CSR method and to the Liouvillian spectral-gap framework are references to established techniques and are not load-bearing in a circular way. The main circularity-adjacent problem is the critical-exponent claim: beta, nu, and z are taken from 2D quantum Ising literature and used to rescale the data, after which the paper reports that it 'extracted' these exponents. This is a consistency check presented as a calculation, and it is the reason the score is not 0. Other concerns about the many-body gap result, such as the absence of NG-convergence checks at the 12x12 critical point and the unstated gap-extraction procedure, are correctness risks rather than circularity and are not counted in the score. Overall, the method is self-contained and externally benchmarked, so the circularity is limited to the exponent-attribution step.

Assumptions & free parameters 1 free parameters · 6 assumptions · 0 invented entities

The method introduces no new physical entities. The free parameters are limited to the critical drive Gc, estimated from the data, and the imported Ising exponents. The main burdens are the expressiveness of the variational ansatz and the identification of finite-time decay with the asymptotic gap.

free parameters (1)
  • Gc = approximately 1.0 gamma
    Critical drive amplitude estimated from the crossing of finite-size parity curves (Fig. 5 inset) and then used to rescale both the steady-state parity and the Liouvillian gap. No uncertainty is reported.
assumptions (6)
  • domain assumption The Lindblad master equation and the Wigner phase-space representation are exact descriptions of the open bosonic system.
    Standard framework used from Eq. (1) onward; the paper assumes Markovian weak coupling and uses the exact Wigner representation, not the truncated Wigner approximation.
  • domain assumption The Dirac-Frenkel variational principle projected onto the finite multi-Gaussian manifold gives a faithful approximation of the exact dynamics for the observables of interest.
    This is the core method assumption in Eq. (8). It is validated for a single-mode benchmark but not for the 12x12 lattice that supports the critical scaling claim.
  • ad hoc to paper The ansatz in Eq. (D1), with even-paired complex Gaussian centers and block-diagonal per-mode squeezed covariances, is sufficiently expressive to capture the relevant critical dynamics.
    The restriction to block-diagonal covariances with per-mode squeezing parameters is a modeling choice that limits multi-mode correlations; the paper argues these are captured by multiple Gaussian components, but this is not rigorously shown.
  • domain assumption The critical exponents beta, nu, and z of the 2D quantum Ising model apply to this driven-dissipative Bose-Hubbard model.
    The scaling collapse in Fig. 5(c) and Fig. 7(c) uses beta = 0.32641871, nu = 0.62997097, and z = 2.0235 from external literature; the conclusion that the model belongs to the Ising universality class depends on this imported hypothesis.
  • domain assumption The asymptotic decay rate obtained from finite-time relaxation equals the Liouvillian spectral gap.
    The paper identifies the late-time decay of observables with the Liouvillian gap, but simulations only reach 20 to 120 gamma^-1, and the extraction method is not fully specified.
  • standard math Gaussian integral identities and Fourier transform manipulations used in Appendix A are correct.
    Standard analytic Gaussian integration is used to derive the generating function Z[J, Jtilde]; no formal proof is given, but the steps are standard.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Variational Multi-Gaussian Phase-Space Dynamics via Automatic Differentiation." pith.science (2026). https://pith.science/paper/X2G7BSI6

@misc{pith2026250714076,
  author       = {Pith},
  title        = {Pith review of: Variational Multi-Gaussian Phase-Space Dynamics via Automatic Differentiation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/X2G7BSI6}},
  note         = {Machine review of arXiv:2507.14076}
}
read the original abstract

We introduce a variational method for simulating the dynamics of interacting open quantum bosonic systems deep in the quantum regime. The method is based on a multi-dimensional Wigner phase-space representation and employs a Variational Multi-Gaussian (VMG) ansatz, whose accuracy is systematically controlled by the number of Gaussian components. The variational equations of motion are derived from the Dirac-Frenkel principle and evaluated efficiently by combining the analytical structure of Gaussian functions with automatic differentiation. As a key first physical application, we study a driven-dissipative two-dimensional Bose-Hubbard lattice with two-boson coherent driving and two-body losses. Using our dynamical approach, we compute the finite-size scaling of the Liouvillian spectral gap, extracted from the relaxation dynamics, which vanishes in the thermodynamic limit. Our results reveal critical slowing down with dynamical exponents of the 2D quantum Ising universality class, demonstrating the power of our method to capture complex quantum dynamics in large open systems.

Figures

Figures reproduced from arXiv: 2507.14076 by the authors.

Figure 1
Figure 1. FIG. 1. Scheme of the variational method demonstrated in this work to calculate accurately the dynamics of a quantum [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (a) Wigner function at time [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Dynamics of observables for the same open quantum system as in Fig. [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Scheme of the considered 12 [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. (a) Comparison between the steady-state expectation value of the global parity operator obtained using both the [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Time evolution of the average boson number [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. (a) Illustrative example of the temporal dynamics of the parity operator expectation value, shown on a log scale, at [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

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

  1. Phase diagram of lasing under correlated pump from GPU-accelerated Truncated Wigner dynamics

    quant-ph 2026-08 conditional novelty 7.0 of 10

    Power-law correlated incoherent pumping of a superradiant laser preserves coherent, ultra-narrow emission down to pump-range exponent alpha≈1 while reducing the required pump rate by up to a factor N^{1-alpha}.

  2. Quantum and Classical Potts Criticality in Driven-Dissipative Bosonic Lattices

    quant-ph 2026-07 conditional novelty 7.0 of 10

    Nonequilibrium steady states of cubically driven dissipative Bose-Hubbard lattices realize classical 2D and quantum 1D three-state Potts criticality according to the presence of three-photon losses.

Reference graph

Works this paper leans on

78 extracted references · 58 canonical work pages · cited by 2 Pith papers

  1. [1]

    (17)) which we also report here for the sake of simplicity: ˆH = NG ∑ j=1 [− ∆ˆa† j ˆaj+ U 2 ˆa†2 j ˆa2 j + G 2 ˆa†2 j + G∗ 2 ˆa2 j] − ∑ ⟨j,j′⟩ J z (ˆa† j ˆaj′+ ˆa† j′ ˆaj)

    Bose-Hubbard Lindbladian in phase-space Let us now consider a M -modes (sites) driven dissipa- tive Bose-Hubbard Hamiltonian (Eq. (17)) which we also report here for the sake of simplicity: ˆH = NG ∑ j=1 [− ∆ˆa† j ˆaj+ U 2 ˆa†2 j ˆa2 j + G 2 ˆa†2 j + G∗ 2 ˆa2 j] − ∑ ⟨j,j′⟩ J z (ˆa† j ˆaj′+ ˆa† j′ ˆaj) . (B8) 13 For clarity, we will denote by LGW the expli...

    work page

  2. [2]

    The parity operator ˆΠ acts on the position and mo- mentum basis as [73]: ˆΠ∣q⟩= ∣− q⟩; ˆΠ∣p⟩= ∣− p⟩

    Parity W eyl symbol In this subsection, we calculate the Weyl symbol cor- responding to the boson number parity operator ˆΠ = exp(iπ ∑i ˆa† i ˆai). The parity operator ˆΠ acts on the position and mo- mentum basis as [73]: ˆΠ∣q⟩= ∣− q⟩; ˆΠ∣p⟩= ∣− p⟩. (B15) The Weyl symbol of the parity operator ˆΠ of a 1-mode bosonic system can be computed directly from th...

    work page

  3. [3]

    Carusotto, A

    I. Carusotto, A. A. Houck, A. J. Koll´ ar, P. Roushan, D. I. Schuster, and J. Simon, Photonic materials in circuit quantum electrodynamics, Nat. Phys. 16, 268 (2020)

    2020

  4. [4]

    Carusotto and C

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

    2013

  5. [5]

    Ozawa, H

    T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Simon, O. Zil- berberg, and I. Carusotto, Topological photonics, Rev. Mod. Phys. 91, 015006 (2019)

    2019

  6. [6]

    To address this technical challenge, we have employed Taylor-mode automatic differentiation [56], a mathemati- cally efficient—though still underutilized—technique [57] based on expanding functions into their Taylor series. This approach achieves at least factorial improvements in computational efficiency with respect to standard auto- matic differentiati...

    work page 2025

  7. [7]

    Mirrahimi, Z

    M. Mirrahimi, Z. Leghtas, V. V. Albert, S. Touzard, R. J. Schoelkopf, L. Jiang, and M. H. Devoret, Dynamically protected cat-qubits: a new paradigm for universal quan- tum computation, New J. Phys. 16, 045014 (2014)

    work page 2014

  8. [8]

    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, Science 347, 853 (2015)

    2015

Show all 78 references
  1. [9]

    D. J. Weigand and B. M. Terhal, Generating grid states from Schr¨ odinger-cat states without postselection, Phys. Rev. A 97, 022341 (2018)

    2018

  2. [10]

    Guillaud and M

    J. Guillaud and M. Mirrahimi, Repetition Cat Qubits for Fault-Tolerant Quantum Computation, Phys. Rev. X 9, 041053 (2019)

    2019

  3. [11]

    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, Nature 584, 205 (2020)

    2020

  4. [12]

    M. V. Larsen, X. Guo, C. R. Breum, J. S. Neergaard-Nielsen, and U. L. Andersen, Deterministic multi-mode gates on a scalable photonic quantum com- puting platform, Nat. Phys. 17, 1018 (2021)

    2021

  5. [13]

    Joshi, K

    A. Joshi, K. Noh, and Y. Y. Gao, Quantum information processing with bosonic qubits in circuit QED, Quantum Science and Technology 6, 033001 (2021)

    2021

  6. [14]

    Gautier, A

    R. Gautier, A. Sarlette, and M. Mirrahimi, Combined Dissipative and Hamiltonian Confinement of Cat Qubits, PRX Quantum 3, 020339 (2022)

    2022

  7. [15]

    Chamberland, K

    C. Chamberland, K. Noh, P. Arrangoiz-Arriola, E. T. Campbell, C. T. Hann, J. Iverson, H. Putterman, T. C. Bohdanowicz, S. T. Flammia, A. Keller, G. Refael, J. Preskill, L. Jiang, A. H. Safavi-Naeini, O. Painter, and F. G. Brand˜ ao, Building a Fault-Tolerant Quantum Computer U...

    2022

  8. [16]

    X. Pan, J. Schwinger, N.-N. Huang, P. Song, W. Chua, F. Hanamura, A. Joshi, F. Valadares, R. Filip, and Y. Y. Gao, Protecting the Quantum Interference of Cat States by Phase-Space Compression, Phys. Rev. X 13, 021004 (2023)

    2023

  9. [17]

    A. Hajr, B. Qing, K. Wang, G. Koolstra, Z. Pedramrazi, Z. Kang, L. Chen, L. B. Nguyen, C. J¨ unger, N. Goss, I. Huang, B. Bhandari, N. E. Frattini, S. Puri, J. Dres- sel, A. N. Jordan, D. I. Santiago, and I. Siddiqi, High- Coherence Kerr-Cat Qubit in 2D Architecture, Phys. Rev...

    2024

  10. [18]

    Marquet, A

    A. Marquet, A. Essig, J. Cohen, N. Cottet, A. Murani, E. Albertinale, S. Dupouy, A. Bienfait, T. Peronnin, S. Jezouin, R. Lescanne, and B. Huard, Autoparametric Resonance Extending the Bit-Flip Time of a Cat Qubit up to 0.3 s, Phys. Rev. X 14, 021019 (2024)

    2024

  11. [19]

    A. Z. Ding, B. L. Brock, A. Eickbusch, A. Koottan- davida, N. E. Frattini, R. G. Corti˜ nas, V. R. Joshi, S. J. de Graaf, B. J. Chapman, S. Ganjam, L. Frunzio, R. J. Schoelkopf, and M. H. Devoret, Quantum control of an oscillator with a Kerr-cat qubit, Nat. Commun. 16, 5279 (2025)

    2025

  12. [20]

    Hutin, P

    H. Hutin, P. Bilous, C. Ye, S. Abdollahi, L. Cros, T. Dvir, T. Shah, Y. Cohen, A. Bienfait, F. Marquardt, and B. Huard, Preparing Schr¨ odinger Cat States in a Mi- crowave Cavity Using a Neural Network, PRX Quantum 6, 010321 (2025)

    2025

  13. [21]

    M. Yuan, A. Seif, A. Lingenfelter, D. I. Schuster, A. A. Clerk, and L. Jiang, Universal control in bosonic systems with weak Kerr nonlinearities, Phys. Rev. A 111, 032606 (2025)

    2025

  14. [22]

    Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Ann

    U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Ann. Phys. 326, 96 (2011)

    2011

  15. [23]

    Or´ us, Tensor networks for complex quantum systems, Nat

    R. Or´ us, Tensor networks for complex quantum systems, Nat. Rev. Phys. 1, 538 (2019)

    2019

  16. [24]

    Kilda, A

    D. Kilda, A. Biella, M. Schir` o, R. Fazio, and J. Keel- ing, On the stability of the infinite Projected Entangled Pair Operator ansatz for driven-dissipative 2D lattices, SciPost Phys. Core 4, 005 (2021)

    2021

  17. [25]

    D. A. Hryniuk and M. H. Szyma´ nska, Tensor-network- based variational Monte Carlo approach to the non- equilibrium steady state of open quantum systems, Quantum 8, 1475 (2024)

    2024

  18. [26]

    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)

    2015

  19. [27]

    R. Rota, F. Minganti, C. Ciuti, and V. Savona, Quan- tum critical regime in a quadratically driven nonlinear photonic lattice, Phys. Rev. Lett. 122, 110405 (2019)

    2019

  20. [28]

    Carleo and M

    G. Carleo and M. Troyer, Solving the quantum many- body problem with artificial neural networks, Science 355, 602 (2017)

    2017

  21. [29]

    Carrasquilla and G

    J. Carrasquilla and G. Torlai, How To Use Neural Networks To Investigate Quantum Many-Body Physics, PRX Quantum 2, 040201 (2021)

    2021

  22. [30]

    Vicentini, A

    F. Vicentini, A. Biella, N. Regnault, and C. Ciuti, Varia- tional Neural-Network Ansatz for Steady States in Open Quantum Systems, Phys. Rev. Lett. 122, 250503 (2019)

    2019

  23. [31]

    M. J. Hartmann and G. Carleo, Neural-Network Ap- proach to Dissipative Quantum Many-Body Dynamics, Phys. Rev. Lett. 122, 250502 (2019)

    2019

  24. [32]

    Yoshioka and R

    N. Yoshioka and R. Hamazaki, Constructing neural sta- tionary states for open quantum many-body systems, Phys. Rev. B 99, 214306 (2019)

    2019

  25. [33]

    Nagy and V

    A. Nagy and V. Savona, Variational Quantum Monte Carlo Method with a Neural-Network Ansatz for Open Quantum Systems, Phys. Rev. Lett. 122, 250501 (2019)

    2019

  26. [34]

    Torlai, G

    G. Torlai, G. Mazzola, J. Carrasquilla, M. Troyer, R. Melko, and G. Carleo, Neural-network quantum state tomography, Nat. Phys. 14, 447 (2018)

    2018

  27. [35]

    Carnazza, F

    F. Carnazza, F. Carollo, D. Zietlow, S. Andergassen, G. Martius, and I. Lesanovsky, Inferring Markovian quantum master equations of few-body observables in in- teracting spin chains, New J. Phys. 24, 073033 (2022). 16

    2022

  28. [36]

    D. Luo, Z. Chen, J. Carrasquilla, and B. K. Clark, Au- toregressive Neural Network for Simulating Open Quan- tum Systems via a Probabilistic Formulation, Phys. Rev. Lett. 128, 090501 (2022)

    2022

  29. [37]

    Denis and G

    Z. Denis and G. Carleo, Accurate neural quantum states for interacting lattice bosons, Quantum 9, 1772 (2025)

    2025

  30. [38]

    Polkovnikov, Phase space representation of quantum dynamics, Ann

    A. Polkovnikov, Phase space representation of quantum dynamics, Ann. Phys. (N. Y.) 325, 1790 (2010)

    2010

  31. [39]

    Deuar, A

    P. Deuar, A. Ferrier, M. Matuszewski, G. Orso, and M. H. Szyma´ nska, Fully Quantum Scalable Description of Driven-Dissipative Lattice Models, PRX Quantum 2, 010319 (2021)

    2021

  32. [40]

    C. C. Margossian, A review of automatic differentia- tion and its efficient implementation, WIREs Data Min. Knowl. Discov. 9, e1305 (2019)

    2019

  33. [41]

    Griewank and A

    A. Griewank and A. Walther, Evaluating Derivatives: Principles and Techniques of Algorithmic Differentia- tion, 2nd ed. (SIAM, 2008)

    2008

  34. [42]

    Breuer and F

    H. Breuer and F. Petruccione, The Theory of Open Quan- tum Systems , 1st ed. (Oxford University Press, Great Clarendon Street, Oxford, UK, 2002) p. 648

    2002

  35. [43]

    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)

    1932

  36. [44]

    R. J. Glauber, Coherent and incoherent states of the ra- diation field, Phys. Rev. 131, 2766 (1963)

    1963

  37. [45]

    E. C. G. Sudarshan, Equivalence of semiclassical and quantum mechanical descriptions of statistical light beams, Phys. Rev. Lett. 10, 277 (1963)

    1963

  38. [46]

    Husimi, Some Formal Properties of the Density Ma- trix, Proc

    K. Husimi, Some Formal Properties of the Density Ma- trix, Proc. Phys.-Math. Soc. Jpn., 3rd Ser. 22, 264 (1940)

    1940

  39. [47]

    A. G. Baydin, B. A. Pearlmutter, A. A. Radul, and J. M. Siskind, Automatic differentiation in machine learning: a survey, J. Mach. Learn. Res. 18, 1 (2018)

    2018

  40. [48]

    J. E. Bourassa, N. Quesada, I. Tzitrin, A. Sz´ ava, T. Isac- sson, J. Izaac, K. K. Sabapathy, G. Dauphinais, and I. Dhand, Fast Simulation of Bosonic Qubits via Gaus- sian Functions in Phase Space, PRX Quantum 2, 040315 (2021)

    2021

  41. [49]

    Walschaers, Non-Gaussian Quantum States and Where to Find Them, PRX Quantum 2, 030204 (2021)

    M. Walschaers, Non-Gaussian Quantum States and Where to Find Them, PRX Quantum 2, 030204 (2021)

    2021

  42. [50]

    Leghtas, G

    Z. Leghtas, G. Kirchmair, B. Vlastakis, R. J. Schoelkopf, M. H. Devoret, and M. Mirrahimi, Hardware-efficient au- tonomous quantum memory protection, Phys. Rev. Lett. 111, 120501 (2013)

    2013

  43. [51]

    Yurke and D

    B. Yurke and D. Stoler, Generating quantum mechanical superpositions of macroscopically distinguishable states via amplitude dispersion, Phys. Rev. Lett. 57, 13 (1986)

    1986

  44. [52]

    Casta˜ nos, R

    O. Casta˜ nos, R. J´ auregui, R. L´ opez-Pe˜ na, J. Recamier, and V. I. Man’ko, Schr¨ odinger-cat states in Paul traps, Phys. Rev. A 55, 1208 (1997)

    1997

  45. [53]

    Vogel and H

    K. Vogel and H. Risken, Determination of quasiproba- bility distributions in terms of probability distributions for the rotated quadrature phase, Phys. Rev. A 40, 2847 (1989)

    1989

  46. [54]

    Y.-F. Qu, M. Stefanini, T. Shi, T. Esslinger, S. Gopalakr- ishnan, J. Marino, and E. Demler, Variational approach to the dynamics of dissipative quantum impurity models, Phys. Rev. B 111, 155113 (2025)

    2025

  47. [55]

    Bettencourt, M

    J. Bettencourt, M. J. Johnson, and D. K. Duvenaud, Taylor-Mode Automatic Differentiation for Higher-Order Derivatives in JAX (2019)

    2019

  48. [56]

    Bradbury, R

    J. Bradbury, R. Frostig, P. Hawkins, M. J. Johnson, C. Leary, D. Maclaurin, G. Necula, A. Paszke, J. Van- derPlas, S. Wanderman-Milne, and Q. Zhang, JAX: com- posable transformations of Python+NumPy programs (2018)

    2018

  49. [57]

    X. Yuan, S. Endo, Q. Zhao, Y. Li, and S. C. Benjamin, Theory of variational quantum simulation, Quantum 3, 191 (2019)

    2019

  50. [58]

    Eeltink, F

    D. Eeltink, F. Vicentini, and V. Savona, Variational dy- namics of open quantum systems in phase space (2023), arXiv:2307.07429 [quant-ph]

    2023 arXiv

  51. [59]

    Griewank, J

    A. Griewank, J. Utke, and A. Walther, Evaluating higher derivative tensors by forward propagation of univariate Taylor series, Math. Comput. 69, 1117 (2000)

    2000

  52. [60]

    Dangel, T

    F. Dangel, T. Siebert, M. Zeinhofer, and A. Walther, Col- lapsing Taylor Mode Automatic Differentiation (2025), arXiv:2505.13644 [cs.LG]

    2025

  53. [61]

    Gravina, F

    L. Gravina, F. Minganti, and V. Savona, Critical Schr¨ odinger Cat Qubit, PRX Quantum4, 020337 (2023)

    2023

  54. [62]

    Johansson, P

    J. Johansson, P. Nation, and F. Nori, QuTiP: An open- source Python framework for the dynamics of open quan- tum systems, Comput. Phys. Commun.183, 1760 (2012)

    2012

  55. [63]

    Bartolo, F

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

    2016

  56. [64]

    Chang, V

    C.-H. Chang, V. Dommes, R. Erramilli, A. Homrich, P. Kravchuk, A. Liu, M. Mitchell, D. Poland, and D. Simmons-Duffin, Bootstrapping the 3d Ising stress tensor, J. High Energy Phys. 2025 (3), 136

    2025

  57. [65]

    Reehorst, J

    M. Reehorst, J. High Energy Phys. 2022 (9), 177

    2022

  58. [66]

    Komargodski and D

    Z. Komargodski and D. Simmons-Duffin, The Random- Bond Ising Model in 2.01 and 3 Dimensions, J. Phys. A: Math. Theor. 50, 154001 (2017)

    2017

  59. [67]

    Aspuru-Guzik and P

    A. Aspuru-Guzik and P. Walther, Photonic quantum sim- ulators, Nat. Phys. 8, 285 (2012)

    2012

  60. [68]

    Tosca, M

    J. Tosca, M. C. Strinati, C. Conti, and C. Ciuti, Emergent Equilibrium in All-Optical Single Quantum- Trajectory Ising Machines, Phys. Rev. Lett. 134, 230404 (2025)

    2025

  61. [69]

    M. P. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, Boson localization and the superfluid-insulator transition, Phys. Rev. B 40, 546 (1989)

    1989

  62. [70]

    Cardy, Scaling and Renormalization in Statistical Physics, Cambridge Lecture Notes in Physics (Cam- bridge University Press, 1996)

    J. Cardy, Scaling and Renormalization in Statistical Physics, Cambridge Lecture Notes in Physics (Cam- bridge University Press, 1996)

    1996

  63. [71]

    Minganti, A

    F. Minganti, A. Biella, N. Bartolo, and C. Ciuti, Spectral theory of Liouvillians for dissipative phase transitions, Phys. Rev. A 98, 042118 (2018)

    2018

  64. [72]

    Vicentini, F

    F. Vicentini, F. Minganti, R. Rota, G. Orso, and C. Ciuti, Critical slowing down in driven-dissipative Bose-Hubbard lattices, Phys. Rev. A 97, 013853 (2018)

    2018

  65. [73]

    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)

    2020

  66. [74]

    Adzhemyan, D

    L. Adzhemyan, D. Evdokimov, M. Hnatiˇ c, E. Ivanova, M. Kompaniets, A. Kudlis, and D. Zakharov, The dy- namic critical exponent z for 2d and 3d Ising models from five-loop ϵ expansion, Phys. Lett. A 425, 127870 (2022)

    2022

  67. [75]

    Vogel and H

    K. Vogel and H. Risken, Quasiprobability distributions in dispersive optical bistability, Phys. Rev. A 39, 4675 (1989)

    1989

  68. [76]

    Haroche and J

    S. Haroche and J. M. Raimond, Exploring the Quantum: Atoms, Cavities, and Photons (Oxford Univ. Press, Ox- 17 ford, 2006)

    2006

  69. [77]

    Simon, N

    R. Simon, N. Mukunda, and B. Dutta, Quantum-noise matrix for multimode systems: U(n) invariance, squeez- ing, and normal forms, Phys. Rev. A 49, 1567 (1994)

    1994

  70. [78]

    D. F. Walls and G. J. Milburn, Quantum Optics , 2nd ed. (Springer, 2008)

    2008

Pith tools

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