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.
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). 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.
In quantum information theory and statistical physics, symmetries of multiple copies, or replicas, of a system play a pivotal role. For unitary ensembles, these symmetries are encoded in the replicated commutant: the algebra of operators commuting with the ensemble across k replicas. Determining the commutant is straightforward for the full unitary group, but remains a major obstacle for structured, computationally relevant circuit families. We solve this problem for matchgate circuits, which prepare fermionic Gaussian states on n qubits. Using a Majorana fermion representation, we show that operators coupling different system copies generate the orthogonal Lie algebra so(k), endowing the space of invariants with rich and tractable structure. This underlying symmetry decomposes the matchgate commutant into irreducible sectors, which we completely resolve via a Gelfand--Tsetlin construction. We provide an explicit orthonormal basis of the matchgate commutant for all k and n, together with a formula for its dimension that grows polynomially in n. Furthermore, we characterize the commutant of the Clifford--matchgate subgroup, showing that restricting to signed permutations of Majorana modes yields a commutant that qualitatively diverges from the matchgate case for k≥4 replicas. Ultimately, our orthonormal basis turns algebraic classification into a working toolbox. Using it, we derive closed-form expressions for matchgate twirling channels and a fermionic analogue of Weingarten calculus, the projector encoding all moments of the Gaussian state orbit, state and unitary frame potentials, the average nonstabilizerness of fermionic Gaussian states, a systematic hierarchy of non-Gaussianity measures, and a fermionic de Finetti theorem.
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.
Lie groups, and therefore Lie algebras, are fundamental structures in quantum physics that determine the space of possible trajectories of evolving systems. However, classification and characterization methods for these structures are often impractical for larger systems. In this work, we provide a comprehensive classification of Lie algebras generated by an arbitrary set of Pauli operators, from which an efficient method to characterize them follows. By mapping the problem to a graph setting, we identify a reduced set of equivalence classes: the free-fermionic Lie algebra, the set of all anti-symmetric Paulis on n qubits, the Lie algebra of symplectic Paulis on n qubits, and the space of all Pauli operators on n qubits, as well as controlled versions thereof. Moreover, out of these, we distinguish 6 Clifford inequivalent cases and find a simple set of canonical operators for each, which allow us to give a physical interpretation of the dynamics of each class. Our findings reveal a no-go result for the existence of small Lie algebras beyond the free-fermionic case in the Pauli setting and offer efficiently computable criteria for universality and extendibility of gate sets. These results bear significant impact in ideas in a number of fields like quantum control, quantum machine learning, or classical simulation of quantum circuits.
In this work, we characterize the t-th order commutants of fermionic Gaussian unitaries and of their particle-preserving subgroup acting on n fermionic modes. These commutants govern Haar averages over the corresponding groups and therefore play a central role in fermionic randomized protocols, invariant theory, and resource quantification. Using Howe dualities, we show that the particle-preserving commutant is generated by generalized copy-hopping operators, while that for general Gaussian commutant is generated by generalized quadratic Majorana bilinears together with parity. We then derive closed formulas for the dimensions of both commutants as functions of t and n, and develop constructive Gelfand--Tsetlin procedures to obtain explicit orthonormal bases, with detailed low-t examples. Our framework also clarifies the structure of replicated fermionic states and connects naturally to measures of fermionic correlations, generalized Plücker-type constraints, and the stabilizer entropy of fermionic Gaussian states. These results provide a unified algebraic description of higher-order invariants for fermionic Gaussian dynamics.
Understanding the dynamical properties of quantum systems is an essential task in quantum computing, quantum control, and many-body physics. Tools such as representation theory and Lie theory provide crucial information on reachability and computational power. However, this information can be difficult to access exactly or compute efficiently for arbitrary generating sets. Here we focus on the setting of Pauli strings, which satisfy numerous exceptional properties that simplify the problem. We find deep connections between Pauli Lie algebras and certain subgroups of the Clifford group generated by transvections, through the symplectic properties of the Pauli strings. This allows us to give an invariant-based perspective on these objects and their reachability, in the language of Pauli orbits, symmetries, and invariant subspaces. The invariant-based approach provides efficient algorithms for identifying Lie algebras and orbits, as well as a simple framework for analyzing structured Pauli generating sets. We also show in an elementary way that Clifford subgroups generated by transvections provide 3-designs for the corresponding Pauli Lie groups. We illustrate the framework through structured examples from variational quantum algorithms, restricted quantum computation, many-body systems, and random circuits.