Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2308.16407 · Journal of the Physical Society of Japan

General Formula for the Green’s Function Approach to the Spin-1/2 Antiferromagnetic Heisenberg Model

Anonymous, Takao Morinari

Published Jan 12, 2024• 3 citations
arXivHTMLPDF

Abstract

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 $J_1$-$J_2$ model. For the cubic lattice case, the Green's function approach provides a good estimation for the transition temperature. Regarding the $J_1$-$J_2$ 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.

DOI
10.7566/jpsj.93.024704

Related papers

Thermodynamic properties of one-dimensional spin-1 antiferromagnetic Heisenberg chains: Green’s-function approach

Shu-Qing Bao, Hu Zhao, Jue-Lian Shen +1 more

The one-dimensional spin-1 Heisenberg antiferromagnetic chains are studied by a generalized two-time Green's-function method at finite temperature. 〈${\mathit{S}}_{\mathit{i}}^{\mathit{z}}$〉 is set to zero at each site to guarantee the lack of the long-range order in the one-dimensional Heisenberg system. The Green's functions are decoupled in terms of the correlation functions of spin operators. A set of self-consistent equations of the correlation functions are derived and solved numerically. Correlation functions of spin operators, thermodynamic properties such as the internal energy, and specific heat in the entire temperature region are obtained. The Haldane gap in the excitation spectrum appears naturally in the analytical result. When k is zero, there is a gap 2\ensuremath{\Delta}, which is 1.0J and close to the theoretical result of White and Huse. And there is a broad maximum in the temperature dependence of the specific heat. The properties of the thermodynamic quantities are consistent with the numerical results on the finite chains. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Jan 1, 1996•24 citations

Green's-Function Formalism of the One-Dimensional Heisenberg Spin System

Jun Kondo, Yamaji K.

The one-dimensional Heisenberg model with S = 1/2 is treated with the use of the two-time Green's functions. The hierarchy of the equations of motion of the Green'f functions is decoupled at a stage one-step further than Tyablikov's decoupling. The thermal average of the spin component, ≪ Sz >, is set to zero, becaure the long-range order does not exist in one dimension. Instead, our Green's functions are expressed in terms of the correlation functions cn ≡4≪ S0zSnz >. The Green's function is essentially of the form representing undamped spin waves, whose spactrum depends on c1, c2 and one more parameter. They are determined by the requirement that c1 and c2 should be self-consistent and that c0 should be unity. The self-consistency equations have been solved analytically at high- and low-temperature limits, and also solved numerically in the whole range of the temperature. Thermodynamic quantities have been calculated using these solutions. It has turned out that the theory gives the correct high-temperature expansion for the thermodynamic quantities and the correlation functions. The latter is expressed by cn = (J / 4 kBT )n. In the case of ferromagnetic coupling, the correlation function cn at T = 0 is equal to 1/3 for all n's. This is what is expected from the correct ground state of the ferromagnetic Heisenberg system. The spin-wave spectrum at T = 0 also agrees with the correct one. At T ≪J/ kB, we find that the specific heat goes as T1/2 and the susceptibility goes as T-2. The gross feature of the temperature-dependence of the thermodynamic quantities agrees with Bonner and Fisher. In the case of the antiferromagnetic Heisenberg model, we find c1 = -0.55407 and c2 = 0.16100 at T = 0, which are fairly close to the exact values. The thermodynamic quantities are also in gross agreement with Bonner and Fisher.

Progress of Theoretical Physics•Mar 1, 1972•246 citations

Thermodynamics of frustrated ferromagnetic spin- 12 Heisenberg chains: Role of interchain coupling

P. Müller, Johannes Richter, Dieter Ihle

The thermodynamics of coupled frustrated ferromagnetic chains is studied within a spin-rotation-invariant Green's-function approach. We consider an isotropic Heisenberg spin-half system with a ferromagnetic in-chain coupling ${J}_{1}<0$ between nearest neighbors and a frustrating antiferromagnetic next-nearest-neighbor in-chain coupling ${J}_{2}>0$. We focus on the moderate strength of frustration ${J}_{2}<|{J}_{1}|/4$ such that the in-chain spin-spin correlations are predominantly ferromagnetic. We consider two interchain couplings (ICs) ${J}_{\ensuremath{\perp},y}$ and ${J}_{\ensuremath{\perp},z}$, corresponding to the two axes perpendicular to the chain, where ferromagnetic as well as antiferromagnetic ICs are taken into account. We discuss the influence of frustration on the ground-state properties for antiferromagnetic ICs, where the ground state is of a quantum nature. The major part of our study is devoted to the finite-temperature properties. We calculate the critical temperature ${T}_{c}$ as a function of the competing exchange couplings ${J}_{2},{J}_{\ensuremath{\perp},y},{J}_{\ensuremath{\perp},z}$. We find that for fixed ICs, ${T}_{c}$ decreases monotonically with increasing frustration ${J}_{2}$, where as ${J}_{2}\ensuremath{\rightarrow}|{J}_{1}|/4$ the ${T}_{c}({J}_{2})$ curve drops down rapidly. To characterize the magnetic ordering below and above ${T}_{c}$, we calculate the spin-spin correlation functions $\ensuremath{\langle}{\mathbf{S}}_{0}{\mathbf{S}}_{\mathbf{R}}\ensuremath{\rangle}$, the magnetic order parameter $M$, the uniform static susceptibility ${\ensuremath{\chi}}_{0}$, as well as the correlation length $\ensuremath{\xi}$. Moreover, we discuss the specific heat ${C}_{V}$ and the temperature dependence of the excitation spectrum ${\ensuremath{\omega}}_{\mathbf{q}}$. As ${J}_{2}\ensuremath{\rightarrow}|{J}_{1}|/4$, some unusual frustration-induced features were found, such as an increase of the in-chain spin stiffness (in the case of ferromagnetic ICs) or of the in-chain spin-wave velocity (in the case of antiferromagnetic ICs) with growing temperature.

arXiv:1703.04573•Physical review. B./Physical review. B•Apr 5, 2017•15 citations

Spin excitations and thermodynamics of the t-J model on the honeycomb lattice

Vladimirov, A., Ihle, D., Plakida, N.

We present a spin-rotation-invariant Green-function theory for the dynamic spin susceptibility in the spin-1/2 antiferromagnetic t-J Heisenberg model on the honeycomb lattice. Employing a generalized mean-field approximation for arbitrary temperatures and hole dopings, the electronic spectrum of excitations, the spin-excitation spectrum and thermodynamic quantities (two-spin correlation functions, staggered magnetization, magnetic susceptibility, correlation length) are calculated by solving a coupled system of self-consistency equations for the correlation functions. The temperature and doping dependence of the magnetic (uniform static) susceptibility is ascribed to antiferromagnetic short-range order. Our results on the doping dependencies of the magnetization and susceptibility are analyzed in comparison with previous results for the t-J model on the square lattice.

•Aug 1, 2018•4 citations

Spin excitations and thermodynamics of the t-J model on the honeycomb lattice

Artem Vladimirov, Dieter Ihle, Nikolay M. Plakida

We present a spin-rotation-invariant Green-function theory for the dynamic spin susceptibility in the spin-1/2 antiferromagnetic t - J Heisenberg model on the honeycomb lattice. Employing a generalized mean-field approximation for arbitrary temperatures and hole dopings, the electronic spectrum of excitations, the spin-excitation spectrum and thermodynamic quantities (two-spin correlation functions, staggered magnetization, magnetic susceptibility, correlation length) are calculated by solving a coupled system of self-consistency equations for the correlation functions. The temperature and doping dependence of the magnetic (uniform static) susceptibility is ascribed to antiferromagnetic short-range order. Our results on the doping dependencies of the magnetization and susceptibility are analyzed in comparison with previous results for the t - J model on the square lattice.

arXiv:1803.03140•The European Physical Journal B•Aug 1, 2018•5 citations

Thermodynamics of the one-dimensional frustrated Heisenberg ferromagnet with arbitrary spin

M. Härtel, Johannes Richter, Dieter Ihle +2 more

The thermodynamic quantities (spin-spin correlation functions $\ensuremath{\langle}{\mathbf{S}}_{0}{\mathbf{S}}_{n}\ensuremath{\rangle}$, correlation length $\ensuremath{\xi}$, spin susceptibility $\ensuremath{\chi}$, and specific heat ${C}_{V}$) of the frustrated one-dimensional ${J}_{1}$-${J}_{2}$ Heisenberg ferromagnet with arbitrary spin quantum number $S$ below the quantum critical point, i.e., for ${J}_{2}<|{J}_{1}|/4$, are calculated using a rotation-invariant Green-function formalism and full diagonalization as well as a finite-temperature Lanczos technique for finite chains of up to $N=20$ sites. The low-temperature behavior of the susceptibility $\ensuremath{\chi}$ and the correlation length $\ensuremath{\xi}$ are well described by $\ensuremath{\chi}=\frac{2}{3}{S}^{4}\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right){T}^{\ensuremath{-}2}+A{S}^{5/2}{\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right)}^{1/2}{T}^{\ensuremath{-}3/2}$ and $\ensuremath{\xi}={S}^{2}\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right){T}^{\ensuremath{-}1}+B{S}^{1/2}{\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right)}^{1/2}{T}^{\ensuremath{-}1/2}$ with $A\ensuremath{\approx}1.1\ensuremath{\cdots}1.2$ and $B\ensuremath{\approx}0.84\ensuremath{\cdots}0.89$. The vanishing of the factors in front of the temperature at ${J}_{2}=|{J}_{1}|/4$ indicates a change of the critical behavior of $\ensuremath{\chi}$ and $\ensuremath{\xi}$ at $T\ensuremath{\rightarrow}0$. The specific heat may exhibit an additional frustration-induced low-temperature maximum when approaching the quantum critical point. Such a maximum appears for $S=1/2$ and $S=1$, but was not found for $S>1$.

arXiv:1106.1013•Physical Review B•Sep 6, 2011•25 citations
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 J1J_1J1​-J2J_2J2​ model. For the cubic lattice case, the Green's function approach provides a good estimation for the transition temperature. Regarding the J1J_1J1​-J2J_2J2​ 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.
The one-dimensional spin-1 Heisenberg antiferromagnetic chains are studied by a generalized two-time Green's-function method at finite temperature. 〈Siz{\mathit{S}}_{\mathit{i}}^{\mathit{z}}Siz​〉 is set to zero at each site to guarantee the lack of the long-range order in the one-dimensional Heisenberg system. The Green's functions are decoupled in terms of the correlation functions of spin operators. A set of self-consistent equations of the correlation functions are derived and solved numerically. Correlation functions of spin operators, thermodynamic properties such as the internal energy, and specific heat in the entire temperature region are obtained. The Haldane gap in the excitation spectrum appears naturally in the analytical result. When k is zero, there is a gap 2\ensuremath{\Delta}, which is 1.0J and close to the theoretical result of White and Huse. And there is a broad maximum in the temperature dependence of the specific heat. The properties of the thermodynamic quantities are consistent with the numerical results on the finite chains. \textcopyright{} 1996 The American Physical Society.
The one-dimensional Heisenberg model with S = 1/2 is treated with the use of the two-time Green's functions. The hierarchy of the equations of motion of the Green'f functions is decoupled at a stage one-step further than Tyablikov's decoupling. The thermal average of the spin component, ≪ Sz >, is set to zero, becaure the long-range order does not exist in one dimension. Instead, our Green's functions are expressed in terms of the correlation functions cn ≡4≪ S0zSnz >. The Green's function is essentially of the form representing undamped spin waves, whose spactrum depends on c1, c2 and one more parameter. They are determined by the requirement that c1 and c2 should be self-consistent and that c0 should be unity. The self-consistency equations have been solved analytically at high- and low-temperature limits, and also solved numerically in the whole range of the temperature. Thermodynamic quantities have been calculated using these solutions. It has turned out that the theory gives the correct high-temperature expansion for the thermodynamic quantities and the correlation functions. The latter is expressed by cn = (J / 4 kBT )n. In the case of ferromagnetic coupling, the correlation function cn at T = 0 is equal to 1/3 for all n's. This is what is expected from the correct ground state of the ferromagnetic Heisenberg system. The spin-wave spectrum at T = 0 also agrees with the correct one. At T ≪J/ kB, we find that the specific heat goes as T1/2 and the susceptibility goes as T-2. The gross feature of the temperature-dependence of the thermodynamic quantities agrees with Bonner and Fisher. In the case of the antiferromagnetic Heisenberg model, we find c1 = -0.55407 and c2 = 0.16100 at T = 0, which are fairly close to the exact values. The thermodynamic quantities are also in gross agreement with Bonner and Fisher.
The thermodynamics of coupled frustrated ferromagnetic chains is studied within a spin-rotation-invariant Green's-function approach. We consider an isotropic Heisenberg spin-half system with a ferromagnetic in-chain coupling J1<0{J}_{1}<0J1​<0 between nearest neighbors and a frustrating antiferromagnetic next-nearest-neighbor in-chain coupling J2>0{J}_{2}>0J2​>0. We focus on the moderate strength of frustration J2<∣J1∣/4{J}_{2}<|{J}_{1}|/4J2​<∣J1​∣/4 such that the in-chain spin-spin correlations are predominantly ferromagnetic. We consider two interchain couplings (ICs) ${J}_{\ensuremath{\perp},y}$ and ${J}_{\ensuremath{\perp},z}$, corresponding to the two axes perpendicular to the chain, where ferromagnetic as well as antiferromagnetic ICs are taken into account. We discuss the influence of frustration on the ground-state properties for antiferromagnetic ICs, where the ground state is of a quantum nature. The major part of our study is devoted to the finite-temperature properties. We calculate the critical temperature Tc{T}_{c}Tc​ as a function of the competing exchange couplings ${J}_{2},{J}_{\ensuremath{\perp},y},{J}_{\ensuremath{\perp},z}$. We find that for fixed ICs, Tc{T}_{c}Tc​ decreases monotonically with increasing frustration J2{J}_{2}J2​, where as ${J}_{2}\ensuremath{\rightarrow}|{J}_{1}|/4$ the Tc(J2){T}_{c}({J}_{2})Tc​(J2​) curve drops down rapidly. To characterize the magnetic ordering below and above Tc{T}_{c}Tc​, we calculate the spin-spin correlation functions $\ensuremath{\langle}{\mathbf{S}}_{0}{\mathbf{S}}_{\mathbf{R}}\ensuremath{\rangle}$, the magnetic order parameter MMM, the uniform static susceptibility ${\ensuremath{\chi}}_{0}$, as well as the correlation length $\ensuremath{\xi}$. Moreover, we discuss the specific heat CV{C}_{V}CV​ and the temperature dependence of the excitation spectrum ${\ensuremath{\omega}}_{\mathbf{q}}$. As ${J}_{2}\ensuremath{\rightarrow}|{J}_{1}|/4$, some unusual frustration-induced features were found, such as an increase of the in-chain spin stiffness (in the case of ferromagnetic ICs) or of the in-chain spin-wave velocity (in the case of antiferromagnetic ICs) with growing temperature.
We present a spin-rotation-invariant Green-function theory for the dynamic spin susceptibility in the spin-1/2 antiferromagnetic t-J Heisenberg model on the honeycomb lattice. Employing a generalized mean-field approximation for arbitrary temperatures and hole dopings, the electronic spectrum of excitations, the spin-excitation spectrum and thermodynamic quantities (two-spin correlation functions, staggered magnetization, magnetic susceptibility, correlation length) are calculated by solving a coupled system of self-consistency equations for the correlation functions. The temperature and doping dependence of the magnetic (uniform static) susceptibility is ascribed to antiferromagnetic short-range order. Our results on the doping dependencies of the magnetization and susceptibility are analyzed in comparison with previous results for the t-J model on the square lattice.
We present a spin-rotation-invariant Green-function theory for the dynamic spin susceptibility in the spin-1/2 antiferromagnetic t - J Heisenberg model on the honeycomb lattice. Employing a generalized mean-field approximation for arbitrary temperatures and hole dopings, the electronic spectrum of excitations, the spin-excitation spectrum and thermodynamic quantities (two-spin correlation functions, staggered magnetization, magnetic susceptibility, correlation length) are calculated by solving a coupled system of self-consistency equations for the correlation functions. The temperature and doping dependence of the magnetic (uniform static) susceptibility is ascribed to antiferromagnetic short-range order. Our results on the doping dependencies of the magnetization and susceptibility are analyzed in comparison with previous results for the t - J model on the square lattice.
The thermodynamic quantities (spin-spin correlation functions $\ensuremath{\langle}{\mathbf{S}}_{0}{\mathbf{S}}_{n}\ensuremath{\rangle}$, correlation length $\ensuremath{\xi}$, spin susceptibility $\ensuremath{\chi}$, and specific heat CV{C}_{V}CV​) of the frustrated one-dimensional J1{J}_{1}J1​-J2{J}_{2}J2​ Heisenberg ferromagnet with arbitrary spin quantum number SSS below the quantum critical point, i.e., for J2<∣J1∣/4{J}_{2}<|{J}_{1}|/4J2​<∣J1​∣/4, are calculated using a rotation-invariant Green-function formalism and full diagonalization as well as a finite-temperature Lanczos technique for finite chains of up to N=20N=20N=20 sites. The low-temperature behavior of the susceptibility $\ensuremath{\chi}$ and the correlation length $\ensuremath{\xi}$ are well described by $\ensuremath{\chi}=\frac{2}{3}{S}^{4}\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right){T}^{\ensuremath{-}2}+A{S}^{5/2}{\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right)}^{1/2}{T}^{\ensuremath{-}3/2}$ and $\ensuremath{\xi}={S}^{2}\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right){T}^{\ensuremath{-}1}+B{S}^{1/2}{\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right)}^{1/2}{T}^{\ensuremath{-}1/2}$ with $A\ensuremath{\approx}1.1\ensuremath{\cdots}1.2$ and $B\ensuremath{\approx}0.84\ensuremath{\cdots}0.89$. The vanishing of the factors in front of the temperature at J2=∣J1∣/4{J}_{2}=|{J}_{1}|/4J2​=∣J1​∣/4 indicates a change of the critical behavior of $\ensuremath{\chi}$ and $\ensuremath{\xi}$ at $T\ensuremath{\rightarrow}0$. The specific heat may exhibit an additional frustration-induced low-temperature maximum when approaching the quantum critical point. Such a maximum appears for S=1/2S=1/2S=1/2 and S=1S=1S=1, but was not found for S>1S>1S>1.