Global Dynamics of Modeling Flocculation of Microorganism

From a biological perspective, a dynamic model describing the cultivation and flocculation of a microorganism that uses two different kinds of nutrients (carbon source and nitrogen source) is proposed. For the proposed model, there always exists a boundary equilibrium, i.e., Rhodopseudomonas palustris -free equilibrium. Furthermore, under additional conditions, the model also has five positive equilibria at most, i.e., the equilibria for which carbon source, nitrogen source, Rhodopseudomonas palustris and flocculants are coexistent. The phenomena of backward and forward bifurcations are extensively discussed by using center manifold theory. The global stability of the boundary equilibrium of the proposed model is deeply investigated. Moreover, the local stability of the positive equilibrium and the uniform persistence of the proposed model are discussed. Under additional conditions, the global stability of the positive equilibrium is studied. Some control strategies are given by the theoretical analysis. Finally, some numerical simulations are performed to confirm the correctness of the theoretical results.


Introduction
Photosynthetic bacteria, which are common microorganisms in the natural environment, have been applied in the field of environmental protection, such as in the treatment of sewage, domestic wastewater and the bioremediation of sediment mud polluted with organic matter (see, for example, [1][2][3][4]).On the other hand, photosynthetic bacteria can produce relatively large amounts of physiologically-active substances, such as vitamin B 12 , ubiquinone (coenzyme Q10 ), 5-aminolevulinic acid (ALA) and RNA (see, for example [5]).In particular, vitamin B 12 has been used in treating anemia and as an eye lotion.Recently, applications as health food supplements have received considerable attention.Coenzyme Q10 has been used in treating heart diseases for many years.Further, coenzyme Q10 has been used not only as a medicine, but also as some food supplements, because of its physiological activities.One of the developments of ALA applications is in the area of photodynamic diagnosis.RNA is an attractive source of 5 -ribonucleotides for use as a flavor enhancer in the food industry.In recent years, the production of RNA has been used as a dietary source of pyrimidine for human immune functions (see, for example [6,7]).
Some photosynthetic bacteria, such as Rhodopseudomonas palustris, are extensively used in the production of lycopene, aquaculture, and so on [8].It can use sunlight, inorganic and organic compounds for energy.Further, Rhodopseudomonas palustris can have practical value for removing microcystin from the water body during algal blooms [9].It can also degrade 2,4,6-trinitrotoluene (TNT), which has negative effects on the human body and aquatic life, resulting in a major threat to drinking and irrigation water supplies, as well as the recreational use of surface waters worldwide.Moreover, Rhodopseudomonas palustris is regarded as the most promising microbial system for the biological production of hydrogen, which has been extensively developed because of its high-energy content and clean product after combustion [10].However, the concentration of Rhodopseudomonas palustris is very low under anaerobic light culture conditions (see, for example, [11,12]).Therefore, cost-efficient harvesting of Rhodopseudomonas palustris is a new challenge.In order to harvest Rhodopseudomonas palustris from the liquid, it is necessary to flocculate the single cells into large cell aggregates.Flocculation is a chemically-based separation process that requires less energy than centrifugation and ultrafiltration and, thus, is regarded as the most promising means for degrading microorganisms.Since algal toxins of blooms have happened occasionally in recent years, the problems of degrading microorganisms have received wide attentions (see, for example, [13][14][15][16]).
Flocculants are a kind of important water treatment reagent, which can be divided into organic flocculants and inorganic flocculants according to the chemical compositions [17].Although organic flocculants, such as polyacrylamide, are frequently used in wastewater treatment and industrial downstream processes because of their high efficiency, some of them are not easily degraded in nature [18,19], and some of the monomers derived from synthetic polymers are harmful to the human body (see, for example, [20,21]).To solve these environmental problems, inorganic flocculants are increasingly being seen as an alternative in the settlement of microorganisms, more specifically in wastewater treatment owing to their inexpensive and nontoxic characteristics.Thus, inorganic flocculants may be used as nontoxic, cost-effective and widely-available flocculants for harvesting Rhodopseudomonas palustris (see, for example, [22][23][24][25][26][27] and the references therein).
Mathematical models have played an important role in better understanding microbiology and population biology (see, for example, [28,29]).In recent years, the dynamics of the chemostat models has received considerable attentions (see, for example, [29]).The article of Smith and Waltman has played an important role in the development of the chemostat models [30].From then on, much research on the chemostat models has been extensively studied by many authors.A model describing two populations of microorganisms competing for one single limiting nutrients was proposed in [31].Later, the model was extended to an arbitrary number of populations in [32,33].These models, which were studied in articles [31][32][33], have proven that they all include a competitive exclusion effect.In the articles [34][35][36][37][38][39], some further developments have been performed on the chemostat models to place the relevant models in a naturally more sensible manner.
Let S(t) and X(t) denote the concentration of the nutrients and Rhodopseudomonas palustris, respectively, in the culture vessel at time t.P(t) denotes the concentration of inorganic flocculants, which are used for harvesting Rhodopseudomonas palustris (see Figure 1).The constants S 0 > 0 and P 0 > 0 denote the input concentration of the nutrients and flocculants, respectively.For simplicity, it is assumed that the input of the nutrients and flocculants is continuous.The constant D > 0 is the dilution rate of the chemostat.The constant τ ≥ 0 denotes the time delay involved in the conversion of nutrients to Rhodopseudomonas palustris.The flocculation rate of microorganisms is assumed to be a bilinear mass-action function response m 1 X(t)P(t), where m 1 ≥ 0 is the per capita contact rate.At the same time, flocculants produce loss or consumption [17], and the loss rate of flocculants is also assumed to be a bilinear mass-action function response m 2 X(t)P(t), where m 2 ≥ 0 is constant.Thus, in [40], the following dynamic model has been proposed: In Model (1), r ≥ 0 and r 1 ≥ 0 are constants, and the bilinear mass-action uptake function S(t)X(t) has been used.
It should be mentioned here that, the analysis reveals that Model (1) proposed in [40] exhibits the phenomenon of backward bifurcation for the existence of positive equilibria.Moreover, the local stability properties of the equilibria have been dealt with in detail.
People found that the influence of different nutrients has played an important role in the culture of microorganisms.In order to take this into consideration, appropriate combinations of nutrients are considered in chemostat models.Models with two competitors and two perfectly-complemented growth-limiting nutrients are studied in [41,42].Local asymptotic conditions for the equilibria are derived.
When there is a microorganism to compete for two or more resources, it may become necessary to consider how the resources, once consumed, interact to promote growth.In [41], the authors employ consumer needs to provide a criterion to classify resources and classify resources as perfectly complementary, perfectly substitutable or imperfectly substitutable.Perfectly-complementary resources are different essential substances that must be taken together.In this case, each substance fulfills different functions with respect to the growth of microorganisms.For example, carbon source and nitrogen source may be complementary for the growth of bacterium.
Motivated by the papers mentioned above, in this paper, we further consider a dynamic model describing the cultivation and flocculation of Rhodopseudomonas palustris, and the nutrients presented in [40] will be divided into carbon source and nitrogen source, which are perfectly complementary in the culture of Rhodopseudomonas palustris (see, for example, [5,[41][42][43][44]).We assume that the growth of Rhodopseudomonas palustris is always co-limited by carbon and nitrogen for all possible nutrient conditions.We do not consider the case where only one of these nutrients limits growth, for example, in environmental scenarios of high carbon, but very low nitrogen loads, the growth of bacteria may be purely nitrogen limited.
Let C(t) and N(t) denote the concentration of carbon source and nitrogen source, respectively, in the culture vessel at time t.The constants C 0 > 0 and N 0 > 0 denote the input concentration of carbon source and nitrogen source, respectively.We assume that the conversion of nutrients to microorganism biomass occurs instantly.That is, the time delay τ in Model (1) equals zero.Hence, we have the following dynamic model describing the cultivation and flocculation of a microorganism: In Model (2), the parameters D, r, m 1 , m 2 and P 0 are the same as Model (1).The term rµ 1 (C(t))µ 2 (N(t)) is the growth rate of Rhodopseudomonas palustris, and the terms r 1 µ 1 (C(t))µ 2 (N(t)) and r 2 µ 1 (C(t))µ 2 (N(t)) represent the quantity of the decreasing of the carbon source and nitrogen source, respectively, where r 1 and r 2 are non-negative constants; the functions µ 1 (C(t)) and µ 2 (N(t)) are nonnegative and continuous for C(t) ≥ 0, N(t) ≥ 0. For the simplicity of the theoretical analysis, in this paper, the functions µ 1 (C(t)) and µ 2 (N(t)) are chosen as Monod-type functions, i.e., , where K 1 > 0 and K 2 > 0 are the half-saturation constants with respect to the carbon source and nitrogen source, respectively.δ i (i = 1, 2) are yield coefficients (see, for example, [29]), which are defined as: Therefore, the dynamic Model (2) can be rewritten in the following form: It is convenient to introduce dimensionless variables.In particular, we define: and still denote C, N, X, P, K 1 , K 2 , t, r 1 , r 2 , r, m 1 and m 2 with C, N, X, P, K 1 , K 2 , t, r 1 , r 2 , r, m 1 and m 2 , then Model (3) becomes: , According to the biological considerations, the initial condition of Model (4) is given as: where the constants C 0 , N 0 , X 0 and P 0 represent the initial concentration of the carbon source, nitrogen source, Rhodopseudomonas palustris and flocculants respectively.The purpose of this paper is to tackle the existence of backward and forward bifurcations by using center manifold theory and to investigate the global stability properties of the two classes of equilibria by constructing the suitable Lyapunov functions.
The organization of the paper is as follows.The global existence, nonnegativity and boundedness of the solutions of Model (4) are investigated in Section 2. In Section 3, the existence of the equilibria and the phenomena of backward and forward bifurcations are extensively discussed.In Section 4, the global stability of the boundary equilibrium of Model ( 4) is discussed by the stability theory of ordinary differential equations.Furthermore, we consider the local stability of positive equilibrium, the uniform persistence of Model (4) and the global asymptotic stability of the positive equilibrium in Section 5.In Section 6, some control strategies are given by the theoretical analysis.Some discussions are given in Section 7.

The Global Existence, Nonnegativity and Boundedness of Solutions
From the biological considerations, it is necessary to show that all of the solutions of Model (4) are nonnegative and bounded for all t ≥ 0. By using the basic theory of ordinary differential equations [45] and some simple calculations, it is not difficult to show the following result.
Theorem 1.The solution (C(t), N(t), X(t), P(t)) of Model (4) with the initial Condition ( 5) is existent, unique and nonnegative for all t ≥ 0 and satisfies: Proof.From the theory of the local existence of solutions for ordinary differential equations, it can be obtained that the solution (C(t), N(t), X(t), P(t)) of Model ( 4) is existent and unique for t ∈ [0, δ).Here, δ is some positive constant [29,45,46].Furthermore, we also have that the solution (C(t), N(t), X(t), P(t)) is nonnegative for t ∈ [0, δ).
We can easily show that C(t), N(t) and P(t) are bounded on t ∈ [0, δ).Let us further show that X(t) is also bounded on t ∈ [0, δ).For t ≥ 0, define: From Model (4), we obtain that, for t ≥ 0, From ( 6) and the well-known comparison principle, we have that M(t) is also bounded for t ∈ [0, δ).Hence, by employing the continuation theorems of the solutions [29,45,46], the solution (C(t), N(t), X(t), P(t)) is existent and unique for any t ≥ 0. Similarly, the solution is nonnegative for any t ≥ 0.
Thus, from the comparison principle, we have that: By the first equation of Model (4), we obtain that, for t ≥ 0, Again, we can conclude that lim inf t→∞ C(t) ≥ C. By using the technique similar above, we can show that lim inf t→∞ N(t) ≥ N, lim inf t→∞ P(t) ≥ P. We complete the proof of Theorem 1.
Theorem 2. The compact set: attracts all of the solutions of Model ( 4) and is positively invariant with respect to Model (4).
Proof.According to Theorem 1, it only needs to be proven that Ω is positively invariant with respect to Model (4).That is, it needs to be shown that In fact, if there exists some . By the Lagrange mean-value theorem, there exists some On the other hand, from ( 6), we have that Ṁ(t 2 ) ≤ α − M(t 2 ) < α − α = 0, which is a contradiction.Thus, M(t) ≤ α for any t ≥ 0. Therefore, from Model (4), we have that for any In fact, if there exists some . By the Lagrange mean-value theorem, there exists some On the other hand, from Model (4), we obtain that: which is a contradiction.Thus, C(t) ≥ C for any t ≥ 0. Therefore, from Model (4), we obtain that, for any t ≥ 0, from which we easily obtain that N ≤ N(t), P ≤ P(t), for any t ≥ 0. This completes the proof of Theorem 2.

The Existence of the Equilibria and Its Classification
Let (C, N, X, P) be any equilibrium of Model (4).Then, (C, N, X, P) satisfies the following nonlinear algebraic equations, Model ( 4) always has the boundary equilibrium E 0 (1, 1, 0, 1).The existence of E 0 indicates that, if there is no Rhodopseudomonas palustris to be added into the culture vessel at the beginning of the culture, the concentrations of the carbon source, nitrogen source and flocculants always maintain the constant values 1, 1 and 1, respectively.The equilibrium E 0 (1, 1, 0, 1) is also called Rhodopseudomonas palustris-free equilibrium.
Define the basic bifurcation parameter as: .
Let (C * , N * , X * , P * ) be any positive equilibrium of Model (4).From ( 7), we have that: Clearly, X * should satisfy the following conditions: Substituting the second, third and forth equations of ( 8) into the first equation gives a fifth order algebraic equation, where: Let us consider the necessary condition for the existence of the positive equilibria of Model (4).
From the first equation in (8), we obtain the following function: which implies that r > 1 is a necessary condition for a positive equilibrium to exist.Using the methods similar to [47], we can give the sufficient conditions of the existence of the positive equilibria of Model (4).The following results (Theorem 3) follow from the various possibilities enumerated in Table A1  Hence, under suitable conditions, there may be at most five different positive roots for the fifth order algebraic equation.Let X = X * be any such positive root, which also satisfies Conditions (9).Thus, from (8), C = C * > 0, N = N * > 0 and P = P * > 0 can be obtained.Therefore, Model (4) at most has five positive equilibria of the type of E * (C * , N * , X * , P * ).The equilibrium E * (C * , N * , X * , P * ) indicates that the carbon source, nitrogen source, Rhodopseudomonas palustris and flocculants may be coexistent for any time t ≥ 0.
Remark 1.The existence of multiple positive equilibria of Model (4) when R 0 < 1 (shown in Table A1; see Appendix A) indicates the possibility of the existence of backward bifurcation (see, for example, [40,47,48]), where the stable boundary equilibrium co-exists with a stable positive equilibrium.This can be explored below via numerical simulations (see Figure 2a,b).A rigorous result can be obtained using center manifold theory [49].The detailed proof is given in Appendix B.  From the analysis given in Appendix B, we have established the following result.
The existence of backward bifurcation implies that stable boundary equilibrium and stable positive equilibrium may be coexistent.In biology, this means that the basic bifurcation parameter R 0 is not the threshold value, which is used to determine whether Rhodopseudomonas palustris can be harvested successfully or not.In this case, more complicated dynamic properties may occur.There may exist a new threshold value, which is less than R 0 and used to determine whether Rhodopseudomonas palustris can be harvested successfully or not.

The Global Stability of the Boundary Equilibrium
Global stability properties of the equilibria E 0 or E * imply that the asymptotic properties of the carbon source, nitrogen source, Rhodopseudomonas palustris and flocculants in the culture vessel are not dependent on the initial values C 0 , N 0 , X 0 and P 0 .For the global stability property of the boundary equilibrium E 0 of Model (4), we have the following result.Theorem 5.If R 0 < 1, then the boundary equilibrium E 0 (1, 1, 0, 1) of Model ( 4) is locally asymptotically stable.Further, if: then the boundary equilibrium E 0 of Model ( 4) is globally asymptotically stable.
Proof.It is easy to show that the boundary equilibrium E 0 (1, 1, 0, 1) of Model ( 4) is locally asymptotically stable by the characteristic equation of the linearization of Model (4).Next, we prove the global asymptotic stability of the boundary equilibrium E 0 of Model (4).
If for some t, C( t) = 1 and N( t) = 1, then we have that Ċ( t) = Ṅ( t) = 0. Hence, from the first or second equation of Model (4), we obtain that X( t) = 0. Thus, for any t ∈ R, we have that X(t) ≡ 0.
By the analysis of Theorem 5, we obtain that the boundary equilibrium E 0 of Model ( 4) is globally asymptotically stable (see Figure 3a).

Stability and Uniform Persistence
In this section, we give the local asymptotic stability of the positive equilibrium and uniform persistence of Model (4).Further, under additional conditions, the global asymptotic stability of the positive equilibrium is obtained by using some techniques of constructing Lyapunov functions.

The Local Asymptotic Stability of the Positive Equilibrium
For the local asymptotic stability of the positive equilibrium of Model (4), we have the following result.Theorem 6.If the positive equilibrium E * (C * , N * , X * , P * ) of Model (4) exists and the condition: holds, then the positive equilibrium E * (C * , N * , X * , P * ) of Model ( 4) is locally asymptotically stable.
Proof.For convenience, we assume that: For the linearized system, the corresponding characteristic equation of Model (4) can be expressed as follows: where: According to the Routh-Hurwitz criterion, we need to show that: For the sake of simplicity, we define: By computation, it can be obtained that: , we obtain that ∆ 1 > 0, ∆ 2 > 0, ∆ 3 > 0, and ∆ 4 > 0. Hence, from the Routh-Hurwitz criterion, we obtain that the positive equilibrium E * (C * , N * , X * , P * ) of Model ( 4) is locally asymptotically stable.This completes the proof of Theorem 6.

Uniform Persistence
As pointed out in [29,45,50], uniform persistence is an important concept in the cultivation of microorganisms.From a biological perspective, one basic question about biological models involves the long-time survival of the species.To this end, we may want to know whether the constructed model is uniformly persistent with respect to one or more species.That is, the number of the species will keep positive and bounded away from zero for any positive time.In recent years, this kind of research has been particularly common in microorganisms (see, for example, [45]), especially in the simulation of harmful algae growth, where we usually want to predict if the harmful algae will go towards extinction or if blooms of the harmful algae will remain, which directly impacts human health and food webs in aquatic ecosystems.From a mathematical point of view, uniform persistence can give sufficient criteria for the existence of a positive equilibrium for the dissipative system [51].
Model ( 4) is said to be uniformly persistent if there are positive constants n 1 , n 2 , n 3 , n 4 , N 1 , N 2 , N 3 , N 4 , such that each positive solution (C(t), N(t), X(t), P(t)) of Model (4) satisfies: Now, we give a result on the uniform persistence of Model (4).To proceed, we introduce the following notation and terminology.Denote by R(t) (t ≥ 0) the family of solution operators corresponding to Model (4).The ω-limit set ω(x) of x consists of y , such that there exists a sequence From Theorem 2, we obtain that G is positively invariant corresponding to R(t).It is easy to see that G 0 is also positively invariant.

Proof.
By Theorem 1, it can be obtained that nonnegative solutions of Model (4) are point dissipative.According to the definitions above, it suffices to show that ∂G 0 repels uniformly nonnegative solutions of Model (4).It is obvious that there is only one equilibrium E 0 in M ∂ .
In the following, let us show that E 0 is factually isolated.That is, there exists some neighborhood U of E 0 , such that E 0 is the largest invariant set in U.
In fact, for sufficiently small positive constant ε, let us choose: We show that E 0 is the largest invariant set of U for some ε.

The Global Asymptotic Stability of the Positive Equilibrium
In this subsection, the global asymptotic stability of the positive equilibrium E * (C * , N * , X * , P * ) of Model ( 4) is studied.Theorem 9. Assume that R 0 > 1 and Model (4) has a unique positive equilibrium E * (C * , N * , X * , P * ).Let: If D 1 ≥ 0, D 2 ≥ 0, D 3 < 0 and D 4 > 0 hold, then the positive equilibrium E * (C * , N * , X * , P * ) of Model ( 4) is globally asymptotically stable.Here: Proof.Let us consider the following Lyapunov function on A 0 , Clearly, V 2 (C * , N * , X * , P * ) = 0, and V 2 is positive definite with respect to E * (C * , N * , X * , P * ).The derivative of V 2 along the solutions of Model ( 4) is: for any t ≥ 0. Noting that: Thus, we have that: for any t ≥ 0. For any t ≥ 0, we have that: .
Hence, it can be obtained that, for any t ≥ 0, Motivated by the papers mentioned in [50,54], the derivative of V 2 along the solutions of Model ( 4) is represented in quadratic form.Then, we obtain that: for any t ≥ 0.
If the condition D 1 ≥ 0 holds, we obtain that, for any t ≥ 0, Furthermore, from the condition D 2 ≥ 0 and C(t) ≤ 1 for any t ≥ 0, we obtain that: These conditions ensure that the coefficients of (C(t) − C * )(X(t) − X * ) and (N(t) − N * )(X(t) − X * ) are nonnegative.Thus, we obtain that, for any t ≥ 0, where W is a symmetric 4 × 4 matrix and W = {b ij } 1≤i,j≤4 with: The other parameters have the same definitions as in Theorem 9.Then, W is negative definite.By Lemma 6.2 provided in [50], we can obtain that the real quadratic form Y T (t)WY(t) is negative definite.Then, by using some classical analysis techniques of differential equations, the positive equilibrium E * (C * , N * , X * , P * ) of Model ( 4) is globally asymptotically stable.Next, we illustrate that the conditions given in Theorem 9 are reasonable.We set r = 5.8, K 1 = 1, K 2 = 0.6, m 1 = 0.0001, m 2 = 1, r 1 = 6, r 2 = 10.It is easy to see that there exists a unique positive equilibrium E * (C * , N * , X * , P * ) ≈ (0.6709, 0.4515, 0.3181, 0.7587), and the conditions D 1 -D 4 hold.By Theorem 9, the positive equilibrium of Model ( 4) is globally asymptotically stable (see Figure 3b).
In general, it is very difficult to obtain the global stability properties of the positive equilibrium of Model (4).In this paper, we have obtained the sufficient conditions to ensure the global stability properties of the positive equilibrium by constructing a suitable Lyapunov function.Because of the complexity of the conditions, it is difficult to grasp biological intuition.However, from the numerical simulations, we have found some interesting biological phenomena.The conditions D 1 -D 4 can be satisfied if the flocculation rate of microorganisms (m 1 ) is small enough (see, Figure 3b).From a biological point of view, these conditions are reasonable.For more detailed biological considerations, we will leave it for further investigation.

Control Strategies
In this section, some control strategies are provided by suitable theoretical analysis.If R 0 < 1 holds, then, we have from Theorem 5 that the boundary equilibrium E 0 of Model ( 4) is locally asymptotically stable.We note that the condition R 0 < 1 is equivalent to the following inequality, All of the parameters r, K 1 , K 2 and m 1 in (12) are defined as in Model ( 4). ( 12) can be further written as the following form, r Here, all of the parameters r, K 1 , K 2 , m 1 , C 0 , N 0 , P 0 and D in ( 13) are defined as in Model (3).
In view of the biological meanings of the parameters in Model (3) and Condition (13), Theorem 5 indicates that the concentration of Rhodopseudomonas palustris in the chemostat tends to zero, and the concentration of the carbon source, nitrogen source and flocculants may tend to the constant values C 0 , N 0 and P 0 , respectively, as time t increases, if one of the following two cases occurs: (a) reducing the absorption of Rhodopseudomonas palustris or the carbon input concentration, or the nitrogen input concentration; (b) improving the velocity or the flocculation effect or flocculant input concentration.These cases are reasonable, since they imply the insufficient sources for Rhodopseudomonas palustris to grow.Hence, in the environmental science field, it can be used to remove algae and heavy metals.
From Theorem 7, we have that Model (3) is uniformly persistent if R 0 > 1.This means that the concentrations of the carbon source, nitrogen source, Rhodopseudomonas palustris and flocculants in the chemostat may be ultimately maintained at some positive constant values, as time t increases, if one of the following two cases occurs: (a) improving the absorption of Rhodopseudomonas palustris, or the carbon input concentration, or the nitrogen input concentration; (b) reducing the velocity or the flocculation effect or flocculant input concentration.
These control strategies can be performed by numerical simulations.
In the following, for convenience, we simulate the extinction or persistence of microorganism (X(t)) numerically by using (12) and Model (4).
If the parameters are chosen as in Table A2 (see Appendix A), Rhodopseudomonas palustris in the chemostat will tend to extinction (see Figure 4a).
If the absorption of Rhodopseudomonas palustris (r) is improved from r = 5.8 to r = 8 and the other parameters are the same as Table A2, Rhodopseudomonas palustris in the chemostat will tend to be constant (see Figure 4b).
If the flocculation effect (m 1 ) is reduced from m 1 = 1 to m 1 = 0.1, Rhodopseudomonas palustris in the chemostat will tend to be constant (see Figure 5a).A1.(b) Rhodopseudomonas palustris in the chemostat will tend to be constant if the absorption of Rhodopseudomonas palustris (r) is improved.The initial conditions are If Michaelis-Menten constant of carbon (K 1 ) is reduced from K 1 = 1 to K 1 = 0.5, Rhodopseudomonas palustris in the chemostat will tend to be constant (see Figure 5b).
Rhodopseudomonas palustris in the chemostat will tend to be constant (see Figure 6).

Discussion and Conclusions
In the paper, based on some biological considerations and chemostat models, a dynamic model governed by ordinary differential equations with four variables (carbon source, nitrogen source, Rhodopseudomonas palustris and flocculants) is presented.There is a boundary equilibrium and at most five positive equilibria for the proposed model.To give a theoretical analysis for the existence of all of the positive equilibria of Model (4), the method of the Descartes rule of signs is applied to the classifications of the positive roots of a fifth order algebraic equation.
The local and global stability properties of the boundary equilibrium of Model (4) have been studied in detail.An interesting phenomenon of backward and forward bifurcations is observed.That is, there may exist two positive equilibria even if the condition R 0 < 1 holds.Hence, sufficient Condition (10) to ensure the global stability of the boundary equilibrium is reasonable in mathematics.
The local stability of the positive equilibrium of Model ( 4) is also carried out.From Condition (11), we have that the positive equilibrium is locally asymptotically stable when the flocculation coefficient m 1 is small enough.Hence, Condition (11) is also reasonable in biology.
Uniform persistence of Model (4) has also been completely studied under the condition R 0 > 1. Uniform persistence has very important significance both in mathematics and biology, and it characterizes the long-term survival of some microorganisms [45].
Finally, some control strategies are provided by simple theoretical analysis.From Theorem 5, we have that Rhodopseudomonas palustris in the chemostat will tend to extinction if R 0 < 1.
In this case, these control strategies can be applied to remove Cyanobacteria, which are well known to produce a variety of toxins and have serious harm on human health.From Theorem 7, we have that Rhodopseudomonas palustris in the chemostat will tend to be positive constant if R 0 > 1.In this case, these control strategies can be widely used for the collection of useful microorganisms.
It is well-known that the existence of time delays is inevitable in biology.For example, in the cultivation of microorganisms, there are always time delays in the process of transferring nutrients and the uptake of nutrients.Hence, chemostat models with time delays that account for the time lapsing between the uptake of nutrients by cells and the incorporation of these nutrients as biomass have been given much attention [40,[55][56][57].Based on Model (3), it may have the following more general form with time delays, In Model ( 14), the constants ρ 1 ≥ 0 and ρ 2 ≥ 0 are the rate constants at which the carbon source and nitrogen source are recycled because of the death of Rhodopseudomonas palustris.The constant σ ≥ 0 is a fixed time during which the carbon source and nitrogen source are released completely from dead Rhodopseudomonas palustris.The constant τ ≥ 0 denotes the time delay involved in the conversion of nutrients to Rhodopseudomonas palustris.The factor e −d 1 τ is the probability constant at which Rhodopseudomonas palustris remains in the culture vessel during the conversion process.The theoretical analysis of Model (14) will be studied separately.

Appendix B
In Section 3, we have discussed the phenomena of backward and forward bifurcations by numerical simulations.In this section, the center manifold theory is used on Model (4) to obtain the rigorous result (see, for example, [40,47,48]).
Let C = x 1 , N = x 2 , X = x 3 , P = x 4 , so that Model (4) can be re-written in the following form: (B1) The Jacobian matrix of Model ( 15) at E 0 (1, 1, 0, 1) is given by: Suppose r is chosen as a bifurcation parameter.Solving R 0 = 1 gives: Eigenvectors of J(E 0 )| r=r * It can be shown that the Jacobian matrix of Model ( 15) at r = r * has a right eigenvector (corresponding to the zero eigenvalue) given by ω = (ω 1 , ω 2 , ω 3 , ω 4 ) T , where: Further, the Jacobian matrix of Model ( 15) at r = r * has a left eigenvector (associated with the zero eigenvalue) given by v = (v 1 , v 2 , v 3 , v 4 ) T , where: Computations of â and b For Model (15), the associated non-zero partial derivatives of g = (g 1 , g 2 , g 3 , g 4 ) T (at E 0 ) are given by: It follows from the above expressions that: from which it can be shown that â > 0 if: For the sign of b, it can be shown that the associated non-vanishing partial derivatives of g are: > 0.
Thus, we have established Theorem 4 in view of [48].The proof is completed.

Figure 1 .
Figure 1.The device for collecting Rhodopseudomonas palustris in the chemostat by inputting inorganic flocculants.

Figure 4 .
Figure 4. (a) Rhodopseudomonas palustris in the chemostat will tend to extinction with the parameters in TableA1.(b) Rhodopseudomonas palustris in the chemostat will tend to be constant if the absorption of Rhodopseudomonas palustris (r) is improved.The initial conditions are C 0 = 2, N 0 = 4, X 0 = 1, P 0 = 3.

Figure 5 .
Figure 5. (a) Rhodopseudomonas palustris in the chemostat will tend to be constant if the flocculation effect (m 1 ) is reduced.(b) Rhodopseudomonas palustris in the chemostat will tend to be constant if the Michaelis-Menten constant of carbon (K 1 ) is reduced.The initial conditions are C 0 = 2, N 0 = 4, X 0 = 1, P 0 = 3.