Abstract
An Innocents-Spreaders-Calmness-Removes (ISCR) rumor propagation model is established with nonlinear incidence and time delay on complex networks in this paper. Based on the mean-field theory, the spreading dynamics of the ISCR model are discussed in detail. Firstly, the basic reproduction number is obtained by the next generation matrix method to ensure the existence of rumor-prevailing equilibrium. Secondly, by utilizing the Routh–Hurwitz criterion and LaSalle’s invariance principle, the local stability and global stability of rumor equilibria are proved. Moreover, the optimal control is presented via Pontryagin’s minimum principle, which is to effectively restrain rumor diffusion. Finally, the theoretical results are verified by numerical simulations.
1. Introduction
Rumors are usually defined as unproven words and may damage personal reputation, affect financial markets, cause social panic and instability, and severely disrupt people’s normal and orderly life. Social networks do build a good platform for people’s communication, but it also provides an opportunity for a dishonest person to spread rumors. In the online virtual social platform, everyone has their own online virtual identity, and the virtual identity is the medium for rumors to diffuse [1]. Compared with traditional rumors, Internet rumors spread faster and wider, so people should pay more attention to them and take certain measures to deal with them when necessary. Therefore, it is vital to study the potential mechanism and control measures of rumor propagation on social networks. Based on mathematical models, the research on the mechanism of rumor propagation has received extensive attention from scholars.
Considering the similarities between the spread of rumors and epidemics, traditional rumor propagation models are mostly based on the dynamic models of infectious diseases, such as SI (Susceptible-Infected), SIS (Susceptible-Infected-Susceptible), SIR (Susceptible-Infected-Removed), and so on. Daley and Kendall first studied the rumor propagation dynamic analysis in 1965 and proposed a rumor spreading model named DK model [2]. After that, Maki and Thomson [3] further improved the DK model to the MT model in 1973. In addition, some scholars also applied the fractal-factional order model to the infectious disease model [4], and others discussed the backward bifurcation and optimal control of the infectious disease model [5, 6]. Based on the work of predecessors, many papers which considered other affecting factors of rumor propagation have been proposed afterwards [7–12], such as incubation [7], the proportion of wiseman in the crowd [8], debunking behavior in emergencies [9], media report [10], psychological factors and forgetting mechanism [11], different attitudes towards rumors [12], superspreaders [13], and so on.
As an effective tool, complex networks have laid a fine foundation for studying the spread of rumors on social networks. Zanette first proposed a rumor propagation model on small-world networks by utilizing complex network theory to the study of rumor propagation [14]. Li et al. had established an I2S2R (Ignorants‐Spreaders 1‐Spreaders 2‐ Stiflers 1‐Stiflers 2) rumor spreading model on homogeneous networks [15]. A novel SIR (Susceptible-Propagating-Recovery) rumor spreading model was proposed in both homogeneous and heterogeneous networks by Zhu et al. [16]. Nowadays, there are still many articles considering the rumor spreading model on complex networks [17–22]. Based on the aforementioned models, a rumor propagation model is proposed on complex networks in this paper.
It is worth emphasizing that the incidence rate plays a significant role in the spread of rumors. The incidence mainly includes bilinear incidence and nonlinear incidence. The characteristic of bilinear incidence is that the number of Spreaders increases linearly. However, the psychological changes of Innocents have a lot of influence on rumor propagation, which make bilinear incidence have some limitations, and the nonlinear incidence is used frequently in [23–25]. For example, the incidence was proposed by Capasso and Serio [23] to represent saturation phenomena for large numbers of infectives. The incidence was utilized to explain the following phenomenon: the incidence of infectious diseases may show a downward trend during the peak of infectious diseases because some individuals take protective steps to reduce contact with others individuals [24]. The incidence was used by Ruan and Wang [25] to describe the inhibition rate that came from an increasing number of susceptible individuals. Although the nonlinear incidence was first used in disease transmission models, it can also be considered in the rumor transmission model because of the similar transmission mechanism. Based on the above discussion, the nonlinear incidence is used in this paper to describe the influence of the psychological changes of Innocents on rumor propagation, where represents the spreading ability of rumors and represents the impact of population crowding or changes on Innocents. Obviously, this is more in line with the reality. Hence, considering the nonlinear incidence is of significance for the study of rumor propagation.
As is well known to all, people sometimes may not timely respond to rumors. When an Innocent receives a rumor, it takes some time to consider whether to spread the rumor. This is the reason why time delay exists. Therefore, this will motivate us to study rumor propagation with time delay. For instance, Jain et al. [26] analyzed the effect of delay to influence thinkers. A delayed rumor spreading model was proposed by Li and Ma [27] in emergencies. Chen et al. [28] established multiple delayed models to explore the new characteristics of rumor spreading process. Meanwhile, timeliness is the important characteristic of rumor propagation, so it is vital to study the rumor propagation model with time delay.
In light of above discussion, our main contributions are reflected as follows:(1)A novel ISCR rumor propagation model is proposed with nonlinear incidence and time delay on complex networks. Based on the traditional SIR model, this paper adds a “Calmness” compartment, which makes the model more realistic.(2)In real social networks, the spread of rumors depends on the degree of nodes in the network, so this paper models the spread process of rumors based on the network structure.(3)Because individuals need a certain reaction time after exposure to rumors and the number of Spreaders is limited, this paper considers the effects of time delay and nonlinear incidence on rumor propagation. Moreover, the nonlinear incidence describes the influence of the psychological changes of Innocents on rumor propagation.(4)By utilizing the Routh–Hurwitz criterion and LaSalle’s invariance principle, the local stability and global stability of rumor equilibria are proved in detail.(5)To reduce the density of Spreaders and control costs, an optimal control strategy with time delay is given and analyzed for the controlled system by Pontryagin’s minimum principle. Meanwhile, this paper also analyzes the influence of time delay on optimal control in numerical simulation.
2. Problem Description
In this section, a novel ISCR rumor propagation model is proposed on complex networks to learn the dynamics of rumor propagation mechanism. Four states are proposed to represent the different states of individuals in the process of rumor propagation. Innocents represent those who do not perceive rumors but may be infected. Spreaders represent those who understand rumors and spread them. Calmness represents those who calm down before they stop spreading rumors. Removes represent those who perceive rumors but do not spread them. , , , and denote the density of Innocents, Spreaders, Calmness, and Removes at time , respectively. Moreover, we assume that .
The process of ISCR model is shown in Figure 1. According to the mean-field theory, the model can be described as follows:where denotes the average degree of complex networks, is the average infectious delay of the infectious rumors, represents the recruitment rate of Innocents, denotes the spread rate of Spreaders, is the saturation coefficient that measures the inhibitory or psychological effect of the general public towards rumors, is the transfer rate from Innocents to Removes due to the immune mechanism, is the forgetting rate from Spreaders to Removes by the forgetting mechanism, is the calmness rate of Spreaders, and is the transfer rate from Calmness to Removes. Suppose that every class has the same emigration rate , and the recruitment rate is equal to the emigration rate, that is, . Assume that above parameters are positive.

The initial conditions of system (1) are taken by the following form:where . The Banach space is nonnegative continuous, mapping the interval into .
Remark 1. It is worth noting that a novel state called “Calmness” is introduced in model (1) which means rumor Spreaders may go through a calm period before becoming Removes. In the early stage of an emergency, if the mainstream media has low credibility, then some Spreaders may first tend to believe the information told by acquaintances rather than the mainstream media. However, when acquaintances know that the previous information is fabricated and tell Spreaders, then Spreaders will calm down and may gradually reduce the spread of rumors. Spreaders confirm that the information is indeed fabricated afterwards and then stop spreading rumors.
Lemma 1. The positive invariant set of system (1) is defined bywith initial conditions (2).
Lemma 2. If , , , and , the solutions , , , and of system (1) with the initial conditions (2) are positive for all .
Proof. If , according to the first equation of system (1), one hasIt can be rewritten asThus,Then, according to the variation of constant formula,Hence,Similarly, one can prove that , , and . So, the solutions , , , and of system (1) with the initial condition (2) are positive for all .
3. Dynamic Analysis of the ISCR Model
In this section, the equilibria and the basic reproduction number of system (1) are calculated via utilizing the next generation matrix method [29]. Then, the stability of equilibria is discussed by using the Routh–Hurwitz criterion [30] and LaSalle’s invariance principle [31].
3.1. Equilibria of Model and the Basic Reproduction Number
Let the right side of system (1) be zero, and one has
Let in system (9), and one can easily get the rumor-free equilibrium of system (1) by the following:
To facilitate the calculation of the basic reproduction number of system (9), let ; then, system (9) can be rewritten aswhere
Then,
By calculation, is obtained as follows:
Assume that rumor-prevailing equilibrium is a solution of system (1), that is,
From (15), one can get
Obviously, rumor-prevailing equilibrium of system (1) exists if , and , , , .
Remark 2. According to the expression of , the basic reproduction number is independent of time delay, that is to say, the time delay is only related to the spread time of rumors and will not affect the spread scale of rumors.
3.2. Stability Analysis
Theorem 1. Rumor-free equilibrium is locally asymptotically stable for all if .
Proof. To simplify the calculation, let , , , and , and the linearized system of (1) is as follows:The characteristic equation of system (17) isThen,If , the characteristic equation of system (19) isOne hasIf , in light of the Routh–Hurwitz criterion [30], is locally asymptotically stable for .
Now, if , one only needs to consider the following formula:Assume that (22) has a purely imaginary root , with . Then, separating real and imaginary parts givesEquation (23) is squared and then added to obtainThen,Therefore, equation (22) has no purely imaginary root if . Hence, rumor-free equilibrium is locally asymptotically stable for any if .
Remark 3. If , rumor-free equilibrium is unstable for any .
Remark 4 (see [30]). Consider the following linear delay differential equation:where , , and . The solution of equation (26) is . Also, the characteristic equation of system (26) is .
Theorem 2. Rumor-prevailing equilibrium of system (1) is locally asymptotically stable for all if .
Proof. Let , , , and , and the linearized system of (1) takes the following form:The characteristic equation of system (27) isThen,whereIf , equation (29) can be rewritten asObviously, . Furthermore, one only need to consider the following equation:whereIf , . Then,If , . Then, one has , . Therefore, rumor-prevailing equilibrium is locally asymptotically stable for if .
LetNow, if , assume that (35) has a purely imaginary root , with . Then,Separating real and imaginary parts givesEquation (37) is squared and then added to obtainLet , and equation (38) can be rewritten aswhereLetIf , . Then,LetIf , . Then, (35) has no purely imaginary roots if . Hence, rumor-prevailing equilibrium is locally asymptotically stable for any if .
Theorem 3. If , rumor-free equilibrium is globally asymptotically stable for all .
Proof. Construct the following Lyapunov function:Differentiate along the solution of system (1) as follows:Therefore, ensures that and if and only if , , , . Adopting LaSalle’s invariance principle [31], we have , , and , that is, is globally asymptotically stable for all if .
Remark 5. Rumor-free equilibrium is globally asymptotically stable for all if , that is to say, rumor will gradually die out with the growth of time, and the density of Spreaders will gradually approach zero.
Theorem 4. If , rumor-prevailing equilibrium is globally asymptotically stable for all .
Proof. Define the Lyapunov function aswhereDifferentiate along the solutions of system (1) as follows:Differentiating along the solutions of system (1), one hasDifferentiate along the solutions of system (1) as follows:Hence,LetBecause ,Obviously, and if and only if , , , . Then, one has , , , and , that is, is globally asymptotically stable for all if .
Remark 6. If for all , rumor-prevailing equilibrium is globally asymptotically stable, that means rumors will spread steadily over time and will not die out in the end.
Remark 7. For the construction of the Lyapunov function of the delayed model in this paper, we must first ensure the positive definiteness of the function. Secondly, uncertain limit integral functions are added to the Lyapunov function to offset the previous delayed term. Finally, the algebraic method is used to ensure that the derivative of the Lyapunov function is negative. Meanwhile, Holling-type II functional response is added in this paper, which makes the analysis of the delayed model more difficult.
4. Optimal Control
In this section, optimized control strategies are proposed so as to control the spread of rumors and reduce the control cost of social platforms, and the optimal control solution is found by utilizing Pontryagin’s minimum principle [32]. Next, a control variable is introduced to represent the function of control strategy for . Then, an admissible control set is defined:
The controlled system can be obtained bywith initial conditions (2). Then, an objective function is considered aswhere is a weight coefficient to keep the density of in balance and control the cost of .
Lemma 3 (see [33]). The controlled system (56) with any initial conditions has a unique solution.
Proof. The controlled system (56) can be rewritten aswhereand . System (57) is a nonlinear system with a bounded coefficient. SetA simple calculation shows thatwhere and are some positive constants, independent of state variables , , , and , andHere , , , and , for . Therefore, it is easy to show thatwhere . Thus, it follows that the function is uniformly Lipschitz continuous. The solution of system (55) exists from (62). And the solution of system (55) takes into account the constraints on the controls and the restrictions on the non-negativeness of the state variables.
In order to find an optimal solution, first we find the Lagrangian function and Hamiltonian function for the optimal control system (55). The Lagrangian function of the problem is taken asand we define the Hamiltonian function aswhere , represent the adjoint variables to be determined appropriately.
Lemma 4. There exists an optimal control such thatfor system (55) under the initial conditions (2).
Proof. In fact, the following conditions are satisfied. (i) The set of control and corresponding state variables is not empty. (ii) The control space is convex and closed by definition. (iii) Each right hand side of the state system is continuous and is bounded by a sum of the bounded control and the state. Furthermore, it can be written as a linear function of the control variate with coefficients depending on time and the state. (iv) is convex on the control set and is bounded below. Thus, according to [34], there exists an optimal control . This completes the proof.
Theorem 5. Let be the optimal state solutions with associated optimal control variable for the optimal control system (55). Then, adjoint variables satisfywith transversality conditionsFurthermore, the optimal controlcan be found.
Proof. Differentiating the Hamiltonian function (64) with respect to , , , and and substituting , , , , , and into equations, one hasBy the optimal conditions, one hasConsidering the range of control variable and the property of Hamiltonian function, one hasSo, the optimal control can be obtained as
5. Numerical Simulations
In this section, the theoretical results are verified by some numerical simulations. We choose the average degree in this paper.
5.1. Stability of Rumor-Free Equilibrium
Case 1. In model (1), let , , , , , , , , and . By simple calculation, the basic reproduction number . From Theorem 1, rumor-free equilibrium of model (1) is locally asymptotically stable which is verified by Figure 2(a). Figure 2(b) describes the asymptotic stability of equilibrium .

(a)

(b)
5.2. Stability of Rumor-Prevailing Equilibrium
Case 2. Choose , , , , , , , , and in model (1). By simple calculation, . The local stability of rumor-prevailing equilibrium is examined as shown in Figure 3(a), and the asymptotic stability of rumor-prevailing equilibrium is examined as shown in Figure 3(b).

(a)

(b)
5.3. The Influence of on Rumor Propagation
Choose , and other parameters are fixed as Case 2. The global stability of rumor-prevailing equilibrium can be obtained when , and rumors always exist. From Figure 4(a), we can find that time delay will suppress the peak of Spreaders’ density. As the time delay increases, the maximum density of Spreaders gradually decreases. Therefore, it is appropriate to consider time delay in the spread of rumors.

(a)

(b)
5.4. The Influence of on Rumor Propagation
Let , and other parameters are fixed as Case 2. The influence of on rumor propagation can be observed from Figure 4(b). With the increase of , when becomes larger, the density of Spreaders will be reduced, which indicates that psychological factors have a positive impact on rumor propagation.
5.5. The Influence of on Rumor Propagation
Choose , and other parameters are fixed as Case 2. From Figure 5(a), the influence of is apparent. Infection rate of rumors will affect the ultimate extent of Spreaders, which provides a novel channel to control rumors. For instance, in the early stage of rumor propagation, forbidding Spreaders in online social networks can reduce the infection rate of rumors.

(a)

(b)
5.6. The Influence of on Rumor Propagation
Let , and other parameters are fixed as Case 2. The positive influence of immune mechanism can be shown by Figure 5(b) on rumor propagation. The density of Removes is affected by different immunization rate. Increasing the density of Removes can be achieved by increasing the immunization rate. For example, more popularization of science videos or articles should be published by some well-known online social networking platforms such as Weibo, WeChat, and so on to improve the level of public scientific knowledge.
5.7. Effect of Optimal Control
In this part, the optimal control can effectively lessen the density of Spreaders and extend the region of rumors by numerical simulations. Next, we discuss the influence of time delay on optimal control and give some suggestions to control rumor propagation.
5.7.1. Without Time Delay
Choose in the objective function , and other parameters are fixed as Case 2 besides . The density of individuals with and without optimal control is shown in Figure 6(a). If time delay is equal to zero, the spread speed of rumors is accelerated, and the control time should also be advanced. In daily life, our country and government should spread the harmfulness of rumors through various channels, strengthen publicity, and strengthen the education of the people. When rumors appear, the government and relevant departments should quickly formulate emergency plans to ensure social security. They will do a good job with the goal of “the most true information, the fastest speed, and the best effect” to reduce the harm to the society caused by the spread of rumors.

(a)

(b)
5.7.2. With Time Delay
Choose in the objective function , and other parameters are fixed as Case 2. The density of individuals with and without optimal control is shown in Figure 6(b). It is obvious that control strategy that we proposed controls rumors successfully. At the beginning of the rumor, a downward trend is shown from the density of the Spreaders, and the density of Removes is increasing rapidly under the control strategy.
Next, the path of optimal control and control cost is shown in Figure 7. Choose , and the optimal control is gradually decreased to 0 with time as shown in Figure 7(a). Moreover, the manifestation of objective function is shown in Figure 7(b). Obviously, with the decrease of control force, the control cost gradually increases for a certain time.

(a)

(b)
5.8. General Comparison
Figure 8 clearly shows the impact of time delay on optimal control. Obviously, time delay only affects the control time and does not affect the final result of control. When the time delay is equal to zero, the density of Spreaders decreases rapidly under control, but the density of Spreaders decreases very slowly when time delay exists. Therefore, in the initial stage of rumor propagation, the government and relevant media industries should refute the rumor in time to shorten the public’s response time when facing the rumor, so as to achieve the purpose of controlling the rumor.

6. Conclusion
In this paper, an ISCR rumor propagation model with nonlinear incidence and time delay is presented on complex networks. According to the mean-field theory, the ISCR model is discussed in detail. Firstly, the basic reproduction number is calculated by utilizing the next generation matrix theory. Secondly, the locally asymptotic stability of rumor-free (prevailing) equilibrium is verified by using the Routh‐Hurwitz criterion and the globally asymptotic stability of equilibria is confirmed by using LaSalle's invariance principle under . Because time delay only affects the spread time of rumors, it does not affect the final spread scale and state of rumors, and the following results are given:(i)If , rumor-free equilibrium is locally asymptotically stable under and rumor-prevailing equilibrium under is globally asymptotically stable.(ii)If , the locally and globally asymptotic stability of rumor-free (prevailing) equilibrium is satisfied under .
Then, according to Pontryagin’s minimum principle, the optimal control is presented to minimize the density of Spreaders and control costs. Finally, the theoretical results of this paper are verified via some numerical simulations.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported in part by the National Natural Science Foundation of China under Grant nos. U1703262 and 62163035 and in part by the Tianshan Innovation Team Program under Grant no. 2020D14017 and Graduate Innovation Program of Xinjiang Uygur Autonomous Region under Grant no. XJ2021G019.