1. Introduction
About 100 years ago, A. J. Lotka [
1] and V. Volterra [
2] independently developed a mathematical model, which nowadays serves as the mathematical background for population dynamics, ecology, chemical reactions, etc. Their model is based on a system of ordinary differential equations (ODEs) with quadratic nonlinearities (typically two equations). The natural generalization of this model in 1D space reads as follows:
where the lower subscripts
t and
x mean differentiation with respect to (w.r.t.) these variables,
and
are two unknown functions, which usually represent densities,
and
are arbitrary constants (some of them can vanish and different types of interactions arise depending on signs of nonvanish constants), and
and
are diffusion coefficients. System (
1) is called the diffusive Lotka–Volterra (DLV) system and is the main object of this work. If the diffusivities are such that
, then (
1) reduces to the classical Lotka–Volterra system.
In contrast to the classical Lotka–Volterra system, the DLV system attracted the attention of scholars much later. Its rigorous study started in the 1970s (see the pioneering works [
3,
4,
5,
6]). At the present time, there are many recent works devoted to qualitative analysis of the DLV system (
1) and its multi-component analogs (see [
7,
8] and the works cited therein). However, the number of the papers devoted to construction of exact solutions of the nonlinear system (
1) is relatively small. Exact solutions in the form of traveling waves were constructed in [
9,
10,
11,
12]. In the case of the three-component DLV system, some traveling waves were found in [
13,
14]. The existence of traveling wave solutions were examined in [
7,
8,
15,
16]. To the best of our knowledge, exact solutions with more complicated structure were derived only in papers [
17,
18] for the two- and three-component DLV systems, respectively. In [
19], nontrivial exact solutions were derived for a natural generalization of system (
1) involving additional linear and/or quadratic terms. We also point out that systems of nonlinear ODEs for finding exact solutions of (
1) in the very special case when
are presented in the handbook [
20]. However, those systems are not solved therein.
Thus, the problem of the construction of exact solutions of the DLV system (
1) and its multi-dimensional analogs, especially those with a biological, physical or chemical interpretation, is a hot topic.
From the very beginning, we point out that the DLV system (
1) is nonlinear; hence, it cannot be integrated in a straightforward way. Here, we examine this system assuming that both equations involve diffusion, are nonlinear and are not autonomous, i.e.,
The most powerful methods for the construction of exact solutions for non-integrable nonlinear partial differential equations (PDEs) are symmetry-based methods. These methods originate from the Lie method, which was created by the prominent Norwegian mathematician Sophus Lie in the 1880s. The Lie method (the terminology ‘the Lie symmetry analysis’ and ‘the group-theoretical analysis’ are also used) still attracts the attention of many researchers, and new results are published on a regular basis (see the recent monographs [
21,
22] and the papers cited therein). However, it is well known that some nonlinear PDEs and systems of PDEs arising in applications have poor Lie symmetry. The Lie method is not efficient for such equations since it enables only those exact solutions to be constructed, which can be easily obtained without using this method. The DLV system (
1) belongs to such systems because one possesses a nontrivial Lie symmetry only under unrealistic restrictions on parameters (see more details in [
10]). As a result, Lie symmetries allow us to construct only traveling wave solutions for (
1).
During recent decades, other symmetry-based methods were developed in order to solve nonlinear PDEs with poor Lie symmetry. The best known among them is the method of nonclassical symmetries proposed by G. Bluman and J. Cole in 1969 [
23]. Notably, following Fushchych’s proposal dating back to the 1980s [
24,
25], we use the terminology ‘
Q-conditional symmetry’ instead of ‘nonclassical symmetry’ (see also a discussion concerning terminology in Chapter 3 of [
22]). Although this method was suggested 50 years ago, its successful applications for solving
nonlinear systems of PDEs were accomplished only in the 2000s, and the majority of such papers were published during the last 10 years (see [
17,
18,
26,
27,
28,
29,
30]). This occurred because application of the nonclassical method (such a terminology was used in [
23] instead of nonclassical symmetries) leads to very complicated nonlinear equations to-be-solved. As a result, one needs to solve a much more complicated PDE system (the so-called system of determining equations (DEs)), comparing with the system in question. In paper [
27], a simpler algorithm was proposed in order to make essential progress in solving systems of DEs and to construct
Q-conditional symmetry. The algorithm is based on the notion of
Q-conditional symmetry of the first type. In paper [
17], we successfully applied the new algorithm for finding new exact solutions of the DLV system (
1).
In this work, using a modification of the algorithm for finding
Q-conditional symmetry of the first type, we make further progress in the construction of new symmetries and exact solutions of the DLV system (
1). Moreover, we demonstrate that some solutions can be useful in population dynamics.
The paper is organized as follows. In
Section 2, we present some definitions and provide a complete description of
Q-conditional symmetries of the first type of the DLV system (
1) in the so-called no-go case. In
Section 3, the symmetries obtained are applied to reduce the DLV system to the systems of ODEs and to construct exact solutions. In
Section 4, the properties of particular exact solutions are examined with the aim to provide their biological interpretation. Finally, we present some conclusions in the last section.
2. Q-Conditional Symmetries of the DLV System
Let us consider the general form of the
Q-conditional symmetry operator of system (
1), namely the first-order operator:
where
and
are smooth functions that can be found, using the well known criterion. Given (
3), one can calculate the second prolongation of the operator
Q:
where the coefficients
and
with relevant subscripts are expressed via the functions
and
by the well-known formulae (see [
21,
31]).
Definition 1. Operator (3) is called the Q-conditional symmetry for the DLV system (1) if the following invariance conditions are satisfied:where the manifoldwhile Since the expressions
and
contain only the derivatives
and
, one needs to consider the following coefficients:
System (
1) is the system of evolution equations. Therefore, the problem of constructing its
Q-conditional symmetries of the form (
3) essentially depends on the value of the function
. Thus, one should consider two different cases:
.
.
In
Case 1, one can set
without loss of generality using the well known property stating that the
Q-conditional symmetry operator can be multiplied by an arbitrary smooth function (see the proof in [
31]). Moreover, in this case, the differential consequences of equations
and
(see the manifold
) w.r.t. the variables
t and
x lead to the second-order PDEs, namely, the following:
We notethat the above equations involve the time derivatives
and the mixed derivatives
, which do not occur in the invariance conditions (
4). As a result, the manifold
can be rewritten as
, i.e., the first-order differential consequences can be omitted.
Case 1 for the DLV system (
1) was investigated in the work [
17] (see also Chapter 3 in [
31]).
Here, we examine
Case 2, for which the terminology
‘no-go case’ is often used. Thus, we are looking for operators of the following form:
In this case, formulae (
5) are essentially simplified and take the following forms:
First of all, we note that the task of constructing the
Q-conditional symmetries with
for scalar evolution equations is equivalent to solving the equation in question [
32]. For this reason, one can obtain only some particular results finding the
Q-conditional symmetry operators of the form (
6) for system (
1).
Our aim is to construct
Q-conditional symmetries of the first type for the DLV system (
1) in
Case 2. The notion of
Q-conditional symmetry of the first type was introduced in the paper [
27] as a special case of
Q-conditional symmetry for
systems of PDEs. Each
Q-conditional symmetry of the first type is automatically a
Q-conditional symmetry (nonclassical symmetry) but not vice versa.
Definition 2. Operator (6) is called the Q-conditional symmetry of the first type for the DLV system (1) if the following invariance conditions are satisfied:where the manifold is either given by or . We point out that the definition was given in [
27] for an arbitrary multi-component system of evolution PDEs, and differential consequences (see above
) were not used therein. In fact, such equations do not play any role if one looks for operators of the form (
3) with
. It was only indicated (see conclusions in [
27]) that differential consequences should be taken into account in the case of arbitrary systems (e.g., involving hyperbolic equations). However, it was not noted in [
27] that one may use differential consequences for searching operators (
3) with
occurring in
Case 2. Here, we show that such an approach leads to new results.
The problem of finding
Q-conditional symmetries of the first type for some reaction–diffusion systems (in particular, two- and three-component DLV systems) are considered in monograph [
31] (see Chapters 3 and 4). In paper [
29], such symmetries were constructed in the no-go case for a wide class of reaction–diffusion systems with nonconstant diffusivities.
Let us apply Definition 2 to construct the system of DEs for finding the
Q-conditional symmetry operators of the form (
6). Firstly, we note that the DLV system (
1) has a symmetric structure and admits the discrete transformation
Thus, to obtain all
Q-conditional symmetries of the first type for the DLV system (
1), it is enough to examine only one manifold from Definition 2, say
. Solving the corresponding system of DEs, the list of inequivalent (up to same specified local transformations) DLV systems and corresponding operators will be derived. In order to obtain a complete classification, the DLV systems from the list derived will be checked to confirm whether they admit additional conditional symmetries satisfying Definition 2 with the manifold
.
Thus, the invariance conditions (
8) corresponding to the manifold
take the forms:
where
and
with indices are calculated by the Formula (
7).
Using the equations generating the manifold
, one can exclude the derivatives
and
:
Note that the derivative
(which can be defined from the equation
) is not presented in conditions (
9). Thus, to construct the system of DEs, one needs to substitute (
10) into (
9) and to split the equations obtained w.r.t.
and
. Omitting straightforward calculations, we present only the following result:
Note that Equations (
13)–(
15) were simplified (using equation
and its differential consequences w.r.t.
x and
u) by excluding the derivatives
and
.
The result of integrating the system of DEs (
11)–(
15) can be formulated as follows.
Theorem 1. The DLV system (1) is invariant under the Q-conditional symmetry operator(s) of the first type (6) if and only if the system and the corresponding operator(s) have the forms listed in Table 1. Any other DLV system (1) admitting a nontrivial Q-conditional symmetry of the first type and the relevant operator(s) are reduced to those listed in Table 1 by local transformations from the following set:where and are some correctly-specified constants. The proof of the theorem is presented in
Appendix A.
Remark 1. In Table 1, the upper indexes u and v mean that the relevant Q-conditional symmetry operators satisfy Definition 2 for the manifolds and , respectively. Remark 2. In Table 1, ; and are arbitrary smooth functions of the corresponding variables, the function is the general solution of the linear ODEthe functions F and G are the general solution of the systemand the function is the general solution of the Burgers equation , whilewhere , , and are arbitrary constants. Remark 3. The general solution of the quasilinear first-order system (17) can be constructed in an implicit form using of computation program, say Maple. In order to avoid cumbersome formulae, here we present its solution in the case :where α is an arbitrary constant, while is an arbitrary smooth function. The above integral is expressed in the terms of elliptic functions, which degenerate to elementary functions in particular cases. For example, if . Remark 4. All the systems arising in Table 1, excepting that in Case 1, are semi-coupled because the second equation is autonomous. Interestingly, the equation in Cases 2, 3, 5 and 7 is nothing else but the Fisher equation [33]. The most interesting system from applicability point of view occurs in Case 1. Using the transformations (see (
16))
the system can be generalized to the following form:
Thus, we conclude that the DLV system (
1) admits exactly two
Q-conditional symmetries of the first type provided,
and
, i.e., has the form (
19). Depending on signs of the parameters, the DLV system (
19) can describe competition or mutualism of two populations of species (cells). However, this system cannot describe the prey–predator interaction because the quadratic terms have the same signs in both equations (see [
34] for the classification of interaction types).
Finally, the following observation should be highlighted. Because each
Q-conditional symmetry of the first type is automatically a usual
Q-conditional (nonclassical) symmetry, all operators listed in
Table 1 are nonclassical symmetries. On the other hand, it can be noted that Cases 7 and 8 of
Table 1 do not present new nonclassical symmetries because the operators
and
are particular cases of those
and
arising in Cases 3 and 4. In fact, setting
and
in
, one obtains exactly the operator
. Similarly,
is a particular case of
. We remind the reader that any
Q-conditional symmetry can be multiplied by an arbitrary smooth function in contrast to the Lie symmetry and
Q-conditional symmetry of the first type. Thus, Cases 7 and 8 of
Table 1 can be skipped if one considers the
Q-conditional (nonclassical) symmetries.
3. Reduction and Exact Solutions
In this section, we present examples of reductions of the DLV system to ODE systems, using the
Q-conditional symmetry from Theorem 1, and solve the ODE systems obtained in order to construct exact solutions of the DLV system. Our aim is to find such exact solutions of the DLV system (
1) that are bounded, nonnegative and satisfy the zero Neumann boundary conditions in some correctly-specified domain (interval).
Consider the DLV system from Case 1 of
Table 1, namely the following:
which is the most interesting from the applicability point of view. Since the
Q-conditional symmetry operators
and
of system (
20) lead to the equivalent solutions (up to discrete transformation
), we use only one of them, namely
. According to the standard procedure, in order to construct the ansatz corresponding to the operator
, one needs to solve the first-order PDE system (see
and
in Definition 1):
Integrating system (
21) for each form of the function
from (
18), one constructs the ansatz as follows:
if
the ansatz
if
and the ansatz is as follows:
if
It can be noted that each ansatz derived above satisfies the simple functional relation
.
Now three reductions of the PDE system in question to the ODE systems can be provided. Substituting the above ansatz into the DLV system (
20), we arrive at the ODE system as follows:
in the case of (
22) and (
23), while the system
is obtained in the case of (
24). Here,
and
are new unknown functions.
It turns out that each of the ODE systems (
25) and (
26) can be integrated by reducing to a single second-order ODE. As a result, the general solution of system (
25) is derived in the following form:
where
and
are arbitrary constants. Substituting the functions
and
from (
27) into formulae (
22) and (
23), one obtains the exact solutions of the DLV system (
20) with
:
if
, and
if
, where
.
Similarly, solving the ODE system (
26), we arrive at the exact solutions
and
of system (
20) with
for
and
, respectively.
To the best of our knowledge, all the exact solutions obtained above are new, although several papers are devoted to finding exact solutions of the two-component DLV system. In fact, a majority of these papers [
9,
10,
11] present the plane wave solutions (traveling waves), which have the following structure:
where
c is the speed of the wave. The solutions derived herein possess a more complicated structure than traveling waves and cannot be presented in the form (
30). So, we compare our results only with paper [
17] (see also Chapter 3 in [
31]), where also nontrivial solutions were constructed for the DLV system as follows:
It can be easily noted that this system is reduced to the form
by the introduction diffusivities
and
(here
and
). Now, one realizes that two nonlinear systems (
20) and (
31) are inequivalent, provided
. It means that any solution derived in [
17] cannot be transformed into a solution of the DLV system (
20). We only point out that the exact solutions (
28) and (
29) with
have the same structure as those of (117)–(118) [
17]. It means that these solutions could be used for description of similar processes. A possible application is presented in the next section.
4. Interpretation of the Solution Obtained
In this section, we present an example that demonstrates remarkable properties of some solutions constructed in the previous section. Obviously, using the transformation
(see Formula (
16)) and introducing the notation
one reduces the DLV system (
20) to the following form:
The nonlinear system (
32) with positive parameters
and
c is widely used for describing the competition of two population of species (or cells) (see, e.g., [
34,
35]). Solution (
28) (we set
just for simplicity) after the above transformation takes the following form:
In order to provide a biological interpretation, the components
u and
v must be bounded and nonnegative in some domain because they describe densities of species. If we consider the domain
, then it can be easily shown that both components are bounded and nonnegative, provided the coefficient restrictions
hold. Moreover, solution (
33) possesses the asymptotical behavior as follows:
Now, one realizes that
and
are steady state points of the competition model (
32) and the asymptotical behavior (
34) is in agreement with the qualitative theory of this model (see [
8] and the papers cited therein).
In real-world applications, the competition takes place in some bounded domain, say,
. Typically, the zero flux conditions are assumed at the boundary of
:
Such boundary conditions reflect a natural assumption that the competing species cannot cross the boundaries (for example, the wide river is a natural obstacle). It can be easily checked that the exact solution (
33) satisfies the boundary conditions only under the following requirement:
where
and
are arbitrary integer parameters and
Thus, we conclude that our solution with correctly-specified parameters describe the competition of two population of species in the bounded domain. An example is presented in
Figure 1.
It should be pointed out that traveling wave solutions, which are widely studied for any nonlinear model and play an important role in qualitative analysis, usually cannot be used for solving the relevant models involving the zero flux boundary conditions in the bounded domains. Let us consider the traveling wave of the DLV system (
32) with
, which was firstly constructed in [
10] (see also Section 3.2.3 in [
31]) and much later rediscovered in [
12] (see formulae (18) and (24) therein)
Clearly, the components
u and
v are bounded and nonnegative, provided that
and the asymptotical behavior is the same as in (
34). However, solution (
36) does not satisfy the boundary condition (
35) for any finite values of
A and
B. It can be done only for
and
. An example of solution (
36) is presented in
Figure 2.
5. Conclusions
In this paper, the two-component DLV system (
1) was examined in order to find
Q-conditional symmetries in the so-called no-go case (see (
6)) and to construct exact solutions and provide their biological meaning.
From the very beginning, we modified the definition of
Q-conditional symmetries of the first type [
27] in the no-go case (see Definition 2). In contrast to the standard definition (see Definition 1), Definition 2 allows us to obtain the integrable system of DEs (
11)–(
15). Solving the system of DEs, the main theoretical result in the form of Theorem 1 was derived. The theorem presents an exhaustive list of
Q-conditional symmetries of the first type, which the DLV system (
1) admits depending on the parameters
and
. All other cases, which are not listed in
Table 1, are reducible to those in
Table 1 by appropriate point transformations of the form (
16).
We used the coefficient restrictions (
2), which are motivated from the mathematical and applied point of view. For example, we excluded from the examination the systems of the form (
1) involving a linear equation, i.e., the following:
It is not plausible that such a system can model any interaction between species (cells, chemicals) because the first equation is linear and autonomous. Interestingly, system (
37) is reduced to the following:
by the local transformation
In the case
, all solutions of system (
38) can be easily derived by substitution of the relevant solutions of the linear diffusion equation into the second equation, which becomes the form of a standard diffusion equation with a linear source (sink), i.e., it is again solvable via classical methods for linear PDEs. In the case
, the situation is more complicated because the second equation is nonlinear w.r.t.
.
The conditional symmetries from
Table 1 allow us to construct exact solutions of the relevant systems. As a result, a variety of new exact solutions of the nonlinear DLV system (
20), which is the most interesting from applicability point of view, were derived. Moreover, it was shown that all the solutions obtained are new.
Finally, we examined a model describing the competition of two populations of species. It was shown that the exact solution (
33) with correctly-specified parameters are bounded, nonnegative and satisfies the zero Neumann boundary conditions at bounded space domains. Moreover, the solution possesses a realistic asymptotic behavior. Thus, we conclude that our solution with correctly-specified parameters describes the competition of two species. Interestingly, the solution is periodical in space (see
Figure 1) in contrast to the known traveling wave (see
Figure 2).
In conclusion, we want to highlight an unsolved problem. The nonlinear system (
32) with positive parameters
and negative
b and
c is the model describing mutualism or cooperation (see [
34,
36]). Obviously, the solutions constructed in this work can be used for these types of interaction in a quite similar way as that in
Section 4. However, these solutions are not applicable for the third most common type of interactions between species (cells), prey–predator models. In the prey–predator model, the parameters should satisfy the conditions
,
and
(see the DLV system (
1)). It can be seen that
Table 1 does not contain such a type of systems; therefore, the relevant exact solutions cannot be found. Moreover, we noted that all the exact solutions derived in other papers [
9,
10,
11,
12,
13,
14,
17,
18] are not applicable for description of the prey–predator interaction as well. Thus, the problem of constructing exact solutions for the DLV system (
1) modeling the interaction between the prey and predator is a hot topic.