Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2610.33613

The Coherent State Projection Eigenvalue Theorem

Peter D. Drummond, Margaret D. Reid

quant-phSubmitted Oct 24, 2026• Announced Oct 27, 2026• 0 citations
arXivHTMLPDF

Abstract

Recent papers have developed a cyclic symmetry projected extension for coherent state representations of bosonic Hilbert space. This extension can greatly reduce sampling errors in phase-space simulations of bosonic quantum systems that have symmetries. Here we prove an eigenvalue theorem for the continuum limit of the projected coherent states. Our intuitively appealing result is that in the continuum limit the projected coherent states are eigenstates of the symmetry generating operator. This result gives a differential identity in phase-space that includes the conserved eigenvalues. Detailed applications of these differential identities to obtaining improved Fokker-Planck and stochastic equations in many-body physics will be given elsewhere.

AI Summary

The paper asks whether cyclic symmetry projections of bosonic coherent states converge, as the cyclic order increases, to eigenstates of the corresponding continuous-symmetry generators. It considers finitely many bosonic modes and mutually commuting generators $\hat g_\nu=\hat{\boldsymbol a}^\dagger g_\nu\hat{\boldsymbol a}$, where the one-particle matrices are Hermitian with integer eigenvalues. For a fixed integer sector label $p$, it proves that $\Pi_p^{(\mathcal M)}\lVert\boldsymbol\alpha\rangle\to E_p\lVert\boldsymbol\alpha\rangle$ in Hilbert-space norm, where $E_p$ is the generator’s spectral projector. It also establishes $\|(\hat g-p)\Pi_p^{(\mathcal M)}\lVert\boldsymbol\alpha\rangle\|\to0$. If the limiting projection is nonzero, normalized projected states converge to its normalization. The proof uses the exact identity $\Pi_p^{(\mathcal M)}=\sum_{\ell\in\mathbb Z}E_{p+\ell\mathcal M}$: unwanted sectors lie at $|m|\ge\mathcal M-|p|$, so their total weight vanishes as a spectral tail. Finite coherent-state generator moments establish the corresponding convergence for the unbounded generator. The argument extends to multiple commuting projections and limits involving only a subset of them. The limiting eigenvalue relations yield first- and second-order differential identities intended to modify drift and diffusion in phase-space simulations. The convergence is strong, not generally in operator norm; generators, amplitudes, and sector labels remain fixed. The construction assumes commuting integer-spectrum one-particle generators, and a finite lattice’s discrete translation symmetry does not itself admit the continuum limit used here. Explicit improved stochastic equations, numerical demonstrations, and sampling-error reductions are deferred to other work.

Little sign of AI-generated textIncremental

Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.

Related papers

Coherent states in projected Hilbert spaces

P. D. Drummond, M. D. Reid

Coherent states in a projected Hilbert space have many useful properties. When there are conserved quantities, a representation of the entire Hilbert space is not necessary. The same issue arises when conditional observations are made with postselected measurement results. In these cases, only a part of the Hilbert space needs to be represented, and one can define this restriction by way of a projection operator. Here coherent state bases and normally ordered phase-space representations are introduced for treating such projected Hilbert spaces, including existence theorems and dynamical equations. These techniques are very useful in studying novel optical or microwave integrated photonic quantum technologies, such as boson sampling or Josephson quantum computers. In these cases, states become strongly restricted due to inputs, nonlinearities, or conditional measurements. This paper focuses on coherent phase states, which have especially simple properties. Practical applications are reported on calculating recurrences in anharmonic oscillators, the effects of arbitrary phase noise on Schr\"odinger cat fringe visibility, and on boson sampling interferometry for large numbers of modes.

arXiv:1610.01261•Physical Review A•Dec 23, 2016•8 citations

Matrix phase-space representations for quantum symmetries

Peter D. Drummond, Alexander S. Dellios, Margaret D. Reid

We introduce a general phase-space representation that includes global quantum symmetries in the basis expansion. This method, called matrix phase-space, projects the basis onto a reduced Hilbert space, which can greatly reduce sampling errors of many-body quantum simulations and unifies several previous phase-space methods. The purpose of this paper is to provide detailed proofs of basic theorems and operator identities. We also treat several different types of symmetries. To illustrate the benefits of matrix phase-space methods, we give a detailed derivation of a recent application to the topical problem of verifying the outputs of Gaussian boson sampling (GBS) quantum computers with photon number resolving detectors. This has exponential complexity, and using parity symmetry reduces sampling errors by very large factors relative to earlier methods.

quant-pharXiv:2606.12769•Jun 12, 2026•0 citations

Angular momentum projection integrals for coherent states

Serdar Kuyucak B, K. Unnikrishnan

Angular momentum projection integrals arising from the coherent state representation of bosonic spectrum generating algebras are evaluated using analytical and algebraic methods. For the vibron model of diatomic molecules, a closed expression is derived by direct integration. For the general case, the integral is evaluated as a 1/N expansion by exploiting the symmetries of the underlying boson system and making use of computer algebra.

Journal of Physics A Mathematical and General•Apr 7, 1995•5 citations

Bosonic Coherent States

Didier Robert, Monique Combescure

In a first part we give a brief presentation of general Fock space setting to describe quantum field theory. Bosons are quantum particles with integer spin and have symmetric wave functions; Fermions are quantum particles with half-integer spin and are represented by anti-symmetric wave functions. The functional setting is given by a symmetric or anti-symmetric tensor product of Hilbert spaces. We describe these spaces and transformations between these spaces. We shall follow the references [23, 42]. Coherent states are defined by translating the vacuum states with the Weyl operators. This is easily done here for bosons. We shall see in the next chapter how to deal with fermions. In a second part we give an interesting application of bosonic coherent states to the study of the classical limit as \(\hbar \searrow 0\) of non-relativistic boson systems with two body interaction in the neighborhood of a solution of the classical system (here the Hartree equation). The classical limit corresponds here to the mean-field limit as the number of particles goes to infinity. As we have done for finite systems, we use here Hepp’s method, which is a linearization procedure of the quantum Hamiltonian around the classical field. The fluctuations around this solution are controlled by a purely quadratic Hamiltonian. In a series of several important papers [114–116] Ginibre and Velo have proven an asymptotic expansion and remainder estimates for these quantum fluctuations. Finally, following the paper [225] one can show that, in the limit \(\hbar \rightarrow 0\), the marginal distribution of the time-evolved coherent states tends in trace-norm to the projector onto the solution of the classical fields equation (Hartree equation) with a uniform remainder estimates in time.

Theoretical and mathematical physics•Jan 1, 2021•0 citations

Path-integral formulation of truncated Wigner approximation for bosonic Markovian open quantum systems

Toma Yoneya, Kazuya Fujimoto, Yuki Kawaguchi

The truncated Wigner approximation (TWA) enables us to investigate bosonic quantum many-body dynamics, including open quantum systems described by the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation. In the TWA, the Weyl-Wigner transformation, a way of mapping from quantum-mechanical operators to c -numbers, of the GKSL equation leads to the Fokker–Planck equation, which we calculate by reducing it to the corresponding stochastic differential equations. However, the Fokker–Planck equation is not always reduced to the stochastic differential equations depending on details of jump operators. In this work, we clarify the condition for obtaining the stochastic differential equations from the Fokker–Planck equation and derive analytical expressions of these equations for a system with an arbitrary Hamiltonian with jump operators that do not couple different states. This result leads to remarkable simplification of the complicated calculations required in the conventional procedures of the TWA. In the course of the derivation, we formulate the GKSL equation by using the path-integral representation based on the Weyl-Wigner transformation, which gives us a clear interpretation of the relation between the TWA and quantum fluctuations and allows us to calculate the non-equal time correlation functions in the TWA. In the benchmark calculations, we numerically confirm that the relaxation dynamics of physical quantities including the non-equal time correlation functions obtained in our formulation agrees well with the exact ones in the numerically solvable models. Since our formulation is widely applicable to bosonic systems such as the Bose-Hubbard model, it facilitates studies of the open quantum many-body phenomenon.

arXiv:2405.11173•Annals of Physics•May 16, 2025•5 citations

The Quantum and Stochastic Toolbox: xSPDE4.2

P. D. Drummond, Teh, Run Yan, Thenabadu, Manushan +5 more

This is the fourth major release of the xSPDE toolbox, which solves stochastic partial and ordinary differential equations, with applications in biology, chemistry, engineering, medicine, physics and quantum technologies. It computes statistical averages, including time-step and sampling error estimation. xSPDE can provide higher order convergence, Fourier spectra and probability densities. The toolbox has graphical output and $χ^{2}$ statistics, as well as weighted, projected, or forward-backward equations. It can generate input-output quantum spectra. The equations can have independent periodic, Dirichlet, and Neumann or Robin boundary conditions in any dimension, for any vector component, and at either end of any interval. xSPDE has functions that can numerically solve both ordinary and partial differential stochastic equations of any type, obtaining correlations, probabilities and averages. The toolbox has a core treating stochastic differential equations, with averages, probability distributions and full error estimates. There are stochastic extensions treating applications to partial differential equations, projected equations, quantum stochastic equations, master equations and quantum phase-space simulations including Gaussian boson sampling experiments.

arXiv:2303.04448•Mar 8, 2023•0 citations
Recent papers have developed a cyclic symmetry projected extension for coherent state representations of bosonic Hilbert space. This extension can greatly reduce sampling errors in phase-space simulations of bosonic quantum systems that have symmetries. Here we prove an eigenvalue theorem for the continuum limit of the projected coherent states. Our intuitively appealing result is that in the continuum limit the projected coherent states are eigenstates of the symmetry generating operator. This result gives a differential identity in phase-space that includes the conserved eigenvalues. Detailed applications of these differential identities to obtaining improved Fokker-Planck and stochastic equations in many-body physics will be given elsewhere.
The paper asks whether cyclic symmetry projections of bosonic coherent states converge, as the cyclic order increases, to eigenstates of the corresponding continuous-symmetry generators. It considers finitely many bosonic modes and mutually commuting generators g^ν=a^†gνa^\hat g_\nu=\hat{\boldsymbol a}^\dagger g_\nu\hat{\boldsymbol a}g^​ν​=a^†gν​a^, where the one-particle matrices are Hermitian with integer eigenvalues.
For a fixed integer sector label ppp, it proves that Πp(M)∥α⟩→Ep∥α⟩\Pi_p^{(\mathcal M)}\lVert\boldsymbol\alpha\rangle\to E_p\lVert\boldsymbol\alpha\rangleΠp(M)​∥α⟩→Ep​∥α⟩ in Hilbert-space norm, where EpE_pEp​ is the generator’s spectral projector. It also establishes ∥(g^−p)Πp(M)∥α⟩∥→0\|(\hat g-p)\Pi_p^{(\mathcal M)}\lVert\boldsymbol\alpha\rangle\|\to0∥(g^​−p)Πp(M)​∥α⟩∥→0. If the limiting projection is nonzero, normalized projected states converge to its normalization. The proof uses the exact identity Πp(M)=∑ℓ∈ZEp+ℓM\Pi_p^{(\mathcal M)}=\sum_{\ell\in\mathbb Z}E_{p+\ell\mathcal M}Πp(M)​=∑ℓ∈Z​Ep+ℓM​: unwanted sectors lie at ∣m∣≥M−∣p∣|m|\ge\mathcal M-|p|∣m∣≥M−∣p∣, so their total weight vanishes as a spectral tail. Finite coherent-state generator moments establish the corresponding convergence for the unbounded generator. The argument extends to multiple commuting projections and limits involving only a subset of them.
The limiting eigenvalue relations yield first- and second-order differential identities intended to modify drift and diffusion in phase-space simulations. The convergence is strong, not generally in operator norm; generators, amplitudes, and sector labels remain fixed. The construction assumes commuting integer-spectrum one-particle generators, and a finite lattice’s discrete translation symmetry does not itself admit the continuum limit used here. Explicit improved stochastic equations, numerical demonstrations, and sampling-error reductions are deferred to other work.
Coherent states in a projected Hilbert space have many useful properties. When there are conserved quantities, a representation of the entire Hilbert space is not necessary. The same issue arises when conditional observations are made with postselected measurement results. In these cases, only a part of the Hilbert space needs to be represented, and one can define this restriction by way of a projection operator. Here coherent state bases and normally ordered phase-space representations are introduced for treating such projected Hilbert spaces, including existence theorems and dynamical equations. These techniques are very useful in studying novel optical or microwave integrated photonic quantum technologies, such as boson sampling or Josephson quantum computers. In these cases, states become strongly restricted due to inputs, nonlinearities, or conditional measurements. This paper focuses on coherent phase states, which have especially simple properties. Practical applications are reported on calculating recurrences in anharmonic oscillators, the effects of arbitrary phase noise on Schr\"odinger cat fringe visibility, and on boson sampling interferometry for large numbers of modes.
We introduce a general phase-space representation that includes global quantum symmetries in the basis expansion. This method, called matrix phase-space, projects the basis onto a reduced Hilbert space, which can greatly reduce sampling errors of many-body quantum simulations and unifies several previous phase-space methods. The purpose of this paper is to provide detailed proofs of basic theorems and operator identities. We also treat several different types of symmetries. To illustrate the benefits of matrix phase-space methods, we give a detailed derivation of a recent application to the topical problem of verifying the outputs of Gaussian boson sampling (GBS) quantum computers with photon number resolving detectors. This has exponential complexity, and using parity symmetry reduces sampling errors by very large factors relative to earlier methods.
In a first part we give a brief presentation of general Fock space setting to describe quantum field theory. Bosons are quantum particles with integer spin and have symmetric wave functions; Fermions are quantum particles with half-integer spin and are represented by anti-symmetric wave functions. The functional setting is given by a symmetric or anti-symmetric tensor product of Hilbert spaces. We describe these spaces and transformations between these spaces. We shall follow the references [23, 42]. Coherent states are defined by translating the vacuum states with the Weyl operators. This is easily done here for bosons. We shall see in the next chapter how to deal with fermions. In a second part we give an interesting application of bosonic coherent states to the study of the classical limit as ℏ↘0\hbar \searrow 0ℏ↘0 of non-relativistic boson systems with two body interaction in the neighborhood of a solution of the classical system (here the Hartree equation). The classical limit corresponds here to the mean-field limit as the number of particles goes to infinity. As we have done for finite systems, we use here Hepp’s method, which is a linearization procedure of the quantum Hamiltonian around the classical field. The fluctuations around this solution are controlled by a purely quadratic Hamiltonian. In a series of several important papers [114–116] Ginibre and Velo have proven an asymptotic expansion and remainder estimates for these quantum fluctuations. Finally, following the paper [225] one can show that, in the limit ℏ→0\hbar \rightarrow 0ℏ→0, the marginal distribution of the time-evolved coherent states tends in trace-norm to the projector onto the solution of the classical fields equation (Hartree equation) with a uniform remainder estimates in time.
The truncated Wigner approximation (TWA) enables us to investigate bosonic quantum many-body dynamics, including open quantum systems described by the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation. In the TWA, the Weyl-Wigner transformation, a way of mapping from quantum-mechanical operators to c -numbers, of the GKSL equation leads to the Fokker–Planck equation, which we calculate by reducing it to the corresponding stochastic differential equations. However, the Fokker–Planck equation is not always reduced to the stochastic differential equations depending on details of jump operators. In this work, we clarify the condition for obtaining the stochastic differential equations from the Fokker–Planck equation and derive analytical expressions of these equations for a system with an arbitrary Hamiltonian with jump operators that do not couple different states. This result leads to remarkable simplification of the complicated calculations required in the conventional procedures of the TWA. In the course of the derivation, we formulate the GKSL equation by using the path-integral representation based on the Weyl-Wigner transformation, which gives us a clear interpretation of the relation between the TWA and quantum fluctuations and allows us to calculate the non-equal time correlation functions in the TWA. In the benchmark calculations, we numerically confirm that the relaxation dynamics of physical quantities including the non-equal time correlation functions obtained in our formulation agrees well with the exact ones in the numerically solvable models. Since our formulation is widely applicable to bosonic systems such as the Bose-Hubbard model, it facilitates studies of the open quantum many-body phenomenon.
This is the fourth major release of the xSPDE toolbox, which solves stochastic partial and ordinary differential equations, with applications in biology, chemistry, engineering, medicine, physics and quantum technologies. It computes statistical averages, including time-step and sampling error estimation. xSPDE can provide higher order convergence, Fourier spectra and probability densities. The toolbox has graphical output and χ2χ^{2}χ2 statistics, as well as weighted, projected, or forward-backward equations. It can generate input-output quantum spectra. The equations can have independent periodic, Dirichlet, and Neumann or Robin boundary conditions in any dimension, for any vector component, and at either end of any interval. xSPDE has functions that can numerically solve both ordinary and partial differential stochastic equations of any type, obtaining correlations, probabilities and averages. The toolbox has a core treating stochastic differential equations, with averages, probability distributions and full error estimates. There are stochastic extensions treating applications to partial differential equations, projected equations, quantum stochastic equations, master equations and quantum phase-space simulations including Gaussian boson sampling experiments.