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:2502.04857 · Physical review. B./Physical review. B

Explicit Pfaffian formula for amplitudes of fermionic Gaussian pure states in arbitrary Pauli bases

M. A. Rajabpour, M. A. Seifi Mirjafarlou, Reyhaneh Khasseh

Published Jun 2, 2025• 2 citations
arXivHTMLPDF

Abstract

The explicit computation of amplitudes for fermionic Gaussian pure states in arbitrary Pauli bases is a long-standing challenge in quantum many-body physics, with significant implications for quantum tomography, experimental studies, and quantum dynamics. These calculations are essential for analyzing complex properties beyond traditional measures, such as formation probabilities, global entanglement, and entropy in nonstandard bases, where exact and computationally efficient methods remain underdeveloped. In addition to these physical applications, having explicit formulas is crucial for optimizing negative log-likelihood functions in quantum tomography, a key task in the NISQ era. In this work, we present an explicit Pfaffian formula (Theorem 1) for determining these amplitudes in arbitrary Pauli bases, utilizing a matrix whose structure reflects the qubit parity. Additionally, we introduce a recursive relation (Theorem 2) that connects amplitudes for systems with varying qubit numbers, enabling scalable computations for large systems. Together, these results provide a versatile framework for studying global entanglement, Shannon-R\'enyi entropies, formation probabilities, and performing efficient quantum tomography, thereby significantly expanding the computational toolkit for analyzing complex quantum systems. Finally, we utilize our formalism to determine the post-measurement entanglement entropy, reflecting how local measurements alter entanglement, and compare the results with conformal field theory predictions.

DOI
10.1103/physrevb.111.235102

Related papers

Matrix Elements of Fermionic Gaussian Operators in Arbitrary Pauli Bases: A Pfaffian Formula

M. A. Rajabpour, M. A. Seifi Mirjafarlou, Reyhaneh Khasseh

Fermionic Gaussian operators are foundational tools in quantum many-body theory, numerical simulation of fermionic dynamics, and fermionic linear optics. While their structure is fully determined by two-point correlations, evaluating their matrix elements in arbitrary local spin bases remains a nontrivial task, especially in applications involving quantum measurements, tomography, and basis-rotated simulations. In this work, we derive a fully explicit and general Pfaffian formula for the matrix elements of fermionic Gaussian operators between arbitrary Pauli product states. Our approach introduces a pair of sign-encoding matrices whose classification leads to a Lie algebra isomorphic to $\mathfrak{so}(2L)$. This algebraic structure not only guarantees consistency of the Pfaffian signs but also reveals deep connections to Clifford algebras. The resulting framework enables scalable computations across diverse fields -- from quantum tomography and entanglement dynamics to algebraic structure in fermionic circuits and matchgate computation. Beyond its practical utility, our construction sheds light on the internal symmetries of Gaussian operators and offers a new lens through which to explore their role in quantum information and computational models.

arXiv:2506.02809•Jun 3, 2025•0 citations

Non-Local Magic Resources for Fermionic Gaussian States

Daniele Iannotti, Beatrice Magni, Riccardo Cioli +2 more

Entanglement and magic are fundamental resources that capture the complexity of quantum many-body systems. Non-local magic isolates the irreducible nonstabilizerness intrinsically tied to entanglement. However, evaluating this quantity generally requires a prohibitive minimization over the full Hilbert space, making it computationally inaccessible beyond a few qubits. Here, we overcome this bottleneck by establishing a closed-form expression for the non-local stabilizer entropies of fermionic Gaussian states over local Gaussian unitaries, which we prove at Rényi index $α=2$ for arbitrary subsystem size, and which can be evaluated in polynomial time directly from the eigenvalues of the reduced Majorana covariance matrix. We apply this framework to characterize fermionic non-local magic across diverse physical regimes: we derive an exact Page-like curve for typical random states, reveal logarithmic scaling at the quantum critical point of the XY model, and establish a quasiparticle picture for magic generation during out-of-equilibrium quantum quenches. Crucially, because our result relies solely on two-point correlation functions, it provides a scalable route for the experimental estimation of fermionic non-local magic in large-scale quantum processors via fermionic shadow tomography.

quant-pharXiv:2604.27049•May 1, 2026•0 citations

Efficient representation of Gaussian fermionic pure states in noncomputational bases

Babak Tarighi, Reyhaneh Khasseh, M. A. Rajabpour

This paper introduces an innovative approach for representing Gaussian fermionic states, pivotal in quantum spin systems and fermionic models, within a range of alternative quantum bases. We focus on transitioning these states from the conventional computational ${\ensuremath{\sigma}}^{z}$ basis to more complex bases, such as $(\ensuremath{\phi},\frac{\ensuremath{\pi}}{2},\ensuremath{\alpha})$, which are essential for accurately calculating critical quantities like formation probabilities and Shannon entropy. Our methodology can be advantageous in quantum quench studies and quantum tomography, where alternative basis representations significantly enhance optimization processes. We present an algorithm that not only simplifies the basis transformation, but also reduces computational complexity, making it feasible to calculate amplitudes of large systems efficiently. Our key contribution is a technique that translates amplitude calculations into the Pfaffian computation of submatrices from an antisymmetric matrix, a process facilitated by understanding domain wall relationships across different bases. As an application, we will determine the formation probabilities for various bases and configurations within the critical transverse field Ising chain, considering both periodic and open boundary conditions. We aim to categorize the configurations and bases by examining the universal constant term that characterizes the scaling of the logarithm of the formation probability in the periodic system, as well as the coefficient of the logarithmic term in the case of open systems. In the open system scenario, this coefficient is influenced by the central charge and the conformal weight of the boundary condition-changing operator. This work is set to expand the toolkit available for researchers in quantum information theory and many-body physics, providing a more efficient and elegant solution for exploring Gaussian fermionic states in nonstandard quantum bases.

arXiv:2403.03289•Physical Review A•Jun 14, 2024•2 citations

Computable fermionic non-Gaussianity from the covariance matrix

Poetri Sonya Tarabunga, Bernhard Jobst, Raúl Morral-Yepes +4 more

Fermionic non-Gaussianity, or fermionic magic, is a key resource underlying the computational complexity of fermionic quantum systems, yet tractable and operationally meaningful ways to quantify it remain limited. We address this challenge by developing a convex resource theory of fermionic non-Gaussianity and introducing two families of computable quantities for pure fermionic states, both derived from the Williamson normal form of the covariance matrix. The first family, occupation number entropies, is defined as the Tsallis-$α$ entropy of the occupation numbers. We prove that one member of this family is monotonic under Gaussian protocols, establishing it as a computable convex resource monotone. It consequently lower bounds the number of non-Gaussian gates needed for state preparation. The second family, natural-orbital participation entropies, is given by the Rényi-$α$ entropy of the squared amplitudes of the state in the natural-orbital basis, defined by the eigenvectors of the covariance matrix. They quantify state compressibility in this basis and thus upper bound the classical simulation cost in an orthonormal Gaussian basis. We analyze both families for stabilizer and translation-invariant states, where they simplify and reveal additional structure. We further study representative examples, including random SWAP-doped matchgate circuits and the bond-modulated XXZ model, highlighting the role of non-Gaussianity in many-body phenomena. Our work establishes a resource-theoretic framework for computable fermionic non-Gaussianity that unifies notions arising across quantum information, condensed-matter physics, and quantum chemistry, opening new directions for studying the complexity of quantum many-body systems and providing practical tools to assess the classical simulability of fermionic states relevant for quantum advantage.

quant-pharXiv:2607.02242•Jul 3, 2026•0 citations

Nonadditivity in fermionic quantum Shannon theory

Farzin Salek, Amir-Reza Negari, Graeme Smith +2 more

We study collective state preparation and channel encoding in fermionic systems under parity superselection and composition governed by the canonical anticommutation relations. For a two-mode Gaussian state $τ$, we obtain the exact entanglement of formation, $E_{\mathrm F}^{\mathrm{phys}}(τ^{\otimes m})=\lceil m/2\rceil\log2$ for every $m\geq1$. Two copies have the same formation value as one, violating additivity and strong superadditivity, while regularization halves the single-copy value. For an associated Gaussian channel, we determine the minimum output entropy and Holevo information at every blocklength and prove nonadditivity of both. Together, these results give exact answers to the fermionic counterparts of the four additivity questions associated with Shor and Hastings. One physical use carries no classical information, whereas two uses transmit one bit perfectly. The unassisted classical and zero-error classical capacities both equal one half bit per mode. A sharp finite-block bound gives a strong converse: the decoding success probability vanishes at every rate above capacity. These optimizations cover all physical pure-state decompositions and channel inputs, with Gaussian ensembles and classical codes attaining the optima. The formation and minimum output entropy formulas hold for every Rényi order $α\in(0,\infty]$. We also determine the unassisted quantum and zero-error quantum capacities, both equal to one half logical qubit per mode, and prove a strong converse: entanglement fidelity vanishes at every rate above capacity. For a noisy Gaussian family, we determine the one- and two-copy formation values and prove strict nonadditivity at full rank. For $τ$, Bell violation under local parity-preserving measurements first appears at three copies.

quant-pharXiv:2610.06758•Oct 6, 2026•0 citations

Explicit Pfaffian formula for amplitudes of fermionic Gaussian pure states in arbitrary Pauli bases

Anonymous

Physical review. B./Physical review. B•May 21, 2025•0 citations
The explicit computation of amplitudes for fermionic Gaussian pure states in arbitrary Pauli bases is a long-standing challenge in quantum many-body physics, with significant implications for quantum tomography, experimental studies, and quantum dynamics. These calculations are essential for analyzing complex properties beyond traditional measures, such as formation probabilities, global entanglement, and entropy in nonstandard bases, where exact and computationally efficient methods remain underdeveloped. In addition to these physical applications, having explicit formulas is crucial for optimizing negative log-likelihood functions in quantum tomography, a key task in the NISQ era. In this work, we present an explicit Pfaffian formula (Theorem 1) for determining these amplitudes in arbitrary Pauli bases, utilizing a matrix whose structure reflects the qubit parity. Additionally, we introduce a recursive relation (Theorem 2) that connects amplitudes for systems with varying qubit numbers, enabling scalable computations for large systems. Together, these results provide a versatile framework for studying global entanglement, Shannon-R'enyi entropies, formation probabilities, and performing efficient quantum tomography, thereby significantly expanding the computational toolkit for analyzing complex quantum systems. Finally, we utilize our formalism to determine the post-measurement entanglement entropy, reflecting how local measurements alter entanglement, and compare the results with conformal field theory predictions.
Fermionic Gaussian operators are foundational tools in quantum many-body theory, numerical simulation of fermionic dynamics, and fermionic linear optics. While their structure is fully determined by two-point correlations, evaluating their matrix elements in arbitrary local spin bases remains a nontrivial task, especially in applications involving quantum measurements, tomography, and basis-rotated simulations. In this work, we derive a fully explicit and general Pfaffian formula for the matrix elements of fermionic Gaussian operators between arbitrary Pauli product states. Our approach introduces a pair of sign-encoding matrices whose classification leads to a Lie algebra isomorphic to so(2L)\mathfrak{so}(2L)so(2L). This algebraic structure not only guarantees consistency of the Pfaffian signs but also reveals deep connections to Clifford algebras. The resulting framework enables scalable computations across diverse fields -- from quantum tomography and entanglement dynamics to algebraic structure in fermionic circuits and matchgate computation. Beyond its practical utility, our construction sheds light on the internal symmetries of Gaussian operators and offers a new lens through which to explore their role in quantum information and computational models.
Entanglement and magic are fundamental resources that capture the complexity of quantum many-body systems. Non-local magic isolates the irreducible nonstabilizerness intrinsically tied to entanglement. However, evaluating this quantity generally requires a prohibitive minimization over the full Hilbert space, making it computationally inaccessible beyond a few qubits. Here, we overcome this bottleneck by establishing a closed-form expression for the non-local stabilizer entropies of fermionic Gaussian states over local Gaussian unitaries, which we prove at Rényi index α=2α=2α=2 for arbitrary subsystem size, and which can be evaluated in polynomial time directly from the eigenvalues of the reduced Majorana covariance matrix. We apply this framework to characterize fermionic non-local magic across diverse physical regimes: we derive an exact Page-like curve for typical random states, reveal logarithmic scaling at the quantum critical point of the XY model, and establish a quasiparticle picture for magic generation during out-of-equilibrium quantum quenches. Crucially, because our result relies solely on two-point correlation functions, it provides a scalable route for the experimental estimation of fermionic non-local magic in large-scale quantum processors via fermionic shadow tomography.
This paper introduces an innovative approach for representing Gaussian fermionic states, pivotal in quantum spin systems and fermionic models, within a range of alternative quantum bases. We focus on transitioning these states from the conventional computational ${\ensuremath{\sigma}}^{z}$ basis to more complex bases, such as $(\ensuremath{\phi},\frac{\ensuremath{\pi}}{2},\ensuremath{\alpha})$, which are essential for accurately calculating critical quantities like formation probabilities and Shannon entropy. Our methodology can be advantageous in quantum quench studies and quantum tomography, where alternative basis representations significantly enhance optimization processes. We present an algorithm that not only simplifies the basis transformation, but also reduces computational complexity, making it feasible to calculate amplitudes of large systems efficiently. Our key contribution is a technique that translates amplitude calculations into the Pfaffian computation of submatrices from an antisymmetric matrix, a process facilitated by understanding domain wall relationships across different bases. As an application, we will determine the formation probabilities for various bases and configurations within the critical transverse field Ising chain, considering both periodic and open boundary conditions. We aim to categorize the configurations and bases by examining the universal constant term that characterizes the scaling of the logarithm of the formation probability in the periodic system, as well as the coefficient of the logarithmic term in the case of open systems. In the open system scenario, this coefficient is influenced by the central charge and the conformal weight of the boundary condition-changing operator. This work is set to expand the toolkit available for researchers in quantum information theory and many-body physics, providing a more efficient and elegant solution for exploring Gaussian fermionic states in nonstandard quantum bases.
Fermionic non-Gaussianity, or fermionic magic, is a key resource underlying the computational complexity of fermionic quantum systems, yet tractable and operationally meaningful ways to quantify it remain limited. We address this challenge by developing a convex resource theory of fermionic non-Gaussianity and introducing two families of computable quantities for pure fermionic states, both derived from the Williamson normal form of the covariance matrix. The first family, occupation number entropies, is defined as the Tsallis-ααα entropy of the occupation numbers. We prove that one member of this family is monotonic under Gaussian protocols, establishing it as a computable convex resource monotone. It consequently lower bounds the number of non-Gaussian gates needed for state preparation. The second family, natural-orbital participation entropies, is given by the Rényi-ααα entropy of the squared amplitudes of the state in the natural-orbital basis, defined by the eigenvectors of the covariance matrix. They quantify state compressibility in this basis and thus upper bound the classical simulation cost in an orthonormal Gaussian basis. We analyze both families for stabilizer and translation-invariant states, where they simplify and reveal additional structure. We further study representative examples, including random SWAP-doped matchgate circuits and the bond-modulated XXZ model, highlighting the role of non-Gaussianity in many-body phenomena. Our work establishes a resource-theoretic framework for computable fermionic non-Gaussianity that unifies notions arising across quantum information, condensed-matter physics, and quantum chemistry, opening new directions for studying the complexity of quantum many-body systems and providing practical tools to assess the classical simulability of fermionic states relevant for quantum advantage.
We study collective state preparation and channel encoding in fermionic systems under parity superselection and composition governed by the canonical anticommutation relations. For a two-mode Gaussian state τττ, we obtain the exact entanglement of formation, EFphys(τ⊗m)=⌈m/2⌉log⁡2E_{\mathrm F}^{\mathrm{phys}}(τ^{\otimes m})=\lceil m/2\rceil\log2EFphys​(τ⊗m)=⌈m/2⌉log2 for every m≥1m\geq1m≥1. Two copies have the same formation value as one, violating additivity and strong superadditivity, while regularization halves the single-copy value. For an associated Gaussian channel, we determine the minimum output entropy and Holevo information at every blocklength and prove nonadditivity of both. Together, these results give exact answers to the fermionic counterparts of the four additivity questions associated with Shor and Hastings. One physical use carries no classical information, whereas two uses transmit one bit perfectly. The unassisted classical and zero-error classical capacities both equal one half bit per mode. A sharp finite-block bound gives a strong converse: the decoding success probability vanishes at every rate above capacity. These optimizations cover all physical pure-state decompositions and channel inputs, with Gaussian ensembles and classical codes attaining the optima. The formation and minimum output entropy formulas hold for every Rényi order α∈(0,∞]α\in(0,\infty]α∈(0,∞]. We also determine the unassisted quantum and zero-error quantum capacities, both equal to one half logical qubit per mode, and prove a strong converse: entanglement fidelity vanishes at every rate above capacity. For a noisy Gaussian family, we determine the one- and two-copy formation values and prove strict nonadditivity at full rank. For τττ, Bell violation under local parity-preserving measurements first appears at three copies.