The Scientific World Journal

Volume 2013, Article ID 274719, 6 pages

http://dx.doi.org/10.1155/2013/274719

## Controllability in Hybrid Kinetic Equations Modeling Nonequilibrium Multicellular Systems

Dipartimento di Scienze Matematiche, Politecnico, Corso Duca degli Abruzzi 24, 10129 Torino, Italy

Received 8 August 2013; Accepted 2 September 2013

Academic Editors: K. Ammari, M. M. Cavalcanti, and S. Sivasundaram

Copyright © 2013 Carlo Bianca. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.

#### 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. Section 2 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 function of cells that at time , possess the triplet of microscopic variables , where is the space variable, is the velocity variable, and is the task-variable (activity). The elementary product is the number of cells whose state, at time , is in the elementary volume of the space of microscopic states . Let be the domain of the all possible microscopic states and the Lebesgue measure on . If , then the following expectation of : represents the global number, at time , of the cells.

Higher order moments of can be defined under suitable integrability assumption on . Let , where denotes the Euclidean distance of and . If, for instance, , the th-order moment of is defined as follows: It is worth stressing that if represents the joint distribution function, marginal distribution functions of 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 in is computed as follows: and represent the energy expressed by the activity variable.

The task dynamics performed by cells is modified, at the time , by *binary interactions*, which occur at the microscopic level and refer to the mutual actions between the cell with microscopic state (test cell) and the cell with microscopic state (field cell), or the test cell and the cell with microscopic state (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 values and the activity variable is continuous . Therefore, the cell microscopic state includes both discrete and continuous variables. Moreover, we assume space homogeneity. Accordingly, the distribution function of the system reads The mathematical framework thus consists of a set of evolution equations for , which represents the distribution function of the test cell. Let be the distribution function of the candidate cell and the distribution function of the field cell.

The time evolution of the distribution function of test cell is modified when the candidate cell, with task-variable , after the interaction with the field cell, with task-variable , can acquire the task-variable . 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 denotes the number of cells, at time , such that , then
is a nonlinear product that refers to the interaction, at time , between the number of candidate cells, with velocity such that and the number of field cells, with velocity , such that ; the possibility of this interaction can be measured by introducing the nonnegative function 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 , which is a probability density with respect to and then
Bearing all above in mind, the (infinitesimal) result of the interaction reads
and summing up with respect to all the candidate and field cells, we obtain the following operator which models the gain of test cell:
Similarly, the loss of test cells is modeled by the following operator:
Finally, if 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
The evolution equation of the distribution function over the microscopic state can be derived by a balance equation of the inlet and outlet flows in the elementary volume of the space of the microscopic states. The hybrid kinetic for active particles framework thus reads:

*Definition 1. *Let , for , be the interaction rate between the -cell distributed according to and the -cell distributed according to . Let be the probability density satisfying the property (7). A function is said to be the solution of (12) if(i);(ii) is differentiable with respect to the variable ;(iii) is an integrable function with respect to the elementary measure ;(iv) is an integrable function with respect to the elementary measure ;(v) is an integrable function with respect to the elementary measure ;(vi) satisfies (12) for all .

Setting , the -order moment of the distribution function , for , is written as follows: In particular, the zero-order (density or mass), first-order (mean activation or linear momentum), and second-order (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 at macroscopic scale. In this case, the kinetic framework reads 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: 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: where 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. *Let , , be an external force field differentiable with respect to ; , for , interaction rate between the -cell distributed according to and the -cell distributed according to ; consider to be the probability density satisfying the property (7). A function is said to be the solution of the model (16) if(i);(ii) is differentiable with respect to the variables and ;(iii) is an integrable function with respect to the elementary measure ;(iv) is an integrable function with respect to the elementary measure ;(v) is an integrable function with respect to the elementary measure ;(vi) is an integrable function with respect to the elementary measure ;(vii) is differentiable with respect to the variable ;(viii) satisfies (16) for all .

*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 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 , the probability density , the net rate of birth/death , and the external force .

#### 3. Differential Equations for the Moments

This section is concerned with the derivation of differential equations for the moments. Let be the following moment: Let be the following function: where

The following result holds true.

Theorem 4. *Let , , and be real constants. If there exists a nonnegative solution of the controlled kinetic framework (16) such that as , then the th-order moment is solution of the following Riccati nonlinear ordinary differential equation:
*

*Proof. *The interaction operator can be written as follows:
Multiplying both sides of by and integrating over , we have
Summing with respect to , multiplying by , and summing with respect to , we obtain
Multiplying by and 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 , we have
and then the proof.

According to Theorem 4, the solution of the Riccati equation (21) can be obtained as follows. The Riccati equation reads Setting if is a solution of (26), the general integral can be written as where is solution of A nonnegative and constant solution of (26) is Therefore, the solution of (26) can be written as follows: where

The next theorem gives the evolution equation for all moments where is an odd number.

Theorem 5. *Let , be an odd number and . Then, the th-order moment of the distribution function satisfies the following Riccati nonlinear ordinary differential equation:
**
where
**
Moreover, if??? is initially bounded, it remains bounded for all . *

*Proof. *The proof follows by multiplying both sides of (14) by 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 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].

#### Acknowledgments

The author acknowledges the support by the FIRB Project RBID08PP3J—Metodi matematici e relativi strumenti per la modellizzazione e la simulazione della formazione di tumori, competizione con il sistema immunitario, e conseguenti suggerimenti terapeutici.

#### References

- S. A. Levin, “Complex adaptive systems: exploring the known, the unknown and the unknowable,”
*Bulletin of the American Mathematical Society*, vol. 40, no. 1, pp. 3–19, 2003. View at Publisher · View at Google Scholar · View at Scopus - D. J. Krause and G. D. Ruxton,
*Living in Groups*, Oxford University Press, Oxford, UK, 2002. - D. L. Abel and J. T. Trevors, “Self-organization vs. self-ordering events in life-origin models,”
*Physics of Life Reviews*, vol. 3, no. 4, pp. 211–228, 2006. View at Publisher · View at Google Scholar · View at Scopus - R. Eftimie, “Hyperbolic and kinetic models for self-organized biological aggregations and movement: a brief review,”
*Journal of Mathematical Biology*, vol. 65, no. 1, pp. 35–75, 2011. View at Publisher · View at Google Scholar · View at Scopus - C. Bianca, “Kinetic theory for active particles modelling coupled to Gaussian thermostats,”
*Applied Mathematical Sciences*, vol. 6, no. 13-16, pp. 651–660, 2012. View at Google Scholar · View at Scopus - C. Bianca, “An existence and uniqueness theorem to the Cauchy problem for thermostatted-KTAP models,”
*International Journal of Mathematical Analysis*, vol. 6, no. 17-20, pp. 813–824, 2012. View at Google Scholar - C. Bianca, “Onset of nonlinearity in thermostatted active particles models for complex systems,”
*Nonlinear Analysis: Real World Applications*, vol. 13, no. 6, pp. 2593–2608, 2012. View at Publisher · View at Google Scholar · View at Scopus - C. Bianca, “Modeling complex systems by functional subsystems representation and thermostatted-KTAP methods,”
*Applied Mathematics & Information Sciences*, vol. 6, pp. 495–499, 2012. View at Google Scholar - C. Bianca, M. Ferrara, and L. Guerrini, “High-order moments conservation in thermostatted kinetic models,”
*Journal of Global Optimization*, 2013. View at Publisher · View at Google Scholar - C. Dogbe, “Nonlinear pedestrian-flow model: uniform well-posedness and global existence,”
*Applied Mathematics & Information Sciences*, vol. 7, no. 1, pp. 29–40, 2013. View at Publisher · View at Google Scholar - C. Dogbe, “On the modelling of crowd dynamics by generalized kinetic models,”
*Journal of Mathematical Analysis and Applications*, vol. 387, no. 2, pp. 512–532, 2012. View at Publisher · View at Google Scholar · View at Scopus - C. Bianca and C. Dogbe, “A mathematical model for crowd dynamics: multiscale analysis, fluctuations and random noise,”
*Nonlinear Studies*, vol. 20, pp. 281–305, 2013. View at Google Scholar - D. Helbing and P. Molnár, “Social force model for pedestrian dynamics,”
*Physical Review E*, vol. 51, no. 5, pp. 4282–4286, 1995. View at Publisher · View at Google Scholar · View at Scopus - V. Bagland, B. Wennberg, and Y. Wondmagegne, “Stationary states for the noncutoff Kac equation with a Gaussian thermostat,”
*Nonlinearity*, vol. 20, no. 3, pp. 583–604, 2007. View at Publisher · View at Google Scholar · View at Scopus - B. Wennberg and Y. Wondmagegne, “Stationary states for the Kac equation with a Gaussian thermostat,”
*Nonlinearity*, vol. 17, no. 2, pp. 633–648, 2004. View at Publisher · View at Google Scholar · View at Scopus - B. Wennberg and Y. Wondmagegne, “The Kac equation with a thermostatted force field,”
*Journal of Statistical Physics*, vol. 124, no. 2-4, pp. 859–880, 2006. View at Publisher · View at Google Scholar · View at Scopus - P. Degond and B. Wennberg, “Mass and energy balance laws derived from high-field limits of thermostated Boltzmann equations,”
*Communications in Mathematical Sciences*, vol. 5, no. 2, pp. 355–382, 2007. View at Publisher · View at Google Scholar · View at Scopus - C. Bianca, “Thermostatted kinetic equations as models for complex systems in physics and life sciences,”
*Physics of Life Reviews*, vol. 9, no. 4, pp. 359–399, 2012. View at Publisher · View at Google Scholar - D. J. Evans, W. G. Hoover, B. H. Failor, B. Moran, and A. J. C. Ladd, “Nonequilibrium molecular dynamics via Gauss's principle of least constraint,”
*Physical Review A*, vol. 28, no. 2, pp. 1016–1021, 1983. View at Publisher · View at Google Scholar · View at Scopus - D. J. Evans and G. P. Morriss,
*Statistical Mechanics of Nonequilibrium Fluids*, Academic Press, New York, NY, USA, 1990. - O. G. Jepps and L. Rondoni, “Deterministic thermostats, theories of nonequilibrium systems and parallels with the ergodic condition,”
*Journal of Physics A*, vol. 43, no. 13, Article ID 133001, 2010. View at Publisher · View at Google Scholar · View at Scopus - G. P. Morriss and C. P. Dettmann, “Thermostats: analysis and application,”
*Chaos*, vol. 8, no. 2, pp. 321–336, 1998. View at Google Scholar · View at Scopus - D. Ruelle, “Smooth dynamics and new theoretical ideas in nonequilibrium statistical mechanics,”
*Journal of Statistical Physics*, vol. 95, no. 1-2, pp. 393–468, 1999. View at Publisher · View at Google Scholar · View at Scopus - S. Sivasundaram and J. Uvah, “Controllability of impulsive hybrid integro-differential systems,”
*Nonlinear Analysis: Hybrid Systems*, vol. 2, no. 4, pp. 1003–1009, 2008. View at Publisher · View at Google Scholar · View at Scopus - K. F. Gauss, “On a New Fundamental Law of Mechanics,”
*Journal für die Reine und Angewandte Mathematik*, vol. 4, pp. 232–235, 1829. View at Google Scholar - M. Nowak,
*Evolutionary Dynamics: Exploring the Equations of Life*, Belknap Press, 2006. - C. Bianca and M. Pennisi, “The triplex vaccine effects in mammary carcinoma: a nonlinear model in tune with SimTriplex,”
*Nonlinear Analysis: Real World Applications*, vol. 13, no. 4, pp. 1913–1940, 2012. View at Publisher · View at Google Scholar · View at Scopus - M. Pennisi, F. Pappalardo, and S. Motta, “Agent based modeling of lung metastasis-immune system competition,”
*Lecture Notes in Computer Science*, vol. 5666, pp. 1–3, 2009. View at Publisher · View at Google Scholar · View at Scopus - R. M. May, “Uses and abuses of mathematics in biology,”
*Science*, vol. 303, no. 5659, pp. 790–793, 2004. View at Publisher · View at Google Scholar · View at Scopus - A. Bellouquid and C. Bianca, “Modelling aggregation-fragmentation phenomena from kinetic to macroscopic scales,”
*Mathematical and Computer Modelling*, vol. 52, no. 5-6, pp. 802–813, 2010. View at Publisher · View at Google Scholar · View at Scopus - C. Bianca, “Existence of stationary solutions in kinetic models with Gaussian thermostats,”
*Mathematical Methods in the Applied Sciences*, vol. 36, no. 13, pp. 1768–1775, 2013. View at Publisher · View at Google Scholar - M. A. Ragusa, “Commutators of fractional integral operators on Vanishing-Morrey spaces,”
*Journal of Global Optimization*, vol. 40, no. 1-3, pp. 361–368, 2008. View at Publisher · View at Google Scholar · View at Scopus