\n\nAlong the line developed in a recent paper by one of the authors (Li W., von Delft J. and Xiang T., Phys. Rev. B, 86 (2012) 195137), a simple update scheme, as an optimal method for optimizing the infinite tree tensor network (iTTN) states, is utilized to study the XXZ model under external field on the Bethe lattice. The ground-state energy, bipartite entanglement, order parameters, and nearest-neighbor correlation functions are calculated in the framework of the iTTN states. With increasing external field, two quantum phase transitions (QPTs) (a first-order and a second-order one) take place successively. Between the antiferromagnetic phase and the saturated ferromagnetic phase, there exists another intermediate (canted Néel) phase, in which both x-axis antiferromagnetic order and unsaturated z-axis ferromagnetic order coexist. An interesting spin-flop transition between antiferromagnetic phase and canted Néel phase is observed. The model-independent bipartite entanglement is found to be capable of describing all the QPTs. It is interesting that the Schmidt gap acts as a local order parameter, and can also be used to describe the QPTs.\n
\n\nAlong the line developed in a recent paper by one of the authors (Li W., von Delft J. and Xiang T., Phys. Rev. B, 86 (2012) 195137), a simple update scheme, as an optimal method for optimizing the infinite tree tensor network (iTTN) states, is utilized to study the XXZ model under external field on the Bethe lattice. The ground-state energy, bipartite entanglement, order parameters, and nearest-neighbor correlation functions are calculated in the framework of the iTTN states. With increasing external field, two quantum phase transitions (QPTs) (a first-order and a second-order one) take place successively. Between the antiferromagnetic phase and the saturated ferromagnetic phase, there exists another intermediate (canted Néel) phase, in which both x-axis antiferromagnetic order and unsaturated z-axis ferromagnetic order coexist. An interesting spin-flop transition between antiferromagnetic phase and canted Néel phase is observed. The model-independent bipartite entanglement is found to be capable of describing all the QPTs. It is interesting that the Schmidt gap acts as a local order parameter, and can also be used to describe the QPTs.\n
We show that the simple update approach proposed by Jiang et al. [H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground-state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper bounded by $1/\mathrm{ln}(q\ensuremath{-}1)$ even at the transition point (q is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar ``critical'' phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.
We have extended the canonical tree tensor network (TTN) method, which was initially introduced to simulate the zero-temperature properties of quantum lattice models on the Bethe lattice, to finite temperature simulations. By representing the thermal density matrix with a canonicalized tree tensor product operator, we optimize the TTN and compute accurately the thermodynamic quantities, including the internal energy, specific heat, spontaneous magnetization, etc., at various temperatures. By varying the anisotropic coupling constant $\mathrm{\ensuremath{\Delta}}$, we obtain the phase diagram of the spin-1/2 Heisenberg XXZ model on the Bethe lattice, where three magnetic ordered phases, namely the ferromagnetic, XY, and antiferromagnetic phases, are found in low temperatures and separated from the high-T paramagnetic phase by a continuous thermal phase transition at Tc. The XY phase is separated from the other two phases by two first-order phase transition lines at the symmetric coupling points $\mathrm{\ensuremath{\Delta}}=\ifmmode\pm\else\textpm\fi{}1$. We have also carried out a linear spin wave calculation on the Bethe lattice, showing that the low-energy magnetic excitations are always gapped, and find the obtained magnon gaps in very good agreement with those estimated from the TTN simulations. Despite the gapped excitation spectrum, Goldstone-like transverse fluctuation modes, as a manifestation of spontaneous continuous symmetry breaking, are observed in the ordered magnetic phases with $|\mathrm{\ensuremath{\Delta}}|\ensuremath{\le}1$. One remarkable feature there is that the transverse correlation length reaches ${\ensuremath{\xi}}_{c}=1/ln(z\ensuremath{-}1)$ for $T\ensuremath{\le}{T}_{c}$, the maximal value allowed on a z-coordinated Bethe lattice.
Bipartite entanglement, entanglement spectrum, and Schmidt gap in S=1 bond-alternative antiferromagnetic Heisenberg chain are investigated by the infinite time-evolving block decimation (iTEBD) method. The quantum phase transition (QPT) from the singlet-dimer phase to the Haldane phase can be detected by the singular behavior of bipartite entanglement, the sudden change of the entanglement spectrum, and the completely vanishing of the Schmidt gap. The critical point is determined to be around r c ≃ 0.587, and the second-order character of the QPT is verified. Doubly degenerate entanglement spectra of both even and odd bonds are observed in the Haldane phase, by which one can distinguish the Haldane phase from the singlet-dimer phase easily. Nearest-neighbor antiferromagnetic correlations and next-nearest-neighbor ferromagnetic correlations are found in the whole parameter region. At the critical massless point, although exponentially decaying antiferromagnetic correlation is observed, it approaches to a constant value finally. Therefore, long-range correlations exist and the correlation length becomes divergent at the critical point.
By Lanczos exact diagonalization and the infinite time-evolving block decimation (iTEBD) technique, the two-site entanglement as well as the bipartite entanglement, the ground state energy, the nearest-neighbor correlations, and the magnetization in the antiferromagnetic Heisenberg (AFH) model under an external field are investigated. With increasing external field, the small size system shows some distinct upward magnetization stairsteps, accompanied synchronously with some downward two-site entanglement stairsteps. In the thermodynamic limit, the two-site entanglement, as well as the bipartite entanglement, the ground state energy, the nearest-neighbor correlations, and the magnetization are calculated, and the critical magnetic field h(c) = 2.0 is determined exactly. Our numerical results show that the quantum entanglement is sensitive to the subtle changing of the ground state, and can be used to describe the magnetization and quantum phase transition. Based on the discontinuous behavior of the first-order derivative of the entanglement entropy and fidelity per site, we think that the quantum phase transition in this model should belong to the second-order category. Furthermore, in the magnon existence region (h < 2.0), a logarithmically divergent behavior of block entanglement which can be described by a free bosonic field theory is observed, and the central charge c is determined to be 1.
The characterization of the quantum phase transition in a lowdimensional system has attracted a considerable amount of attention in quantum manybody systems. As one of the fundamental models in quantum magnetism, spin-1 models have richer phase diagrams and show more complex physical phenomena. In the spin-1 antiferromagnetic XXZ model, the Haldane phase and the Nel phase are the gapped topologic phases which cannot be characterized by the local order parameters. To characterize the nature in such phases, one has to calculate the non-local long range order parameters. Normally, the non-local order parameter in the topological phase is obtained from the extrapolation of finite-sized system in numerical study. However, it is difficult to extract the critical exponents with such an extrapolated non-local order parameter due to the numerical accuracy. In a recently developed tensor network representation, i.e., the infinite matrix product state (iMPS) algorithm, it was shown that the non-local order can be directly calculated from a very large lattice distance in an infinite-sized system rather than an extrapolated order parameter in a finite-sized system. Therefore, it is worthwhile using this convenient technique to study the non-local orders in the topological phases and characterize the quantum criticalities in the topological quantum phase transitions. In this paper, by utilizing the infinite matrix product state algorithm based on the tensor network representation and infinite time evolving block decimation method, the quantum entanglement, fidelity, and critical exponents of the topological phase transition are investigated in the one-dimensional infinite spin-1 bond-alternating XXZ Heisenberg model. It is found that there is always a local dimerization order existing in the whole parameter range when the bond-alternative strength parameter changes from 0 to 1. Also, due to the effect of the bond-alternating, there appears a quantum phase transition from the long-rang ordering topological Nel phase to the local ordering dimerization phase. The von Neumann entropy, fidelity per lattice site, and order parameters all give the same phase transition point at c = 0.638. To identify the type of quantum phase transition, the central charge c 0.5 is extracted from the ground state von Neumann entropy and the finite correlation length, which indicates that the phase transition belongs to the two-dimensional Ising universality class. Furthermore, it is found that the Nel order and the susceptibility of Nel order have power-law relations to |-c|. From the numerical fitting of the Nel order and its susceptibility, we obtain the characteristic critical exponents ' = 0.236 and ' = 0.838. It indicates that such critical exponents from our method characterize the nature of the quantum phase transition. Our critical exponents from the iMPS method can provide guidance for studying the properties of the phase transition in quantum spin systems.
The recent progress in the optimization of two-dimensional tensor networks [H.-J. Liao, J.-G. Liu, L. Wang, and T. Xiang, Phys. Rev. X 9, 031041 (2019)] based on automatic differentiation opened the way towards precise and fast optimization of such states and, in particular, infinite projected entangled-pair states (iPEPS) that constitute a generic-purpose Ansatz for lattice problems governed by local Hamiltonians. In this work, we perform an extensive study of a paradigmatic model of frustrated magnetism, the J_1-J_2 J 1 − J 2 Heisenberg antiferromagnet on the square lattice. By using advances in both optimization and subsequent data analysis, through finite correlation-length scaling, we report accurate estimations of the magnetization curve in the N'eel phase for J_2/J_1 \le 0.45 J 2 / J 1 ≤ 0.45 . The unrestricted iPEPS simulations reveal an U(1) U ( 1 ) symmetric structure, which we identify and impose on tensors, resulting in a clean and consistent picture of antiferromagnetic order vanishing at the phase transition with a quantum paramagnet at J_2/J_1 \approx 0.46(1) J 2 / J 1 ≈ 0.46 ( 1 ) . The present methodology can be extended beyond this model to study generic order-to-disorder transitions in magnetic systems.