By combining the Grassmann algebra with multi-scale entanglement renormalization ansatz (MERA), we introduce a new unbiased and effective numerical method for simulating 2D strongly correlated electronic systems. The new GMERA method inherits all the advantages of MERA, which constructs the variational wave function based on complicated tensor network. Besides it can deal with fermionic properties of the system due to Grassmann algebra through local tensor contractions. This general method can treat different tensor network structures in a universal way. We show several benchmark calculations of the GMERA method, including the free fermion model, tight binding model, as well as the t-J model with hole doping.
By combining the Grassmann algebra with multi-scale entanglement renormalization ansatz (MERA), we introduce a new unbiased and effective numerical method for simulating 2D strongly correlated electronic systems. The new GMERA method inherits all the advantages of MERA, which constructs the variational wave function based on complicated tensor network. Besides it can deal with fermionic properties of the system due to Grassmann algebra through local tensor contractions. This general method can treat different tensor network structures in a universal way. We show several benchmark calculations of the GMERA method, including the free fermion model, tight binding model, as well as the t-J model with hole doping.
In this thesis, we use the new numerical method multiscale entanglement
renormalization ansatz (MERA) to study the one dimensional spin system. This
method, MERA, utilizes a network of tensors to represent quantum many-body
systems on a lattice. Other examples of tensor network state are matrix product
states (MPSs) for 1D systems, tree tensor networks (TTNs) for systems with a
tree shape, and projected entangled-pair states (PEPSs) for 2D systems and be-
yond. The three structures dier in the graph that denes how the tensors are
interconnected into a network: The graphs for MPSs, TTNs, and 2D PEPSs are,
respectively, a chain, a tree, and a 2D lattice. Importantly, from these tensor
networks the expectation value of local observables can be computed eciently.
Developing algorithms to simulate quantum many-body systems is important
for understanding and studying the physical principles. And these may become
powerful tools for the analysis of quantum computation and quantum information
theory. But in general ecient simulation is dicult, due to the fundamental
quantum physical laws.
The main method we have studied and worked on is multiscale entangle-
ment renormalization ansatz (MERA) which is a variational ansatz for many-
body states. MERA is a tree tensor network for 1D, 2D systems and beyond.
We have developed code of ternary-MERA to compute the ground state energy,
expected values of local observables, and correlators for 1D lattice systems. In
a nutshell, MERA is an algorithm for entanglement renormalization which is a
numerical technique based on locally reorganizing the Hilbert space of a quantum
many-body system with the aim to reduce the amount of entanglement in its wave
function. We use this method to study 1D Ising model, XXZ model, and J1 - J2
model, and then comparing the numerical results with exact diagonalization (ED)
and density matrix renormalization group (DMRG).
Strongly correlated quantum matter, such as the insufficiently understood high- temperature superconductors, are difficult to simulate numerically. The multiscale entanglement renormalization ansatz (MERA) is a hierarchical class of tensor networks that efficiently captures the typical ground-state entanglement structures. Motivated by recent advances to employ MERA in quantum-classical algorithms, the authors take a fresh look here, analyzing new variations of the approach including tensor Trotterization and importance sampling. Optimal contraction sequences and algorithmic phase diagrams show that these quantum-inspired methods can substantially reduce computational costs.
The projective construction (the slave-particle approach) has played an very important role in understanding strongly correlated systems, such as the emergence of fermions, anyons, and gauge theory in quantum spin liquids and quantum Hall states. Recently, fermionic Projected Entangled Pair States (fPEPS) have been introduced to effciently represent many-body fermionic states. In this paper, we show that the strongly correlated bosonic/fermionic states obtained both from the projective construction and the fPEPS approach can be represented systematically as Grassmann tensor product states. This construction can also be applied to all other tensor network states approaches. The Grassmann tensor product states allow us to encode many-body bosonic/fermionic states effciently with a polynomial number of parameters. We also generalize the tensor-entanglement renormalization group (TERG) method for complex tensor networks to Grassmann tensor networks. This allows us to approximate the norm and average local operators of Grassmann tensor product states in polynomial time, and hence leads to a variational approach for describing strongly correlated bosonic/fermionic systems in higher dimensions.
Tensor network quantum states are powerful tools for strongly correlated systems, tailored to capture local correlations such as in ground states with entanglement area laws. When applying tensor network states to interacting fermionic systems, a proper choice of the basis or orbitals can reduce the bond dimension of tensors and provide physically relevant orbitals. We introduce such a change of basis with unitary gates obtained from compressing fermionic Gaussian states into quantum circuits corresponding to various tensor networks. These circuits can reduce the ground-state entanglement entropy and improve the performance of algorithms such as the density matrix renormalization group. We study the Anderson impurity model with one and two impurities to show the potential of the method for improving computational efficiency and interpreting impurity physics. Furthermore, fermionic Gaussian circuits can also suppress entanglement during the time evolution out of low-energy state. Last, we consider Gaussian multiscale entanglement renormalization ansatz (GMERA) circuits which compress fermionic Gaussian states hierarchically. The emergent coarse-grained physical models from these GMERA circuits are studied in terms of their entanglement properties and suitability for performing time evolution.
In a recent contribution [P. Corboz, G. Evenbly, F. Verstraete, and G. Vidal, arXiv:0904.4151 (unpublished)] entanglement renormalization was generalized to fermionic lattice systems in two spatial dimensions. Entanglement renormalization is a real-space coarse-graining transformation for lattice systems that produces a variational ansatz, the multiscale entanglement renormalization ansatz (MERA), for the ground states of local Hamiltonians. In this paper we describe in detail the fermionic version of the MERA formalism and algorithm. Starting from the bosonic MERA, which can be regarded both as a quantum circuit or in relation to a coarse-graining transformation, we indicate how the scheme needs to be modified to simulate fermions. To confirm the validity of the approach, we present benchmark results for free and interacting fermions on a square lattice with sizes between $6\ifmmode\times\else\texttimes\fi{}6$ and $162\ifmmode\times\else\texttimes\fi{}162$ and with periodic boundary conditions. The present formulation of the approach applies to generic tensor network algorithms.
Abstract The strongly correlated fermions play a vital role in modern physics. For a given fermionic Hamiltonian system, the most widely used approach to exploring the underlying physics is to study the wave function that incorporates Fermi–Dirac statistics, which can be obtained variationally by energy minimization or by imaginary-time evolution. In this work, we develop an accurate tensor network method for one-dimensional interacting fermionic models based on the coherent-state path-integral representation of the fermionic partition function. Employing the coherent-state representation, the partition function is effectively represented as a (1+1)-dimensional anisotropic Grassmann-valued tensor network, and the Grassmann version of the corner transfer-matrix renormalization group algorithm is developed to contract the tensor network and evaluate physical quantities. We validate our method on the one-dimensional fermionic Hubbard model with a magnetic field, where the essential features of the phase diagram in the ( μ , B ) plane are quantitatively captured. Our work offers a promising approach to interacting fermionic models within the framework of tensor networks.