Modeling the transmission dynamics of racism propagation with community resilience

Racism spreading can have a vital influence on people’s lives, declining adherence, pretending political views, and recruiters’ socio-economical crisis. Besides, Web 2.0 technologies have democratized the creation and propagation of racist information, which facilitated the rapid spreading of racist messages. In this research work, the impact of community resilience on the spread dynamics of racism was assessed. To investigate the effect of resilience-building, new SERDC mathematical model was formulated and analyzed. The racism spread is under control where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0<1$$\end{document}R0<1 , whereas persist in the community whenever \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0>1$$\end{document}R0>1. Sensitivity analysis of the parameters value of the model are conducted. The rising of transmission and racial extremeness rate provides the prevalence of racism spread. Effective community resilience decline the damages, mitigate, and eradicate racism propagation. Theoretical analysis of the model are backed up by numerical results. Despite the evidence of numerical simulations, reducing the transmission and racial extremeness rate by improving social bonds and solidarity through community resilience could control the spread of racism.

social aspect. Moreover, racism also fragments the groups themselves, allowing individuals quarantined from their peers. Further likely to act in destructive ways either against themselves or against those they see as their enemies. Racism provides public health consequences, and it can be considered a public health crisis [7,22]. Resilience concerns the strength of a system to sustain itself through reform via friendship and casual transformation [1,23].
Community resilience is the continued ability of a community to utilize available resources to respond to, confront, and recover from adverse situations. Moreover, a resilient community would decide what it requires to conquer damage and then sustain to strengthen health, social and economic systems. The community would learn from the past and enhance mobility, primarily through adequate social wealth advancement [8,9]. It requires community resilience building and the brace of identity within convinced narratives and shared responses to threat. A resilience program identifies public affairs, needs systematic and transparent modes of solidarity. Community resilience admits communities to sense granted [3,14]. Resilient communities, consequently, acquire to struggle beside, accommodate, and shape change, and such community resilience is a symbol of social sustainability [2,19,24].
Mathematical modeling has a vital role in system dynamics understanding of realworld phenomena in different aspects. It has been widely applicable in physical, biological as well as social sciences; see reference [5,6,10,13,16,17,20]. An effective mathematical model for spread control of racism will have principal theoretical significance in inhibiting the spreading processes, and then offers practical measures for minimizing losses induced by racism. The study of its spreading follows either model of opinion dynamics or a model of biological epidemics. Recently, [21] studied racism spreading in cyberspace and proposed the SEID model, in which population was separated into four categories: Susceptible (people who are unaware of the racism), Exposed (people who are aware of the racism, yet not decided to be racist or denier), Infected (people who spread the racism), Denier (people who are aware of the racism, but choose not to spread it).
This study is focusing on the assessment of community resilience's effect on the spread of racism. The new model incorporates community resilience as a control strategy to minimize the harm induced by racism. The theoretical analysis, sensitivity analysis, and numerical simulation will be conducted.

Model formulation
In this section, a racism spreading model, which considers the impact of community resilience (anti-racism) campaign, is proposed and referred as SERDC model. Based on the classical SIR model and hesitating mechanism, the total population is divided into four classes: namely S(t)-susceptible population, E(t)-exposed population, R(t) -racist population, and D(t)-denier population. The class S(t) stands for people who do not know the racism. The class E(t) represents people who are aware of the racism, but hesitate to be racist and do not spread it. The class R(t) represents people who are racially racist and actively spreading racism. The class D(t) represents people who are aware of the racism, but choose to deny to spread the racism; individuals who have either stop spreading it or never spread it due to community resilience or anti-racism campaign. Furthermore, let C(t) be the cumulative number of community resilience building campaigns in the region at time t. Each individual must belong to only one of the four classes. Let Λ be recruitment rate of individuals in the susceptible class. Racism increases through the contact of susceptible individuals with racial racists. When an individual in class S comes into contact with an individual in class R get to aware the racism, he or she will switch to E with a probability β , where β is the spreading probability, which depicts the probability that a susceptible individual knows the racism via contact with an individual in class R . If the individuals are well-informed about the threat of racism, then they may have abstained from racist group propagandas. An anti-racism campaign mitigates the propagation of racism by creating awareness about racism threats. We consider the growth rate of the community resilience building is proportional to the number of racist individuals. The constant C 0 is the baseline number of community resilience, which is regularly operating in the community. The resilience campaign is targeted to induce behavioral change in susceptible individuals. Some portions of susceptible class individuals will become racism denier class by the rate θ , where θ is the rate of effective community resilience campaign on susceptible individuals ( β + θ ≤ 1 ). Individuals in class E can switch to class R or D depending on trust and interest about racism. We assume that individuals in class E switch to class R with the probability of φ r when they decide to spread the racism and switch to class D with the probability φ d when they choose to deny racism spreading ( φ r + φ d ≤ 1 ). Racially racist individuals convert to be a racist-denier class due to various reasons with the probability of γ per unit time. Further, it assumed that some resilience community building campaigns fade or lose their impact on people. Moreover, the parameters α and µ 0 are constants denoting the rate of implementation coverage and depletion of community resilience campaigns, respectively. The constant µ denotes the natural death rate.
By considering the above assumptions and the schematic diagram in Fig. 1, the spread dynamics of racism with community resilience is described by using the following system of differential equations: The total population N(t) would be a constant all the time in this model, i.e.,

Mathematical analysis
Well-posedness

From the first equation of the system (1), we obtain
Let H (t) = (βR(t) + θC(t) + µ) and multiply both sides the last inequality of (2) by Integrating this inequality from 0 to t gives Then It can be shown in a similar way that E ≥ 0, R ≥ 0, D ≥ 0 and C ≥ 0 . Therefore, all the solution sets are positive for t ≥ 0 .
Here note that By combining the first four equations of model (1), we have So, the feasible domain of the system (1) is For the well-posedness of the model, we have the following lemma.

Lemma 1
The set Γ is positively invariant to system (1).

Proof
From the equation (3), we have and Now integrating the above equation from 0 to t , we have Therefore, by using the theory of differential inequalities [18], we obtain lim t→∞ sup N (t) ≤ Λ µ , i.e., the human population can not be more than Further, from the last equation of model system (1) and using the fact, R(t) ≤ Λ µ for large t > 0 , we obtain From the theory of differential inequality, we obtain This implies that 0 ≤ C(t) ≤ αΛ µµ 0 for large t ≥ 0 . The proof is complete.

Equilibria and their stability analysis
In this sub-section, we show the feasibility of all equilibria by setting the rate of change with respect to time t of all dynamical variables to zero. Model (1) has two feasible equilibria, which are listed as follows: • The racism-free equilibrium (RFE) of racism prevalence equilibrium point is computed after we have the threshold value R 0 .

Spreading threshold of the model
Here, we will find the basic reproduction number (R 0 ) of model (1) using next-generation matrix approach [12]. We have the matrix of new racist F(X) and the matrix of transfer V(X) . Let X = (E, R, S, D, C) T , model (1) can be rewritten as: where The Jacobian matrices of F(X) and V(X) at The inverse of V is computed as The next-generation matrix K L = FV −1 is given by Therefore, basic reproduction number is R 0 = ρ(K L ) = max(|µ| : µ ∈ ρ(K L )) is spectral radius of matrix K L and it obtained as follows:

Stability of the RFE
Here, we determine the result of linear stability of the equilibrium E 1 of model (1).

Theorem 2
If R 0 < 1 , the racism-free equilibrium E 1 of system (1) is locally asymptotically stable, and it is unstable if R 0 > 1.

and consider a Lyapunov function:
Differentiating L in the solutions of system (1), we get .
Therefore, L < 0 whenever R 0 < 1 . Furthermore, L = 0 if and only if R 0 = 1 . Thus the largest invariant set in W ∈ Γ |L(E, R) = 0 is the singleton. By LaSalle's Invariance Principle the racism-free equilibrium is globally asymptotically stable in Γ , completing the proof.
Theorem 3 completely determines the global dynamics of model (1) in when R 0 ≤ 1 . It establishes the basic reproduction number R 0 as a sharp threshold parameter. Namely, if R 0 < 1 , all solutions in the feasible region converge to the RFE E 1 , and the racism will die out from the community irrespective of the initial conditions. If R 0 > 1 , E 1 is unstable and the system is uniformly persistent, and a racism spread will always exist.

Racism prevalence equilibrium stability
Here, we show the feasibility of racism prevalence equilibrium E 2 . The values of S * , E * , R * , D * , and C * are obtained by solving the following set of algebraic equations: After some algebraic calculations, we get the value of E 2 as: Therefore, there exists a unique positive solution only when R 0 > 1 . Thus, it has a unique racism prevalence equilibrium, E 2 .

Theorem 4
If R 0 > 1 , then the racism prevalence equilibrium point E 2 of system (1) is locally asymptotically stable.

Proof
The Jacobian matrix of the model at E 2 iṡ L ≤ κξ φ r β κξ S * − 1 R = κξ (R 0 − 1)R, since S ≤ S(0) and S ∈ Γ . (4) From the Jacobian matrix (5), we get 1 = −µ . Furthermore, the stability of the system could prove using the Routh-Hurwitz criteria. The characteristic polynomial of (5) is where The polynomial (6) has negative roots (eigenvalues) if all its coefficients terms are positive, or it satisfies Routh-Hurwitz criteria of stability [11]. From (7) we can verify that a 1 > 0, a 4 > 0, a 1 a 2 − a 3 > 0 and a 3 (a 1 a 2 − a 3 ) − a 2 1 a 4 > 0 , when R 0 > 1 . Hence, according to the Routh-Hurwitz criterion, all roots have negative real parts. Thus, the racism equilibrium is asymptotically stable. The proof is complete.
The result tells as, if the initial values of any trajectory are near the equilibrium E 2 , then the solution trajectories approach to the E 2 under the condition R 0 > 1.

Sensitivity analysis
We explore R 0 sensitivity analysis of system (1) to determine the model robustness to parameter values. This is a strategy to identify the most significance parameters of the model dynamics. The normalized sensitivity index Υ is given by: Thus normalized sensitivity indices for parameters are obtained as From the sensitivity indices calculation results, we can identify some parameters that strongly influence the model dynamics. Parameters β and φ r have a positive influence on the basic reproduction number R 0 , that is, an increase in these parameters implies . , an increase in R 0 . While parameters φ d , θ , and γ have a negative influence on the basic reproduction number R 0 , that is, an increase in these parameters implies a decrease in R 0 . Here, we illustrate graphically demonstration of model (1) parameters values sensitivity. Figure 2 illustrates that the rise in contact rate β facilitates the rapid spread of racism. Similarly, Fig. 3, despite racial extremeness rate φ r also increase the propagation dynamics of racism. Further, applying an appropriate intervention measure reduces the contact and racial extremeness rate.
Furthermore, the parameter θ is highly sensitive. The result provides community resilience is an effective control measure to eradicate the racism spread. From

Numerical simulations and analysis
In this section, we carry out numerical simulation of model (1) by using Matlab standard ordinary differential equations (ODEs) solver function ode45.

General dynamics
Here, numerically illustrate the asymptotic behavior of model (1). Consider the initial conditions S(0) = 0.9, E(0) = 0.06, R(0) = 0.04, D(0) = 0, and C(0) = 0.1.   Figure 8 shows the simulation of several solutions of system (1) with different randomly given initial values. Figure 8 shows that all of these solutions eventually converge to the racismfree equilibrium E 1 . Furthermore, social stability is substantially continuous with the strengthening of community cohesion and unity through diversity. Moreover, the damages caused by the existence of racism will be removed. Finally, we have to get a racismfree community by eliminating racist activists from the environment. . We can see that even for a small fraction of the racist state at the beginning the racism persists in the community and stabilizes in time. This means that the trajectories converge to the racism prevalence equilibrium point. Figure 10 shows the simulation of its four solutions with varying initial values. The simulation has shown that all of these solution trajectories finally tend to the racism prevalence equilibrium E 2 . Thus, as established in Theorem 4, the racism persists in the community whenever R 0 > 1 . Further, incompetent community resilience and less recovery of racist individuals lead to the prevalence of racism spread.

Impact of effective contact rate
To evaluate the impact of effective contact in the system dynamics, we carry out some numerical simulations. We set the contact rate β as 0.01, 0.05, 0.15 and β = 0.5 . From Fig. 12, we can see that racist individuals reach a higher peak level as β increases. Figure 12 depicts the influence effective contact rate as shown in the sensitivity analysis. The result implying that containing contact between racist and susceptible individuals leads to reduce the racism spreading. Figure 11 shows the susceptible population also decreases whenever β increases (almost 90% of susceptible population will be racist where β = 0.5).

Impact of racial extremeness rate
Exposed individuals being racist with a rate of φ r , we conduct numerical simulations to shows the effect of φ r on system (1).
We set the rate φ r as 0.0, 0.1, 0.25 and 0.5. In Fig. 14, we can see that the racist individuals increase as φ r increases. This figure shows the notable impact of φ r as shown in Fig. 3, in the sensitivity analysis. The result implies that if the rate φ r of hesitating state individual becomes large, then the racially racist individuals increase and the spread becomes persistent. Similarly, Fig. 13 provides the result susceptible individuals rapidly become the racist group whenever the rate φ r increases.

Impacts of community resilience
To evaluate the effects of community resilience, we carry out simulations by varying effect resilience campaign rate θ. Figure 15 shows the role of resilience campaign on the racist class. We set θ = 0, 0.01, 0.05 and θ = 0.2 , it can observed that R solution trajectories converges to zero when θ values are increases. Thus, to intensify the coverage of community resilience, provide conquering the size of racist individuals. From the simulations, we learn efficient community resilience is a relevant control measure against racism propagation. This simulation result shows that more than 10% effective resilience-building is required to conquer the propagation of racism.

Discussion and conclusions
In this paper, we investigated the dynamics of racism spreading with control measures. A new SERDC racism spreading model with community resilience has been exhibited and simulated. We have studied the existence and stability of the equilibrium points of the offered model. Using the next-generation matrix approach, we have derived an explicit formula for the basic reproduction number, R 0 . Furthermore, we have shown that the RFE is globally asymptotically stable if R 0 ≤ 1 and unstable whenever R 0 > 1 . For the case where R 0 > 1 , we have proven that there exists a unique racism prevalence equilibrium, which is locally asymptotically stable. The sensitivity analysis exhibits the significant control parameter of the model. The key parameters are the contact rate ( β ), the racial extremeness rate ( φ r ), recovery rate ( γ ), and effective community resilience campaign ( θ ). Some numerical simulations were conducted to verify the theoretical analysis. Further, we demonstrate the impact of the model parameter value on the system dynamic. The general dynamics simulation result implies the racism spread is under control where R 0 < 1 , but persists in the population whenever R 0 > 1.
Moreover, a resilient community would determine what it demands to conquer damage and then sustain to improve health, social, and economic systems. Similarly, multicultural communities are strengthening cohesion through diversity, and they enhance socio-economical benefits equally. Truly robust community resilience can promote the power of integrity and minimize damages by removing racism from the living environment.