Matrix Product Operators Research Articles

Characterization of topological phase transitions from a non-Abelian topological state and its Galois conjugate through condensation and confinement order parameters

Topological phases exhibit unconventional order that cannot be detected by any local order parameter. In the framework of Projected Entangled Pair States(PEPS), topological order is characterized by an entanglement symmetry of the local tensor which describes the model. This symmetry can take the form of a tensor product of group representations, or in the more general case a correlated symmetry action in the form of a Matrix Product Operator(MPO), which encompasses all string-net models. Among other things, these entanglement symmetries allow for the description of ground states and anyon excitations. Recently, the idea has been put forward to use those symmetries and the anyonic objects they describe as order parameters for probing topological phase transitions, and the applicability of this idea has been demonstrated for Abelian groups. In this paper, we extend this construction to the domain of non-Abelian models with MPO symmetries, and use it to study the breakdown of topological order in the double Fibonacci (DFib) string-net and its Galois conjugate, the non-hermitian double Yang-Lee (DYL) string-net. We start by showing how to construct topological order parameters for condensation and deconfinement of anyons using the MPO symmetries. Subsequently, we set up interpolations from the DFib and the DYL model to the trivial phase, and show that these can be mapped to certain restricted solid on solid(RSOS) models, which are equivalent to the $((5\pm\sqrt{5})/2)$-state Potts model, respectively. The known exact solutions of the statistical models allow us to locate the critical points, and to predict the critical exponents for the order parameters. We complement this by numerical study of the phase transitions, which fully confirms our theoretical predictions; remarkably, we find that both models exhibit a duality between the order parameters for condensation and deconfinement.

Computable Rényi mutual information: Area laws and correlations

The mutual information is a measure of classical and quantum correlations of great interest in quantum information. It is also relevant in quantum many-body physics, by virtue of satisfying an area law for thermal states and bounding all correlation functions. However, calculating it exactly or approximately is often challenging in practice. Here, we consider alternative definitions based on Rényi divergences. Their main advantage over their von Neumann counterpart is that they can be expressed as a variational problem whose cost function can be efficiently evaluated for families of states like matrix product operators while preserving all desirable properties of a measure of correlations. In particular, we show that they obey a thermal area law in great generality, and that they upper bound all correlation functions. We also investigate their behavior on certain tensor network states and on classical thermal distributions.

Operator front broadening in chaotic and integrable quantum chains

Operator spreading under unitary time evolution has attracted a lot of attention recently as a way to probe many-body quantum chaos. While quantities such as out-of-time-ordered correlators (OTOCs) do distinguish interacting from noninteracting systems, it has remained unclear to what extent they can truly diagnose chaotic vs integrable dynamics in many-body quantum systems. Here, we analyze operator spreading in generic one-dimensional many-body quantum systems using a combination of matrix product operator (MPO) and analytical techniques, focusing on the operator right weight. First, we show that while small bond dimension MPOs allow one to capture the exponentially decaying tail of the operator front, in agreement with earlier results, they lead to significant quantitative and qualitative errors for the actual front---defined by the maximum of the right weight. We find that while the operator front broadens diffusively in both integrable and chaotic interacting spin chains, the precise shape and scaling of the height of the front in integrable systems is anomalous for all accessible times. We interpret these results using a quasiparticle picture. This provides a sharp, although rather subtle, signature of many-body quantum chaos in the operator front.

Classical restrictions of generic matrix product states are quasi-locally Gibbsian

We show that norm squared amplitudes with respect to a local orthonormal basis (the classical restriction) of finite quantum systems on one-dimensional lattices can be exponentially well approximated by Gibbs states of local Hamiltonians (i.e., they are quasi-locally Gibbsian) if the classical conditional mutual information (CMI) of any connected tripartition of the lattice is rapidly decaying in the width of the middle region. For injective matrix product states, we, moreover, show that the classical CMI decays exponentially whenever the collection of matrix product operators satisfies a “purity condition,” a notion previously established in the theory of random matrix products. We, furthermore, show that violation of the purity condition enables a generalized notion of error correction on the virtual space, thus indicating the non-generic nature of such violations. We make this intuition more concrete by constructing a probabilistic model where purity is a typical property. The proof of our main result makes extensive use of the theory of random matrix products and may find applications elsewhere.

Classical simulation of lossy boson sampling using matrix product operators

Characterizing the computational advantage from noisy intermediate-scale quantum (NISQ) devices is an important task from theoretical and practical perspectives. Here, we numerically investigate the computational power of NISQ devices focusing on boson sampling, one of the well-known promising problems which can exhibit quantum supremacy. We study hardness of lossy boson sampling using matrix product operator (MPO) simulation to address the effect of photon-loss on classical simulability using MPO entanglement entropy (EE), which characterizes a running time of an MPO algorithm. An advantage of MPO simulation over other classical algorithms proposed to date is that its simulation accuracy can be efficiently controlled by increasing an MPO's bond dimension. Notably, we show by simulating lossy boson sampling using an MPO that as an input photon number grows, its computational cost, or MPO EE, behaves differently depending on a loss-scaling, exhibiting a different feature from that of lossless boson sampling. Especially when an output photon number scales faster than the square root of an input photon number, our study shows an exponential scaling of time complexity for MPO simulation. On the contrary, when an output photon number scales slower than the square root of an input photon number, MPO EE may decrease, indicating that an exponential time complexity might not be necessary.

Entanglement barriers in dual-unitary circuits

After quantum quenches in many-body systems, finite subsystems evolve nontrivially in time, eventually approaching a stationary state. In typical situations, the reduced density matrix of a given subsystem begins and ends this endeavor as a low-entangled vector in the space of operators. This means that if its entanglement in the operator space initially grows (which is generically the case), it must eventually decrease, describing a barrier-shaped curve. Understanding the shape of this ``entanglement barrier'' is interesting for three main reasons: (i) it quantifies the dynamics of entanglement in the (open) subsystem; (ii) it gives information on the approximability of the reduced density matrix by means of matrix product operators; and (iii) it shows qualitative differences depending on the type of dynamics undergone by the system, signaling quantum chaos. Here we compute exactly the shape of the entanglement barriers described by different R\'enyi entropies after quantum quenches in dual-unitary circuits initialized in a class of solvable matrix product states (MPS)s. We show that, for free (SWAP-like) circuits, the entanglement entropy behaves as in rational conformal field theories (CFT)s. On the other hand, for completely chaotic dual-unitary circuits, it behaves as in holographic CFTs, exhibiting a longer entanglement barrier that drops rapidly when the subsystem thermalizes. Interestingly, the entanglement spectrum is nontrivial in the completely chaotic case. Higher R\'enyi entropies behave in an increasingly similar way to rational CFTs, such that the free and completely chaotic barriers are identical in the limit of infinite replicas (i.e., for the so called min-entropy). We also show that, upon increasing the bond dimension of the MPSs, the barrier maintains the same shape. It simply shifts to the left to accommodate for the larger initial entanglement.

Quantum Heat Statistics with Time-Evolving Matrix Product Operators

We present a numerically exact method to compute the full counting statistics of heat transfer in non-Markovian open quantum systems, which is based on the time-evolving matrix product operator (TEMPO) algorithm. This approach is applied to the paradigmatic spin-boson model in order to calculate the mean and fluctuations of the heat transferred to the environment during thermal equilibration. We show that system-reservoir correlations make a significant contribution to the heat statistics at low temperature and present a variational theory that quantitatively explains our numerical results. We also demonstrate a fluctuation-dissipation relation connecting the mean and variance of the heat distribution at high temperature. Our results reveal that system-bath interactions make a significant contribution to heat transfer even when the dynamics of the open system is effectively Markovian. The method presented here provides a flexible and general tool to predict the fluctuations of heat transfer in open quantum systems in non-perturbative regimes.

Anomalous quantum information scrambling for Z3 parafermion chains

Parafermions are exotic quasiparticles with non-Abelian fractional statistics that could be exploited to realize universal topological quantum computing. Here, we study the scrambling of quantum information in one-dimensional parafermionic chains, with a focus on $\mathbb{Z}_3$ parafermions in particular. We use the generalized out-of-time-ordered correlators (OTOCs) as a measure of the information scrambling and introduce an efficient method based on matrix product operators to compute them. With this method, we compute the OTOCs for $\mathbb{Z}_3$ parafermions chains up to $200$ sites for the entire early growth region. We find that, in stark contrast to the dynamics of conventional fermions or bosons, the information scrambling light cones for parafermions can be both symmetric and asymmetric, even for inversion-invariant Hamiltonians involving only hopping terms. In addition, we find a deformed light cone structure with a sharp peak at the boundary of the parafermion chains in the topological regime, which gives a unambiguous evidence of the strong zero modes at infinite temperature.

Efficient Exploration of Hamiltonian Parameter Space for Optimal Control of Non-Markovian Open Quantum Systems.

We present a general method to efficiently design optimal control sequences for non-Markovian open quantum systems, and illustrate it by optimizing the shape of a laser pulse to prepare a quantum dot in a specific state. The optimization of control procedures for quantum systems with strong coupling to structured environments-where time-local descriptions fail-is a computationally challenging task. Wemodify the numerically exact time evolving matrix product operator (TEMPO) method, such that it allows the repeated computation of the time evolution of the reduced system density matrix for various sets of control parameters at very low computational cost. This method is potentially useful for studying numerous optimal control problems, in particular in solid state quantum devices where the coupling to vibrational modes is typically strong.

Multipartite nonlocality in one-dimensional quantum spin chains at finite temperatures

Multipartite nonlocality is an important measure of multipartite quantum correlations. In this paper, we show that the nonlocal $n$-site Mermin-Klyshko operator ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{M}}_{n}$ can be exactly expressed as a matrix product operator with a bond dimension $D=2$, and then the calculation of nonlocality measure $\mathcal{S}$ can be simplified into standard one-dimensional (1D) tensor networks. With the help of this technique, we analyze finite-temperature multipartite nonlocality in several typical 1D spin chains, including an $XX$ model, an $XXZ$ model, and a Kitaev model. For the $XX$ model and the $XXZ$ model, in a finite-temperature region, the logarithm measure of nonlocality (${log}_{2}\mathcal{S}$) is a linear function of the temperature $T$, i.e., ${log}_{2}\mathcal{S}\ensuremath{\sim}\ensuremath{-}aT+b$. It provides us with an intuitive picture about how thermodynamic fluctuations destroy multipartite nonlocality in 1D quantum chains. Moreover, in the $XX$ model $\mathcal{S}$ presents a magnetic-field-induced oscillation at low temperatures. This behavior has a nonlocal nature and cannot be captured by local properties such as the magnetization. Finally, for the Kitaev model, we find that in the limit $T\ensuremath{\rightarrow}0$ and $N\ensuremath{\rightarrow}\ensuremath{\infty}$ the nonlocality measure may be used as an alternative order parameter for the topological-type quantum phase transition in the model.

Tensor network approach to thermalization in open quantum many-body systems.

We investigate the relaxation dynamics of open nonintegrable quantum many-body systems in the thermodynamic limit by using a tensor-network formalism. We simulate the Lindblad quantum master equation (LQME) of infinite systems by making use of the uniform matrix product operators (MPO) as the ansatz of their density matrices. Furthermore, we establish a method to measure the thermodynamic equivalence between two states described by the uniform MPOs. We numerically show that when an initial state of the LQME is a thermal Gibbs state, a time evolved state is always indistinguishable from a Gibbs state with a time-dependent effective temperature in the weak-dissipation and thermodynamic limit.

Simulating noisy quantum circuits with matrix product density operators

Simulating quantum circuits with classical computers requires resources growing exponentially in terms of system size. Real quantum computer with noise, however, may be simulated polynomially with various methods considering different noise models. In this work, we simulate random quantum circuits in 1D with Matrix Product Density Operators (MPDO), for different noise models such as dephasing, depolarizing, and amplitude damping. We show that the method based on Matrix Product States (MPS) fails to approximate the noisy output quantum states for any of the noise models considered, while the MPDO method approximates them well. Compared with the method of Matrix Product Operators (MPO), the MPDO method reflects a clear physical picture of noise (with inner indices taking care of the noise simulation) and quantum entanglement (with bond indices taking care of two-qubit gate simulation). Consequently, in case of weak system noise, the resource cost of MPDO will be significantly less than that of the MPO due to a relatively small inner dimension needed for the simulation. In case of strong system noise, a relatively small bond dimension may be sufficient to simulate the noisy circuits, indicating a regime that the noise is large enough for an `easy' classical simulation. Moreover, we propose a more effective tensor updates scheme with optimal truncations for both the inner and the bond dimensions, performed after each layer of the circuit, which enjoys a canonical form of the MPDO for improving simulation accuracy. With truncated inner dimension to a maximum value $\kappa$ and bond dimension to a maximum value $\chi$, the cost of our simulation scales as $\sim ND\kappa^3\chi^3$, for an $N$-qubit circuit with depth $D$.

Two-dimensional topological order and operator algebras

We review recent interactions between mathematical theory of two-dimensional topological order and operator algebras, particularly the Jones theory of subfactors. The role of representation theory in terms of tensor categories is emphasized. Connections to two-dimensional conformal field theory are also presented. In particular, we discuss anyon condensation, gapped domain walls and matrix product operators in terms of operator algebras.

Temperature-dependent energy diffusion in chaotic spin chains

We study the temperature dependence of energy diffusion in two chaotic gapped quantum spin chains, a tilted-field Ising model and an XZ model, using an open system approach. We introduce an energy imbalance by coupling the chain to thermal baths at its boundary and study the non-equilibrium steady states of the resulting Lindblad dynamics using a matrix product operator ansatz for the density matrix. We define an effective local temperature profile by comparing local reduced density matrices in the steady state to those of a uniform thermal state. We then measure the energy current for a variety of driving temperatures and extract the temperature dependence of the energy diffusion constant. For the Ising model, we are able to study temperatures well below the energy gap and find a regime of dilute excitations where rare three-body collisions control energy diffusion. A kinetic model correctly predicts the observed exponential increase of the energy diffusion constant at low temperatures. For the XZ model, we are only able to access intermediate to high temperatures relative to the energy gap and we show that the data is well described by an expansion around the infinite temperature limit. We also discuss the limitations of the particular driving scheme and suggest that lower temperatures can be accessed using larger baths.

Improved Thermal Area Law and Quasilinear Time Algorithm for Quantum Gibbs States

One of the most fundamental problems in quantum many-body physics is the characterization of correlations among thermal states. Of particular relevance is the thermal area law, which justifies the tensor network approximations to thermal states with a bond dimension growing polynomially with the system size. In the regime of sufficiently low temperatures, which is particularly important for practical applications, the existing techniques do not yield optimal bounds. Here, we propose a new thermal area law that holds for generic many-body systems on lattices. We improve the temperature dependence from the original $\mathcal{O}(\beta)$ to $\tilde{\mathcal{O}}(\beta^{2/3})$, thereby suggesting diffusive propagation of entanglement by imaginary time evolution. This qualitatively differs from the real-time evolution which usually induces linear growth of entanglement. We also prove analogous bounds for the R\'enyi entanglement of purification and the entanglement of formation. Our analysis is based on a polynomial approximation to the exponential function which provides a relationship between the imaginary-time evolution and random walks. Moreover, for one-dimensional (1D) systems with $n$ spins, we prove that the Gibbs state is well-approximated by a matrix product operator with a sublinear bond dimension of $e^{\sqrt{\tilde{\mathcal{O}}(\beta \log(n))}}$. This proof allows us to rigorously establish, for the first time, a quasi-linear time classical algorithm for constructing an MPS representation of 1D quantum Gibbs states at arbitrary temperatures of $\beta = o(\log(n))$. Our new technical ingredient is a block decomposition of the Gibbs state, that bears resemblance to the decomposition of real-time evolution given by Haah et al., FOCS'18.

Matrix product operator symmetries and intertwiners in string-nets with domain walls

We provide a description of virtual non-local matrix product operator (MPO) symmetries in projected entangled pair state (PEPS) representations of string-net models. Given such a PEPS representation, we show that the consistency conditions of its MPO symmetries amount to a set of six coupled equations that can be identified with the pentagon equations of a bimodule category. This allows us to classify all equivalent PEPS representations and build MPO intertwiners between them, synthesising and generalising the wide variety of tensor network representations of topological phases. Furthermore, we use this generalisation to build explicit PEPS realisations of domain walls between different topological phases as constructed by Kitaev and Kong [Commun. Math. Phys. 313 (2012) 351-373]. While the prevailing abstract categorical approach is sufficient to describe the structure of topological phases, explicit tensor network representations are required to simulate these systems on a computer, such as needed for calculating thresholds of quantum error-correcting codes based on string-nets with boundaries. Finally, we show that all these string-net PEPS representations can be understood as specific instances of Turaev-Viro state-sum models of topological field theory on three-manifolds with a physical boundary, thereby putting these tensor network constructions on a mathematically rigorous footing.

Hydrodynamics of nonintegrable systems from a relaxation-time approximation

We develop a general kinetic theory framework to describe the hydrodynamics of strongly interacting, nonequilibrium quantum systems in which integrability is weakly broken, leaving a few residual conserved quantities. This framework is based on a generalized relaxation-time approximation; it gives a simple, but surprisingly accurate, prescription for computing nonequilibrium transport even in strongly interacting systems. This approximation reproduces the crossover from generalized to conventional hydrodynamics in interacting one-dimensional Bose gases with integrability-breaking perturbations, both with and without momentum conservation. It also predicts the hydrodynamics of chaotic quantum spin chains, in good agreement with matrix product operator calculations.

Entanglement entropy of two-species hard-core bosons in one dimension

We study the von Neumann entropy of a model for two-species hard-core bosons in one dimension. In this model, the same-species bosons satisfy hard-core conditions, while the different-species bosons are allowed to occupy the same site with a local interaction U . At half-filling, by Jordan-Wigner transformation, the model is exactly mapped to a fermionic Hubbard model. The phase transition from superfluid U = 0 to Mott insulator U > 0 can be explained by simple one-band theory at half-filling. We measure the von Neumann entanglement entropy of the ground states for the half-filled case and away from half-filling to understand quantum phase transitions. To achieve this goal, we use a time-evolution-block decimation method with infinite-size matrix product state and also use a density matrix renormalization group with matrix product operators with large bond dimensions up to 300. We found strong evidence that the local minimum point of the von Neumann entanglement entropy is the quantum critical point for finite-bond-dimension matrix product states.

Tangent-space methods for truncating uniform MPS

An essential primitive in quantum tensor network simulations is the problem of approximating a matrix product state with one of a smaller bond dimension. This problem forms the central bottleneck in algorithms for time evolution and for contracting projected entangled pair states. We formulate a tangent-space based variational algorithm to achieve this goal for uniform (infinite) matrix product states. The algorithm exhibits a favourable scaling of the computational cost, and we demonstrate its usefulness by several examples involving the multiplication of a matrix product state with a matrix product operator.

Phase transitions in the one-dimensional ionic Hubbard model

We study quantum phase transitions by measuring the bond energy, the number density, and the half-chain entanglement entropy in the one-dimensional ionic Hubbard model. By using the matrix product operator to perform the infinite density matrix renormalization group, we obtain the ground states in the canonical form of matrix product states. Depending on the chemical potential and the staggered potential, the number density and the half-chain entanglement entropy shows clear signatures of a Mott transition. Our results confirm the success of using the matrix product operator method to investigate itinerant fermion systems.