No Arabic abstract
We show a way to perform the canonical renormalization group (RG) prescription in tensor space: write down the tensor RG equation, linearize it around a fixed-point tensor, and diagonalize the resulting linearized RG equation to obtain scaling dimensions. The tensor RG methods have had a great success in producing accurate free energy compared with the conventional real-space RG schemes. However, the above-mentioned canonical procedure has not been implemented for general tensor-network-based RG schemes. We extend the success of the tensor methods further to extraction of scaling dimensions through the canonical RG prescription, without explicitly using the conformal field theory. This approach is benchmarked in the context of the Ising models in 1D and 2D. Based on a pure RG argument, the proposed method has potential applications to 3D systems, where the existing bread-and-butter method is inapplicable.
We introduce a general method for optimizing real-space renormalization-group transformations to study the critical properties of a classical system. The scheme is based on minimizing the Kullback-Leibler divergence between the distribution of the system and the normalized normalizing factor of the transformation parametrized by a restricted Boltzmann machine. We compute the thermal critical exponent of the two-dimensional Ising model using the trained optimal projector and obtain a very accurate thermal critical exponent $y_t=1.0001(11)$ after the first step of the transformation.
An algorithm of the tensor renormalization group is proposed based on a randomized algorithm for singular value decomposition. Our algorithm is applicable to a broad range of two-dimensional classical models. In the case of a square lattice, its computational complexity and memory usage are proportional to the fifth and the third power of the bond dimension, respectively, whereas those of the conventional implementation are of the sixth and the fourth power. The oversampling parameter larger than the bond dimension is sufficient to reproduce the same result as full singular value decomposition even at the critical point of the two-dimensional Ising model.
We show that the recent renormalization-group analysis of Lifshitz critical behavior presented by Leite [Phys. Rev. B {bf 67}, 104415 (2003)] suffers from a number of severe deficiencies. In particular, we show that his approach does not give an ultraviolet finite renormalized theory, is plagued by inconsistencies, misses the existence of a nontrivial anisotropy exponent $theta e 1/2$, and therefore yields incorrect hyperscaling relations. His $epsilon$-expansion results to order $epsilon^2$ for the critical exponents of $m$-axial Lifshitz points are incorrect both in the anisotropic ($0<m<d$) and the isotropic cases ($m=d$). The inherent inconsistencies and the lack of a sound basis of the approach makes its results unacceptable even if they are interpreted in the sense of approximations.
We present a new strategy for contracting tensor networks in arbitrary geometries. This method is designed to follow as strictly as possible the renormalization group philosophy, by first contracting tensors in an exact way and, then, performing a controlled truncation of the resulting tensor. We benchmark this approximation procedure in two dimensions against an exact contraction. We then apply the same idea to a three dimensional system. The underlying rational for emphasizing the exact coarse graining renormalization group step prior to truncation is related to monogamy of entanglement.
Analyzing in detail the first corrections to the scaling hypothesis, we develop accelerated methods for the determination of critical points from finite size data. The output of these procedures are sequences of pseudo-critical points which rapidly converge towards the true critical points. In fact more rapidly than previously existing methods like the Phenomenological Renormalization Group approach. Our methods are valid in any spatial dimensionality and both for quantum or classical statistical systems. Having at disposal fast converging sequences, allows to draw conclusions on the basis of shorter system sizes, and can be extremely important in particularly hard cases like two-dimensional quantum systems with frustrations or when the sign problem occurs. We test the effectiveness of our methods both analytically on the basis of the one-dimensional XY model, and numerically at phase transitions occurring in non integrable spin models. In particular, we show how a new Homogeneity Condition Method is able to locate the onset of the Berezinskii-Kosterlitz-Thouless transition making only use of ground-state quantities on relatively small systems.