A wide range of analytical and numerical methods are available to study quantum spin systems. However, the complexity of spin correlations and interactions limits their applicability to specific temperature ranges. The analytical approach utilizing Green's function has proved advantageous, as it allows for formulation without restrictions on the presence of long-range order and facilitates estimation of the spin excitation spectrum and thermodynamic quantities across the entire temperature range. In this work, we present a generalized formulation of the Green's function method that can be applied to diverse spin systems. As specific applications, we consider the hypercubic lattice and the
J1-
J2 model. For the cubic lattice case, the Green's function approach provides a good estimation for the transition temperature. Regarding the
J1-
J2 model, we include nematic correlations in the analysis and find no signature of such correlations, though accurate numerical calculations are required in the presence of strong frustration. Although our focus is on the spin one-half antiferromagnetic Heisenberg model on an arbitrary lattice, the Green's function approach can be generalized to incorporate other interactions and higher spin values.
Since the breakthrough of twistronics a plethora of topological phenomena in correlated systems has appeared. These devices can be typically analyzed in terms of lattice models using Green's function techniques. In this work we introduce a general method to obtain the boundary Green's function of such models taking advantage of the numerical Faddeev-LeVerrier algorithm to circumvent some analytical constraints of previous works. We illustrate our formalism analyzing the edge features of a Chern insulator, the Kitaev square lattice model for a topological superconductor and the Checkerboard lattice hosting topological flat bands. The efficiency and accuracy of the method is demonstrated by comparison to recursive Green's function calculations.
The linked-cluster expansion technique for the high-temperature expansion of spin model is reviewed. A new algorithm for the computation of three-point and higher Green’s functions is presented. Series are computed for all components of two-point Green’s functions for a generalized 3D Ising model, to 25th order on the bcc lattice and to 23rd order on the sc lattice. Series for zero-momentum four-, six-, and eight-point functions are computed to 21st, 19th, and 17th order respectively on the bcc lattice. 1