The κ-Deformed Calogero–Leyvraz Lagrangians and Applications to Integrable Dynamical Systems

The Calogero–Leyvraz Lagrangian framework, associated with the dynamics of a charged particle moving in a plane under the combined influence of a magnetic field as well as a frictional force, proposed by Calogero and Leyvraz, has some special features. It is endowed with a Shannon “entropic” type kinetic energy term. In this paper, we carry out the constructions of the 2D Lotka–Volterra replicator equations and the N=2 Relativistic Toda lattice systems using this class of Lagrangians. We take advantage of the special structure of the kinetic term and deform the kinetic energy term of the Calogero–Leyvraz Lagrangians using the κ-deformed logarithm as proposed by Kaniadakis and Tsallis. This method yields the new construction of the κ-deformed Lotka–Volterra replicator and relativistic Toda lattice equations.


Introduction
Recently, Calogero and Leyvraz [1,2] demonstrated a time-independent Hamiltonian description of the motion of a charged particle moving in a two dimensional space under the influence of a magnetic field perpendicular to the plane of motion and a frictional force proportional to the velocity. This motion may be viewed as a dynamics of cyclotron motion with friction; this model arises through the coupling of a particle to a large number of external degrees of freedom. The most interesting feature of this Lagrangian is the kinetic energy term-entropic type kinetic energy. This makes the Lagrangian a nonstandard one. We have explored the applications of this class of Lagrangians in our earlier papers [3,4].
It is worth noting that some physical systems cannot be described by the Boltzmann-Gibbs statistical mechanics, for example, systems such as long-range interactions, long-time memory and multifractal or hierarchical structures are some of them. To overcome at least some of these difficulties, Tsallis [5,6] proposed a generalized entropic form based on a κ-deformed logarithm. Later, an example of self-dual κ-deformed logarithmic functions is found in the work of Kaniadakis [7][8][9][10][11]. Over the last decade or so, scientists have observed that many physical and social phenomena often follow the so-called power law distributions (see for example, [12][13][14][15]). We demonstrated that many (generalized) power law distribution equations can be derived from Calogero-Leyvraz Lagrangian formalism using κ-deformation theory.
These are the two popular ways to deform logarithmic and exponential terms in physics. We propose a new Lagrangian where the logarithm term appearing in Calogero-Leyvraz is replaced by a deformed logarithmic function, and study the dynamics. The algebraic structures arising in this κ-deformed framework have been carefully grafted by Scarfone in [16]; in fact, the concept of generalized algebras has been employed constructively to study entropic forms in [16,17]. It is worth noting that the generalized entropies [18] play an important role in generalized distribution theory in complex systems [19] and they have been studied extensively from the information geometric point of view in [20,21]. The aim of this work is to carry out the formulation of the κ-deformed well-known dynamical systems, namely, the Lotka-Volterra replicator equation and N = 2 relativistic Toda lattice system.
Among ecological models, the Lotka-Volterra equation for predator-prey systems [22][23][24] has played a significant role in dynamical systems. It is known that the solutions to this conservative system in phase space are level curves of the energy function. In the mathematical investigations of ecological models, conservative dynamics are often very useful from the (geometrical) mechanics point of view. It is known that the replicator equation in the evolutionary game theory [25] is closely related to the Lotka-Volterra equation. The replicator equation is the first and most important game dynamics studied in connection with evolutionary game theory. It was originally developed for symmetric games with finitely many strategies. Evolutionary game theory [26] studies the behavior of large populations of agents who repeatedly engage in strategic interactions. Note that changes in behavior in these populations are driven either by natural selection via differences in birth and death rates, or by the application of myopic decision rules by individual agents.
The Hamiltonian of the two-dimensional motion of electrons in the presence of the periodic potential and the magnetic field perpendicular to the two-dimensional plane is described by: The spectrum of this system yields a butterfly like structure, known as Hofstadter's butterfly [27]. In a completely independent line of research, the string community investigated a system associated to Hamiltonian, when q and p are restricted to be purely imaginary; this equation reduces to Hofstadter's Hamiltonian (1). In general, q and p are complex coordinates, hence the equation determines a real two-dimensional Riemann surface, or equivalently a complex one-dimensional curve, whose shape is parameterized by the value of H [28,29]. This appears when mirror symmetry is applied to a non-compact Calabi-Yau geometry known as the local P 1 × P 1 geometry. This curve is connected to the Seiberg-Witten curve, encoding the information on instantons in N = 2 supersymmetric pure SU(2) gauge theory [30][31][32].
The system known nowadays as the relativistic Toda lattice (RTL) was invented by S.N.M. Ruijsenaars-the Hamiltonian of the periodic relativistic Toda lattice with just N = 2 particle, after removing the center-of-mass mode. We can illustrate this as follows [28,29]. The Ruijsenaars Hamiltonian [33] of N = 2 quantum Toda lattice is given by: Let us consider center of mass frame p 1 + p 2 = 0 and define which yields an one parameter family of (2), Results of the paper The Calogero-Leyvraz Lagrangian has some interesting features; it is endowed with the Shannon entropic [34] kinetic energy term. The Legendre transform of this Lagrangian yields a Hamiltonian with the exponential momentum term. We have seen that the Calogero-Leyvraz Lagrangian/Hamiltonian allows us to formulate several features related to deformed dynamical systems, balanced loss-gain systems and generalized rate equations [3,4]. In particular, most of the power law distributions and rate equations can be manufactured from this class of Lagrangian. It has been explored that the Calogero-Leyvraz theory of cyclotron-friction motion is closely related to the "curl force" theory as proposed by Berry and Shukla [35,36], although the latter is a totally position-dependent nonconservative force with a nonvanishing curl, whereas the former is totally velocity dependent.
In this paper, we present a different formulation of the celebrated Lotka-Volterra equation [22,24] using the Calogero-Leyvraz Lagrangian. We also give a new derivation of the replicator equation using the Calogero-Leyvraz type Lagrangian. The final example is related to the N = 2 relativistic Toda lattice system. We formulate the latter equation using a different type of entropic Lagrangian; this entropic kinetic energy is described via a cross-entropy term, which yields a new formulation of the N = 2 relativistic Toda lattice system.
Since the kinetic energy term of the Calogero-Leyvraz Lagrangian involves a logarithmic term, we deform this logarithmic term using the Kaniadakis method [17] and obtain κ-deformed Lotka-Volterra, replicator and N = 2 relativistic Toda Lattice equations. In fact, the entire reason to shift the usual formalism to Calogero-Leyvraz formalism is to formulate κ-deformed integrable models. We also formulate these deformed equations using the Tsallis logarithm.
In the introductory section, we review a κ -deformed Liénard equation which satisfies the Chiellini integrability condition. This condition allows us to integrate the Liénard type equation using the Abel equation of the first kind. In general, the Liénard equation does not give a Lagrangian formulation, but with the imposition of the Chiellini condition, it yields a Lagrangian formulation.
This paper is organized as follows. In Section 2, we review the Calogero-Leyvraz Lagrangians and Hamiltonians and their applications. In particular, we also describe the κ-deformed Liénard type equation using Kaniadakis and Tsallis type deformation of the kinetic energy term and demonstrate that this Liénard equation admits the Chiellini integrability condition [37]. This integrability condition plays an important role in the formulation of Lagragian and solutions of the integrable class of the Liénard equation. A nonexhaustive list of applications includes, among others, those in [38][39][40]. After giving a gentle introduction to the Calogero-Leyvraz method, we apply this scheme for the construction of the Lotka-Volterra, replicator and N = 2 relativistic Toda lattice systems in Section 3. Our Section 4 is dedicated to the construction of κ-deformed equations. We give a formulation of the deformed Lotka-Volterra, replicator and N = 2 relativistic Toda lattice equations using the Kaniadakis and Tsallis methods.

Review of Calogero-Leyvraz's Lagrangian and Hamiltonian Formulation of the Dynamics of Cyclotron with Friction System
The Hamiltonian of the free particle moving against friction is given by: according to the Newtonian equation of motionz = −ż. The corresponding Lagrangian description of this system is given by: A minor modification of the Hamiltonian H(p, z) = e p + λp c + cz yields the dynamics of a particle moving against friction in a constant force field λ.
Legendre transformation of Calogero-Leyvraz Lagrangian: Let us recall the Calogero-Leyvraz method first; the Lagrangian is given as L = −γq + v ln v, where v =q. The equation of motionq + γq = 0 yields a constant of motion The corresponding momentum Substituting this in the Legendre transformation, thus we obtain the Calogero-Leyvraz Hamiltonian, where we have scaled the momentum to ignore the constant term. Hence we establish the connection between the Calogero-Leyvraz Lagrangian and Hamiltonian via Legendre transformation. This construction can be extended to Lagrangian involving a timedependent coefficient.

Calogero-Leyvraz Hamiltonian and Planar Systems
Calogero and Leyvraz straight-forwardly generalized this motion by complexification to describe motions taking place in a plane. Physically, this is connected to the motion against friction of a charged particle in the presence of a perpendicular constant magnetic field, or a constant electric field lying in that plane, or both these forces. If we set c = γ + iω and go to the complex plane, the following pair of Poisson commuting Hamiltonians are obtained: At first we consider a minor change; coefficients are considered to be time dependent. The two-dimensional Calogero-Leyvraz model is given by the following Hamiltonian: where γ(t), ω(t) are parameters, (x, y) are coordinates and (p x , p y ) are corresponding momenta. Here we note that the potential energy is a linear function of coordinates while the kinetic energy Ψ = e p x cos p y . The Hamiltonian (10) yields the following equations of motion: Calogero and Leyvraz reformulated (11) in a 3-dimensional context by introducing the 3-vector r = (x, y, 0) in the xy-Cartesian plane and the unit vectorẑ = (0, 0, 1) orthogonal to that plane; this yields: We obtain the sister (or mirror) equations of (11) if we consider a different K.E., viz., Φ(p x , p y ) = e p x sin p y , with the same potential energy γx − ωy, this is given by: It is easy to check both (11) and (13). The linear equation can be generalized to a nonlinear equation from the Calogero-Leyvraz Hamiltonian using generalized potential energy. Suppose we consider where φ and ψ are some functions of x and y. This yields The complex kinetic energy is given by: If we consider the (complex) potential energy U c = Γ(t)z, then the equation of motion resulting from the Hamiltonian H = e P + Γz is Suppose Re(P) = p and Re(z) = a(t)x, then the Hamiltonian becomes H R = e p + q and this yieldsẍ = −a(t)ẋ. Note that a nonlinear potential energy with the same K.E. yields nonlinear equations, for example, H 1 = e p + a(t) ln x and H 2 = e p + a(t)x n yield the following equations of motion: viz respectively. Equation (17) admits, for constant Γ, a Lagrangian: with (17) following from the associated Euler-Lagrange equation.

Illustration: Generalized Liénard Equation and the Calogero-Leyvraz Lagrangian
We have seen that the Calogero and Leyvraz construction yields interesting sets of dynamical equations. In this section, we formulate a nonlinear ODE belonging to the Liénard class of equations.
The Liénard type ordinary second order nonlinear differential equation is given by: f and g are two continuously differentiable functions on R. Since the Liénard equation itself is also an autonomous differential equation, the substitution, y = dq dt or q = y(t)dt, leads the Liénard equation to become a first order differential, which belongs to the Abel equation of the second kind. This can also be expressed in terms of the Abel equation of the first kind, which we use later. Let us define the following generalization of the Calogero-Leyvraz Lagrangian: The Euler-Lagrange equation yields: For small values of (q + f (q)), the above equation reduces to the well-known nonstandard Lagrangian,q For different choices of f we get different types of equations. Let f (q) = λq n , then n = 1; this becomes a damped oscillator equation, n = 2; this maps to a second Riccati or modified Emden-Fowler equation. The corresponding Lagrangian is given by:

Deformations of Calogero-Leyvraz Lagrangians and κ-Deformed Oscillator Equations
The most attractive feature of the Calogero-Leyvraz Lagrangian is the involvement of a logarithmic term. We grab this opportunity and deform the (entropic) kinetic energy term. We use primarily Kaniadakis and Tsallis logarithms. At first we deform (22) using the Tsallis logarithm.
Let us introduce the Tsallis logarithm. We assume q > 0 for practical purposes. The Tsallis q-logarithm and q-exponential functions are defined by: where q = 1 and 1
Proof. By direct computation we obtain: Here, we have tacitly used the formula of ln κ x.
yields theq Let us demonstrate this with a couple of examples. Suppose we take κ = 1 2 , then (28) reduces to If we take f (q) = a q 2 , then Equation (30) yields the second-order Riccati ( also known as the modified Emden ) equation,q
This equation is a fractional damped system except for κ = ± n 2 , where n ∈ Z.
The Liénard type ordinary second order nonlinear differential equation can be mapped to the first kind first order Abel differential equation, This Abel equation allows us to find some exact general solutions of the Liénard type equations by using the integrability conditions of the Abel equation.

Lemma 1.
A first kind Abel type differential equation of the form (36) can be exactly integrated if the functions q(x) and p(x) satisfy the condition:

Claim 1. Equation (35) reduces to the Liénard equation for
which satisfies the Chiellini condition.

Entropic Lagrangian and Integrable Class of Systems
In this section, at first we give a new derivation of Lotka-Volterra and replicator equations using Calogero-Leyvraz Lagrangians endowed with the "Shannon entropic" [34] type kinetic energy terms. Then, using the cross entropy type kinetic energy term, we derive the N = 2 relativistic Toda lattice equation.

Calogero-Leyvraz Lagrangian and Lotka-Volterra Equation
We start with the derivation of the Lotka-Volterra equation. Consider the following logarithmic Lagrangian endowed with an entropic kinetic term: The Euler-Lagrange equation yields:q Let us write this equation as a system of first-order equations. Define: Then (41) equation can recasted as: This is a standard form of the celebrated 2D Lotka-Volterra equation in non-dimensionalized form [23].
In the standard formalism of the Lotka-Volterra equation (43), the Hamiltonian is given by: The nonstandard Hamiltonian form, yields Equation (43). The two Hamiltonians can be connected easily through exponential mapping.

Replicator Equation
In 1978, Taylor and Jonker [41] introduced a system of differential equations which were designated later on as the replicator equation. This equation plays an important role in evolutionary game theory. The replicator equation models the frequency evolution of certain strategic behaviors within a biological population. Hofbauer [25] unveiled an equivalence relation between the Lotka-Volterra equation and the replicator equation.
Consider the first population where individuals interact with each other according to a set of n + 1 pure strategies E 0 , · · · , E n with relative frequencies x 0 , · · · , x n , and the second population plays different m + 1 pure strategies F 0 , · · · , F m with frequencies y 0 , · · · , y m . After a contest E i versus F j , the payoff for the first player is a ij whereas for the second player it is b ji . Let A = (a ij ) be the matrix consisting of these a ij and so also B, then for such games the evolutionary dynamics is given by: Adding or multiplying a (positive) constant to each column of A or B does not alter the dynamics.
In the case of n = m = 1, the above equations simplify to: We will use the Calogero-Leyvraz type Lagrangian to derive the planar replicator type equation.

Calogero-Leyvraz Lagrangian and Replicator Equation
We define the Lagrangian of the coupled system as: It is straight forward to see that the Lagrangian (48) yields: Let us define:q We obtain the replicator equation from (49):
The Ruijsenaars Hamiltonian is given by: where q = e iRh .
For N = 2 case this equation reduces to: Consider the centre of mass frame p 1 + p 2 = 0. Let us define: We express H 2 as:Ĥ for practical purposes we scaled R = 1. Let us express Hamiltonian (57) in terms of cosine hyperbolic function where we drop the factor 2. The Hamiltonian equation yields: which reduces to Equation (53).

Connection to Calabi-Yau Manifold
We must note that the energy function E = e p + e −p + e q + e −q can be expressed as: This defines a genus one Riemann surface. The complex 3D space V = X + X −1 + Y + Y −1 − E describes a Calabi-Yau manifold. This sets up a connection with the Calabi-Yau manifold. The Riemann surface has enough information to describe this Calabi-Yau manifold. The energy function E(p, q) considered to be Hamiltonian appears in the string theory.
Mirror symmetry states that a CY manifold has its mirror dual. The Kähler structure of the original CY is mapped to the complex structure of the mirror CY, and vice versa. In our case, the mirror curve is given by: where µ 1 and µ 2 are the complex moduli of the mirror CY.
The new equation is in the same form as the Lagrangian of the N = 2 relativistic Toda lattice equation. Our case is similar to the case of the quantized mirror curve for the local P 1 × P 1 , which is related to the quantum eigenvalue problem of the relativistic Toda lattice with just two particles.

Calogero-Leyvraz Type Lagrangian with Coupling Constant and Mirror Map
Consider the following map: where γ −1 is a coefficient ( or coupling ) parameter. The corresponding Hamiltonian or energy function is given as: This yields: Let us change the variable,

The κ-Deformed 2D Lotka-Volterra, Replicator and Relativistic Toda Lattice Equations
By deforming the natural logarithm and exponential functions we present κ-deformed 2D Lotka-Volterra and relativistic Toda lattice equations in this section. We will use both Kaniadakis and Tsallis deformations to derive new sets of κ-deformed systems.

The κ Deformation of 2D Lotka-Volterra Equation
We describe two types of κ-deformed systems, semi-deformation and full deformation. In the first case we only deformed the kinetic part, whereas in the second case we consider both the kinetic energy (K.E.) and potential energy (P.E.) parts of the deformations.
Case 1: (deforming only K.E) Consider the following κ-deformed Lagrangian: Proposition 4. The Euler-Lagrange equation of the deformed Lagrangian (65) reads: Proof. It is easy to see that: Thus we obtain: where the generalized κ deformed exponential is given by The final result follows from the remaining part of the calculation.
It is clear that when κ → 0 we recover the ordinary Lotka-Volterra equation.
Case 2: (deforming both K.E. and P.E.) In this case we also deformed the exponential term e −q in the potential. We take the following Lagrangian: Thus we obtain the following result from the straightforward computation.
Proposition 5. The Euler-Lagrange equation of the deformed Lagrangian L d κ yields:

Expressing κ-Deformed Equation
The inverse of the generalized κ deformed exponential Exp κ (x) is given by: We now express Equation (66) in a standard form. Let us define: Thus we obtain:ẏ This can be expressed as:ẏ We recover the original equation when κ → 0. Let w = ln z or z = e w . Then the second equation becomeṡ A further change of variable p = 1 2 w 2 yields a modified set of deformed Lotka-Volterra equations:q Hence we express the deformed Lotka-Volterra equation in a standard form using the generalized κ deformed exponential function.

Tsallis Logarithm and Deformed Lotka-Volterra System
In section we express the Tsallis logarithm and exponential in terms of κ, which are given as: We now deform the Lotka-Volterra Lagrangian using the κ-deformed Tsallis logarithm and exponential. It is defined as: with where 1 − κq > 0.
Proof. Using the properties of the Tsallis logarithm and exponential functions we arrive at our desired result.
One can readily check that, when κ → 0, Equation (79) reduces to the usual Lotka-Volterra equation. If we assume only the deformation of the kinetic term using the Tsallis logarithm, then Equation (79) reduces to:

The κ-Deformed Replicator Equation
We can deform the Lagrangian of the replicator Equation (48) using the Kaniadakis or Tsallis deformation of logarithm term. Using Kaniadakis deformation we obtain: (80) Proposition 7. With the Euler-Lagrange equation related to the deformed Lagrangian (80), we obtain the Kaniadakis deformed coupled equation: Corollary 3. The Kaniadakis κ deformed replicator equation is given by: where x =q 1 and y =q 2 .
One can readily see when κ → 0 the deformed replicator Equations (83) and (84) reduces to the original replicator equation.
We can repeat the same procedure using Tsallis deformation of the logarithm. The coupled equations are given by: which leads to the Tsallis deformed replicator equation, This again reduces to the original replicator equation for κ → 0.

The κ-Deformed N = 2 Relativistic Toda Lattice System
In this section, at first we also modify the entropic kinetic energy term, keeping all other terms unchanged. The Kaniadakis κ-deformed Lagrangian for the N = 2 relativistic Toda lattice system is defined as: Proposition 8. The Euler-Lagrange equation corresponding to the κ-deformed Lagrangian L 1κ yields: q 1 +q 2 cosh κ ln (q + 1 +q 2 ) + d dt Proof. This proof follows from the direct computation.
One can readily check that for κ → 0 (89) reduces to the ordinary N = 2 relativistic Toda lattice equation.
For the most general case, we deform the potential term cosh q too, which yields the following Lagrangian: We obtain κ-deformed N = 2 relativistic Toda lattice equation, q 1 +q 2 cosh κ ln (q + 1 +q 2 ) + d dt q 1 +q 2 cosh κ ln (q + 1 +q 2 ) − 1 + 1 This yields the most general Kaniadakis κ-deformation of the N = 2 relativistic Toda lattice equation which reduces to the original one when κ goes to zero.
Tsallis Deformed N = 2 Relativistic Toda Lattice Equation In this section we present the deformation of the the N = 2 relativistic Toda lattice equation using Tsallis deformation. Let the entropic part of the kinetic term be given by: We now compute the equation of motion using the Tsallis deformed kinetic energy.

Outlook
In this paper we considered a special class of Lagrangians proposed by Calogero and Leyvraz with an "exotic" kinetic energy term. This term has a close resemblance to the Shannon entropy function,q lnq. Using this new type of Lagrangian, we derived the celebrated Lotka-Volterra and replicator equation. We then generalized the construction of Calogero and Leyvraz and considered a different type of kinetic energy term based on cross entropy. We then manufactured an N = 2 relativistic Toda lattice system. We also discussed the significance of this equation in modern physics. Different avatars of this equation appeared in string theory and theoretical high energy physics-purely imaginary position and momentum coordinates lead to the Hofstadter model.
The main goal is to express all these celebrated equations in terms of logarithmic kinetic energy using the deformation of the entropic kinetic energy term. We used the Kaniadakis κ-deformed logarithm and exponential functions to deform these Calogero-Leyvraz type Lagrangians to give a new formulation of κ -deformed Lotka-Volterra, replicator equation and N = 2 relativistic Toda lattice system. We also extended this deformation to the Tsallis class and derived Tsallis-deformed equations. All the original equations can be recovered from the deformed systems when κ → 0.
In a nutshell, this paper elucidated the strength of the Calogero-Leyvraz formalism based on entropic kinetic terms. It would be interesting to derive more known and not so well known systems using this method. The predator-prey models are one of the best places to apply our scheme. We may try to apply this scheme to planar generalized Lotka-Volterra (GLV) equations; for example, consider two interacting populations with densities x > O and y > O with the simplest formal description of interaction with the linear dependence of the growth ratesẋ/x andẏ/y. This yields the following GLV equation: x = x a + bx + cy ,ẏ = y d + ex + f y).
We can generalize this construction and check whether we can manufacture this new equation using the Calogero-Leyvraz formalism. We then implement the κ -deformation of such equations and study their dynamics.
It would be worth investigating the Calabi-Yau manifold connected to the κ-deformed Hamiltonian of the N = 2 relativistic Toda lattice equation.