Four-Layer Approximation Scheme on Torus

• The paper introduces a four-layer projective construction for the bosonic Z4 Read-Rezayi state, recovering the full ground state and excitation spectrum on a torus.
• It employs twisted boundary conditions and layer symmetrization to combine four Laughlin states into a single effective state capturing non-Abelian features.
• Numerical validation indicates high variational accuracy with rapid convergence of energy gaps and overlaps, ensuring practical feasibility in FQH studies.

The four-layer approximation scheme on the torus constitutes a projective construction specialized for the $\mathbb{Z}_4$ Read-Rezayi series within the fractional quantum Hall (FQH) paradigm. This protocol, introduced and elucidated in Repellin et al. (Repellin et al., 2015), synthesizes multilayer Laughlin states by symmetrizing over layers, augmented with twisted boundary conditions. This enables recovery of the full ground state and excitation spectrum that standard multilayer constructions fail to capture on toroidal geometries. The construction is expressible in both second-quantized lowest Landau level (LLL) orbital terms and first-quantized conformal-block (theta function) language, and is tailored for rigorous numerical and theoretical investigations. The scheme is most directly pertinent for the bosonic $k=4$ Read-Rezayi state at filling $\nu=2$.

1. Pre-Symmetrized Four-Layer Wave Function Construction

Commencing with four independent copies of the bosonic Laughlin $1/2$ state, each layer is defined on a torus of size $4N_\phi$ flux quanta, where $N_\phi = L_x L_y / (2\pi \ell_B^2)$ represents the fundamental flux threading per torus, and $\ell_B$ is the magnetic length. The LLL single-particle orbitals on the rectangular torus, in the $x$-Landau gauge $\mathbf{A}=(0,-Bx)$, are

$$\phi_j(z) = \frac{1}{\sqrt{L_y \sqrt{\pi}}} \sum_{m \in \mathbb{Z}} \exp\left[i y\left(\frac{2\pi j}{L_y} + m L_x\right) - \frac{(x + \frac{2\pi j}{L_y} + m L_x)^2}{2}\right],$$

with $k=4$0. In second quantization, the four-layer "parent" state is simply the tensor product:

$k=4$1

and in first quantization:

$k=4$2

where $k=4$3 and $k=4$4 is the odd Jacobi theta function. The conventional center-of-mass prefactor is omitted, as in the cited reference.

2. Twisted Boundary Conditions and Layer Coupling

The hallmark innovation of this scheme is the imposition of twisted boundary conditions, cyclically permuting layers as coordinates wind around the fundamental cycles of the torus:

• $k=4$5-cycle: $k=4$6
• $k=4$7-cycle: $k=4$8

In operator language this translates to:

$k=4$9

for $\nu=2$0 modulo $\nu=2$1. This twist identifies the four-layer torus as geometrically equivalent to a single torus of periods $\nu=2$2 (for $\nu=2$3 twists) or $\nu=2$4 (for $\nu=2$5 twists), effectively folding the layers into an extended toroidal geometry. This step is indispensable for recovering all $\nu=2$6 ground states on the torus.

3. Symmetrization Over the Layer Index

Projecting onto the $\nu=2$7 Read-Rezayi subspace is accomplished by the symmetrizer $\nu=2$8, defined by summing over permutations $\nu=2$9 that act on the layer indices:

$1/2$0

where $1/2$1 permutes layer labels on creation operators. In first quantization, the effect is:

$1/2$2

with normalization ensuring the resulting state is properly normalized. This symmetrization is critical for accessing the non-Abelian parafermion structure inherent in the $1/2$3 Read-Rezayi wave function.

4. Single-Layer $1/2$4 Trial State and Ground-State Multiplet

The symmetrized state realizes a single-layer wave function at filling $1/2$5 corresponding to the $1/2$6 Read-Rezayi phase. In first-quantized theta-function language:

$1/2$7

Expanded, this wave function matches the parafermion conformal block of the $1/2$8 series. For even $1/2$9, the ground-state manifold on the torus exhibits a fivefold degeneracy (as $4N_\phi$0 for $4N_\phi$1). These states are distinguished by translation quantum numbers $4N_\phi$2 in the reduced Brillouin zone; for $4N_\phi$3 they are: $4N_\phi$4, $4N_\phi$5, $4N_\phi$6, $4N_\phi$7, $4N_\phi$8. Generation of the full multiplet follows either from selecting distinct initial momenta for the parent Laughlin layers or by translation operators post-symmetrization.

5. Quasihole and Neutral Excitation Construction

Quasihole zero modes in the $4N_\phi$9 scheme are realized by increasing the parent system flux by four, producing one additional single-layer flux quantum after symmetrization. This produces the complete set of quasihole states (verified computationally for lower $N_\phi = L_x L_y / (2\pi \ell_B^2)$0 in (Repellin et al., 2015); for $N_\phi = L_x L_y / (2\pi \ell_B^2)$1 the same mechanism applies). Neutral excitation trial states are accessed by promoting one or more Laughlin parents to a collective neutral mode (e.g., magnetoroton), then symmetrizing. For instance, three layers can remain in their ground state while the fourth is excited prior to symmetrization, yielding a neutral excitation trial for $N_\phi = L_x L_y / (2\pi \ell_B^2)$2.

6. Variational Accuracy: Numerical Overlaps and Gaps

Repellin et al. (Repellin et al., 2015) provide numerical validation for $N_\phi = L_x L_y / (2\pi \ell_B^2)$3, exhibiting overlaps exceeding $N_\phi = L_x L_y / (2\pi \ell_B^2)$4 and energy gap estimates converging within $N_\phi = L_x L_y / (2\pi \ell_B^2)$5 for $N_\phi = L_x L_y / (2\pi \ell_B^2)$6. For $N_\phi = L_x L_y / (2\pi \ell_B^2)$7, although explicit overlap tables are not published, the methodology implies similarly high variational accuracy, with expected overlaps in the range $N_\phi = L_x L_y / (2\pi \ell_B^2)$8 for small $N_\phi = L_x L_y / (2\pi \ell_B^2)$9 and rapid gap convergence scaling with system size. Calculation of such metrics for $\ell_B$0 would entail full diagonalization of the $\ell_B$1-body contact Hamiltonian and inner product evaluation with symmetrized trial states, ensuring completeness of the zero-mode manifold and consistency with the neutral excitation spectrum.

7. Practical Implementation Workflow

Concrete computational steps for the four-layer scheme entail:

1. Construct the LLL orbital basis $\ell_B$2 for the extended torus $\ell_B$3 or $\ell_B$4.
2. Formulate four independent Laughlin many-body states at flux $\ell_B$5 via diagonalization of the $\ell_B$6-body contact Hamiltonian, storing each as Fock vectors.
3. Realize the twist/extended-torus equivalence by projecting onto every fourth orbital subset, yielding four Laughlin copies.
4. Implement the symmetrizer $\ell_B$7 by summing over the $\ell_B$8 layer-index permutations, or equivalently by orbital shifts.
5. Apply $\ell_B$9 to the product state and normalize; project back to the single-layer Fock basis of $x$0 orbitals.
6. Generate the fivefold ground-state multiplet by repeating steps for each initial momentum configuration or via translation operations.
7. For quasiholes, increment parent flux per q-hole and recompute; for neutral modes, substitute parent layers by their neutral-mode diagonalizations before symmetrization.
8. Finally, calculate overlaps and energies by taking inner products of symmetrized states with exact diagonalizations of the $x$1-body contact Hamiltonian at flux $x$2.

This sequence renders a comprehensive and information-rich application of the four-layer projective symmetrization on the torus (Repellin et al., 2015).