The Green's-function method is employed to study the two-dimensional system of classical
O(N)-vector spins. The second-order double-time retarded Green's functions of the corresponding quantum Hamiltonian system on a one-dimensional lattice are decoupled approximately in a selfconsistent manner. The expressions for the mass gap and the spin-correlation function thus obtained reproduce the known exact results in the large-
N limit and in the strong-coupling region. Some directions are discussed to improve the present method in the region of weak coupling with
N not large.
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.