Fermion discretization effects in the two-flavor lattice Schwinger model: A study with matrix product states

We present a comprehensive tensor network study of staggered, Wilson, and twisted mass fermions in the Hamiltonian formulation, using the massive two-flavor Schwinger model as a benchmark. Particular emphasis is placed on twisted mass fermions, whose properties in this context have not been systematically explored before. We confirm the expected <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mi mathvariant="script">O</a:mi> <a:mo stretchy="false">(</a:mo> <a:mi>a</a:mi> <a:mo stretchy="false">)</a:mo> </a:math> improvement in the free theory and observe that this improvement persists in the interacting case. By leveraging an electric-field-based method for mass renormalization, we reliably tune to maximal twist and establish the method’s applicability in the two-flavor model. Once mass renormalization is included, the pion mass exhibits rapid convergence to the continuum limit. Finite-volume effects are addressed using two complementary approaches: dispersion relation fits and finite-volume scaling. Our results show excellent agreement with semiclassical predictions and reveal a milder volume dependence for twisted mass fermions compared to staggered and Wilson discretizations. In addition, we observe clear isospin-breaking effects, suggesting intriguing parallels with lattice QCD. These findings highlight the advantages of twisted mass fermions for Hamiltonian simulations and motivate their further exploration—particularly in view of future applications to higher-dimensional lattice gauge theories.