Applications of the nonlinear Klein/Sinh-Gordon equations in modern physics: a numerical study

1.   Introduction

Malaria is one of the world's most significant infectious diseases ^[1]. Malaria is a life-threatening disease caused by parasites that is usually transmitted to persons through the bites of female Anopheles mosquitoes ^[2]. Malaria gives rise to great pressure for the global prevention and control of infectious diseases ^[3]. World Health Organization reported ^[2] that there were an estimated 247 million malaria cases, including 619,000 deaths worldwide in 2021, and the majority of cases and deaths occurred in sub-Saharan Africa. The African region accounted for a disproportionate share of the global malaria burdens ^[4,5]. In 2021, the African region was home to 95% of global malaria cases and 96% of global malaria deaths, and children under 5 years old accounted for about 80% of all malaria deaths there ^[2]. There are 5 kinds of parasite species that cause malaria in humans, and two of these species P. falciparum and P. vivax pose the greatest threat ^[6]. The first malaria symptoms such as headache, fever and chills usually appear 10–15 days after the bite of a malaria mosquito and may be mild and difficult to be recognized as malaria, which implies that malaria exists the incubation period ^[2]. It was reported that asymptomatic infections were more prevalent in sub-Saharan Africa, where an estimated 24 million people had asymptomatic malaria infections ^[7]. Thus asymptomatic infections can occur during malaria transmission.

Asymptomatic infected people have no clinical symptoms, but they are contagious and the impact of asymptomatic infections on malaria transmission is enormous ^[8,9]. Bousema et al. ^[10] pointed out that asymptomatic carriers contributed to sustained transmission of malaria in local populations, and there was substantial evidence that an increase in the number of asymptomatic carriers at specific time intervals affected the dynamics of malaria transmission. Laishram et al. ^[4] concluded that asymptomatic malaria infections was a challenge for malaria control programs.

Since the emergence of malaria, scholars at home and abroad have been studying the pathogenesis and transmission dynamics of malaria. In all research methods, mathematical modeling is undoubtedly one of the most intuitive and effective methods. Many researchers have studied the dynamic evolution of malaria transmission by applying some mathematical models of malaria. In 1911, Ross ^[11] put forward a basic ordinary differential equations (ODEs) malaria model. Afterwards, MacDonald ^[12] extended Ross's model, and gave first the definition of the basic reproduction number. The extended Ross's model was said to be the Ross-Macdonald model. Subsequently, the Ross-Macdonald model has been extended to higher dimensions and more factors affecting malaria transmission have been taken into account (see, e.g., ^[1,6,13]). For example, Kingsolver ^[14] extended the Ross-Macdonald model and explained the greater attraction of infectious humans to mosquitoes in 1987. Safan and Ghazi ^[1] developed a 4D ODEs malaria transmission model with standard incidence rates, and analyzed the dynamic properties of equilibria of the malaria model. In 2020, Aguilar and Gutierrez ^[6] established a high-dimensional ODEs malaria model with asymptomatic carriers and standard incidence rate, and dealt with local dynamics of the disease-free equilibrium of the malaria model.

Over the years, considering the incubation period of malaria, lots of researchers established some time-delayed malaria models (see, e.g., ^{[3,5,15,16,17,18]}). For instance, in 2008, Ruan et al. ^[5] first established a class of Ross-Macdonald model with two time delays, and investigated the stability of equilibria of the model and the impact of time delays on the basic reproduction number. In 2019, Ding et al. ^[3] proposed a malaria model with time delay, and investigated the global stability of the uninfected equilibrium of the model as well as its uniform persistence. For the moment, there are few theoretical analysis of the model of malaria with standard incidence rate. Recently, Guo et al. ^[13] established a malaria transmission model with time delay and standard incidence rate, and they studied the global dynamic properties of equilibria of the model. Based on this, we extend and improve the model in ^[13], namely, we establish a malaria transmission model with asymptomatic infections, standard incidence rate and time delay, and then study the global dynamic properties of equilibria of the malaria model.

The remainder of this paper is organized as follows. In Section 2, we put forward a time-delayed dynamic model of malaria with asymptomatic infections and standard incidence rate, and prove the well-posedness as well as dissipativeness of the system. In Section 3, we obtain the existence conditions of malaria-free and malaria-infected equilibria of the system, and verify the local dynamic properties of equilibria in terms of the basic reproduction number $R_{0}$. In Section 4, to obtain the global dynamic property of the malaria-infected equilibrium for $R_{0} > 1$, we acquire the weak persistence of the system through some analysis techniques. In Section 5, by utilizing the Lyapunov functional method and the limiting system of the model combining stability of partial variables, we obtain the global stability results of malaria-free and malaria-infected equilibria in terms of $R_{0}$, respectively.

2.   Model formulation

In order to delve into the details of malaria transmission, we develop a time-delayed model with asymptomatic infections and standard incidence rate. The population is classified into four compartments, which are denoted by ${S}_{h}$: susceptible individuals, ${A}_{h}$: asymptomatic infected individuals, ${I}_{h}$: symptomatic infected individuals, ${R}_{h}$: recovered individuals, respectively. The mosquitoes are classified into two compartments, which are denoted by ${S}_{m}$: susceptible mosquitoes and ${I}_{m}$: infected mosquitoes, respectively. Then the model of malaria transmission is proposed as follows:

where ${N}_{m}(t) = {S}_{m}(t)+{I}_{m}(t).$ Here, time delay $\tau\geq0$, and all other parameters of system (2.1) are assumed to be positive and $p\in(0, 1)$. The description of parameters are listed in Table 1.

The phase space of system (2.1) is

where ${C}$ is the Banach space of continuous functions mapping from $[-\tau, 0]$ to $\mathbb{R}_{+}^{6}$ with $\mathbb{R}_{+} = [0, \infty)$ and the supremum norm. In the following, the well-posedness as well as dissipativeness of system (2.1) will be investigated in ${C}_{+}$.

Theorem 2.1. The solution $u(t) = {({S}_{m}(t), {I}_{m}(t), {S}_{h}(t), {A} _{h}(t), {I}_{h}(t), {R}_{h}(t))}^{T}$ of system (2.1) with any $\phi \in{C}_{+}$ exists uniquely, and is non-negative and ultimately bounded on $\mathbb{R}_{+}$. In particular, ${({S}_{m}(t), {S}_{h}(t))}^{T}\gg{\bf{0}}$ on $(0, \infty)$, and ${C_{+}}$ is positively invariant for system (2.1).

Proof. In view of the basic theory of delay differential equations (DDEs) ^[19,20], the solution $u(t)$ of system (2.1) with any $\phi\in C_{+}$ is unique on its maximum interval $[0, T_{\phi})$ of existence. Firstly, we will prove that the solution $u(t)$ is non-negative on $[0, T_{\phi})$. According to the continuous dependence of solutions of DDEs on parameters ^[19,20], then for any $b\in(0, T_{\phi})$ and a sufficiently small $\varepsilon > 0,$ the solution $u(t, \varepsilon) = (u_{1}(t, \varepsilon), u_{2}(t, \varepsilon), u_{3} (t, \varepsilon), u_{4}(t, \varepsilon), u_{5}(t, \varepsilon), u_{6}(t, \varepsilon))^{T}$ through $\phi$ of the following model:

uniformly exists on $[0, b]$. Consequently, we claim $u(t, \varepsilon)\gg{\bf{0}}$ on $[0, b)$. It is clear that $\dot{u}_{i}(0, \varepsilon) > 0, i\in I_{6} = \{1, 2, 3, 4, 5, 6\}$ whenever $u_{i}(0, \varepsilon) = 0.$ Next, we prove the claim by contradiction. Suppose that there exists $\bar{t}\in(0, b)$ such that $u_{i}(\bar{t}, \varepsilon) = 0$ for some $i\in I_{6}$ and $u(t, \varepsilon)\gg{\bf{0}}$ for $t\in(0, \bar{t})$, where

As a result, it holds

Since

it follows from (2.2) that $\dot{u}_{i}(\bar{t}, \varepsilon) > 0$, which yields a contradiction to (2.3). Thus, we have $u(t, \varepsilon)\gg{\bf{0}}$ for $t\in(0, b)$.

Letting $\varepsilon\rightarrow0^{+}$ gives that $u(t, 0) = u(t)\geq0$ for any $t\in\lbrack0, b)$. Note that $b\in(0, T_{\phi})$ is chosen arbitrarily, so that $u(t)\geq0$ on $[0, T_{\phi})$. It is obvious that $T_{\phi} > \tau$. Therefore, from system (2.1), we have that for any $t\geq\tau$,

As a consequence, by the comparison principle, we can obtain that $u(t)$ is bounded. Accordingly, from the continuation theorem of solutions of DDEs ^[19], it follows $T_{\phi} = \infty$. Consequently, we have

Therefore, the solution $u(t)$ with any $\phi\in{C}_{+}$ uniquely exists, and is non-negative and ultimately bounded on $\mathbb{R}_{+}$. Moreover, it is not difficult to get that ${({S}_{m}(t), {S}_{h}(t))}^{T}\gg{\bf{0}}$ on $\mathbb{R}_{+}\backslash\{0\}$, and ${C_{+}}$ is a positive invariant set for system (2.1).

3.   Local stability

To begin with, it follows easily the malaria-free equilibrium $E_{0} = ({S} _{m}^{0}, 0, {S}_{h}^{0}, 0, 0, 0)^{T}$, where ${S}_{m}^{0} = {\lambda}_{m}/{\mu} _{m}$ and ${S}_{h}^{0} = {\lambda}_{h}/{\mu}_{h}$. By using the similar method in ^[21,22], we can calculate the basic reproduction number

To get a malaria-infected equilibrium (i.e., positive equilibrium) $E^{\ast } = \left({S}_{m}^{\ast}, {I}_{m}^{\ast}, {S}_{h}^{\ast}, {A}_{h}^{\ast}, {I} _{h}^{\ast}, {R}_{h}^{\ast}\right) ^{T}$, we have the following lemma.

Lemma 3.1. System (2.1) exists a unique $E^{\ast}\gg{\bf{0}}$ when and only when ${R}_{0} > 1$.

Proof. First of all, the malaria-infected equilibrium equations can be obtained as follows:

Note that ${S}_{m}^{\ast}+{I}_{m}^{\ast} = \lambda_{m}/\mu_{m}$, it follows from (3.1) that

Substituting the third and the fourth equations in (3.2) and ${S} _{m}^{\ast} = \lambda_{m}/\mu_{m}-{I}_{m}^{\ast}$ into the first equation in (3.2), there holds

In consequence, (3.3) possesses a unique positive root

if and only if ${R}_{0} > 1.$ Thus, we can conclude that $E^{\ast}$ is a unique malaria-infected equilibrium of system (2.1) if and only if $R_{0} > 1$, where $E^{\ast}$ satisfies

Next, by adopting similar techniques in ^[23,24,25], we will discuss the local dynamic properties of the malaria-free equilibrium $E_{0}$ and the the malaria-infected equilibrium $E^{\ast}$ with respect to $R_{0}$. First of all, for the local stability of the equilibrium $E_{0}$, we have the theorem as follows.

Theorem 3.1. For any $\tau\geq0$, the malaria-free equilibrium $E_{0}$ is locally asymptotically stable (LAS) when $R_{0} < 1$, and unstable when $R_{0} > 1$.

Proof. With some calculations, the characteristic equation of the linear system of system (2.1) at $E_{0}$ can be obtained as follows:

where

Clearly, Eq (3.5) possesses three negative real roots: $-\mu_{h}$ (double) and $-\mu_{m}$. The other roots of Eq (3.5) satisty $\mathcal{H}(\lambda) = 0$. Next, we will prove that any root $\lambda$ of $\mathcal{H}(\lambda) = 0$ has negative real part. Suppose, by contradiction, $\lambda$ has the nonegative real part. Then it follows from $\mathcal{H}(\lambda) = 0$ that

Taking the modulus of both sides in (3.7), we have

and

for $R_{0} < 1$ and $\tau\geq0,$ which leads to a contradiction. Therefore, the real part of each root of the Eq (3.5) is negative. Accordingly, $E_{0}$ is LAS for $R_{0} < 1$ and $\tau\geq0.$

Now, we prove that the $E_{0}$ is unstable for $R_{0} > 1$ and $\tau\geq0$ by the zero theorem. Clearly, for $R_{0} > 1$ and $\tau\geq0$, we can get

According to the zero theorem, there must exsit a positive real root in Eq (3.6). Thus, $E_{0}$ is unstable for $R_{0} > 1$ and $\tau\geq0$.

For the local stability of the equilibrium $E^{\ast}$, we can obtain the theorem as follows.

Theorem 3.2. For any $\tau\geq0$, the malaria-infected equilibrium $E^{\ast}$ is LAS if and only if $R_{0} > 1$.

Proof. By Lemma 3.1, we just require to demonstrate the sufficiency. Let

With direct calculation, the characteristic equation of the linear system of system (2.1) at $E^{\ast}$ can be got as follows:

where

Clearly, Eq (3.8) has two negative roots: $-\mu_{h}$ and $-\mu_{m}$. The other roots of Eq (3.8) satisfy $g(\lambda) = 0$. Then, we will prove that any root $\lambda$ of $g(\lambda) = 0$ has negative real part by contradiction. Assume that $\lambda$ has the non-negative real part. By $g(\lambda) = 0$, we can get

Taking the modulus of both sides in (3.9), for $R_{0} > 1$ any $\tau \geq0,$ it follows

and

Obviously, this is a contradiction. Hence, the real part of each root of the Eq (3.8) is negative for $R_{0} > 1$ and $\tau\geq0$, which ensures the local stability of the equilibrium $E^{\ast}$.

4.   Weak persistence

Generally, to obtain the global stability of the equilibrium $E^{\ast}$, we need to prove the strong persistence or uniform persistence of system (2.1). However, we study the weak persistence of system (2.1), which can ensure the global stability of the equilibrium $E^{\ast}.$ Of course, the weak persistence of system (2.1) is more accessible than its strong or uniform persistence. Now, we define

and let

be the solution of system (2.1) with any $\phi\in\digamma.$ It follows easily that $\digamma$ is a positive invariant set of system (2.1), and $u(t)\gg{\bf{0}}$ for $t > 0$. Hence, we discuss the weak persistence of system (2.1) in $\digamma$.

According to ^[26], system (2.1) is said to be weakly persistent if

We define $u_{t} = \left(S_{mt}, I_{mt}, S_{ht}, A_{ht}, I_{ht}, R_{ht}\right) ^{T}\in C_{+}$ to be $u_{t} (\theta) = u(t+\theta)$, $\theta\in\lbrack-\tau, 0]$ for $t\geq0$, and $u_{t}$ is the solution of system (2.1) with $\phi$. Inspired by the work in ^[13], we study the weak persistence of system (2.1). First, we have the following lemma.

Lemma 4.1. Assume that ${R}_{0} > 1$, $\theta\in(0, 1)$ and $\limsup _{t\rightarrow \infty}I_{m}(t)\leq\theta I_{m}^{\ast}$. Then

Proof. It follows from (2.4) that

For any $\epsilon > 1$, there can be found ${\varrho} = {\varrho}(\phi, \epsilon)\geq0$ such that

and then

Consequently,

Letting $\epsilon\rightarrow1^{+}$, it holds

The malaria-infected equilibrium equations imply that

Theorem 4.1. Let $R_{0} > 1$. Then $\limsup_{t\rightarrow \infty}I_{m}(t)\geq I_{m}^{\ast}$.

Proof. We will use the proof by contradiction to verify this result. Provided that $\limsup _{t\rightarrow \infty}I_{m}(t) < I_{m}^{\ast}$. Whereupon, one can find a $\theta \in(0, 1)$ such that $\limsup_{t\rightarrow \infty}I_{m}(t)\leq\theta I_{m}^{\ast}.$ Using Lemma 4.1, we can get that there is an $\epsilon_{0} > 0$ such that for any $\epsilon\in(0, \epsilon_{0}),$

Thanks to Lemma 4.1, it follows that for any $\epsilon \in(0, \epsilon_{0})$, there can be found $\mathscr{T}\equiv \mathscr{T}(\epsilon, \phi) > 0$ such that

Now, we define the functional on $\digamma$ as follows,

Obviously, $L(u_{t})$ is bounded since $L$ is continous on $\digamma$. Then for $t\geq\mathscr{T}$, the derivative of $L$ along the solution $u_{t}$ is given by

Denote

Next, we will prove that $I_{m}(t)\geq c$ for $t\geq\mathscr{T}$. If not, there exists a $\mathscr{T}_{0}\geq0$ such that $I_{m}(t)\geq c$ for $t\in\lbrack\mathscr{T}, \mathscr{T}+\tau+\mathscr{T}_{0}]$, $I_{m} (\mathscr{T}+\tau+\mathscr{T}_{0}) = c$ and $\dot{I}_{m}(\mathscr{T}+\tau +\mathscr{T}_{0})\leq0$. Then it follows that for $t\in\lbrack\mathscr{T}, \mathscr{T}+\tau +\mathscr{T}_{0}]$,

It easily follows from Eq (4.2) that for $t\in\lbrack \mathscr{T}, \mathscr{T}+\tau+\mathscr{T}_{0}]$,

Analogously, one can get

for $t\in\lbrack\mathscr{T}, \mathscr{T}+\tau+\mathscr{T}_{0}]$. By $R_{0} ^{2} = S_{h}^{0}S_{m}^{0}/S_{h}^{\ast}S_{m}^{\ast}$ and (4.1), we have

Accordingly, we conclude that

Clearly, this contradicts $\dot{I}_{m}(\mathscr{T}+\tau+\mathscr{T}_{0})\leq 0$. As a result, $I_{m}(\mathscr{T})\geq c$ for $t\geq\mathscr{T}$. Hence, for $t\geq\tau,$

which hints $L(u_{t})\rightarrow \infty$ as $t\rightarrow \infty$. Accordingly, this contradicts the boundedness of $L(u_{t})$.

According to Theorem 4.1, we have the following result.

Corollary 4.1. If $R_{0} > 1$, then for any $\tau\geq0$, system (2.1) is weakly persistent.

5.   Global stability

We will study the global asymptotic stability of the equilibria $E_{0}$ and $E^{\ast}$ with respect to $R_{0}$. For this purpose, we get from (2.4) the following limiting system of system (2.1):

Adopting a similar argument as in the proof of Theorem 2.1, it follows that the solution

of system (5.1) through any $\varphi = (\varphi_{1}, \varphi_{2}, \varphi_{3}, \varphi_{4}, \varphi_{5}, \varphi_{6})^{T}\in C_{+}$ uniquely exists, and is non-negative and ultimately bounded on $[0, \infty)$. Setting

gives that $z_{t} = (S_{mt}, I_{mt}, S_{ht}, A_{ht}, I_{ht}, R_{ht})^{T}\in C_{+}$ is also the solution of system (5.1) through $\varphi$ for $t\geq0$. We can find easily that $E^{\ast}$ and $E_{0}$ are also the equilibria of system (5.1), and $C_{+}$ is a positive invariant set of system (5.1). By the way, ${(S_{m}(t), S_{h}(t)})^{T}\gg{\bf{0}}$ for $t > 0$. Define $H(v) = v-1-\ln v, \; v > 0.$ Thereupon, for the global dynamic property of the equilibrium $E_{0}$ of system (2.1), we have the theorem as follows.

Theorem 5.1. For any $\tau\geq0$, the malaria-free equilibrium $E_{0}$ is GAS when $R_{0} < 1$ and GA when $R_{0} = 1$ in $C_{+}$.

Proof. By Theorem 3.1, it follows that for $R_{0} < 1$, $E_{0}$ is LAS. Thus, we only need to prove that for $R_{0}\leq1$, $E_{0}$ is $GA$. Let $u_{t}$ be the solution of system (2.1) with any $\phi\in C_{+}$ and $z_{t}$ be the solution of system (5.1) though any $\varphi\in C_{+}$. Let $\omega(\phi)$ be the $\omega$-limit set of $\phi$ with respect to system (2.1). In order to prove the global attractivity of $E_{0}$, we just need to show that $\omega (\phi) = \left\{ E_{0}\right\}.$ By Theorem 2.1, we know that $u_{t}$ is bounded on $C_{+}$. Hence, it follows from (2.4) that $\omega (\phi)$ is a compact set, and is also a subset of $C_{+}$.

Let us define the following functional ${\mathcal{V}}$ on $L_{1} = \left\{ \varphi\in {C}_{+}:\varphi_{1}(0) > 0, \varphi_{3}(0) > 0\right\} \subseteq C_{+}$

where

Obviously, ${\mathcal{V}}_{1}$ is continuous on $L_{1}$. Since $z_{t}\in L_{1}$ for $t\geq1,$ the derivative of ${\mathcal{V}}$ along $z_{t}$ ($t\geq1$) is given by

where $h = \beta_{h}\mu_{m}/\lambda_{m}$. Considering (5.2) and (5.3), we can conclude that both $S_{h}(t)$ and $S_{m}(t)$ are persistent. In other words, there exists a $\sigma = \sigma(\varphi) > 0$ such that $\liminf_{t\rightarrow \infty}S_{h}(t) > \sigma$ and $\liminf_{t\rightarrow \infty}S_{m}(t) > \sigma$. As a result, $\omega(\varphi)\subseteq L_{1}$, where $\omega(\varphi)$ is the $\omega$-limit set of $\varphi$ with respect to system (5.1). It is evident that ${\mathcal{V}}$ is a Lyapunov functional on $\left\{ z_{t}:t\geq1\right\} \subseteq L_{1}$. Then it follows from [27,Corollry 2.1] that $\dot{\mathcal{V}}(\psi) = 0$, $\forall\psi\in\omega(\varphi)$.

Assume that $z_{t}$ is the solution of system (5.1) through any $\psi\in \omega(\varphi)$. Then the invariance of $\omega(\varphi)$ gives that $z_{t}\in\omega(\varphi)$ for $t\in\mathbb{R}$. According to (5.3), we have ${S}_{m}(t) = {S}_{m}^{0}$ and ${S}_{h}(t) = {S}_{h}^{0}$ for $t\in \mathbb{R}$. From system (5.1) and the invariance of $\omega(\varphi),$ it follows that $I_{m}(t) = A_{h}(t) = I_{h}(t) = R_{h}(t) = 0$ for $t\in\mathbb{R}$. Thus for $R_{0}\leq1$, it holds that $\omega(\varphi) = \left\{ E_{0}\right\}$, which implies that $W^{s}(E_{0}) = C_{+}$, where $W^{s}(E_{0})$ is the stable set of $E_{0}$ with respect to system (5.1).

Now, we prove that the equilibrium $E_{0}$ of system (5.1) is uniformly stable for ${R}_{0}\leq1$ by using the similar approach in ^[28,29]. Observe that the first five equations of system (5.1) can constitute an independent subsystem

It is not difficult to find that

is a positive invariant set with respect to system (5.4). Clearly, system (5.4) has a malaria-free equilibrium $X^{0} = \left(S_{m}^{0}, 0, S_{h}^{0}, 0, 0\right) ^{T}$. According to (5.2), (5.3), [27,Corollary 3.3] hints that $X^{0}$ is uniformly stable. Define

By the definition of uniform stability of $X^{0}$, it follows that for any $\epsilon > 0$, there is $\delta\leq2\epsilon/3$ such that for any $\xi \in\mathcal{C}_{+}$ and $\Vert\xi-X^{0}\Vert < \delta$, there holds

where $X_{t}$ is the solution of system (5.4) with $\xi$. Then considering the sixth equation of system (5.1), we can get

Consequently, for any $\phi\in C_{+}$ and $\Vert\phi-E_{0}\Vert < \delta$, it follows that for any $t\geq0$,

and then

Thus, the equilibrium $E_{0}$ is uniformly stable for system (5.1).

Next, we claim that $\omega(\phi) = \left\{ E_{0}\right\}$ for ${R}_{0} \leq1.$ We first have $E_{0}\in\omega(\phi)$ since $\omega(\phi)\subseteq C_{+} = W^{s}(E_{0}).$ Assume that there exists $\psi\in\omega(\phi)$ such that $\psi\neq E_{0}.$ Let $\alpha(\psi)$ be the $\alpha$-limit set of $\psi$ for system (5.1). Then it follows from the invariance and the compactness of $\omega(\phi)$ that $\alpha(\psi)\subseteq\omega(\phi)$. The invariance of $\alpha(\psi)$ and the stable set $C_{+}$ of $E_{0}$ yield that $E_{0} \in\alpha(\psi)$. Obviously, this contradicts to the stability of $E_{0}$ for system (5.1). Therefore, $\omega(\phi) = \left\{ E_{0}\right\}.$

Remark 5.1. In fact, the stability of the malaria-free equilibrium ${E}_{0}$ of system (5.1) can be acquired for $R_{0} = 1$ in the proof of Theorem 5.1. But using the proof of Theorem 3.1, we can not obtain the stability of the equilibrium ${E}_{0}$ for $R_{0} = 1$.

For the global dynamic property of the equilibrium $E^{\ast}$ of system (2.1), we can draw the following theorem.

Theorem 5.2. For any $\tau\geq0$, the malaria-infected equilibrium $E^{\ast}$ is globally asymptotically stable if and only if $R_{0} > 1$ in $\digamma$.

Proof. From Lemma 3.1 and Theorem 3.2, we just require to prove that $E^{\ast}$ is GA for $R_{0} > 1$. Let $u_{t}$ be the solution of system (2.1) with any $\phi\in\digamma$ and $z_{t}$ be the solution of system (5.1) through any $\varphi\in\digamma$. We can obtain that $\digamma$ is positively invariant for system (5.1), and $z(t)\gg{\bf{0}}$ for $t\geq0$. In order to show that $E^{\ast}$ is GA, we only need to show that $\omega(\phi) = \{E^{\ast}\}$. It follows from Theorem 2.1 that $u_{t}$ is bounded on $\digamma$. Thus, it holds that $\omega(\phi)$ is compact.

Let us define a functional ${\mathcal{V}}$ on $L_{2} = \left\{ \varphi\in C_{+}:\varphi_{i}(0) > 0, i = 1, 2, 3, 4, 5\right\} \subseteq\digamma$ as follows

where

Clearly, ${\mathcal{V}}_{2}$ is continuous on $L_{2}$. In as much as $z_{t}\in L_{2}$ for $t\geq\tau+1$, the derivative of the functional ${\mathcal{V}}$ along $z_{t}$ in $t\geq\tau+1$ is given by

where $h = \beta_{h}\mu_{m}/\lambda_{m}$. Further, it follows from the equilibrium equations that

By (5.6) and (5.7), it follows that $\omega(\varphi)\subseteq L_{2}$. It is clear that ${\mathcal{V}}$ is a Lyapunov functional on $\left\{z_{t}: t\geq\tau+1\right\} \subseteq L_{2}$. As a consequence, [27,Corollary2.1] implies that $\dot{\mathcal{V}}(\psi) = 0$ for any $\psi\in\omega(\varphi)$.

Let $z_{t}$ be the solution of system (5.1) for any $\psi\in \omega(\varphi)$. Then the invariance of $\omega(\varphi)$ indicates that $z_{t}\in\omega(\varphi)$ for any $t\in\mathbb{R}$. Thus, from (5.7), it follows that for any $t\in\mathbb{R}$,

By (5.8) and the third equation of system (5.1), we get that for any $t\in\mathbb{R}, {I}_{m}(t) = {I}_{m}^{\ast},$ ${A}_{h}(t) = {A}_{h}^{\ast}$ and ${I}_{h}(t) = {I}_{h}^{\ast}.$ Consequently, by the invariance of $\omega(\varphi)$ and system (5.1), it holds that ${R}_{h}(t) = {R}_{h}^{\ast}$ for any $t\in\mathbb{R}$. Therefore, it follows that $z_{0} = \psi = E^{\ast},$ and then $\omega(\varphi) = \left\{ E^{\ast }\right\}$, which implies $W^{s}(E^{\ast}) = \digamma$, where $W^{s}(E^{\ast })$ is the stable set of $E^{\ast}$ with respect to system (5.1).

Now, we prove that the equilibrium $E^{\ast}$ of system (5.1) is uniformly stable by using the similar argument in ^[28,29]. Note that system (5.4) has a unique malaria-infected equilibrium $X^{\ast} = \left(S_{m}^{\ast}, I_{m}^{\ast}, S_{h}^{\ast}, A_{h}^{\ast}, I_{h}^{\ast}\right)$. It follows from (5.6), (5.7) and [27,Corollary 3.3] that $X^{\ast}$ is uniformly stable. By the definition of uniform stability of $X^{\ast}$, it follows that for any $\epsilon > 0$, there is $\delta\leq2\epsilon/3$ such that for any $\xi \in\mathcal{C}_{+}$ and $\Vert\xi-X^{\ast}\Vert < \delta$, there holds

where $X_{t}$ is the solution of system (5.4) through $\xi$. Hence, for any $\phi\in L_{2}$ and $\Vert\phi-E^{\ast}\Vert < \delta$, it follow from (5.5) that

for any $t\geq0$, where

is used. Thus, we have

which gives that the equilibrium $E^{\ast}$ is uniformly stable with respect to system (5.1).

Next, we claim that $\omega(\phi) = \{E^{\ast}\}$. From Theorem 4.1, it follows that $\omega(\phi)\cap\digamma\neq\emptyset,$ which gives that $E^{\ast}\in\omega(\phi)$. Assume that there is $\psi\in\omega(\phi)$ such that $\psi\neq E^{\ast}.$ Then the invariance and the compactness of $\omega(\phi)$ implies that $\alpha(\psi)\subseteq\omega(\phi)$. By Theorem 3.2, we have that $E^{\ast}\in\alpha(\psi)$. Obviously, this contradicts to the stability of $E^{\ast}$ with respect to system (5.1). Therefore, $\omega(\phi) = \left\{ E^{\ast}\right\}.$

Remark 5.2. Indeed, the proof of Theorem 5.2 can be simplified, i.e., the stability of system (5.1) is not a must, because $E^{\ast}\in\omega(\phi)$ and Theorem 3.2 can imply that $\omega(\phi) = \{E^{\ast}\}.$ But if we use [30,Theorem 4.1] to prove that $\omega(\phi) = \{E^{\ast}\},$ then the stability of system (5.1) is required.

Acknowledgements

This work is partially supported by the National NSF of China (No. 11901027), the Major Program of the National NSF of China (No. 12090014), the State Key Program of the National NSF of China (No. 12031020), the China Postdoctoral Science Foundation (No. 2021M703426), the Pyramid Talent Training Project of BUCEA (No. JDYC20200327), and the BUCEA Post Graduate Innovation Project (No. PG2022143).

Conflict of interest

The authors declare there is no conflict of interest.