We study brick wall quantum circuits enjoying a global fermionic symmetry. The constituent 2-qubit gate, and its fermionic symmetry, derive from a 2-particle scattering matrix in integrable, supersymmetric quantum field theory in 1+1 dimensions. Our 2-qubit gate, as a function of three free parameters, is of so-called free fermionic or matchgate form, allowing us to derive the spectral structure of both the brick wall unitary U_F U F and its, non-trivial, Hamiltonian limit H_{\gamma} H γ in closed form. We find that the fermionic symmetry pins H_{\gamma} H γ to a surface of critical points, whereas breaking that symmetry leads to non-trivial topological phases. We briefly explore quench dynamics for this class of circuits.
We present a classical simulation method for fermionic quantum systems which, without loss of generality, can be represented by parity-preserving circuits made of two-qubit gates in a brick-wall structure. We map such circuits to a fermionic tensor network and introduce a novel decomposition of non-Matchgate gates into a Gaussian fermionic tensor and a residual quartic term, inspired by interacting fermionic systems. The quartic term is independent of the specific gate, which allows us to precompute intermediate results independently of the exact circuit structure and leads to significant speedups when compared to other methods. Our decomposition suggests a natural perturbative expansion which can be turned into an algorithm to compute measurement outcomes and observables to finite accuracy when truncating at some order of the expansion. For particle number conserving gates, our decomposition features a unique truncation cutoff reducing the computational effort for high precision calculations. Our algorithm significantly lowers resource requirements for simulating parity-preserving circuits while retaining high accuracy, making it suitable for simulations of interacting systems in quantum chemistry and material science. Lastly, we discuss how our algorithm compares to other classical simulation methods for fermionic quantum systems.
We present quantum circuits with a brick wall structure using the optimal number of parameters and two-qubit gates to parametrize
SU(2n), and provide evidence that these circuits are universal for
n≤5. For this, we successfully compile random matrices to the presented circuits and show that their Jacobian has full rank almost everywhere in the domain. Our method provides a new state of the art for synthesizing typical unitary matrices from
SU(2n) for
n=3,4,5, and we extend it to the subgroups
SO(2n) and
Sp∗(2n). We complement this numerical method by a partial proof, which hinges on an open conjecture that relates universality of an ansatz to it having full Jacobian rank almost everywhere.
Matchgate circuits realize free-fermion dynamics: they are efficiently classically simulable, yet cannot on their own generate the generic randomness required for universal computation or unitary design formation. We study a controlled route beyond this integrable limit by doping matchgate circuits with non-Gaussian gates-physically, the injection of fermionic interactions into an otherwise free system. Using the matchgate commutant framework, we obtain analytic control over unitary
2-design formation. For globally scrambled dynamics, the design problem maps exactly onto a classical birth-death Markov chain with an Ornstein-Uhlenbeck continuum limit, recasting the emergence of quantum randomness in terms of spectral gaps and mixing times and yielding rigorous bounds on the number of non-Gaussian gates needed for approximate
2-designs. These bounds hold for a broad class of parity-preserving non-Gaussian gates, independently of microscopic details, with numerics indicating that the same mechanism governs higher-order designs. Used as local building blocks in a glued-circuit architecture, they yield approximate parity-preserving
2-designs in polylogarithmic depth with a sparse non-Gaussian gate count, with implications for Page-like entanglement growth and fermionic classical-shadow protocols. Finally, locality reshapes this picture: in local brickwork dynamics, design formation is diffusion-limited and far slower. Our results establish doped matchgate circuits as a controlled, analytically tractable route from free fermions to interaction-generated quantum designs.
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.
Random circuits giving rise to unitary designs are key tools in quantum information science and many-body physics. In this work, we investigate a class of random quantum circuits with a specific gate structure. Within this framework, we prove that one-dimensional structured random circuits with non-Haar-random local gates can exhibit substantially more global randomness compared to Haar-random circuits with the same underlying circuit architecture. In particular, under a solvability condition, we exactly diagonalize the second-moment operators for these structured random circuits by establishing a link to free-fermionic models, and we show that their spectral gaps can exceed those of Haar-random circuits. While solvable structured random circuits exhibit free-fermionic behavior, the underlying fixed gate need not be a matchgate. Our findings have applications in improving circuit depth bounds for randomized benchmarking, in generating approximate unitary 2-designs from shallow random circuits, and in enabling exact analyses of quantum many-body dynamics.
In this paper, we address the problem of Yang-Baxter integrability of doubled quantum circuit of qubits (spins 1/2) with open boundary conditions where the two circuit replicas are only coupled at the left or right boundary. We investigate the cases where the bulk is given by elementary six vertex unitary gates of either the free fermionic XX type or interacting XXZ type. By using the Sklyanin’s construction of reflection algebra, we obtain the most general solutions of the boundary Yang-Baxter equation for such a setup. We use this solution to build, from the transfer matrix formalism, integrable circuits with two step discrete time Floquet (aka brickwork) dynamics. We prove that, only if the bulk is a free-model, the boundary matrices are in general non-factorizable, and for particular choice of free parameters yield non-trivial unitary dynamics with boundary interaction between the two chains. Then, we consider the limit of continuous time evolution and we give the interpretation of a restricted set of the boundary terms in the Lindbladian setting. Specifically, for a particular choice of free parameters, the solutions correspond to an open quantum system dynamics with the source terms representing injecting or removing particles from the boundary of the spin chain.