Search

We perform the state-of-the-art tensor network simulations directly in the thermodynamic limit to clarify the critical properties of the $q$-state clock model on the square lattice. We determine accurately the two phase transition temperatures through the singularity of the classical analog of the entanglement entropy, and provide extensive numerical evidences to show that both transitions are of the Berezinskii-Kosterlitz-Thouless (BKT) type for $q\ge 5$ and that the low-energy physics of this model is well described by the $\mathbb{Z}_q$-deformed sine-Gordon theory. We also determine the characteristic conformal parameters, especially the compactification radius, that govern the critical properties of the intermediate BKT phase., 5 pages, 5 figures plus supplemental material; accepted for publication in Phys. Rev. E (Rapid Com.), minor typo corrections

Tensor renormalization group (TRG) constitutes an important methodology for accurate simulations of strongly correlated lattice models. Facilitated by the automatic differentiation technique widely used in deep learning, we propose a uniform framework of differentiable TRG ($\partial$TRG) that can be applied to improve various TRG methods, in an automatic fashion. Essentially, $\partial$TRG systematically extends the concept of second renormalization [PRL 103, 160601 (2009)] where the tensor environment is computed recursively in the backward iteration, in the sense that given the forward process of TRG, $\partial$TRG automatically finds the gradient through backpropagation, with which one can deeply "train" the tensor networks. We benchmark $\partial$TRG in solving the square-lattice Ising model, and demonstrate its power by simulating one- and two-dimensional quantum systems at finite temperature. The deep optimization as well as GPU acceleration renders $\partial$TRG manybody simulations with high efficiency and accuracy., 13 pages, 9 figures

The Chebyshev expansion offers a numerically efficient and easy-implement algorithm for evaluating dynamic correlation functions using matrix product states (MPS). In this approach, each recursively generated Chebyshev vector is approximately represented by an MPS. However, the recurrence relations of Chebyshev polynomials are broken by the approximation, leading to an error which is accumulated with the increase of the order of expansion. Here we propose a reorthonormalization approach to remove this error introduced in the loss of orthogonality of the Chebyshev polynomials. Our approach, as illustrated by comparison with the exact results for the one-dimensional XY and Heisenberg models, improves significantly the accuracy in the calculation of dynamical correlation functions., 9 pages, 8 figures

Using the tensor renormalization group method based on the higher-order singular value decomposition,we have studied the phase transitions of the five-state clock model on the square lattice.The temperature dependence of the specific heat indicates the system has two phase transitions, as verified clearly in the correlation function. By investigating the magnetic susceptibility, we can only obtain the upper critical temperature as Tc2 = 0.9565(7). From the fixed-point tensor, we locate the transition temperatures at Tc1 = 0.9029(1) and Tc2 = 0.9520(1), consistent with the MC and the DMRG results., Comment: 5 pages, 7 figures

In the tensor-network framework, the expectation values of two-dimensional quantum states are evaluated by contracting a double-layer tensor network constructed from initial and final tensor-network states. The computational cost of carrying out this contraction is generally very high, which limits the largest bond dimension of tensor-network states that can be accurately studied to a relatively small value. We propose an optimized contraction scheme to solve this problem by mapping the double-layer tensor network onto an intersected single-layer tensor network. This reduces greatly the bond dimensions of local tensors to be contracted and improves dramatically the efficiency and accuracy of the evaluation of expectation values of tensor-network states. It almost doubles the largest bond dimension of tensor-network states whose physical properties can be efficiently and reliably calculated, and it extends significantly the application scope of tensor-network methods., 7 pages, 6 figures, accepted version for Phys. Rev. B. The title is changed for easier understanding

Frustrated quantum magnetism has moved to the forefront of physics research, posing fundamental questions concerning quantum disordered states, entanglement, topology and the nature of the quantum wavefunction. The defining problem in the field is one of the simplest, the ground state of the nearest-neighbour $S = 1/2$ antiferromagnetic Heisenberg model on the kagome lattice, but has defied all theoretical and numerical methods employed to date. We apply the formalism of tensor-network states (TNS), specifically the method of projected entangled simplex states (PESS), whose combination of a correct accounting for multipartite entanglement and infinite system size provides qualitatively new insight. By studying the ground-state energy, the staggered magnetization we find at all finite tensor bond dimensions and the effects of a second-neighbour coupling, we demonstrate that the ground state is a gapless spin liquid. We discuss the comparison with other numerical studies and the physical interpretation of the gapless ground state., Comment: 8 pages, 4 figures, revised version accepted for publication in PRL

We propose a generalized Lanczos method to generate the many-body basis states of quantum lattice models using tensor-network states (TNS). The ground-state wave function is represented as a linear superposition composed from a set of TNS generated by Lanczos iteration. This method improves significantly both the accuracy and the efficiency of the tensor-network algorithm and allows the ground state to be determined accurately using TNS with very small virtual bond dimensions. This state contains significantly more entanglement than each individual TNS, reproducing correctly the logarithmic size dependence of the entanglement entropy in a critical system. The method can be generalized to non-Hamiltonian systems and to the calculation of low-lying excited states, dynamical correlation functions, and other physical properties of strongly correlated systems., 5 pages, 5 figures

We perform a systematic study of the antiferromagnetic Heisenberg model on the Husimi lattice using numerical tensor-network methods based on Projected Entangled Simplex States (PESS). The nature of the ground state varies strongly with the spin quantum number, $S$. For $S = 1/2$, it is an algebraic (gapless) quantum spin liquid. For $S = 1$, it is a gapped, non-magnetic state with spontaneous breaking of triangle symmetry (a trimerized simplex-solid state). For $S = 2$, it is a simplex-solid state with a spin gap and no symmetry-breaking; both integer-spin simplex-solid states are characterized by specific degeneracies in the entanglement spectrum. For $S = 3/2$, and indeed for all spin values $S \ge 5/2$, the ground states have $120$-degree antiferromagnetic order. In a finite magnetic field, we find that, irrespective of the value of $S$, there is always a plateau in the magnetization at $m = 1/3$., 18 pages, 25 figures;minor changes, refs updated

We develop coarse-graining tensor renormalization group algorithms to compute physical properties of two-dimensional lattice models on finite periodic lattices. Two different coarse-graining strategies, one based on the tensor renormalization group and the other based on the higher-order tensor renormalization group, are introduced. In order to optimize the tensor-network model globally, a sweeping scheme is proposed to account for the renormalization effect from the environment tensors under the framework of second renormalization group. We demonstrate the algorithms by the classical Ising model on the square lattice and the Kitaev model on the honeycomb lattice, and show that the finite-size algorithms achieve substantially more accurate results than the corresponding infinite-size ones., 14 pages, 14 figures

We propose an eigen-operator scheme to study the lattice model of interacting spinless fermions at half filling and show that this model possesses a hidden form of reflection positivity in its Majorana fermion representation. Based on this observation, we prove rigourously that the ground state of this model is either unique or doubly degenerate if the lattice size $N$ is even, and is always doubly degenerate if $N$ is odd. This proof holds in all dimensions with arbitrary lattice structures.

The Potts model plays an essential role in classical statistical mechanics, illustrating many fundamental phenomena. One example is the existence of partially long-range-ordered states, in which some degrees of freedom remain disordered. This situation may arise from frustration of the interactions, but also from an irregular but unfrustrated lattice structure. We study partial long-range order in a range of antiferromagnetic $q$-state Potts models on different two-dimensional lattices and for all relevant values of $q$. We exploit the power of tensor-based numerical methods to evaluate the partition function of these models and hence to extract the key thermodynamic properties -- entropy, specific heat, magnetization, and susceptibility -- giving deep insight into the phase transitions and ordered states of each system. Our calculations reveal a range of phenomena related to partial ordering, including different types of entropy-driven phase transition, the role of lattice irregularity, very large values of the critical $q_c$, and double phase transitions., 23 pages, 27 figures

Using the tensor renormalization group method based on the higher-order singular value decom- position, we have studied the thermodynamic properties of the continuous XY model on the square lattice. The temperature dependence of the free energy, the internal energy and the specific heat agree with the Monte Carlo calculations. From the field dependence of the magnetic susceptibility, we find the Kosterlitz-Thouless transition temperature to be 0.8921 \pm 0.0019, consistent with the Monte Carlo as well as the high temperature series expansion results. At the transition temperature, the critical exponent \delta is estimated as 14.5, close to the analytic value by Kosterlitz., Comment: 5 pages, 6 figures

We propose a new class of tensor-network states, which we name projected entangled simplex states (PESS), for studying the ground-state properties of quantum lattice models. These states extend the pair-correlation basis of projected entangled pair states (PEPS) to a simplex. PESS are an exact representation of the simplex solid states and provide an efficient trial wave function that satisfies the area law of entanglement entropy. We introduce a simple update method for evaluating the PESS wave function based on imaginary-time evolution and the higher-order singular-value decomposition of tensors. By applying this method to the spin-1/2 antiferromagnetic Heisenberg model on the kagome lattice, we obtain accurate and systematic results for the ground-state energy, which approach the lowest upper bounds yet estimated for this quantity., 12 pages, 8 figures

We consider the sign problem for classical spin models at complex $\beta =1/g_0^2$ on $L\times L$ lattices. We show that the tensor renormalization group method allows reliable calculations for larger Im$\beta$ than the reweighting Monte Carlo method. For the Ising model with complex $\beta$ we compare our results with the exact Onsager-Kaufman solution at finite volume. The Fisher zeros can be determined precisely with the TRG method. We check the convergence of the TRG method for the O(2) model on $L\times L$ lattices when the number of states $D_s$ increases. We show that the finite size scaling of the calculated Fisher zeros agrees very well with the Kosterlitz-Thouless transition assumption and predict the locations for larger volume. The location of these zeros agree with Monte Carlo reweighting calculation for small volume. The application of the method for the O(2) model with a chemical potential is briefly discussed., Comment: 8 pages, 9 figures

We introduce a slave-fermion formulation in which to study the charge dynamics of the half-filled Hubbard model on the square lattice. In this description, the charge degrees of freedom are represented by fermionic holons and doublons and the Mott-insulating characteristics of the ground state are the consequence of holon–doublon bound-state formation. The bosonic spin degrees of freedom are described by the antiferromagnetic Heisenberg model, yielding long-ranged (Neel) magnetic order at zero temperature. Within this framework and in the self-consistent Born approximation, we perform systematic calculations of the average double occupancy, the electronic density of states, the spectral function and the optical conductivity. Qualitatively, our method reproduces the lower and upper Hubbard bands, the spectral-weight transfer into a coherent quasiparticle band at their lower edges and the renormalisation of the Mott gap, which is associated with holon–doublon binding, due to the interactions of both quasiparticle species with the magnons. The zeros of the Green function at the chemical potential give the Luttinger volume, the poles of the self-energy reflect the underlying quasiparticle dispersion with a spin-renormalised hopping parameter and the optical gap is directly related to the Mott gap. Quantitatively, the square-lattice Hubbard model is one of the best-characterised problems in correlated condensed matter and many numerical calculations, all with different strengths and weaknesses, exist with which to benchmark our approach. From the semi-quantitative accuracy of our results for all but the weakest interaction strengths, we conclude that a self-consistent treatment of the spin-fluctuation effects on the charge degrees of freedom captures all the essential physics of the antiferromagnetic Mott–Hubbard insulator. We remark in addition that an analytical approximation with these properties serves a vital function in developing a full understanding of the fundamental physics of the Mott state, both in the antiferromagnetic insulator and at finite temperatures and dopings.

We propose a modified spin-wave theory to study the 1/3 magnetization plateau of the antiferromagnetic Heisenberg model on the kagome lattice. By the self-consistent inclusion of quantum corrections, the 1/3 plateau is stabilized over a broad range of magnetic fields for all spin quantum numbers, S. The values of the critical magnetic fields and the widths of the magnetization plateaus are fully consistent with recent numerical results from exact diagonalization and infinite projected entangled paired states., Comment: 5 pages, 4 figures, accepted in Chinese Physics Letters

We propose a novel coarse graining tensor renormalization group method based on the higher-order singular value decomposition. This method provides an accurate but low computational cost technique for studying both classical and quantum lattice models in two- or three-dimensions. We have demonstrated this method using the Ising model on the square and cubic lattices. By keeping up to 16 bond basis states, we obtain by far the most accurate numerical renormalization group results for the 3D Ising model. We have also applied the method to study the ground state as well as finite temperature properties for the two-dimensional quantum transverse Ising model and obtain the results which are consistent with published data., Comment: 11 pages, 16 figures

We show that the matrix (or more generally tensor) product states in a finite translation invariant system can be accurately constructed from the same set of local matrices (or tensors) that are determined from an infinite lattice system in one or higher dimensions. This provides an efficient approach for studying translation invariant tensor product states in finite lattice systems. Two methods are introduced to determine these size-independent local tensors.

We have discussed the tensor-network representation of classical statistical or interacting quantum lattice models, and given a comprehensive introduction to the numerical methods we recently proposed for studying the tensor-network states/models in two dimensions. A second renormalization scheme is introduced to take into account the environment contribution in the calculation of the partition function of classical tensor network models or the expectation values of quantum tensor network states. It improves significantly the accuracy of the coarse grained tensor renormalization group method. In the study of the quantum tensor-network states, we point out that the renormalization effect of the environment can be efficiently and accurately described by the bond vector. This, combined with the imaginary time evolution of the wavefunction, provides an accurate projection method to determine the tensor-network wavfunction. It reduces significantly the truncation error and enable a tensor-network state with a large bond dimension, which is difficult to be accessed by other methods, to be accurately determined., 18 pages 23 figures, minor changes, references added

We propose a second renormalization group method to handle the tensor-network states or models. This method reduces dramatically the truncation error of the tensor renormalization group. It allows physical quantities of classical tensor-network models or tensor-network ground states of quantum systems to be accurately and efficiently determined., Comment: 5 figures, 4 pages