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

Disentangling interacting systems with fermionic Gaussian circuits: Application to quantum impurity models

Ang-Kun Wu, Benedikt Kloss, Wladislaw Krinitsin, Matthew Fishman, J. H. Pixley, Edwin Miles Stoudenmire

Published Jan 9, 2025• 11 citations
arXivHTMLPDF

Abstract

Tensor network quantum states are powerful tools for strongly correlated systems, tailored to capture local correlations such as in ground states with entanglement area laws. When applying tensor network states to interacting fermionic systems, a proper choice of the basis or orbitals can reduce the bond dimension of tensors and provide physically relevant orbitals. We introduce such a change of basis with unitary gates obtained from compressing fermionic Gaussian states into quantum circuits corresponding to various tensor networks. These circuits can reduce the ground-state entanglement entropy and improve the performance of algorithms such as the density matrix renormalization group. We study the Anderson impurity model with one and two impurities to show the potential of the method for improving computational efficiency and interpreting impurity physics. Furthermore, fermionic Gaussian circuits can also suppress entanglement during the time evolution out of low-energy state. Last, we consider Gaussian multiscale entanglement renormalization ansatz (GMERA) circuits which compress fermionic Gaussian states hierarchically. The emergent coarse-grained physical models from these GMERA circuits are studied in terms of their entanglement properties and suitability for performing time evolution.

DOI
10.1103/physrevb.111.035119

Related papers

Entanglement renormalization circuits for $2d$ Gaussian Fermion States

Sing Lam Wong, Andrew C. Potter

The simulation of entangled ground-states of quantum materials remains challenging for classical computational methods in more than one spatial dimension, and is a prime target for quantum computational advantage. To this end, an important goal is to identify efficient quantum state preparation protocols that minimize the physical qubit number and circuit depth resources required to capture higher-dimensional quantum correlations. This work introduces a quantum circuit compression algorithm for Gaussian fermion states based on the multi-scale entanglement renormalization ansatz (MERA), which provides an exponential reduction in the circuit depth required to approximate highly-entangled ground-states relevant for quantum materials simulations. The algorithm, termed two-dimensional Gaussian MERA ($2d$ GMERA), extends MERA techniques to compress higher-dimensional Gaussian states. Through numerical simulations of the Haldane model on a honeycomb lattice, the method is shown to accurately capture area-law entangled states including topologically trivial insulators, Chern insulators, and critical Dirac semimetals. While Gaussian states alone are classically simulable, this approach establishes empirical upper bounds on quantum resources needed to prepare free fermion states that are adiabatically connected to correlated ground states, providing guidance for implementing these protocols on near-term quantum devices and offering a foundation for simulating more complex quantum materials. Finally, we develop a novel fermion-to-qubit encoding scheme, based on an expanding $2d$ topological order, that enables implementing fermionic rotations via qubit Pauli rotations with constant Pauli weight independent of system size.

arXiv:2506.04200•Jun 4, 2025•0 citations

Fermionic Orbital Optimization in Tensor Network States

Christian Krumnow, Libor Veis, Örs Legeza +1 more

Tensor network states and specifically matrix-product states have proven to be a powerful tool for simulating ground states of strongly correlated spin models. Recently, they have also been applied to interacting fermionic problems, specifically in the context of quantum chemistry. A new freedom arising in such nonlocal fermionic systems is the choice of orbitals, it being far from clear what choice of fermionic orbitals to make. In this Letter, we propose a way to overcome this challenge. We suggest a method intertwining the optimization over matrix product states with suitable fermionic Gaussian mode transformations. The described algorithm generalizes basis changes in the spirit of the Hartree-Fock method to matrix-product states, and provides a black box tool for basis optimization in tensor network methods.

arXiv:1504.00042•Physical Review Letters•Nov 18, 2016•84 citations

Clifford Circuits Augmented Grassmann Matrix Product States

Atis Yosprakob, Tu, Wei-Lin, Tsuyoshi Okubo +2 more

Recent progress in combining Clifford circuits with tensor-network (TN) methods has shown that local Clifford disentanglers can reduce bipartite entanglement across TN bonds prior to tensor compression, thereby improving the efficiency of TN simulations. In this work, we embed local Clifford disentanglers in the Grassmann-tensor language to define a Clifford-augmented Grassmann matrix product state (CAGMPS) ansatz, and develop a density-matrix renormalization group (DMRG) framework based on this ansatz while preserving locality and fermion-parity structure. We benchmark the resulting CAGMPS--DMRG method on representative fermionic lattice systems, including the tight-binding, $t$-$V$, and $t$-$V$-$V'$ models. In all cases, Clifford augmentation systematically suppresses bipartite entanglement and improves the accuracy of the ground-state energy at a fixed bond dimension. We further show that the Grassmann-evenness condition, together with equivalence under entangling action, restricts the relevant two-site Clifford candidates to 12 inequivalent representatives, enabling a more economical disentangling search than approaches based on the standard two-qubit Clifford gate set. Our results suggest that the CAGMPS--DMRG method provides a scalable and efficient variational tool for strongly correlated fermionic systems.

arXiv:2510.04164•Oct 5, 2025•0 citations

Quantum impurity models using superpositions of fermionic Gaussian states: Practical methods and applications

Samuel Boutin, Bela Bauer

The coherent superposition of nonorthogonal fermionic Gaussian states has been shown to be an efficient approximation to the ground states of quantum impurity problems [Bravyi and Gosset, Commun. Math. Phys. 356, 451 (2017)]. We propose a practical approach for performing a variational calculation based on such states. Our method is based on approximate imaginary-time equations of motion that decouple the dynamics of each Gaussian state forming the ansatz. It is independent of the lattice connectivity of the model and the implementation is highly parallelizable. To benchmark our variational method, we calculate the spin-spin correlation function and R\'enyi entanglement entropy of an Anderson impurity, allowing us to identify the screening cloud and compare to density matrix renormalization group calculations. Secondly, we study the screening cloud of the two-channel Kondo model, a problem difficult to tackle using existing numerical tools.

arXiv:2105.01088•Physical Review Research•Aug 24, 2021•9 citations

Combining Grassmann algebra with entanglement renormalization method

Jie Lou, Yan Chen

By combining the Grassmann algebra with multi-scale entanglement renormalization ansatz (MERA), we introduce a new unbiased and effective numerical method for simulating 2D strongly correlated electronic systems. The new GMERA method inherits all the advantages of MERA, which constructs the variational wave function based on complicated tensor network. Besides it can deal with fermionic properties of the system due to Grassmann algebra through local tensor contractions. This general method can treat different tensor network structures in a universal way. We show several benchmark calculations of the GMERA method, including the free fermion model, tight binding model, as well as the t-J model with hole doping.

arXiv:1506.03716•Jun 11, 2015•0 citations

Variational Power of Quantum Circuit Tensor Networks

Reza Haghshenas, Johnnie Gray, Andrew C. Potter +1 more

We characterize the variational power of quantum circuit tensor networks in the representation of physical many-body ground states.Such tensor networks are formed by replacing the dense block unitaries and isometries in standard tensor networks by local quantum circuits.We explore both quantum circuit matrix product states and the quantum circuit multiscale entanglement renormalization Ansatz, and introduce an adaptive method to optimize the resulting circuits to high fidelity with more than 10 4 parameters.We benchmark their expressiveness against standard tensor networks, as well as other common circuit architectures, for the 1D and 2D Heisenberg and 1D Fermi-Hubbard models.We find quantum circuit tensor networks to be substantially more expressive than other quantum circuits for these problems, and that they can even be more compact than standard tensor networks.Extrapolating to circuit depths which can no longer be emulated classically, this suggests a region of advantage in quantum expressiveness in the representation of physical ground states.

arXiv:2107.01307•Physical Review X•Mar 11, 2022•80 citations
Tensor network quantum states are powerful tools for strongly correlated systems, tailored to capture local correlations such as in ground states with entanglement area laws. When applying tensor network states to interacting fermionic systems, a proper choice of the basis or orbitals can reduce the bond dimension of tensors and provide physically relevant orbitals. We introduce such a change of basis with unitary gates obtained from compressing fermionic Gaussian states into quantum circuits corresponding to various tensor networks. These circuits can reduce the ground-state entanglement entropy and improve the performance of algorithms such as the density matrix renormalization group. We study the Anderson impurity model with one and two impurities to show the potential of the method for improving computational efficiency and interpreting impurity physics. Furthermore, fermionic Gaussian circuits can also suppress entanglement during the time evolution out of low-energy state. Last, we consider Gaussian multiscale entanglement renormalization ansatz (GMERA) circuits which compress fermionic Gaussian states hierarchically. The emergent coarse-grained physical models from these GMERA circuits are studied in terms of their entanglement properties and suitability for performing time evolution.
Entanglement renormalization circuits for 2d2d2d Gaussian Fermion States
The simulation of entangled ground-states of quantum materials remains challenging for classical computational methods in more than one spatial dimension, and is a prime target for quantum computational advantage. To this end, an important goal is to identify efficient quantum state preparation protocols that minimize the physical qubit number and circuit depth resources required to capture higher-dimensional quantum correlations. This work introduces a quantum circuit compression algorithm for Gaussian fermion states based on the multi-scale entanglement renormalization ansatz (MERA), which provides an exponential reduction in the circuit depth required to approximate highly-entangled ground-states relevant for quantum materials simulations. The algorithm, termed two-dimensional Gaussian MERA (2d2d2d GMERA), extends MERA techniques to compress higher-dimensional Gaussian states. Through numerical simulations of the Haldane model on a honeycomb lattice, the method is shown to accurately capture area-law entangled states including topologically trivial insulators, Chern insulators, and critical Dirac semimetals. While Gaussian states alone are classically simulable, this approach establishes empirical upper bounds on quantum resources needed to prepare free fermion states that are adiabatically connected to correlated ground states, providing guidance for implementing these protocols on near-term quantum devices and offering a foundation for simulating more complex quantum materials. Finally, we develop a novel fermion-to-qubit encoding scheme, based on an expanding 2d2d2d topological order, that enables implementing fermionic rotations via qubit Pauli rotations with constant Pauli weight independent of system size.
Tensor network states and specifically matrix-product states have proven to be a powerful tool for simulating ground states of strongly correlated spin models. Recently, they have also been applied to interacting fermionic problems, specifically in the context of quantum chemistry. A new freedom arising in such nonlocal fermionic systems is the choice of orbitals, it being far from clear what choice of fermionic orbitals to make. In this Letter, we propose a way to overcome this challenge. We suggest a method intertwining the optimization over matrix product states with suitable fermionic Gaussian mode transformations. The described algorithm generalizes basis changes in the spirit of the Hartree-Fock method to matrix-product states, and provides a black box tool for basis optimization in tensor network methods.
Recent progress in combining Clifford circuits with tensor-network (TN) methods has shown that local Clifford disentanglers can reduce bipartite entanglement across TN bonds prior to tensor compression, thereby improving the efficiency of TN simulations. In this work, we embed local Clifford disentanglers in the Grassmann-tensor language to define a Clifford-augmented Grassmann matrix product state (CAGMPS) ansatz, and develop a density-matrix renormalization group (DMRG) framework based on this ansatz while preserving locality and fermion-parity structure. We benchmark the resulting CAGMPS--DMRG method on representative fermionic lattice systems, including the tight-binding, ttt-VVV, and ttt-VVV-V′V'V′ models. In all cases, Clifford augmentation systematically suppresses bipartite entanglement and improves the accuracy of the ground-state energy at a fixed bond dimension. We further show that the Grassmann-evenness condition, together with equivalence under entangling action, restricts the relevant two-site Clifford candidates to 12 inequivalent representatives, enabling a more economical disentangling search than approaches based on the standard two-qubit Clifford gate set. Our results suggest that the CAGMPS--DMRG method provides a scalable and efficient variational tool for strongly correlated fermionic systems.
The coherent superposition of nonorthogonal fermionic Gaussian states has been shown to be an efficient approximation to the ground states of quantum impurity problems [Bravyi and Gosset, Commun. Math. Phys. 356, 451 (2017)]. We propose a practical approach for performing a variational calculation based on such states. Our method is based on approximate imaginary-time equations of motion that decouple the dynamics of each Gaussian state forming the ansatz. It is independent of the lattice connectivity of the model and the implementation is highly parallelizable. To benchmark our variational method, we calculate the spin-spin correlation function and R\'enyi entanglement entropy of an Anderson impurity, allowing us to identify the screening cloud and compare to density matrix renormalization group calculations. Secondly, we study the screening cloud of the two-channel Kondo model, a problem difficult to tackle using existing numerical tools.
By combining the Grassmann algebra with multi-scale entanglement renormalization ansatz (MERA), we introduce a new unbiased and effective numerical method for simulating 2D strongly correlated electronic systems. The new GMERA method inherits all the advantages of MERA, which constructs the variational wave function based on complicated tensor network. Besides it can deal with fermionic properties of the system due to Grassmann algebra through local tensor contractions. This general method can treat different tensor network structures in a universal way. We show several benchmark calculations of the GMERA method, including the free fermion model, tight binding model, as well as the t-J model with hole doping.
We characterize the variational power of quantum circuit tensor networks in the representation of physical many-body ground states.Such tensor networks are formed by replacing the dense block unitaries and isometries in standard tensor networks by local quantum circuits.We explore both quantum circuit matrix product states and the quantum circuit multiscale entanglement renormalization Ansatz, and introduce an adaptive method to optimize the resulting circuits to high fidelity with more than 10 4 parameters.We benchmark their expressiveness against standard tensor networks, as well as other common circuit architectures, for the 1D and 2D Heisenberg and 1D Fermi-Hubbard models.We find quantum circuit tensor networks to be substantially more expressive than other quantum circuits for these problems, and that they can even be more compact than standard tensor networks.Extrapolating to circuit depths which can no longer be emulated classically, this suggests a region of advantage in quantum expressiveness in the representation of physical ground states.