Next Article in Journal
Optimal Control Analysis of Cholera Dynamics in the Presence of Asymptotic Transmission
Next Article in Special Issue
Approximations of an Equilibrium Problem without Prior Knowledge of Lipschitz Constants in Hilbert Spaces with Applications
Previous Article in Journal
Qualitative Analyses of Integro-Fractional Differential Equations with Caputo Derivatives and Retardations via the Lyapunov–Razumikhin Method
Previous Article in Special Issue
Approximation Results for Equilibrium Problems Involving Strongly Pseudomonotone Bifunction in Real Hilbert Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Towards the Dependence on Parameters for the Solution of the Thermostatted Kinetic Framework

Dipartimento di Matematica e Fisica, Università degli Studi della Campania “L. Vanvitelli”, Viale Lincoln 5, I-81100 Caserta, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Axioms 2021, 10(2), 59; https://doi.org/10.3390/axioms10020059
Submission received: 13 February 2021 / Revised: 29 March 2021 / Accepted: 31 March 2021 / Published: 12 April 2021
(This article belongs to the Special Issue Numerical Analysis and Computational Mathematics)

Abstract

:
A complex system is a system involving particles whose pairwise interactions cannot be composed in the same way as in classical Mechanics, i.e., the result of interaction of each particle with all the remaining ones cannot be expressed as a sum of its interactions with each of them (we cannot even know the functional dependence of the total interaction on the single interactions). Moreover, in view of the wide range of its applications to biologic, social, and economic problems, the variables describing the state of the system (i.e., the states of all of its particles) are not always (only) the usual mechanical variables (position and velocity), but (also) many additional variables describing e.g., health, wealth, social condition, social rôle ⋯, and so on. Thus, in order to achieve a mathematical description of the problems of everyday’s life of any human society, either at a microscopic or at a macroscpoic scale, a new mathematical theory (or, more precisely, a scheme of mathematical models), called KTAP, has been devised, which provides an equation which is a generalized version of the Boltzmann equation, to describe in terms of probability distributions the evolution of a non-mechanical complex system. In connection with applications, the classical problems about existence, uniqueness, continuous dependence, and stability of its solutions turn out to be particularly relevant. As far as we are aware, however, the problem of continuous dependence and stability of solutions with respect to perturbations of the parameters expressing the interaction rates of particles and the transition probability densities (see Section The Basic Equations has not been tackled yet). Accordingly, the present paper aims to give some initial results concerning these two basic problems. In particular, Theorem 2 reveals to be stable with respect to small perturbations of parameters, and, as far as instability of solutions with respect to perturbations of parameters is concerned, Theorem 3 shows that solutions are unstable with respect to “large” perturbations of interaction rates; these hints are illustrated by numerical simulations that point out how much solutions corresponding to different values of parameters stay away from each other as t + .
MSC:
82B40; 37F05; 45M10; 35B30; 34A12

1. Introduction

The present paper deals with the system of equations governing the behavior of so-called complex systems (see Section 2 for details). Roughly speaking, a complex system is a set of a large number of individuals (particles) whose behavior is strongly influenced by their mutual interactions, in addition to external forces and possibly to a thermostat [1,2,3,4], so that the evolution of the system cannot be by no means deterministic, but must be described in terms of the probability distribution fuction on the set of possible values of a suitable variable describing the state of each individual. In this connection, it must be carefully noted that—though the notion of a complex system was originated in a purely mechanical framework and could be traced back to Boltzmann’s Kinetic Theory of Gases [5,6,7,8]—a complex system is not nowadays considered as simply consisting of material particles, whose state is completely described by the two variables position and velocity. The notion of complex system has been exported in several different contexts (biology [9,10], medicine [11,12,13,14,15], economy [16,17], psychology [18,19], social dynamics [20,21,22,23] ⋯), in which the state variables are non-mechanical and, in at least one case, vectorial [24].
Though the behavior of each particle of the system is of course deterministic, i.e., it is uniquely determined by its interactions with other particles; nevertheless, the number of particles and interactions is so large as to prevent us from following the evolution of the state of each particle. Accordingly, the model is based on the choice to describe the evolution of the system as a whole, turning the attention to the probability distribution on the states of particles; thus, the evolution equation takes the form (1) (see Section 2). As usual for nonlinear differential and integral equations, also in this case we have to tackle the classical problems about existence, uniqueness, continuous dependence, and stability with respect to initial values (and boundary values, when required). These problems have been tackled in [25,26,27,28,29,30].
In many cases of interest, the solutions to the problems about stability and continuous dependence of solutions depend on the coefficients of the equations, especially in the cases in which they are not constants but functions of the independent variables. In this last case, the question of whether two solutions, corresponding to the same assigned data but to two different systems of coefficients, are close when such are the coefficients spontaneously arises. This question seems to be of special relevance for Equation (1). As we shall see in more detail in Section 2, in Equation (1), denoting by u the state variable and by D u the state space, that is the set of all possible values of u, we find two kinds of coefficients:
1.
the coefficient η ( u * , u * ) , a function defined on D u 2 , expressing the interaction rate of the particles whose state is u * with the particles whose state is u * , i.e.,— roughly speaking—the number of their interactions per unit time;
2.
the coefficient A ( u * , u * , u ) , a function defined on D u 3 , expressing the transition probability, i.e., the probability (density) that any individual in the state u * , when interacting with a particle in the state u * , falls in the state u.
In any context, it is obvious that, to different prescriptions on the form of function η ( u * , u * ) or of function A ( u * , u * , u ) , there will correspond different probability distributions on D u (or, in a strictly statistical interpretation, different distributions of relative frequencies on D u over the system). However, what should we expect about the dependence of the difference of distributions on the difference between prescribed coefficients? Should an accordingly small difference between the corresponding distributions correspond to small perturbations to the interaction rate or to the transition distribution?
These questions are quite similar to those posed in all the classical problems associated with differential and integro-differential equations, but—as far as we are aware—have not been tackled yet for Equation (1). Nevertheless, in view of the large number of applications of complex systems (and of Equation (1), which describes their evolution) to so many basic problems of collective life of the whole mankind (for instance, let us mention the prediction of the evolution of epidemic diseases, or of the emergence of unsustainable economic inequalities), these questions are of special relevance in the framework of KTAP. In addition, the present paper is the first attempt to tackle them and to give some initial results about both the continuous dependence of solutions of Equation (1) on the coefficients and their instability. In connection with this last topic, the paper also offers some numerical simulations that show the separation between solutions corresponding to different values of parameters.
The contents of the paper are distributed as follows: in Section 2, we recall the structure of KTAP theory and Equation (1), and report the Cauchy problem associated with it, in the case in which the activity variable is assumed to be continuous (Section 2.1), as well as in the case in which it is assumed to be discrete (Section 2.2); Section 3 will be devoted to draw the notion of dependence of solutions on the parameters, and we state and prove a result concerning the continuous dependence of solutions on parameters, again in both the continuous case (Section 3.1) and the discrete case (Section 3.2), and, in Section 4, we give first results about instability of solutions (in both cases); in this connection, special attention should be paid to the numerical simulations presented in Section 4.3 for the discrete case that offers a clear perception of the fact that the solutions to Equation (4) depend continuously on the parameters, but stay apart from each other when the perturbation of the parameters is greater than a well-defined threshold value. Finally, in Section 5, we outline some research perspectives based on some general and meaningful conclusions that can be drawn from the results found in the previous sections.

2. The Basic Equations

2.1. The Continuous Activity Framework

Let D u R and F > 0 . According to what has been laid out in the Introduction, this paper is devoted to the analysis of properties of solutions f ( t , u ) : [ 0 , + [ × D u R + of the following nonlinear integro-differential equation, with quadratic nonlinearity:
t f ( t , u ) + F u ( 1 u E 1 [ f ] ( t ) ) f ( t , u ) = J [ f , f ] ( t , u ) ,
where the operator J [ f , f ] ( t , u ) is defined as follows:
J [ f , f ] ( t , u ) = G [ f , f ] ( t , u ) L [ f , f ] ( t , u ) = D u × D u η ( u * , u * ) A ( u * , u * , u ) f ( t , u * ) f ( t , u * ) d u * d u * + f ( t , u ) D u η ( u , u * ) f ( t , u * ) d u * ,
and
  • η ( u * , u * ) : D u × D u R + ;
  • A ( u * , u * , u ) : D u × D u × D u R + with the property:
    D u A ( u * , u * , u ) d u = 1 , u * , u * D u ;
  • E 1 [ f ] ( t ) = D u u f ( t , u ) d u .
The Cauchy problem associated with Equation (1) reads
t f ( t , u ) + F u ( 1 u E 1 [ f ] ( t ) ) f ( t , u ) = J [ f , f ] ( t , u ) ( t , u ) [ 0 , + [ × D u f ( 0 , u ) = f 0 ( u ) u D u .
Let
E 0 [ f ] ( t ) = D u f ( t , u ) d u
and
E 2 [ f ] ( t ) = D u u 2 f ( t , u ) d u .
Consider the function space K ( D u ) defined as
K ( D u ) : = f ( t , u ) [ 0 , + [ × D u R + : E 0 [ f ] ( t ) = E 2 [ f ] ( t ) = 1 .
The existence and uniqueness of solutions
f ( t , u ) C ( 0 , + ) × D u ; L 1 ( D u ) K ( D u )
of the Cauchy problem (3) are proved in [25], under the condition
f ( t , u ) = 0 , u D u .
The existence of solutions of the nonequilibrium stationary problem related to (1) is proved in [31]. A proof of the convergence of the solution of (3) to the nonequilibrium stationary solution as time goes to infinity is given in [32].
In many cases of interest, as for example in the description of the diffusion of epidemics η ( u * , u * ) can be supposed to be constant, i.e., there exists η > 0 such that η ( u * , u * ) = η , for all u * , u * D u .
Remark 1.
If C is a complex system, homogeneous with respect to the mechanical variables, i.e., space and velocity, (1) describes the evolution of the distribution function f ( t , u ) of C , and is called thermostatted kinetic framework [2].
The microscopic state is described by a scalar variable u, called activity, which attains its values in a real continuous subset D u . In this frame:
  • η ( u * , u * ) is the interaction rate between the particles that are in the state u * and the particles in the state u * ;
  • A ( u * , u * , u ) is the transition probability density i.e., the probability (density) that a particle in the state u * falls into the state u after interacting with a particle in the state u * ;
  • F > 0 is the value of the external force field acting on the system C ;
  • E 0 [ f ] ( t ) is the density, E 1 [ f ] ( t ) is the linear momentum and E 2 [ f ] ( t ) is the global energy;
  • G [ f , f ] ( t , u ) is the gain-term operator and L [ f , f ] ( t , u ) is the loss-term operator.
Equation (1) and the related problem (3) describe the evolution of a system C such that the global activation energy, E 2 [ f ] ( t ) , is kept constant by means of a thermostat [33].

2.2. The Discrete Activity Framework

Let I u = { u 1 , u 2 , , u n } be a discrete subset of R . The operator J i [ f ] ( t ) , for i { 1 , 2 , , n } is defined as:
J i [ f ] ( t ) = G i [ f ] ( t ) L i [ f ] ( t ) = h = 1 n k = 1 n η h k B h k i f h ( t ) f k ( t ) f i ( t ) k = 1 n η i k f k ( t ) ,
where η h k : I u × I u R + , for h , k { 1 , 2 , , n } , and the functions B h k i : I u × I u × I u R + (where i , h , k { 1 , 2 , , n } ) obey the condition
i = 1 n B h k i = 1 , h , k { 1 , 2 , , n } .
Let f ( t ) = f 1 ( t ) , f 2 ( t ) , , f n ( t ) , where, for any i { 1 , 2 , n } ,
f i ( t ) : = f ( t , u i ) : [ 0 , + [ × I u R +
is a solution of the nonlinear ordinary differential equation
d f i d t ( t ) = J i [ f ] ( t ) + F i ( t ) i = 1 n u i 2 J i [ f ] + F i E 2 [ f ] f i ( t ) ,
for F ( t ) = F 1 ( t ) , F 2 ( t ) , , F n ( t ) with F i ( t ) > 0 . The 2-nd order moment function E 2 [ f ] ( t ) of f takes now the form
E 2 [ f ] ( t ) = i = 1 n u i 2 f i ( t ) .
Consider the function space:
R f 2 = R f 2 R + ; E 2 = f C [ 0 , + ] ; R + n : E 2 [ f ] = E 2
where E 2 R + . The existence and uniqueness of solutions to the Cauchy problem associated with Equation (4), with initial data f 0 such that i = 1 n u i 2 f 0 = 1 , has been proved in [34] under the following assumption:
H1 
There exist η , F > 0 , such that F i ( t ) F , for t > 0 , and η h k η , for h , k { 1 , 2 , , n } .
A nonequilibrium stationary solution of Equation (4), for i { 1 , 2 , , n } , is a function f i satisfying the equation
J i [ f ] + F i i = 1 n u i 2 J i [ f ] + F i E 2 f i = 0 .
Let R ˜ f 2 denote the function space:
R ˜ f 2 R + ; E 2 = f R + n : E 2 [ f ] = E 2 .
The existence of nonequilibrium stationary solutions g ( u ) R ˜ f 2 has been proved in [28], under the assumption H1.
In particular, under the further assumptions:
H2 
i = 1 n u i B h k i = 0 , for all h , k { 1 , 2 , , n } ,
H3 
i = 1 n u i 2 B h k i = u h 2 , for all h , k { 1 , 2 , , n } ,
it has been proved in [28] that any nonequilibrium stationary solution is unique if the force field verifies the constraint
F > 2 η E 2 2 1 + 1 u 2 2 .
Proposition 1
([28]). If assumptions H1H3 are met, together with the assumption
H4 
E 0 [ f ] = E 2 [ f ] = 1
then
1.
The evolution equation of E 1 [ f ] ( t ) = i = 1 n u i f i ( t ) takes the form
E 1 [ f ] ( t ) + η + i = 1 n u i 2 f i E 1 [ f ] ( t ) i = 1 n u i F i = 0 ;
2.
as t + ,
E 1 [ f ] ( t ) K : = i = 1 n u i F i η + i = 1 n u i 2 F i ;
3.
Denoting by f 0 the initial data of the Cauchy problem related to (4), one has
E 1 [ f ] ( t ) K c exp η + i = 1 n u i 2 F i t ,
where c is a constant depending on the system.
Remark 2.
In [35], the existence of solutions of Equation (4) and of the related nonequilibrium stationary problem has been proved for more general values of the real discrete variable u i .
Remark 3.
Equation (4) has been proposed in [34] in order to model a complex system which partitioned in n subsystems, called functional subsystems. In particular:
  • the function f i ( t ) , for i { 1 , 2 , , n } , denotes the distribution function of the i-th functional subsystem;
  • the function F ( t ) = ( F 1 ( t ) , F 2 ( t ) , , F n ( t ) ) is the external force field acting on the whole system;
  • The term
    α : = i = 1 n u i 2 J i [ f ] + F i E 2
    represents the thermostat term, which allows for keeping constant the quantity E 2 [ f ] ( t ) :
  • the term η h k is the interaction rate related to the encounters between the functional subsystem h and the functional subsystem k, for h , k { 1 , 2 , , n } ;
  • the function B h k i denotes the transition probability density that the functional subsystem h falls into the i after interacting with the functional subsystem k, for i , h , k { 1 , 2 , , n } ;
  • the operator J i [ f ] ( t ) , for i { 1 , 2 , , n } , models the net flux to the i-th functional subsystem; G i [ f ] ( t ) denotes the gain term operator (incoming flux) and L i [ f ] ( t ) the loss term operator (outgoing flux).
Remark 4.
Let p N . Equation (4) can be further generalized as follows:
d f i d t ( t ) = J i [ f ] ( t ) + F i ( t ) i = 1 n u i p J i [ f ] + F i E p [ f ] f i ( t ) .
This framework allows for keeping thep-th order moment
E p [ f ] ( t ) = i = 1 n u i p f i ( t ) .
of the distribution f constant.
Remark 5.
The convergence of any solution of (4) to a corresponding nonequilibrium stationary state (solution to (5)), as time goes to infinity, has been proved in [27].

3. The Continuous Dependence on the Parameters

3.1. The Continuous Activity Framework

Let A ( u * , u * , u ) , A ˜ ( u * , u * , u ) , η and η ˜ two classes of parameters for Equation (1). Let J [ f , f ] ( t , u ) = G [ f , f ] ( t , u ) L [ f , f ] ( t , u ) be the operator related to the parameters A ( u * , u * , u ) and η , and J ˜ [ f , f ] ( t , u ) = G ˜ [ f , f ] ( t , u ) L ˜ [ f , f ] ( t , u ) the operator related to the parameters A ˜ ( u * , u * , u ) and η ˜ .
The related Cauchy problems, with the same initial data f 0 ( u ) , are defined as follows:
t f ( t , u ) + F u ( 1 u E 1 [ f ] ( t ) ) f ( t , u ) = J [ f , f ] ( t , u ) f ( 0 , u ) = f 0 ( u ) ,
t f ( t , u ) + F u ( 1 u E 1 [ f ] ( t ) ) f ( t , u ) = J ˜ [ f , f ] ( t , u ) f ( 0 , u ) = f 0 ( u ) .
By [25], there exist two functions f ( t , u ) C ( 0 , + ) × D u ; L 1 ( D u ) K ( D u ) and f ˜ ( t , u ) C ( 0 , + ) × D u ; L 1 ( D u ) K ( D u ) that are solutions to problem (9) and to problem (10), respectively.
The present paper aims to give a contribution in two directions:
(1)
the continuous dependence of the solutions of Equation (1) on the parameters A ( u * , u * , u ) and η ;
(2)
a first attempt towards the instability of the solutions of Equation (1) for certain values of the two classes of parameters ( A ( u * , u * , u ) , η ) and ( A ˜ ( u * , u * , u ) , η ˜ ) .
The first result will be a proof of the continuous dependence of solutions to Equation (1) with respect to the parameters.
Let Θ ( u * , u * , u ) be the function defined as:
Θ ( u * , u * , u ) : = | η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) | .
Theorem 1.
Let f ( t , u ) , f ˜ ( t , u ) C ( 0 , + ) × D u ; L 1 ( D u ) K ( D u ) the solutions to problems (9) and (10), respectively. Assume that Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) . If there exist δ , δ ^ > 0 such that | η η ˜ | < δ and Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) δ ^ then, for all T > 0 :
f ( t , u ) f ˜ ( t , u ) C ( 0 , T ) × D u ; L 1 ( D u ) δ + δ ^ T e ( 2 η + η ˜ ) T .
Proof. 
Integrating Equation (1) from 0 to t and recalling that f ( 0 , u ) = f 0 ( u ) in mind, one has
f ( t , u ) = f 0 ( u ) + 0 t J [ f , f ] ( τ , u ) d τ F 0 t u ( 1 u E 1 [ f ] ( τ ) ) f ( τ , u ) d τ ,
where
0 t J [ f , f ] ( τ , u ) d τ = 0 t G [ f , f ] ( τ , u ) d τ 0 t L [ f , f ] ( τ , u ) d τ = 0 t D u × D u η A ( u * , u * , u ) f ( t , u * ) f ( t , u * ) d u * d u * d τ η 0 t f ( τ , u ) d τ .
The integral expression (12) for the solutions f ( t , u ) and f ˜ ( t , u ) becomes
f ( t , u ) = f 0 ( u ) + 0 t D u × D u η A ( u * , u * , u ) f ( t , u * ) f ( t , u * ) d u * d u * d τ η 0 t f ( τ , u ) d τ F 0 t u ( 1 u E 1 [ f ] ( τ ) ) f ( τ , u ) d τ ,
and
f ˜ ( t , u ) = f 0 ( u ) + 0 t D u × D u η ˜ A ˜ ( u * , u * , u ) f ˜ ( t , u * ) f ˜ ( t , u * ) d u * d u * d τ η ˜ 0 t f ˜ ( τ , u ) d τ F 0 t u ( 1 u E 1 [ f ˜ ] ( τ ) ) f ˜ ( τ , u ) d τ
respectively.
Subtracting side by side Equation (14) from Equation (13), we readily get
f ( t , u ) f ˜ ( t , u ) = = 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ + 0 t η ˜ f ˜ ( τ , u ) η f ( τ , u ) d τ + F 0 t u f ˜ ( τ , u ) f ( τ , u ) + f ( τ , u ) u E 1 [ f ] ( τ ) f ˜ ( τ , u ) u E 1 [ f ˜ ] ( τ ) d τ
which leads immediately to the estimate
| f ( t , u ) f ˜ ( t , u ) | 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ + η ˜ 0 t f ˜ ( τ , u ) d τ η 0 t f ( τ , u ) d τ + F 0 t u f ˜ ( τ , u ) f ( τ , u ) + f ( τ , u ) u E 1 [ f ] ( τ ) f ˜ ( τ , u ) u E 1 [ f ˜ ] ( τ ) d τ .
Since f ( t , u ) = f ˜ ( t , u ) = 0 for u D u , the third term on the right-hand side of inequality (16) vanishes, so that we obtain the relation
| f ( t , u ) f ˜ ( t , u ) | 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ + 0 t η ˜ f ˜ ( τ , u ) η f ( τ , u ) d τ
and, by straightforward calculations, one finds that the first term on the right-hand side of (17) can be estimated as follows:
0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ 0 t D u × D u η f ( τ , u * ) A ( u * , u * , u ) | f ( τ , u * ) f ˜ ( τ , u * ) | d u * d u * d τ + 0 t D u × D u η ˜ f ˜ ( τ , u * ) A ˜ ( u * , u * , u ) | f ˜ ( τ , u * ) f ( τ , u * ) | d u * d u * d τ + 0 t D u × D u f ˜ ( τ , u * ) f ( τ , u * ) η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) d u * d u * d τ .
Now, using inequality (18) and integrating both sides of relation (17) on D u , we get
f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) η 0 t f ( τ , u ) f ˜ ( τ , u ) L 1 ( D u ) d τ + η ˜ 0 t f ( τ , u ) f ˜ ( τ , u ) L 1 ( D u ) d τ + Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) + 0 t D u | η ˜ f ˜ ( τ , u ) η f ( τ , u ) | d u d τ .
Since:
| η ˜ f ˜ ( τ , u ) η f ( τ , u ) | = | η ˜ f ˜ ( τ , u ) η f ˜ ( τ , u ) + η f ˜ ( τ , u ) η f ( τ , u ) | f ˜ ( τ , u ) | η ˜ η | + η | f ˜ ( τ , u ) f ( τ , u ) | ,
relation (19), bearing in mind that | η η ˜ | < δ and Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) δ ^ , may be rewritten in the form
f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) 0 t ( η + η ˜ ) f ( τ , u ) f ˜ ( τ , u ) L 1 ( D u ) d τ + Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) + | η ˜ η | t + η 0 t f ( τ , u ) f ˜ ( τ , u ) L 1 ( D u ) d τ 0 t ( 2 η + η ˜ ) f ( τ , u ) f ˜ ( τ , u ) L 1 ( D U ) d τ + δ + δ ^ t .
By Grönwall’s inequaility [36],
f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) ( δ + δ ^ ) t e ( 2 η + η ˜ ) t .
By (20) and (21), relation (11) is proved, i.e., for T > 0 ,
f ( t , u ) f ˜ ( t , u ) C ( 0 , T ) × D u ; L 1 ( D u ) ( δ + δ ^ ) T e ( 2 η + η ˜ ) T .
Remark 6.
It is worth pointing out that
1.
the assumption | η η ˜ | < δ is an estimate of the distance between the interaction rates;
2.
the assumption Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) δ ^ is an estimate on the distance between the transition probability densities, “weighted” by the interaction rates.
Remark 7.
The conclusion (11) of Theorem 1 ensures the continuous dependence of solutions on the parameters A ( u * , u * , u ) and η. Indeed,
f ( t , u ) f ˜ ( t , u ) C ( 0 , + ) × D u ; L 1 ( D u ) δ , δ ^ 0 0 .

3.2. The Discrete Activity Framework

This section aims to prove the continuous dependence of the solutions of Equation (4) on the parameters η h k and B h k i when f 0 , F ( t ) = F , constant in time, and T > 0 are fixed.
Let f ( t ) = f 1 ( t ) , f 2 ( t ) , , f n ( t ) , f ^ ( t ) = f ^ 1 ( t ) , f ^ 2 ( t ) , , f ^ n ( t ) be the solutions of the systems
d f i d t ( t ) = J i [ f ] ( t ) + F i ( t ) i = 1 n u i 2 J i [ f ] + F i E 2 [ f ] f i ( t ) t [ 0 , T ] f ( 0 ) = f 0 ,
and
d f i d t ( t ) = J ^ i f ] ( t ) + F i ( t ) i = 1 n u i 2 J ^ i f ] + F i E 2 [ f ] f i ( t ) t [ 0 , T ] f ( 0 ) = f 0 ,
respectively; the operators J [ f ] and J ^ [ f ] are defined by the parameters η h k , B h k i , and η ^ h k , B ^ h k i , respectively.
The following stability result holds.
Theorem 2.
Let f ( t ) = f 1 ( t ) , f 2 ( t ) , , f n ( t ) , f ^ ( t ) = ( f ^ 1 ( t ) , f ^ 2 ( t ) , , f ^ n ( t ) ) be the solutions of (12) and (13), respectively. Assume η h k η , η ^ h k η ^ , for h , k { 1 , 2 , , n } , and F i F , for i { 1 , 2 , , n } , for η , η ^ , F > 0 . If Λ : = i = 1 n h , k = 1 n η h k B h k i η ^ h k B ^ h k i , then
max t [ 0 , T ] f ( t ) f ^ ( t ) 1 Λ T e η + η ^ + i = 1 n u i 2 F i T
where
f ( t ) f ^ ( t ) 1 : = i = 1 n | f i ( t ) f ^ i ( t ) | .
Proof. 
Bearing assumption H4 in mind, and integrating Equations (22) and (23) on [ 0 , t ] , we get
f i ( t ) = f i 0 + 0 t J i [ f ] ( t ) + F i i = 1 n u i 2 J i [ f ] ( t ) + F i f i ( t ) d t ,
and
f ^ i ( t ) = f i 0 + 0 t J ^ i f ] ( t ) + F i i = 1 n u i 2 J ^ i [ f ^ ] ( t ) + F i f ^ i ( t ) d t
for i { 1 , 2 , , n } . Now, subtracting (26) from (25), we find
f i ( t ) f ^ i ( t ) = = 0 t J i [ f ] ( t ) J ^ i [ f ^ ] ( t ) d t 0 t i = 1 n u i 2 J i [ f ] ( t ) + F i f i ( t ) i = 1 n u i 2 J ^ i [ f ^ ] ( t ) + F i f ^ i ( t ) d t .
By taking the side-by-side sum on i { 1 , 2 , , n } of these last relations, we arrive at
i = 1 n f i ( t ) f ^ i ( t ) 0 t i = 1 n J i [ f ] ( t ) J ^ i [ f ^ ] ( t ) + 0 t i = 1 n i = 1 n u i 2 J i [ f ] ( t ) + F i f i ( t ) i = 1 n u i 2 J ^ i [ f ^ ] ( t ) + F i f ^ i ( t ) d t .
Now, first of all, observe that
0 t i = 1 n | J i [ f ] ( t ) J ^ i [ f ^ ] ( t ) | 0 t i = 1 n h , k = 1 n η h k B h k i f h ( t ) f k ( t ) η ^ h k B ^ h k i f ^ h ( t ) f ^ k ( t ) d t + 0 t i = 1 n f i ( t ) k = 1 n η i k f k ( t ) f ^ i ( t ) k = 1 n η ^ i k f ^ k ( t ) .
The integrand in the first term on the right-hand side of (29) can be estimated as follows:
η h k B h k i f h ( t ) f k ( t ) η ^ h k B ^ h k i f ^ h ( t ) f ^ k ( t ) = = η h k B h k i f h ( t ) f k ( t ) B h k i η h k f h ( t ) f ^ k ( t ) + B h k i η h k f h ( t ) f ^ k ( t ) η ^ h k B ^ h k i f ^ h ( t ) f ^ k ( t ) η h k B h k i f h ( t ) f k ( t ) f ^ k ( t ) + f ^ k ( t ) η h k B h k i f h ( t ) η ^ h k B ^ h k i f ^ h ( t ) η h k B h k i f h ( t ) f k ( t ) f ^ k ( t ) + f ^ k ( t ) η h k B h k i f h ( t ) η ^ h k B ^ h k i f h ( t ) + η ^ h k B ^ h k i f h ( t ) η ^ h k B ^ h k i f ^ h ( t ) η h k B h k i f h ( t ) f k ( t ) f ^ k ( t ) + f ^ k ( t ) f h ( t ) η h k B h k i η ^ h k B ^ h k i + η ^ h k B ^ h k i f h ( t ) f ^ h ( t ) .
By using this estimate, the first integral on the right-hand side of (29) turns out to be majorized as follows:
0 t i = 1 n h , k = 1 n η h k B h k i f h ( t ) f k ( t ) η ^ h k B ^ h k i f ^ h ( t ) f ^ k ( t ) d t 0 t i = 1 n h , k = 1 n η h k B h k i f h ( t ) f k ( t ) f ^ k ( t ) d t + 0 t i = 1 n h , k = 1 n η ^ h k B ^ h k i f ^ k ( t ) f h ( t ) f ^ h ( t ) d t + 0 t h , k = 1 n f ^ k ( t ) f h ( t ) i = 1 n η h k B h k i η ^ h k B ^ h k i d t η 0 t f ( t ) f ^ ( t ) 1 d t + η ^ 0 t f ( t ) f ^ ( t ) 1 d t + Λ t .
As far as the second term on the right-hand side of (28) is concerned, one has
i = 1 n u i 2 J i [ f ] ( t ) = i = 1 n u i 2 G i [ f ] ( t ) L i [ f ( t ) = i = 1 n u i 2 h , k = 1 n η h k B h k i f h ( t ) f k ( t ) f i ( t ) k = 1 n η i k f k ( t ) = h , k = 1 n i = 1 n u i 2 B h k i η h k f h ( t ) f k ( t ) i = 1 n u i 2 f i ( t ) k = 1 n η i k f k ( t ) = 0
which in turn implies
0 t i = 1 n i = 1 n u i 2 J i [ f ] ( t ) + F i f i ( t ) i = 1 n u i 2 J ^ i [ f ^ ] ( t ) + F i f ^ i ( t ) d t 0 t i = 1 n u i 2 F i f ( t ) f ^ ( t ) 1 d t .
Finally, by using relations (31) and (33), inequality (28) becomes
f ( t ) f ^ ( t ) 1 0 t η h + η ^ k + i = 1 n u i 2 F i f ( t ) f ^ ( t ) 1 d t + Λ t
and now Grönwall’s inequaility [36] yields
f ( t ) f ^ ( t ) 1 Λ T e η h + η ^ k + i = 1 n u i 2 F i t
leading at once to (24). □
Remark 8.
The conclusion of Theorem 2 is the continuous dependence of solution of Equation (4) on the parameters of the system, i.e., the interaction rate η h k and the transition probability density B h k i . In fact,
max t [ 0 , T ] f ( t ) f ^ ( t ) 1 Λ 0 0 .
Remark 9.
The coefficient Λ defined in Theorem 2 is a first estimate of the distance between the two classes of parameters, i.e., ( η h k , B h k i ) and ( η ^ h k , B ^ h k i ) .

4. A First Attempt towards the Instability with Respect to the Parameters

4.1. The Continuous Activity Framework

This section aims to give a first result about instability of solutions of Equation (1) with respect to the parameters, interaction rate η , and transition probability density A ( u * , u * , u ) .
Theorem 3.
Let f ( t , u ) , f ˜ ( t , u ) C ( 0 , + ) × D u ; L 1 ( D u ) K ( D u ) be the solutions to problems (9) and (10), respectively. Assume that Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) . If there exist two constants M 1 and M 1 ^ , with M 1 > M 1 ^ such that | η η ˜ | > M 1 and Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) M 1 ^ , then, for all T > 0 :
f ( t , u ) f ˜ ( t , u ) C ( 0 , T ) × D u ; L 1 ( D u ) ( M 1 M 1 ^ ) 1 + ( η + η ˜ ) T T > 0 .
Proof. 
As in Theorem 1, by using the integral formulation of (1) and by straightforward calculations, one has
| f ( t , u ) f ˜ ( t , u ) | = | 0 t η ˜ f ˜ ( τ , u ) η f ( τ , u ) d τ + 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ + F 0 t u ( f ˜ ( τ , u ) f ( τ , u ) + f ( τ , u ) u E 1 [ f ] ( τ ) f ˜ ( τ , u ) u E 1 [ f ˜ ] ( τ ) ) d τ | 0 t η ˜ f ˜ ( τ , u ) η f ( τ , u ) d τ 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ F 0 t u ( f ˜ ( τ , u ) f ( τ , u ) + f ( τ , u ) u E 1 [ f ] ( τ ) f ˜ ( τ , u ) u E 1 [ f ˜ ] ( τ ) ) d τ
whence, by integrating the (36) on D u , we obtain
f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) D u 0 t η ˜ f ˜ ( τ , u ) η f ( τ , u ) d τ d u D u 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ d u F D u 0 t u ( f ˜ ( τ , u ) f ( τ , u ) + f ( τ , u ) u E 1 [ f ] ( τ ) f ˜ ( τ , u ) u E 1 [ f ˜ ] ( τ ) ) d τ .
Since f ( t , u ) = f ˜ ( t , u ) = 0 for u D u , the third term on the right-hand side of the (37) vanishes.
The first term at the right-hand side of inequality (37) is estimated as follows:
D u 0 t η ˜ f ˜ ( τ , u ) η f ( τ , u ) d τ d u 0 t D u η ˜ f ˜ ( τ , u ) η f ( τ , u ) d u d τ = | η ˜ η | t .
Consider now the second term of the right-hand side of inequality (37). First of all,
D u 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ d u D u 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ d u .
By straightforward calculations,
η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) = = | η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η A ( u * , u * , u ) f ( τ , u * ) f ˜ ( τ , u * ) + η A ( u * , u * , u ) f ( τ , u * ) f ˜ ( τ , u * ) + η ˜ A ˜ ( u * , u * , u ) f ( τ , u * ) f ˜ ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ( τ , u * ) f ˜ ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) | = | η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) f ˜ ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ( τ , u * f ( τ , u * ) ) + f ( τ , u * ) f ˜ ( τ , u * ) η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) | η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) f ˜ ( τ , u * ) + η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ( τ , u * ) f ˜ ( τ , u * ) + f ( τ , u * ) f ˜ ( τ , u * ) η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) .
In virtue of inequalities (39) and (40),
D u 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ d u η 0 t D u f ( τ , u * ) D u | f ( τ , u * ) f ˜ ( τ , u * ) | D u A ( u * , u * , u ) d u d u * d u * d τ η ˜ 0 t D u f ˜ ( τ , u * ) D u | f ( τ , u * ) f ˜ ( τ , u * ) | D u A ( u * , u * , u ) d u d u * d u * d τ 0 t D u D u × D u f ( τ , u * ) f ˜ ( τ , u * ) η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) d u * d u * d u d τ = η t f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) η ˜ t f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) 0 t D u D u × D u f ( τ , u * ) f ˜ ( τ , u * ) η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) d u * d u * d u d τ
and, using Hölder’s inequality,
0 t D u D u × D u f ˜ ( τ , u * ) f ( τ , u * ) η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) d u * d u * d u d τ 0 t D u D u f ˜ ( τ , u * ) max u * D u f ( τ , u * ) D u η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) d u * d u * d u d τ 0 t D u D u max u * D u f ( τ , u * ) D u η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) d u * d u * d u d τ 0 t D u × D u × D u η A ( u * , u * , u ) η ˜ A ˜ ( u * , u * , u ) d u * d u * d u = Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) t .
Thanks to relations (41) and (42), inequality (39) becomes
D u 0 t D u × D u η A ( u * , u * , u ) f ( τ , u * ) f ( τ , u * ) η ˜ A ˜ ( u * , u * , u ) f ˜ ( τ , u * ) f ˜ ( τ , u * ) d u * d u * d τ d u ( η + η ˜ ) t f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) t .
Finally, by (38) and (43), inequality (37) yields
1 + ( η + η ˜ ) t f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) | η ˜ η | t Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) t .
Then:
f ( t , u ) f ˜ ( t , u ) L 1 ( D u ) | η ˜ η | Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) 1 + ( η + η ˜ ) t t .
Relation (35) is then proved by using the (44), and keeping in mind the fact that | η η ˜ | > M 1 and Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) M 1 ^ , with M 1 > M ^ 1 :
f ( t , u ) f ˜ ( t , u ) C ( 0 , T ) × D u ; L 1 ( D u ) | η ˜ η | Θ ( u * , u * , u ) L 1 ( D u × D u × D u ) 1 + ( η + η ˜ ) T T ( M 1 M 1 ^ ) 1 + ( η + η ˜ ) T T .
Remark 10.
In Theorem 3, the instability is related to the variation of the interaction rate.
Remark 11.
For instance, if D u = [ 0 , 1 2 ] is taken into account with A = A ˜ , then the right-hand side of relation (35) is strictly positive, so that the instability of the solutions follows at once.

4.2. The Discrete Activity Framework

In this section, we want to outline a first step of a study of instability in the discrete framework (4). This is an important issue in view of future numerical analysis.
Theorem 4.
Let f ( t ) = f 1 ( t ) , f 2 ( t ) , , f n ( t ) , f ^ ( t ) = ( f ^ 1 ( t ) , f ^ 2 ( t ) , , f ^ n ( t ) ) be the solutions of Equations (22) and (23), respectively. Let η , η ^ , F > = 0 such that η h k η , η ^ h k η ^ , for h , k { 1 , 2 , , n } , and F i F , for i { 1 , 2 , , n } . Furthermore, let
Γ ( t ) : = min f , f ^ ( C ( [ 0 , T ] ) ) n i = 1 n 0 t h , k = 1 n η h k B h k i f h ( t ) f k ( t ) η ^ h k B ^ h k i f ^ h ( t ) f ^ k ( t ) .
Then,
max [ 0 , T ] f ( t ) f ^ ( t ) 1 max [ 0 , T ] Γ ( t ) ( η + η ^ ) t 1 + i = 1 n u i 2 F i t .
Proof. 
Bearing the (27) in mind, and by using the (32), straightforward calculations show, for i { 1 , 2 , , n } :
f i ( t ) f ^ i ( t ) 0 t h , k = 1 n η h k B h k i f h ( s ) f k ( s ) η ^ h k B ^ h k i f ^ h ( s ) f ^ k ( s ) d s 0 t f ^ i ( s ) k = 1 n η ^ i k f ^ k ( s ) f i ( s ) k = 1 n η i k f k ( s ) d s i = 1 n u i 2 F i 0 t f ^ i ( s ) f i ( s ) d s .
By taking the sum on i 1 , 2 , , n of relations (46), we find
f ( t ) f ^ ( t ) 1 i = 1 n 0 t h , k = 1 n η h k B h k i f h ( s ) f k ( s ) η ^ h k B ^ h k i f ^ h ( s ) f ^ k ( s ) d s i = 1 n 0 t f ^ i ( s ) k = 1 n η ^ i k f ^ k ( s ) f i ( s ) k = 1 n η i k f k ( s ) d s i = 1 n i = 1 n u i 2 F i 0 t f ^ i ( s ) f i ( s ) d s .
Now, observe that
i = 1 n 0 t f ^ i ( s ) k = 1 n η ^ i k f ^ k ( s ) f i ( s ) k = 1 n η i k f k ( s ) d s ( s ) 0 t i = 1 n f ^ i ( s ) k = 1 n η ^ i k f k ^ ( s ) d s + 0 t i = 1 n f i ( s ) k = 1 n η i k f k ( s ) d s ( η + η ^ ) t ,
and
i = 1 n i = 1 n u i 2 F i 0 t f ^ i ( s ) f i ( s ) d t i = 1 n u i 2 F i 0 t f ( t ) f ^ ( t ) 1 .
Using these two last relations, inequality (47) may be rewritten in the form
f ( t ) f ^ ( t ) 1 Γ ( t ) ( η + η ^ ) t i = 1 n u i 2 F i 0 t f ( t ) f ^ ( t ) 1 ,
so that
f ( t ) f ^ ( t ) 1 Γ ( t ) ( η + η ^ ) t 1 + i = 1 n u i 2 F i t ,
and inequality (45) is achieved. □
Remark 12.
By using inequality (45), we see that, if Γ ( t ) > ( η + η ^ ) t , then an instability appears in the framework (4). In addition, it is important to note that this is a condition involving the parameters of the system, i.e., interaction rate and transition probability density.

4.3. Numerical Simulations

This section aims to present some numerical simulations in the framework described by (22). Specifically, the parameters of the system, i.e., interaction rate and transition probability, acquire different values. All the simulations that follow have been performed by using the routine Ode45 of MatLab.
Let n = 3 , which is three functional subsystems that are taken into account. The initial data are the vector:
f 0 = 3 / 8 , 1 / 2 , 1 / 8 .
The interaction rate parameter has the following form:
η h k = exp η | h k | , η > 0 .
In the first set of simulations, the transition probability is constant, whether the interaction rate varies. Specifically, three cases are considered:
  • η = 1 ;
  • η = 3 ;
  • η = 6 .
The solution f ( t ) = f 1 ( t ) , f 2 ( t ) , f 3 ( t ) is of course different from value to value of the interaction rate. Specifically, Figure 1 shows the three plots of the solution f ( t ) respectively corresponding to the three different values of η listed above. In addition, Figure 2 offers a comparison between the solutions corresponding to the values η = 3 and η = 3 , 2 , respectively.
In the second set of simulations, the interaction rate is constant, while the transition probability density acquires different real values. Precisely, Figure 3 shows the three plots of the solution f ( t ) respectively corresponding to three different values of B h k i . The considered cases are:
  • B h k i = c i h k 1 s g | h i | , i , h , k { 1 , 2 , , n } ,
    where g is a non-increasing function of | h i | and s, and the parameters c i h k , for i , h , k { 1 , 2 , , n } , are positive real numbers, depending on the particular system taken into account;
  • B ^ h k i that differs from B h k i only for h = 3 and k = 2 ;
  • B h k i uniform.
Furthermore, Figure 4 shows the solutions corresponding to the values B h k i and B ^ h k i in the same plot in order to compare their behaviors in time.
It is worth being stressed that the shape of solution strictly depends on the value of the parameters of the system (see Figure 1 and Figure 3) as the results reported in Section 4.1 and Section 4.2 show for both the continuous and the discrete framework. Moreover, bearing the Figure 2 and Figure 4 in mind, a small perturbation of a parameter may determine that the related solution has the same shape, but they are not so “close” to each other.

5. Conclusions and Research Perspectives

The results proved in Section 3 and Section 4 are, in some sense, complementary. The former shows that the difference between two arbitrary solutions of Equation (1), corresponding to two different sets of parameters, i.e., different systems of interaction rates and different probability distributions on the results of interactions, varies with a suitable measure of the difference between the systems of parameters, and the variation is continuous; the latter shows that, if the difference between the interaction rates is sufficiently large in a suitable sense, then the corresponding solutions—though starting from the same initial value—will move at once away from each other and will stay apart at any future time, i.e., their distance has a constant positive lower bound. In addition, the numerical simulations plotted in the figures shown in Section 4.3 seem to give a good visual counterpart of this result.
As far as we are aware, no similar results have been previously reported in the literature about KTAP, perhaps because the study of the dependence of solutions on the perturbations of parameters (interaction rates and transition probabilities) seems to be too difficult in relation to its relevance for applications, so that tackling it is considered as an almost useless effort. However, on the contrary, results like the ones found and reported in the present work are probably intended for becoming of the greatest relevance for applications, with special concern with social and economic sciences. In this connection, we can observe that economic interactions in any human society are ruled by the government: in a country in which some commercial transactions are allowed, they will produce exchanges of goods and money, with a subsequent modification of the distribution of wealth; but, in another country, where the same transactions are forbidden, the interaction rate referred to them is zero, and we must expect that the distribution of wealth could not be modified by these transactions, regardless of the values of transition probabilities that are allowed to be the same in both cases. This remark has worked as a suggestion of a search for instability results of the kind of Theorems 3 and 4. Of course, these Theorems cannot be considered as more than a first step on the way towards much more general instability results, for at least the good reason that they only refer to the very special case in which the interaction rates are constant with respect to the couples of states. Accordingly, this research about instability requires to be deepened along at least three lines, which will be the object of future work.
As laid out in the Introduction, KTAP is not a theory or simply a model, but a whole scheme of models to describe and—above all— predict the behavior of complex systems. In addition, as a matter of fact, our prime scope is its application to human collectivities, in order to suggest some ways to solve the problems raised by many and diffused bad mental habits that control not only human behaviors but also the criteria according to which legislators decide the (inter-)actions that can be allowed and the (inter-)actions that must be forbidden. Laws can modify both interaction rates and transition probabilities, so a complete and detailed view of the behaviors they produce could avoid that past mistakes from being repeated in the future.
In this line of thought, first of all, one should find possible conditions of instability in the quite general case in which interaction rates are arbitrary functions defined on D u × D u : from a purely mathematical viewpoint, this will require in turn a suitable definition of their distance.
Next, one has to find possible instability conditions on the transition probabilities, also in the case in which the interaction rates are left unchanged.
Finally, one has to study the reciprocal influence between the perturbation of interaction rates and the perturbation of transition probabilities. In this connection, it should be noted that Theorem 3 already furnishes a first hint in this direction.
These three lines of search give good and hopefully—on a pragmatic ground—useful perspectives for the development of the study started and reported in the present paper.

Author Contributions

Conceptualization, B.C. and M.M.; methodology, B.C. and M.M.; software, B.C. and M.M.; validation, B.C. and M.M.; formal analysis, B.C. and M.M.; writing—original draft preparation, M.M.; writing—review and editing, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Marco Menale is supported by the Research Project “ANDROIDS” (AutoNomous DiscoveRy Of depressIve Disorder Signs) for VALERE (VAnviteLli pEr la RicErca), developed by Università degli Studi della Campania “L. Vanvitelli”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bar-Yam, Y. Dynamics of Complex Systems; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  2. Bianca, C. Modeling complex systems by functional subsystems representation and thermostatted-KTAP methods. Appl. Math. Inf. Sci. 2021, 6, 495–499. [Google Scholar]
  3. Cilliers, P. Complexity and Postmodernism: Understanding Complex Systems. S. Afr. J. Philos. 1999, 18, 275–278. [Google Scholar] [CrossRef]
  4. Morriss, G.P.; Dettmann, C.P. Thermostats: Analysis and application. Chaos Interdiscip. J. Nonlinear Sci. 1998, 8, 321–336. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bobylev, A.V.; Cercignani, C. Exact eternal solutions of the Boltzmann equation. J. Stat. Phys. 2002, 106, 1019–1038. [Google Scholar] [CrossRef]
  6. Cercignani, C. The boltzmann equation. In The Boltzmann Equation and Its Applications; Springer: New York, NY, USA, 1988; pp. 40–103. [Google Scholar]
  7. Cercignani, C.; Gabetta, E. Transport Phenomena and Kinetic Theory: Applications to Gases, Semiconductors, Photons, and Biological Systems; Springer Science & Business Media: Basel, Switzerland, 2007. [Google Scholar]
  8. Cercignani, C.; Illner, R.; Pulvirenti, M. The Mathematical Theory of Dilute Gases; Springer Science & Business Media: Basel, Switzerland, 2013; Volume 106. [Google Scholar]
  9. Bonabeau, E.; Theraulaz, G.; Deneubourg, J.L. Mathematical model of self-organizing hierarchies in animal societies. Bull. Math. Biol. 1996, 58, 661–717. [Google Scholar] [CrossRef]
  10. Thieme, H.R. Mathematics in Population Biology; Princeton University Press: Princeton, NJ, USA, 2018; Volume 12. [Google Scholar]
  11. Giorno, V.; Román-Román, P.; Spina, S.; Torres-Ruiz, F. Estimating a non-homogeneous Gompertz process with jumps as model of tumor dynamics. Comput. Stat. Data Anal. 2017, 107, 18–31. [Google Scholar] [CrossRef]
  12. Masurel, L.; Bianca, C.; Lemarchand, A. On the learning control effects in the cancer-immune system competition. Phys. A Stat. Mech. Its Appl. 2018, 506, 462–475. [Google Scholar] [CrossRef] [Green Version]
  13. Pappalardo, F.; Palladini, A.; Pennisi, M.; Castiglione, F.; Motta, S. Mathematical and computational models in tumor immunology. Math. Model. Nat. Phenom. 2012, 7, 186–203. [Google Scholar] [CrossRef] [Green Version]
  14. Poleszczuk, J.; Macklin, P.; Enderling, H. Agent-based modeling of cancer stem cell driven solid tumor growth. In Stem Cell Heterogeneity; Human Press: New York, NY, USA, 2016; pp. 335–346. [Google Scholar]
  15. Spina, S.; Giorno, V.; Román-Román, P.; Torres-Ruiz, F. A stochastic model of cancer growth subject to an intermittent treatment with combined effects: Reduction in tumor size and rise in growth rate. Bull. Math. Biol. 2014, 76, 2711–2736. [Google Scholar] [CrossRef]
  16. Bianca, C.; Kombargi, A. On the inverse problem for thermostatted kinetic models with application to the financial market. Appl. Math. Inf. Sci. 2017, 11, 1463–1471. [Google Scholar] [CrossRef]
  17. Bisi, M.; Spiga, G.; Toscani, G. Kinetic models of conservative economies with wealth redistribution. Commun. Math. Sci. 2009, 7, 901–916. [Google Scholar] [CrossRef] [Green Version]
  18. Carbonaro, B.; Serra, N. Towards mathematical models in psychology: A stochastic description of human feelings. Math. Model. Methods Appl. Sci. 2002, 12, 1453–1490. [Google Scholar] [CrossRef]
  19. Carbonaro, B.; Giordano, C. A second step towards a stochastic mathematical description of human feelings. Math. Comput. Model. 2005, 41, 587–614. [Google Scholar] [CrossRef]
  20. Bronson, R.; Jacobson, C. Modeling the dynamics of social systems. Comput. Math. Appl. 1990, 19, 35–42. [Google Scholar] [CrossRef] [Green Version]
  21. Buonomo, B.; Della Marca, R. Modelling information-dependent social behaviors in response to lockdowns: The case of COVID-19 epidemic in Italy. medRxiv 2020. [Google Scholar] [CrossRef]
  22. Fryer, R.G., Jr.; Roland, G. A model of social interactions and endogenous poverty traps. Ration. Soc. 2007, 19, 335–366. [Google Scholar] [CrossRef]
  23. Kacperski, K. Opinion formation model with strong leader and external impact: A mean field approach. Phys. A Stat. Mech. Its Appl. 1999, 269, 511–526. [Google Scholar] [CrossRef]
  24. Bianca, C.; Carbonaro, B.; Menale, M. On the Cauchy Problem of Vectorial Thermostatted Kinetic Frameworks. Symmetry 2020, 12, 517. [Google Scholar] [CrossRef] [Green Version]
  25. Bianca, C. An existence and uniqueness theorem to the Cauchy problem for thermostatted-KTAP models. Int. J. Math. Anal. 2012, 6, 813–824. [Google Scholar]
  26. Bianca, C.; Mogno, C. Qualitative analysis of a discrete thermostatted kinetic framework modeling complex adaptive systems. Commun. Nonlinear Sci. Numer. Simul. 2018, 54, 221–232. [Google Scholar] [CrossRef]
  27. Bianca, C.; Menale, M. A Convergence Theorem for the Nonequilibrium States in the Discrete Thermostatted Kinetic Theory. Mathematics 2019, 7, 673. [Google Scholar] [CrossRef] [Green Version]
  28. Bianca, C.; Menale, M. Existence and uniqueness of nonequilibrium stationary solutions in discrete thermostatted models. Commun. Nonlinear Sci. Numer. Simul. 2019, 73, 25–34. [Google Scholar] [CrossRef]
  29. Carbonaro, B.; Menale, M. Dependence on the Initial Data for the Continuous Thermostatted Framework. Mathematics 2019, 7, 612. [Google Scholar] [CrossRef] [Green Version]
  30. Carbonaro, B.; Menale, M. The mathematical analysis towards the dependence on the initial data for a discrete thermostatted kinetic framework for biological systems composed of interacting entities. AIMS Biophys. 2020, 7, 204. [Google Scholar]
  31. Bianca, C. Existence of stationary solutions in kinetic models with Gaussian thermostats. Math. Methods Appl. Sci. 2013, 2013, 1768–1775. [Google Scholar] [CrossRef]
  32. Bianca, C.; Menale, M. On the convergence towards nonequilibrium stationary states in thermostatted kinetic models. Math. Methods Appl. Sci. 2019, 42, 6624–6634. [Google Scholar] [CrossRef]
  33. Bianca, C. Thermostated kinetic equations as models for complex systems in physics and life sciences. Phys. Life Rev. 2012, 9, 359–399. [Google Scholar] [CrossRef] [PubMed]
  34. Bianca, C.; Mogno, C. Modelling pedestrian dynamics into a metro station by thermostatted kinetic theory methods. Math. Comput. Model. Dyn. Syst. 2018, 24, 207–235. [Google Scholar] [CrossRef]
  35. Bianca, C.; Menale, M. Mathematical Analysis of a Thermostatted Equation with a Discrete Real Activity Variable. Mathematics 2020, 8, 57. [Google Scholar] [CrossRef] [Green Version]
  36. Walter, W. Differential and Integral Inequalities; Springer-Verlag: Berlin/Heidelberg, Germany, 2012; Volume 55. [Google Scholar]
Figure 1. From top left to bottom η = 1 , η = 3 , γ = 6 . f 1 dot-dashed, f 2 dashed, f 3 full.
Figure 1. From top left to bottom η = 1 , η = 3 , γ = 6 . f 1 dot-dashed, f 2 dashed, f 3 full.
Axioms 10 00059 g001
Figure 2. In black the solution for η = 3 , in red the solution for η = 3 , 2 .
Figure 2. In black the solution for η = 3 , in red the solution for η = 3 , 2 .
Axioms 10 00059 g002
Figure 3. From top left to bottom B h k i , B ^ h k i perturbed and uniform distribution. f 1 dot-dashed, f 2 dashed, f 3 full.
Figure 3. From top left to bottom B h k i , B ^ h k i perturbed and uniform distribution. f 1 dot-dashed, f 2 dashed, f 3 full.
Axioms 10 00059 g003
Figure 4. In black the solution for B h k i , in red the solution for B ^ h k i .
Figure 4. In black the solution for B h k i , in red the solution for B ^ h k i .
Axioms 10 00059 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Carbonaro, B.; Menale, M. Towards the Dependence on Parameters for the Solution of the Thermostatted Kinetic Framework. Axioms 2021, 10, 59. https://doi.org/10.3390/axioms10020059

AMA Style

Carbonaro B, Menale M. Towards the Dependence on Parameters for the Solution of the Thermostatted Kinetic Framework. Axioms. 2021; 10(2):59. https://doi.org/10.3390/axioms10020059

Chicago/Turabian Style

Carbonaro, Bruno, and Marco Menale. 2021. "Towards the Dependence on Parameters for the Solution of the Thermostatted Kinetic Framework" Axioms 10, no. 2: 59. https://doi.org/10.3390/axioms10020059

APA Style

Carbonaro, B., & Menale, M. (2021). Towards the Dependence on Parameters for the Solution of the Thermostatted Kinetic Framework. Axioms, 10(2), 59. https://doi.org/10.3390/axioms10020059

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop