# Controllability in Hybrid Kinetic Equations Modeling Nonequilibrium Multicellular Systems

Authors: Carlo Bianca
Journal: The Scientific World Journal
 (2013)

---

## Abstract

This paper is concerned with the derivation of hybrid kinetic partial integrodifferential equations that can be proposed for the mathematical modeling of multicellular systems subjected to external force fields and characterized by nonconservative interactions. In order to prevent an uncontrolled time evolution of the moments of the solution, a control operator is introduced which is based on the Gaussian thermostat. Specifically, the analysis shows that the moments are solution of a Riccati-type differential equation.


---

## Body

## 1. Introduction

The derivation of mathematical models for the evolution of multicellular systems is a challenge of the last and also this century. Several mathematical frameworks have been proposed with the aim of capturing the main properties of these systems, which have also the feature to be complex [1]. Indeed different emerging behaviors appear as result of the interactions: competition, pattern formation, and organization (see the book [2] and the review paper [3]).The mathematical frameworks present in the pertinent literature principally pertain to multicellular (biological) systems that are not subjected to external forces (systems at equilibrium); the interested reader is referred to the recent review [4] (and the references cited therein) for what concern the hyperbolic and kinetic models for living systems at equilibrium.However, dealing with multicellular systems whose time evolution depends on the presence of external fields requires a different treatment. An external force field moves the system away from equilibrium and the blow-up of some observable (such as the kinetic energy) occurs.Recently a new mathematical framework has been proposed and analyzed: the thermostated kinetic theory for active particles [5–9]. This framework, which refers to general complex systems (including vehicular traffic, and crowd and swarm dynamics; see, among others, [10–13]), is based on kinetic theory for active particles and deterministic thermostat. Kinetic equations coupled with the Gaussian isokinetic thermostat are already present in the literature: Kac equation [14–16] and Boltzmann equation with the one-dimensional BGK operator [17] (see also the recent review paper [18]). The solution of these equations is a distribution (probability) function defined in the microscopic state of the particles, which can include mechanical variables (space and velocity) and/or biological variables.Roughly speaking, the deterministic thermostat is a mathematical framework that has been proposed in nonequilibrium molecular dynamics in order to achieve equilibration to nonequilibrium situation where there is a flux of energy through the system induced by external fields. The reader interested in further details is referred to papers [19, 20] and the review paper [21]. Accordingly, a control term (a damping term) is inserted into the equations of motion in order to ensure the conservation of the kinetic energy (or a different moment); see [22, 23].The thermostated kinetic for active particles framework proposed in [7] is concerned with nonequilibrium complex systems characterized by conservative interactions only (namely, interactions that only modify the microscopic state of the particles). Therefore, the mathematical modeling of nonequilibrium multicellular systems, where proliferative and destructive events may occur, cannot directly performed within this framework.This paper is concerned with a further generalization of the thermostated kinetic framework [7] that takes into account the modeling of multicellular systems with nonconservative interactions. Moreover, the velocity variable is discretized, generating an hybrid framework with discrete and continuous variables. Therefore, the new thermostated kinetic framework consists of nonlinear hybrid partial integrodifferential equations with quadratic nonlinearity.It is worth stressing that, to the best of our knowledge, this is the first time that nonconservative interactions are taken into account in the thermostated kinetic for active particles framework.Specifically, in the present paper the time evolution of the moments is analyzed. The control term, which is based on the mathematical thermostats, allows the derivation of nonlinear ordinary differential equations (Riccati-type equations) fulfilled by the moments.The framework proposed here is certainly worth of future research activity concerning both its qualitative analysis and the application to modelling complex biological systems; see [7] and the references cited therein. Mathematical control in integrodifferential equations is already presented in the pertinent literatur; see, among others, [24].The contents of the present paper are divided into three sections which follow this introduction. Section2 highlights the essential hybrid mathematical settings. Section 3 deals with the analysis of the moments evolution. Finally, Section 4 is meant to future research perspectives.


## 2. The Underlying Hybrid Kinetic Framework

This section is meant to a concise description of an hybrid kinetic for active particle framework that constitutes a new paradigm for the modelling of complex multicellular systems characterized by nonconservative particle interactions. Specifically, these systems are composed by cells (active particles), with different genotypes and phenotypes, which are able to perform a task; for example, the immune system is constituted by cells that protect the human body from pathogens, and the connective tissue is made up of fibroblast cells that synthesize the extracellular matrix and collagen.The mathematical structure consists of autonomous partial integrodifferential equations system with quadratic nonlinearity constituted by three different operators: the operator which models conservative interactions, the operator that takes care of nonconservative interactions, and the operator that ensures the control of the evolution of the moments. Each equation is a kinetic-type equation whose solution is the distribution functionf(t,x,v,u) of cells that at time t∈[0,∞[, possess the triplet of microscopic variables (x,v,u), where x∈Dx⊂ℝ3 is the space variable, v∈Dv⊂ℝ3 is the velocity variable, and u∈Du⊂ℝ is the task-variable (activity). The elementary product f(t,x,v,u)dxdvdu is the number of cells whose state, at time t, is in the elementary volume of the space of microscopic states dΩ=[x,x+dx]×[v,v+dv]×[u,u+du]. Let Ω=Dx×Dv×Du be the domain of the all possible microscopic states and
(1)dμΩ=dxdvdu
the Lebesgue measure on Ω. If f∈L1[Ω,dμΩ], then the following expectation 𝔼0[f] of f:
(2)𝔼0[f](t)=∫Ω‍f(t,x,v,u)dμΩ
represents the global number, at time t, of the cells.Higher order moments off can be defined under suitable integrability assumption on f. Let ωm,n,p=|x|m|v|n|u|p, where |·| denotes the Euclidean distance of · and m,n,p∈ℕ. If, for instance, ωm,n,pf∈L1[Ω,dμΩ], the (m,n,p)th-order moment of f is defined as follows:
(3)𝔼m,n,p[f](t)=∫Ω‍ωm,n,pf(t,x,v,u)dμΩ.
It is worth stressing that if f represents the joint distribution function, marginal distribution functions of f refer either to the distribution over the mechanical state or to distribution over the microscopic activity. These marginal distribution functions define the local quantities. For instance, the local quadratic activity moment at time t in x is computed as follows:
(4)𝔸2[f](t,x)=∫Dv×Du‍u2f(t,x,v,u)dvdu
and represent the energy expressed by the activity variable.The task dynamics performed by cells is modified, at the timet, by binary interactions, which occur at the microscopic level and refer to the mutual actions between the cell with microscopic state (x,v,u) (test cell) and the cell with microscopic state (x*,v*,u*) (field cell), or the test cell and the cell with microscopic state (x*,v*,u*) (candidate cell). The candidate cell is the cell that reaches the microscopic state of the test cell after the interaction with the field cell.


### 2.1. The Hybrid Kinetic Setting at Equilibrium

Looking at the motion of cells, it is nonrestrictive to assume that velocity variables can attain discrete valuesv∈{v1,v2,…,vn} and the activity variable is continuous u∈Du. Therefore, the cell microscopic state includes both discrete and continuous variables. Moreover, we assume space homogeneity. Accordingly, the distribution function f of the system reads
(5)f(t,v,u)=∑i=1n‍f(t,vi,u)δ(v-vi)=∑i=1n‍fi(t,u)δ(v-vi).
The mathematical framework thus consists of a set of evolution equations for fi(t,u), which represents the distribution function of the test cell. Let fi(t,u*) be the distribution function of the candidate cell and fi(t,u*) the distribution function of the field cell.The time evolution of the distribution function of test cellfi(t,u) is modified when the candidate cell, with task-variable u*, after the interaction with the field cell, with task-variable u*, can acquire the task-variable u. The interactions in multicellular systems are conservative (modify the magnitude activity variable of the cells) and nonconservative (that modify the number of cells as a consequence of proliferation and destruction events occurring for natural birth/death and competition). Since fi(t,u)du denotes the number of cells, at time t, such that u∈[u,u+du], then
(6)fi(t,u*)du*fj(t,u*)du*
is a nonlinear product that refers to the interaction, at time t, between the number of candidate cells, with velocity vi such that u*∈[u*,u*+du*] and the number of field cells, with velocity vj, such that u*∈[u*,u*+du*]; the possibility of this interaction can be measured by introducing the nonnegative function ιij(u*,u*) which represents the interaction rate; the probability that after this interaction the candidate cell undergoes a change in its microscopic state (that of test particle) is measured by the nonnegative function 𝒜ij(u*,u*,u), which is a probability density with respect to u and then
(7)∫Du‍𝒜ij(u*,u*,u)du=1,∀u*,u*∈Du.
Bearing all above in mind, the (infinitesimal) result of the interaction reads
(8)𝒜ij(u*,u*,u)ιij(u*,u*)fi(t,u*)du*fj(t,u*)du*,
and summing up with respect to all the candidate and field cells, we obtain the following operator which models the gain of test cell:
(9)𝒢ij[fi,fj](t,u)=∫Du×Du‍ιij(u*,u*)𝒜ij(u*,u*,u)×fi(t,u*)fj(t,u*)du*du*.
Similarly, the loss of test cells is modeled by the following operator:
(10)ℒij[fi,fj](t,u)=fi(t,u)∫Du‍ιij(u,u*)fj(t,u*)du*.
Finally, if πij(u,u*) is the net birth/death rate of the test cell due to the interaction with the field cell, then the operator which models nonconservative interactions reads
(11)𝒩ij[fi,fj](t,u)=fi(t,u)∫Du‍ιij(u,u*)×πij(u,u*)fj(t,u*)du*.
The evolution equation of the distribution function fi over the microscopic state can be derived by a balance equation of the inlet and outlet flows in the elementary volume [u,u+du] of the space of the microscopic states. The hybrid kinetic for active particles framework thus reads:
(12)∂tfi(t,u)=∑j=1n‍(𝒢ij[fi,fj](t,u)-ℒij[fi,fj](t,u)+𝒩ij[fi,fj](t,u)).Definition 1.
Letιij(u1,u2):Du×Du→ℝ+, for i,j∈{1,2,…,n}, be the interaction rate between the u1-cell distributed according to fi(t,u1) and the u2-cell distributed according to f2(t,u2). Let 𝒜ij(u1,u2,u):Du×Du×Du→ℝ+ be the probability density satisfying the property (7). A function fi=fi(t,u):(0,∞)×Du→ℝ+ is said to be the solution of (12) if(i)
fi(t,u)∈C((0,∞),L1(Du));(ii)
fi is differentiable with respect to the variable t;(iii)
ιij(u1,u2)𝒜ij(u1,u2,u)fi(t,u1)fj(t,u2) is an integrable function with respect to the elementary measure du1du2;(iv)
ιij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(v)
ιij(u1,u2)πij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(vi)
fi satisfies (12) for all (t,u)∈(0,∞)×Du.Settingf=(f1(t,u),f2(t,u),…,fn(t,u))∈ℝn, the (p,q)-order moment of the distribution function f(t,v,u), for p,q∈ℕ, is written as follows:
(13)𝔼p,q[f](t)=∑i=1n‍vip∫Du‍uqfi(t,u)du.
In particular, the zero-order 𝔼0,0 (density or mass), first-order 𝔼1,1 (mean activation or linear momentum), and second-order 𝔼2,2 (activation energy or kinetic energy) moments fulfill an important role depending on the system under consideration.


### 2.2. The Controlled Hybrid Kinetic Setting at Nonequilibrium

The mathematical framework (12) is concerned with multicellular systems at equilibrium. Nonequilibrium conditions occur when the system is subjected to external fields Fi(u):Du→ℝ+ at macroscopic scale. In this case, the kinetic framework reads
(14)∂tfi(t,u)+∂u(Fi(u)fi(t,u))=∑j=1n‍(𝒢ij[fi,fj](t,u)-ℒij[fi,fj](t,u)+𝒩ij[fi,fj](t,u)).
The external field does work on the system thereby moving it away from equilibrium. Therefore, it follows uncontrolled increases of the activation energy (kinetic energy). The action of the Gaussian isokinetic thermostat is modeled, according to the Gauss’ principle of least constrain [20, 25], by the following operator:
(15)𝒯i[Fi,f](t,u)=uFi(u)fi(t,u)∑i=1n‍vi∫Du‍ufi(t,u)du
which is a damping operator (thermostat operator) that is adjusted so as to control the activation energy. The introduction of the thermostat operator modifies the mathematical framework as follows:
(16)∂tfi(t,u)+∂u(Fi(u)fi(t,u)-𝒯i[Fi,f](t,u))=∑j=1n‍(𝒞ij[fi,fj](t,u)+𝒩ij[fi,fj](t,u)),
where 𝒞ij[fi,fj](t,u)=𝒢ij[fi,fj](t,u)-ℒij[fi,fj](t,u) is the operator for the conservative interactions. In what follows, we refer to framework (16) as the controlled kinetic framework with conservative and nonconservative interactions.Definition 2.
LetFi=Fi(u), u∈Du, be an external force field differentiable with respect to u; ιij(u1,u2):Du×Du→ℝ+, for i,j∈{1,2,…,n}, interaction rate between the u1-cell distributed according to fi(t,u1) and the u2-cell distributed according to f2(t,u2); consider 𝒜ij(u1,u2,u):Du×Du×Du→ℝ+ to be the probability density satisfying the property (7). A function fi=fi(t,u):(0,∞)×Du→ℝ+ is said to be the solution of the model (16) if(i)
fi(t,u)∈C((0,∞),L1(Du));(ii)
fi is differentiable with respect to the variables t and u;(iii)
ufi is an integrable function with respect to the elementary measure du;(iv)
ιij(u1,u2)𝒜ij(u1,u2,u)fi(t,u1)fj(t,u2) is an integrable function with respect to the elementary measure du1du2;(v)
ιij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(vi)
ιij(u1,u2)πij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(vii)
𝒯i[Fi,f] is differentiable with respect to the variable u;(viii)
fi satisfies (16) for all (t,u)∈(0,∞)×Du.Remark 3.
The theorem of existence and uniqueness of the solution for the controlled kinetic framework (16) has been obtained in [7] when the nonconservative operator 𝒩ij is equal to zero (conservative interactions only). The proof of the theorem can be adapted in order to obtain existence and uniqueness of the solution also for the nonconservative interactions case. Nevertheless, global existence may not occur. This is a work in progress and results will be reported in due course.
The depicted hybrid controlled kinetic framework (16) is quite general and can be exploited to originate specific models for multicellular systems by acting on the specific forms of the grid velocity, interaction rate ιij, the probability density 𝒜ij, the net rate of birth/death πij, and the external force Fi.


## 2.1. The Hybrid Kinetic Setting at Equilibrium

Looking at the motion of cells, it is nonrestrictive to assume that velocity variables can attain discrete valuesv∈{v1,v2,…,vn} and the activity variable is continuous u∈Du. Therefore, the cell microscopic state includes both discrete and continuous variables. Moreover, we assume space homogeneity. Accordingly, the distribution function f of the system reads
(5)f(t,v,u)=∑i=1n‍f(t,vi,u)δ(v-vi)=∑i=1n‍fi(t,u)δ(v-vi).
The mathematical framework thus consists of a set of evolution equations for fi(t,u), which represents the distribution function of the test cell. Let fi(t,u*) be the distribution function of the candidate cell and fi(t,u*) the distribution function of the field cell.The time evolution of the distribution function of test cellfi(t,u) is modified when the candidate cell, with task-variable u*, after the interaction with the field cell, with task-variable u*, can acquire the task-variable u. The interactions in multicellular systems are conservative (modify the magnitude activity variable of the cells) and nonconservative (that modify the number of cells as a consequence of proliferation and destruction events occurring for natural birth/death and competition). Since fi(t,u)du denotes the number of cells, at time t, such that u∈[u,u+du], then
(6)fi(t,u*)du*fj(t,u*)du*
is a nonlinear product that refers to the interaction, at time t, between the number of candidate cells, with velocity vi such that u*∈[u*,u*+du*] and the number of field cells, with velocity vj, such that u*∈[u*,u*+du*]; the possibility of this interaction can be measured by introducing the nonnegative function ιij(u*,u*) which represents the interaction rate; the probability that after this interaction the candidate cell undergoes a change in its microscopic state (that of test particle) is measured by the nonnegative function 𝒜ij(u*,u*,u), which is a probability density with respect to u and then
(7)∫Du‍𝒜ij(u*,u*,u)du=1,∀u*,u*∈Du.
Bearing all above in mind, the (infinitesimal) result of the interaction reads
(8)𝒜ij(u*,u*,u)ιij(u*,u*)fi(t,u*)du*fj(t,u*)du*,
and summing up with respect to all the candidate and field cells, we obtain the following operator which models the gain of test cell:
(9)𝒢ij[fi,fj](t,u)=∫Du×Du‍ιij(u*,u*)𝒜ij(u*,u*,u)×fi(t,u*)fj(t,u*)du*du*.
Similarly, the loss of test cells is modeled by the following operator:
(10)ℒij[fi,fj](t,u)=fi(t,u)∫Du‍ιij(u,u*)fj(t,u*)du*.
Finally, if πij(u,u*) is the net birth/death rate of the test cell due to the interaction with the field cell, then the operator which models nonconservative interactions reads
(11)𝒩ij[fi,fj](t,u)=fi(t,u)∫Du‍ιij(u,u*)×πij(u,u*)fj(t,u*)du*.
The evolution equation of the distribution function fi over the microscopic state can be derived by a balance equation of the inlet and outlet flows in the elementary volume [u,u+du] of the space of the microscopic states. The hybrid kinetic for active particles framework thus reads:
(12)∂tfi(t,u)=∑j=1n‍(𝒢ij[fi,fj](t,u)-ℒij[fi,fj](t,u)+𝒩ij[fi,fj](t,u)).Definition 1.
Letιij(u1,u2):Du×Du→ℝ+, for i,j∈{1,2,…,n}, be the interaction rate between the u1-cell distributed according to fi(t,u1) and the u2-cell distributed according to f2(t,u2). Let 𝒜ij(u1,u2,u):Du×Du×Du→ℝ+ be the probability density satisfying the property (7). A function fi=fi(t,u):(0,∞)×Du→ℝ+ is said to be the solution of (12) if(i)
fi(t,u)∈C((0,∞),L1(Du));(ii)
fi is differentiable with respect to the variable t;(iii)
ιij(u1,u2)𝒜ij(u1,u2,u)fi(t,u1)fj(t,u2) is an integrable function with respect to the elementary measure du1du2;(iv)
ιij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(v)
ιij(u1,u2)πij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(vi)
fi satisfies (12) for all (t,u)∈(0,∞)×Du.Settingf=(f1(t,u),f2(t,u),…,fn(t,u))∈ℝn, the (p,q)-order moment of the distribution function f(t,v,u), for p,q∈ℕ, is written as follows:
(13)𝔼p,q[f](t)=∑i=1n‍vip∫Du‍uqfi(t,u)du.
In particular, the zero-order 𝔼0,0 (density or mass), first-order 𝔼1,1 (mean activation or linear momentum), and second-order 𝔼2,2 (activation energy or kinetic energy) moments fulfill an important role depending on the system under consideration.


## 2.2. The Controlled Hybrid Kinetic Setting at Nonequilibrium

The mathematical framework (12) is concerned with multicellular systems at equilibrium. Nonequilibrium conditions occur when the system is subjected to external fields Fi(u):Du→ℝ+ at macroscopic scale. In this case, the kinetic framework reads
(14)∂tfi(t,u)+∂u(Fi(u)fi(t,u))=∑j=1n‍(𝒢ij[fi,fj](t,u)-ℒij[fi,fj](t,u)+𝒩ij[fi,fj](t,u)).
The external field does work on the system thereby moving it away from equilibrium. Therefore, it follows uncontrolled increases of the activation energy (kinetic energy). The action of the Gaussian isokinetic thermostat is modeled, according to the Gauss’ principle of least constrain [20, 25], by the following operator:
(15)𝒯i[Fi,f](t,u)=uFi(u)fi(t,u)∑i=1n‍vi∫Du‍ufi(t,u)du
which is a damping operator (thermostat operator) that is adjusted so as to control the activation energy. The introduction of the thermostat operator modifies the mathematical framework as follows:
(16)∂tfi(t,u)+∂u(Fi(u)fi(t,u)-𝒯i[Fi,f](t,u))=∑j=1n‍(𝒞ij[fi,fj](t,u)+𝒩ij[fi,fj](t,u)),
where 𝒞ij[fi,fj](t,u)=𝒢ij[fi,fj](t,u)-ℒij[fi,fj](t,u) is the operator for the conservative interactions. In what follows, we refer to framework (16) as the controlled kinetic framework with conservative and nonconservative interactions.Definition 2.
LetFi=Fi(u), u∈Du, be an external force field differentiable with respect to u; ιij(u1,u2):Du×Du→ℝ+, for i,j∈{1,2,…,n}, interaction rate between the u1-cell distributed according to fi(t,u1) and the u2-cell distributed according to f2(t,u2); consider 𝒜ij(u1,u2,u):Du×Du×Du→ℝ+ to be the probability density satisfying the property (7). A function fi=fi(t,u):(0,∞)×Du→ℝ+ is said to be the solution of the model (16) if(i)
fi(t,u)∈C((0,∞),L1(Du));(ii)
fi is differentiable with respect to the variables t and u;(iii)
ufi is an integrable function with respect to the elementary measure du;(iv)
ιij(u1,u2)𝒜ij(u1,u2,u)fi(t,u1)fj(t,u2) is an integrable function with respect to the elementary measure du1du2;(v)
ιij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(vi)
ιij(u1,u2)πij(u1,u2)fj(t,u2) is an integrable function with respect to the elementary measure du2;(vii)
𝒯i[Fi,f] is differentiable with respect to the variable u;(viii)
fi satisfies (16) for all (t,u)∈(0,∞)×Du.Remark 3.
The theorem of existence and uniqueness of the solution for the controlled kinetic framework (16) has been obtained in [7] when the nonconservative operator 𝒩ij is equal to zero (conservative interactions only). The proof of the theorem can be adapted in order to obtain existence and uniqueness of the solution also for the nonconservative interactions case. Nevertheless, global existence may not occur. This is a work in progress and results will be reported in due course.
The depicted hybrid controlled kinetic framework (16) is quite general and can be exploited to originate specific models for multicellular systems by acting on the specific forms of the grid velocity, interaction rate ιij, the probability density 𝒜ij, the net rate of birth/death πij, and the external force Fi.


## 3. Differential Equations for the Moments

This section is concerned with the derivation of differential equations for the moments. Let𝔼1[f] be the following moment:
(17)𝔼1[f](t)≔𝔼1,1[f](t)=∑i=1n‍vi∫Du‍ufi(t,u)du.
Let μ(t) be the following function:
(18)μ(t)≔∑i=1n‍μi(t),
where
(19)μi(t)≔∫Du‍fi(t,u)du,(20)μ¯(t)≔∑i=1n‍viμi(t).The following result holds true.Theorem 4.
Letιij=ι, πij=π, and Fi=F be real constants. If there exists a nonnegative solution f of the controlled kinetic framework (16) such that fi(t,u)=0 as u∈∂Du, then the 1th-order moment 𝔼1[f](t) is solution of the following Riccati nonlinear ordinary differential equation:
(21)ddt𝔼1[f](t)=F(μ¯(t)-(𝔼1[f](t))2)-ι(1-π)μ(t)𝔼1[f](t).Proof.
The interaction operator𝒥[fi,fj] can be written as follows:
(22)𝒥ij[fi,fj](t,u)=𝒢ij[fi,fj](t,u)-ιμj(t)fi(t,u)+ιπμj(t)fi(t,u).
Multiplying both sides of 𝒥ij[fi,fj] by u and integrating over Du, we have
(23)∫Du‍u𝒥ij[fi,fj](t,u)du=-ιμj(t)(1-π)∫Du‍ufi(t,u)du.
Summing with respect to j, multiplying by vi, and summing with respect to i, we obtain
(24)∑i=1n‍vi∑j=1n‍∫Du‍u𝒥ij[fi,fj](t,u)du=-ιμ(t)(1-π)𝔼1[f](t).
Multiplying by u and vi the second term of the left hand side of (16), integrating with respect to the activity variable, performing integration by parts and summing with respect to i, we have
(25)∑i=1n‍vi∫Du‍u∂u((1-u𝔼1[f](t))fi(t,u))du=(𝔼1[f](t))2-μ¯(t)
and then the proof.According to Theorem4, the solution of the Riccati equation (21) can be obtained as follows. The Riccati equation reads
(26)ddt𝔼1[f](t)+F(𝔼1[f](t))2+ι(1-π)μ(t)𝔼1[f](t)-Fμ¯(t)=0.
Setting
(27)γ1(t)=F,γ2(t)=ι(1-π)μ(t),γ3(t)=-Fμ¯(t),
if 𝔼1[f](t)¯ is a solution of (26), the general integral can be written as
(28)𝔼1[f](t)=𝔼1[f](t)¯+1λ(t),
where λ(t) is solution of
(29)λ′-[γ2(t)+2γ1(t)𝔼1[f](t)¯]λ=γ1(t).
A nonnegative and constant solution of (26) is
(30)𝔼1[f](t)¯=ι2(1-π)2μ2(t)+4F2μ¯(t)-ι(1-π)μ(t)2F.
Therefore, the solution of (26) can be written as follows:
(31)𝔼1[f](t)=𝔼1[f](t)¯+e∫0t‍β(τ)dτ(𝔼1[f](0)-𝔼1[f](t)¯)-1+F∫0t‍eβ(τ)dτ,
where
(32)β(τ)=ι2(1-π)2μ2(τ)+4F2μ¯(τ).The next theorem gives the evolution equation for all moments wherep is an odd number.Theorem 5.
Letp∈ℕ, q be an odd number and t≥0. Then, the (p,q)th-order moment of the distribution function f satisfies the following Riccati nonlinear ordinary differential equation:
(33)ddt𝔼p,q[f](t)=pF𝔼p,q-1[f](t)(μ~(t)-𝔼p,q[f](t))+ιμ(t)(1-π)𝔼p,q[f](t),
where
(34)μ~(t)=∑i=1n‍vipμi(t).
Moreover, if   𝔼p,q[f](t) is initially bounded, it remains bounded for all t>0.Proof.
The proof follows by multiplying both sides of (14) by up and performing integration by parts on the control term.


## 4. Research Perspectives

The controlled kinetic framework proposed in this paper allows the derivation of specific models for multicellular systems characterized by nonconservative interactions. This framework belongs to the class of thermostated kinetic for active particles models.The mathematical framework (16) can be further generalized in order to include the role of mutations; see Nowak [26]. This is an important issue in the cancer modeling [27, 28].A future research perspective is the generalization of the mathematical framework (16) to open systems subjected to external actions at the microscopic scale, for example, the role that the outer environment has in the whole dynamics [29].Perspectives include also the introduction of stochastic terms that model jump processes in the activity or in the velocity variable; see paper [30] and the references section.Moreover, the proof of the existence of solutions to the stationary problem(35)∂u(Fi(u)fi(u)-𝒯i[Fi,f](u))=∑j=1n‍(𝒞ij[fi,fj](u)+𝒩ij[fi,fj](u))
is missing. This proof has been gained for the conservative interactions case; see [31].An important research perspective is the development of suitable asymptotic limits for deriving macroscopic equations for the evolution of moments; see [17] and the references section. The derivation of these equations is based on suitable assumptions on the operators in order to obtain convergence results. Optimization investigation can also be performed on the solutions of these macroscopic equations, such as regularization and embedding results [32].It is worth stressing that the analysis of the moments performed in this paper can be straightforwardly applied also for the thermostated framework proposed in [9].



---

*Source: 274719-2013-09-26.xml*