A Critical Phenomenon in Solitonic Ising Chains
We discuss a phase transition of the second order taking place in non-local 1D Ising chains generated by specific infinite soliton solutions of the KdV and BKP equations.
Gespeichert in:
| Datum: | 2007 |
|---|---|
| Hauptverfasser: | , |
| Format: | Artikel |
| Sprache: | English |
| Veröffentlicht: |
Інститут математики НАН України
2007
|
| Schriftenreihe: | Symmetry, Integrability and Geometry: Methods and Applications |
| Online Zugang: | https://nasplib.isofts.kiev.ua/handle/123456789/147379 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Zitieren: | A Critical Phenomenon in Solitonic Ising Chains / I.M. Loutsenko, V.P. Spiridonov // Symmetry, Integrability and Geometry: Methods and Applications. — 2007. — Т. 3. — Бібліогр.: 13 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
nasplib_isofts_kiev_ua-123456789-147379 |
|---|---|
| record_format |
dspace |
| spelling |
nasplib_isofts_kiev_ua-123456789-1473792025-02-09T16:35:10Z A Critical Phenomenon in Solitonic Ising Chains Loutsenko, I.M. Spiridonov, V.P. We discuss a phase transition of the second order taking place in non-local 1D Ising chains generated by specific infinite soliton solutions of the KdV and BKP equations. The work of I.M. Loutsenko has been supported by European Community grant MIFI-CT-2005-007323 and V.P. Spiridonov is partially supported by the Russian Foundation for Basic Research (grant no. 05-01-01086). The authors are grateful to V.B. Priezzhev for useful remarks. 2007 Article A Critical Phenomenon in Solitonic Ising Chains / I.M. Loutsenko, V.P. Spiridonov // Symmetry, Integrability and Geometry: Methods and Applications. — 2007. — Т. 3. — Бібліогр.: 13 назв. — англ. 1815-0659 2000 Mathematics Subject Classification: 70H06; 82B20 https://nasplib.isofts.kiev.ua/handle/123456789/147379 en Symmetry, Integrability and Geometry: Methods and Applications application/pdf Інститут математики НАН України |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| language |
English |
| description |
We discuss a phase transition of the second order taking place in non-local 1D Ising chains generated by specific infinite soliton solutions of the KdV and BKP equations. |
| format |
Article |
| author |
Loutsenko, I.M. Spiridonov, V.P. |
| spellingShingle |
Loutsenko, I.M. Spiridonov, V.P. A Critical Phenomenon in Solitonic Ising Chains Symmetry, Integrability and Geometry: Methods and Applications |
| author_facet |
Loutsenko, I.M. Spiridonov, V.P. |
| author_sort |
Loutsenko, I.M. |
| title |
A Critical Phenomenon in Solitonic Ising Chains |
| title_short |
A Critical Phenomenon in Solitonic Ising Chains |
| title_full |
A Critical Phenomenon in Solitonic Ising Chains |
| title_fullStr |
A Critical Phenomenon in Solitonic Ising Chains |
| title_full_unstemmed |
A Critical Phenomenon in Solitonic Ising Chains |
| title_sort |
critical phenomenon in solitonic ising chains |
| publisher |
Інститут математики НАН України |
| publishDate |
2007 |
| url |
https://nasplib.isofts.kiev.ua/handle/123456789/147379 |
| citation_txt |
A Critical Phenomenon in Solitonic Ising Chains / I.M. Loutsenko, V.P. Spiridonov // Symmetry, Integrability and Geometry: Methods and Applications. — 2007. — Т. 3. — Бібліогр.: 13 назв. — англ. |
| series |
Symmetry, Integrability and Geometry: Methods and Applications |
| work_keys_str_mv |
AT loutsenkoim acriticalphenomenoninsolitonicisingchains AT spiridonovvp acriticalphenomenoninsolitonicisingchains AT loutsenkoim criticalphenomenoninsolitonicisingchains AT spiridonovvp criticalphenomenoninsolitonicisingchains |
| first_indexed |
2025-11-28T00:02:20Z |
| last_indexed |
2025-11-28T00:02:20Z |
| _version_ |
1849990204888186880 |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 3 (2007), 059, 11 pages
A Critical Phenomenon in Solitonic Ising Chains?
Igor M. LOUTSENKO † and Vyacheslav P. SPIRIDONOV ‡
† Mathematical Institute, University of Oxford, 24-29 St. Gilles’, Oxford, OX1 3LB, UK
E-mail: loutsenk@maths.ox.ac.uk
‡ Bogoliubov Laboratory of Theoretical Physics, JINR, Dubna, Moscow Region, 141980 Russia
E-mail: spiridon@theor.jinr.ru
Received December 04, 2006, in final form April 17, 2007; Published online April 24, 2007
Original article is available at http://www.emis.de/journals/SIGMA/2007/059/
Abstract. We discuss a phase transition of the second order taking place in non-local 1D
Ising chains generated by specific infinite soliton solutions of the KdV and BKP equations.
Key words: Ising chain; solitons; phase transition
2000 Mathematics Subject Classification: 70H06; 82B20
To the memory of Vadim B. Kuznetsov
1 The Korteweg–de Vries solitonic spin chain
In a series of papers [9, 10], we described a direct relation between soliton solutions of integrable
hierarchies and lattice gas systems (e.g., Coulomb gases on two dimensional lattices). The latter
models can be reformulated also as some Ising spin systems with a non-local exchange. In
particular, the grand partition functions of specific N -site Ising chains for some fixed values
of the temperature were shown to coincide with the tau-functions of N -soliton solutions of the
Korteweg–de Vries (KdV) and Kadomtsev–Petviashvili (KP) equations.
We would like to complete here the consideration of [9] and investigate a critical phenomenon
appearing in these models in the zero temperature limit. The N -soliton solution of the KdV
equation ut + uxxx − 6uux = 0 has the form [1, 11]
u(x, t) = −2∂2
x log τN (x, t),
where τN is the determinant of a N ×N matrix M ,
τN = detM, Mij = δij +
2
√
kikj
ki + kj
e(θi+θj)/2,
θi = kix− k3
i t + θ
(0)
i , i, j = 1, 2, . . . , N.
The parameters ki describe amplitudes of solitons, θ
(0)
i /ki are the zero time phases of solitons,
and k2
i are their velocities. The tau-function τN admits the following Hirota type representa-
tion [7]:
τN =
∑
µi=0,1
exp
∑
1≤i<j≤N
Aijµiµj +
N∑
i=1
θiµi
, (1)
?This paper is a contribution to the Vadim Kuznetsov Memorial Issue ‘Integrable Systems and Related Topics’.
The full collection is available at http://www.emis.de/journals/SIGMA/kuznetsov.html
mailto:loutsenk@maths.ox.ac.uk
mailto:spiridon@theor.jinr.ru
http://www.emis.de/journals/SIGMA/2007/059/
http://www.emis.de/journals/SIGMA/kuznetsov.html
2 I.M. Loutsenko and V.P. Spiridonov
where the soliton phase shifts Aij are expressed in terms of the spectral variables ki as
eAij =
(ki − kj)2
(ki + kj)2
.
As remarked in [9] for θi = θ(0), this τN defines the grand partition function of a lattice gas
model with the chemical potential θ(0) and µi being the filling factors of the lattice sites by
molecules. The constants Aij describe interaction energy of the molecules.
Substituting in (1) µi = (σi + 1)/2, where σi = ±1 are other discrete variables, one can pass
from the lattice gases to Ising spin chains [3]:
τN = eΦZN , Φ =
1
4
∑
i<j
Aij +
1
2
N∑
j=1
θj ,
where
ZN =
∑
σi=±1
e−βE , β =
1
kT
, (2)
E =
∑
1≤i<j≤N
Jijσiσj −
N∑
i=1
Hiσi. (3)
Here, Jij are the exchange constants, Hi is an external magnetic field, T is the temperature,
and k is the Boltzmann constant:
βJij = −1
4
Aij , βHi =
1
2
θi +
1
4
N∑
j=1,i6=j
Aij .
τN coincides thus with the partition function of a one-dimensional Ising chain with the specific
non-local exchange. A similar situation holds for the KP hierarchy and some other partial
differential or difference nonlinear integrable equations.
From the thermodynamic point of view, it is interesting to understand the N →∞ behaviour
of these “solitonic” statistical mechanics models. In general there are infinitely many free
parameters, and it is difficult to classify qualitatively different cases. An interesting class of
models is related to the so-called self-similar potentials [12]. These potentials are characterized
by the q-periodicity constraints kj+M = qkj and θ
(0)
j+M = θ
(0)
j , where q is an arbitrary parameter,
0 < q < 1, and M is a positive integer. Exchange constants satisfy in this case the constraints
Ji+M,j+M = Jij . For M = 1, this translational invariance takes the simplest form Jij = J(i− j)
with ki forming one geometric progression, ki = k1q
i−1, q = e−2α, where k1 and α > 0 are free
parameters. More precisely,
Aij = 2 log | tanh α(i− j)|.
The KdV coordinate x and time t describe a part of the magnetic field Hi decaying exponen-
tially fast for i →∞ because q < 1. Only the values of constants θ
(0)
i are therefore relevant for
the leading asymptotics of the partition function in the N →∞ thermodynamic limit. Neglect-
ing this (x, t)-dependence, we come to the constraint Hi+M = Hi, which is just the homogeneity
condition for M = 1. In principle, it is possible to compute the N → ∞ asymptotics for the
partition function for arbitrary M -periodic magnetic fields, but in [9] only the M = 1, 2 cases
were considered.
Since 0 < | tanh α(i− j)| < 1, we have Jij = −Aij/4β > 0, which corresponds to an antifer-
romagnetic Ising chain (a similar picture holds for M > 1). Although we have a long distance
A Critical Phenomenon in Solitonic Ising Chains 3
interaction, its intensity falls off exponentially fast, and the absence of phase transitions in such
systems at non-zero temperatures is well known [5]. It appears, however, that there exists a spe-
cial limit leading to a nontrivial critical phenomenon in this model. Indeed, we consider the limit
α → 0+ or q → 1−. The phase shifts Aij ∝ Jij/kT then diverge. We can take nevertheless as
the true exchange constants J ren
ij = Jij(q−1−q) and as the true temperature kTren = kT (q−1−q).
For the self-consistency of the temperature definition, we should renormalize the magnetic field
as well, H ∝ h/(q−1 − q), and assume that h is finite.
As a result, the interaction energy of any spin “in the bulk” with all others,
Ei =
∞∑
j=−∞, 6=i
J ren
ij σj ,
is finite for q → 1− (or α → 0+). Indeed, the maximal value of this interaction energy is
Emax = lim
α→0
4α
∞∑
j=−∞, 6=0
J0j = − 2
β
lim
α→0
α
∞∑
j=−∞, 6=0
log
∣∣∣∣1− qj
1 + qj
∣∣∣∣ = − 4
β
lim
α→0
α log
(q; q)∞
(−q; q)∞
.
We use the notation (a; q)∞ =
∞∏
k=0
(1− aqk) and (a1, . . . , am; q)∞ =
m∏
j=1
(aj ; q)∞.
For the following considerations, we need theta functions [2]
θ1(ν, q) = −i
∑
n∈Z
(−1)nq(n+1/2)2e(2n+1)iν = iq1/4e−iν
(
q2, e2iν , q2e−2iν ; q2
)
∞,
θ2(ν, q) =
∑
n∈Z
q(n+1/2)2e(2n+1)iν = q1/4e−iν
(
q2,−e2iν ,−q2e−2iν ; q2
)
∞,
θ3(ν, q) =
∑
n∈Z
qn2
e2niν =
(
q2,−qe2iν ,−qe−2iν ; q2
)
∞,
θ4(ν, q) =
∑
n∈Z
(−1)nqn2
e2niν =
(
q2, qe2iν , qe−2iν ; q2
)
∞,
where q = eπiτ , Im(τ) > 0, and their modular transformations
θ1(ν/τ, q̃) = −i
√
−iτeiν2/πτθ1(ν, q), θ2(ν/τ, q̃) =
√
−iτeiν2/πτθ4(ν, q),
θ3(ν/τ, q̃) =
√
−iτeiν2/πτθ3(ν, q), θ4(ν/τ, q̃) =
√
−iτeiν2/πτθ2(ν, q),
where q̃ = e−πi/τ and
√
−iτ is positive for purely imaginary τ . Using these formulas, we obtain
(q; q)2∞
(−q; q)2∞
=
θ′1(0, q1/2)
θ2(0, q1/2)
=
θ′1(0, q̃1/2)
(−iτ)θ4(0, q̃1/2)
= q̃1/8 2
(−iτ)
(q̃; q̃)2∞
(q̃1/2; q̃)2∞
,
where q = e2πiτ = e−2α and q̃ = e−2πi/τ = e−2π2/α. As a result,
Emax =
π2
2β
< ∞,
and Ei ≤ Emax. The limit α → 0 corresponds thus to an infinitely small and infinitely long-range
nonlocal interaction model at a low value (zero) effective temperature.
There are some other interesting limits. For instance, for q → 0 and finite h, we obtain the
high temperature nearest neighbor interaction spin chain, J ren
ij ∝ δi+1,j , Tren →∞. For finite H,
this limit corresponds to the non-interacting spins. The solitonic interpretation describes thus
4 I.M. Loutsenko and V.P. Spiridonov
only a two-dimensional subspace of parameters (T,H, q). For fixed q, the temperature T is also
fixed, and we can set the “KdV temperature” equal to β = 1.
Using the Wronskian representation for τN , the leading asymptotics of ZN for N → ∞ was
determined in [9] for the M = 1 translationally invariant model and a homogeneous magnetic
field. Namely, ZN → exp(−NβfI), where the free energy per site fI has the form
−βfI(q, H) = log
2(q4; q4)∞ coshβH
(q2; q2)1/2
∞
+
1
4π
∫ 2π
0
dν log
(
|ρ(ν)|2 − q tanh2 βH
)
,
|ρ(ν)|2 =
(q2eiν , q2e−iν ; q4)2∞
(q4eiν , q4e−iν ; q4)2∞
1
4 sin2(ν/2)
= q
θ2
4(ν/2, q2)
θ2
1(ν/2, q2)
.
The total magnetization of the lattice takes the form:
m(H) = −∂fI
∂H
= lim
N→∞
1
N
N∑
i=1
〈σi〉 =
(
1− 1
π
∫ π
0
dν
1 + d(ν) cosh2 βH
)
tanh βH, (4)
where
d(ν) =
θ2
4(ν, q2)
θ2
1(ν, q2)
− 1.
The function m(H) grows monotonically with H and reflects qualitative predictions of the
general theory of 1D systems with the fast decaying interactions [5]. However, the limit α → 0
with the renormalized exchange and magnetic field breaks down the corresponding necessary
conditions, and we obtain a non-trivial critical phenomenon.
We substitute in (4) βH = h/(q−1 − q), h > 0, and take the limit α → 0. Since
θ2
4(ν, q2)
θ2
1(ν, q2)
= −θ2
2(ν/τ, q̃2)
θ2
1(ν/τ, q̃2)
=
(−e2iν/τ ,−q̃4e−2iν/τ ; q̃4)2∞
(e2iν/τ , q̃4e−2iν/τ ; q̃4)2∞
,
where q = eπiτ/2 = e−2α, q̃ = e−πi/2τ = q−π2/8α, we have
χ(ν) ≡ lim
α→0
d(ν) cosh2 h
4α
= lim
α→0
(
(1 + e2iν/τ )2(1 + e2i(ν−π)/τ )2
(1− e2iν/τ )2(1− e2i(ν−π)/τ )2
− 1
)
eh/2α
4
= lim
α→0
(
e−2iν/τ + e2i(ν−π)/τ
)
eh/2α, 0 < ν < π.
Substituting this result in (4) and using relation lim
α→0
tanh(h/4α) = 1, we obtain
m(h) = 1− 2
π
∫ π/2
0
dν
1 + χ(ν)
.
For 0 < ν < π/2, we have
χ(ν) =
{
0, if h/π < ν < π/2,
∞, if 0 < ν < h/π,
for h < π2/2, and χ(ν) = ∞, for h > π2/2. The final result therefore can be represented in the
form
m(h) =
1− 2
π
∫ π/2
h/π
dν =
2
π2
h, if |h| < π2/2,
1, if |h| ≥ π2/2.
A Critical Phenomenon in Solitonic Ising Chains 5
We have an obvious point of non-analyticity of m(h) or of the free energy fI(h) at the cri-
tical value of the magnetic field hcrit = π2/2, such that the magnetic susceptibility χ(h) =
β−1dm(h)/dh has a jump at it (i.e., we have the phase transition of the second order). This is
a typical phenomenon in the systems with long-range interaction, where the mean field approxi-
mation gives exact values for the one-point correlation functions (see, e.g., [3]).
2 The mean field approximation
In the mean field theory, one considers a few degrees of freedom (usually, just one) of a taken
system in an effective mean field of the remaining part of the system. This effective or mean
field depends itself on the analysis of the one-body dynamics. As an example, we consider the
general spin chain with energy (3) and the mean magnetization at the i-th site of the lattice
〈σi〉 =
∑
σ1,σ2,...
σie
−βE∑
σ1,σ2,...
e−βE
.
Instead of calculating the above sums, we stick to the i-th spin and evaluate its contribution to
the energy as
Ei(σi) = −σiH̃i, (5)
where
H̃i = −
∑
j 6=i
Jij〈σj〉+ Hi (6)
is an effective mean magnetic field at the i-th site created by the external field Hi and the rest
of the system ∝
∑
j 6=i
Jij〈σj〉.
In the one body problem (5), the configuration space consists of two states σi = ±1, and
therefore
〈σi〉 =
eβH̃i − e−βH̃i
eβH̃i + e−βH̃i
= tanh(βH̃i).
Substituting the values of effective fields (6) in the last equation, we obtain a system of tran-
scendental equations for mean values of all spins
〈σi〉 = tanhβ
−∑
j 6=i
Jij〈σj〉+ Hi
. (7)
We consider now the translationally invariant system Jij = J(i− j), Hi = H. In such a system,
all mean values of the spins are the same in the thermodynamic limit, 〈σi〉 = 〈σ〉, and from (7),
we obtain
〈σ〉 = tanhβ
−∑
j 6=i
J(i− j)〈σ〉+ H
= tanhβ(−J〈σ〉+ H), J =
∑
j 6=0
J(j).
The solution of this equation is an intersection of graphs of two functions: y = x and y =
tanh β(−Jx + H). We consider two different cases.
6 I.M. Loutsenko and V.P. Spiridonov
1) A ferromagnet in the zero magnetic field: J < 0, H = 0. Our system of equations
has only the trivial solution x = y = 0 for β|J | ≤ 1 and three solutions x = y = 0 and
x = y = ±m for β|J | > 1 and some 0 < m < 1. The latter nontrivial solutions describe the
spontaneous magnetization m at the temperatures smaller than the critical value 1/βcrit = |J |
(the tanh(−βJx)-function becomes steeper at the origin as β increases and starts to intersect
the line x = y in two additional points as its slope exceeds the critical value).
2) An antiferromagnet, J > 0. The only solution at H = 0 is the trivial solution x = y = 0.
We consider now the zero-temperature limit β →∞. In this case, the tanh-function transforms
to the sign-function, and we obtain{
y = x,
y = sgn(−Jx + H),
sgn(x) =
{
1, x > 0,
−1, x < 0.
Solving these equations is rather easy. The function sgn(x) is shifted by H/J from the origin
along the x-axis. When |H/J | < 1, it intersects with the line y = x by its vertical part, and the
magnetization equals to H/J . When |H| exceeds J , the line y = x intersects with one of the
horizontal branches of sgn(x), and the magnetization becomes equal to ±1. The magnetization
is thus a continuous piecewise linear function of H consisting of three parts: two constant
〈σ〉 = ±1 for |H| > J and the linear piece 〈σ〉 = H/J connecting them through the origin. In
our KdV-solitonic model, we denoted βH = h/(q−1 − q) and Emax = (q−1 − q)J . In the limit
q → 1, we have Emax = π2/2β and 〈σ〉 = H/J = 2h/π2, which leads to the exact value of
the critical magnetic field hcrit = π2/2. Such a qualitative behaviour of the system is obvious:
when the external field exceeds the interaction energy between spins, they all flip in the field
direction.
The mean field approximation is known to give exact one-point correlation functions (e.g.,
the magnetization) for systems with the long-range interaction (as our q → 1 limit). It might
be non-suitable, however, for the two point correlators (e.g., 〈σiσj〉).
We consider now the M -periodic chain. For general (not necessarily solitonic) M -periodic
chain, it is reasonable to introduce multi-index exchange Jnm(i− j), where i− j is the distance
between the cells and 1 ≤ n, m ≤ M are the respective internal cell indices for the n-th and
m-th sublattices. In our particular solitonic KdV case, we have
βJnm(i) = −1
2
log
∣∣∣∣knqi − km
knqi + km
∣∣∣∣ , Jnm(i) = Jmn(−i), (8)
and the energy
E =
∑
i,j∈Z, n6=m
Jnm(i− j)σ(n)
i σ
(m)
j +
∑
i,j∈Z, i6=j, n
Jnn(i− j)σ(n)
i σ
(n)
j +
∑
i∈Z, n
H(n)σ
(n)
i .
In the latter sum the magnetic field is also M -periodic, H
(n)
i = H(n), inhomogeneous only inside
the cells. The analysis similar to the M = 1 case yields from (7) the following system of M
equations for M unknowns 〈σ(n)〉:
〈σ(n)〉 = tanhβ
(
−
∑
m
Jnm〈σ(m)〉+ H(n)
)
, n,m = 1, . . . ,M.
In these equations,
Jnn =
∑
i∈Z,6=0
Jnn(i), Jnm =
∑
i∈Z
Jnm(i), n 6= m.
A Critical Phenomenon in Solitonic Ising Chains 7
In the zero-temperature limit, we have
〈σ(n)〉 = sgn
(
−
∑
m
Jnm〈σ(m)〉+ H(n)
)
, n,m = 1, . . . ,M, (9)
where sgn(x) is the sign-function.
In the solitonic case (8), the M×M matrix J = Jnm is symmetric with the constant diagonal:
Jnn = − 1
2β
∑
i∈Z, i6=0
log
∣∣∣∣1− qi
1 + qi
∣∣∣∣ ≡ −A, Jnm = Jmn = − 1
2β
∑
i∈Z
log
∣∣∣∣knqi − km
knqi + km
∣∣∣∣ . (10)
For M = 2, we have J12 = J21 ≡ −B and
〈σ(1)〉 = sgn(A〈σ(1)〉+ B〈σ(2)〉+ H(1)),
〈σ(2)〉 = sgn(B〈σ(1)〉+ A〈σ(2)〉+ H(2)).
If we take the uniform magnetic field H(1) = H(2) = H, these equations become symmetric
in σ(1) and σ(2) and have the solution 〈σ〉 = 〈σ(1)〉 = 〈σ(2)〉 stemming from one equation
〈σ〉 = sgn((A + B)〈σ〉+ H).
For the completely uniform magnetic field, there exists thus a solution when the spins flip
simultaneously for both sublattices for sufficiently large magnetic fields.
Simultaneous phase transition exists for all sublattices in the uniform field H(1) = · · · =
H(M) = H, when 〈σ(1)〉 = · · · = 〈σ(M)〉 and all spins in all sublattices are aligned simultaneously
for a sufficiently large H. As seen from (9), such a solution exists, if
∑
m Jnm are equal, which
is certainly true for M = 2 because of the permutational symmetry. But it may be not so for
M > 2. For instance, in the solitonic case (10) for M = 3, we have∑
m
J1m ∝
∑
i∈Z,6=0
log
∣∣∣∣1− qi
1 + qi
∣∣∣∣+∑
i∈Z
log
∣∣∣∣k1q
i − k2
k1qi + k2
∣∣∣∣+∑
i∈Z
log
∣∣∣∣k1q
i − k3
k1qi + k3
∣∣∣∣ ,
∑
m
J2m ∝
∑
i∈Z,6=0
log
∣∣∣∣1− qi
1 + qi
∣∣∣∣+∑
i∈Z
log
∣∣∣∣k2q
i − k1
k2qi + k1
∣∣∣∣+∑
i∈Z
log
∣∣∣∣k2q
i − k3
k2qi + k3
∣∣∣∣ ,
∑
m
J3m ∝
∑
i∈Z,6=0
log
∣∣∣∣1− qi
1 + qi
∣∣∣∣+∑
i∈Z
log
∣∣∣∣k3q
i − k1
k3qi + k1
∣∣∣∣+∑
i∈Z
log
∣∣∣∣k3q
i − k2
k3qi + k2
∣∣∣∣ .
These three sums are certainly different for 0 < q < 1, and in the limit β → ∞ (which we
cannot reach within the solitonic interpretation for q < 1), the magnetization would become
a piecewise linear function of H of a more complicated form than in the M = 1, 2 cases. In our
model, however, the zero temperature is reached by multiplication of the above sums by q−1− q
and taking the limit q → 1− (or α → 0+). All three sums become then equal yielding the same
magnetization as in the M = 1 and M = 2 cases.
3 The BKP solitonic spin chain
Another Ising chain model solved in [9] appears from the multisoliton solution of the KP equation
of B type, i.e. the BKP equation [4]. The corresponding partition function has the same form (2),
where the exchange constants are
βJij = −1
4
Aij , eAij =
(ai − aj)(bi − bj)(ai − bj)(bi − aj)
(ai + aj)(bi + bj)(ai + bj)(bi + aj)
.
8 I.M. Loutsenko and V.P. Spiridonov
For ai = bi = ki/2, this model coincides with the KdV-inspired model at the twice lower value
of the temperature obtained after the change β → 2β.
The translational invariance of this spin chain, Jij = J(i− j), yields
ai = qi−1, bi = bqi−1, q = e−2α,
where we normalize a1 = 1 and assume that 0 < q < 1 as before. This gives the exchange
βJij = −1
4
log
tanh2 α(i− j)− (b− 1)2/(b + 1)2
coth2 α(i− j)− (b− 1)2/(b + 1)2
,
where the parameter b is restricted to three regions (because of the b → 1/b invariance): either
−1 < b < −q (the ferromagnetic chain, Jij < 0), or q < b ≤ 1 or |b| = 1, b 6= −1 (the
antiferromagnetic chain, Jij > 0).
In the thermodynamic limit N →∞, the free energy per site for the homogeneous magnetic
field Hi = H takes the form [9]:
−βfI(H) =
1
4
log
(q, q, bq, q/b; q)∞
(−q,−q,−bq,−q/b; q)∞
+
1
4π
∫ 2π
0
dν log |2ρ(ν)|,
where
ρ(ν) = cosh 2βH +
(−q; q)2∞
(−eiν ,−qe−iν ; q)∞
(
(b−1eiν , qbe−iν ; q)∞
(b−1, qb; q)∞
+
(beiν , qb−1e−iν ; q)∞
(b, qb−1; q)∞
)
.
Taking the derivative with respect to H, we find the magnetization
m(H) =
(
1− 1
π
∫ π
0
dν
1 + d(ν) cosh 2βH
)
tanh 2βH, (11)
where
d(ν) = 2
θ1(φ/2, q1/2)
θ2(0, q1/2)
θ2(ν, q1/2)
θ1(ν + φ/2, q1/2)− θ1(ν − φ/2, q1/2)
(12)
with b = eiφ. For real φ we have |b| = 1, the choice φ = iγ, 0 < γ < 2α, yields q < b < 1, and
for φ = π + iγ, we have −1 < b < −q. The limit b → 1 describes the magnetization for the
“KdV-spin chain” at the twice lower value of the temperature. A simple test of this expression
consists in the choice b = −1 corresponding to the non-interacting spins, Jij = 0. In this case
d(ν) = 1, and we obtain m(H) = tanhβH as it should be for the free system.
We substitute now q = e−2α and H = 2h/(q−1 − q) in (11) and consider the limit α → 0+.
The factor 2 in front of h was chosen for coincidence of this model with the the KdV spin chain
with the effective replacement β → 2β (i.e., the twice lower value of the temperature). We apply
the modular transformation to theta functions in (12) and obtain
d(ν) = 2
θ1(φ/2τ, q̃1/2)
θ4(0, q̃1/2)
θ2(ν/τ, q̃1/2)
e−i νφ
πτ θ1
(
ν+φ/2
τ , q̃1/2
)
− ei νφ
πτ θ1
(
ν−φ/2
τ , q̃1/2
) ,
where τ = iα/π and q̃ = e−2πi/τ = e−2π2/α. Denoting
χ(ν) = lim
α→0
d(ν) cosh
4qh
1− q2
,
A Critical Phenomenon in Solitonic Ising Chains 9
we therefore obtain
χ(ν) = lim
α→0
1
2e−
πφ
2α
(
1− e
πφ
α
)(
1− e
π(2ν−π)
α
)
e
h
α
e−
πν
α
(1+ φ
π
)−πφ
2α
(
1− e2π
ν+φ/2
α
)
− e−
πν
α
(1−φ
π
)+πφ
2α
(
1− e2π
ν−φ/2
α
) .
In the region φ = iγ, 0 < γ < 2α, we have
χ(ν) = lim
α→0
sin(πγ/2α)
(
1− e
π
α
(2ν−π)
)
e
h
α
−πν
α
2 sin((π/2− ν)γ/α)
.
Since 0 < γ/α < 2, the sin-factors do not influence the asymptotic behaviour, and for 0 < ν <
π/2, we find χ(ν) = 0 for h/π < ν < π/2 and χ(ν) = ∞ for ν < h/π. As a result, we obtain
m(h) = 2h/π2 for |h| < π2/2 and m(h) = 1 for |h| ≥ π2/2. This is the same picture as for the
“KdV-chain”, as it should be because the limit α → 0 assumes the limit γ → 0 or b → 1.
In a similar way, for φ = π + iγ, 0 ≤ γ < 2α, and α → 0, we find b → −1, i.e. the trivial
situation of free spins. The most interesting behaviour appears in the region 0 < φ < π, for
which we find
m(h) = 1− 2
π
∫ π/2
0
dν
1 + χ(ν)
,
where
χ(ν) =
{
0, if h/(π − φ) < ν,
∞, if ν < h/(π − φ).
As a result,
m(h) =
2h
π(π − φ)
, if |h| < π(π − φ)
2
,
and m(h) = 1, if |h| ≥ π(π − φ)/2. The critical value of the magnetic field, for which we have
the phase transition,
hcrit = ±π(π − φ)
2
,
depends explicitly on the parameter of the model φ, and for φ = 0, we obtain the previous result.
The Ising spin systems associated with the multisoliton solutions of integrable nonlinear equa-
tions provide thus the models with phase transitions already in their simplest one-dimensional
spin chain realizations.
It is interesting to analyze consequences of the antiferromagnetic nature of the exchange. We
take for this the KdV-inspired Ising chain and apply different magnetic fields to the odd, H1,
and even, H2, sites. The corresponding magnetization for the odd sites sublattice was derived
in [9]:
modd(H1,H2) = −2
dfI
dH1
= lim
p→∞
1
p
p∑
j=1
〈σ2j−1〉
= tanhβH1 −
tanh βH2
cosh2 βH1
1
π
∫ π
0
dν
θ2
4(ν,q2)
θ2
1(ν,q2)
− tanh βH1 tanh βH2
.
10 I.M. Loutsenko and V.P. Spiridonov
The magnetization for the even sites sublattice is obtained after permuting H1 and H2 in this
expression. Obviously, if we take the alternating magnetic field H1 = −H2 = H, then the total
magnetization is equal to zero, though the sublattice magnetizations remain non-trivial:
modd(H) =
1 +
1
π
∫ π
0
dν
θ2
4(ν,q2)
θ2
1(ν,q2)
cosh2 βH + sinh2 βH
tanh βH.
However, after substituting βH = h/(q−1 − q) and taking the limit α → 0, we see that our zero
temperature critical phenomenon disappears: modd(h) = 1 for h > 0 and modd(h) = −1 for
h < 0, similar to the free spins system.
4 Conclusion
As shown in [10], soliton solutions of integrable hierarchies with the complex values of spectral
variables are connected to the intrinsic Coulomb gases on two dimensional lattices with some
nontrivial dielectric or conductor boundaries. In this picture, the Coulomb interaction energy
between two charges and their effective images created by the boundary conditions plays the
role of the soliton phase shifts, the coordinates of charges coincide with the spectral parameters
of solitons, and the external electrostatic field is expanded in some series with the coefficients
playing the role of integrable hierarchy times. This transparent relation serves as a clue for
building new Coulomb lattice gas models exactly solvable at some fixed temperatures. The
latter temperatures are also related to the random matrix models [6, 9, 10], but we do not
discuss here applications of the described phase transition within these interpretations.
Our phase transitions are of a rather simple nature. In the lattice gas language, the transition
in the M = 1 periodic case describes the situation when the lattice is filled to its limit, i.e. the
number of particles equals to the number of sites, and no more particles can be added to the
system. In the M -periodic case, such transitions may happen separately, when each sublattice is
filled completely one by one, or simultaneously, for an appropriate choice of parameters. Their
qualitative features can be found from the mean field theory. There are several interesting
questions which would be interesting to analyze in the future, like influences of the hierarchy
times on the thermodynamical quantities, understanding of our systems beyond the “solitonic”
temperature values, investigation of the higher correlation functions, and so on.
Acknowledgements
The work of I.M. Loutsenko has been supported by European Community grant MIFI-CT-2005-
007323 and V.P. Spiridonov is partially supported by the Russian Foundation for Basic Research
(grant no. 05-01-01086). The authors are grateful to V.B. Priezzhev for useful remarks.
I vaguely remember a tall man with glasses vigorously explaining to me something during my
poster presentation at the IGTMP colloquium in Moscow in 1990. Probably that was Vadim –
I never asked him about that later on. We got acquainted at the first SIDE meeting near Montréal
in 1994 and had sufficiently long discussions during his few days visit to CRM after that confe-
rence. I remember telling him that by the work on separation of variables [8], which impressed me
much, he closed to me that field, and it is necessary to think about other directions of research.
We became closer during Vadim’s stay at the CRM in 1994–1995. It was very nice time from
many points of view. I visited his house and have known his family during one of the parties
he was gathering. It appeared that he likes Russian “bards” singing, which I was bond to as
well. Once he even sent to me a web-link to some new mp3-recordings coming from his native
Saint-Petersburg.
A Critical Phenomenon in Solitonic Ising Chains 11
In June 1999, he chaired my talk at the Hong Kong meeting on special functions, where
I reported results of the paper [13]. He was interested by that much, and we discussed possible
intersections with integrable systems. I saw Vadim last time at the conference in Edinburgh in
September 2003, where he was the main organizer. During the preparation of the corresponding
proceedings, I actively communicated with him, and it was clear that he is extremely busy by all
kinds of obligations. In May 2005, I suggested to him to form a team in order to try to get an
INTAS grant. His first reaction was positive, but after he has known the rules and procedures,
he rejected this idea by saying that there is too much bureaucracy and he has too many other
commitments for the next few months. I totally agreed with his critics and accepted his excuse.
Later on he listed to me a number of other possibilities to get research funding from the UK
sources which sounded quite reasonable. We exchanged by about ten e-mails with him over the
May–October 2005 period, and it was devastating to know that he has passed away. We have
another deeply regrettable loss in the FSU scientific community, which was possible, probably, to
prevent in other circumstances.
V.P. Spiridonov
I have met Vadim first when I have been pursuing my PhD studies at the Centre de Recherches
Mathématiques in Montréal. At that time, Vadim was a postdoctoral fellow there. I remember
him as an open-heart person, frank and honest, very enthusiastic and completely devoted to the
problems he did and had in mind. Always full of energy, he showed the keenest interest for many
questions of science and life. To meet Vadim was very interesting to me.
I.M. Loutsenko
References
[1] Ablowitz M.J., Segur H., Solitons and the inverse scattering transform, SIAM, Philadelphia, 1981.
[2] Andrews G.E., Askey R., Roy R., Special functions, Encyclopedia of Math. Appl., Vol. 71, Cambridge Univ.
Press, Cambridge, 1999.
[3] Baxter R.J., Exactly solved models in statistical mechanics, Academic Press, London, 1982.
[4] Date E., Jimbo M., Kashiwara M., Miwa T., Transformation groups for soliton equations, in Nonlinear
Integrable Systems, World Scientific, Singapore, 1983, 41–119.
[5] Evans M. R., Phase transitions in one-dimensional nonequilibrium systems, cond-mat/0007293.
[6] Gaudin M., Une famille à une paramètre d’ensembles unitaires, Nucl. Phys. 85 (1966), 545–575.
Gaudin M., Gaz coulombien discret à une dimension, J. Phys. (France) 34 (1973), 511–522.
[7] Hirota R., Exact solution of the Korteweg–de Vries equation for multiple collisions of solitons, Phys. Rev.
Lett. 27 (1971), 1192–1194.
[8] Kuznetsov V.B., Quadrics on real Riemannian spaces of constant curvature: separation of variables and
connection with Gaudin magnet, J. Math. Phys. 33 (1992), 3240–3254.
[9] Loutsenko I.M., Spiridonov V.P., Self-similar potentials and Ising models, Pis’ma v ZhETF (JETP Letters)
66 (1997), 747–753.
Loutsenko I.M., Spiridonov V.P., Spectral self-similarity, one-dimensional Ising chains and random matrices,
Nucl. Phys. B 538 (1999), 731–758.
[10] Loutsenko I.M., Spiridonov V.P., Soliton solutions of integrable hierarchies and Coulomb plasmas, J. Stat.
Phys. 99 (2000), 751–767, cond-mat/9909308.
[11] Matveev V.B., Salle M.A., Darboux transformations and solitons, Springer Series in Nonlinear Dynamics,
Springer-Verlag, 1991.
[12] Spiridonov V.P., Universal superpositions of coherent states and self-similar potentials, Phys. Rev. A 52
(1995), 1909–1935, quant-ph/9601030.
[13] Spiridonov V.P., Zhedanov A.S., Spectral transformation chains and some new biorthogonal rational func-
tions, Comm. Math. Phys. 210 (2000), 49–83.
http://arxiv.org/abs/cond-mat/0007293
http://arxiv.org/abs/cond-mat/9909308
http://arxiv.org/abs/quant-ph/9601030
1 The Korteweg-de Vries solitonic spin chain
2 The mean field approximation
3 The BKP solitonic spin chain
4 Conclusion
References
|