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:2507.16100

Derivation of the Loop Hafnian Generating Function for Arbitrary Symmetric Matrices via Gaussian Integration

Sergey V. Tarasov

Published Jul 21, 2025• 0 citations
arXivHTMLPDF

Abstract

This short note shows that the recently proposed generating function for loop hafnians -- originally derived using quantum-optical methods for a restricted class of matrices -- is in fact valid for arbitrary symmetric matrices. The proof relies solely on Gaussian integration and does not assume any additional properties inherited from the covariance matrices of quantum Gaussian states.

DOI
10.48550/arxiv.2507.16100

Related papers

Klauder-Bargmann Integral Representation of Gaussian Symmetries and Generating Functions of Gaussian States

Tiju Cherian John, Kalyanapuram Rangachari Parthasarathy

Let $\Gamma(\mathcal{H})$ be the boson Fock space over a finite dimensional Hilbert space $\mathcal{H}$. It is shown that every gaussian symmetry admits a Klauder-Bargmann integral representation in terms of coherent states. Furthermore, gaussian symmetries, gaussian states and second quantization contractions, all of these operators belong to a weakly closed, selfadjoint semigroup $\mathcal{E}_2(\mathcal{H})$ of bounded operators in $\Gamma(\mathcal{H})$. This yields, a new parametrization of gaussian states, which is a very fruitful alternative to the customary parametrization by position-momentum mean vectors and covariance matrices. This leads to a rich harvest of corollaries: (i) every gaussian state $\rho$ admits a factorization $\rho = Z_{1}^{\dagger}Z_{1}$, where $Z_{1}$ is an element of $\mathcal{E}_2(\mathcal{H})$ and has the form $Z_{1} = \sqrt{c}\Gamma(\sqrt{\Lambda})\exp{\sum_{r=1}^{n} \lambda_ra_r+\sum_{r,s=1}^{n} \alpha_{rs}a_{r}a_{s}}$ on the dense linear manifold generated by all exponential vectors, $\Lambda$ being a positive operator in $\mathcal{H}$, $a_{r}, 1\leq r \leq n$ are the annihilation operators corresponding to the $n$ different modes in $\Gamma(\mathcal{H})$, $\lambda_r\in \mathbb{C}$ and $[\alpha_{rs}]$ is a symmetric matrix in $M_n(\mathbb{C})$; (ii) an explicit particle basis expansion of an arbitrary mean zero pure gaussian state vector along with a density matrix formula for a general gaussian state in terms of its $\mathcal{E}_2(\mathcal{H})$-parameters; (iii) an easy test for the entanglement of pure gaussian states and a class of examples of pure $n$-mode gaussian states which are completely entangled; (iv) tomography of an unknown gaussian state in $\Gamma(\mathbb{C}^n)$ by the estimation of its $\mathcal{E}_2(\mathbb{C}^n)$-parameters using $O(n^2)$ measurements with a finite number of outcomes.

•Nov 15, 2019•0 citations

Tutorial Notes on One-Party and Two-Party Gaussian States

Berthold‐Georg Englert, Krzysztof Wódkiewicz

Gaussian states — or, more generally, Gaussian operators — play an important role in Quantum Optics and Quantum Information Science, both in discussions about conceptual issues and in practical applications. We describe, in a tutorial manner, a systematic operator method for first characterizing such states and then investigating their properties. The central numerical quantities are the covariance matrix that specifies the characteristic function of the state, and the closely related matrices associated with Wigner's and Glauber's phase space functions. For pedagogical reasons, we restrict the discussion to one-dimensional and two-dimensional Gaussian states, for which we provide illustrating and instructive examples.

arXiv:quant-ph/0307196•International Journal of Quantum Information•Jun 1, 2003•44 citations

Combinatorial Matrix Functions in Ising Models and Quantum Circuit Representation

Minhyeok Kang, Gwonhak Lee, Youngrong Lim +1 more

Permanents, hafnians, and loop-hafnians are combinatorial matrix functions defined on the perfect matchings of graphs. These matrix functions arise in the quantum amplitudes for output photon configurations in linear optical networks, which has led both to sampling hardness results and to a variety of applications. The same matrix functions also appear in Ising models, since Ising models naturally encode graphs, and the matrix functions emerge from quantum interference. However, previous results have been restricted to balanced bipartite graphs, and hence to the permanent. In this paper, we extend the graph structure of the Ising model so that all three matrix functions can be obtained. In particular, self-loops can be encoded by introducing additional spins, and the loop-hafnian is obtained from a superposition of spin states, which we call the loop-hafnian state. We show how to recover the hafnian and permanent cases as the number of self-loops decreases. Moreover, we show that the Ising dynamics we consider can be simulated on a quantum circuit with \(O(n^2)\) gates for a \(2n \times 2n\) real symmetric matrix. Specifically, we propose a quantum circuit that maps the loop-hafnian state to a basis state with \(O(n^2)\) gates, and show that, by an adaptive measurement, the loop-hafnian states for all subsets of vertices can be measured simultaneously.

quant-pharXiv:2602.01868•Feb 3, 2026•0 citations

Generation of Matrices with Specified Eigenvalues and their Diagonalization

Habatwa Vincent Mweene

We present a prescription for forming matrices with specified eigenvalues and known eigenvectors. With this method, we can form Hermitian, anti-Hermitian, symmetric and general matrices with arbitrary eigenvalues. In addition we propose an algorithm for diagonalizing such matrices. The functions required for the realization of this are probability amplitudes connecting observables with discrete eigenvalue spectra and can be obtained from spin theory. For the example case of $5\times 5$ matrices, these functions are given.

arXiv:quant-ph/0506074•Jun 9, 2005•0 citations

Representation of the three-dimensional rotation operator in the anisotropic Gaussian basis

R. G. Nazmitdinov, Luis Miguel Robledo

A method based on the generating function technique is proposed to calculate matrix elements of the three-dimensional rotation operator for anisotropic Gaussian wavefunctions.

Journal of Physics A Mathematical and General•Feb 21, 1997•0 citations

Generating function of product of bivariate Hermite polynomials and their applications in studying quantum optical states

Hong-Yi Fan, Pengfei Zhang, Zhen Wang

By virtue of the operator-Hermite-polynomial method, we derive some new generating function formulas of the product of two bivariate Hermite polynomials. Their applications in studying quantum optical states are presented.

Chinese Physics B•Apr 30, 2015•6 citations
This short note shows that the recently proposed generating function for loop hafnians -- originally derived using quantum-optical methods for a restricted class of matrices -- is in fact valid for arbitrary symmetric matrices. The proof relies solely on Gaussian integration and does not assume any additional properties inherited from the covariance matrices of quantum Gaussian states.
Let Γ(H)\Gamma(\mathcal{H})Γ(H) be the boson Fock space over a finite dimensional Hilbert space H\mathcal{H}H. It is shown that every gaussian symmetry admits a Klauder-Bargmann integral representation in terms of coherent states. Furthermore, gaussian symmetries, gaussian states and second quantization contractions, all of these operators belong to a weakly closed, selfadjoint semigroup E2(H)\mathcal{E}_2(\mathcal{H})E2​(H) of bounded operators in Γ(H)\Gamma(\mathcal{H})Γ(H). This yields, a new parametrization of gaussian states, which is a very fruitful alternative to the customary parametrization by position-momentum mean vectors and covariance matrices. This leads to a rich harvest of corollaries: (i) every gaussian state ρ\rhoρ admits a factorization ρ=Z1†Z1\rho = Z_{1}^{\dagger}Z_{1}ρ=Z1†​Z1​, where Z1Z_{1}Z1​ is an element of E2(H)\mathcal{E}_2(\mathcal{H})E2​(H) and has the form Z1=cΓ(Λ)exp⁡∑r=1nλrar+∑r,s=1nαrsarasZ_{1} = \sqrt{c}\Gamma(\sqrt{\Lambda})\exp{\sum_{r=1}^{n} \lambda_ra_r+\sum_{r,s=1}^{n} \alpha_{rs}a_{r}a_{s}}Z1​=c​Γ(Λ​)exp∑r=1n​λr​ar​+∑r,s=1n​αrs​ar​as​ on the dense linear manifold generated by all exponential vectors, Λ\LambdaΛ being a positive operator in H\mathcal{H}H, ar,1≤r≤na_{r}, 1\leq r \leq nar​,1≤r≤n are the annihilation operators corresponding to the nnn different modes in Γ(H)\Gamma(\mathcal{H})Γ(H), λr∈C\lambda_r\in \mathbb{C}λr​∈C and [αrs][\alpha_{rs}][αrs​] is a symmetric matrix in Mn(C)M_n(\mathbb{C})Mn​(C); (ii) an explicit particle basis expansion of an arbitrary mean zero pure gaussian state vector along with a density matrix formula for a general gaussian state in terms of its E2(H)\mathcal{E}_2(\mathcal{H})E2​(H)-parameters; (iii) an easy test for the entanglement of pure gaussian states and a class of examples of pure nnn-mode gaussian states which are completely entangled; (iv) tomography of an unknown gaussian state in Γ(Cn)\Gamma(\mathbb{C}^n)Γ(Cn) by the estimation of its E2(Cn)\mathcal{E}_2(\mathbb{C}^n)E2​(Cn)-parameters using O(n2)O(n^2)O(n2) measurements with a finite number of outcomes.
Gaussian states — or, more generally, Gaussian operators — play an important role in Quantum Optics and Quantum Information Science, both in discussions about conceptual issues and in practical applications. We describe, in a tutorial manner, a systematic operator method for first characterizing such states and then investigating their properties. The central numerical quantities are the covariance matrix that specifies the characteristic function of the state, and the closely related matrices associated with Wigner's and Glauber's phase space functions. For pedagogical reasons, we restrict the discussion to one-dimensional and two-dimensional Gaussian states, for which we provide illustrating and instructive examples.
Permanents, hafnians, and loop-hafnians are combinatorial matrix functions defined on the perfect matchings of graphs. These matrix functions arise in the quantum amplitudes for output photon configurations in linear optical networks, which has led both to sampling hardness results and to a variety of applications. The same matrix functions also appear in Ising models, since Ising models naturally encode graphs, and the matrix functions emerge from quantum interference. However, previous results have been restricted to balanced bipartite graphs, and hence to the permanent. In this paper, we extend the graph structure of the Ising model so that all three matrix functions can be obtained. In particular, self-loops can be encoded by introducing additional spins, and the loop-hafnian is obtained from a superposition of spin states, which we call the loop-hafnian state. We show how to recover the hafnian and permanent cases as the number of self-loops decreases. Moreover, we show that the Ising dynamics we consider can be simulated on a quantum circuit with O(n2)O(n^2)O(n2) gates for a 2n×2n2n \times 2n2n×2n real symmetric matrix. Specifically, we propose a quantum circuit that maps the loop-hafnian state to a basis state with O(n2)O(n^2)O(n2) gates, and show that, by an adaptive measurement, the loop-hafnian states for all subsets of vertices can be measured simultaneously.
We present a prescription for forming matrices with specified eigenvalues and known eigenvectors. With this method, we can form Hermitian, anti-Hermitian, symmetric and general matrices with arbitrary eigenvalues. In addition we propose an algorithm for diagonalizing such matrices. The functions required for the realization of this are probability amplitudes connecting observables with discrete eigenvalue spectra and can be obtained from spin theory. For the example case of 5×55\times 55×5 matrices, these functions are given.