Active and Purely Dissipative Nambu Systems in General Thermostatistical Settings Described by Nonlinear Partial Differential Equations Involving Generalized Entropy Measures

In physics, several attempts have been made to apply the concepts and tools of physics to the life sciences. In this context, a thermostatistic framework for active Nambu systems is proposed. The so-called free energy Fokker–Planck equation approach is used to describe stochastic aspects of active Nambu systems. Different thermostatistic settings are considered that are characterized by appropriately-defined entropy measures, such as the Boltzmann–Gibbs–Shannon entropy and the Tsallis entropy. In general, the free energy Fokker–Planck equations associated with these generalized entropy measures correspond to nonlinear partial differential equations. Irrespective of the entropy-related nonlinearities occurring in these nonlinear partial differential equations, it is shown that semi-analytical solutions for the stationary probability densities of the active Nambu systems can be obtained provided that the pumping mechanisms of the active systems assume the so-called canonical-dissipative form and depend explicitly only on Nambu invariants. Applications are presented both for purely-dissipative and for active systems illustrating that the proposed framework includes as a special case stochastic equilibrium systems.


Introduction
In a seminal work, Nambu generalized classical Hamiltonian mechanics to what is nowadays known as Nambu mechanics [1].Nambu mechanics describes n-dimensional dynamical systems that exhibit n − 1 invariants (i.e., integrals of motion).Applications of Nambu mechanics can be found in various fields of physics (classical mechanics, astrophysics, electrodynamics, solid state physics, hydrodynamics, nonlinear physics) and even in disciplines beyond physics.
In electrodynamics, Nambu mechanics has found several applications [5,11,12,20,21].In particular, the motion of a charged particle in a constant magnetic field can be studied from the perspective of Nambu mechanics [22].In solid state physics, the Calogero-Moser system has been studied from the perspective of Nambu mechanics [23,24].Moreover, Nambu mechanics has been used to calculate energy levels by means of a particular perturbation theoretical method [5].Nambu mechanics has been used to address certain hydrodynamic problems [25,26].In nonlinear physics, the Lorenz system has turned out to be an interesting model for applications of Nambu mechanics.Typically, a modified version of the Lorenz system that is known to be integrable is studied rather than the original chaotic system [11,[27][28][29].Furthermore, it has been shown that Nambu mechanics provides a theoretical framework to examine the dynamics of chiral models [16].
Beyond physics, it has been suggested that Nambu mechanics may be used to describe biochemical reaction equations [30,31].
Only recently, Nambu mechanics has been used to address active systems, such as limit cycle oscillators [19,28,29,[32][33][34].The notion of an active system is closely related to pumping and negative friction [35,36].Pumping and negative friction are dissipative components.Systems exhibiting such dissipative components are likely to be affected by fluctuating forces [37].Therefore, at issue is to discuss active Nambu systems in the context of stochastic processes.A first step in this direction has been conducted within the framework of Boltzmann-Gibbs-Shannon thermostatistics and Fokker-Planck equations that are linear with respect to their probability densities [7].
In the following sections, a more comprehensive picture of stochastic, active Nambu systems will be developed.First, general entropy measures will be considered rather than the special case of the Boltzmann-Gibbs-Shannon entropy.This will lead to Fokker-Planck equations that are nonlinear with respect to their probability densities.From a mathematical perspective, these nonlinear Fokker-Planck equations correspond to nonlinear partial differential equations.Second, non-equilibrium thermodynamic state variables will be defined that are consistent with the aforementioned generalized entropy measures.The thermodynamic state variables can be used to characterize the stochastic, active Nambu systems.Furthermore, it is pointed out that the proposed framework includes purely dissipative equilibrium systems as special cases.In Section 2, the general theoretical framework will be developed.In Section 3, examples of purely dissipative and active Nambu systems will be presented.

Nambu Dynamics: Deterministic Case
We consider an n-dimensional state space described by the state vector r = (x 1 , . . ., x n ) involving n components x k .The state vector is assumed to evolve in time t.Our focus is on Nambu systems [1].Therefore, our departure point is the deterministic evolution of r as defined in Nambu mechanics.Accordingly, in the deterministic case, there are n − 1 functions H 1 , . . ., H n−1 of r, and r(t) satisfies where I = (I 1 , . . ., I n ) denotes a so-called conservative force vector given by: In Equation ( 2), the symbol j 1 ,...,j n denotes the n-dimensional Levi-Civita tensor, which equals one for 1,2,3,4,...,n and changes the sign when two indices are switched.For any other cases (i.e., if two or more than two indices assume the same integer values), the tensor equals zero.Note that ∇I holds, where ∇ = (∂/∂x 1 , • • • , ∂/∂x n ) is the nabla operator.That is, I is divergence free.Importantly, we have ∑ k I k ∂H j /∂r k = 0 for any j, which implies that dH j /dt = 0 holds for all j = 1, . . ., n − 1.That is, the functions H j are invariants of the Nambu dynamics defined by Equation (1).In the context of active, stochastic Nambu systems, we will refer to the functions H j as pseudo-invariants because in the active, stochastic case, they can vary as functions of time.

Nambu Dynamics: Stochastic Case
In earlier work, a stochastic version of Nambu mechanics has been proposed that is suitable to address both active and purely dissipative systems [7,19,32,33].These studies focused on the Boltzmann-Gibbs-Shannon thermostatistics.In order to address this stochastic approach in general thermostatistic settings, in what follows, the free energy Fokker-Planck equation approach is used that has been proposed in the literature [38].
Accordingly, the time-dependent probability density P(r, t) = δ(r − r(t)) of the state r is considered.Here and in what follows, • denotes ensemble averaging, and δ(•) is the Dirac delta function.Following the free energy Fokker-Planck equation approach, three thermodynamic state variables are defined: the thermostatistic entropy S, the non-equilibrium internal energy U NL and the non-equilibrium free energy F NL .The free energy Fokker-Planck equation approach applies to a non-equilibrium active system that operates close to thermal equilibrium, where linear non-equilibrium thermodynamics holds.The approach also includes purely dissipative equilibrium systems as a special case.In this special case, the thermodynamic functions U NL and F NL become equilibrium functions.In order to derive explicit solutions, it is useful to define S by means of the outer function B(y) and the entropy kernel s(P) like [38]: with d n x = ∏ n k=1 d x k , y = s(P) d n x and concave entropy kernels: In the special case B(y) = y and s(P) = −P ln P, we have the Boltzmann-Gibbs-Shannon entropy: Let us put B(y) = y again.For s(P) = (P q − P)/(1 − q) with q = 1, q > 0 and s(P) = −P ln P with q = 1, we have the Tsallis entropy of non-extensive thermostatistics [39][40][41] that reads explicitly: for q = 1 and S T = S BGS for q = 1.Obviously, the Tsallis entropy involves a parameter q and includes the Boltzmann-Gibbs-Shannon entropy as a special case for q = 1.Other entropy measures, such as the Renyi entropy or the Sharma-Mittal entropy, can be obtained as special cases of Equation (3) for appropriate choices of B and s [38].
As such, the internal energy functional U NL can depend in various ways on r and P. In order to arrive at analytical solutions, it has been suggested to follow the so-called canonical-dissipative approach [36,[42][43][44].Accordingly, U NL satisfies two features: it is linear in P, and it depends explicitly only on the functions H j .That is, there is only an implicit dependency on the state r.In summary, we have: For example, in the special case of a Brownian particle evolving in a potential (purely dissipative case) with n = 2, H 1 = H and g(H) = H, we find that U NL reduces to the equilibrium internal energy given by U = H = HPd 2 x, where H is the ordinary Hamiltonian function of the particle (for more details, see below).
Finally, the non-equilibrium free energy F NL is defined in analogy to the equilibrium free energy like: where θ ≥ 0 is a weight parameter that is considered as the counterpart to the temperature T in equilibrium systems.
With the definitions made above, we can propose a model for stochastic Nambu systems in general thermostatistic settings.As we will show below, the model applies both to active and purely-dissipative systems.The model is defined by the free energy Fokker-Planck equation [38]: In Equation ( 9), we have introduced the incomplete nabla operator ∇ D = (a 1 ∂/∂x 1 , • • • , a n ∂/∂x n ) with a j = 0 or a j = 1.If a j = 1 for all j, then ∇ D = ∇.The incomplete nabla operator plays a crucial role for the dissipative part of the dynamics.For this reason, it is denoted by a subindex "D".
In Equation ( 9), the symbol M stands for a semi-positive definite n × n matrix, the so-called mobility matrix [38].Semi-positive definite means that for any vector Z = (Z 1 , . . ., Z n ), we have ZMZ ≥ 0. For example, a diagonal matrix with semi-positive diagonal elements is a semi-positive matrix.Finally, in Equation ( 9), the expression δF NL /δP denotes the variational derivative of the functional F NL [P] with respect to P.
For the special case of the Boltzmann-Gibbs-Shannon entropy (5), by substituting Equations ( 5), ( 7) and (8) into Equation ( 9), we obtain an ordinary Fokker-Planck equation of the form: which is a partial differential equation that is linear with respect to P. However, for entropy measures different from S BGS , Equation ( 9) yields partial differential equations that are nonlinear with respect to P. For example, for stochastic Nambu systems in a non-extensive thermostatistic setting related to the Tsallis entropy (6), by substituting Equations ( 6)-( 8) into Equation ( 9), we obtain: with q > 0. The diffusion term, i.e., the third expression on the right-hand side of the equal sign in Equation (11), is nonlinear with respect to P. This type of nonlinear diffusion coefficient has been introduced by Plastino and Plastino in the context of Fokker-Planck equations associated with the non-extensive Tsallis entropy [45] and is a benchmark nonlinearity of nonlinear diffusion equations in material physics [46][47][48].
Due to the fact that in general, the model ( 9) is nonlinear with respect to P (e.g., see Equation (11)) and in view of the fact that the structure of the model ( 9) is similar to a Fokker-Planck equation, it has been suggested to refer to Equation (9) as the nonlinear Fokker-Planck equation.That is, from a mathematical point of view, the proposed model ( 9) can be regarded as a nonlinear partial differential equation or a nonlinear Fokker-Planck equation.Taking a thermodynamics point of view with a focus on thermodynamic state variables, S, U NL , F NL and, in particular, the role of F NL for the approach to stationarity (see below), Equation ( 9) is regarded as a free energy Fokker-Planck equation.

Approach to Stationarity and Stationary Solutions
Solutions of the ordinary Fokker-Planck equation (i.e., the Fokker-Planck equation that is linear with respect to P) are known to approach under certain conditions stationary solutions [49].Several studies have generalized this result for Fokker-Planck equations that are nonlinear with respect to P and associated with generalized statistical entropy measures different from the Boltzmann-Gibbs-Shannon measure [38,[50][51][52][53][54].In the context of the stochastic Nambu equation model described by the free energy Fokker-Planck Equation (9), the approach to stationarity can be shown conveniently using concepts from the theory of linear non-equilibrium thermodynamics.To this end, following [38], let us introduce the thermodynamic force X th defined by: Note that for generalized entropy measures S, the force X th can depend on P as indicted.Having introduced X th , Equation (9) becomes: Differentiating the free energy F NL with respect to time, from Equations ( 8) and ( 13), we obtain: The proof of Equation ( 14) can be found in [38] and is reviewed in Appendix A. As indicated in Equation ( 14), due to the fact that M is semi-positive definite, it follows that the free energy decays as a function of time or is constant.The decay of the free energy is related to the increase of the statistical entropy due to irreversible processes.The latter will be denoted by d i S. It can be shown (see [38] and Appendix A) that: If F NL is bounded from below [38,55], then Equation ( 14) yields a so-called H-theorem that states that any transient solution in terms of P(r, t) of the stochastic Nambu mechanics system eventually converges to a stationary one.More precisely, from F NL [P] > C for any P and dF NL /dt ≤ 0, we obtain: From Equation ( 14), it then follows that X th (P) = 0. From Equation ( 13), we conclude that: holds for t → ∞.This implies we have either ∂P/∂t = 0 and ∇IP = 0 or there is a time-dependent probability density P(r, t) that satisfies Equation (17).However, Equation ( 17) only allows for the so-called weak solution of the form: These weak solutions are composed of deterministic trajectories.This leads to a contradiction because for θ > 0, the model ( 9) exhibits a diffusion term.That is, trajectories of the model ( 9) exhibit a random component even in the stationary case.Therefore, we conclude that the transient probabilities P(r, t) become stationary in the long-term limit.
Let us derive analytical expressions for the stationary probability densities P of stochastic Nambu models defined by Equation (9).From our previous discussion, it follows that stationary solutions satisfy the conditions X th = 0 and ∇IP = 0.The latter condition implies that: because the conservative force I has zero divergence.We will return to Equation ( 19) below.Let us dwell on the condition X th = 0, which reads: Here, µ is an integration constant.From Equations ( 8) and ( 20), it follows that: Substituting the definitions ( 3) and ( 7) into Equation ( 21) and rearranging terms, we obtain: where f denotes the slope of the entropy kernel, and we have: Since the entropy kernel is assumed to be concave (see Equation ( 4)), the function f is monotonically decaying, and the inverse f −1 of f exists.To obtain a more concise formulation involving a smaller number of symbols, let us denote this inverse f −1 alternatively by [ds/dP] −1 .Then, Equation ( 22) reads: and taking the inverse yields: Equation ( 26) is an implicit definition for the stationary probability density because in general, the expression dB/dy depends on P, as well; see Equation (24).However, if the entropy measure does not exhibit an outer function, that is if B(y) = y holds, then Equation ( 26) becomes an explicit definition for the stationary probability density and reads: Let us illustrate applications of Equation ( 27) for the Boltzmann-Gibbs-Shannon entropy (5) and the Tsallis entropy (6).For S BGS , we have s(P) = −P ln P, ds/dP = −1 − ln P and [ds/dP] −1 (ξ) = exp{−1 − ξ}.Consequently, Equation ( 27) reduces to: with Z = exp{1 − (µ/θ)}.The stationary probability density has the form of a generalized Boltzmann distribution and includes the ordinary Boltzmann distribution for the special case n = 2, H 1 = H, g(H) = H and θ = T.For S T , we have s(P) = (P q − P)/(1 − q), ds/dP = (qP (q−1) − 1)/(1 − q) and [ds 1) for q > 1.Consequently, Equation ( 27) becomes: for q ∈ (0, 1) and: for q > 1.In the expressions above, we have used the operator [•] + , which is defined by [z] + = max(z, 0), where z is a real number.The operator makes sure that probability densities are semi-positive definite.
For q > 1, stationary probability densities of non-extensive systems described by Tsallis entropy frequently correspond to cutoff distributions, that is they decay to zero at certain boundaries and equal zero outside these boundaries.
Finally, let us return to the condition (19).By definition of the canonical-dissipative approach, Equation ( 19) is satisfied for any canonical-dissipative distribution.That is, the canonical-dissipative approach involves kernels g of internal energy functionals that do not depend explicitly on the state vector r, but only implicitly via the invariants H j , see Equation (7).As a consequence of this choice, the stationary distributions only depend on the invariants H j (see Equation ( 26)) and do not explicitly depend on the state vector r.Therefore, the product I∇P vanishes as required by Equation ( 19), which might be shown as follows:

Active Nambu Systems Exhibit Attractors Defined by Classical Nambu Systems
In the "zero temperature" case θ = 0, Equation ( 9) with U NL given by Equation (7) reads: and exhibits a "weak solution" P(u, t) = δ(u − r(t)) , where r(t) are trajectories satisfying a canonical-dissipative Nambu dynamics without noise [19]: Let us show next that the kernel function g of the internal energy U NL can act as a Lyapunov function.Differentiating g with respect to time, we obtain: The first term on the right-hand side of the equal sign vanishes because we have: and ∇H k I = 0 for any k as mentioned in Section 2.1.Consequently, Equation (34) becomes: Let us focus on semi-positive definite mobility matrices that assume a diagonal form: M ij = m j δ ij with m j ≥ 0, where δ ij is the Kronecker symbol.Then, ∇gM∇ D g = ∇ D gM∇ D g holds, which implies that the expression ∇gM∇ D g is of the form ZMZ with Z = ∇ D g, and therefore, we have ∇gM∇ D g = ZMZ ≥ 0. In summary, Equation (36) leads to the inequality: Next, we assume that g is bounded from below.For example, the quadratic function: with parameters A k ≥ 0 and B k is bounded from below with ∀ r : g ≥ C ≥ 0, where C is a constant.From dg/dt ≤ 0 and g ≥ C, it follows that g becomes stationary in the long time limit: Assuming again that M ij = m j δ ij and that ∇ D "matches" the mobility matrix (as in all examples of Section 3), such that m j = 0 ⇔ a j = 0 and m j > 0 ⇔ a j = 1, then Equation (39) implies that ∇ D g = 0; see Equation (36).Substituting this result into Equation ( 33), we see that r satisfies the original Nambu dynamics (1).That is, any active Nambu system converges to a trajectory that is described by a classical Nambu system (1).This trajectory acts as an attractor.
The question arises whether or not the values of the pseudo-invariants H j are fixed by the impact of the function g.If g assumes the quadratic form (38) for a single pseudo-invariant j, with A j > 0, then this question can be answered affirmatively under particular circumstances.
Let us assume that H j is bounded from below and exhibits a minimum value C for a point r or a trajectory r(t), such that H j ≥ C ∧ ∃r : H j (r) = C.Let us assume that Properties (i) and (ii) hold.Then, from Equation (41), it follows that if B j ≥ C, then: That is, if at an initial time point t 0 we have H j = B j , then H j converges to B j in the long time limit.In contrast, if B j < C, then by definition, H j cannot converge to B j .However, in the long time limit, g becomes stationary, and we have ∇ D g = 0. Therefore, from Equation (41), it follows that if B j < C, then: Typically, B j = C denotes a bifurcation point that separates two qualitatively different dynamical domains.If so, from a physics point of view, the parameter B j would be considered as a pumping parameter, whereas from a dynamical systems perspective, B j would be regarded as a bifurcation parameter.

Examples of Active and Purely-Dissipative Systems
In what follows, we will illustrate the general theoretical framework by means of three examples.The first two examples are about purely dissipate systems satisfying Boltzmann-Gibbs-Shannon thermostatistics.The third example describes an active stochastic Nambu system.The system is discussed from the perspectives of classical Boltzmann-Gibbs-Shannon thermostatistics and non-extensive thermostatistics.

Brownian Motion in a Potential Field
Nambu mechanics reduces to Hamiltonian mechanics for n = 2. Therefore, the Brownian motion of a particle of mass m moving in one direction x with velocity v and momentum p = mv is a special case of Nambu mechanics.We put: where V(x) is a globally-attractive potential with min{V} = C and V → ∞ for x → ±∞.In terms of Nambu mechanical systems, we have n = 2, r = (x, p), ∇ = (∂/∂x, ∂/∂p) and H 1 = H.The particle not affected by friction and fluctuating forces satisfies Equation (1) with the conservative force: From Equations ( 1) and ( 45), we obtain the Hamiltonian equations: of a classical particle moving in a potential V(x).The classical particle is assumed to move through a medium.Due to the interaction with the medium, the particle is subjected to friction, on the one hand, and a fluctuating force, on the other hand [49].To account for these two effects, we add dissipative forces in the evolution equation for the moment.Accordingly, the model ( 9) is used with: where γ ≥ 0 denotes the friction coefficient.Likewise, we use ∇ D = (0, ∂/∂p).As far as the thermodynamic state variables are concerned, we put: Since we use S BGS , Equation (10) applies.Substituting Equations ( 47) and (48) into Equation (10) yields: where D denotes the diffusion coefficient defined by D = γθ.The Fokker-Planck Equation ( 49) is also known as Kramers equation.Trajectories of the Brownian particle can be computed from the Langevin equation associated with the Fokker-Planck Equation (49).The Langevin equation reads [49]: where Γ(t) is a fluctuating force defined in terms of a Langevin force normalized with respect to the Dirac delta function like Γ(t)Γ(t ) = 2δ(t − t ) [49].The Langevin equation for the Brownian particle (50) includes Equation ( 46) as special case for γ = 0 (⇒ D = 0).Substituting Equation (48) into Equation ( 28), we obtain the stationary probability density in form of the Boltzmann distribution: for θ = T and Z = exp{−H(x, p)/T} dxdp.

Charged Particle in a Magnetic Field
As pointed out by Pletnev [22] (p.291/292), a charged particle in a magnetic field is an example of a Nambu system.To this end, the particle is described only in the velocity state space with n = 3.For the sake of simplicity, let us put the particle mass equal to unity.The state vector is given by the three-dimensional velocity vector: r = v = (v 1 , v 2 , v 3 ).Accordingly, we have ∇ = (∂/∂v 1 , ∂/∂v 2 , ∂/∂v 3 ).The invariants are: Here, H 1 is the kinetic energy.The parameter q denotes the charge of the particle and should not be confused with the parameter q introduced above in the context of the Tsallis entropy.The vector B denotes the magnetic field vector in three dimensions and is assumed to be constant.The second invariant is the scalar product (dot product) of the particle velocity and the magnetic field vector.This invariant is a consequence of the Lorentz force (produced by the magnetic field), which only acts perpendicular to the particle velocity.Therefore, the Lorentz force cannot change the scalar product vB.In three dimensions, for arbitrary invariants H 1 and H 2 , the conservative force (2) can be expressed in terms of the following cross product [1]: Substituting Equations ( 52) and (53) into Equation (1), we obtain: Since r = v, we re-obtain the classical result from electromagnetodynamics: where F L is the aforementioned Lorentz force.Following our previous example (see Section 3.1), we consider next the motion of a charged particle subject to a magnetic field in a medium (that does not affect the magnetic field).That is, we add two dissipative components to the model ( 55): friction and a fluctuating force.We consider again classical Boltzmann-Gibbs-Shannon statistics.Accordingly, we put ∇ D = ∇ and: where E is the diagonal matrix and γ ≥ 0 is the friction coefficient.Furthermore, we put: Substituting Equations ( 56) and ( 57) into Equation (10) yields: where ∆ = ∇ 2 denotes the three-dimensional Laplace operator and D = γθ again.The Fokker-Planck Equation ( 58) defines the evolution of the probability density P(v, t).Stochastic trajectories v(t) can be computed from the corresponding Langevin equation [49] that reads: where Γ(t) = (Γ 1 (t), Γ 2 (t), Γ 3 (t)) denotes a three-dimensional Langevin force with Γ j (t), Γ k (t ) = 2δ jk δ(t − t ).Here, δ jk denotes the Kronecker symbol.Finally, the stationary probability density can be obtained by substituting Equation (57) into Equation ( 28) and corresponds to a Maxwell distribution: for θ = T.

Active Spinning Top Featuring Non-Extensive Statistics: An Approach Involving Thermodynamic State Variables
So far, we have discussed purely-dissipative systems and stochastic systems within the framework of the classical Boltzmann-Gibbs-Shannon thermostatistics.Let us turn next to an example of an active stochastic Nambu system.Nambu mentioned in his seminal study on Nambu mechanics that the Euler equations of the spinning top can be cast into the form of Equation ( 1).The freely rotating spinning top can be considered as a Nambu system.In two previous studies, a stochastic, active version of a spinning top has been introduced exploiting the Nambu mechanics approach [7,19].This proposed stochastic active spinning top has been discussed from the perspective of classical Boltzmann-Gibbs-Shannon thermostatistics.In what follows, we will introduce a non-extensive variant based on the Tsallis entropy of the stochastic active spinning top proposed in the aforementioned earlier studies [7,19].Moreover, in this context, the thermodynamic state variables F NL , U NL and S will be introduced, and the corresponding thermodynamic perspective will be taken as a departure point.
The dynamics of a spinning top can conveniently be described by means of the three-dimensional angular momentum vector L = (L 1 , L 2 , L 3 ).Let I 1 , I 2 and I 3 denote the principle inertia moments of the top.Then, the evolution equation for L is given by the Euler equation: where we have introduced the diagonal matrix: of "inverse" inertia coefficients.Let us put n = 3 and r = L.Moreover, we introduce H 1 defined by half of the squared amount of the angular momentum and H 2 given by the kinetic energy as invariants: Then, substituting Equation ( 63) together with Equation ( 53) into the Nambu Equation ( 1), we re-obtain the Euler Equation (61).In doing so, it can be shown that the spinning top is a Nambu system.
The Fokker-Planck Equation (65) defines the evolution of the probability density P(L, t).For q = 1, the model reduces to the stochastic version proposed earlier in the literature [7] describing a linear partial differential equation, that is, an ordinary Fokker-Planck equation.In contrast, for q = 1, the model is described by a nonlinear partial differential equation.
In order to highlight the active nature of the model (65), we would like to point out that in the "zero temperature limit" given by θ = 0, we have an H-theorem for g; see Section 2.4.In particular, for A 2 = 0, this implies that in the limiting case t → ∞, the angular momentum L approaches in the amount the value √ 2B 1 , such that H 1 = B 1 .Likewise, for A 1 = 0, it follows that the kinetic energy approaches in the long time limit the parameter value B 2 .In this sense, the model features a pumping mechanisms (i.e., negative damping) that leads to an increase of angular momentum (A 1 > 0, A 2 = 0) or kinetic energy (A 1 = 0, A 2 > 0).
Stochastic trajectories L(t) can be computed by means of [38]: where Γ(t) denotes the three-dimensional Langevin force mentioned in the previous example of Section 3.2.The stochastic evolution Equation (66) can be interpreted in two alternative ways.First, we may consider Equation (66) as part of a two-tiered (or two-layered [38]) description of the stochastic trajectories L(t).Accordingly, first, the Fokker-Planck Equation (65) is solved in the time-interval [t 1 , t 2 ] of interest (tier one).In doing so, P(L, t) is obtained.Subsequently, the solution P(L, t) is substituted into Equation (66).In this case, Equation (66) corresponds to a non-autonomous, multiplicative noise Langevin equation (tier two) interpreted according to Ito calculus.Alternative to this two-tiered interpretation, Equation (66) can be regarded as a self-consistent Langevin equation [38].In line with this alternative interpretation, the probability density P(L, t) is computed from trajectories L(t) like P(u, t) = δ(u − L(t)) .However, the evolution of the trajectories L(t) in turn depends on P(L, t).
For q > 1, we obtain: In order to arrive at a more precise interpretation of the stationary distributions, we may discuss the stationary distribution in appropriately-defined one-dimensional spaces related to the invariants H 1 and H 2 , rather than in the original state space [57][58][59].That is, we define the distributions: for j = 1, 2. In order to reduce the number of symbols, we may consider H j both as a coordinate and a random variable, depending on its context.Accordingly, we replace ξ by H j in Equation (70) and define the distributions of the invariants like: For illustration purposes, let us focus on an active system that features a pumping mechanism that makes the angular moment become asymptotically stable at the amount √ 2B 1 (in the deterministic case θ = 0).That is, we put A 2 = 0. From Equation (71) and Equations ( 67)-(69), we obtain the distributions: for q ∈ (0, 1), for q = 1 with Z = exp{−A 1 (H 1 − B 1 ) 2 /(2θ)}d 3 L and: for q > 1.In Equations ( 72)-( 74), the pre-factor √ H 1 /Z is a re-normalization factor that occurs due to the fact that we transform distributions from a three-dimensional state space to a one-dimensional one.Roughly speaking, we put P 1 (H 1 )dH 1 = P(L)d 3 L and express the three-dimensional state space element d 3 L like d 3  74) can be characterized as follows.For q ∈ (0, 1), the energy distribution is a power law distribution.The probability density of the pseudo-invariant H 1 decays like for sufficiently large values.Since we have |L| 2 = 2H 1 , the probability density of the squared angular moment |L| behaves qualitatively in the same way.For q = 1, the distribution would be a truncated normal distribution (truncated at H 1 = 0) with a maximum at H 1 = B 1 if we ignore the pre-factor √ H 1 /Z .This truncated normal distribution has been studied in detail (both in theoretical and experimental works) in the context of the standard canonical-dissipative limit-cycle oscillator [36,44,57,58,60,61].Taking the pre-factor √ H 1 /Z into account, the truncated normal distribution becomes distorted.For q > 1, the distribution is a cutoff distribution.It is useful to introduce the parameter µ eff = µ + θ/(q − 1), such that Equation (74) reduces to: From Equation (75), it is clear that the shape of the distribution follows an inverted U with a maximum at H 1 = B 1 .The distribution exhibits two cutoff points at then the cutoff points are located at H 1+ = B 1 + 2µ eff /A 1 and H 1− = 0.The distribution assumes finite values in the interval (H 1− , H 1+ ) (which is the so-called "support" of the distribution) and is zero otherwise.

Numerics
In order to illustrate the active nature of the canonical-dissipative spinning top model (65), we solved numerically the corresponding Langevin Equation (66).
We focused on the Boltzmann-Gibbs-Shannon case q = 1 and a pumping mechanism that stabilized the amount of angular moment (A 1 > 0, A 2 = 0).Figures 1 and 2 show one component of L, as well as the pseudo-invariant H 1 as functions of time for two different parameters γ 1 in the deterministic case (zero temperature limit θ = 0 ⇒ D = 0).As expected from our general considerations in Section 2.4, the pseudo-invariant H 1 converged to the stationary value B 1 .The relaxation time towards stationarity became shorter when γ 1 was increased.That is, γ 1 determined the characteristic time scale of the transient period.This interpretation of γ 1 can also be seen from Equation (36).Substituting g = A 1 δH 2 1 /2 with δH 1 = H 1 − B 1 into Equation (36) (together with ∇ D = ∇), we obtain: Accordingly, the difference δH 1 between H 1 and its fixed point value B 1 decays faster when γ 1 is larger, which is what we found in the numerical simulations.
The stationary probability density (73) is shown in Figures 3 and 4 for two different diffusion constants.For a relatively small diffusion constant, the pre-factor √ H 1 has relatively little impact, and the distribution resembles a normal distribution; see Figure 3.In contrast, if the diffusion constant is sufficiently larger, the distorting impact of the pre-factor √ H 1 becomes dominant, and the shape of the distribution qualitatively differs from a normal distribution; see Figure 4. Irrespective of D, the distribution of the pseudo-invariant H 1 corresponding to the squared angular moment is unimodal and exhibits its maximum at a value H 1 > 0. This peak at a value different from the "ground state value" H 1 = 0 indicates that the system is an active system rather than a purely dissipative one [62].As in Figure 1, but for a larger parameter γ 1 : γ 1 = 3/100.Other parameters as in Figure 1: H  2: Other parameters as in Figure 3:

Invariants of Nambu Dynamics as Pseudo-Invariants
We proposed a stochastic framework to discuss active Nambu systems that allows one to address thermodynamic state variables, on the one hand, and generalized entropic measures, on the other hand.The approach is based on classical Nambu mechanics.Consequently, for an n-dimensional system, there are n − 1 functions H 1 , . . ., H n−1 that act as invariants in the classical Nambu system.However, in the stochastic, active variant of that system, the functions H 1 , . . ., H n−1 vary over time.There are two reasons.
in line with the so-called canonical-dissipative modeling effort for active systems, the systems feature a pumping mechanism that makes sure that at least one of the functions H 1 , . . ., H n−1 approaches a fixed point value.That is, while in classical Nambu mechanics, the function values H 1 , . . ., H n−1 are determined by the initial conditions (i.e., initial values of the state vector (r 1 , . . ., r n )); for active Nambu systems, at least one function H j converges to a fixed point value B j and in doing so assumes in the long time limit a value that does not depend on the initial conditions.Figures 1 and 2 illustrate this transient dynamics.After the transient period, we have |H j − B j | < C, where C > 0 can be made arbitrarily small.In this case, variation in time of H j may be considered as being negligibly small, and we may refer to this domain as the stationary domain.We may say that in the stationary domain, the system has approached an attractor defined by the original Nambu system that does not feature the pumping mechanism.In this context, note that the phrase "pumping mechanisms" actually refers to a pumping and damping mechanism.For example, if the initial conditions are such that H j > B j , then H j decays towards B j .
Second, taking fluctuating forces into account associated with the entropic form S, due to the impact of these fluctuating forces, the functions H 1 , . . ., H n−1 vary over time in an erratic fashion.This holds both in the transient and stationary domain.In the stationary domain, the fluctuations can be illustrated by means of the stationary probability density P j (H j ).We demonstrated the erratic variations of the functions H 1 , . . ., H n−1 for a stochastic, active spinning top model involving the functions H 1 and H 2 , where H 1 corresponds to half of the squared amount of the angular momentum.The probability density P 1 (H 1 ) is shown in Figures 3 and 4 for two systems that differ with respect to the "strength" of the impact of the fluctuating force (as quantified by the diffusion coefficient).
In view of these consideration, it might be useful to refer to the functions H 1 , . . ., H n−1 occurring in stochastic, active Nambu systems described by Equation (9) as pseudo-invariants rather than invariants.The functions H 1 , . . ., H n−1 of Equation ( 9) are pseudo-invariants in the sense that in the limit t → ∞ and in the "zero temperature limit" θ → 0, the variants over time of these functions become negligibly small.

Active, Stochastic Systems and Generalized, Non-Extensive Entropic Measures
Let us motivate the utilization of generalized entropy measures of the form (3) and in particular the Tsallis entropy (6) for active systems.To this end, we would like to focus on active systems that describe human and animal behavior and in particular the motor control system of humans and animals.
The Boltzmann-Gibbs-Shannon entropy is an extensive measure.That is, when applied to many body systems, it does not account for long-range interactions between individual subsystems of many body systems.While for many applications (for example in the statistical mechanics description of ideal gases), such long-range interaction can be neglected, human and animal behavior may results from a plenitude of interacting components that feature long-range interactions.For example, on the neuronal level, neurons in the motor cortex may need to synchronize their activity in order to produce a descending signal that controls limb movement.The axonal connections between these neurons may be considered as representing long-range connections rather than local connections.The Tsallis entropy (6) has been proposed with this goal in mind to account for non-extensivity effect as a result of long-range interactions between components [39].
Second, the control of human and animal motor behavior involves several subsystems, such as the neural systems, the muscular system and the human and animal body as such (i.e., the system of the limbs).These subsystem feature again subsystems.For example, within the neural system, we may distinguish between processes at the synapses, the membrane and in the axons.The point that we would like to make here is that the human and animal motor control system from a mechanistic (microscopic, bottom-up) perspective involves various components that are likely to evolve on different characteristic time scales.To model explicitly all of these components is a challenging task.Alternatively, one may think of using stochastic models that are able to capture processes that take place at different levels of consideration at the same time.Scale independent models as described for example by the Levy distribution are models that exhibit this feature.This kind of model is known to exhibit power law distributions rather than normal distributions or Boltzmann distributions.As we have shown in Section 3, when taking general entropic forms into account, the shape of the stationary distributions is determined by the entropic form.In particular, using a non-extensive thermostatistic approach based on the Tsallis entropy (6), for non-extensivity parameters q < 1, the power law distribution can be modeled; see Section 3.2.In other words, using a generalized thermostatistic approach that leads to a stationary power law distribution allows us to capture in a crude way that active biological systems and in particular human and animal motor control systems involve various subcomponents that are likely to act on different characteristic time scales.
Let us reiterate that in view of the considerations made above, the proposed stochastic Nambu mechanics models are promising models to capture on appropriately-defined descriptive, macroscopic levels characteristic properties of biological systems (e.g., that their observables satisfy power law distributions).In fact, for the migration (cell motion) of certain bacteria, experimental evidence has been found that the bacteria velocity scores are distributed in a certain range according to a power law [63].It has further been suggested that the absolute velocity scores satisfy a Tsallis distribution similar to the power law distribution (29).In short, while a stochastic model of the form (9) may not be able to provide a mechanistic account of the many interacting components of a biological system under consideration, the model might be used as a phenomenological model to address certain macroscopic aspects of the system.If so, the benefit of the approach outlined above is that thermodynamic variables can be introduced in a well-defined manner such that the performance of the system under consideration can be studied from a thermodynamical perspective.
Having argued that the case q < 1 in non-extensive thermostatistics can help us to build more realistic models of active systems, one can also argue that the alternative case q > 1 has relevance for our understanding of active systems [64].Active systems as frequently been found in biology are in general limited in their performance measures.For example, humans and animals can run just so fast, and translational and rotational limb movement can be performed only up to certain maximal values of speed and angular velocity.Cutoff distributions can take these limitations into account.Modeling human and animal behavior by means of cutoff distributions implies that we acknowledge that human behavior is subjected to limits.Accordingly, the observable related to the behavior of interest will not assume values beyond certain limits.The probability to find observable values beyond these limits equals zero.
The cases q < 1 and q > 1 are conflicting alternatives.When modeling a certain phenomenon of interest, one may give priority to either one of the two aforementioned aspects motivating the utilization of non-extensive thermostatistics with q < 1 or q > 1. Future efforts may be directed to combine both aspects in a single model.

Conclusions
Stochastic Nambu mechanics models provide a relatively novel theoretical framework to address systems in the life sciences from a physics perspective.This perspective includes both dynamic and thermodynamic aspects.Interestingly, the approach does not aim to establish a new physics of the animate world that is separated from the physics concerned with the inanimate world.Rather, as has been demonstrated above, the stochastic Nambu mechanics perspective is sufficiently broad to address both active systems as typically observed in the life sciences and equilibrium or close to equilibrium systems as frequently studied in classical physics.From a mathematical point of view, the stochastic models under consideration can be described by nonlinear partial differential equations for probability densities.Although in general, it is difficult to determine the properties of equations of this kind, the underlying thermodynamic aspects of the theory turns out to provide tools to analyze solutions.For example, we were able to determine certain properties of transient and stationary probability densities of the stochastic Nambu mechanics models.Importantly, we illustrated that the type of the nonlinearity of a partial differential equation under consideration corresponds to a certain type of an entropy measure.It is this link between the type of the entropy measure and the type of the nonlinearity that endows the proposed theoretical framework with powerful tools that can be used in applications to experimental data and theory-building.
Stationary probability density of the stochastic active spinning top model.Analytical (solid line) and numerical results (symbols) are shown.The diffusion constant D was relatively small.The full stochastic model defined by Equation (66) was solved numerically using a stochastic Euler forward scheme with a single time step of 0.005 time units in the time interval [0, 1000].From L(t) thus obtained, the invariant H 1 (t) was calculated.The numerical results show the probability density estimated from H 1 (t) in [500, 1000] (neglecting the transient period) using kernel density estimation with positive support.The analytical results were drawn from Equation (73).The effective integration factor Z Z was determined numerically.Model parameters: γ 1 = 0.3, D = 0.1.Other parameters as in Figure