On solutions to nonlinear integral equation of the Hammerstein type and its applications to Gibbs measures for continuous spin systems Текст научной статьи по специальности «Математика»

NANOSYSTEMS:

PHYSICS, CHEMISTRY, MATHEMATICS Original article

Mavlonov I.M., et al. Nanosystems: Phys. Chem. Math., 2024,15 (1), 23-30.

http://nanojournal.ifmo.ru DOI 10.17586/2220-8054-2024-15-1-23-30

On solutions to nonlinear integral equation of the Hammerstein type and its applications to Gibbs measures for continuous spin systems

IsmoilM. Mavlonov3 ", AloberdiM. Sattarov4,b, Sevinchbonu A. Karimova3, FarhodH. Haydarov1,2,3,c

1 Institute of Mathematics, Tashkent, Uzbekistan 2New Uzbekistan University, Tashkent, Uzbekistan 3National University of Uzbekistan, Tashkent, Uzbekistan 4University of Business and Science, Namangan, Uzbekistan

[email protected], [email protected],[email protected],

Corresponding author: I.M. Mavlonov, [email protected]

Abstract The paper deals with the problem of constructing kernels of Hammerstein-type equations whose positive solutions are not unique. This problem arises from the theory of Gibbs measures, and each positive solution of the equation corresponds to one translation-invariant Gibbs measure. Also, the problem of finding kernels for which the number of positive solutions to the equation is greater than one is equivalent to the problem of finding models which has phase transition. In these articles, the number of solutions corresponding to the constructed kernels does not exceed 3, and in turn, it only gives us a chance to check the existence of phase transitions. The constructed kernels in this paper are more general than the kernels in the above-mentioned papers and except for checking phase transitions, it allows us to classify the set of Gibbs measures.

Keywords Generalized SOS model, spin values, Cayley tree, gradient Gibbs measure, periodic boundary law

Acknowledgements We thank the referee for careful reading of the manuscript; in particular, for some comments, useful discussions, and suggestions which have improved the paper.

For citation Mavlonov I.M., Sattarov A.M., Karimova S.A., Haydarov F.H. On solutions to nonlinear integral equation of the Hammerstein type and its applications to Gibbs measures for continuous spin systems. Nanosystems: Phys. Chem. Math., 2024, 15 (1), 23-30.

1. Introduction

The Hammerstein equation covers a large variety of areas and is of much interest to a wide audience due to the fact that it has applications in numerous areas. Several problems arise in differential equations (ordinary and partial), for instance, elliptic boundary value problems whose linear parts possess Green's function can be transformed into the Hammerstein integral equations. Equations of the Hammerstein type play a crucial role in the theory of optimal control systems and in automation and network theory (see e.g., [1]).

The Hammerstein equation is a mathematical model that describes the response of a system to an input signal. It is commonly used in control systems engineering and signal processing.

Let G(u(t)) represent the static nonlinearity component of the system, which is a function of the input signal, and H(u(t)) represent the dynamic linear component of the system, which is a linear function of the input signal. Then the equation is typically represented as: y(t) = G(u(t)) + H(u(t)), where y(t) and u(t) represent the output and input signal to the system at time t.

The Hammerstein equation allows one modeling of systems that exhibit both linear and nonlinear behavior, making it useful in various applications.

In statistical mechanics, the Gibbs measure is a probability distribution that describes the equilibrium state of a system. It is commonly used to describe systems with countable spin values, where the spins can take on a finite or countably infinite number of discrete values. However, there are also models with uncountable spin values, where the spins can take uncountably many values, such as in continuous spin systems. In these cases, the Gibbs measure needs to be adapted to handle the uncountable nature of the spin values.

The definition of the Gibbs measure for models with uncountable spin values often involves considering the Hamil-tonian or energy function that governs the system. The Gibbs measure assigns a probability to each spin configuration based on its energy, such that spin configurations with lower energy have higher probabilities.

Specifically, for a given spin configuration, the probability assigned by the Gibbs measure is proportional to the exponential of the negative energy of that configuration. The proportionality constant is chosen such that the total probability

© Mavlonov I.M., Sattarov A.M., Karimova S.A., Haydarov F.H., 2024

assigned by the Gibbs measure sums to 1 over all possible spin configurations. Note that all the above-mentioned and additional information can be found in References [2-4].

In recent years, increasing attention was given to models with a uncountable many spin values on particle systems (especially, the Cayley tree). By this time, models with four interactions with the set of spin values [0,1] were considered, and the following results were achieved: In [5] splitting Gibbs measures on the Cayley tree of order k are described by solutions to a special nonlinear integral equation of Hammerstain's type. The existence of positive solutions to the nonlinear integral equation of Hammerstain's type and a sufficient condition of the uniqueness of the equation given in [6]. Also, in [7-10], results on non-uniqueness of positive solutions to the equation with non-degenerate kernels are given and for the case of degenerate kernel can be seen in [11,12]. Finally, in [13,14] translation-invariant Gibbs measures for models with four competing interactions are described by solutions to the Lyapunov integral equation. In this paper, we construct more general kernels which gives us a chance to check the existence of phase transitions and to classify the set of Gibbs measures.

2. Special equation of Hammerstein type

Let Q be a set. A nonlinear integral equation of Hammerstein type on Q is one of the form

u(x)+ f K(x,y)f (l(y))dy = h(x) (2.1)

Jo

where dy stands for a a- finite measure on the measure space Q; the mappings l and K are measurable on Q and Q x Q respectively, f is a real-valued function defined on R and is in general nonlinear and h is given function on Q. If we define the operator A by

Av(x)= K (x,y)v(y)dy, (2.2)

Jo

and Nf to be the Nemystkii operator associated with f :

Nf l(x) = f (l(x)) (2.3)

then integral equation (2.1) can be transformed to functional equation form as follows:

l + ANf l = 0, (2.4)

where without the loss of generality, we have taken h = 0 (see detail in [15]). We consider a special equation of Hammerstein type, i.e.

i K(x,y)lk(y)dy = l(x), k e N. (2.5)

J 0

Here and below, K : [0,1]2 ^ (0,

Actually, the last equation is very sufficient in the Gibbs measure theory. In the present paper, we give short information about relation between the last equation and the Gibbs measures. Detailed information can be found in references [5,8].

Let rk = (V, L) be a Cayley tree and we consider models where the spin takes values in the set [0,1], and is assigned to the vertices of the tree. For A c V a configuration aA on A is an arbitrary function aA : A ^ [0,1]. For J e R \ {0} and £ : (u, v) e [0,1]2 ^ £uv e R we define a Hamiltonian

H(a) = -J £a(x)a(v), (2.6)

(x,y)£L

which is the generalization of some classic models such as Ising, Potts,.... Here and below, the notation (x, y) means the nearest neighbor vertices.

Let h(t, x) := ht,x, (t, x) e [0,1] x V be real mapping of x e V \ {x0}, where x0 is a root of the tree. Now, we define the family of probability measures {^(n)}n&N on QVn, i.e.,

V(n)(an) = Z- exp ( -pH(an) + ^ K{x),x) , (2.7)

V xeWn J

Here, as before, an : x e Vn ^ a(x) and Zn is the partition function.

If the family of probability measures {^(n)}neN is compatible then by Kolmogorov's extension theorem, there exists a unique measure p on QV such that p ({a|Vn = an}) = p(n)(an) for all n e N. The measure p is called splitting Gibbs measure corresponding to Hamiltonian (2.6)

The following result describes a new condition which is equivalent to the condition of compatibility.

Proposition 1. [5] The probability distributions ^(n)(&n), n = 1, 2,..., in (2.7) are compatible iff for any x G V \ {x0} the following equation holds:

Z01 exp( J^tu)/ (u,y)du

f (t,x) = JJ ±

yes(x) Jo exP(J0£o«)f (u,y)du

Here, and below f (t, x) = exp(ht,x - h0,x), t G [0,1] and du = A(dw) is the Lebesgue measure.

(2.8)

Let / (t, x) do not depend on the vertices of the Cayley tree. Then the last equation (2.8) has strongly positive solution if and only if the equation (2.5) has strongly positive solution in M0 = {/ G C+[0,1] : / (0) = 1}, where C +[0,1] is the set of all positive continuous functions on [0,1] (see [6]).

Let / ( 1 ) = g ( 2 ) = c and denote that

= <

f (t) if t G

c if t G g(t) if t G

»■3

1 2 3 ' 3

2, i 3

and (t) = <

g(1 - t) if t G

c if t G

f (1 - t) if t G

o, 1

3 1 2 3 ' 3

2 ■ 1

Also, for f 1 ( 3 j = g1 ( 2 ) = c1, we define functions:

^i(u)= <

fi(u) if u G ci if u G gi(u) if u G

»■ 3

1 2

3 ' 3

2, 1 3

and ^2 (u) = <

gi(1 — u), if u G ci , if u G

fi(1 — u), if u G

», 1 ' 3 1 2 3 ' 3

'2 • 1

By using these functions, we define a degenerate kernel:

KT(t, u) = (t)^i (u) + y>2(t)^2(u).

Then equation (2.5) can be written as

(2.9)

f(t) = J tei(t)Vi(u) + ^(t)^(«))fk(u)du.

0

1 1 f (t) = Vi(«)f k(u)d(u) + ^(t) y v2(u)fk(u)d(u) = f (t).

00

1 1 Ci = J V>i(u)fk(u)d(u), C2 = J ^(«)fk(u)du. 00 Consequently, by taking into account f (t) = C1^1 (t) + C2y2(t), we obtain

(2.10)

Namely,

Put

i

Ci = J ^¿(u)(Cm(u) + C2^2(u))fcdu, i G {1, 2}.

The last equality is equivalent to

Ci =

C2

k o k o

Cf ai + C2k ai +

CÎ iC2a2 +.... +

C| iCia2

+ .... +

C|afc+i

Ci afc+i.

(2.11)

l

Here and below,

1 1 1

ai = j h (u)<fk (u)du, a2 = J h(u)^kL~1 (u)<fi2(u)du, ... a.k+1 = j h(u)<fk+1(u)du, 0 0 0 1 1 1

P1 = J Mu)^ (u)du, P2 = j Mu)¥k-1(u)¥2 (u)du, ... Pk+1 = J Mu)fk+1(u)du. (2.12)

C

Let x = —, then the system of equations (2.11) can be written as

C2

)atxk + I | a2xk-1 + I | a3xk-2 + ... + ( | akx + ( I ak+i

W V2 J_vfc-w vfc

k I I k I I k I I k I I k I

I ak+ixk + I I akxk 1 + I I ak-ixk 2 + ... + I I a2x + I I ai

0 \ 1 \ 2 \ k - 1 V k

Namely,

Qk (x) := | k J ak+1Xk+1 + II k J ak - I ^ | a1 | xk + ...

... + ( ( k J a-1 - ( k 1 J ak J x - I k I ak+1. (2.13)

Now, we can conclude the following result:

Proposition 2. Finding positive solutions of equation (2.5) with the kernel (2.9) is equivalent to finding positive roots of the polynomial Qk (x).

In the next sections, we consider the positive roots of the polynomial Qk (x) in special cases.

3. Solutions to the equation (2.5) with the kernel (2.9) for the case k = 2.

In this section, we study positive solutions to equation (2.5) with the kernel (2.9) for the case k = 2. Let k = 2 then the system of equations (2.11) can be rewritten as

C1 = of a1 + 2C1 C2a2 + C2a3

Cf = Cfp1 + 2C1C2P2 + Cfp3. By (2.12) and the construction of kernel one obtains:

(3.1)

a.1 = j h(u)<f1(u)d(u) = J f1(u)f2(u)d(u) + J C1 c2d(u) + j g1(u)g2(u)du = 0 0 12

3 3

1 2 . .

3 31 1

= J g 1(1 - u)g2(l - u)d(u) + J C1C2d(u) + J f1(l - u)f2(l - u)du = j 1f(u)^f(u)d(u) = P3. 0 12 0

3 3

Similarly,

1 1 2 <

1 3 31

a2 = j I1(u)f1(u)<f2(u)d(u) = j f1 (u)f (u)g(l - u)d(u) + J C1C2d(u) + J g1 (u)g(u)f (1 - u)du = 0 0 12

3 3

1 1

= j g 1(1 - u)f (u)g(l - u)d(u)+ J C1C2d(u)+ J f1(1 - u)g(u)f (1 - u)d(u) = J ^2(u)^1 (u)<p2(u)d(u) = P2

1

1

3

3

3

3

3

and

1 1 2 -i

J- 3 3 i

«3 = j (u)d(u) = j /i(u)g2(1 — u)d(u) + j cic2d(u) + j gi(u)/2(1 — u)d(u) =

0 0 12

3 3

1 2 1 1 3 3 1 1

= j gi(1 — u)/2(u)d(u) + j cjc2d(u) + j /-(1 — u)g2(u)d(u) = j —2(u)^>2(u)d(u) = ßi.

Hence,

ai = ß3, "2 = ß2, «3 = ßi. (3.2)

C

By taking into account x = — and the equality (2.13), we have

C2

Q2(x) = «3x3 + (2«2 — a-)x2 + (a- — 2«2 )x — «3 = 0.

Thus, we obtain

Q2(x) = (x — 1)(a3x2 + (03 + 2a2 — a-)x + 0,3) = 0.

Put

D := (202 — 0- — 03)(202 — 0- + 303). It is easy to check 0- + 03 > 202. Therefore, the sign of D is the same as the sign of the expression 0- — 202 — 303.

Proposition 3. Let k = 2 and KK(t, u) be the kernel which defined in (2.9). Then the following statements hold:

(1) If 0- < 202 + 303 then there is unique positive solution of (2.10);

(2) If 0- = 202 + 303 then there are exactly two positive solutions of (2.10);

(3) If 0- > 202 + 303 then there are exactly three positive solutions of (2.10). In the language of Gibbs measure theory, we obtain:

Theorem 1. Let k = 2 and Kf(t, u) be the function of the Hamiltonian (2.6). Then the following assignments hold:

(1) If 0- < 202 + 303 then there is unique translation-invariant Gibbs measure of the model (2.6);

(2) If 0- = 202 + 303 then there are exactly two translation-invariant Gibbs measures of the model (2.6);

(3) If 0- > 202 + 303 then there are exactly three translation-invariant Gibbs measures of the model (2.6).

4. Solutions to the equation (2.5) with the kernel (2.9) for the case k = 3.

Now, we consider positive solutions to equation (2.5) with the kernel (2.9) for the case k = 3. Then the system of equations (2.11) can be written as

C- = C3«- + 3C2C2a2 + 3C-C22a3 + C|a4

(4.1)

C = C3ßi + 3C2C2ß2 + 3CiC2ß3 + C23 ß4

and

ai = j — i (w)<£>3(w)d(w) = j /i(u)/3(u)du + j c ic3du + j gi (u)g3(u)du = ß4.

0 0 12

3 3

Analogously,

i 1 2 i

i 3 3 i

«2 = j — i(u)^2(u)^2(u)du = j /i(u)/2(u)g(1 — u)du + j c ic3du + j g i(u)g2(u)/(1 — u)du = ß3,

i 3 3 i

3

«3 = j — i (u)<^> i (u)^2(u)du = j /i (u)/(u)g2(1 — u)du + j c ic3du + j g i (u)g(u)/2(1 — u)du = ß2,

0 0 12

3 3

i 1 2 i

i 3 3 i

«4 = j — i(u)^|(u)du = j /i(u)g3(1 — u)du + j c ic3du + j g i(u)/3(1 — u)du = ßi.

0 0 12

3

3

1

3

3

As a result, one obtains

ai = $4, «2 = $3, «3 = $2, «4 = (4.2)

Consequently, the polynomial, which is defined in (2.13), is

Q3(x) = «4x4 + (3«3 — «i)x3 + («i — 3«3 )x — «4 = 0.

Namely,

Q3(x) = (x — 1)(x + 1)( «4x2 + (3«3 — «i)x + «4) = 0. The discriminant of the polynomial «4x2 + (3«3 — «1)x + «4 is

D := (3«3 — «i — 2«4)(3«3 — «i + 2«4)-

Proposition 4. Let k = 3 and KK(t, u) be the kernel which defined in (2.9). Then the following statements hold:

(1) If |3«3 — «i| > 2«4 then there is unique positive solution of (2.10);

(2) If |3«3 — «i| = 2«4 then there are exactly two positive solutions of (2.10);

(3) If |3«3 — «i| < 2«4 then there are exactly three positive solutions of (2.10). In the language of Gibbs measure theory, we can conclude the following theorem.

Theorem 2. Put k = 3 and Kf(t, u) is the function of the Hamiltonian (2.6). Then the following assignments hold:

(1) If |3«3 — «i| > 2«4 then there is unique translation-invariant Gibbs measure of the model (2.6);

(2) If |3«3 — «i | = 2«4 then there are exactly two translation-invariant Gibbs measures of the model (2.6);

(3) If |3«3 — «i | < 2«4 then there are exactly three translation-invariant Gibbs measures of the model (2.6).

5. Solutions to the equation (2.5) with the kernel (2.9) for the case k = 4.

In the Gibbs measure theory, it is more interesting to construct Hamiltonian which has more than three translationinvariant Gibbs measures. In this section, we show that the Hamiltonian which corresponds to our kernel has more than three Gibbs measures. For the case k = 4, the system of equations (2.11) can be written as

Ci = C4 «i + 4C3C2«2 + 6C12C22«3 + 4CiC|«4 + C^s C2 = C4 + 4C3C2$2 + 6C2C22,03 + 4CiC23$4 + C24$s

and

i 3 3 i

ai = j фi(u)<^>4(u)du = j /-(u)/4(u)du + j c -c4du + j gi(u)g4(u)du = вб,

«2 = j Ф i (u)^2(u)du = j /i (u)/3(u)g(1 — u)du + j c ic4du + j g i (u)g3(u)/(1 — u)du = в4,

0 0 12

3 3

1 1 2-, i 3 3 i

a3 = j ф i(м)^2(м)^2(м)йм = j /i(u)/2(u)g2(1 — u)du + j c -c4du + j g i(u)g2(u)/2(1 — u)du = вз,

i 3 3 i

«4 = j ф - (u)<^> - (u)y>2(u)du = j /i (u)/(u)g3(1 — u)du + j c -c4du + j g - (u)g(u)/3(1 — u)du = в2, 0 0 12

i 3

«5 = j фi(u)du = j / (u)g4(1 — u)du + j c ic4du + j g i(u)/4(1 — u)du = вь

0 0 12

3 3

Hence,

a l = в5, «2 = в4, «з = вз, «4 = в2, «5 = в 1- (5.1)

Ci

By taking into account x = — and the equality (5.1), one obtains

C2

Q4(x) = «5Х5 + (4«4 — a i )x4 + (6«з — 4«2)x3 + (4«2 — 6«з)х2 + (a i — 4«4 )x — «5 =0.

3

3

1

1

3

3

3

3

1

3

Namely,

Q4(x) = (x — 1) («5x4 + («5 + 4«4 — a-)x3 + («5 + 4«4 — a- + 6«3 — 4«2)x2 + («5 + 4«4 — a-)x + «5) = 0. After short calculations, we have (for x = 1)

Q4(x)

x2(x — 1)

= a5y + (a5 + 4a4 — ai)y + (4a4 — a5 — ai + 6a3 — 4a2),

where y = x + 1. The discriminant of the right hand side of the last equality is

x

D = («5 + 4«4 — ai)2 — 4«5(4«4 — «5 — ai + 6a3 — 4«2). Thus, we can give the following result:

Proposition 5. Let k = 4 and K(t, u) be the kernel which defined in (2.9). Then the following statements hold:

(1) If (a5 + 4a4 — a1 )2 < 4a5(4a4 — a5 — a1 + 6a3 — 4a2) then there is unique positive solution of (2.10);

(2) If (a5 +4a4 — a1)2 = 4a5(4a4 — a5 — a1 +6a3 — 4a2) then there are at most three positive solutions of (2.10);

(3) If (a5 + 4a4 — a1)2 > 4a5(4a4 — a5 — a1 + 6a3 — 4a2) then there are at most five positive solutions of (2.10).

In the Gibbs measure theory, we obtain the following theorem.

Theorem 3. Put k = 4 and let K (t, u) be the function of the Hamiltonian (2.6). Then the following assignments hold:

(1) If (a5 +4a4 — a1)2 < 4a5(4a4 — a5 — a1 + 6a3 — 4a2) then there is unique translation-invariant Gibbs measure of the model (2.6);

(2) If (a5 +4a4 — a1)2 = 4a5(4a4 — a5 — a1 + 6a3 — 4a2) then there are at most three translation-invariant Gibbs measures of the model (2.6);

(3) If (a5 + 4a4 — a1)2 > 4a5(4a4 — a5 — a1 + 6a3 — 4a2) then there are at most five translation-invariant Gibbs measures of the model (2.6).

6. Conclusion

The central concept in understanding phase transitions on the Cayley trees is the notion of cluster decomposition. In the absence of phase transitions, there is a unique Gibbs measure that describes the equilibrium behavior of the system. However, in the presence of a phase transition, multiple Gibbs measures can coexist. The existence of multiple Gibbs measures is often associated with the breaking of symmetry in the system. This symmetry breaking can manifest itself in different ways, such as the appearance of multiple stable states or the coexistence of different phases (e.g. [3]).

Mathematically, the presence of multiple Gibbs measures can be established by studying the behavior of the system under different boundary conditions or by considering the behavior of relevant observables. The occurrence of non-uniqueness in the limit of large system size indicates the existence of phase transitions and the coexistence of multiple Gibbs measures. Also, the properties of the different Gibbs measures can also provide insights into the nature of the phase transition. For example, the entropy of the Gibbs measures can exhibit non-analytic behavior or the correlation length may diverge at the critical point (see [2] and [4]).

In short, the existence of multiple Gibbs measures is closely related to the presence of phase transitions on the Cayley trees. These phase transitions are characterized by breaking of symmetry and the coexistence of different equilibrium states. The identification and characterization of different Gibbs measures play a crucial role in understanding the critical phenomenon and the behavior of the system near the phase transition point. In this paper, we showed that under certain conditions (see Theorems 1,2,3) there exist multiple Gibbs measures of the model (2.6) (this model is a generalization of Ising, Potts, SOS,...) on the Cayley trees of order two, three and four.

Now, we will give short information about the novelty of our paper. Firstly, in [5], authors consider the model (2.6) on the Cayley trees and prove that there is not any phase transition on the Cayley tree of order k = 1. For the case k > 2, models in which phase transitions exist are constructed in [8] and the generalization of the constructed model is studied in [7,9]. In [6], the problem of finding translation-invariant (periodic with period two) Gibbs measures of the model (2.6) was reduced to finding positive fixed points of the nonlinear operator of Hammerstein type and the problem of the existence of fixed points of this operator considered in [10,15,16] ([13,14]). But, from the Gibbs measure theory, it is interesting to study fixed points of the nonlinear operator of the Hammerstein type with degenerate kernel and all the above works corresponding to non-degenerate kernels. Then works for degenerate kernels considered in [11] and [12] but in these papers, the fixed point of the operator with the degenerate kernel does not correspond to any Gibbs measure. In the present paper, each fixed point of the operator with a degenerate kernel is corresponding to one Gibbs measure.

References

[1] Dolezale V. Monotone operators and its applications in automation and network theory, in: Studies in Automation and Control. Elsevier Science Publ, New York, 1979.

[2] Friedli S., Velenik Y. Statistical Mechanics of Lattice Systems. Cambridge University Press, 2017.

[3] Georgii H.O. Gibbs Measures and Phase Transitions. de Gruyter Studies in Mathematics, 2011.

[4] Rozikov U.A. Gibbs measures on Cayley trees, World Sci. Pub, Singapore, 2013.

[5] Rozikov U.A. and Eshkabilov Yu.Kh. On models with uncountable set of spin values on a Cayley tree: Integral equations. Math. Phys. Anal. Geom., 2010,13, P. 275-286.

[6] Eshkabilov Yu.Kh., Haydarov F.H., Rozikov U.A. Uniqueness of Gibbs measure for models with uncountable set of spin values on a Cayley tree. Math. Phys. Anal. Geom., 2013, 16(1), P. 1-17.

[7] Botirov G., Jahnel B. Phase transitions for a model with uncountable spin space on the Cayley tree: the general case. Positivity, 2019, 23, P. 291301.

[8] Eshkabilov Yu.Kh., Haydarov F.H., Rozikov U.A. Non-uniqueness of Gibbs measure for models with uncountable set of spin values on a Cayley Tree. J.Stat.Phys., 2012,147, P. 779-794.

[9] Jahnel B., Christof K., Botirov G. Phase transition and critical value of nearest-neighbor system with uncountable local state space on Cayley tree. Math. Phys. Anal. Geom., 2014, 17, P. 323-331.

[10] Eshkabilov Yu.Kh., Haydarov F.H. On positive solutions of the homogeneous Hammerstein integral equation. Nanosystems: Phys. Chem. Math., 2015, 6(5), P. 618-627.

[11] Eshkabilov Yu.Kh., Nodirov Sh.D., Haydarov F.H. Positive fixed points of quadratic operators and Gibbs Measures. Positivity, 2016.

[12] Eshkabilov Yu.Kh., Nodirov Sh.D. Positive Fixed Points of Cubic Operators on R2 and Gibbs Measures, Jour.Sib.Fed. Univer.Math.Phys., 2019, 12(6), P. 663-673.

[13] Haydarov F.H. Fixed points of Lyapunov integral operators and Gibbs measures. Positivity, 2018, 22(4), P. 1165-1172.

[14] Eshkabilov Yu.Kh., Haydarov F.H. Lyapunov operator L with degenerate kernel and Gibbs measures. Nanosystems: Phys. Chem. Math., 2017, 8(5), P. 553-558.

[15] Haydarov F.H. Existence and uniqueness of fixed points of an integral operator of Hammerstein type. Theor. Math. Phys., 2021, 208, P. 1228-1238.

[16] Ganikhodjaev N.N. Exact solution of an Ising model on the Cayley tree with competing ternary and binary interactions. Theor. Math. Phys., 2002, 130, P. 419-424.

Submitted 7 November 2023; revised 12 January 2024; accepted 14 January 2024

Information about the authors:

IsmoilM. Mavlonov - National University of Uzbekistan, University str., 4 Olmazor district, Tashkent, 100174, Uzbekistan; ORCID 0009-0009-4889-1889; [email protected]

AloberdiM. Sattarov - University of Business and Science, Namangan region, Namangan city, 111 Beshkapa str., Namangan, 160100, Uzbekistan; ORCID 0009-0008-0115-3937; [email protected]

Sevinchbonu A. Karimova - National University of Uzbekistan, University str., 4 Olmazor district, Tashkent, 100174, Uzbekistan; ORCID 0009-0006-3128-6605;

FarhodH. Haydarov - Institute of Mathematics, 9, University str., Tashkent, 100174, Uzbekistan; New Uzbekistan University, 54, Mustaqillik Ave., Tashkent, 100007, Uzbekistan; National University of Uzbekistan, University str., 4 Olmazor district, Tashkent, 100174, Uzbekistan; ORCID 0000-0001-9388-122X; haydarov [email protected]

Conflict of interest: the authors declare no conflict of interest.