Noninteracting tight-binding models for Fock parafermions

This paper demonstrates that noninteracting $p$-state Fock parafermions on a one-dimensional lattice can be mapped to simpler $q_m$-state parafermions (or fermions when $p$ is a power of two), enabling the construction of solvable parafermionic Hamiltonians with single-particle spectra identical to their fermionic counterparts.

The Gist

Imagine you are trying to organize a very strict, very crowded party.

In the world of physics, there are two main types of guests: Bosons (the party animals who love to crowd into the same spot, like a mosh pit) and Fermions (the introverts who refuse to share a seat; if one person is sitting there, no one else can). This rule is called the "Exclusion Principle."

But what if you wanted a guest who is somewhere in between? A guest who is polite enough to let a few people in, but not too many? This is the world of Parafermions. They are the "Goldilocks" particles: not too crowded, not too empty. Specifically, this paper focuses on 4-state Parafermions, which means each "seat" (or orbital) can hold 0, 1, 2, or 3 guests, but never 4.

The Problem: The Party is Too Complicated

Usually, when physicists try to calculate how these "Goldilocks" particles behave in a chain (like a line of seats), the math gets incredibly messy. Because they have these weird "in-between" rules, you can't just use the standard, simple formulas that work for regular fermions. It's like trying to solve a Rubik's cube while wearing oven mitts.

The Solution: The "Fermion Translator"

The author of this paper, Edward McCann, discovered a clever trick. He realized that for these specific 4-state particles, you don't need to invent new math. Instead, you can translate them into a language we already understand perfectly: regular fermions.

Here is the analogy:
Imagine the 4-state Parafermion seat is a special VIP booth.

• 0 guests = The booth is empty.
• 1 guest = One person is sitting there.
• 2 guests = Two people are sitting there.
• 3 guests = Three people are sitting there.

McCann showed that you can model this VIP booth by using two separate, regular fermion booths right next to each other:

1. Booth A (Spin Up): Can hold 0 or 1 person.
2. Booth B (Spin Down): Can hold 0 or 1 person.

Now, here is the magic translation:

• 0 Parafermions = 0 in Booth A + 0 in Booth B.
• 1 Parafermion = 1 in Booth A + 0 in Booth B.
• 2 Parafermions = 0 in Booth A + 1 in Booth B. (Wait, why 0 and 1? Because the math treats the "Down" booth as having double the weight! So 1 person in the "Down" booth counts as 2 guests in the VIP booth).
• 3 Parafermions = 1 in Booth A + 1 in Booth B.

By splitting the complex 4-state problem into two simple 2-state (regular fermion) problems, the math suddenly becomes easy. It's like realizing that a complicated 4-digit combination lock is actually just two simple 2-digit locks stuck together.

Why This Matters: The "Thermometer" Test

The paper doesn't just stop at the math trick. It asks: "Does this translation actually work in the real world?"

To test this, the author looked at temperature. If you heat up a system of these particles, how does their energy change?

• For regular fermions, there is a standard formula (Fermi-Dirac statistics).
• For these 4-state particles, there is a middle-ground formula (Gentile statistics).

The author found that the "middle-ground" behavior is exactly what you get if you take your two regular fermion booths and treat one of them as if it were twice as cold as the other.

• The "Up" fermions act normally.
• The "Down" fermions act like they are at half the temperature.

When you combine the behavior of these two "temperatures," you get the exact same result as the complex 4-state particle. It's like saying, "The average temperature of a room with a heater and an AC unit is the same as a room with a single thermostat set to a weird middle value."

The Big Picture: Building Blocks

The paper also suggests this isn't just a one-time trick. If you have a particle that can hold 6 guests, or 8 guests, you can break them down into smaller, simpler pieces (like prime numbers).

• 4-state breaks into 2 fermions.
• 8-state breaks into 3 different types of fermions.
• 6-state breaks into a 3-state particle and a 2-state fermion.

The "Why Should I Care?"

Why do we care about these weird particles?

1. Quantum Computing: These particles are the "holy grail" for building super-stable quantum computers. They are naturally resistant to errors (noise) that usually destroy quantum information.
2. Simplicity: Before this paper, simulating these particles on a computer was a nightmare. Now, scientists can use standard, fast computer programs designed for regular electrons to simulate these exotic particles. It turns a super-computer problem into a laptop problem.

Summary

Think of this paper as a Rosetta Stone for a strange new language.

• The Problem: 4-state particles are weird and hard to calculate.
• The Discovery: They are secretly just two regular particles playing a game of "double or nothing."
• The Result: We can now easily predict how they behave, how they store energy, and how they might help us build the quantum computers of the future.

The author has essentially handed us a map, showing us that the mysterious "Goldilocks" particles aren't aliens after all; they are just regular fermions wearing a disguise.

Technical Summary

Here is a detailed technical summary of the paper "Noninteracting tight-binding models for Fock parafermions" by Edward McCann.

1. Problem Statement

Fock parafermions are quantum particles with $p$-state occupation numbers ($0, 1, \dots, p-1$) that exhibit statistics intermediate between fermions and bosons (Gentile statistics). While they are promising candidates for topological quantum computing, they typically arise in strongly correlated systems and lack a general single-particle description.

• The Challenge: For arbitrary $p$, the relationship between the number operator and creation/annihilation operators is nonlinear, making standard tight-binding models (which rely on bilinear Hamiltonians) unsolvable for Fock parafermions.
• The Goal: To construct a solvable, non-interacting tight-binding model for Fock parafermions (specifically for $p=4$) that possesses a real single-particle spectrum, allowing for the calculation of thermodynamic properties and energy levels via standard diagonalization techniques.

2. Methodology

The author employs a mapping strategy to transform the complex parafermionic problem into a solvable fermionic one.

• Mapping to Spinful Fermions:
  The core innovation is mapping $p=4$ Fock parafermions to a system of spin-1/2 itinerant fermions (spin-up $\uparrow$ and spin-down $\downarrow$).

  • The parafermion occupation number $n_j \in \{0, 1, 2, 3\}$ is decomposed into fermion occupation numbers:
    $$\hat{N}_j = \hat{n}_{j\uparrow} + 2\hat{n}_{j\downarrow}$$
    where $\hat{n}_{j\sigma} \in \{0, 1\}$. This linear combination allows the four states ($0,1,2,3$) to be represented by the four combinations of two spin species ($00, 10, 01, 11$).
  • Explicit non-local transformations (involving string phase factors) are defined to map the parafermion operators ($c_j, c_j^\dagger$) to fermion operators ($f_{j\sigma}, f_{j\sigma}^\dagger$).

• Hamiltonian Construction:
  A tight-binding Hamiltonian is constructed for the parafermions. Through the mapping, this Hamiltonian transforms into a sum of two bilinear (quadratic) fermionic Hamiltonians:
  $$H = \Phi_\uparrow^\dagger \mathbf{H} \Phi_\uparrow + 2 \Phi_\downarrow^\dagger \mathbf{H} \Phi_\downarrow$$
  where $\mathbf{H}$ is an $L \times L$ matrix representing the single-particle tight-binding model (with onsite energies $u_j$ and hopping $t_j$).

  • Crucially, the "string" phases that usually make parafermion models non-local or non-Hermitian are shown to be gauged away for nearest-neighbor hopping, resulting in a standard Hermitian matrix $\mathbf{H}$.

• Generalization:
  The method is generalized to arbitrary composite $p$. If $p = q_1 q_2 \dots q_k$ (prime factors), the $p$-state parafermions decompose into a linear combination of $q_m$-state parafermions. Specifically, if $p$ is a power of two, the system decomposes entirely into fermions.

3. Key Contributions

1. Solvable Model for $p=4$: The paper provides the first explicit construction of a non-interacting tight-binding Hamiltonian for four-state Fock parafermions that yields a real, single-particle spectrum.
2. Exact Mapping: It establishes a rigorous mathematical mapping between $p=4$ parafermions and two species of spin-1/2 fermions, proving that the many-body energy spectrum of the parafermions is a linear combination of the fermionic spectra.
3. Thermodynamic Consistency: The work demonstrates that the thermodynamic distribution function for these parafermions (Gentile statistics) is mathematically consistent with the fermionic mapping.
4. Extension to Topological Phases: The methodology is extended to the Kitaev superconducting chain, showing that the parafermionic counterpart exhibits a four-fold degenerate ground state in the topological phase, distinguished by the occupation of Majorana edge modes.

4. Key Results

• Energy Spectrum: The many-body energy $E$ for a system with $L$ orbitals is given by:
  $$E = \sum_{\ell=1}^L n_\ell \epsilon_\ell$$
  where $\epsilon_\ell$ are the eigenvalues of the single-particle matrix $\mathbf{H}$ (identical to the fermionic case), and $n_\ell \in \{0, 1, 2, 3\}$.
• Thermodynamics (Gentile Statistics):
  • The mean occupation number $\langle n_\ell \rangle$ follows Gentile statistics.
  • The distribution is shown to be the sum of two Fermi-Dirac distributions:
    $$\langle n_\ell \rangle = \langle n_{\ell\uparrow} \rangle + 2\langle n_{\ell\downarrow} \rangle$$
  • The spin-down component behaves as a Fermi-Dirac gas at an effective temperature $T_{eff} = T/2$.
• Specific Heat and Internal Energy:
  • The internal energy $u_E(T)$ and specific heat $c_V(T)$ are derived as sums of fermionic contributions: $u_E(T) = u_0(T) + 2u_0(T/2)$.
  • At low temperatures, the specific heat of the parafermions is 1.5 times that of non-degenerate fermions ($c_V \approx \frac{3}{2} c_{fermion}$).
• Numerical Verification: Exact diagonalization of a 2-site system confirms that the 16 many-body energy levels match the linear combinations of the single-particle eigenvalues derived from the matrix $\mathbf{H}$.

5. Significance

• Bridging Theory and Experiment: By reducing a complex parafermionic problem to a solvable fermionic tight-binding model, the paper provides a concrete platform for simulating parafermions using existing fermionic systems (e.g., quantum dots or cold atoms).
• Topological Quantum Computing: The demonstration of a four-fold degenerate ground state in the parafermionic Kitaev chain reinforces the potential of parafermions for topological quantum computation, offering a new route to realize non-Abelian anyons.
• Statistical Mechanics: The work clarifies the nature of Gentile statistics, showing how intermediate statistics can emerge from the superposition of fermionic species with different effective temperatures.
• Scalability: The decomposition principle for composite $p$ suggests that a wide class of parafermionic models can be solved analytically, provided $p$ is a composite number (specifically a power of two for full fermionic decomposition).

In summary, McCann successfully constructs a "free" parafermion model by exploiting the algebraic structure of composite occupation numbers, transforming a non-linear, interacting-looking problem into a linear, solvable fermionic one. This opens the door to calculating exact thermodynamic and topological properties of parafermionic systems.