Fully Discrete Massive Thirring Model

• Fully discrete Massive Thirring Model is a lattice analogue of the continuous Thirring system that discretizes both space and time while maintaining integrability.
• Distinct formulations such as bilinear/KP-Toda and Lax-pair methods illustrate how coordinate choices and reductions yield various discrete integrable systems.
• The approach enables explicit soliton solutions, Yang–Baxter maps, and links to lattice Hamiltonian formulations, bridging continuous and discrete theories.

Searching arXiv for papers on fully discrete Massive Thirring model and related discretizations. The fully discrete Massive Thirring Model denotes lattice analogues of the classical massive Thirring system in which both independent variables are discretized, most commonly in light-cone coordinates, while the integrable structure is preserved through a discrete zero-curvature condition, Hirota bilinearization, or hierarchy reduction. In the recent literature, the subject is represented most explicitly by a bilinear/KP-Toda construction of a fully discrete light-cone model and by a Lax-pair construction that also produces an associated Yang–Baxter map (Chen et al., 13 Jul 2025, Tsuchida, 2024). Closely related work on semi-discrete systems, Darboux transformations, and lattice Hamiltonian regularizations clarifies that the phrase “fully discrete” is used in more than one technical sense across integrable-systems and lattice-field-theory communities (Joshi et al., 2018, Bañuls et al., 2017).

1. Continuous origin and coordinate dependence

The continuous massive Thirring model appears in several equivalent normalizations and coordinate systems. In one light-cone normalization, the model is

$$\mathrm{i}u_x+v+u|v|^2=0,\qquad \mathrm{i}v_t+u+v|u|^2=0,$$

which is the continuous target of the fully discrete bilinear construction (Chen et al., 13 Jul 2025). In another Lax-pair normalization, the light-cone system is written for fields $q,r,u,v$ as

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$

$iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$

and, under the scalar reduction $r=q^*$, $v=u^*$, becomes

$$iq_t+u-\frac{u}{2}|q|^2=0,\qquad iu_x-q+\frac{q}{2}|u|^2=0.$$

This is the continuous system recovered from the fully discrete Lax construction (Tsuchida, 2024).

The coordinate choice is not a superficial matter. A separate line of work emphasizes laboratory coordinates,

$$\mathrm{i}u_t + \mathrm{i}u_x + v = |v|^2 u, \qquad \mathrm{i}v_t - \mathrm{i}v_x + u = |u|^2 v,$$

and shows that previously known integrable discretizations were mainly in characteristic variables rather than in the coordinates natural for the Cauchy problem (Joshi et al., 2018). This establishes that “fully discrete Massive Thirring Model” does not refer to a unique canonical lattice equation, but to a family of integrable discretizations whose precise form depends on normalization, reduction, and coordinate system.

2. Bilinear and KP-Toda construction of a fully discrete light-cone model

A direct fully discrete light-cone model is given by the pair

$$\frac{\rm i}{a}(u^l_{k+1}-u^l_k)+v^l_{k+1}+u^l_{k+1}v^l_{k+1}\tilde v^l_k=0,$$

$$\frac{\rm i}{b}(v^{l+1}_k-v^l_k)+u^{l+1}_k+v^{l+1}_k u^{l+1}_k\tilde u^l_k=0,$$

where $q,r,u,v$0 is the discrete spatial index, $q,r,u,v$1 is the discrete temporal index, and $q,r,u,v$2 are the lattice spacings in the two directions (Chen et al., 13 Jul 2025). The dependent variables satisfy

$q,r,u,v$3

and the field equations bilinearize into

$q,r,u,v$4

$q,r,u,v$5

$q,r,u,v$6

$q,r,u,v$7

The derivation proceeds from determinant tau functions of the two-component discrete KP-Toda hierarchy, followed by a period-2 reduction

$q,r,u,v$8

and then by a complex-conjugate reduction with

$q,r,u,v$9

The final identification of hierarchy indices,

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$0

produces the fully discrete massive Thirring equations themselves (Chen et al., 13 Jul 2025).

This hierarchy-based route is conceptually distinct from the use of auxiliary discrete indices in continuous KP reductions. In the rogue-wave literature, the tau functions $iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$1 and the shift $iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$2 appear only as internal reduction machinery, and no discrete or fully discrete MT equation is obtained (Chen et al., 2022). The fully discrete model above is therefore a genuine lattice equation rather than a continuous MT system written with auxiliary hierarchy labels.

3. Lax-pair formulation, discrete zero curvature, and Yang–Baxter structure

A second explicit construction starts from a Lax-pair representation in light-cone coordinates and discretizes both directions through discrete spatial and temporal Lax matrices of the same structural type (Tsuchida, 2024). The discrete compatibility condition is

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$3

which is the fully discrete zero-curvature equation. From this identity, the paper derives four rational matrix equations, equations (2.29)–(2.32), for the discrete evolution of $iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$4.

The model admits a Hermitian reduction

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$5

with parameter constraints

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$6

For the scalar case, the parameter choice

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$7

yields the fully discrete scalar MTM, written explicitly in the paper as the rational lattice equations (2.38)–(2.39). The same framework also supplies a binary Bäcklund–Darboux transformation and a one-soliton solution in the scalar reduction (Tsuchida, 2024).

A distinctive feature of this formulation is the associated Yang–Baxter map. Using the factorization identities

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$8

and the analogous identity for the temporal variables, the discrete dynamics is recast as a parameter-dependent Yang–Baxter map with explicit rational formulas, equations (4.3)–(4.6). The map satisfies the set-theoretic Yang–Baxter relation and has a continuous limit, so the fully discrete MTM is tied not only to a lattice zero-curvature representation but also to a factorization-based Yang–Baxter structure (Tsuchida, 2024).

4. Semi-discrete precursors and the role of coordinates

The fully discrete theory is closely connected to semi-discrete precursor models. In laboratory coordinates, an integrable spatial semi-discretization was obtained from a gauge-modified MTM Lax pair together with a Bäcklund–Darboux transformation for the Ablowitz–Ladik lattice, producing the system

$iq_t+u-quu=0,\qquad ir_t-v-vqr=0,$9

supplemented by two difference relations for $iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$0 and $iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$1 (Joshi et al., 2018). That paper explicitly identifies the result as the first integrable semi-discretization of the MTM in laboratory coordinates and emphasizes its relevance to the Cauchy problem.

A related semi-discretization in non-characteristic coordinates is built from the Lax pair

$iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$2

with compatibility

$iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$3

and, after the reduction $iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$4, $iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$5, recovers

$iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$6

in the continuum limit (Tsuchida, 12 May 2025). The appendix identifies the discrete spatial Lax problem with a binary Bäcklund–Darboux transformation, showing why the semi-discrete model is integrable.

The fully discrete bilinear model reduces to the semi-discrete one as

$iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$7

and then to the continuous MT model as

$iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$8

while the fully discrete Lax construction passes first to a semi-discrete MTM in one limit and then to the continuous light-cone system (Chen et al., 13 Jul 2025, Tsuchida, 2024). At the same time, the laboratory-coordinate literature stresses that one cannot simply transfer a characteristic-coordinate discretization back to laboratory coordinates, because the coordinate rotation mixes positive and negative powers of the spectral parameter in the Lax operators (Joshi et al., 2018). This is a central reason why distinct discrete MTM families coexist.

5. Lattice Hamiltonian formulations and the broader use of “fully discrete”

In lattice field theory and quantum many-body work, the phrase “fully discrete” is used in a different but related sense. One study starts from the continuum $iu_x-q+qru=0,\qquad iv_x+r-vqr=0,$9-dimensional Thirring action,

$r=q^*$0

discretizes space with staggered fermions,

$r=q^*$1

maps it to a spin-$r=q^*$2 chain by Jordan–Wigner, and studies the resulting discrete Hamiltonian with MPS/DMRG in the zero-charge sector (Bañuls et al., 2017). The paper explicitly states that the computational study is performed on a fully discrete lattice Hamiltonian formulation.

Two quantum-computation papers make the same distinction from a Hamiltonian perspective. One uses a periodic 1D spatial lattice with Wilson fermions, continuous time, Jordan–Wigner encoding onto $r=q^*$3 qubits, and a three-site benchmark to compute the mass gap by a hybrid classical-quantum method (Mishra et al., 2019). Another uses staggered fermions, Jordan–Wigner mapping, and an explicit $r=q^*$4 Pauli Hamiltonian,

$r=q^*$5

to study finite-temperature chiral and topological transitions with QMETTS (Gong et al., 2024). These works discretize the many-body Hilbert space and spatial coordinate explicitly, but retain continuous time; they are therefore lattice-Hamiltonian realizations of the massive Thirring model rather than fully discrete integrable evolution equations.

6. Solitons, reductions, and structural issues

The fully discrete MTM is primarily valued for preserving integrability under discretization. In the bilinear/KP-Toda construction, integrability is established through Hirota bilinear form, reduction from an integrable hierarchy, and Gram-determinant tau functions that give explicit multi-bright soliton solutions (Chen et al., 13 Jul 2025). The one-soliton solution is written in terms of

$r=q^*$6

and the $r=q^*$7-soliton sector is given by determinant formulas of the same type. The paper also emphasizes that its light-cone discretization is different from the older discretizations of Nijhoff et al. and has no square singularity (Chen et al., 13 Jul 2025).

A technically important point is that, in the discrete bright-soliton solutions, the auxiliary fields $r=q^*$8 are not automatically the complex conjugates of $r=q^*$9. The conjugate structure is recovered only in the continuum limit $v=u^*$0 (Chen et al., 13 Jul 2025). This is one of the most common sources of confusion when discrete MTM formulas are compared directly with the continuous model.

The Lax-pair construction adds a different integrability package: discrete zero curvature, binary Bäcklund–Darboux transformation, one-soliton solutions, and a Yang–Baxter map with continuous limit (Tsuchida, 2024). In parallel, a continuous Darboux-matrix analysis of the classical MTM develops polynomial Darboux matrices constrained by a dihedral $v=u^*$1 reduction, with root quadruplets

$v=u^*$2

and explicit lowest-degree dressing formulas (Degasperis, 2014). That work does not derive a semi-discrete or fully discrete MTM, but it identifies the algebraic ingredients usually required for one: symmetry-preserving Darboux steps, automorphic determinant polynomials, and an iterative dressing mechanism. Within the integrable-systems literature, this suggests that fully discrete MTM equations are best understood not as arbitrary finite-difference replacements, but as reductions of larger discrete Lax or hierarchy structures constrained by the same reduction symmetries as the continuous theory.