Accurate simulation of the finite density lattice Thirring model
Abstract
We present a study of the finite density lattice Thirring model in 1+1 dimensions using the worldline/fermionbag algorithm. The model has features similar to QCD and provides a test case for exploring the accuracy of various methods of solving sign problems. In the massless limit and with open boundary conditions we show that the sign problem is an artifact of the auxiliary field approach and is completely eliminated in the fermion bag approach. With periodic boundary conditions the sign problem is mild in the fermion bag method. We present accurate results for various quantities in the model that can be used as a benchmark for comparison with other methods of solving sign problems.
Accurate simulation of the finite density lattice Thirring model
Jarno Rantaharju^{}^{}email: jmr108@phy.duke.edu
1Duke University,

Abstract.
1 Introduction
In lattice field theories with finite chemical potential, Monte Carlo methods are often hindered by the sign problem. The partition function may not be strictly positive and cannot be interpreted as a probability density. The problem is encountered in QCD, among other models. Several general methods for solving sign problems are currently being studied, including the complex Langevin method and Monte Carlo simulations on Lefschetz thimbles.
In order to verify the accuracy of these methods, it is necessary to compare them to accurate data in models where the sign problem can be solved. Several exactly solvable one dimensional models have been used for this purpose. [1, 2, 3, 4, 5, 6, 7] Two dimensional models offer more complex test cases that can more accurately represent the complex physics of target models, but may not be easily solvable.
The Thirring model in 1+1 dimensions provides a physically significant test case with features similar to QCD. The model was studied in [8] using a holomorphic flow to produce an integration contour similar to Lefschetz thimbles. It is asymptotically free and has a well defined continuum limit at the Gaussian fixed point. It includes a massless boson and a fermion with a dynamically generated mass, reminiscent of the pions and baryons in QCD. However, these features do not result from a spontaneous symmetry breaking as in QCD. Spontaneus breaking of continuous symmetries is forbidden in two dimensional systems. Instead the model is critical at the massless limit [9].
We use the fermion bag algorithm [10] to simulate the massless 1+1 dimensional Thirring model and provide benchmark data on the fermion number at finite density. We show that the sign problem is absent with open spatial boundary conditions using a world line representation. The world line representation also provides an efficient worm algorithm for producing fermion bag configurations. Using open boundary conditions, we provide data with large lattice sizes and study critical scaling in the model.
2 Numerical Studies
The massless staggered Thirring model is defined by the action
(1)  
(2) 
where is the chemical potential. The action can be expressed equivalently using an auxiliary gauge field
(3)  
(4) 
This form of the action was used in [8]. The coupling is related to as
(5) 
where and are modified Bessel functions.
2.1 Fermion Bag Representation
The weight of a configuration of fermion fields is given by the partition function
(6) 
The fermion bag [10] representation is found be expanding the partition function
(7)  
(8) 
where is a new dimer variable that takes values and , [d] is the set of possible dimer configurations and is the set of sites not occupied by a dimer. An example configuration is shown in Figure 1. Each dimer occupies two sites, and .
It is now straightforward to integrate over the fermion fields on occupied sites
(9)  
(10) 
where the matrix is constructed out of the fermion matrix by projecting to the space of unoccupied sites . The weight of a given configuration of dimers is then
(11) 
Configurations of can be generated using heatbath updates. There are two possible updates, adding and removing a dimer. The determinant of the suggested configuration is calculated and the new configuration is accepted with the probability , where is the partition function of the new configuration. If negative weights are encountered, the observables must be reweighed to the actual probability density
(12) 
In order to measure the severity of the sign problem at nonzero chemical potential, we measure the average sign of the partition function. The results are shown in Figure 2 in the left panel. The right panel shows the average sign in the auxiliary field representation for comparison. In both cases we use both periodic and antiperiodic boundary conditions in spatial directions and and . In the fermion bag model we use the coupling and in the auxiliary field model we use the comparable . The temporal boundary condition are always antiperiodic.
The sign problem is mild with periodic boundary conditions. The average sign differs from around two values of the chemical potential, and . With antiperiodic boundary conditions, the sign problem is severe in a wider region between and . It can still be overcome using large statistics in these small volumes.
In a finite volume system the energy levels of possible states are quantized. As the chemical potential is increased the fermion number increases in a step like fashion as the chemical potential exceeds these energy levels. An adequate solution to the sign problem should reproduce this behavior. In Figure 3 we show the fermion number
(13) 
As a function of the chemical potential. The steplike behavior is reproduced in both cases and we find a shift in the energy of the lowest state due to the interaction term, consistent with a dynamically generated mass. In the antiperiodic case we find a region of nontrivial behavior corresponding to the region of small average sign. We also notice that the regions with average sign different from correspond to the jumps in the fermion number in both cases.
2.2 World Line Representation
The sign problem is notably absent when the spatial boundary conditions are open,
(14) 
and the temporal boundary conditions are antiperiodic. This result can be shown exactly using a world line representation of the partition function. The method also allows a fast worm algorithm for producing fermion bag configurations.
The world line representation is an expansion of the exponential form of the fermion determinant
(15)  
(16)  
(17) 
where are new directed link variables, with representing a term with and a term with . Since the Grassmann integral is only nonzero if we have exactly one instance of and at all sites, the links must form closed loops and cannot overlap. In the last step we performed the Grassmann integral over all fermion fields and arranged the sum over the link variables into loops.
The weight of a configuration is then a product of the weights of the loops in the configuration and the weights of the dimers. We will show that all allowed loops have positive weight given open boundary conditions. Consider first loops that do not wrap around the time direction. We can enumerate all such loops by starting from a simple loop and using the two deformations shown in Figure 4.
Deformation (a) replaces a link with a staple or a staple with a link and does not change the sign of the loop. Deformation (b) replaces a corner with an inverted corner and does change the sign of the loop. It also also changes the volume enclosed by one. Repeating these two step we find that all loops that enclose an even number of sites have the same sign. Since the simple twosite loop has a positive sign, any loop that encloses an even volume also has a positive sign. Since odd volumes cannot be filled with dimers or loops, only even volumes are allowed.
The same argument applies for time wrapping loops that wrap around the time direction. These loops split the lattice into two volumes, one on each side. Starting from a loop pointing directly up or directly down, we can enumerate all other loops. The weight of this starting loop is positive and since only even volumes are allowed, the weights of all allowed loops are positive. The argument does not apply if a mass term in included in the action or when the spatial boundary conditions are periodic or antiperiodic. In these cases it is trivial to construct allowed loops with a negative weight.
Since all possible configurations of loops have positive weight, it is straightforward to generate configurations within the world line representation using a worm update. The update begins by changing a random link from to with the probability
(18) 
If accepted, the update creates two defects. There is now a site with no link pointing to it and a site with two links pointing to it. The first is the tail and the second the head of the worm. The new link is kept and a new direction is suggested for the old link pointing at the head. With this update the head propagates to a new site. The process is repeated until a link pointing at the tail is accepted.
The set of dimers may be updated using a simple heatbath step. Loops consisting of two site can be turned into dimers with the probability and dimers into loops with the probability . An example of the worm update is shown in Figure 5.
It is now straightforward to produce data even on large lattices. As before, we reproduce the steplike behavior of the fermion number on finite lattices. The fermion number in the world line formalism is simply
(19) 
Figure 6 shows the fermion number at with several lattice sizes up to and . The steplike behavior of the fermion number is clearly visible. At nonzero coupling , the steps are shifted to a larger , providing evidence of the dynamically generated fermion mass.
The model is asymptotically free and at small , in the continuum limit, it should be possible to connect with perturbation theory. This requires large lattice sizes to reduce finite volume effects. At infinite volume but small coupling , the fermion mass scales as
(20) 
where is the oneloop coefficient of the function.
We measure the fermion mass using large square lattices at small . In this case the fermion number density increases linearly with the chemical potential
(21) 
We show measurements of the fermion number density at in Figure 7 on the left. The fermion mass is in this case. On the right in Figure 7 we show the fermion mass at small coupling and large lattice size, up to , together with the asymptotic scaling with and . We find qualitative agreement with the perturbative expectation.
3 Conclusions
We study the massless 1+1 dimensional lattice Thirring model with finite chemical potential using the fermion bag formulation. With periodic and antiperiodic boundary conditions we find a milder sign problem than in the auxiliary field formalism. We are able extract accurate measurements with small lattice sizes.
We show that the sign problem is removed with open boundary conditions in the spatial direction using a world line formalism. The world line expansion naturally leads to an efficient worm update for generating fermion bag configurations and allows accurate measurements with very large lattice sizes. We then study the model at multiple lattice sizes, robustly reproducing the steplike behavior of the fermion number and showing that a fermion mass is generated dynamically by the interaction. At small coupling and large lattice size we find a qualitative match to the expected asymptotic behavior of the fermion mass.
4 Acknowledgements
This study is done in collaboration with S. Chandrasekharan and V. Ayyar. We thank A. Alexandru and P.F. Bedaque for useful discussion. This work was supported by the U.S. Department of Energy, Office of Science, Nuclear Physics program under Award Number DEFG0205ER41368.
References
 [1] A. Alexandru, G. Basar, P.F. Bedaque, G.W. Ridgway, N.C. Warrington, JHEP 05, 053 (2016), 1512.08764
 [2] A. Mukherjee, M. Cristoforetti, Phys. Rev. B90, 035134 (2014), 1403.5680
 [3] H. Fujii, S. Kamata, Y. Kikukawa, JHEP 11, 078 (2015), [Erratum: JHEP02,036(2016)], 1509.08176
 [4] Y. Tanizaki, Y. Hidaka, T. Hayata, New J. Phys. 18, 033002 (2016), 1509.07146
 [5] H. Fujii, S. Kamata, Y. Kikukawa, JHEP 12, 125 (2015), [Erratum: JHEP09,172(2016)], 1509.09141
 [6] S. Kamata, Y. Kikukawa, H. Fujii, PoS LATTICE2015, 190 (2016), 1511.00160
 [7] A. Alexandru, G. Basar, P. Bedaque, Phys. Rev. D93, 014504 (2016), 1510.03258
 [8] A. Alexandru, G. Basar, P.F. Bedaque, G.W. Ridgway, N.C. Warrington, Phys. Rev. D95, 014502 (2017), 1609.01730
 [9] E. Witten, Nucl. Phys. B145, 110 (1978)
 [10] S. Chandrasekharan, Phys. Rev. D82, 025007 (2010), 0910.5736