Tangent-Fermion Lattice Formulation

• Tangent-Fermion Lattice is defined by replacing the sine dispersion with a tangent function in the Dirac operator to achieve a single Dirac cone and eliminate fermion doubling.
• It employs a generalized eigenvalue problem with local, Hermitian operators that maintain exact chiral and time-reversal symmetries on the lattice.
• The approach has been extended to handle boundary conditions, magnetic fields, interacting one-dimensional systems, and quantum-circuit implementations, demonstrating its versatility.

Searching arXiv for the cited tangent-fermion papers to ground the article in the current literature. {"query":"tangent fermions Dirac lattice Beenakker Pacholski Donis Vela arXiv", "max_results": 10} {"query":"(Vela et al., 2024)", "max_results": 5} A tangent-fermion lattice is a lattice formulation of the Dirac equation in which the continuum linear dispersion is replaced by a tangent dispersion, while the dynamics are recast as a generalized eigenvalue problem with local Hermitian operators. In the two-dimensional setting, this construction was introduced to place a massless Dirac or Majorana fermion on a lattice with only a single Dirac cone in the Brillouin zone, preserving chiral and time-reversal symmetries and avoiding fermion doubling (Beenakker et al., 2023). Subsequent work extended the same framework to boundary-value problems for Dirac edge states (Vela et al., 2024), to the magnetic-field problem and the zeroth Landau level (Vela et al., 19 May 2025), to interacting one-dimensional chiral and helical systems (Zakharov et al., 14 Jan 2026, Zakharov et al., 23 Jun 2026), and to quantum-circuit block encodings of the generalized Dirac operator (Beenakker, 17 Jun 2026).

1. Definition and formal construction

In the two-dimensional massless case, the continuum starting point is the Dirac equation

$$-i\hbar v(\sigma_x \partial_x + \sigma_y \partial_y)\psi(x,y) = E \psi(x,y),$$

with $\psi=(\psi_1,\psi_2)^T$. A naive nearest-neighbor discretization replaces $\partial_x$ by a symmetric finite difference and produces a sine dispersion $E(k)\propto \sin(ak)$, which has zeros both at $k=0$ and $k=\pi/a$. In two dimensions this generates extra Dirac cones at Brillouin-zone boundaries and hence fermion doubling (Vela et al., 2024, Beenakker et al., 2023).

The tangent-fermion construction replaces the sine dispersion by a tangent dispersion. In the static two-dimensional formulation, the Bloch Hamiltonian is

$$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$

with dispersion

$$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$

Near $\mathbf{k}=0$, $\tan(ak_\alpha/2)\approx ak_\alpha/2$, so the continuum cone is recovered. At the Brillouin-zone boundary the tangent diverges, giving a pole rather than an additional zero (Beenakker et al., 2023).

The central technical step is to write the Dirac problem as a generalized eigenvalue problem

$\psi=(\psi_1,\psi_2)^T$0

with

$\psi=(\psi_1,\psi_2)^T$1

$\psi=(\psi_1,\psi_2)^T$2

Both $\psi=(\psi_1,\psi_2)^T$3 and $\psi=(\psi_1,\psi_2)^T$4 are local, Hermitian, and sparse; $\psi=(\psi_1,\psi_2)^T$5 is positive definite. For Bloch states this yields

$\psi=(\psi_1,\psi_2)^T$6

This is the characteristic tangent-fermion lattice dispersion: one Dirac cone per Brillouin zone, with the nonlinearity pushed to the zone edge as a tangent pole (Vela et al., 2024).

2. Fermion doubling, locality, and symmetry

The tangent-fermion lattice was developed against the background of the Nielsen–Ninomiya obstruction. In the standard lattice-Hamiltonian setting, a local discretization preserving the relevant symmetries cannot realize a single massless Dirac cone without doublers. Tangent fermions evade this obstruction in two linked ways: the underlying tangent derivative is nonlocal, and the practical formulation uses a generalized eigenvalue problem with two local operators rather than a single local Hamiltonian (Beenakker et al., 2023, Beenakker, 17 Jun 2026).

In one dimension, Stacey’s derivative is

$\psi=(\psi_1,\psi_2)^T$7

with Fourier symbol $\psi=(\psi_1,\psi_2)^T$8. This has only one zero in the Brillouin zone and a pole at $\psi=(\psi_1,\psi_2)^T$9, so there is no second low-energy cone. Pacholski et al. showed that the same spectrum can be implemented through a local operator pencil $\partial_x$0, with

$\partial_x$1

so that

$\partial_x$2

This factorization retains the tangent dispersion while restoring locality at the level of the generalized eigenproblem (Beenakker, 17 Jun 2026).

A defining feature of the two-dimensional tangent-fermion lattice is exact preservation of chiral symmetry and time-reversal symmetry. In the square-lattice formulation, $\partial_x$3, and under $\partial_x$4 with $\partial_x$5, both $\partial_x$6 and $\partial_x$7 remain invariant (Vela et al., 2024). The earlier tangent-fermion analysis identifies this symmetry content as the reason the Dirac cone is topologically protected against disorder and staggered perturbations: attempts to couple the Dirac point at $\partial_x$8 to Brillouin-zone-boundary states encounter the tangent divergence, so a gap cannot be opened by the same mechanism that destabilizes sine-dispersion fermions (Beenakker et al., 2023).

The literature distinguishes this construction sharply from Wilson and staggered formulations. Wilson fermions remove doublers by adding a momentum-dependent mass term, but this breaks chiral and time-reversal symmetries. Staggered schemes reduce doubling but retain multiple Dirac points in two dimensions and do not preserve the full symmetry structure of a single topological-insulator surface cone (Beenakker et al., 2023, Vela et al., 2024).

3. Boundary conditions and single-cone Dirac edge states

A major extension of the tangent-fermion lattice is the treatment of boundary conditions for confined Dirac fermions. In the continuum, current conservation and self-adjointness do not permit simply setting $\partial_x$9 on the boundary. For a straight edge at $E(k)\propto \sin(ak)$0, the admissible one-parameter family is

$E(k)\propto \sin(ak)$1

Special cases include the infinite-mass boundary condition at $E(k)\propto \sin(ak)$2 and the zigzag boundary condition at $E(k)\propto \sin(ak)$3 (Vela et al., 2024).

On the surface of a three-dimensional topological insulator, these boundary conditions arise from a magnetic insulator. A magnetization

$E(k)\propto \sin(ak)$4

adds a large term $E(k)\propto \sin(ak)$5, and in the limit $E(k)\propto \sin(ak)$6 the surface Dirac fermion acquires exactly the boundary condition above. Changing the magnetization direction sweeps out the full one-parameter family (Vela et al., 2024).

The lattice implementation proceeds by restricting the infinite-lattice operators $E(k)\propto \sin(ak)$7 and $E(k)\propto \sin(ak)$8 to a finite domain, rotating the spinor basis locally at each boundary site with a unitary $E(k)\propto \sin(ak)$9 chosen so that

$k=0$0

and then removing the spin-down rows and columns at the boundary sites. The resulting matrices $k=0$1 and $k=0$2 satisfy

$k=0$3

Because only principal submatrices are taken and the intermediate transformation is unitary, Hermiticity and positive definiteness are preserved (Vela et al., 2024).

In channel geometry, with width $k=0$4 and conserved longitudinal momentum $k=0$5, the continuum quantization condition is

$k=0$6

with $k=0$7 and $k=0$8. The tangent-fermion strip spectrum reproduces the expected cases: no edge states for infinite-mass boundaries, a flat zero-energy edge band for zigzag boundaries, and dispersive edge states for intermediate $k=0$9. The low-energy lattice spectra agree almost exactly with the analytic continuum results, while remaining built from a single Dirac cone in the Brillouin zone (Vela et al., 2024).

The same study also identifies limitations. For oblique boundaries, the simple sharp-boundary prescription can couple low-energy states to the tangent pole in the folded Brillouin zone and generate spurious oscillations; a finite-mass boundary layer remedies this. For perfectly zigzag boundaries, the flat edge band exhibits a doubled degeneracy, and a small perturbation of $k=\pi/a$0 splits it into one physical and one spurious branch (Vela et al., 2024).

4. Magnetic fields and the zeroth Landau level

The tangent-fermion lattice also has a gauge-invariant magnetic-field formulation. On a square lattice, the field is introduced through gauge-covariant translation operators

$k=\pi/a$1

with

$k=\pi/a$2

The lattice generalized eigenproblem uses

$k=\pi/a$3

$k=\pi/a$4

This formulation preserves exact lattice chiral symmetry, $k=\pi/a$5 and $k=\pi/a$6 (Vela et al., 19 May 2025).

In the continuum, the massless Dirac equation in a perpendicular field $k=\pi/a$7 has Landau levels

$k=\pi/a$8

and the zeroth Landau level is exactly at zero energy, independent of $k=\pi/a$9, with definite chirality

$$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$0

On an infinite lattice, however, Stacey’s theorem enforces equal numbers of zero modes of opposite chirality, so the zeroth level becomes doubly degenerate and loses its continuum-style topological protection (Vela et al., 19 May 2025).

The resolution proposed in the 2025 work is boundary-condition-assisted chirality selection. On a finite domain, the boundary condition

$$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$1

removes the spin-down degree of freedom at boundary sites and changes the zero-mode count. In channel geometry at $$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$2, the finite tangent-fermion lattice acquires a doubly degenerate zero-energy mode consisting of two spin-up edge states, giving

$$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$3

Turning on a uniform magnetic field then produces a zero-energy sector composed of a spin-polarized edge state and a spin-polarized bulk Landau state localized around $$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$4 (Vela et al., 19 May 2025).

The central robustness result is that, in a non-uniform magnetic field preserving chiral symmetry, higher Landau levels broaden but the zero mode remains exactly flat. In the fully two-dimensional rectangular geometry with $$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$5 and $$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$6, all $$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$7 Landau levels broaden into smeared steps in the integrated density of states, while the zero-mode plateau remains perfectly sharp. The zero-energy band has total degeneracy $$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$8 and consists of bulk zeroth-Landau-level states together with additional spin-polarized edge states, all of the same chirality (Vela et al., 19 May 2025).

5. One-dimensional interacting tangent fermions

The tangent-fermion idea was extended to one-dimensional helical and chiral systems by replacing the usual sine dispersion with

$$H_{\rm Stacey} = \frac{2\hbar v}{a}\bigl[\sigma_x\tan(ak_x/2)+\sigma_y\tan(ak_y/2)\bigr],$$9

or, in condensed notation,

$$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$0

In real space this corresponds to long-range hopping with alternating sign. For the helical Luttinger liquid,

$$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$1

which preserves time-reversal symmetry while eliminating the doubler at the Brillouin-zone edge (Zakharov et al., 14 Jan 2026).

That work uses the tangent-fermion kinetic term together with forward-scattering and Umklapp interactions to study spontaneous time-reversal symmetry breaking in a helical Luttinger liquid. Bosonization gives the Luttinger parameter

$$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$2

and the two-particle backscattering term becomes relevant for $$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$3. Density-matrix renormalization group calculations on finite tangent-fermion lattices confirm the expected transition: at half filling and with nonzero two-particle Umklapp, the propagator crosses over from Luttinger scaling to exponential decay and transverse spin correlators saturate, indicating a gapped phase with spontaneous time-reversal symmetry breaking (Zakharov et al., 14 Jan 2026).

A distinct one-dimensional application is the anomaly-free 3–4–5–0 chiral model. There the free Hamiltonian is

$$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$4

again with

$$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$5

Because $$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$6 is strictly monotonic on $$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$7, each species has a single chiral branch and no mirror node at $$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$8 (Zakharov et al., 23 Jun 2026).

The interacting terms include the six-fermion 3–4–5–0 gapping operators and a Hubbard-type density-density interaction that tunes the Luttinger parameter to

$$E_{\rm Stacey}^2(\mathbf{k}) = \left(\frac{2\hbar v}{a}\right)^2 \bigl[\tan^2(ak_x/2)+\tan^2(ak_y/2)\bigr].$$9

The scaling dimension of the six-fermion interaction becomes

$\mathbf{k}=0$0

so in the symmetric case the interaction is relevant for

$\mathbf{k}=0$1

The DMRG results show the opening of an excitation gap in this regime without the appearance of a degenerate ground state, which the paper identifies as the hallmark of symmetric mass generation (Zakharov et al., 23 Jun 2026).

In both one-dimensional studies, the computational practicality of tangent fermions relies on the fact that the nonlocal hopping admits an exact matrix-product-operator representation with finite bond dimension independent of system size (Zakharov et al., 14 Jan 2026, Zakharov et al., 23 Jun 2026). This suggests that the generalized-eigenproblem logic of tangent fermions is not limited to single-particle lattice Dirac equations.

6. Quantum algorithms, computational representation, and scope

A recent development concerns quantum-circuit decomposition of the tangent-fermion Dirac operator. Direct linear-combination-of-unitaries representations of nonlocal discretizations such as the tangent derivative require a number of terms and a subnormalization factor that grow with lattice size. The generalized-eigenproblem formulation avoids this by block-encoding each member of the operator pencil separately (Beenakker, 17 Jun 2026).

In one dimension, with antiperiodic translation operator $\mathbf{k}=0$2, the local operators are

$\mathbf{k}=0$3

These have exact LCU decompositions with

$\mathbf{k}=0$4

$\mathbf{k}=0$5

In $\mathbf{k}=0$6 dimensions,

$\mathbf{k}=0$7

and

$\mathbf{k}=0$8

all independent of lattice size (Beenakker, 17 Jun 2026).

This provides an efficient block-encoding primitive for generalized eigenvalue solvers and quantum linear-system methods applied to Dirac spectra and Green functions without fermion doubling. The paper’s characterization is that the tangent-fermion pencil has complexity “on a par with elliptic operators” in the LCU sense (Beenakker, 17 Jun 2026).

Across these developments, the term “tangent-fermion lattice” has a precise and consistent meaning. It does not refer to a curved-space tangent bundle; in the one-dimensional chiral-fermion literature, the phrase means a lattice whose discretized Dirac operator has a tangent dispersion instead of the standard sine dispersion (Zakharov et al., 23 Jun 2026). In the foundational two-dimensional work, it denotes a symmetry-preserving lattice regularization of a single Dirac cone, implemented through a local generalized eigenproblem and extendable to magnetic fields, boundaries, edge modes, Landau levels, interacting one-dimensional phases, and quantum algorithms (Beenakker et al., 2023, Vela et al., 2024, Vela et al., 19 May 2025, Beenakker, 17 Jun 2026).

A plausible implication is that tangent fermions occupy a distinct position among lattice-fermion regularizations: they trade the standard single-operator eigenproblem for an operator pencil, and in return realize a single Dirac cone without Wilson mass terms, preserve exact lattice chiral symmetry, and remain compatible with sparse-matrix, tensor-network, and block-encoding methods.