Journal of Nonlinear Mathematical Physics

Volume 26, Issue 3, May 2019, Pages 468 - 482

Quasiperiodic Solutions of the Heisenberg Ferromagnet Hierarchy

Authors
Peng Zhao
1College of Arts and Sciences, Shanghai Maritime University, Shanghai, 201306, People’s Republic of China
2College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao, 266590, People’s Republic of China
3School of Mathematical Sciences, Qufu Normal University, Qufu, 273165, People’s Republic of China,pengzhao@shmtu.edu.cn
Engui Fan*
School of Mathematical Sciences, Institute of Mathematics and Key Laboratory of Mathematics for Nonlinear Science, Fudan University, Shanghai, 200433, People’s Republic of China,faneg@fudan.edu.cn
Temuerchaolu
College of Arts and Sciences, Shanghai Maritime University, Shanghai, 201306, People’s Republic of China,tmchaolu@shmtu.edu.cn
*Corresponding author
Corresponding Author
Engui Fan
Received 17 October 2018, Accepted 25 March 2019, Available Online 6 January 2021.
DOI
10.1080/14029251.2019.1613053How to use a DOI?
Keywords
Heisenberg ferromagnet hierarchy; spectral curve; Riemann theta function; quasiperiodic solution
Abstract

We present quasi-periodic solutions in terms of Riemann theta functions of the Heisenberg ferromagnet hierarchy by using algebrogeometric method. Our main tools include algebraic curve and Riemann surface, polynomial recursive formulation and a special meromorphic function.

Copyright
© 2019 The Authors. Published by Atlantis and Taylor & Francis
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

1. Introduction

Integrable discretizations of soliton equations have attracted much attention during the past few years. The celebrated Landau-Lifshits equation [2]

st=[s,sxx+Js],s3,J=diag(J1,J2,J3),J1+J2+J3=0,(1.1)
has two well-known integrable discretizations. One of them is the Sklyanin lattice which is determined by the Poisson brackets algebra (see (13) in [7] or [1])
{sα(n),s0(n)}=Jβγsβ(n)sγ(n),{sα(n),sβ(n)}=s0(n)sγ(n),(1.2)
and the discrete variant of the Landau-Lifshits Hamiltonian (see (23) in [7])
H(0)=nln(s0(n+1)s0(n)+α=13(K1K0Jα)sα(n+1)sα(n)).(1.3)

Another discretization of the Landau-Lifshits equation (1.1) is the so-called Shabat-Yamilov lattice [5,6]

ut=2hu+v+hv,vt=2huvhu,(1.4)
where u = u(n, t), v = v(n, t), ut=du(n,t)dt, vt=dv(n,t)dt, u+ = E+u = u(n + 1, t), v = Eb = b(n − 1, t) (E± are the shift operators and (n, t) ∈ 𝕑 × ℝ), h = h(u,v) is a symmetric polynomial of u, v degree no higher than two with respect to each of the variables and huuu = 0. By introducing the complexified stereographic projection
S=S(u,v)=1uv(1uv,i+iuv,u+v),i=1,(1.5)
one may relate the Shabat-Yamilov lattice (1.4) to the Sklyanin lattice by a special transformation of variables [1]. The linear combination of the flows of (1.4) and its shift
ut=2huv+hv,vt=2huv+hu,(1.6)
give rise to [5]
ut=ρ1(2hu+v+hv)+ρ2(2huv+hv),vt=ρ1(2huvhu)+ρ2(2huv+hu),(1.7)
where ρ1, ρ2 are real constants. By specifying ρ1 = 1, ρ2 = 0, h=12(uv)2, the lattice (1.7) reduces to the so-called Heisenberg ferromagnet (HF) lattice
{ut=(uv)(uu+)((u+v)1,vt=(uv)(vv)(uv)1.(1.8)

In the spin variables S (see (1.5)), the system (1.7) reads [1]

St=ρ1S,KS,([S,S+]1+[S,S+]+[S,S]1+[S,S])2ρ1[S,KS]+ρ2S,KS(S+S+1+[S,S+]+S+S1+[S,S]),(1.9)
where K = diag(K1, K2, K3) and |S| = 1. It is well-known that the lattice (1.9) is integrable with a zero-curvature representation and can be reduced to the well-known Heisenberg lattice (see [1] and references therein).

In this paper, we study quasi-periodic solutions of the whole Heisenberg ferromagnet hierarchy by using algebrogeometric method. In section 2, we construct the stationary and time-dependent HF hierarchy from its zero-curvature representation. Then spectral curves and an auxiliary function φ are introduced in section 3. In section 4, based on analytic and asymptotic properties of φ, we derive theta function representations for the entire HF hierarchy.

2. Heisenberg ferromagnet hierarchy

In this section, we construct the Heisenberg ferromagnet hierarchy by developing zero-curvature formulism of the HF lattice. For later use we denote by E± the shift operators acting on ψ={ψ(n,t)}n=+𝕑 according to (E± ψ)(n, t) = ψ(n±1, t), or E± ψ = ψ± for convenience. These notations are followed by [4].

It is well-known that the HF lattice has symplectic operator, Hamiltonian structure, recursion operator, nontrivial generalised symmetry and Lax representation [5]. To construct the HF hierarchy we start from the following two 2 × 2 Lax matrices

U=U(λ,u,v)=(λ2u(uv)12(uv)12uv(uv)1λ+2v(uv)1),(2.1)
V=V(λ,u,v)=λ1(uv)1(u+v22uv(u+v)),(2.2)
where λ ∈ ℂ is a spectral parameter and u, v are functions of the lattice variable n and the time variable t.

In the following we temporarily view u, v as functions of only n and define the sequence {a, b, c}∈ℕ0

2δ(a+a+ubvb+)=b1+b1,0,(2.3)
2uvδ(a+a++c/u+c+/v)=c1+c1,0,(2.4)
2vδ(aa+)2uvδb2δc+=a1+a1+,0,(2.5)
where a, b, c are polynomials of u, v and their shifts, and δ = (uv)−1. In matrix form, the relations (2.3)(2.5) can be written as
(I+E+uvE+0I+E+0u1+v1E+IE+uv1E+)(abc)=(0(2δ)1(E+I)000(2uvδ)1(E+I)E+I00)(a1b1c1),(2.6)
where I is the identity operator from ℂ𝕑 to ℂ𝕑. Apparently, once the initial value (a0, b0, c0) is given, {a, b, c}∈ℕ can be recursively determined by the relations (2.3)(2.5) or (2.6). However, the calculations involved are rather big and therefore it is uneasy to obtain explicit form of {a, b, c}∈ℕ. In spite of difficulties, we can still get the following result.

Theorem 2.1.

Solutions of the system (2.3)(2.5) are explicitly given by the following recursion relations

a=u+vuv(E+I)1S1,1,0,(2.7)
b=S2,1+2uv(E+I)1S1,1,0,(2.8)
c=S3,12uvuv(E+I)1S1,1,0, (2.9)
where S1,−1, S2,−1, S3,−1 are defined by
S1,1=uu2(v)2(v(b1+b1)+a1a1+2δ)v(u+)2v2(u(b1+b1)+a1a1+2δ)+uu2(v)2u(b1+b1)+a1a1+2δ+v(u+)2v2v(b1+b1)+a1a1+2δ+b1+b12δ,
S2,1=1u2(v)2((v(b1+b1)+a1a1+2δ)u(b1+b1)+a1a1+2δ),
S3,1=uvu2(v)2(uv(v(b1+b1)+a1a1+2δ)vuv(b1+b1)+a1a1+2δ).

Proof.

First, from (2.3)(2.5) it follows

2uδ(aa+)=2uvδb+2δc+a1a1+.(2.10)

Eliminating b+ from (2.3) and (2.10), we have

2uδ(2aub+cu)=u(b1+b1)+a1a1+.(2.11)

Then insertion of (2.3) into (2.5) yields

4vδa+2v2δb++2δc+=v(b1+b1)+a1a1+,(2.12)
or equivalently,
2avb+cv=(v(b1+b1)+a1a1+2δ).(2.13)

Combining (2.11) with (2.13), we obtain

b+1uvc=(v(b1+b1)+a1a1+2vδ)u(b1+b1)+a1a1+2uδ(2.14)
b2au+v=1u2(v)2(v(b1+b1)+a1a1+2vδ)1u2(v)2×u(b1+b1)+a1a1+2δ.(2.15)

Thus, inserting (2.15) into (2.3), we conclude (2.7) and (2.8) hold. Finally, the formula (2.9) can be derived from the following two equalities

2uauvb+uvc=u(v(b1+b1)+a1a1+2vδ),(2.16)
2vauvb+vuc=v(v(b1+b1)+a1a1+2uδ).(2.17)
by eliminating b.

Next we compute the stationary HF hierarchy. To this end we start from the stationary zero-curvature equation

0=UW(EW)U,W=(ABCA),(2.18)
where A, B, C, D are polynomials of λ:
A==1N+1aN+1λ,B==1N+1bN+1λ,C==1N+1cN+1λ,N.(2.19)

Using (2.1), one may rewrite Eq. (2.18) as

(λ2uδ)A2δC=(λ2uδ)A++2uvδB+,(2.20)
(λ2uδ)B+2δA=2δA++(λ+2vδ)B+,(2.21)
2uvδA+(λ+2vδ)C=(λ2uδ)C+2uvδA+,(2.22)
2uvδB(λ+2vδ)A=2δC+(λ+2uv)δA+.(2.23)

Inserting (2.19) into (2.20)(2.23) and taking into account (2.3)(2.5), one derives the stationary HF hierarchy

aNaN+=0,bNbN+=0,cNcN+=0.(2.24)

In the case N = 0, the system (2.24) is reduced to the stationary HF lattice

(uv)(uu+)(u+v)1=0,(uv)(vv)(uv)1=0.(2.25)

The time-dependent HF hierarchy can be derived as follows. First, a, b, c are considered as functions of both n and t. To distinguish different higher order flows that are related to different time variables, we replace t by tr. Then inserting (2.1), (2.3)(2.5) and (2.19) into the time-dependent zero-curvature equation

0=Utr+UW(E+W)U,(2.26)
we obtain the r-th equation in the HF hierarchy
(2uδ)tr=ar+ar+,(2.27)
(2δ)tr=br+br+,(2.28)
(2uvδ)tr=crcr+,(2.29)
(2vδ)tr=ar+ar+.(2.30)

The system (2.27)(2.30) is overdetermined and can be simplified into the r-th HF lattice.

Theorem 2.2.

The r-th HF lattice has the form of

utr=uv2(u(brbr+)ar+ar+),(2.31)
vtr=uv2(v(brbr+)ar+ar+),r0.(2.32)

Proof.

First, using (2.27), we obtain

2utrδ+2uδtr=ar+ar+.(2.33)

Then the relations (2.28) and (2.33) give rise to

2utrδ=u(brbr+)ar+ar+,(2.34)
and hence (2.31) holds. Using (2.28) and (2.30), we arrive at (2.32). To complete the proof, it remains to show that (2.28), (2.29) are compatible with (2.31), (2.32). Actually, using (2.31) and (2.32), we obtain
(2δ)tr=(2uv)tr=2δ2(utrvtr)=δ(u(brbr+)ar+ar+)+δ(v(brbr+)ar+ar+)=br+br+.(2.35)

Then by (2.31), (2.32) and (2.35), it follows

2utrvδ=uv(brbr+)v(arar+),2uvtrδ=uv(brbr+)u(arar+),2uvδtr=uv(brbr+),
and consequently,
(2uvδ)tr=2utrvδ+2uvtrδ+2uvδtr=uv(brbr+)(u+v)(arar+)=crcr+.(2.36)

Here we use the recursion relations

2vδ(aa+)2uvδb2δc+=a1+a1+,0,2uδ(aa+)+2uvδb+2δc+=a1a1+,0
in the last equality of (2.36).

3. Spectral curves and a basic meromorphic function

In this section, we first introduce the spectral curves associated with the HF hierarchy. Then we introduce a basic meromorphic function φ and study its analytic and asymptotic properties.

Let (ψ1, ψ2)T and (ϕ1, ϕ2)T be two fundamental solutions of auxiliary linear problem

E+Φ(λ,n,tr)=U(λ,n,tr)Φ(λ,n,tr),Φtr(λ,n,tr)=W(λ,n,tr)Φ(λ,n,tr),(3.1)
where Φ(λ, n, tr) = (Φ1(λ, n, tr),Φ2(λ, n, tr))T . Moreover, we define
f(λ,n,tr)=(λ22λ)nψ1(λ,n,tr)φ1(λ,n,tr),(3.2)
h(λ,n,tr)=(λ22λ)nψ2(λ,n,tr)φ2(λ,n,tr),(3.3)
g(λ,n,tr)=21(λ22λ)n(ψ1(λ,n,tr)φ2(λ,n,tr)+φ1(λ,n,tr)ψ2(λ,n,tr)).(3.4)

Theorem 3.1.

The functions f, g, h defined in (3.2)(3.4) satisfy

(gfhg)+=U(gfhg)U1,(3.5)
(gfhg)tr=[(ABCA),(gfhg)],(3.6)
where [·,·] is the Lie bracket of two matrices P, Q
[P,Q]=PQQP
and g2 + fh is independent of x and tr.

Proof.

The relation (3.5) can be derived from (3.1), (3.2), (3.3) and (3.4). Indeed, we have

g+=21(λ22λ)(n+1)ψ1+φ2++φ1+ψ2+=21(λ22λ)(n+1)(((λ2uδ)ψ12δψ2)(2uvδφ1+(λ+2vδ)φ2)+((λ2uδ)φ12δφ2)(2uvδψ1+(λ+2vδ)ψ2))=(λ22λ)1(2uvδ(λ2uδ)f(λ2uδ)(λ+2vδ)g+4uvδ2g+2δ(λ+2vδ)h,
f+=(λ22λ)(n+1)ψ1+φ1+=(λ22λ)(n+1)((λ2uδ)ψ12δψ2)((λ2uδ)φ12δφ2)=(λ22λ)1((λ2uδ)(λ2uδ)f+2δ(λ2uδ)2g(2δ)2h),
h+=(λ22λ)(n+1)ψ2+φ2+,=(λ22λ)(n+1)(2uvδψ1+(λ+2vδ)ψ2)(2uvδφ1+(λ+2vδ)φ2)=(λ22λ)1((2uvδ)2f+2uvδ(λ+2vδ)2g+(λ+2vδ)2h).

Similarly, according to the definition of f, g, h, ϕ1, ϕ2, ψ1, ψ1, it follows that

gtr=21(λ22λ)n(ψ1,trφ2+ψ1φ2,tr+φ1,trψ2+φ1ψ2,tN)=21(λ22λ)n((Aψ1+Bψ2)φ2+ψ1(Cφ1Aφ2)+(Aφ1+Bφ2)ψ2+φ1(Cψ1Aψ2))+BhCf,
ftr=(λ22λ)n(ψ1φ1)tr=(λ22λ)n(ψ1,trφ1+ψ1φ1,tr)=(λ22λ)n((Aψ1+Bψ2)φ1+ψ1(Aφ1+Bφ2))=2Af2Bg
htr=(λ22λ)n(ψ2φ2)tr=(λ22λ)n(ψ2,trφ2+ψ2φ2,tr)=(λ22λ)n((Cψ1Aψ2)φ2+ψ2(Cφ1Aφ2))=2Cg2Ah,
which indicates (3.6) holds. Finally, using
(det(gfhg))tr=det(gfhg)tr((gfhg)tr(gfhg)1),(3.7)
det(gfhg)=g2fh,(3.8)
we conclude that g2 + fh is independent of n and tr.

In what follows it is convenient to introduce

gM==1M+1aM+1λλ,fM==1M+1bM+1λ,hM==1M+1cM+1λ,M0.(3.9)

Then one may obtain the following result.

Theorem 3.2.

Assume u, v are solutions of the r-th HF lattice and the M-th stationary HF lattice (2.24). Then for fixed M ∈ ℕ0,

g=gM,f=fM,h=hM,(3.10)
solve the system (3.6).

Proof.

First we show that (3.10) satisfy the first equation in (3.6). To this end we make the ansatz

g==1M+1aM+1λλ,f==1M+1bM+1λ,h==1M+1cM+1λ,(3.11)
where aj, bj, cj, j = 1,...,M +1 are polynomials with respect to u, v and their shift. Inserting (3.11) into (3.6), one derives
aj=aj,bj=bj,cj=cj,j=1,,M+1.(3.12)
and hence (3.10) and (3.9) hold. To complete the proof, one has to prove that g, f, h also satisfy (3.6). Since u, v are both solutions of the M-th time-independent lattice, there exists a common eigenfunction χ = (χ1, χ2) for the following two linear problems
χ+=Uχ,VMχ=λ(M+1)yχ,(3.13)
where the 2 × 2 matrix
VM=(gMfMhMgM).(3.14)

The relation (3.13) implies

λM+1(gM2+fMhM)1/2=y.

Next we introduce a new function

ϕ=χ2χ1=λ(M+1)ygMfM=hMλ(M+1)y+gM.(3.15)

A direct computation shows that φ′ satisfies

(λ2uδ)ϕ+2δϕϕ=2uvδ+(λ+2vδ)ϕ.(3.16)

Differentiating (3.16) with respect to tN then yields

((λ2uδ2δϕ)E+2δϕ+(λ+2vδ))ϕtr=(2uvδ)tr+(2vδ)trϕ+(λ+2vδ)ϕtr.(3.17)

On the other hand, we have

((λ2uδ2δϕ)E+2δϕ+(λ+2vδ))(C2AϕB(ϕ)2)=(λ2uδ2δϕ)C+2δϕ+C(λ+2vδ)C+(2uvδ+(λ+2vδ)ϕϕ+E++2uvδ(λ2uδ)ϕ+ϕ)(2AϕB(ϕ+)2)=(2uvδ)tr+(2vδ)trϕ+(λ+2vδ)ϕtr.(3.18)

Combining (3.17) with (3.18), one obtains

(λ2uδ2δϕ)E+2δϕ+(λ+2vδ))(ϕtrC+2AϕB(ϕ)2)=0,(3.19)
and consequently,
ϕtr=C2AϕB(ϕ)2.(3.20)

On the other hand, from the definition of ψi, ϕi, one infers that the function φ (see (3.32)) satisfies

(λ2uδ)ϕ+2δϕ+ϕ2uvδ(λ+2vδ)ϕ=0.(3.21)

Then by (3.1), we have

(ψ2ψ1)tr=ψ2,tNψ1ψ2ψ1,tNψ12=Cψ1Aψ2ψ1ψ2ψ1Aψ1+Bψ2ψ1=C2Aψ2ψ1B(ψ2ψ1)2.(3.22)

Similarly one may derive

(φ2φ1)tr=C2Aφ2φ1B(φ2φ1)2.(3.23)

Thus one arrives at

ϕtr=C2AϕBϕ2.(3.24)

A comparison of (3.17), (3.20), (3.21) and (3.24) yields

ϕ=F(λ)ϕ,(3.25)
where F(λ) is an arbitrary function of λ. Without loss of generality, one can take F(λ) ≡ 1. Therefore, one infers from (3.25) that there exists common eigenfunctions of the following linear system
χ+=Uχ,(3.26)
χtr=Vχ,(3.27)
VMχ=λ(M+1)yχ.(3.28)

From the compatibility condition of (3.27) and (3.28), we finally obtain

VM,tr=[V,VM],(3.29)
which completes the proof.

Next let us define

Δ±(λ,n,tr)=±21(λ22λ)n(ψ1(λ,n,tr)φ2(λ,n,tr)φ1(λ,n,tr)ψ2(λ,n,tr)).(3.30)

It is not difficult to check that Δ± satisfy

Δ±2(λ,n,tr)=g2(λ,n,tr)+f(λ,n,tr)h(λ,n,tr),Δ+(λ,n,tr)+g(λ,n,tr)=Δ(λ,n,tr)+g(λ,n,tr)=(λ22λ)nψ1(λ,n,tr)φ2(λ,n,tr),φ2(λ,n,tr)/φ1(λ,n,tr)=(Δ(λ,n,tr)g(λ,n,tr))/f(λ,n,tr),ψ2/ψ1(λ,n,tr)=h(λ,n,tr)/(Δ+(λ,n,tr)+g(λ,n,tr))=(Δ+(λ,n,tr)g(λ,n,tr))/f(λ,n,tr).

For fixed M ∈ ℕ, the spectral curve 𝒦M−1 of genus M − 1 can be introduced as follows

𝒦M1:y2=λ2M+2(f2+gh)=j=02M1(λEj),Ej.

Throughout this paper, one assumes that

EjE,j,j=0,1,,2M1.(3.31)

The next step is crucial; we lift the functions ϕ21 and ψ21 on the Riemann surface 𝒦M−1 by

ϕ(P,x,tr)={ψ2(λ,n,tr)/ψ1(λ,n,tr),P𝒦M1+φ2(λ,n,tr)/φ1(λ,n,tr),P𝒦M1,(3.32)
where 𝒦M1± are two sheets of 𝒦M−1. The spectral curve 𝒦M−1 is compactified by the point at infinity P. A point on 𝒦M−1 is denoted by P = (λ,y(P)), where λ ∈ ℂ and y(P) is a holomorphic function defined on the two sheets 𝒦M1± of 𝒦M−1:
y(P)={j=02M1(λEj),P𝒦M1+j=02M1(λEj),P𝒦M1.(3.33)

In particular, one introduces P0±=(0,±(j=02M1Ej)1/2).

Based on above preparations, we turn to study the analytic and asymptotic properties of φ.

Lemma 3.3.

The function φ defined in (3.32) is meromorphic. Its divisor is given by a

(ϕ(,n,tr))=𝒟ν^_(n,tr)ν^M(n,tr)𝒟μ^_(n,tr)μ^M(n,tr).(3.34)

Here the notations ν^_(n,tr) and ν^_(n,tr) represent two vectors

ν^_(n,tr)=(ν^1(n,tr),,ν^M1(n,tr)),μ^_(n,tr)=(μ^1(n,tr),,μ^M1(n,tr)).

Lemma 3.4.

Near the points P, the function φ has the following expansions

ϕ(,n,tr)=c++O(ζ),PP±,ζ=λ1(3.35)
and
ϕ(,n,tr)={vE1(vv)(uv)2(uv)ζ+O(ζ2),PP0+,ζ=λ,u+(u+u)(uv)2(u+v)ζ+O(ζ2).PP0,ζ=λ,(3.36)
where c∞± are lattice constants.

Proof.

First one observes that the expansion (3.35) can be easily derived from the definition of φ. Thus it suffices to compute the expansions of φ near P. Using the relation

Δ±2=g2+fh=λ2(M+1)((i=0Ngiλi)2(i=0Nfiλi)(i=0Nhiλi))=λ2(M+1)i=0N(k=0i(gkgikfkhik))λi,
one infers
Δ±=ζ0±ζ(M+1)g02+f0h0(1+122g0g1+f0h1+f1h0g02+f0h0ζ+O(ζ2)),(3.37)
as PP. Then from the relation g02+f0h0=1 and (3.32), it follows
ϕ(,n,tr)={1u+vuv2uv+O(ζ),PP0+,λ=ζ,1u+vuv2uv+O(ζ),PP0,λ=ζ,={v+O(ζ),PP0+,λ=ζ,u+O(ζ),PP0,λ=ζ.(3.38)

To get the coefficient of ζ in (3.36), we use the equation (see (3.16))

(λ2uδ)ϕ+2δϕ+ϕ=2uvδ+(λ+2vδ)ϕ.(3.39)

Inserting the ansatz

ϕ=ϕ0±+ϕ1±ζ+ϕ2±ζ2+O(ζ3),asPP0±(3.40)
into (3.39) and comparing the coefficients of ζ lead to
ϕ0±+2uδϕ1±+2δϕ0±+ϕ1±2δϕ0±ϕ1±+=ϕ0±+2vδϕ1±.(3.41)

Explicitly, Eq. (3.41) can be equivalently written as

v2uδϕ1±++2vδϕ1++2vδϕ1++=v+2vδϕ1+,(3.42)
u+2uδϕ1++2u+δϕ1+2uδϕ1+=u+2vδϕ1,(3.43)
which implies
ϕ1+=E1(vv)(uv)2(vu),ϕ1=(u+u)(uv)2(u+v).(3.44)

Using similar method we can prove (3.35).

4. Quasi-periodic solutions

In this section, we obtain theta function representations for the meromorphic function φ and quasiperiodic solutions u, v.

First, we choose a convenient base point Q0KM−1\{P∞±, P}. The Abel maps A_Q0(), α_Q0() are defined by

A_Q0:XJ(X)=M1/LM1,
PA_Q0(P)=(AQ0,1(P),,QQ0,M1(P))=(Q0Pω1,,Q0PωM1)(modLM1),
and
α_Q0:Div(X)J(X),𝒟α_Q0(𝒟)=P𝒦M1𝒟(P)A_Q0(P),
where LM1={z_M1|z_=N_+ΓM_,N_,M_𝕑M1}, and Γ, Ξ_Q0 are the Riemann matrix and the vector of Riemann constants, respectively. Moreover, we choose a homology basis {aj,bj}j=1M1 on X in such a way that the intersection matrix of the cycles satisfies
ajbk=δj,k,ajak=0,bjbk=0,j,k=1,,n.(4.1)

For brevity, we introduce

z_(P,Q_)=Ξ_Q0A_Q0(P)+α_Q0(𝒟Q_),P𝒦M1,Q_=(Q1,,QM1)σM1𝒦M1,(4.2)
where z_(,Q_) is independent of the choice of Q0. The Riemann theta function θ(z_) associated with X and the homology homology basis {aj,bj}j=1M1 is defined by
θ(z_)=n_𝕑exp(2πin_,z_+πin_,n_Γ),z_M1,
where B_,C_=j=1M1B¯jCj is the scalar product in ℂM−1.

Next, we introduce the differential of the third kind with simple zero and pole respectively at ν^M(n,tr) and μ^M(n,tr) by

ων^M(n,tr)μ^M(n,tr)(3)(P)=(y+y(ν^M(n,tr))λνM(n,tr)y+y(μ^M(n,tr))λμM(n,tr))dλ2y+λ0yi=1M2(λλj)dλ,
where λj, j = 0,...,M − 2, are constants that are uniquely determined by the requirement of vanishing a-period, i.e.
ajων^M(n,tr)μ^M(n,tr)(3)(P)=0,j=1,,M1.(4.3)

A simple computation shows

Q0Pων^M(n,tr)μ^M(n,tr)(3)(P)={lnζ+d0(n,tr)+O(ζ),Pν^M(n,tr),lnζ+d1(n,tr)+O(ζ),Pμ^M(n,tr),dP0±+O(ζ)PP0±,dP±+d2(n,tr)ζ+O(ζ2),PP±.

Here d0(n,tr), d1(n,tr), d2(n,tr) are functions of variables n, tr and dP∞±, dP are integration constants.

Theta function representations for quasiperiodic solutions of the rth HF lattice can be obtained as follows.

Theorem 4.1.

The function φ defined in (3.32) is meromorphic on 𝒦M−1 and has the following theta function representation

ϕ(P,n,tr)=cP+θ(λ_(P+,μ^_(n,tr)))θ(λ_(P+,ν^_(n,tr)))θ(λ_(P,ν^_(n,tr)))θ(λ_(P,μ^_(n,tr)))exp(Q0Pων^M(n,tr)μ^M(n,tr)(3)dP+).(4.4)

Proof.

To derive the formula (4.4), one has to use Riemann vanishing theorem and Riemann-Roch theorem in standard literature [3]. From Lemma 2.3, one infers that the function φ(P,n,tN) take the form

ϕ(P,n,tr)=C(n,tr)θ(λ_(P,ν^_(n,tr)))θ(λ_(P,μ^_(n,tr)))exp(Q0Pων^M(n,tr)μ^M(n,tr)(3)),(4.5)
where C(n,tr) is an undetermined function. Then taking the limit PP∞+ on both sides of (4.5), one derives
cP+=C(n,tN)θ(λ_(P+,ν^_(n,tr)))θ(λ_(P+,μ^_(n,tr)))exp(dP+).(4.6)

Inserting(4.6) into (4.5) and eliminating C(n,tr), one finally obtains (4.4).

Next, we obtain theta function representations for solutions u,v with the help of (4.4).

Theorem 4.2.

Quasiperiodic solutions of the r-th HF lattice have the following Riemann theta function representation

dP0+-dP+,(4.7)
v(n,tr)=-cP+θ(λ_(P+,μ^_+(n,tr)))θ(λ_(P+,ν^_+(n,tr)))θ(λ_(P0+,ν^_+(n,tr)))θ(λ_(P0+,μ^_+(n,tr)))exp(dP0+-dP+).(4.8)

Proof.

Using the coordinate λ = ζ near P, we have

ϕ(n,x,tr)=cP+θ(λ_(P+,μ^_(n,tr)))θ(λ_(P+,ν^_(n,tr)))θ(λ_(P,ν^_(x,tr)))θ(λ_(P,μ^_(x,tr)))×exp(Q0Pων^M(x,tr)μ^M(x,tr)(3)dP+)=ζ0cP+θ(λ_(P+,μ^_(n,tr)))θ(λ_(P+,ν^_(n,tr)))θ(λ_(P0±,ν^_(x,tr)))θ(λ_(P0±,μ^_(x,tr)))×exp(dP0±dP+)(1+O(ζ)).(4.9)

Comparing the leading coefficients of (3.36) and (4.9) naturally gives rise to (4.7) and (4.8).

Acknowledgements

Our deepest gratitude goes to the anonymous reviewers for their careful work and thoughtful suggestions that have helped improve this paper. The work of Fan in this paper was supported by grants from the National Science Foundation of China (Project Nos. 11671095, 51879045 and 11271079). Zhao was supported by grants from the National Science Foundation of China (Project Nos.11526137, 11547199 and 11601321). Temuerchaolu was supported by grants from the National Science Foundation of China (Project No.11571008).

Footnotes

a

For the meaning of these notations one may refer to [4].

Journal
Journal of Nonlinear Mathematical Physics
Volume-Issue
26 - 3
Pages
468 - 482
Publication Date
2021/01/06
ISSN (Online)
1776-0852
ISSN (Print)
1402-9251
DOI
10.1080/14029251.2019.1613053How to use a DOI?
Copyright
© 2019 The Authors. Published by Atlantis and Taylor & Francis
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

Cite this article

TY  - JOUR
AU  - Peng Zhao
AU  - Engui Fan
AU  - Temuerchaolu
PY  - 2021
DA  - 2021/01/06
TI  - Quasiperiodic Solutions of the Heisenberg Ferromagnet Hierarchy
JO  - Journal of Nonlinear Mathematical Physics
SP  - 468
EP  - 482
VL  - 26
IS  - 3
SN  - 1776-0852
UR  - https://doi.org/10.1080/14029251.2019.1613053
DO  - 10.1080/14029251.2019.1613053
ID  - Zhao2021
ER  -