-
Variational adiabatic transport of tensor networks
Authors:
Hyeongjin Kim,
Matthew T. Fishman,
Dries Sels
Abstract:
We discuss a tensor network method for constructing the adiabatic gauge potential -- the generator of adiabatic transformations -- as a matrix product operator, which allows us to adiabatically transport matrix product states. Adiabatic evolution of tensor networks offers a wide range of applications, of which two are explored in this paper: improving tensor network optimization and scanning phase…
▽ More
We discuss a tensor network method for constructing the adiabatic gauge potential -- the generator of adiabatic transformations -- as a matrix product operator, which allows us to adiabatically transport matrix product states. Adiabatic evolution of tensor networks offers a wide range of applications, of which two are explored in this paper: improving tensor network optimization and scanning phase diagrams. By efficiently transporting eigenstates to quantum criticality and performing intermediary density matrix renormalization group (DMRG) optimizations along the way, we demonstrate that we can compute ground and low-lying excited states faster and more reliably than a standard DMRG method at or near quantum criticality. We demonstrate a simple automated step size adjustment and detection of the critical point based on the norm of the adiabatic gauge potential. Remarkably, we are able to reliably transport states through the critical point of the models we study.
△ Less
Submitted 28 November, 2023; v1 submitted 1 November, 2023;
originally announced November 2023.
-
Fragmented superconductivity in the Hubbard model as solitons in Ginzburg-Landau theory
Authors:
Niccolò Baldelli,
Benedikt Kloss,
Matthew Fishman,
Alexander Wietek
Abstract:
The phenomena of superconductivity and charge density waves are observed in close vicinity in many strongly correlated materials. Increasing evidence from experiments and numerical simulations suggests both phenomena can also occur in an intertwined manner, where the superconducting order parameter is coupled to the electronic density. Employing density matrix renormalization group simulations, we…
▽ More
The phenomena of superconductivity and charge density waves are observed in close vicinity in many strongly correlated materials. Increasing evidence from experiments and numerical simulations suggests both phenomena can also occur in an intertwined manner, where the superconducting order parameter is coupled to the electronic density. Employing density matrix renormalization group simulations, we investigate the nature of such an intertwined state of matter stabilized in the phase diagram of the elementary $t$-$t^\prime$-$U$ Hubbard model in the strong coupling regime. Remarkably, the condensate of Cooper pairs is shown to be fragmented in the presence of a charge density wave where more than one pairing wave function is macroscopically occupied. Moreover, we provide conclusive evidence that the macroscopic wave functions of the superconducting fragments are well-described by soliton solutions of a Ginzburg-Landau equation in a periodic potential constituted by the charge density wave. In the presence of an orbital magnetic field, the order parameters are gauge invariant, and superconducting vortices are pinned between the stripes. This intertwined Ginzburg-Landau theory is proposed as an effective low-energy description of the stripe fragmented superconductor.
△ Less
Submitted 21 July, 2023;
originally announced July 2023.
-
Equilibrium Quantum Impurity Problems via Matrix Product State Encoding of the Retarded Action
Authors:
Benedikt Kloss,
Julian Thoenniss,
Michael Sonner,
Alessio Lerose,
Matthew T. Fishman,
E. M. Stoudenmire,
Olivier Parcollet,
Antoine Georges,
Dmitry A. Abanin
Abstract:
In the $0+1$ dimensional imaginary-time path integral formulation of quantum impurity problems, the retarded action encodes the hybridization of the impurity with the bath. In this Article, we explore the computational power of representing the retarded action as matrix product state (RAMPS). We focus on the challenging Kondo regime of the single-impurity Anderson model, where non-perturbative str…
▽ More
In the $0+1$ dimensional imaginary-time path integral formulation of quantum impurity problems, the retarded action encodes the hybridization of the impurity with the bath. In this Article, we explore the computational power of representing the retarded action as matrix product state (RAMPS). We focus on the challenging Kondo regime of the single-impurity Anderson model, where non-perturbative strong-correlation effects arise at very low energy scales. We demonstrate that the RAMPS approach reliably reaches the Kondo regime for a range of interaction strengths $U$, with a numerical error scaling as a weak power law with inverse temperature. We investigate the convergence behavior of the method with respect to bond dimension and time discretization by analyzing the error of local observables in the full interacting problem and find polynomial scaling in both parameters. Our results show that the RAMPS approach offers promise as an alternative tool for studying quantum impurity problems in regimes that challenge established methods, such as multi-orbital systems. Overall, our study contributes to the development of efficient and accurate non-wavefunction-based tensor-network methods for quantum impurity problems.
△ Less
Submitted 29 June, 2023;
originally announced June 2023.
-
Disentangling Interacting Systems with Fermionic Gaussian Circuits: Application to the Single Impurity Anderson Model
Authors:
Ang-Kun Wu,
Matthew T. Fishman,
J. H. Pixley,
E. M. Stoudenmire
Abstract:
Tensor network representations of quantum many-body states provide powerful tools for strongly correlated systems, tailored to capture local correlations such as ground states exhibiting entanglement area laws. When applying tensor network states to interacting fermionic systems, a proper choice of basis or orbitals can greatly reduce the bond dimension of tensors and speed up calculations. We int…
▽ More
Tensor network representations of quantum many-body states provide powerful tools for strongly correlated systems, tailored to capture local correlations such as ground states exhibiting entanglement area laws. When applying tensor network states to interacting fermionic systems, a proper choice of basis or orbitals can greatly reduce the bond dimension of tensors and speed up calculations. We introduce such a change of basis with unitary gates obtained via 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 1D single impurity Anderson model to show the power of the method in improving computational efficiency and interpreting impurity physics. Furthermore, fermionic Gaussian circuits also show potential for suppressing entanglement during the time evolution of a low-lying excited state that is used to compute the impurity Green's function. Lastly, we consider Gaussian multi-scale 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.
△ Less
Submitted 19 December, 2022;
originally announced December 2022.
-
Extracting Off-Diagonal Order from Diagonal Basis Measurements
Authors:
Bo Xiao,
Javier Robledo Moreno,
Matthew Fishman,
Dries Sels,
Ehsan Khatami,
Richard Scalettar
Abstract:
Quantum gas microscopy has developed into a powerful tool to explore strongly correlated quantum systems. However, discerning phases with topological or off-diagonal long range order requires the ability to extract these correlations from site-resolved measurements. Here, we show that a multi-scale complexity measure can pinpoint the transition to and from the bond ordered wave phase of the one-di…
▽ More
Quantum gas microscopy has developed into a powerful tool to explore strongly correlated quantum systems. However, discerning phases with topological or off-diagonal long range order requires the ability to extract these correlations from site-resolved measurements. Here, we show that a multi-scale complexity measure can pinpoint the transition to and from the bond ordered wave phase of the one-dimensional extended Hubbard model with an off-diagonal order parameter, sandwiched between diagonal charge and spin density wave phases, using only diagonal descriptors. We study the model directly in the thermodynamic limit using the recently developed variational uniform matrix product states algorithm, and draw our samples from degenerate ground states related by global spin rotations, emulating the projective measurements that are accessible in experiments. Our results will have important implications for the study of exotic phases using optical lattice experiments.
△ Less
Submitted 21 September, 2022;
originally announced September 2022.
-
Field-Tunable Topological Phase Transitions and Spin-Hall Effects in 2D Crystals
Authors:
Maxwell Fishman,
Debdeep Jena
Abstract:
As recent additions to the catalog of 2D crystals, silicene and other silicene-class crystals have numerous unique properties currently being investigated and considered for use in novel device applications. In this paper, we investigate electronic and transport properties of silicene in a field effect transistor geometry. We find that the Berry curvature of silicene-class crystals can be continuo…
▽ More
As recent additions to the catalog of 2D crystals, silicene and other silicene-class crystals have numerous unique properties currently being investigated and considered for use in novel device applications. In this paper, we investigate electronic and transport properties of silicene in a field effect transistor geometry. We find that the Berry curvature of silicene-class crystals can be continuously tuned by a perpendicular electric field. By direct calculation of the $Z_2$ invariant, we confirm that an electronic phase transition from a topological insulator to a band insulator occurs when the electric field passes a critical value. In a device setting with asymmetric gate voltages, this field-tunable Berry curvature generates a large spin current transverse to the charge current. When the electric field strength surpasses the critical value, the bulk spin current is found to change direction and greatly decrease in magnitude. This finding of a large magnitude, switchable spin current suggests that the silicene family of 2D crystals could be an attractive candidate for field-tunable charge-spin conversion. Such field-tunable phase transitions between topologically distinct phases could be useful for robust qubits as well.
△ Less
Submitted 17 December, 2020;
originally announced December 2020.
-
The Hubbard model on the Bethe lattice via variational uniform tree states: metal-insulator transition and a Fermi liquid
Authors:
Peter Lunts,
Antoine Georges,
E. Miles Stoudenmire,
Matthew Fishman
Abstract:
We numerically solve the Hubbard model on the Bethe lattice with finite coordination number $z=3$, and determine its zero-temperature phase diagram. For this purpose, we introduce and develop the `variational uniform tree state' (VUTS) algorithm, a tensor network algorithm which generalizes the variational uniform matrix product state algorithm to tree tensor networks. Our results reveal an antife…
▽ More
We numerically solve the Hubbard model on the Bethe lattice with finite coordination number $z=3$, and determine its zero-temperature phase diagram. For this purpose, we introduce and develop the `variational uniform tree state' (VUTS) algorithm, a tensor network algorithm which generalizes the variational uniform matrix product state algorithm to tree tensor networks. Our results reveal an antiferromagnetic insulating phase and a paramagnetic metallic phase, separated by a first-order doping-driven metal-insulator transition. We show that the metallic state is a Fermi liquid with coherent quasiparticle excitations for all values of the interaction strength $U$, and we obtain the finite quasiparticle weight $Z$ from the single-particle occupation function of a generalized "momentum" variable. We find that $Z$ decreases with increasing $U$, ultimately saturating to a non-zero, doping-dependent value. Our work demonstrates that tensor-network calculations on tree lattices, and the VUTS algorithm in particular, are a platform for obtaining controlled results for phenomena absent in one dimension, such as Fermi liquids, while avoiding computational difficulties associated with tensor networks in two dimensions. We envision that future studies could observe non-Fermi liquids, interaction-driven metal-insulator transitions, and doped spin liquids using this platform.
△ Less
Submitted 2 November, 2020; v1 submitted 13 October, 2020;
originally announced October 2020.
-
The ITensor Software Library for Tensor Network Calculations
Authors:
Matthew Fishman,
Steven R. White,
E. Miles Stoudenmire
Abstract:
ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagram notation, which allows users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of tensor network algorithms. After discussing the philosophy behind the ITensor…
▽ More
ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagram notation, which allows users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of tensor network algorithms. After discussing the philosophy behind the ITensor approach, we show examples of each part of the interface including Index objects, the ITensor product operator, tensor factorizations, tensor storage types, algorithms for matrix product state (MPS) and matrix product operator (MPO) tensor networks, quantum number conserving block-sparse tensors, and the NDTensors library. We also review publications that have used ITensor for quantum many-body physics and for other areas where tensor networks are increasingly applied. To conclude we discuss promising features and optimizations to be added in the future.
△ Less
Submitted 20 December, 2021; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Wavefunction positivization via automatic differentiation
Authors:
Giacomo Torlai,
Juan Carrasquilla,
Matthew T. Fishman,
Roger G. Melko,
Matthew P. A. Fisher
Abstract:
We introduce a procedure to systematically search for a local unitary transformation that maps a wavefunction with a non-trivial sign structure into a positive-real form. The transformation is parametrized as a quantum circuit compiled into a set of one and two qubit gates. We design a cost function that maximizes the average sign of the output state and removes its complex phases. The optimizatio…
▽ More
We introduce a procedure to systematically search for a local unitary transformation that maps a wavefunction with a non-trivial sign structure into a positive-real form. The transformation is parametrized as a quantum circuit compiled into a set of one and two qubit gates. We design a cost function that maximizes the average sign of the output state and removes its complex phases. The optimization of the gates is performed through automatic differentiation algorithms, widely used in the machine learning community. We provide numerical evidence for significant improvements in the average sign for a two-leg triangular Heisenberg ladder with next-to-nearest neighbour and ring-exchange interactions. This model exhibits phases where the sign structure can be removed by simple local one-qubit unitaries, but also an exotic Bose-metal phase whose sign structure induces "Bose surfaces" with a fermionic character and a higher entanglement that requires deeper circuits.
△ Less
Submitted 10 February, 2020; v1 submitted 11 June, 2019;
originally announced June 2019.
-
Faster Methods for Contracting Infinite 2D Tensor Networks
Authors:
M. T. Fishman,
L. Vanderstraeten,
V. Zauner-Stauber,
J. Haegeman,
F. Verstraete
Abstract:
We revisit the corner transfer matrix renormalization group (CTMRG) method of Nishino and Okunishi for contracting two-dimensional (2D) tensor networks and demonstrate that its performance can be substantially improved by determining the tensors using an eigenvalue solver as opposed to the power method used in CTMRG. We also generalize the variational uniform matrix product state (VUMPS) ansatz fo…
▽ More
We revisit the corner transfer matrix renormalization group (CTMRG) method of Nishino and Okunishi for contracting two-dimensional (2D) tensor networks and demonstrate that its performance can be substantially improved by determining the tensors using an eigenvalue solver as opposed to the power method used in CTMRG. We also generalize the variational uniform matrix product state (VUMPS) ansatz for diagonalizing 1D quantum Hamiltonians to the case of 2D transfer matrices and discuss similarities with the corner methods. These two new algorithms will be crucial to improving the performance of variational infinite projected entangled pair state (PEPS) methods.
△ Less
Submitted 2 December, 2018; v1 submitted 15 November, 2017;
originally announced November 2017.
-
Variational optimization algorithms for uniform matrix product states
Authors:
V. Zauner-Stauber,
L. Vanderstraeten,
M. T. Fishman,
F. Verstraete,
J. Haegeman
Abstract:
We combine the Density Matrix Renormalization Group (DMRG) with Matrix Product State tangent space concepts to construct a variational algorithm for finding ground states of one dimensional quantum lattices in the thermodynamic limit. A careful comparison of this variational uniform Matrix Product State algorithm (VUMPS) with infinite Density Matrix Renormalization Group (IDMRG) and with infinite…
▽ More
We combine the Density Matrix Renormalization Group (DMRG) with Matrix Product State tangent space concepts to construct a variational algorithm for finding ground states of one dimensional quantum lattices in the thermodynamic limit. A careful comparison of this variational uniform Matrix Product State algorithm (VUMPS) with infinite Density Matrix Renormalization Group (IDMRG) and with infinite Time Evolving Block Decimation (ITEBD) reveals substantial gains in convergence speed and precision. We also demonstrate that VUMPS works very efficiently for Hamiltonians with long range interactions and also for the simulation of two dimensional models on infinite cylinders. The new algorithm can be conveniently implemented as an extension of an already existing DMRG implementation.
△ Less
Submitted 19 April, 2019; v1 submitted 24 January, 2017;
originally announced January 2017.
-
Compression of Correlation Matrices and an Efficient Method for Forming Matrix Product States of Fermionic Gaussian States
Authors:
Matthew T. Fishman,
Steven R. White
Abstract:
Here we present an efficient and numerically stable procedure for compressing a correlation matrix into a set of local unitary single-particle gates, which leads to a very efficient way of forming the matrix product state (MPS) approximation of a pure fermionic Gaussian state, such as the ground state of a quadratic Hamiltonian. The procedure involves successively diagonalizing subblocks of the co…
▽ More
Here we present an efficient and numerically stable procedure for compressing a correlation matrix into a set of local unitary single-particle gates, which leads to a very efficient way of forming the matrix product state (MPS) approximation of a pure fermionic Gaussian state, such as the ground state of a quadratic Hamiltonian. The procedure involves successively diagonalizing subblocks of the correlation matrix to isolate local states which are purely occupied or unoccupied. A small number of nearest neighbor unitary gates isolates each local state. The MPS of this state is formed by applying the many-body version of these gates to a product state. We treat the simple case of compressing the correlation matrix of spinless free fermions with definite particle number in detail, though the procedure is easily extended to fermions with spin and more general BCS states (utilizing the formalism of Majorana modes). We also present a DMRG-like algorithm to obtain the compressed correlation matrix directly from a hopping Hamiltonian. In addition, we discuss a slight variation of the procedure which leads to a simple construction of the multiscale entanglement renormalization ansatz (MERA) of a fermionic Gaussian state, and present a simple picture of orthogonal wavelet transforms in terms of the gate structure we present in this paper. As a simple demonstration we analyze the Su-Schrieffer-Heeger model (free fermions on a 1D lattice with staggered hopping amplitudes).
△ Less
Submitted 28 April, 2015;
originally announced April 2015.