The Collective Behaviors of Self-Excitation InformationDiffusion Processes for a Large Number of Individuals

)e opinion dynamics is a complex and interesting process, especially for the systems with a large number of individuals. It is usually hard to describe the evolutionary features of these systems. In some previous works, it has been shown that the selfexcitation type model has superior performance in learning and predicting opinions. Following this line, we consider the selfexcitation opinion model and study the collective behaviors of the self-excitation model. We propose a Mckean–Vlasov-type integrodifferential equation to describe the asymptotic behaviors of the model and show that the introduced equation, by coupling with the initial distribution, has the ability of capturing the influence of the self-excitation process, which describes the mutually exciting and recurrent nature of individuals. We also find that the steady-state distribution is a “contraction” of the initial distribution in the linear and bounded confidence (DWmodel) interaction cases, which is different from the results of the model with nonself-excitation interaction.


Introduction
Expressing opinions and then influencing others are the primary forms of social behavior of people. e dynamics of public opinion has long been a major concern of social science and computational social science. ere has been an increasing interest on the perception and prediction of social opinion, such as quantitative investment firms which measure investor sentiment and trade using social media Karppi and Crawford [1] and prediction of election results [2].
Based on sentiment analysis with deep learning, there are many frameworks to perceive opinion on the social network. For the prediction of social opinion, there are some noticeable models proposed at the beginning of the 21 st century, such as the Sznajd model [3], Hegselmann-Krause model [4,5], Deffuant-Weisbuch model [6], and so on, which had been established to understand the co-ordinated movements of opinion as a group. However, these works held some limitations: (i) most of these models on opinion dynamics are theoretical and have not proved their effectiveness quantificationally. (ii) ese models are represented by a cellular automaton, which makes it difficult to analyze the global behavior of the dynamics.
More recently, some researchers began to study the models, which can provide more accurate predictions [7][8][9]. De et al. [9] proposed a framework of opinion dynamics named SLANT, providing accurate predictions of users' opinions. In the early models, the opinion and the states of the next moment only depend on the states of the last moment, and thus have nothing to do with more previous ones. It means that opinion formation is seen as a Markov process. However, in the real world, this assumption may not be precise. e model given in the study by DE et al. [9] indicates that self-excitation opinion models with exponential function weight have superior performance in learning and predicting opinions.
In many works [9][10][11][12][13], self-excitation models have shown efficient abilities to capture users' behavior pattern. ese models focus on the microbehaviors basing on cellular automata, which would not express the global dynamics of opinion evolution analytically.
In this paper, we study the self-excitation opinion dynamics on the large homogeneous network. We show that when the number of the individuals is very large, the opinion distribution of this system evolves according to a nonlinear partial differential equation (PDE) of Mckean-Vlasov type. Due to the dependence on the initial value in the local rules, the individuals with different initial values will have different dynamic equations; thus, the mean-field equation deriving process is different from the cases of the previous works [14,15]. We show that this PDE can be decomposed into two parts where the first part is the same as the Fokker-Planck model in the study by Toscani [16]. And the second part can capture the mutually exciting and recurrent nature of individual behaviors by an integral term related to the initial value.
is term leads to different steady state behaviors from the Markov model, such as the model given in the study by Toscani [16]. Generally, the stationary solutions of the Toscani type model will trend to a single-point support Dirac measure when noise parameters tend to 0, but the model with additional term makes the steady state to be a "contraction" of the initial distribution. e rest of this paper is organized as follows. Section 2 is for related work. In Section 3, we introduce the mathematical model of opinion dynamics we study and the main mean-field theorem in this paper. In Section 4, using this theorem, we consider the linear interaction model and the bounded confidence (DW) model and study the steady-state distribution of these models. Section 5 is the proof outline of this theorem. Section 6 is a summary. e detailed proof is in the appendix.

Information Diffusion Model.
e early theory of information diffusion in mathematical sociology was proposed by Sznajd in the study by Sznajd-Weron and Sznajd [3]. Sznajd's model is designed to explain the features of opinion dynamics, in which every individual is on a lattice and have a state 0 or 1 to express the opinion. Individuals update their opinion by the opinions of their neighbors. e numerical simulations of this model are also investigated by many researchers [17][18][19]. On the complete graph, Slanina and Lavicka [20] show that when the number of individuals is very large, the probability densities evolve according to the partial differential equations where P(m) is the probability distribution of opinion. Equation (1) [6] and Hegselmann-Krause model [4,5] deal with two different cases. Both these models are so-called bounded confidence models, i.e., the individuals can only be influenced by some others whose opinions are close enough to theirs. e characteristics of bounded confidence lead to the formation of communities, which coincides with the phenomenon in the real world. A very beautiful analysis of the convergence of the Deffuant-Weisbuch model is given in the study by Zhang et al. [21], which reveals the community structure of the model. Bhattacharyya [22] proves the convergence of the Hegselmann-Krause model. On the other hand, numerical simulations on complex networks [23] show that "the more heterogeneous the complex network is, the weaker the ability of polarization and consensus of the complex network will be." Models with both continuous opinion and discrete states are also considered, such as the SHIR model [24], where the individuals have a four-state variable to represent susceptible, hesitated, infected, and removed. Opinion dynamics of the social network and the impact on the recommendation are studied in the studies by Xiong et al. [25] and Li et al. [26].
Como and Fagnani [27] gave a general framework to prove that the scaling limits of the many-body continuous opinion system can be described by a measured value ordinary differential equation (ODE). In such models, every individual will change its opinion to ωX + (1 − ω)X′, a convex combination with the neighbor with probability k(X, X ′ ). When k(X, X ′ ) � 1 [0,ε] |(|X − X ′ |), it becomes the Deffuant-Weisbuch model. e key ideas of their work are the large deviation principle and some estimations based on the Lipschitz property, which are the same as the general method used to prove the chaotic propagation properties [28].
Based on the Boltzmann equation with granular gas like interactions, Toscani [16] introduced a collision model, in which the communication between individuals was considered as a collision with the following form: where function D and P describe the local relevance of the compromise and diffusion for a given opinion. Toscani also considered the quasi-invariant opinion limit for the Boltzmann-type equation. e main idea of that is to scale the interaction frequency, strength, and the diffusion in the integral equation, and then, the equation reduces to a Fokker-Planck type equation.
en, on the base of this model, considering the influence of the structure of the social network, Albi et al. [29] consider the kinetic opinion dynamics on the large scale networks evolving over time. e evolutions of both the network and the opinion were involved.

Self-Excitation Point Process in Social
Systems. Point processes type models are generally used to analyze the impact of events on the system such as CSMA/CD [30], financial data [31], and the dynamics of book sales [32]. A very succinct and effective analysis for the time series of daily views for videos on YouTube is introduced [33]. e authors showed that most of the video viewing record data could be described statistically by a Poisson process, but still, about 10% data show a point process with the instantaneous rate λ(t), where μ i is the number of potential viewers influenced directly by person i who views the video through the social network, and V(t) captures all spontaneous views that are not triggered by epidemic effects on the network. For the self-excitation social opinion models, De et al. [9] proposed an opinion model on a social network, in which every user has a latent opinion about the given topic and can post messages about the topic. e experiments on real data show that the self-excitation model, with the intensity depending on the messages sent from a neighbor in the past and the weight function having the form like e − ωt , performs much more better than the Poisson process (ω � 0) model.
Wang et al. [13] used point process to predicting user activities.
ey proposed a generic framework for point process prediction problem, and they used a mass transport equation to update the transport rate and compute the conditional mass function.

Mean Field eory.
Mean-field theory (MFT) is introduced to study many-body problem by using a single averaged effect to approximate the effect of all individuals. e core of MFT is to estimate the error of the mean-field approximation. If there is no interaction, all individuals are independent of each other, and then, the law of large numbers works and gives the mean-field approximation. In the general case, we still want to be able to use the law of large numbers. at is to say, if one picks a chaotic (i.i.d) initial distribution of particles, we hope that this distribution is still chaotic as time-evolving, which is the so-called propagation of chaos. e theory of propagation of chaos in the cases of Wiener noise is summed in the study by Sznitman [28]. If N particles with initial "chaotic" distribution u ⊗ N 0 satisfy the stochastic differential equation (SDE), the Mckean-Vlasov mean-field equation is where u is the law of X. e method to show this is to use the Lipschitz character of b and Gronwall's lemma. Mckean-Vlasov equation for SDE with Poisson jumps is in the study by Andreis et al. [15]. e authors used Burkholder-Davis-Gundy inequality for martingales. e main step in the study by Andreis et al. [15] is to use Doob-Meyer decomposition and estimate the compensated Poisson process.
Using the large deviation principle [34], Arous and Guionnet [35] studied the mean-field simplification dynamics for Langevin spin glass. e framework of their proof makes use of Sanov's theorem and Varadhan's lemma in the path space C[0, T] to get the rate function for the interaction system. e main obstacle to use these methods to prove MFT is the unboundedness of coefficients in SDE. Because of Varadhan's lemma, if we use a stop time to get the localization of SDE, we can get the so-called local large deviation principle. Dawsont and Gartner [36] gave some compactness criteria to convert the "local" result into "global" one.
ey also show that a Lyapunov function for the system of weakly interacting diffusions will let the compactness condition be satisfied. Puhalskii [37] set up the whole framework of local to global LDP and introduce the C-exponential tightness conditions.

Self-Excitation Opinion Model.
In this paper, we consider a simplified version of the model [9] and study the large number limit of the self-excitation opinion model with the following form: ere are N individuals with opinion x i , which is a function of the time t, such that where h(x j , h i ) is a interaction function, * is the convolution operator, k(t) � e − ωt is the triggering kernel, and N(t) is a Poisson process. is model has been previously studied [11,12]. e original model proposed in the study by De et al. [9] has the following form: where m j is the message sent by user j, α ij is the interaction intensity between i and j, b i � x i (0) is the initial opinion of user i, H j (t) is the history of events up to time t for user j, and k(t) is a triggering kernel. e values of message m come from a sentiment distribution related to the opinion x and As the models have been previously studied [11,12], we only consider the opinion part, the simplified model (7). In the case that h(x i , x j ) � h(x j ), where x j is the opinion of the individual j, De et al. and Wang et al. [9,11] showed the validity of the model in the study by DE et al. [9] on the real world data. And when h(x j , x i ) � x j − x i , the form of this model will degenerate into the form of the DW model [6] and Toscani model [16] when k(t) � 1.

Complexity 3
In studies by De et al. and Wang et al. [9,11], the selfexcitation opinion model with exponential function type weight was shown to be very effective in the opinion dynamics; therefore, we set k(t) � e − ωt , where ω is a constant. Due to the fact that differential of the convolution of two where there is a jump part in g, and dk(t) � − ωk(t)dt. We can use a stochastic differential equation (SDE) to describe x(t) [11].
It is also necessary to add σdW into the model to represent noise, where W is the Wiener noise with e above two equations are given in the study by Wang et al. [11]. Because of the dependence on the path of the history, it is apparently not a Markov process. We consider the limit that the population size tends to infinity. Let X N � (X 1 , X 2 , . . . , X N ) be the opinion variable of the N-particles system. X satisfies If we only consider the case ω � 0 and replace the Poisson Noise dN with dt and consider the form e mean-field theory of this equation is well studied [28]. e Mckean-Vlasov mean-field equation is In this paper, we first use Doob-Meyer decomposition for dN. en, we estimate the compensated Poisson process by Burkholder-Davis-Gundy inequality for martingales. is method helps us to reduce our equation to the form in (12). As for the dependence on the initial values, we regard the initial as a random media. Using the methods in the study by Pra and Hollander [38].
e Mckean-Vlasov equation can be obtained using Sanov's theorem and Varadhan's lemma in the path space C(0, T). In this paper, we use a direct estimation to prove the convergence to the Mckean-Vlasov equation.

Main Results.
e main result in this paper is the following theorem: Theorem 1. When N tends to infinity, there is a limit process X for X N in equation (11), with SDE: where P is the law of X, and b � X[0] is a random variable. And given the initial distribution μ, the law of X is where δ b is a Dirac measure with support on b.
From this theorem, we can study the asymptotic behavior of the self-excitation opinion model. In Section 4, we will study the steady state of this system by using Fokker-Planck equation and show the effects of self-excitation, which makes it impossible to give a perfect consensus, such that it can avoid Abelson's diversity puzzle [39].
In order to prove this theorem, we need some assumptions. Stemming from the physical meaning, we can assume that the support of initial distribution μ(b) is bounded, such as the uniform distribution on [0, 1]. Also, for the proof of Mckean-Vlasov limit, we need the function h(x) to have good enough properties to make where P is the solution of the Mckean-Vlasov equation.
If h is bounded, the first condition is satisfied trivially. In the appendix, we will give other examples of h that can make these conditions be satisfied. Also, we assume that h(x, y) is Lipschitz for both of the variables. However, this is not the necessary condition [14].
Most of these conditions can be relaxed, but we hope to prove our results without excessive technical details on tightness; therefore, these restrictions are imposed.

Mckean-Vlasov Equation.
For the equation we consider a simple case 4 Complexity en, we can see that where m � xdP(x), the mean value of x.
Since h is an odd function, it is easy to see that en, we have the Fokker-Planck equation associated with the Mckean-Vlasov process: Another example is the bounded confidence type model [4,6] where every agent interacts only within a certain level of confidence: where k(y − x) is a continuous function, and if Since, in this case, h(x) is bounded and Lipschitz, our assumptions are satisfied.
Note that for the case that h(y, x) � αy, which is used in the SLANT system [9], where As a sanity check, this formula matches (4) in the study by De et al. [9]. Solving this equation, we get [27] showed that when the population size tends to infinity, the limit behavior of the Deffuant-Weisbuch model can be described by a measurevalued ODE. And also, Toscani introduced a kinetic model of opinion formation [16]. In this model, the opinion is changed by the binary interaction of collision, and the dynamics of the opinion distribution is modeled by Boltzmann-type integrodifferential equation.

Discussion. Como and Fagnani
In these models, the dynamics of opinion distribution is depicted as follows: where μ t is the measure of opinion distribution, and ϕ is an arbitrary test function. e right side of (29) has the following form [27]: and in Boltzmann case [16], Intuitively, both models can be derived from the stochastic differential equation with Poisson jump interactions by the method mentioned in the study by Méléard [14]. Events that interact with other particles are regarded as a Poisson process.
For the system with h(y, x) � y − x, in the case that ω � 0, which is the nonself-excitation case, it degenerates into the general opinion dynamics model, or equally, the quasi-invariant opinion limit for Toscani's Boltzmann equation [16]: e ω term is a modification to consider the influence of the self-excitation process.
In this paper, we only consider the homogeneous network. However, using the method in the study by Pra and Hollander and Albi et al. [38,40], we can study the large system of agents interacting through a network with a given distribution of the number of connections by the following way: (1) Consider a random media c ∈ 0, 1, 2, . . . , c max , which is a discrete variable describing the number of connections. c should be a random variable to Complexity 5 represent the complex network. For example, c obeys power law distribution. (2) Construct a new interaction equation as in the study by Albi et al. [40].
en, it is the same as eorem 1 in the study by Pra and Hollander [38] such that in this random media, we can also prove a mean-field theorem. And there is a distribution P(c, x, t) about the variable c, opinion x, and time t. Using a similar argument as Albi et al. [40], we have obtained the stationary solutions.

Steady-State Distribution of the Linear Case.
Considering a simple case of initial distribution P 0 � 1/2δ − 10 + 1/2δ 10 . We will use this simple form to explain how the initial distribution is coupled to the equation. en, the equation will become e overall distribution P is equal to 1/2P − 10 { } + 1/2P 10 { } . e evolution process of this distribution is shown in in Figure 1(a) . And Figures 1(b) and 1(c) show the differences between two kinds of steady states.
As shown in Figure 1(b), for the initial distribution P 0 , the steady state distribution is normal. In fact, this is true for any initial distribution, and we have the following theorem.

Theorem 2.
In the case that ω � 0, for all type of initial distribution P 0 , we have the steady-state distribution P(∞) � 1/Ze (α(x− m)/σ) 2 , where Z is the normalization constant. This is because the form of PDE for P has the form The steady-state equation and this steady state is asymptotically stable, such that when t ⟶ ∞, the solution of the equation will trend to the steady state. It is easy to see that the steady state is a normal distribution as shown in Figure 1(b).
From Theorem 2, this ω � 0 model cannot describe the actual situation since, in this model, the final state of the opinion distribution has no community structures but a perfect consensus (in fact, this is true for any monotonic function h, since the steady state is the Gibbs distribution of a potential energy V(x), such that ∇V � h), which is so-called Abelson's diversity puzzle [39], a persistent research puzzle in the social sciences. Generally, this puzzle is solved by the bounded-confidence mode, such as HK [4] and DW [6] models. It is proved in the study by Como and Fagnani [27] for a large number of individuals case and by Zhang et al. [21] for the general case that the steady state of the DW model has the form However, in the case ω > 0, the equation has the steady state as in Figure 1(c) when the initial distribution has the form 1/2(δ − 10 + δ 10 ). In the general case, we have Theorem 3. In the case that ω > 0, given an initial distribution μ(x), we have the steady-state distribution The steady state of P b is easy to calculate: where Z is the normalization constant, so we have P(∞) � dμ(b)P b (∞).
Considering the limit behavior σ ⟶ 0, the steady state of Mckean-Vlasov equation when h(y, x) � y − x has the form So that even in the case that the variance of the noise trends to 0, this equation will not converge to a single Dirac measure, which is different from the ω � 0 case. Intuitively, the term ω(b − x) makes the opinion try not to deviate from the initial position too far, which can be considered as soft-bounded confidence.
In the case that h(y, x) � αy, which is used in the SLANT system in the study by De et al. [9], supposing ω > α, the steady state will have the similar form since the equation of which is where m � m(∞) � (α/ω − α)m(0) (this is from (25)), the final mean value of x. We have P(∞) � dμ(b)δ(x− (αm + ωb)/ω).

Steady-State Distribution of Bounded Confidence Models.
In the case that ω � 0 and h(y, x) � (y − x)k(y − x), the steady-state distribution is studied [27]. We can follow their methods to study the weak solution of our equation.
6 Complexity Proof. Considering equation (19) with ω � 0, let where k(y − x) is a bounded confidence kernel. In (a) of (42), we use the symmetry of P(x) and P(y) and the integration by parts. Considering the steady state and σ ⟶ 0 limit, we have d/dt(Q t ) � 0, yet if in the support of P(∞), there are two points x and y such that ‖x − y‖ < R and x ≠ y, and we have − (y − x) 2 k(y − x)P(y)P(x) < 0, which is a contradiction, so that our claim follows. □ Supposing ω > 0, it is hard to obtain a steady state directly. However, we can use Q t function to give an estimation. We can show that in the bounded confidence case, the steady-state solution is still a "contraction" of initial distribution.
And the degree of offset is depending on the value of (y − x) 2 k(y − x)P(dy)P(dx). We can see that (y − x) 2 k(y − x)P(dy)P(dx) will increase as the radius R of the kernel k(y − x) grows. When R ⟶ ∞ and k(y − x) � 1, the model will degenerate into the linear case in eorem 3.

SDE Related to Initial Values.
Since the process we study is not a Markov process, we need to consider SDE with the form It has no Markov generator, but the non-Markovian of it is "not so bad". e law of X can be solved by the following method. For a given X[0], we can solve the law by e tA P 0 , where A is the generator of the Fokker-Planck equation, and P 0 � δ x 0 is the initial Dirac distribution. en, the solution of the SDE has the form where A(X 0 ) is the generator with given X[0], δ X 0 is the Dirac measure with the support on X 0 , and P 0 is the initial distribution of X 0 . e essence of this method is to regard X[0] as an additional random variable. When SDE has the above form, we need the double layer empirical measure, which will regain the symmetry.
We will consider the SDEs where b i is a stochastic variable such that X i [0] � b i ; then, we will prove the large number law for the double layer empirical measure 1/N i δ X i ,b i , which will give the Mckean-Vlasov limit. is is equal to consider SDEs where X b is the stochastic process with given b, and E(X) � E( dμ(b)X b ).

Intermediate Process.
Andreis et al. [15] proved that the Mckean-Vlasov limit of the equation has the form where P is the law of X. We use the same way as Andreis et al. [15] to reduce the jump-SDE to the "averaging dynamics." eir method is to decompose the jump terms into a martingale and a continuous part. en, the martingale can be easily estimated by Burkholder-Davis-Gundy inequality.
We will show that there is an intermediate process Y i with the following SDE: 8 Complexity such that By comparing these two process directly, we can notice that Using Doob-Meyer decomposition, where N is a martingale. Intuitively, the supremum of a martingale would not be very large. is can be proved by Burkholder-Davis-Gundy inequality and our assumption on h. en, we can obtain the approximation (52). See the second part of the appendix for a detailed proof of (52).

Mckean-Vlasov Process.
e intermediate process has the form like where ω i is a random variable. e asymptotic behavior of the double layer empirical measure, with the case that f and g are bounded, has been studied [38]. ey use Varadhan's lemma to reduce the system to the case without interaction. Since f and g are not bounded, we use the method in the study by Sznitman [28] to compare the two SDEs directly. Using the Lipschitz condition, we can control E[sup r ‖X(r) − X(r)‖] by the first Wasserstein distance ρ. e first Wasserstein distance is defined by where H(μ 1 , μ 2 ) is the set of all probability measures on R d × R d with marginals μ 1 and μ 2 . By Kantorovich-Rubinstein duality theorem, it has a dual representation: where ‖f‖ L is the minimal Lipschitz constant for f. Using this representation, since we have assumed that the coefficients are Lipschitz, we can see that where μ N X is the empirical measure, and μ t is the solution of Mckean-Vlasov equation. For where μ N X is the empirical measure for N independent Mckean-Vlasov processes. en, we can derive our result by considering the Wasserstein distance between the solution of Mckean-Vlasov equation and empirical measure. Using eorem 1 in the study by Zanella and Guillin [41], if the two order moment is finite, when the size tends to infinity and the Wasserstein distance trends to zero, so our claim follows. e detailed proof is given in the appendix.

Conclusion
In this paper, we studied the asymptotic behaviors of the exponential function weight of the self-excitation interaction system X N on the large homogeneous network when N, the number of the individuals, trends to infinity. We proved that there is a Mckean-Vlasov process X, such that lim N⟶+∞ E|X N − X| is zero, whose opinion distribution evolve according to a Mckean-Vlasov type integrodifferential equation, which couples with the initial distribution. e steady state of this equation is also studied. We investigated the steady state distribution of linear and bounded confidence models and showed that the steady distribution of these models is a "contraction" of initial distribution. By virtue of the self-excitation interaction, even if we do not consider bounded confidence, it will give not a perfect consensus, but a distribution depending on the initial distribution, such as this model, can avoid Abelson's diversity puzzle.
and it is shown in the study by Situ [42] that since both b − X and h has linear growth, this SDE has a strong solution, and also, we can estimate the moment.
As for the second condition (17), the above method can be also applicable.

B. Reduce to the Averaging Dynamics
e SDE for stochastic process X i is Following the study by Andreis et al. [15], it is useful to introduce an intermediate process Since the Poisson process is a semimartingale, N(t) can be decomposed into N(t) � N(t) + λt by Doob-Meyer decomposition, where N(t) is a martingale. e λdt in Y is the second part of Doob-Meyer decomposition of N. To simplify the symbol, we set λ � 1.
We will show that for large N, so, X can be approximated by Y. Let In our model (9), g( As in the study by Andreis et al. [15], we can use the Burkholder-Davis-Gundy inequality for martingales. where we use the symmetry, assumption (16), and the Lipschitz condition for h. Combine the two results; then, we get (B.8) Complexity en, using the standard technology of Gronwall's lemma, we get the approximation.

C. Mckean-Vlasov process
As shown in the previous section, we can turn to analyze the stochastic process Y. Let X be a stochastic process with SDE: dX � dP(y)αh(y, X)dt + ω(b − X)dt + σdW, (C. 1) where P is the law of X, and b is a random variable, such that b � X[0].
Using the trick in (47), we can show that P(y) � dμ(b)P b (y), (C.2) where P b is the law of SDE with given b, the initial distribution is δ b , and μ(b) is the given initial distribution of X. We will show that Let Q N be the empirical measures j (1/N)δ X j , we have

(C.4)
Since h is Lipschitz, we can see that Summing over i, So, we only need to estimate the right side. Using the Lipschitz condition for h again, this can be controlled by E[ρ[(Q N X (t), Q X (t))]]where ρ is the Wasserstein distance between the two measures, and Q X is the solution of Mckean-Vlasov equation (72). Q N is the empirical measures j (1/N)δ X j . It follows from [41] eorem 1 that if the moments is bounded (this is satisfied by (17)) us, our claim follows.

Data Availability
is paper is not data-driving, and the model we study is from other data-driving ones. e data used in this article are only the data of the simulation of PDE, which are contained in Figure 1.

Conflicts of Interest
e authors declare that they have no conflicts of interest.