The analysis of the error of perturbative series by comparing it to the exact solution is an important tool to understand the non-perturbative physics of statistical models. For some toy models, a new method can be used to calculate higher order weak coupling expansion and modified perturbation theory can be constructed. However, it is nontrivial to generalize the new method to understand the critical behavior of high dimensional spin and gauge models. Actually, it is a big challenge in both high energy physics and condensed matter physics to develop accurate and efficient numerical algorithms to solve these problems. In this thesis, one systematic way named tensor renormalization group method is discussed. The applications of the method to several spin and gauge models on a lattice are investigated. theoretically, the new method allows one to write an exact representation of the partition function of models with local interactions. E.g. O(N) models, Z2 gauge models and U(1) gauge models. Practically, by using controllable approximations, results in both finite volume and the thermodynamic limit can be obtained. Another advantage of the new method is that it is insensitive to sign problems for models with complex coupling and chemical potential. Through the new approach, the Fisher's zeros of the 2D O(2) model in the complex coupling plane can be calculated and the finite size scaling of the results agrees well with the Kosterlitz-Thouless assumption. Applying the method to the O(2) model with a chemical potential, new phase diagram of the models can be obtained. The structure of the tensor language may provide a new tool to understand phase transition properties in general.
We show that the Tensor Renormalization Group (TRG) method can be applied to O(N) spin models, principal chiral models and pure gauge theories (Z2, U(1) and SU(2)) on (hyper) cubic lattices.We explain that contrarily to some common belief, it is very difficult to write compact formulas expressing the blockspinning of lattice models.We show that in contrast to other approaches, the TRG formulation allows us to write exact blocking formulas with numerically controllable truncations.The basic reason is that the TRG blocking separates neatly the degrees of freedom inside the block and which are integrated over, from those kept to communicate with the neighboring blocks.We argue that the TRG is a method that can handle large volumes, which is crucial to approach quasi-conformal systems.The method can also get rid of some sign problems.We discuss recent results regarding the critical properties of the 2D O(2) nonlinear sigma model with complex β and chemical potential.As some of these results appeared in a recently published paper (PRD 88, 056005) and two recent preprints (
Abstract The tensor renormalization group method is a promising approach to lattice field theories, which is free from the sign problem unlike standard Monte Carlo methods. One of the remaining issues is the application to gauge theories, which is so far limited to U(1) and SU(2) gauge groups. In the case of higher rank, it becomes highly nontrivial to restrict the number of representations in the character expansion to be used in constructing the fundamental tensor. We propose a practical strategy to accomplish this and demonstrate it in 2D U(N) and SU(N) gauge theories, which are exactly solvable. Using this strategy, we obtain the singular-value spectrum of the fundamental tensor, which turns out to have a definite profile in the large-N limit. For the U(N) case, in particular, we show that the large-N behavior of the singular-value spectrum changes qualitatively at the critical coupling of the Gross-Witten-Wadia phase transition. As an interesting consequence, we find a new type of volume independence in the large-N limit of the 2D U(N) gauge theory with the θ term in the strong coupling phase, which goes beyond the Eguchi-Kawai reduction.
We report recent progress on the application of the tensor renormalization group (TRG) to quantum field theories pursued by the Tsukuba group. We explain how to treat the scalar, fermion, and gauge theories with the TRG method presenting the results for the phase transitions in the (3+1)-dimensional ((3+1)
d) complex
ϕ4 theory at finite density, (1+1)
d pure U(1) lattice gauge theory with a
θ term, (3+1)
d Nambu--Jona-Lasinio model at finite density and (1+1)
d and (2+1)
d Hubbard models at an arbitrary chemical potential. It is demonstrated that the TRG method is free from the sign problem in practical calculations and applicable to the four-dimensional models.
Using the example of the two-dimensional (2D) Ising model, we show that in contrast to what can be done in configuration space, the tensor renormalization group formulation allows one to write exact, compact, and manifestly local blocking formulas and exact coarse-grained expressions for the partition function. We argue that similar results should hold for most models studied by lattice gauge theorists. We provide exact blocking formulas for several 2D spin models [the
O(2) and
O(3) sigma models and the SU(2) principal chiral model] and for the three-dimensional gauge theories with groups
Z2, U(1) and SU(2). We briefly discuss generalizations to other groups, higher dimensions and practical implementations.
The phase diagram of QCD is a central research topic in modern particle physics. Its calculation, however, is notoriously difficult especially when not only nonzero temperatures
T but also non-vanishing baryon densities are considered – as is relevant for heavy‑ion collisions or inside large neutron stars. At the theoretical level, nonzero baryon densities are achieved by introducing a nonzero chemical potential
μ, which creates a particle-antiparticle asymmetry. It turns out that for nonzero chemical potential, the fermion determinant is no longer real but complex, and therefore standard approaches based on importance sampling break down. In the literature, this is commonly referred to as the sign problem. In order to circumvent this sign problem, a variety of methods has been proposed. Examples include a Taylor expansion in
μ, reweighting, complex Langevin, thimbles or path optimization. Nevertheless, none of the proposed remedies can successfully reach regimes in which
μ/T>1, and therefore even today the full QCD phase diagram in the
μ−T plane cannot be calculated. There is a promising further approach, called tensor-network approach, that differs crucially from all of the aforementioned approaches. In this approach, the partition function is written as a contraction of a large number of local tensors. Of special relevance is the higher order tensor renormalization group (HOTRG) approach, which can be applied for arbitrary number of dimensions, and in particular its generalization, the Grassmann HOTRG (GHOTRG) approach, which is applicable to theories including fermions. Recently, GHOTRG was tailored to strong-coupling QCD using one flavor of staggered quarks in two dimensions. It was observed that, for a non‑zero chemical potential, no numerical instabilities arise. In this thesis, we go beyond both the strong-coupling limit and the single-flavor case. We express the partition function of lattice QCD with staggered fermions in the strong-coupling expansion in terms of a tensor network. This is done for an arbitrary number of dimensions, flavors, and orders of the strong-coupling expansion. Further, the GHOTRG approach is extended such that every coefficient
Zn of the partition function
Z=∑nZnβn can be calculated separately. Finally, numerical results in two dimensions are studied.
We present a tensor-network approach for two-dimensional strong-coupling QCD with staggered quarks at nonzero chemical potential. After integrating out the gauge fields at infinite coupling, the partition function can be written as a full contraction of a tensor network consisting of coupled local numeric and Grassmann tensors. To evaluate the partition function and to compute observables, we develop a Grassmann higher-order tensor renormalization group method, specifically tailored for this model. During the coarsening procedure, the blocking of adjacent Grassmann tensors is performed analytically, and the total number of Grassmann variables in the tensor network is reduced by a factor of two at each coarsening step. The coarse-site numeric tensors are truncated using higher-order singular value decompositions. The method is validated by comparing the partition function, the chiral condensate and the baryon density computed with the tensor method with exact analytical results on small lattices up to volumes of 4×4. For larger volumes, we present first tensor results for the chiral condensate as a function of the mass and volume, and observe that the chiral symmetry is not broken dynamically in two dimensions. We also present tensor results for the number density as a function of the chemical potential, which hint at a first-order phase transition.