Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain
We calculate the low temperature asymptotics of a function γ that generates the temperature dependence of all static correlation functions of the isotropic Heisenberg chain.
Saved in:
| Published in: | Symmetry, Integrability and Geometry: Methods and Applications |
|---|---|
| Date: | 2010 |
| Main Authors: | , , |
| Format: | Article |
| Language: | English |
| Published: |
Інститут математики НАН України
2010
|
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/146513 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain / N. Crampé, F. Göhmann, A. Klümper // Symmetry, Integrability and Geometry: Methods and Applications. — 2010. — Т. 6. — Бібліогр.: 14 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859771140545708032 |
|---|---|
| author | Crampé, N. Göhmann, F. Klümper, A. |
| author_facet | Crampé, N. Göhmann, F. Klümper, A. |
| citation_txt | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain / N. Crampé, F. Göhmann, A. Klümper // Symmetry, Integrability and Geometry: Methods and Applications. — 2010. — Т. 6. — Бібліогр.: 14 назв. — англ. |
| collection | DSpace DC |
| container_title | Symmetry, Integrability and Geometry: Methods and Applications |
| description | We calculate the low temperature asymptotics of a function γ that generates
the temperature dependence of all static correlation functions of the isotropic Heisenberg
chain.
|
| first_indexed | 2025-12-02T07:01:51Z |
| format | Article |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 6 (2010), 082, 5 pages
Universal Low Temperature Asymptotics of
the Correlation Functions of the Heisenberg Chain?
Nicolas CRAMPÉ †, Frank GÖHMANN ‡ and Andreas KLÜMPER ‡
† LPTA, UMR 5207 CNRS-UM2, Place Eugène Bataillon, 34095 Montpellier Cedex 5, France
E-mail: ncrampe@um2.fr
‡ Fachbereich C – Physik, Bergische Universität Wuppertal, 42097 Wuppertal, Germany
E-mail: goehmann@physik.uni-wuppertal.de, kluemper@uni-wuppertal.de
Received August 18, 2010, in final form October 04, 2010; Published online October 09, 2010
doi:10.3842/SIGMA.2010.082
Abstract. We calculate the low temperature asymptotics of a function γ that generates
the temperature dependence of all static correlation functions of the isotropic Heisenberg
chain.
Key words: correlation functions; quantum spin chains; thermodynamic Bethe ansatz
2010 Mathematics Subject Classification: 81Q80; 82B23
1 Introduction
Over the past few years the mathematical structure of the static correlation functions of the
XXZ chain was largely resolved. After an appropriate regularization by a disorder parameter
they all factorize into polynomials in only two functions ρ and ω [8]. These are the one-point
function and a special neighbor two-point function which, in turn, can be represented as integrals
over solutions of certain linear and non-linear integral equations [2]. This resembles much the
situation with free fermions, and what is behind is indeed a remarkable fermionic structure on
the space of quasi-local operators acting on the spin chain [5]. It allows us, for instance, to
calculate short-range correlators with high numerical precision directly in the thermodynamic
limit [1, 12].
The low temperature asymptotics of ρ and ω universally determines the low temperature
properties of all static correlation functions. In this short note we obtain the low temperature
asymptotics in the special case of the isotropic Hamiltonian
H = J
∑
j
(
σx
j−1σ
x
j + σy
j−1σ
y
j + σz
j−1σ
z
j
)
(1.1)
with no magnetic field applied and vanishing disorder parameter. Then ρ = 1 and we are left
with only one function (and its derivatives) which, up to a trivial factor, is the function γ defined
in [3].
2 Definition of the basic function γ
For our purpose here it is convenient to introduce the function γ within the context of a special
realization of a six-vertex model (see e.g. [4]) and its associated quantum transfer matrix [10].
?This paper is a contribution to the Proceedings of the International Workshop “Recent Advances in Quantum
Integrable Systems”. The full collection is available at http://www.emis.de/journals/SIGMA/RAQIS2010.html
mailto:ncrampe@um2.fr
mailto:goehmann@physik.uni-wuppertal.de
mailto:kluemper@uni-wuppertal.de
http://dx.doi.org/10.3842/SIGMA.2010.082
http://www.emis.de/journals/SIGMA/RAQIS2010.html
2 N. Crampé, F. Göhmann and A. Klümper
By definition the latter has 2(N+M) vertical lines alternating in direction and carrying spectral
parameters
u,−u, u,−u, . . . ,−u︸ ︷︷ ︸
2N
, u′ + µ1, µ1 − u′, u′ + µ1, µ1 − u′, . . . , µ1 − u′︸ ︷︷ ︸
2M
.
The spectral parameter on the horizontal line will be denoted µ2. We consider this system in
the limit N ,M → +∞ with the fine tuning uN = i J
T and u′M = i δ
T . With an appropriate
overall normalization the largest eigenvalue Λ(µ2, µ1) is given by
ln(Λ(µ2, µ1)) =
4πJ
T
K(µ2) +
4πδ
T
K(µ2 − µ1)
+
∫ ∞
−∞
dt
ln
[
(1 + b(t, µ1))(1 + b(t, µ1))
]
2 cosh(π(µ2 − t))
. (2.1)
Let us note that we recover the familiar system of equations, allowing us to study the thermo-
dynamical properties of the Hamiltonian (1.1), by setting δ = 0. The function K(x) is defined
as
K(x) =
1
2π
∫ ∞
−∞
dk
e−ikx
1 + e|k|
=
1
4π
(
ψ
(
1− ix
2
)
− ψ
(
1 + ix
2
)
− ψ
(
1− ix
2
)
+ ψ
(
1 + i
x
2
))
,
where ψ is the digamma function. The auxiliary functions b(x, µ) and b(x, µ) are solutions of
a pair of non-linear integral equations given by
ln(b(x, µ1)) = − 2πJ
T cosh(πx)
− 2πδ
T cosh(π(x− µ1))
+
∫ ∞
−∞
dtK(x− t) ln(1 + b(t, µ1))
−
∫ ∞
−∞
dtK(x− t+ i) ln(1 + b(t, µ1)) (2.2)
and a similar equation obtained by exchanging b ↔ b and i↔ −i in (2.2). The function γ can
now be introduced as
γ(µ1, µ2) = −1 +
(
1 + (µ1 − µ2)2
)
T
∂
∂δ
ln(Λ(µ2, µ1))
∣∣∣
δ=0
. (2.3)
It has been conjectured [3] that the correlation functions of the isotropic Heisenberg chain at
any finite temperature (for vanishing magnetic field) are polynomials in γ and its derivatives
evaluated at (0, 0). A similar statement (involving a function ω and its derivative with respect
to the disorder parameter) was proved for the anisotropic XXZ chain [5, 8, 2]. Amazingly the
isotropic limit seems non-trivial and is still a subject of ongoing work. Here we would only like
to mention that the nearest- and next-to-nearest-neighbor two-point functions were expressed
in terms of γ in [3] starting from the multiple integral representation for the density matrix of
the Heisenberg chain obtained in [7]. The formulae for the longitudinal two-point functions are,
for instance,
〈σz
1σ
z
2〉T = −1
3
γ(0, 0), (2.4)
〈σz
1σ
z
3〉T = −1
3
γ(0, 0)− 1
6
γxx(0, 0) +
1
3
γxy(0, 0). (2.5)
They will be used below to test our results for the low-temperature expansion. We denoted
derivatives with respect to the first (resp. second) argument by the subscript x (resp. y). Similar
results for four sites can be obtained from [1] in the isotropic limit. In previous work [13] the
high-temperature expansion (up to order 25) of the two-point functions was obtained analytically
based on (2.4) and (2.5).
Universal Low Temperature Asymptotics of the Correlation Functions 3
3 Low-temperature expansion
To compute the low-temperature expansion of γ, we follow the line of reasoning of the article [9],
where a similar task was performed for the free energy. There are, however, two differences
between the usual equations and the ones used in this note: the additional driving term in (2.2)
proportional to δ and the shift µ2 in the kernel of the integration in (2.1).
The computation is based on the introduction of a shift L = 1
π ln
(
π J
T
)
in the auxiliary
functions:
bL(x) = b(x+ L) and b̃L(x) = b(−x− L).
In the low-temperature limit these functions satisfy
ln(bL(x, µ1)) ∼ −4e−πx − 4
δ
J
e−π(x−µ1) +DL(x)
+
∫ ∞
−L
dt
[
K(x− t) ln(1 + bL(t, µ1))−K(x− t+ i) ln(1 + bL(t, µ1))
]
, (3.1)
where DL(x) is the rest of the integral which does not contribute to the low-temperature limit,
when the magnetic field vanishes (see [9]). A similar relation holds with b ↔ b and i ↔ −i
exchanged.
In terms of the shifted functions the largest eigenvalue becomes
ln(Λ(µ2, µ1)) ∼
4πJ
T
K(µ2) +
4πδ
T
K(µ2 − µ1)
+
T
Jπ
∫ ∞
−L
dteπ(µ2−t) ln
[
(1 + bL(t, µ1))(1 + bL(t, µ1))
]
+
T
Jπ
∫ ∞
−L
dte−π(µ2+t) ln
[
(1 + b̃L(t, µ1))(1 + b̃L(t, µ1))
]
.
To evaluate these integrals we compute
I =
∫ ∞
−L
dt
[
ln(1 + bL(t, µ1)) ln(bL(t, µ1))′ + ln(1 + bL(t, µ1)) ln(bL(t, µ1))′
]
using two different methods. Here the prime stands for the derivative with respect to t. First,
we compute it explicitly using the change of variables z = ln(bL) or z = ln(bL), respectively,
which results in
I = 2
∫ 0
−∞
ln
(
1 + ez
)
dz =
π2
6
.
Second, we replace ln(bL(t, µ1)) and ln(bL(t, µ1)) by their scaling limits (3.1) and simplify the
resulting expression by taking into account that the derivative of K(x) is odd and contributions
by double integrals cancel pairwise. This way we obtain
I = 4π
(
1 +
δ
J
eπµ1
)∫ ∞
−L
dte−πt ln
[
(1 + b(t, µ1))(1 + b(t, µ1))
]
.
The same type of manipulation can be performed for the functions b̃, and a similar result is
obtained with µ1 replaced by −µ1.
Gathering these findings we obtain the asymptotic form of the largest eigenvalue,
ln(Λ(µ2, µ1)) ∼
4πJ
T
K(µ2) +
4πδ
T
K(µ2 − µ1) +
T
24J
(
eπµ2
1 + δ
J e
πµ1
+
e−πµ2
1 + δ
J e
−πµ1
)
.
4 N. Crampé, F. Göhmann and A. Klümper
0
0.05
0.1
0.15
0.2
0.25
1e-02 1e-01 1e+00 1e+01 1e+02
ne
xt
-t
o-
ne
ar
es
t n
ei
gh
bo
ur
z
z-
co
rr
el
at
io
n
fu
nc
tio
n
temperature T/J
NLIE
LTE
HTE
QMC
Figure 1. Comparison of the high- and low-temperature expansions (HTE, LTE) of 〈σz
1σ
z
3〉 with the
full numerical solution obtained from the integral equations (NLIE) and with Monte-Carlo data (QMC).
Thus, using (2.3), the function γ behaves asymptotically for small temperatures as
γ(µ1, µ2) ∼ −1 +
(
1 + (µ1 − µ2)2
)(
4πK(µ2 − µ1)−
T 2
12J2
cosh(π(µ1 + µ2))
)
.
This is our main result.
Using (2.4) and (2.5), we obtain the low-temperature expansion of the longitudinal correlation
functions
〈σz
1σ
z
2〉T ∼
1
3
− 4
3
ln(2) +
T 2
J2
1
36
,
〈σz
1σ
z
3〉T ∼
1
3
− 16
3
ln(2) + 3ζ(3)− T 2
J2
1
36
(
π2
2
− 4
)
.
The constant terms (independent of the temperature) in these expansions are in agreement
with those originally found in [11, 6]. In the figure we compare the combined low- and high-
temperature results for the next-to-nearest neighbor zz-correlation functions with the full nu-
merical curve obtained by implementing the linear and non-linear integral equations that deter-
mine γ and its derivatives [3] on a computer. The high-temperature data and some additional
Monte-Carlo data are taken from [14]. We find that the numerical curves (NLIE, QMC) are
amazingly well approximated by its low- and high-temperature approximations.
Acknowledgments
The authors are grateful to Jens Damerau for providing his computer program and to Christian
Trippe for producing the figure. They wish to express their gratitude to Z. Tsuboi and M. Shi-
roishi for providing their data. N.C. thanks the department of physics of Wuppertal University,
where this work was initiated, for hospitality.
Universal Low Temperature Asymptotics of the Correlation Functions 5
References
[1] Boos H., Damerau J., Göhmann F., Klümper A., Suzuki J., Weiße A., Short-distance thermal correlations
in the XXZ chain, J. Stat. Mech. Theory Exp. 2008 (2008), P08010, 23 pages, arXiv:0806.3953.
[2] Boos H., Göhmann F., On the physical part of the factorized correlation functions of the XXZ chain,
J. Phys. A: Math. Theor. 42 (2009), 315001, 27 pages, arXiv:0903.5043.
[3] Boos H., Göhmann F., Klümper A., Suzuki J., Factorization of multiple integrals representing the density
matrix of a finite segment of the Heisenberg spin chain, J. Stat. Mech. Theory Exp. 2006 (2006), P04001,
13 pages, hep-th/0603064.
[4] Boos H., Göhmann F., Klümper A., Suzuki J., Factorization of the finite temperature correlation functions
of the XXZ chain in a magnetic field, J. Phys. A: Math. Theor. 40 (2007), 10699–10728, arXiv:0705.2716.
[5] Boos H., Jimbo M., Miwa T., Smirnov F., Takeyama Y., Hidden Grassmann structure in the XXZ model.
II. Creation operators, Comm. Math. Phys. 286 (2009), 875–932, arXiv:0801.1176.
[6] Boos H.E., Korepin V.E., Quantum spin chains and Riemann zeta function with odd arguments, J. Phys. A:
Math. Gen. 34 (2001), 5311–5316, hep-th/0104008.
[7] Göhmann F., Klümper A., Seel A., Integral representation of the density matrix of the XXZ chain at finite
temperature, J. Phys. A: Math. Gen. 38 (2005), 1833–1841, cond-mat/0412062.
[8] Jimbo M., Miwa T., Smirnov F., Hidden Grassmann structure in the XXZ model. III. Introducing Matsubara
direction, J. Phys. A: Math. Theor. 42 (2009), 304018, 31 pages, arXiv:0811.0439.
[9] Klümper A., The spin-1/2 Heisenberg chain: thermodynamics, quantum criticality and spin-Peierls expo-
nents, Eur. Phys. J. B 5 (1998), 677–685, cond-mat/9803225.
[10] Suzuki M., Transfer-matrix method and Monte-Carlo simulation in quantum spin systems, Phys. Rev. B 31
(1985), 2957–2965.
[11] Takahashi M., Half-filled Hubbard model at low temperature, J. Phys. C 10 (1977), 1289–1301.
[12] Trippe C., Göhmann F., Klümper A., Short-distance thermal correlations in the massive XXZ chain, Eur.
Phys. J. B 73 (2010), 253–264, arXiv:0908.2232.
[13] Tsuboi Z., A note on the high temperature expansion of the density matrix for the isotropic Heisenberg
chain, Phys. A 377 (2007), 95–101, cond-mat/0611454.
[14] Tsuboi Z., Shiroishi M., High temperature expansion of the emptiness formation probability for the isotropic
Heisenberg chain, J. Phys. A: Math. Gen. 38 (2005), L363–L370, cond-mat/0502569.
http://dx.doi.org/10.1088/1742-5468/2008/08/P08010
http://arxiv.org/abs/0806.3953
http://dx.doi.org/10.1088/1751-8113/42/31/315001
http://arxiv.org/abs/0903.5043
http://dx.doi.org/10.1088/1742-5468/2006/04/P04001
http://arxiv.org/abs/hep-th/0603064
http://dx.doi.org/10.1088/1751-8113/40/35/001
http://arxiv.org/abs/0705.2716
http://dx.doi.org/10.1007/s00220-008-0617-z
http://arxiv.org/abs/0801.1176
http://dx.doi.org/10.1088/0305-4470/34/26/301
http://dx.doi.org/10.1088/0305-4470/34/26/301
http://arxiv.org/abs/hep-th/0104008
http://dx.doi.org/10.1088/0305-4470/38/9/001
http://arxiv.org/abs/cond-mat/0412062
http://dx.doi.org/10.1088/1751-8113/42/30/304018
http://arxiv.org/abs/0811.0439
http://dx.doi.org/10.1007/s100510050491
http://arxiv.org/abs/cond-mat/9803225
http://dx.doi.org/10.1103/PhysRevB.31.2957
http://dx.doi.org/10.1088/0022-3719/10/8/031
http://dx.doi.org/10.1140/epjb/e2009-00417-7
http://dx.doi.org/10.1140/epjb/e2009-00417-7
http://arxiv.org/abs/0908.2232
http://dx.doi.org/10.1016/j.physa.2006.11.018
http://arxiv.org/abs/cond-mat/0611454
http://dx.doi.org/10.1088/0305-4470/38/20/L05
http://arxiv.org/abs/cond-mat/0502569
1 Introduction
2 Definition of the basic function
3 Low-temperature expansion
References
|
| id | nasplib_isofts_kiev_ua-123456789-146513 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1815-0659 |
| language | English |
| last_indexed | 2025-12-02T07:01:51Z |
| publishDate | 2010 |
| publisher | Інститут математики НАН України |
| record_format | dspace |
| spelling | Crampé, N. Göhmann, F. Klümper, A. 2019-02-09T19:40:09Z 2019-02-09T19:40:09Z 2010 Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain / N. Crampé, F. Göhmann, A. Klümper // Symmetry, Integrability and Geometry: Methods and Applications. — 2010. — Т. 6. — Бібліогр.: 14 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 81Q80; 82B23 DOI:10.3842/SIGMA.2010.082 https://nasplib.isofts.kiev.ua/handle/123456789/146513 We calculate the low temperature asymptotics of a function γ that generates the temperature dependence of all static correlation functions of the isotropic Heisenberg chain. This paper is a contribution to the Proceedings of the International Workshop “Recent Advances in Quantum Integrable Systems”. The full collection is available at http://www.emis.de/journals/SIGMA/RAQIS2010.html. The authors are grateful to Jens Damerau for providing his computer program and to Christian Trippe for producing the figure. They wish to express their gratitude to Z. Tsuboi and M. Shiroishi for providing their data. N.C. thanks to the department of physics of Wuppertal University, where this work was initiated, for hospitality. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain Article published earlier |
| spellingShingle | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain Crampé, N. Göhmann, F. Klümper, A. |
| title | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain |
| title_full | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain |
| title_fullStr | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain |
| title_full_unstemmed | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain |
| title_short | Universal Low Temperature Asymptotics of the Correlation Functions of the Heisenberg Chain |
| title_sort | universal low temperature asymptotics of the correlation functions of the heisenberg chain |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/146513 |
| work_keys_str_mv | AT crampen universallowtemperatureasymptoticsofthecorrelationfunctionsoftheheisenbergchain AT gohmannf universallowtemperatureasymptoticsofthecorrelationfunctionsoftheheisenbergchain AT klumpera universallowtemperatureasymptoticsofthecorrelationfunctionsoftheheisenbergchain |