Modelling and Analysis of Bacteria Dependent Infectious Diseases with Variable C

J.B.Shukla,Shikha Singh,Jitendra Singh and Sunil Kumar Sharma

1Innovative Internet University for Research,Kanpur,208017,India

2P.P.N.(P.G.)College,CSJM University,Kanpur,208001,India

3College of Computer and Information Sciences,Majmaah University,Majmaah,11952,Saudi Arabia

Abstract:In this research, we proposed a non-linear SIS model to study the effect of variable interaction rates and non-emigrating population of the human habitat on the spread of bacteria-infected diseases.It assumed that the growth of bacteria is logistic with an intrinsic growth rate is a linear function of infectives.In this model, we assume that contact rates between susceptibles and infectives as well as between susceptibles and bacteria depend on the density of the non-emigrating population and the total population of the habitat.The stability theory has been analyzed to analyzed to study the crucial role played by bacteria in the increased spread of an infectious disease.It is shown that as the density of non-emigrating population increases,the spread of an infectious disease increases.It is shown further that as the emigration increases, the spread of the disease decreases in both the cases of contact mentioned above rates, but this spread increases as these contact rates increase.It suggested that the control of bacteria in the human habitat is very useful to decrease the spread of an infectious disease.These results are confirmed by numerical simulation.

Keywords: Mathematical modelling; density dependent contact rates;stability analysis

1 Introduction

Most of the deaths worldwide are caused due to infectious diseases.A severe threat to the wellbeing and public health is caused not only due to the new infectious diseases but also due to the increasing prevalence of drug-related diseases as well as the resurgence of chronic infectious diseases.Recently, considerable evidence found to suggest that common strategies are adopted by different pathogens to cause disease and infection.The infectious organism causes Infectious diseases; these are bacteria, viruses, fungi, etc.Under normal circumstances, disease symptoms may not develop when the immune system of the host is fully functional, but an infectious disease ensues if the immune system of the host is compromised.Bacteria, protozoa, viruses, etc.cause most of the infections in living organisms.Bacteria is a unicellular prokaryotic micro-organism.In the human habitat, due to household emission, various kinds of carriers and vectors grow and survive.They become agents to carry bacteria to food and water of Susceptibles, causing them to be infected by various diseases.Notably, in habitats, which are not clean enough, various types of bacteria such as mycobacterium tuberculosis, vibrio cholera etc.are present, the growth rate of which depends upon the following factors:

• The Bacteria discharges by infective in the environment.

• The characteristics of household discharges such as nutrients, salts, amino acids, and vitamins etc.

• The characteristics of bacteria related to particular diseases such as TB, Typhoid etc.

• The natural conditions such as the climate of the habitat.

Therefore, we assume that the bacteria population density is proportional to the density of carriers, such as house flies, and hence its growth rate is assumed to follow the logistic model [1].These bacteria get transported by carriers to susceptible and their food as well as water, making them infectives indirectly [2].Mathematical models have been used in the study of the spread and control of infectious diseases for a long time.In classical models, the contact rates between susceptible and infectives have been assumed to be constant [3-6], but recently the effect of density-dependent contact rate, death rate etc.on the spread of infectious diseases have been studied.However, the effect of the non-emigrating population of the habitat on the spread of the infectious disease has never been considered in various studies.However, in a realistic situation, the contact rates between susceptible and infectives as well as between susceptible and bacteria depends upon the density of non-emigrating population [7-9].The massive consumption of alcohol affects almost all parts of the body, which in turn responsible for gonorrhea, which is a bacterial disease.This disease infects the parts of the body like urine, eyes, throat, vagina anus, and the female’s fallopian tubes, uterus etc.[10].Recently researchers work in the field of application of mathematical modeling multiscale fast correlation filtering tracking algorithm [11],grammar model [12], grammar mapping and the integer modulo Arithmetic [13] and decision model of knowledge [14] transfer The spread of an infectious disease can be controlled by using awareness programs.However, the disease remains endemic due to immigration [11,12,15,16].

Therefore, we proposed a non-linear mathematical model to study the effects of the following aspects of the spread of bacteria-infected diseases:

• Effect of the non-emigrating population of the habitat.

• Effect of emigration dependent contact rate between susceptible and infective population as well as between susceptible and bacteria population.

• Effects of discharge of bacteria by infectives in the habitat.

• Effect of control of bacteria by using chemicals and pesticides.

2 Mathematical Model with Emigration Dependent Contact Rates

Let the total human population of the habitat at timetbeN(t), which consists of the susceptible population densityX(t)and infective population densityY(t).We assume that the susceptible population of the habitat may be affected by the bacterial populations with densityB(t)(Fig.1).The bacteria population grows logistically with its intrinsic growth rate coefficient as a linear function of infective populationY(t).The contact rateβ(N)between susceptibles and infectives is assumed to be emigration dependent, and the contact rateλ(N)between susceptibles and bacteria depends upon the densities of non-emigrating population and the total populationN.It is assumed as follows,where,N0is the non-emigrating population density of the habitat which is a fraction ofN,βandβ0are constant contact rate coefficients.It is noted from Eq.(1),β(N)increases asN0increases,but it decreases asβ0increases.

Figure 1:Interaction phenomena between susceptibles and infectives

Similarly, the emigration dependent contact between susceptibles and bacteria population is assumed to be a non-negative linear function of the total population as,

where,λandλ0denotes the contact rate coefficients.From Eq.(2) asN0increases,λ(N)increases butλ(N)decreases asλ0increases.The primary purpose of this paper is to study the effect of variable contact rate between susceptibles and infectives, variable contact rate between susceptibles and bacteria, and non-emigrating population on the spread of bacteria-infected disease and the growth of bacteria due to infectives.

In view of the above situation, the diseases dynamics model given below:

where,X(0)=X0>0,Y(0)=Y0≥0,N(0)=N0>0,B(0)=B0≥0 andN=X+Y.

The model system Eq.(3) has the following parameters:

•A:Immigration rate of the human population from outside.

•β:Interaction coefficient due to infective human population density.

•β0:Emigration dependent transmission coefficient due to infectives density.

•λ:Constant transmission coefficient due to bacteria population density.

•λ0:Emigration dependent transmission coefficient due to bacteria population density.

•d:Natural death rate coefficient.

•α:Disease related death rate coefficient.

•ν:Recovery rate coefficient.

•θ:Growth rate coefficient of bacteria population due to natural conducive factors in the habitat.

•θ0:Natural depletion rate coefficient of bacteria population.

•φ:Growth rate coefficient of bacteria by the discharge of the infective population.

•φ0:The control rate of bacteria population by using chemicals and pesticides etc.

The model system Eq.(3) is reduced as follows by usingX=N−Yandθs=θ−φ0>0,we get

With the initial conditions:Y(0)=Y0≥0,N(0)=N0>0,B(0)=B0≥0.

Region of attraction:The following set gives the domain region of the model system Eqs.(4)-(6) as

where,

3 Basic Reproduction Number and Equilibrium Analysis of the Model

The basic reproduction numberR0is used to measure the transmission potential of a disease.It is the average number of secondary infections produced by an infection in a human habitat where the living population is susceptible [17,18].

The non-linear mathematical model Eqs.(4)-(6) has three non-negative equilibria as:

1.E0(0,,0), this is the trivial equilibrium point which exists always.

2.this is the bacteria-free equilibrium point, which exists if the reproduction number

3.E2(Y∗,N∗,B∗), this is the non-trivial equilibrium point, which exists if

It is noted that whenλ=0,λ0=0 thenR1=R0.Proof.The equilibrium pointE0exists obviously.For the model system Eqs.(4)-(6), we prove the existence ofLetY/= 0 andB= 0 thenare the solution, which is obtained from the following two equations:

By using Eqs.(10) and (11) we define a polynomial function inYas,

We note from Eq.(12),

Thus, the equationF(Y)=0 has at least one rootTo show the root is unique,we show thatF′(Y)<0.By differentiating Eq.(12) with respect toY, we get

then, by using Eq.(12) in Eq.(13), we get

HenceF′(Y)<0 forR0>1.Thus a unique root ofF(Y)=0,exists.Now by using the value of, the value ofcan be uniquely determined from Eq.(11).

The other equilibrium pointE2(Y∗,N∗,B∗)is obtained by solving the following equations:

By using the equations Eqs.(16) and (17) in Eq.(15) and setting

We define the following function

It is noted from Eq.(18)

Hence, the equationH(Y)=0 has at least one root inTo show the uniqueness of rootWe need to show thatH′(Y)<0.For this, we differentiate Eq.(18) with respect toY, we get

By using Eq.(18) in Eq.(19), we have

ProvidedR1>1, whereR1is the reproduction number.Hence the equation Eq.(18) have unique root in the interval

4 Stability Analysis

The stability of the equilibrium pointsandE2(Y∗,N∗,B∗)are stated in the following two theorems [19].The equilibrium pointsare unstable andE2(Y∗,N∗,B∗)is locally asymptotically stable provided the following conditions are satisfied

Proof.See Appendix A.

The equilibrium pointE2(Y∗,N∗,B∗)is non-linearly stable in Ω provided the following conditions are satisfied:

where,Bmis defined in Section 2 Eq.(7).

Proof.See Appendix B:It is noted from Eq.(4), the inequalities Eqs.(23) and (24) are stronger than inequalities Eqs.(21) and (22), as expected.Further, the inequalities Eqs.(21) and(23) are satisfied automatically whenα=0.Also, whenφ=0, then inequalities Eqs.(22) and (24)are satisfied automatically.Hence the death rate related coefficient and rates of growth of bacteria caused by infective population density have destabilizing effects on the system.

5 Numerical Simulation of the Model

By using MAPLE, we show the existence and stability of the equilibrium pointE2(Y∗,N∗,B∗)for the following set parameters as given in Tab.1 value of the set of parameters

Table 1:Value of parameters

For these values of parameters, the value of the non-trivial equilibrium pointE2(Y∗,N∗,B∗)corresponding to Eqs.(4)-(6) is obtained as follows:

The Jacobian matrix atE2(Y∗,N∗,B∗)is

The eigenvalues of the above matrix atE2(Y∗,N∗,B∗)are:

Since the sign of one eigenvalue is negative and two of them having negative real part,therefore the non-trivial equilibrium pointE2(Y∗,N∗,B∗)is stable asymptotically.It is noted here that for the assumed set of positive parameters, the conditions Eqs.(21)-(24) are satisfied.By using MAPLE software, the graph ofY vs.Nfor the model Eqs.(4)-(6) is shown in Fig.2,which indicated the non-linear stability of the equilibrium point in theYNplane.

Figure 2:Phase plots between infected human population density Y(t) and human population density N(t)

Figure 3:Effect of transmission coefficient β due to infectives on infective population density Y

Figure 4:Effect of emigration coefficient β0 due to infectives on infective population density Y

Figure 5:Effect of transmission coefficient λ due to bacteria population density on infective population density Y

It depicted in Fig.3, asβincreases, the number of infectives increases.From Fig.4, we can see that as the emigration dependent contact rateβ0between susceptibles and infectives increases, the number of infectives decreases.From Fig.5, as constant contact rateλbetween susceptible and bacterial population increases, the number of infectives increases and from Fig.6,as emigration dependent contact rateλ0between susceptibles and Bacteria population density of habitat increases, then the number of infectives decreases.From Fig.7, it is also observed that as the growth rate coefficient of bacteria due to infective population increases, the infectives increases.It observed from Fig.8, that as the constant depletion rate coefficientφ0increases, the spread of the disease decreases, and from Fig.9, as the natural growth rateθof bacteria increases,the spread of bacteria-infected disease increases.From Fig.10, we observed that immigration increases, then the number of infectives increases, and from Fig.11, as non-emigrating population increases, then the number of infective populations increases.

Figure 6:Effect of emigration coefficient λ0 due to bacteria population density on infective population density Y

Figure 7:Effect of growth rate coefficient φ of bacteria due to infective population density on infective population density Y

Figure 8:Effect of the depletion rate coefficient φ0 of bacteria due to chemicals and pesticides etc.on infective population density Y

Figure 9:Effect of the natural growth rate coefficient of bacteria population θ on infective population density Y

Figure 10:Effect of constant immigration rate A on infective population density Y

Figure 11:Effect of non-emigrating population N0 on infective population density Y

6 Conclusions

In this study, an SIS non-linear model with emigration has been proposed and analyzed to study the effects of the following factors on the spread of bacteria-infected diseases.

• Effect of non-emigrating population.

• Effect of the emigration dependent contact rate between susceptibles and infectives, which is dependent on non-emigrating population and the total human population in the habitat.

• Effect of the emigration dependent contact rate between susceptibles and bacteria population, which depends on non-emigrating population and the total population in the habitat.

• Effect of the discharge rate of bacteria by the infectives in the habitat.

• Effect of the growth rate of bacteria, which is assumed to follow the logistic model, the growth rate of which is a linear function of the infective population.

Using the stability theory, the analysis of the model has shown the following results.

• As the discharge rate of bacteria by infectives increases, the spread of bacteria-infected disease increases.

• As the non-emigrating population density increases, the spread of bacteria-infected diseases increases.

• As the direct contact rate between susceptibles and infectives increases, the spread of a bacteria infected disease increases.

• As the contact rate between susceptibles and bacteria population increases, the spread of bacteria infected disease increases.

• As the emigration rate increases, the spread of bacteria-infected disease decreases.

• As the natural growth rate of bacteria increases, the spread of an infectious disease increases.

• As the immigration rate increases, the spread of bacteria-infected disease increases.

• As the control rate of bacteria in the habitat increases, the spread of the disease decreases.

The simulation study of the non-linear model confirms the above outcomes.It has been concluded that if the bacteria population in the habitat controlled by using pesticides, the spread of bacteria-infected disease can be reduced considerably.

Acknowledgement:The authors would like to express their heartfelt thanks to the editors and anonymous referees for their most valuable comments and constructive suggestions which leads to the significant improvement of the earlier version of the manuscript.The authors are thankful to the Deanship of Scientific Research at Majmaah University.

Funding Statement:Dr.Sunil Kumar Sharma would like to thank Deanship of Scientific Research at Majmaah University for supporting this work under the Project No.R-2021-8.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

Appendix A.Proof of Theorem 4.1

Proof.The model system Eqs.(4)-(6) as:

For the above model system Eqs.(4)-(6), the Jacobian matrix is defined for the system aswhere,

and

The Jacobian matrix for the first equilibrium point

From the above Jacobian matrix, it is clear that one of the eigenvalues isθs>0.Hence the disease and the bacteria-free equilibrium pointis unstable.Now the Jacobian matrix for the second equilibrium pointis

It is clear from the above Jacobian matrix, one eigenvalueis positive.Hence the bacteria-free equilibrium pointis unstable.We study the local stability behavior ofE2(Y∗,N∗,B∗)by direct Lyapunov’s method.Here we use the transformation of variables to linearize the model system as

Y−Y∗=y,N−N∗=n,B−B∗=band by considering a positive definite functionVas:

Differentiating Eq.(27) with respect to t, we get

The linearization of the functions gives,

Then by Eq.(28), we have

Here we choose the constantk1such thatwe get

By combining conditions Eqs.(35), (36) and on substituting the value ofk1in Eq.(34), we get the conditions as stated in Theorem 4.

Appendix B.Proof of Theorem 4.2

Proof.We consider the Lyapunov function,

On differentiating Eq.(39) with respect to ‘t’, we get

By choosingwe can write the equation Eq.(41) as:

and,

By combining the inequalities Eqs.(45)-(46)

Hence the non-trivial equilibrium point is non-linearly stable in Ω provided the following inequalities holds.