---
title: Gaussian Fermionic PEPS (GFPEPS)
url: https://www.emergentmind.com/topics/gaussian-fermionic-peps-gfpeps
type: topic
---

# Gaussian Fermionic PEPS (GFPEPS)

Gaussian Fermionic Projected Entangled Pair States (GFPEPS) are a class of tensor network states tailored for simulating quantum many-body systems with fermionic degrees of freedom. Defined by local tensors that are constrained to generate fermionic Gaussian states, GFPEPS provide an efficient, sign-problem-free representation for both free-fermion and certain interacting gapped and gapless phases in two and higher spatial dimensions. Their gauged extensions (GGFPEPS) further enable sign-free, variational Monte Carlo simulations of lattice gauge theories with dynamical fermions coupled to compact gauge groups, providing a numerically tractable ansatz for non-perturbative regimes of lattice gauge theories in two and three dimensions [2512.13812, 2404.13123, 2211.00023, 2412.01737].

## 1. Formal Definition and Covariance Structure

A GFPEPS is constructed on a lattice (typically d-dimensional), where each site carries physical fermionic modes $\psi_m(x)$ and 4 "legs" of virtual Majorana modes $a_{i\mu m}(x)$, with $i$ labeling the leg (e.g., corresponding to directions on the lattice), and $\mu, m$ labeling internal indices such as flavor or color. The local site operator is a fermionic Gaussian creator,
\[
A(x) = \exp\Big[ T_{\alpha\beta} \Psi^\dagger_\alpha(x) \Psi^\dagger_\beta(x) \Big],
\]
where $\Psi^\dagger_\alpha$ collects all physical and virtual creation operators [2512.13812].

Neighboring sites are entangled by applying Gaussian projectors $w(x,k)$ on each link, forming maximally entangled bond states between the virtual degrees of freedom [2008.11176, 2208.04623]. The global GFPEPS state is represented as
\[
|\Psi\rangle = \langle \Omega_v | \prod_{x,i} w^\dagger(x,i) \prod_x A(x) | \Omega_p \rangle.
\]
All physical properties of the GFPEPS, as a pure Gaussian state, are encoded in its antisymmetric covariance matrix,
\[
\Gamma_{ij} = \frac{i}{2} \langle [\gamma_i, \gamma_j] \rangle,
\]
where $\gamma_i$ are Majorana operators spanning all physical and virtual fermion spaces. For translationally invariant systems, covariance matrices reduce to momentum-space blocks, greatly facilitating analytic and numerical studies [2008.11176, 2304.06744].

In practical terms, the full PEPS tensor reduces, in the Gaussian case, to three real matrices per site, typically denoted $\{A,B,D\}$, encoding the local Gaussian map, with the output covariance constructed via a Schur complement:
\[
\Gamma_{\rm out} = A + B (D + \Gamma_{\rm in}^{-1})^{-1} B^T.
\]

## 2. Gauging Procedure and Extension to Lattice Gauge Theories

To construct gauged GFPEPS (GGFPEPS), a global symmetry group $G$ (such as $U(1)$, $\mathbb Z_N$, or $SU(N)$) is promoted to a local gauge symmetry. On each lattice link, a Hilbert space for the gauge field is introduced, and controlled-unitary operators $\mathcal{U}_G(x,i)$ act to entangle the virtual fermions with the physical gauge-field registers. The total gauged state on a 2D lattice is written as
\[
|\Psi\rangle = \langle \Omega_v | \prod_{x,i} w^\dagger(x,i) \mathcal{U}_G(x,i) \prod_x A(x) | \Omega_p \rangle \otimes |0\rangle,
\]
where $|0\rangle$ is the bare gauge-field vacuum [2512.13812, 2404.13123].

The physical content of the GGFPEPS is found by integrating over all gauge configurations, decomposing the state as
\[
|\Psi\rangle = \int \mathcal{D}\mathcal{G}\, | \mathcal{G} \rangle \otimes |\psi(\mathcal{G}) \rangle,
\]
where $|\psi(\mathcal{G}) \rangle$ is a pure fermionic Gaussian state for fixed gauge-field background $\mathcal{G}$. Expectation values of gauge-invariant observables reduce to Gaussian integrals over the covariance matrices, enabling the elimination of the sign problem inherent to conventional Monte Carlo algorithms for fermions [2211.00023, 2412.01737].

## 3. Variational Monte Carlo Framework and Algorithmic Aspects

The ground state of a gauge theory Hamiltonian is variationally optimized within the GGFPEPS manifold using Markov chain Monte Carlo over the gauge field configurations. The weight for each configuration $\mathcal{G}$ is the normalized squared norm of the (Gaussian) fermion wavefunction $|\psi(\mathcal{G})|^2$:
\[
p(\mathcal{G}) = \| \psi(\mathcal{G}) \|^2/\mathcal{Z},
\]
where $\mathcal{Z}$ is the normalization. Local updates are proposed and accepted with Metropolis probability $P_{\rm acc} = \min\{1, p(\mathcal{G}')/p(\mathcal{G}) \}$ [2512.13812].

Key computational bottlenecks—the evaluation of overlaps and observables—are polynomial in system size, as the Gaussian overlap (and hence the PEPS contraction) reduces to the computation of submatrices' Pfaffians or determinants. For local updates,
\[
p_0(\mathcal{G}) \propto \sqrt{\det \left[ (1 - \Gamma_{\rm in}(\mathcal{G}) D)/2 \right]},
\]
and Woodbury/matrix-determinant lemma techniques bring the update cost to $O(N^2)$ per step for local changes [2512.13812, 2211.00023, 2404.13123].

Algorithmic studies find that updating $\simeq$25–50% of the links at each Monte Carlo step yields optimal error convergence per wall clock time. Gauge fixing, while reducing the configuration space, can slow the convergence rate due to induced nonlocality in updates; partial “chessboard” gauge fixing avoids this problem, providing a favorable trade-off [2512.13812].

A summary of key algorithmic strategies is shown below:

| Algorithmic Parameter     | Optimal Choice/Effect                  | Reference               |
|--------------------------|----------------------------------------|-------------------------|
| Update size              | $k\sim$25–50% links per step           | [2512.13812]            |
| Gauge fixing             | Avoid maximal trees; chessboard/isolate | [2512.13812]            |
| Averaging (energy terms) | Magnetic: average all; Electric: one   | [2512.13812]            |
| Cost per MC step         | $O(N^2)$                               | [2512.13812,2211.00023] |

## 4. Physical Realizations and Applications

GFPEPS and their

Source: https://www.emergentmind.com/topics/gaussian-fermionic-peps-gfpeps