No Arabic abstract
Using a recently introduced tensor network method, we study the density of states of the lattice Schwinger model, a standard testbench for lattice gauge theory numerical techniques, but also the object of recent experimental quantum simulations. We identify regimes of parameters where the spectrum appears to be symmetric and displays the expected continuum properties even for finite lattice spacing and number of sites. However, we find that for moderate system sizes and lattice spacing of $gasim O(1)$, the spectral density can exhibit very different properties with a highly asymmetric form. We also explore how the method can be exploited to extract thermodynamic quantities.
Numerical investigations of the XY model, the Heisenberg model and the J-J Heisenberg model are conducted, using the exact diagonalisation, the numerical renormalisation and the density matrix renormalisation group approach. The low-lying energy levels are obtained and finite size scaling is performed to estimate the bulk limit values. The results are found to be consistent with the exact values. The DMRG results are found to be most promising. The Schwinger model is also studied using the exact diagonalisation and the strong coupling expansion. The massless, the massive model and the model with a background electric field are explored. Ground state energy, scalar and vector particle masses and order parameters are examined. The achieved values are observed to be consistent with previous results and theoretical predictions. Path to the future studies is outlined.
We numerically study the single-flavor Schwinger model with a topological $theta$-term, which is practically inaccessible by standard lattice Monte Carlo simulations due to the sign problem. By using numerical methods based on tensor networks, especially the one-dimensional matrix product states, we explore the non-trivial $theta$-dependence of several lattice and continuum quantities in the Hamiltonian formulation. In particular, we compute the ground-state energy, the electric field, the chiral fermion condensate, and the topological vacuum susceptibility for positive, zero, and even negative fermion mass. In the chiral limit, we demonstrate that the continuum model becomes independent of the vacuum angle $theta$, thus respecting CP invariance, while lattice artifacts still depend on $theta$. We also confirm that negative masses can be mapped to positive masses by shifting $thetarightarrow theta +pi$ due to the axial anomaly in the continuum, while lattice artifacts non-trivially distort this mapping. This mass regime is particularly interesting for the (3+1)-dimensional QCD analog of the Schwinger model, the sign problem of which requires the development and testing of new numerical techniques beyond the conventional Monte Carlo approach.
The massive Schwinger model is studied, using a density matrix renormalization group approach to the staggered lattice Hamiltonian version of the model. Lattice sizes up to 256 sites are calculated, and the estimates in the continuum limit are almost two orders of magnitude more accurate than previous calculations. Colemans picture of `half-asymptotic particles at background field (theta = pi) is confirmed. The predicted phase transition at finite fermion mass (m/g) is accurately located, and demonstrated to belong in the 2D Ising universality class.
The massive Schwinger model is studied, using a density matrix renormalisation group approach to the staggered lattice Hamiltonian version of the model. Lattice sizes up to 256 sites are calculated, and the estimates in the continuum limit are almost two orders of magnitude more accurate than previous calculations. Colemans picture of `half-asymptotic particles at background field theta = pi is confirmed. The predicted phase transition at finite fermion mass (m/g) is accurately located, and demonstrated to belong in the 2D Ising universality class.
A trial application of the method of Feynman-Kleinert approximants is made to perturbation series arising in connection with the lattice Schwinger model. In extrapolating the lattice strong-coupling series to the weak-coupling continuum limit, the approximants do not converge well. In interpolating between the continuum perturbation series at large fermion mass and small fermion mass, however, the approximants do give good results. In the course of the calculations, we picked up and rectified an error in an earlier derivation of the continuum series coefficients.