We analyze the tensor network loop cluster expansion, introduced in [G. Park, J. Gray, and G. K.-L. Chan, Phys. Rev. B 112, 174310 (2025)] as a systematic correction to belief propagation, in the context of general quantum many-body problems. We provide numerical examples of the accuracy and practical applicability of the approach for the computation of ground-state observables for high bond dimension tensor networks, in two- and three-dimensions, with open and periodic boundary conditions, and for spin and fermion problems. We find that the contraction error converges approximately exponentially with cluster size, enabling accurate local observable and energy estimates for many systems where standard contraction methods are otherwise impractical.
The
d-dimensional Ising transfer matrix is expressed exactly in terms of Fermi-type operators. In two dimensions this leads to the standard solution of the model. In three or more dimensions the transfer matrix is no longer a quadratic form in the Fermi variables, but rather contains quartic and higher-order terms. However, it can be expressed as a well-defined expansion about a quadratic form in Fermi operators. The zeroth-order transfer-matrix-called herein the "fermion model"---can be exactly diagonalized, and is well approximated in the critical region by ${V}^{(0)}\ensuremath{\approx}{[2sinh(2K)]}^{\frac{N}{2}}\mathrm{exp}\left[\ensuremath{-}N\ensuremath{\int}d\stackrel{\ensuremath{\rightarrow}}{\mathrm{q}}{\left({\ensuremath{\kappa}}^{2}+\ensuremath{\Sigma}{n=1}^{d\ensuremath{-}1}\ensuremath{\Sigma}{m=1}^{d\ensuremath{-}1}{q}_{n}{q}_{m}\right)}^{\frac{1}{2}}t(\stackrel{\ensuremath{\rightarrow}}{\mathrm{q}})\right]$where
N is the number of sites per layer, $\ensuremath{\mu}$ is the inverse correlation length, $t(\stackrel{\ensuremath{\rightarrow}}{\mathrm{q}})$ takes values \ifmmode\pm\else\textpm\fi{}1/2,
K is the reduced inverse temperature, and the integral extends over a Brillouin zone of the ($d\ensuremath{-}1$)-dimensional layer. The critical point of the model is located by the equation ${x}_{c}^{d}+{x}_{c}^{d\ensuremath{-}1}+{x}_{c}=1$, where ${x}_{c}={e}^{\ensuremath{-}2{K}_{c}}$. The critical exponent $\ensuremath{\nu}$ describing the divergence of the correlation length at
Kc is found to take on its two-dimensional value, $\ensuremath{\nu}=1$, in all dimensions
d. The possibility of using the fermion model as the basis for an expansion of the Ising critical exponents as a function of $\ensuremath{\delta}(=d\ensuremath{-}2)$ is pointed out.
The generalisation of the quantum inverse scattering method to the study of direct and inverse problems for the multicomponent non-linear Schrodinger model of bosons or fermions with repulsive coupling is made. Two sets of Yang-Baxter equations are solved to obtain commutation relations between the scattering state operators. The eigenfunctions have been constructed for the infinite number of conserved quantities and the eigenvalues of the first three conserved quantities-number of particles, momentum and energy-are obtained. The global Izergin-Korepin relations and the relations between the quantum Jost functions are derived, and from them, the quantum Gel'fand-Levitan equations are established. Finally, the series expansion for the field operators in terms of the scattering state operators is written out explicitly.