---
title: Hidden Fermion Pfaffian State (HFPS)
url: https://www.emergentmind.com/topics/hidden-fermion-pfaffian-state-hfps
type: topic
---

# Hidden Fermion Pfaffian State (HFPS)

Searching arXiv for recent and foundational papers on Hidden Fermion Pfaffian State and related Pfaffian literature.
First, locating the dedicated HFPS paper and closely related Pfaffian references.
Hidden Fermion Pfaffian State (HFPS) denotes a neural quantum state architecture in which hidden fermions and neural-network backflow ideas are extended from Slater determinants to Pfaffians, thereby yielding a variational class that can represent both unpaired and superconducting fermionic phases within scalable variational Monte Carlo [2507.10705]. In a broader but less formal sense, the phrase also resonates with several strands of Pfaffian literature in which paired composite-fermion states, Pfaffian daughter hierarchies, particle-hole-symmetric paired states, and geometric response of Pfaffian phases are central, although those works do not explicitly use the term HFPS [2309.00111], [1704.05107], [1902.09563]. The narrow, formal meaning is therefore architectural and computational; the broader significance lies in its relation to Pfaffian pairing as a recurring organizing principle in strongly correlated many-body physics.

## 1. Formal definition and variational setting

HFPS is introduced as a response to a specific limitation of earlier fermionic neural quantum states. Hidden fermion determinant state (HFDS) and neural network backflow (NNBF) are built around Slater determinants, so they are naturally suited to unpaired fermions and Fermi-liquid-like states. The motivating observation of HFPS is that strongly correlated fermion systems frequently involve pairing, including superconductivity, and that the natural number-projected mean-field language for pairing is the Pfaffian, or equivalently Thouless or BCS form, rather than a determinant [2507.10705].

The starting point is the standard bilinear mean-field Hamiltonian
$$
\hat{\mathcal H}_0=\sum_{p,q}^{M}\left(t_{pq}\hat c_p^\dagger \hat c_q+\frac{1}{2}\Delta_{pq}\hat c_p^\dagger \hat c_q^\dagger+\frac{1}{2}\Delta_{qp}^*\hat c_p\hat c_q\right),
$$
whose ground state is a generalized Gaussian state,
$$
\ket{\psi_0}=\exp\left(\frac{1}{2}\sum_{p,q}^{M}F_{pq}\hat c_p^\dagger \hat c_q^\dagger\right)\ket{0},
$$
with antisymmetric pairing matrix $\mathbf F$. After projection to fixed particle number $N$, the state becomes
$$
\ket{\psi_{\mathrm{pf}}}=\hat{\mathcal P}_N\ket{\psi_0}
=\frac{1}{(N/2)!}\left(\sum_{p<q}F_{pq}\hat c_p^\dagger \hat c_q^\dagger\right)^{N/2}\ket{0},
$$
and for a Fock configuration $\mathbf n$,
$$
\psi_{\mathrm{pf}}(\mathbf n)=\mathrm{pf}(\mathbf n\star \mathbf F\star \mathbf n).
$$
This is the basic Pfaffian wave-function that HFPS generalizes [2507.10705].

The formal role of HFPS is therefore not merely to add neural corrections to an existing paired state, but to replace the determinant-centric hidden-fermion architecture by a Pfaffian-centric one. That substitution is substantive because projected BCS states are a special case of Pfaffians, and determinants are contained inside the Pfaffian class rather than the reverse. A plausible implication is that HFPS is aimed at regimes where variational bias toward unpaired reference states would otherwise be problematic.

## 2. Hidden-sector construction and neural parametrization

HFPS enlarges the physical Hilbert space by adding hidden orbitals and hidden fermions. If the visible sector has $M$ orbitals and $N$ fermions, the hidden sector has $\tilde M$ orbitals and $\tilde N$ hidden fermions. In the enlarged space the wave-function is written as
$$
\ket{\psi}=
\frac{1}{(N_{\mathrm{tot}}/2)!}
\Big(
\sum_{p<q}F^{vv}_{pq}\hat c_p^\dagger \hat c_q^\dagger
+\sum_{\tilde p<\tilde q}F^{hh}_{\tilde p\tilde q}\hat d_{\tilde p}^\dagger \hat d_{\tilde q}^\dagger
+\sum_{p,\tilde p}F^{vh}_{p\tilde p}\hat c_p^\dagger \hat d_{\tilde p}^\dagger
\Big)^{N_{\mathrm{tot}}/2}\ket{0},
$$
with $N_{\mathrm{tot}}=N+\tilde N$ [2507.10705].

For visible configuration $\mathbf n$ and hidden configuration $\tilde{\mathbf n}$,
$$
\psi(\mathbf n,\tilde{\mathbf n})=
\mathrm{pf}
\begin{pmatrix}
\mathbf n\star \mathbf F^{vv}\star \mathbf n &
\mathbf n\star \mathbf F^{vh}\star \tilde{\mathbf n} \\
-\tilde{\mathbf n}\star (\mathbf F^{vh})^T\star \mathbf n &
\tilde{\mathbf n}\star \mathbf F^{hh}\star \tilde{\mathbf n}
\end{pmatrix}.
$$
To obtain a physical wave-function, the hidden occupations are tied to the visible configuration, yielding the configuration-dependent form
$$
\psi(\mathbf n)=
\mathrm{pf}
\begin{pmatrix}
\mathbf n\star \mathbf F^{vv}\star \mathbf n &
\mathbf n\star \mathbf F^{vh}(\mathbf n) \\
-\mathbf F^{vh}(\mathbf n)^T\star \mathbf n &
\mathbf F^{hh}(\mathbf n)
\end{pmatrix}.
$$

A practical issue is that $\mathbf F^{hh}(\mathbf n)$ is not directly compatible with standard network outputs. The construction therefore factors it through a unitary matrix $\mathbf V(\mathbf n)$,
$$
\mathbf F^{hh}(\mathbf n)=\mathbf V^T(\mathbf n)\,\tilde{\mathbf F}^{hh}\,\mathbf V(\mathbf n),
$$
and defines
$$
\tilde{\mathbf F}^{vh}(\mathbf n)=\mathbf F^{vh}(\mathbf n)\mathbf V^\dagger(\mathbf n),
\qquad
\mathbf U(\mathbf n)=\mathrm{diag}(\mathbf 1,\mathbf V(\mathbf n)).
$$
Using $\mathrm{pf}(BAB^T)=\det(B)\mathrm{pf}(A)$, the state becomes
$$
\psi(\mathbf n)=J(\mathbf n)\times
\mathrm{pf}
\begin{pmatrix}
\mathbf n\star \mathbf F^{vv}\star \mathbf n &
\mathbf n\star \tilde{\mathbf F}^{vh}(\mathbf n) \\
-\tilde{\mathbf F}^{vh}(\mathbf n)^T\star \mathbf n &
\tilde{\mathbf F}^{hh}
\end{pmatrix},
$$
where
$$
J(\mathbf n)=\det(\mathbf U(\mathbf n))
$$
is an additional generalized Jastrow factor produced by the ANN. In implementation, $\mathbf F^{vv}$ and $\tilde{\mathbf F}^{hh}$ are direct variational parameters, while a CNN, GCNN, or other ANN outputs $\tilde{\mathbf F}^{vh}(\mathbf n)$ and $J(\mathbf n)$ [2507.10705].

This formulation is also recast as a Pfaffian backflow ansatz,
$$
\psi(\mathbf n)=\tilde J(\mathbf n)\,
\mathrm{pf}\!\left[\mathbf n\star(\mathbf F^{vv}+\mathbf F'(\mathbf n))\star \mathbf n\right],
$$
with
$$
\tilde J(\mathbf n)=J(\mathbf n)\,\mathrm{pf}(\tilde{\mathbf F}^{hh}),
\qquad
\mathbf F'(\mathbf n)=
\tilde{\mathbf F}^{vh}(\mathbf n)(\tilde{\mathbf F}^{hh})^{-1}\tilde{\mathbf F}^{vh}(\mathbf n)^T.
$$
The backflow correction has rank bounded by the number of hidden fermions $\tilde N$, which is identified as a key advantage of the formulation [2507.10705].

## 3. Expressive hierarchy and relation to determinant-based ansätze

HFPS is designed to represent both unpaired states and paired superconducting states. The argument is structural. A Pfaffian can reduce to a Slater determinant, so determinant states are a subset of Pfaffian states. The paper states this explicitly through
$$
\psi_{\mathrm{det}}(\mathbf n)=\det(\mathbf n\star \mathbf u)
=
\mathrm{pf}(\mathbf n\star \mathbf u J \mathbf u^T\star \mathbf n),
$$
for any antisymmetric $\mathbf J$ with $\mathrm{pf}(\mathbf J)=1$ [2507.10705].

The relation to HFDS is equally explicit. The earlier hidden fermion determinant state is
$$
\psi_{\mathrm{HFDS}}(\mathbf n)=
\det
\begin{pmatrix}
\mathbf n\star \mathbf u^v \\
\mathbf u^h(\mathbf n)
\end{pmatrix}.
$$
HFPS can reproduce HFDS by choosing
$$
\mathbf F^{vv}=\mathbf u^v\mathbf J(\mathbf u^v)^T,\quad
\mathbf F^{vh}(\mathbf n)=\mathbf u^v\mathbf J(\mathbf u^h(\mathbf n))^T,\quad
\mathbf F^{hh}(\mathbf n)=\mathbf u^h(\mathbf n)\mathbf J(\mathbf u^h(\mathbf n))^T,
$$
but not vice versa. The formal conclusion drawn in the source is that HFPS is strictly more general because some Pfaffians cannot be represented as HFDS, whereas HFDS imposes a determinant structure on the enlarged space [2507.10705].

Relative to NNBF, HFPS is mapped to a Pfaffian backflow form with a controlled low-rank correction. The emphasis on controlled rank is not incidental: it underwrites both expressivity and computational tractability. This suggests that the architectural novelty of HFPS is not only the replacement of determinants by Pfaffians, but also the replacement of generic backflow by a rank-bounded correction tied directly to hidden-fermion number.

## 4. Computational scaling and large-system evaluation

The computational claim of HFPS is that it remains compatible with large-scale VMC despite the use of Pfaffians. A naive Pfaffian evaluation costs $\mathcal O(n^3)$ for an $n\times n$ antisymmetric matrix. If the full enlarged Pfaffian were recomputed directly at every Monte Carlo step, the cost would be substantial. The paper identifies two accelerations that change this conclusion [2507.10705].

The first is low-rank updates. Since VMC changes only a few occupations at each Monte Carlo move, the sliced Pfaffian matrix changes by a low-rank modification. With stored inverse information, a Pfaffian update is reduced from
$$
\mathcal O((N+\tilde N)^3)\to \mathcal O(N^2\tilde N).
$$
The second is sublattice or translation symmetry. By choosing a unit-cell structure in $\mathbf F^{vv}$, translations that differ only by a sublattice permutation do not require recomputing the Pfaffian from scratch. The effect is to reduce the cost of symmetry projection from scaling with the full system size to scaling with the unit cell [2507.10705].

The summary scaling given in the source is that, without acceleration, a full VMC step scales like
$$
\mathcal O(NM(N+\tilde N)^3),
$$
which for fixed density becomes roughly $\mathcal O(N^5)$. With low-rank updates and sublattice symmetry, the HFPS forward pass scales as $\mathcal O(N^2)$ and a full VMC step as $\mathcal O(N^3)$ [2507.10705]. In the terminology of the abstract, HFPS thus scales to large systems with favorable asymptotic complexity.

## 5. Numerical performance in Hubbard-model benchmarks

The benchmark model is the square-lattice Hubbard Hamiltonian
$$
\hat{\mathcal H}=
-t\sum_{\langle ij\rangle,\sigma}\hat c_{i,\sigma}^\dagger \hat c_{j,\sigma}
+U\sum_i \hat n_{i,\uparrow}\hat n_{i,\downarrow},
$$
and performance is measured using the relative variational energy error
$$
\epsilon_{\mathrm{rel}}=\frac{E-E_0}{E_\infty-E_0},
$$
where $E_0$ is the exact or benchmark ground-state energy and $E_\infty$ is the infinite-temperature reference [2507.10705].

In the weakly correlated Fermi-liquid regime on a $4\times4$ lattice with PBC and $N=10$ at filling $n=5/8$, HFPS is compared to determinant+RBM and HFDS. The reported result is that HFPS improves the variational energy significantly over HFDS, reduces the relative error by roughly two orders of magnitude compared with the best previous NQS result, and yields raw energies extremely close to exact diagonalization [2507.10705].

For the attractive Hubbard model on an $8\times8$ lattice at $n=7/8$, where the system exhibits s-wave superconductivity interpolating from BCS at weak coupling to BEC-like local pairs at stronger coupling, the diagnostic is the pair-pair correlator
$$
C_p(\mathbf r)=\langle \hat\Delta^\dagger(\mathbf r)\hat\Delta(\mathbf 0)\rangle,
$$
with
$$
|\langle \hat\Delta\rangle|
\sim
\left|\lim_{|\mathbf r|\to\infty}C_p(\mathbf r)\right|^{1/2}
$$
used as the ODLRO diagnostic. HFPS, initialized from a BCS mean-field state, produces pair correlations that match DQMC within error bars and reaches energy errors of order
$$
\epsilon_{\mathrm{rel}}\sim 10^{-4},
$$
while the BCS wave-function remains around $10^{-2}$. This is presented as a clear demonstration that HFPS accurately captures s-wave superconductivity across the BCS–BEC crossover [2507.10705].

At half-filling on the repulsive $8\times8$ Hubbard model, the benchmarks are against AFQMC, with comparisons to PP+RBM + Lanczos, HFDS, fPEPS ($D=16$), and Transformer-NNBF. The reported conclusion is that HFPS is more accurate than these methods, with relative energy error at the $10^{-4}$ to $10^{-3}$ scale and often better than the competing NQSs. The same source notes that benchmark uncertainty from AFQMC becomes non-negligible when the variational error becomes this small [2507.10705].

The most demanding test is the stripe phase on $L\times4$ systems at $U=8$ and $n=7/8$. HFPS is reported to outperform MLP-NNBF, HFDS, TBF, TBF+Lanczos, and Transformer-NNBF; to achieve a nearly exact energy on $4\times4$ with
$$
\epsilon_{\mathrm{rel}}\approx 1.2\times 10^{-5};
$$
and to reach
$$
E/N=-0.76413
$$
on $16\times4$, better than the reported Transformer-NNBF result
$$
E/N=-0.76298.
$$
In that regime the diagnostics show spin density wave order with period 16, charge density wave order with period 8, and finite short-range d-wave pairing correlations. The paper does not claim a d-wave superconductor in the $1/8$-doped stripe phase; rather, it argues that d-wave pairing is an important short-range ingredient of the energetics and that HFPS represents it accurately [2507.10705].

The d-wave pairing operator used in that analysis is
$$
\hat\Delta(\mathbf r)=
\frac{1}{4}\sum_{\boldsymbol\delta}
\frac{1}{\sqrt 2}\,\mathrm{sign}(\boldsymbol\delta)
\left(
\hat c_{\mathbf r,\uparrow}\hat c_{\mathbf r+\boldsymbol\delta,\downarrow}
-
\hat c_{\mathbf r,\downarrow}\hat c_{\mathbf r+\boldsymbol\delta,\uparrow}
\right),
$$
with $\mathrm{sign}(\boldsymbol\delta)=+1$ for horizontal bonds and $-1$ for vertical bonds. The mean-field Hamiltonian used to seed that structure is
$$
\hat{\mathcal H}'=
\hat{\mathcal H}
+
\Delta\sum_{\langle ij\rangle}\mathrm{sign}(i,j)
\left(
\hat c^\dagger_{i,\uparrow}\hat c^\dagger_{j,\downarrow}
-
\hat c^\dagger_{i,\downarrow}\hat c^\dagger_{j,\uparrow}
+\mathrm{H.c.}
\right).
$$
Taken together, these experiments support the claim that HFPS captures both s-wave and d-wave pairing channels and provides state-of-the-art variational accuracy in different regimes of both the attractive and repulsive Hubbard models [2507.10705].

## 6. Relation to Pfaffian topological phases and broader usage

The formal HFPS architecture should be distinguished from the broader Pfaffian literature in quantum Hall and synthetic platforms, where the phrase “Hidden Fermion Pfaffian State” is usually absent even when the phenomenology is close. The following works are central points of contact.

| Paper | System | Relevance to HFPS |
|---|---|---|
| [2309.00111] | Wide GaAs quantum wells near $\nu=1/2$ | One-component Pfaffian parent and Pfaffian daughter fractions |
| [1902.09563] | Pfaffian state on curved Riemann surfaces | Universal geometric and anomaly response |
| [1704.05107] | ZnO with strong Landau level mixing | PH-symmetric Pfaffian-like paired state |
| [2207.00834] | AB-stacked bilayer graphene | Tunable half-filled Pfaffian regime |
| [1603.02348] | Coupled atom-cavity arrays | Bosonic Pfaffian-like analogue with effective three-body interaction |

In wide GaAs quantum wells, the observed $\nu=1/2$ FQHS is argued to be a one-component Pfaffian rather than a competing two-component Jain/331 state. The core evidence is the abrupt emergence of unusually strong $\nu=8/17$ and $\nu=7/13$ states at the same tuning point where $\nu=1/2$ becomes especially robust. Since $8/17$ and $7/13$ are identified as precisely the theoretically predicted, simplest daughter states of a one-component Pfaffian $1/2$ FQHS, the data are interpreted as showing a topological phase transition from ordinary Jain-sequence fractions to Pfaffian hierarchical daughter states [2309.00111]. That work does not use the term HFPS, but it is explicitly described as closely aligned with the idea of a robust $\nu=1/2$ parent whose side fractions behave like Pfaffian daughters rather than standard Jain states.

In the geometric-response literature, the Pfaffian is formulated on arbitrary compact Riemann surfaces with curved metric and inhomogeneous magnetic field through an Ising-CFT correlator composed of a free compact boson and a free Majorana fermion. The resulting generating functional encodes linear response to metric and magnetic-field variations, while the effective action contains Aubin–Yau, Mabuchi, and Liouville terms. The density expansion,
$$
\langle \rho(\mathbf z)\rangle_{\rm Pfaffian}
=
\frac{\nu}{2\pi}N_\phi
+\frac{\sqrt\nu\,Q}{8\pi}R
+\frac{1}{8\pi N_\phi}
\left[
\frac18+\frac{\sqrt\nu\,Q}{4}\left(1-\frac1\nu\right)
\right]\Delta_g R
+\mathcal O(N_\phi^{-2}),
$$
and the long-wavelength static structure factor,
$$
S_{\rm Pfaffian}(k)=
\frac{k^2}{2}
+\frac{1-\nu}{8\nu}k^4
+\frac{(1-2\nu)(2-\nu)}{64\nu^2}k^6
+\mathcal O(k^8),
$$
exhibit the gravitational-anomaly contribution associated with central charge $c=3/2$ [1902.09563]. If HFPS is interpreted as a physical Pfaffian state rather than only a variational architecture, this suggests inheritance of the Pfaffian’s shift or spin-response structure, Hall-viscosity-like curvature coupling, and $c=3/2$ anomaly signature.

In ZnO, where Landau level mixing is non-perturbatively strong, the screened Coulomb interaction approach yields evidence for a particle-hole-symmetric Pfaffian-like incompressible state. The overlap with the Pfaffian trial state is strongly system-size dependent and never very high, reaching at most about $0.65$ in the calculations summarized, while the range $2<\kappa<3$ is identified as a possible topological phase-transition regime [1704.05107]. This work also does not use HFPS terminology, but it is described as conceptually close because it moves beyond a simple Pfaffian versus anti-Pfaffian dichotomy toward a PH-symmetric paired phase.

In AB-stacked bilayer graphene, the Pfaffian problem is recast in a tunable Dirac-fermion setting. The special bilayer Landau levels $0^{(+)}_{-1}$ and $0^{(-)}_1$ can host a $\nu=1/2$ Pfaffian state, and the overlap with the Pfaffian trial state together with the collective gap shows a nonmonotonic dependence on magnetic field, with a representative maximum near $B\approx 10$ T for $\gamma_1=400$ meV and $\gamma_1/\epsilon_B\approx 4.9$ [2207.00834]. The same source emphasizes quasiparticles of fractional charge $e/4$ and non-Abelian statistics. Here again, the paper does not label the phase HFPS, but it does describe a paired composite-fermion Pfaffian state of the type to which the phrase may be applied informally.

Finally, a synthetic analogue appears in coupled atom-cavity arrays. By replacing the two-level atom in the Jaynes–Cummings–Hubbard model with a three-level atom and tuning
$$
\epsilon_1=\omega-\Delta,\qquad
\epsilon_2=\omega-2\Delta,\qquad
g_2/g_1=\sqrt{2},
$$
the effective two-body interaction can be suppressed while a three-body repulsion remains comparatively dominant. On a torus, the reported signatures are a quasi-gap, a three-dimensional quasi-degenerate ground-state manifold on a $4\times4$ lattice with 4 particles and 4 flux quanta, a total Chern number of 3, and significant overlap with the Pfaffian trial wavefunction [1603.02348]. The authors describe the state as “Pfaffian-like,” not as a rigorous proof of the exact Moore–Read phase. The connection to HFPS is therefore analogical rather than terminological.

A recurrent misconception is that HFPS is simply another name for the Moore–Read Pfaffian. The literature provided here does not support that equivalence. The explicit formal term “Hidden Fermion Pfaffian State” names a neural-network-augmented Pfaffian variational ansatz for interacting lattice fermions [2507.10705]. By contrast, the fractional quantum Hall, graphene, ZnO, curved-space, and cavity-QED papers discuss physical Pfaffian or Pfaffian-like phases, and several state explicitly that they do not use the term HFPS [2309.00111], [1704.05107], [2207.00834], [1603.02348]. The more precise statement is that HFPS, in its formal sense, is a computational architecture whose conceptual center is Pfaffian pairing, while several neighboring literatures supply the topological, geometric, and experimental contexts in which Pfaffian structure acquires physical meaning.

Source: https://www.emergentmind.com/topics/hidden-fermion-pfaffian-state-hfps