No Arabic abstract
We derive an exact matrix product state representation of the Haldane-Rezayi state on both the cylinder and torus geometry. Our derivation is based on the description of the Haldane-Rezayi state as a correlator in a non-unitary logarithmic conformal field theory. This construction faithfully captures the ten degenerate ground states of this model state on the torus. Using the cylinder geometry, we probe the gapless nature of the phase by extracting the correlation length, which diverges in the thermodynamic limit. The numerically extracted topological entanglement entropies seem to only probe the Abelian part of the theory, which is reminiscent of the Gaffnian state, another model state deriving from a non-unitary conformal field theory.
We study the thin torus limit of the Haldane-Rezayi state. Eight of the ten ground states are found to assume a simple product form in this limit, as is known to be the case for many other quantum Hall trial wave functions. The two remaining states have a somewhat unusual thin torus limit, where a broken pair of defects forming a singlet is completely delocalized. We derive these limits from the wave functions on the cylinder, and deduce the dominant matrix elements of the thin torus hollow-core Hamiltonians. We find that there are gapless excitations in the thin torus limit. This is in agreement with the expectation that local Hamiltonians stabilizing wave functions associated with non-unitary conformal field theories are gapless. We also use the thin torus analysis to obtain explicit counting formulas for the zero modes of the hollow-core Hamiltonian on the torus, as well as for the parent Hamiltonians of several other paired and related quantum Hall states.
The density-matrix renormalization group method has become a standard computational approach to the low-energy physics as well as dynamics of low-dimensional quantum systems. In this paper, we present a new set of applications, available as part of the ALPS package, that provide an efficient and flexible implementation of these methods based on a matrix-product state (MPS) representation. Our applications implement, within the same framework, algorithms to variationally find the ground state and low-lying excited states as well as simulate the time evolution of arbitrary one-dimensional and two-dimensional models. Implementing the conservation of quantum numbers for generic Abelian symmetries, we achieve performance competitive with the best codes in the community. Example results are provided for (i) a model of itinerant fermions in one dimension and (ii) a model of quantum magnetism.
We prove that ground states of gapped local Hamiltonians on an infinite spin chain can be efficiently approximated by matrix product states with a bond dimension which scales as D~(L-1)/epsilon, where any local quantity on L consecutive spins is approximated to accuracy epsilon.
While general quantum many-body systems require exponential resources to be simulated on a classical computer, systems of non-interacting fermions can be simulated exactly using polynomially scaling resources. Such systems may be of interest in their own right, but also occur as effective models in numerical methods for interacting systems, such as Hartree-Fock, density functional theory, and many others. Often it is desirable to solve systems of many thousand constituent particles, rendering these simulations computationally costly despite their polynomial scaling. We demonstrate how this scaling can be improved by adapting methods based on matrix product states, which have been enormously successful for low-dimensional interacting quantum systems, to the case of free fermions. Compared to the case of interacting systems, our methods achieve an exponential speedup in the entanglement entropy of the state. We demonstrate their use to solve systems of up to one million sites with an effective MPS bond dimension of 10^15.
We present a unified framework for renormalization group methods, including Wilsons numerical renormalization group (NRG) and Whites density-matrix renormalization group (DMRG), within the language of matrix product states. This allows improvements over Wilsons NRG for quantum impurity models, as we illustrate for the one-channel Kondo model. Moreover, we use a variational method for evaluating Greens functions. The proposed method is more flexible in its description of spectral properties at finite frequencies, opening the way to time-dependent, out-of-equilibrium impurity problems. It also substantially improves computational efficiency for one-channel impurity problems, suggesting potentially emph{linear} scaling of complexity for $n$-channel problems.