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

Physical review. B, Solid state

Correlation function and susceptibility of site- and bond-diluted Heisenberg paramagnets

和子 川崎, Raza A. Tahir-Kheli

Published Mar 1, 1977• 9 citations
Publisher

Abstract

Generalizing the recent work of Collins, the wave-vector-dependent susceptibility and the static correlation function for a randomly diluted Heisenberg paramagnet on a Bravais lattice with quenched-site, or nearest-neighbor exchange-bond, dilution is studied by high-temperature series-expansion techniques. The first five coefficients for the spin correlation function $S(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}})$ and susceptibility $\ensuremath{\chi}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}})$ are calculated for arbitrary $\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}$, the sign of the exchange and the spin magnitude and the magnetic site $c$ and bond $p$ concentrations. Numerical results for these functions are presented for various spins and for $\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}$ along the symmetry directions (1,1,1) and (1,0,0). Because of the generality of the present series, their length is rather limited, and as such they are not suitable for a careful analysis of the critical region. Nonetheless, to get a qualitative feel for the accuracy of the present results near the transition region, we have used Pad\'e approximants for determining the correlation length $\ensuremath{\xi}(T)$ as a function of the reduced temperature $\ensuremath{\epsilon}=\frac{[T\ensuremath{-}{T}_{C(N)}(\mathrm{random})]}{{T}_{C(N)}(\mathrm{random})}$. In this manner, we estimate that the numerical accuracy of our results should be reasonable as long as the system is moderately above the transition temperature, e.g., $\ensuremath{\epsilon}\ensuremath{\gtrsim}0.5$.

DOI
10.1103/physrevb.15.2741

Related papers

The Randomly Diluted Quantum Heisenberg Antiferromagnet at Zero Temperature

Chuncheng Wan

We use linear spin wave theory and series expansion method to study the diluted quantum antiferromagnetic Heisenberg model at zero temperature near the p = 1 and p = $p\\sb{Q}$ respectively, where p is the concentration of the spins or bonds, and $p\\sb{Q}$ is the transition point from order to disorder of the quantum spin system. For the first case, our complete calculation to the first order in x $\\equiv$ 1 $-$ p gives a dispersion relation for long-wavelength spin waves in the form $\\omega(\\vec q)$ = $c(x)q + i\\gamma(x)q\\sp\\tau$. Our results for $c(x)$ agrees with previous work, where $c(x)$ = $c(0)\\lbrack 1 - \\alpha x + O(x\\sp2)\\rbrack$. In particular the coefficient $\\alpha$, which can be related to the helicity modulus and the uniform perpendicular susceptibility, diverges in the limit $d\\to 2$. One major new results is that $\\tau$ = $d - 1$, which means that d = 2 is the lower critical dimension at which infrared divergences affect the dynamical response due to vacancies. For the second case, the critical behavior of the quantum system near $p\\sb{Q}$ is found to be different from that of the percolation and pure systems. Thus there exists a new fixed point in the phase diagram, which specifies a new universality class. The correlation length exponents $\\nu$ (which comes from the relation $\\xi\\sim(p\\sb{Q} - p)\\sp\\nu$, where $\\xi$ is the correlation length) are estimated, which are 0.65 $\\pm$ 0.2 and 0.55 $\\pm$ 0.2 in 2 and 3 dimensions respectively. The exponents of crossover with temperature and time are determined. We also find evidence that $p\\sb{Q} > p\\sb{c}$.

•Jan 1, 1991•0 citations

Critical temperature of the Heisenberg model with random bond dilution

E Brown, J W Essam, C. M. Place

The method of high-temperature series expansions is used to study the dependence of T c on the concentration p of interactions in a randomly diluted S= 1 / 2 Heisenberg model. It is found that, for all three cubic lattices, the interval of p over which the Pade approximants to the series converge extends further towards the critical probability for percolation than for the corresponding site diluted model. In this interval the critical curve is linear and extrapolates to intersect the T=0 boundary at a value of p marginally above the critical probability for bond percolation. The anomalous dependence of the susceptibility exponent gamma on p which has recently been reported for the site diluted model is not observed for the model studied here.

Journal of Physics C Solid State Physics•Feb 7, 1975•19 citations

Quantum Phase Transition of Two-Dimensional Diluted Heisenberg Antiferromagnet

Synge Todo, Chitoshi Yasuda, Kiyoshi Katō +4 more

Quantum phase transition of site-diluted and bond-diluted Heisenberg antiferromagnets on square lattices is studied. By using the novel continuous-time loop algorithm, we perform quantum Monte Carlo simulations on quite larger lattices at extremely lower temperatures than the previous numerical studies. It is found that the antiferromagnetic long-range order at T = 0 persists so long as a cluster of magnetic sites percolates, that is, the critical concentration is equal to the classical percolation threshold, in both of the site-diluted and bond-diluted cases. Furthermore, we find that some critical exponents, such as the magnetization exponent β, are non-classical and strongly depend on the spin size S. On the other hand, we show that the correlation-length exponent ν is universal and is equal to the classical value (ν=4/3).

Progress of Theoretical Physics Supplement•Jan 1, 2000•2 citations

Theory of Critical-Point Scattering and Correlations. II. Heisenberg Models

Douglas S. Ritchie, Michael Ellis Fisher

The spin-spin correlation functions and the critical-scattering intensity for Heisenberg models of general spin, $S=\frac{1}{2}$ to $\ensuremath{\infty}$, on the sc, bcc, and fcc lattices are studied on the basis of high-temperature series expansions along the lines developed in Paper I [M. E. Fisher and R. J. Burford, Phys. Rev. 156, 583 (1967)]. Subject to increased uncertainties for low spin, it is concluded that the exponents $\ensuremath{\gamma}={1.375}_{\ensuremath{-}0.01}^{+0.02}$, $2\ensuremath{\nu}={1.405}_{\ensuremath{-}0.01}^{+0.02}$, and $\ensuremath{\eta}=0.043\ifmmode\pm\else\textpm\fi{}0.014$ describe all lattices and all spin. Explicit formulas are presented for the susceptibility/zero-angle scattering ${\ensuremath{\chi}}_{0}(T)$, for the inverse correlation length ${\ensuremath{\kappa}}_{1}(T)$, for the effective interaction range ${r}_{1}(T)$, and using the Fisher-Burford approximant, for the total scattering $\stackrel{^}{\ensuremath{\chi}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}},T)$. The shape parameter ${\ensuremath{\varphi}}_{c}$ attains the "universal" value ${\ensuremath{\varphi}}_{c}\ensuremath{\simeq}0$. 11 for large spin but shows signs of spin dependence (and lattice dependence) for low spin. At fixed k the scattering is predicted to display a maximum above ${T}_{c}$ determined by $\frac{{\ensuremath{\kappa}}_{1}({T}_{max})}{k}\ensuremath{\simeq}0.10 (\mathrm{for} S\ensuremath{\gtrsim}2)$ to 0.15. A detailed study is made of the structure dependence of the critical-point correlations ${〈{S}_{0}^{z}{S}_{\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}^{z}〉}_{c}$ for various models. This leads to the revised, universal estimate ${\ensuremath{\varphi}}_{c}\ensuremath{\simeq}0$. 15 for all three cubic lattice, spin-\textonehalf{} Ising models. The results are compared d briefly with various experiments which support $\ensuremath{\eta}\ensuremath{\gtrsim}0.05$.

Physical review. B, Solid state•Apr 1, 1972•196 citations

Magnetic Properties of the Random Mixture of the Classical Heisenberg Spins

Shigetoshi Katsura

The specific heat, the susceptibility, and the correlation function at zero field above the critical temperature of the random mixture (quenched site and bond problems) of the classical Heisenberg spins with nearest neighbor interaction were obtained exactly for the linear chain and for an infinite Bethe lattice (Bethe approximation of the two and three dimensional lattices) above the critical temperature. The results are simply expressed by the replacements of 2 cosh K → 4π (sinh K)/K and tanh K → L(K) (L(K) = Langevin function) for K = KAA, KAB, KBA, and KBB in the corresponding expressions of the random mixture of the Ising spins, and qualitative properties of both models are similar.

Canadian Journal of Physics•May 1, 1975•16 citations

Partial Test of the Universality Hypothesis: The Case of Different Coupling Strengths in Different Lattice Directions

Gerald Paul, Eugene Stanley

High-temperature series expansions are used to examine the dependence of critical-point exponents upon lattice anisotropy (different interaction strengths in different directions in the lattice). The two-spin correlation function ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ is calculated to tenth order in ($\frac{1}{{k}_{B}T}$) for the Ising Hamiltonian ${\mathcal{H}}_{l \mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}\stackrel{\mathrm{xy}}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}\stackrel{z}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}$ for a wide range of anisotropy parameters $R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}$ and for both the sc and fcc lattices; here the first summation is over all pairs of nearest-neighbor sites whose relative displacement vector ${\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{\mathrm{ij}}$ has no $z$ component, while the second summation is over all other pairs of nearest-neighbor sites. Hence for $R=0$, both the sc and fcc lattices reduce to the two-dimensional square lattice, while in the limit $R\ensuremath{\rightarrow}\ensuremath{\infty}$, the sc becomes a one-dimensional linear chain and the fcc becomes two noninteracting three-dimensional bcc lattices. The series for ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ are then used to obtain series of corresponding lengths for the specific heat, susceptibility, and second moment. Analysis of these series yields results consistent with the universality hypothesis of critical-point exponents. Specifically, it is found that when lattice anisotrophy is introduced, the critical-point exponents studied (the susceptibility exponent $\ensuremath{\gamma}$ and the correlation length exponent $\ensuremath{\nu}$) do not appear to change from their values for an isotropic lattice. The problem of next-nearest-neighbor interactions is treated using similar methods in Paper II of this series (and briefly discussed in this paper).

Physical review. B, Solid state•Apr 1, 1972•56 citations
Generalizing the recent work of Collins, the wave-vector-dependent susceptibility and the static correlation function for a randomly diluted Heisenberg paramagnet on a Bravais lattice with quenched-site, or nearest-neighbor exchange-bond, dilution is studied by high-temperature series-expansion techniques. The first five coefficients for the spin correlation function S(k\ensuremath→)S(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}})S(k\ensuremath→​) and susceptibility \ensuremathχ(k\ensuremath→)\ensuremath{\chi}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}})\ensuremathχ(k\ensuremath→​) are calculated for arbitrary k\ensuremath→\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}k\ensuremath→​, the sign of the exchange and the spin magnitude and the magnetic site ccc and bond ppp concentrations. Numerical results for these functions are presented for various spins and for k\ensuremath→\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}k\ensuremath→​ along the symmetry directions (1,1,1) and (1,0,0). Because of the generality of the present series, their length is rather limited, and as such they are not suitable for a careful analysis of the critical region. Nonetheless, to get a qualitative feel for the accuracy of the present results near the transition region, we have used Pad'e approximants for determining the correlation length \ensuremathξ(T)\ensuremath{\xi}(T)\ensuremathξ(T) as a function of the reduced temperature \ensuremathϵ=[T\ensuremath−TC(N)(random)]TC(N)(random)\ensuremath{\epsilon}=\frac{[T\ensuremath{-}{T}_{C(N)}(\mathrm{random})]}{{T}_{C(N)}(\mathrm{random})}\ensuremathϵ=TC(N)​(random)[T\ensuremath−TC(N)​(random)]​. In this manner, we estimate that the numerical accuracy of our results should be reasonable as long as the system is moderately above the transition temperature, e.g., \ensuremathϵ\ensuremath≳0.5\ensuremath{\epsilon}\ensuremath{\gtrsim}0.5\ensuremathϵ\ensuremath≳0.5.
We use linear spin wave theory and series expansion method to study the diluted quantum antiferromagnetic Heisenberg model at zero temperature near the p = 1 and p = psbQp\\sb{Q}psbQ respectively, where p is the concentration of the spins or bonds, and psbQp\\sb{Q}psbQ is the transition point from order to disorder of the quantum spin system. For the first case, our complete calculation to the first order in x equiv\\equivequiv 1 −-− p gives a dispersion relation for long-wavelength spin waves in the form omega(vecq)\\omega(\\vec q)omega(vecq) = c(x)q+igamma(x)qsptauc(x)q + i\\gamma(x)q\\sp\\tauc(x)q+igamma(x)qsptau. Our results for c(x)c(x)c(x) agrees with previous work, where c(x)c(x)c(x) = c(0)lbrack1−alphax+O(xsp2)rbrackc(0)\\lbrack 1 - \\alpha x + O(x\\sp2)\\rbrackc(0)lbrack1−alphax+O(xsp2)rbrack. In particular the coefficient alpha\\alphaalpha, which can be related to the helicity modulus and the uniform perpendicular susceptibility, diverges in the limit dto2d\\to 2dto2. One major new results is that tau\\tautau = d−1d - 1d−1, which means that d = 2 is the lower critical dimension at which infrared divergences affect the dynamical response due to vacancies. For the second case, the critical behavior of the quantum system near psbQp\\sb{Q}psbQ is found to be different from that of the percolation and pure systems. Thus there exists a new fixed point in the phase diagram, which specifies a new universality class. The correlation length exponents nu\\nunu (which comes from the relation xisim(psbQ−p)spnu\\xi\\sim(p\\sb{Q} - p)\\sp\\nuxisim(psbQ−p)spnu, where xi\\xixi is the correlation length) are estimated, which are 0.65 pm\\pmpm 0.2 and 0.55 pm\\pmpm 0.2 in 2 and 3 dimensions respectively. The exponents of crossover with temperature and time are determined. We also find evidence that psbQ>psbcp\\sb{Q} > p\\sb{c}psbQ>psbc.
The method of high-temperature series expansions is used to study the dependence of T c on the concentration p of interactions in a randomly diluted S= 1 / 2 Heisenberg model. It is found that, for all three cubic lattices, the interval of p over which the Pade approximants to the series converge extends further towards the critical probability for percolation than for the corresponding site diluted model. In this interval the critical curve is linear and extrapolates to intersect the T=0 boundary at a value of p marginally above the critical probability for bond percolation. The anomalous dependence of the susceptibility exponent gamma on p which has recently been reported for the site diluted model is not observed for the model studied here.
Quantum phase transition of site-diluted and bond-diluted Heisenberg antiferromagnets on square lattices is studied. By using the novel continuous-time loop algorithm, we perform quantum Monte Carlo simulations on quite larger lattices at extremely lower temperatures than the previous numerical studies. It is found that the antiferromagnetic long-range order at T = 0 persists so long as a cluster of magnetic sites percolates, that is, the critical concentration is equal to the classical percolation threshold, in both of the site-diluted and bond-diluted cases. Furthermore, we find that some critical exponents, such as the magnetization exponent β, are non-classical and strongly depend on the spin size S. On the other hand, we show that the correlation-length exponent ν is universal and is equal to the classical value (ν=4/3).
The spin-spin correlation functions and the critical-scattering intensity for Heisenberg models of general spin, S=12S=\frac{1}{2}S=21​ to $\ensuremath{\infty}$, on the sc, bcc, and fcc lattices are studied on the basis of high-temperature series expansions along the lines developed in Paper I [M. E. Fisher and R. J. Burford, Phys. Rev. 156, 583 (1967)]. Subject to increased uncertainties for low spin, it is concluded that the exponents $\ensuremath{\gamma}={1.375}_{\ensuremath{-}0.01}^{+0.02}$, $2\ensuremath{\nu}={1.405}_{\ensuremath{-}0.01}^{+0.02}$, and $\ensuremath{\eta}=0.043\ifmmode\pm\else\textpm\fi{}0.014$ describe all lattices and all spin. Explicit formulas are presented for the susceptibility/zero-angle scattering ${\ensuremath{\chi}}_{0}(T)$, for the inverse correlation length ${\ensuremath{\kappa}}_{1}(T)$, for the effective interaction range r1(T){r}_{1}(T)r1​(T), and using the Fisher-Burford approximant, for the total scattering $\stackrel{^}{\ensuremath{\chi}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}},T)$. The shape parameter ${\ensuremath{\varphi}}_{c}$ attains the "universal" value ${\ensuremath{\varphi}}_{c}\ensuremath{\simeq}0$. 11 for large spin but shows signs of spin dependence (and lattice dependence) for low spin. At fixed k the scattering is predicted to display a maximum above Tc{T}_{c}Tc​ determined by $\frac{{\ensuremath{\kappa}}_{1}({T}_{max})}{k}\ensuremath{\simeq}0.10 (\mathrm{for} S\ensuremath{\gtrsim}2)$ to 0.15. A detailed study is made of the structure dependence of the critical-point correlations ${〈{S}_{0}^{z}{S}_{\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}^{z}〉}_{c}$ for various models. This leads to the revised, universal estimate ${\ensuremath{\varphi}}_{c}\ensuremath{\simeq}0$. 15 for all three cubic lattice, spin-\textonehalf{} Ising models. The results are compared d briefly with various experiments which support $\ensuremath{\eta}\ensuremath{\gtrsim}0.05$.
The specific heat, the susceptibility, and the correlation function at zero field above the critical temperature of the random mixture (quenched site and bond problems) of the classical Heisenberg spins with nearest neighbor interaction were obtained exactly for the linear chain and for an infinite Bethe lattice (Bethe approximation of the two and three dimensional lattices) above the critical temperature. The results are simply expressed by the replacements of 2 cosh K → 4π (sinh K)/K and tanh K → L(K) (L(K) = Langevin function) for K = KAA, KAB, KBA, and KBB in the corresponding expressions of the random mixture of the Ising spins, and qualitative properties of both models are similar.
High-temperature series expansions are used to examine the dependence of critical-point exponents upon lattice anisotropy (different interaction strengths in different directions in the lattice). The two-spin correlation function ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ is calculated to tenth order in (1kBT\frac{1}{{k}_{B}T}kB​T1​) for the Ising Hamiltonian ${\mathcal{H}}_{l \mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}\stackrel{\mathrm{xy}}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}\stackrel{z}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}$ for a wide range of anisotropy parameters $R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}$ and for both the sc and fcc lattices; here the first summation is over all pairs of nearest-neighbor sites whose relative displacement vector ${\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{\mathrm{ij}}$ has no zzz component, while the second summation is over all other pairs of nearest-neighbor sites. Hence for R=0R=0R=0, both the sc and fcc lattices reduce to the two-dimensional square lattice, while in the limit $R\ensuremath{\rightarrow}\ensuremath{\infty}$, the sc becomes a one-dimensional linear chain and the fcc becomes two noninteracting three-dimensional bcc lattices. The series for ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ are then used to obtain series of corresponding lengths for the specific heat, susceptibility, and second moment. Analysis of these series yields results consistent with the universality hypothesis of critical-point exponents. Specifically, it is found that when lattice anisotrophy is introduced, the critical-point exponents studied (the susceptibility exponent $\ensuremath{\gamma}$ and the correlation length exponent $\ensuremath{\nu}$) do not appear to change from their values for an isotropic lattice. The problem of next-nearest-neighbor interactions is treated using similar methods in Paper II of this series (and briefly discussed in this paper).