We devise a method based on the tensor-network formalism to calculate genuine multisite entanglement in ground states of infinite spin chains containing spin-1/2 or spin-1 quantum particles. The ground state is obtained by employing an infinite time-evolving block decimation method acting upon an initial matrix product state for the infinite spin system. We explicitly show how such infinite matrix product states with translational invariance provide a natural framework to derive the generalized geometric measure, a computable measure of genuine multisite entanglement, in the thermodynamic limit of quantum many-body systems with both spin-1/2 and higher-spin particles.
We devise a method based on the tensor-network formalism to calculate genuine multisite entanglement in ground states of infinite spin chains containing spin-1/2 or spin-1 quantum particles. The ground state is obtained by employing an infinite time-evolving block decimation method acting upon an initial matrix product state for the infinite spin system. We explicitly show how such infinite matrix product states with translational invariance provide a natural framework to derive the generalized geometric measure, a computable measure of genuine multisite entanglement, in the thermodynamic limit of quantum many-body systems with both spin-1/2 and higher-spin particles.
We employ a genuine multipartite entanglement measure, the generalized geometric measure, for investigating the quantum phase transition in an infinite quantum spin-1/2 chain with two-spin as well as three-spin interactions. We show that in contrast to bipartite entanglement of the ground state, the genuine multiparty one unerringly indicates the quantum phase transition in the system. The system lends itself to a complementary behavior between bipartite and genuine multipartite entanglements. We also investigate the trend of the multipartite entanglement in finite chains, in which we perform a finite-size scaling to confirm the signature of the quantum phase transition already in finite systems. The finite-size investigation renders the phenomenon to be observable in real systems.
We compute the multipartite entanglement measures such as the global entanglement of various one- and two-dimensional quantum systems to probe the quantum criticality based on the matrix and tensor product states (MPSs and TPSs). We use the infinite time-evolving block decimation (iTEBD) method to find the ground states numerically in the form of MPSs and TPSs, and then evaluate their entanglement measures by the method of tensor renormalization group (TRG). We find that these entanglement measures can characterize the quantum phase transitions by their derivative discontinuity right at the critical points in all models considered here. We also comment on the scaling behaviors of the entanglement measures by the ideas of quantum-state renormalization group transformations.
The study of interacting quantum many-body systems is one of the central challenges in modern physics. Interestingly, by understanding the physical properties of such systems from a microscopic description, one can explain many complex quantum phenomena such as the high-temperature superconductivity, fractional quantum Hall effect, etc. However, due to high complexity of these systems, theoretical investigations such as analytical and perturbation theories are limited to studying only a small subset of problems. In the past few decades, numerical simulation has emerged as an indispensable method for studying many phenomena, and has led to immense success and progress in recent years. The main obstacle in addressing quantum many-body systems is that the number of parameters required to represent a state of the system scales exponentially with the system size. However, one hardly needs the full Hilbert space but only a small subspace to describe physical properties of the system. From the point of view in quantum information theory, quantum entanglement is one of the main factors that is responsible for the complexity of studying quantum many-body systems. In particular, the amount of quantum entanglement in such system quantifies the correlation between constituents and therefore it determines how big the subspace is required to represent the system. In this thesis, we are interested in using tensor network states, a new emerging numerical technique to study dynamical and ground state properties of the system. The key ingredient is that a tensor network state which consists of a fairly small number of parameters can represent very well the complete wave function of a system. In one spatial dimension, Matrix Product State (MPS) has become a standard tool, combined with related algorithms such as Time Evolving Block Decimation algorithm (TEBD) to study dynamical properties of a locally-perturbed quantum lattice system. In this thesis, we introduce a technique called ``infinite boundary conditions'' that allows us to use a finite-size MPS to represent the state of an infinite-size system. As a result, one can simulate dynamics of a system for a longer time with a significant computational cost reduction compared to traditional methods. In two spatial dimensions, the tensor network known as Projected Entangled Pair States (PEPS) is emerging as a useful approach. However algorithms for optimizing the matrix elements of these tensors are complex and not very efficient. In particular, the efficiency of the optimization depends critically on the quality of the so-called ``environment'', which is expensive to compute. We have introduced a new method to accelerate the convergence of calculating the PEPS ground state by recycling the environment tensors for many optimization steps, vastly reducing the number of times that the full environment needs to be calculated.
We investigate finite-size scaling of genuine multisite entanglement in the ground state of quantum spin-1/2 Heisenberg ladders. We obtain the ground states of odd- and even-legged Heisenberg ladder Hamiltonians and compute genuine multisite entanglement, the generalized geometric measure (GGM), which shows that for even rungs, GGM increases for odd-legged ladder while it decreases for even ones. Interestingly, the ground state obtained by short-range dimer coverings, under the resonating valence bond ansatz, encapsulates the qualitative features of GGM for both the ladders. We find that while the quantity converges to a single value for higher legged odd- and even-ladders, in the asymptotic limit of a large number of rungs, the finite-size scaling exponents of the same tend to diverge. The scaling exponent of GGM is therefore capable to distinguish the odd–even dichotomy in Heisenberg ladders, even when the corresponding multisite entanglements merge. Introduction and background . The sequences of odd and even numbers mathematically reach the same infinity. The question is whether limits of functions of odd/even numbers reach the same limiting function as the respective sequences tend to infinity. In many-body physics, the dichotomy between physical properties of odd and even quantum ladders is known, and a typical result is exemplified by the seminal work of White, Noack, and Scalapino (PRL 73,886 (1994)), wherein the spin gap was shown to exhibit different behaviors in odd and even ladders. The values of the spin gaps were shown as converging to a single value with the increase in the number of legs, irrespective of whether we are with odd or even ladders. It may now be asked whether one can identify a physical quantity that would have different limiting values depending on whether odd or even legged ladders are followed to reach the infinite 2D square lattice. We answer this query in the affirmative by identifying the scaling of multisite entanglement as the physical quantity that is capable of doing the job. Main results . We evaluate the asymptotic genuine multisite entanglement of states of quantum spin models on ladder lattices with odd and even legs, and perform finite size scaling analysis of the same. We find that scaling exponents of the multisite entanglement for large lattices tend to diverge for odd and even ladders, as the number of legs are increased, even though the amount of the multisite entanglement converges in the asymptotic limit. Wider implications . We believe that the results presented in the manuscript will be of interest to the broad physics community working on stronglycorrelated systems, quantum information in many-body physics, and related quantum technologies. The methodology presented in the manuscript will, we believe, also appeal to the general theoretical physics community.
In this paper, we study global quantum discord (GQD) in infinite-size spin chains. For this purpose, in the framework of matrix product states (MPSs), we propose an effective procedure to calculate GQD (denoted as Gn) for consecutive n-site subchains in infinite chains. For a spin-1/2 three-body interaction model, whose ground state can be exactly expressed as MPSs, We use the procedure to study Gn with n up to 24. Then for a spin-1/2 XXZ chain, we firstly use infinite time-evolving block decimation (iTEBD) algorithm to obtain the approximate wavefunction in the from of MPSs, and then figure out Gn with n up to 18. In both models, Gn shows an interesting linear growth as the increase of n, that is, Gn = k*n+b. Moreover, in non-critical regions the slope k of Gn converges very fast, while in critical regions it converges relatively slow, and the behaviors are explained in a clear physical picture with the short-range and long-range correlations. Based on these results, we propose to use Gn/n to describe the global correlations in infinite chains. Gn/n has twofold physical meanings. Firstly, it can be regarded as "global discord per site", very similar to "energy per site" or "magnetization per site" in quantum magnetic systems. Secondly, Gn/n (when n is large enough) describes the quantum correlation between a single site and an (n-1)-site block. Then we successfully apply our theory to an exactly soluble infinite-size spin XY chain which is beyond the matrix product formula, and the Hamiltonian can reduce to the transverse-field Ising model and the XX model. The relation between GQD and quantum phase transitions in these models is discussed.
We study quantum phase transitions of a XXZ spin model with spin S = 1/2 and 1 in one dimension. The XXZ spin chain is one of basic models in understanding various one-dimensional magnetic materials. To study this model, we construct infinite-lattice matrix product state (iMPS), which is a tensor product form for a one-dimensional many-body quantum wave function. By using timeevolution- block-decimation method (TEBD) on iMPS, we obtain the ground states of the XXZ model at zero temperature. This method is very delicate in calculating ground states so that we developed a reliable method of finding the ground state with the dimension of entanglement coefficients up to 300, which is beyond the previous works. By analyzing ground-state energies, half-chain entanglement entropies, and entanglement spectrum, we found the signatures of quantum phase transitions between ferromagnetic phase, XY phase, Haldane phase, and antiferromagnetic phase.