1. Introduction
The Poisson-Nernst-Planck model is a well-known model of ion transport, which plays a crucial role in the study of many biological and physical problems, such as ion channels in cell membranes [1] [2] and semiconductor devices [3].
A steady-state Poisson-Nernst-Planck model [4] [5] [6] is
(1.1)
where
is the electric potential;
is the concentration for the ith ion species;
is the valence;
is the permanent charge of the channel;
is the electrochemical potential;
is the area of the cross-section of the channel;
is the flux density;
is the diffusion coefficient;
is the relative dielectric coefficient;
is the vacuum permittivity; k is the Boltzmann constant; T is the absolute temperature; and e is the elementary charge.
The boundary conditions are, for
,
(1.2)
in the classical Poisson-Nernst-Planck model takes the following form
(1.3)
with
is a constant.
Migration of charges for ionic flow through ion channels is often described mathematically by the Poisson-Nernst-Planck model (1.1), which can be viewed as a simplified version of the Maxwell-Boltzmann equations [7] [8] and the Langevin-Poisson equations [9] [10] by focusing on the key features of biological functions. Recently, the model (1.1) has been greatly studied [11] - [17]. In [18], the author obtained the existence and uniqueness of solutions for systems (1.1) and (1.2) under the assumption that
and
. In [19], the authors completely solved the existence and uniqueness of solutions for the boundary value problems (1.1) and (1.2) under the assumption that
and n ions with the different valences are involved. However, the analysis of the dynamics of Poisson-Nernst-Planck model with nonzero permanent charges
is much more difficult. In [20], the authors justified the existence of multiple solutions for the boundary value problems (1.1) and (1.2) under the assumption that
is a piecewise constant function and
. For the case that n ions with the arbitrary valences are involved and
is a piecewise constant function, the general geometric framework for analyzing the dynamics of systems (1.1) and (1.2) is provided in [21] based on the geometric singular perturbation theory [22] [23] [24]. In this paper, we intend to study the dynamics of the classical Poisson-Nernst-Planck model under the following hypotheses.
(H1).
and
.
(H2).
for
,
for
and
for
, where Q is a constant.
By re-scaling,
The model (1) is reduced to a standard singularly perturbed system of the following
(1.4)
with the boundary condition, for
.
(1.5)
Actually, for the case
, system (1.4) and (1.5) corresponds to the equations studied in [20]. Additionally, due to the above hypothesis (H1), system (1.4) and (1.5) is a special case of the equations studied in [21], but in this special case, the explicit formulae for the solutions of its limit slow system can be obtained, which is crucial for the analysis of the existence of solutions for systems (1.4) and (1.5) in this paper by combining the technique of the geometric singular perturbation theory.
2. Limiting Fast Orbits and Limiting Slow Orbits on [0, 1]
Let
,
. System (1.4) becomes
(2.6)
By using the rescaling
, one has
(2.7)
Define
(2.8)
Due to the fact that
is a piecewise constant function, therefore, we identify the limiting fast and slow orbits on three intervals
,
and
respectively.
2.1. Limiting Fast Orbits on [0, a]
Let
,
,
,
, where
are unknowns to be determined later. Let
We will identify the limiting fast and limiting slow orbits connecting
to
on the interval
, where
. Letting
in (2.7), one gets the limiting fast system
(2.9)
By letting
in (2.6), we obtain the critical manifold
(2.10)
which is normally hyperbolic.
The flow of
under system (2.9) in forward time is denoted by
, and the flow of
under system (2.9) in backward time is denoted by
. Then the following results can be established [20] [21].
Lemma 2.1. System (2.9) has the following
nontrivial first integrals:
Proposition 2.2. (i) The stable manifold
intersects
transversally at points
and the
-limit set of
is
where
for
are arbitrary, and
(ii) The unstable manifold
intersects
transversally at points
and the
-limit set of
is
where
for
are arbitrary, and
Remark 2.3. At
, the limiting fast orbits
are a segment connecting
to
, and at
the limiting fast orbits
are a segment connecting
to
.
2.2. Limiting Slow Orbits on [0, a]
Now we identify the limiting slow orbits
on the critical manifold
. By using a rescaling
System (2.6) becomes
(2.11)
where
, and its limiting slow system is
(2.12)
For system (2.12), the critical manifold is
where
.
It follows that the limiting slow system on
is
(2.13)
For convenience, we denote
and
Lemma 2.4. There is a unique solution
of (2.13) such that
and
where
are given in Proposition 2.2. It is given by
(2.14)
Proof. By system (2.13), it follows that
Therefore,
By inserting the above formula for
into system (2.13) and using the variation of constants formula, the formulas for
in the statement can be obtained. ¨
2.3. Limiting Fast Orbits on [a, b]
Let
,
,
,
, where
are unknowns to be determined later. Let
We will identify the limiting fast and limiting slow orbits connecting
to
on the interval
, where
. The limiting fast system is obtained by letting
in (2.7):
(2.15)
By letting
in (2.6), we obtain the critical manifold
(2.16)
which is normally hyperbolic.
The flow of
under system (2.15) in forward time is denoted by
, and the flow of
under system (2.15) in backward time is denoted by
. Then the following results can be established [20] [21].
Lemma 2.5. System (2.15) has the following
nontrivial first integrals:
Proposition 2.6. (i) Letting
be the unique solution of
and letting
The stable manifold
intersects
transversally at points
and the
-limit set of
is
where
for
are arbitrary.
(ii) Letting
be the unique solution of
and letting
The unstable manifold
intersects
transversally at points
and the
-limit set of
is
where
for
are arbitrary.
Remark 2.7. At
, the limiting fast orbits
are a segment connecting
to
, and at
the limiting fast orbits
are a segment connecting
to
.
2.4. Limiting Slow Orbits on [a, b]
Now we identify the limiting slow orbits
on the critical manifold
. By using a rescaling
System (2.6) becomes
(2.17)
where
, and its limiting slow system is
(2.18)
For system (2.18), the critical manifold is
where
.
It follows that the limiting slow system on
is
(2.19)
Following the idea in [20] [21], system (2.19) can be transformed to
(2.20)
The solution of (2.20) is
(2.21)
where
It can be seen that there exists a
such that
, which implies that
and
,
,
, then it follows that
(2.22)
2.5. Limiting Fast Orbits on [b, 1]
In this section, we will identify the limiting fast and limiting slow orbits connecting
to
on the interval
, where
. The limiting fast system is obtained by letting
in (2.7):
(2.23)
By letting
in (2.6), we obtain the critical manifold
(2.24)
which is normally hyperbolic.
The flow of
under system (2.23) in forward time is denoted by
, and the flow of
under system (2.23) in backward time is denoted by
. Then the following results can be established [20] [21].
Lemma 2.8. System (2.23) has the following
nontrivial first integrals:
Proposition 2.9. (i) The stable manifold
intersects
transversally at points
and the
-limit set of
is
where
for
are arbitrary, and
(ii) The unstable manifold
intersects
transversally at points
and the
-limit set of
is
where
for
are arbitrary, and
Remark 2.10. At
, the limiting fast orbits
are a segment connecting
to
, and at
the limiting fast orbits
are a segment connecting
to
.
2.6. Limiting Slow Orbits on [b, 1]
Now we identify the limiting slow orbits
on the critical manifold
. Just as in sections 2.1 and 2.2, it can be shown that the limiting slow system is
(2.25)
For convenience, we denote
and
Lemma 2.11. There is a unique solution
of (2.25) such that
and
where
are given in Proposition 2.9. It is given by
(2.26)
3. Results
Based on Propositions 2.2, 2.6, 2.9 and formulas (2.14), (2.22), (2.26), we get
(3.27)
and
(3.28)
where
Now, we consider a special case that
, and
. Following the idea in [20], let
and
Then system (3.27) reduces to
(3.29)
where
and system (3.28) reduces to
(3.30)
Additionally, it is demonstrated in [20] that
has solutions, therefore, it follows that the following unknowns
can be determined. The remaining unknown
will be determined by Equation (3.30).
By solving (3.30), we get
(3.31)
where
and
Remark 3.1. Once
are determined, then
are also determined. Therefore, all unknowns involved in Equations (3.27) and (3.28) are determined.
Therefore, a limiting fast and limiting slow orbit is identified as follows, see Figure 1 for an illustration.
By employing the Exchange Lemma [20] [21] [22] [23] [24], it can be verified that
Theorem 3.2. For
sufficiently small, there exists a unique solution of
Figure 1. A limiting fast and limiting slow orbit connecting
to
, where
are limiting fast orbits and
are limiting slow orbits.
system (1.4) and (1.5) near the limiting fast and limiting slow orbits.
4. Conclusion
In this paper, a steady-state Poisson-Nernst-Planck model with n ion species is studied under the assumption that
positively charged ion species have the same valence and there is only one negatively charged ion species. By using the geometric theory for singularly perturbed system, the existence of solutions for systems (1.4) and (1.5) is justified. As we know, it can be seen that the results in [20] correspond to those in this paper in the case that
. Also, this paper shows that the number of solutions for systems (1.4) and (1.5) essentially is determined by the number of solutions for the algebraic Equation (3.29), that is, an increase in the number of positively charged ions with the same valence does not change the number of solutions for systems (1.4) and (1.5). Generally, for the case that there are more than two species ions involved in the Poisson-Nernst-Planck model, the dynamics become more subtle and complicated. Moreover, the mixture of multi-species ions, such as sodium (Na+), potassium (K+), calcium (Ca2+), chloride (Cl−), plays the very important role in many vital biological functions, for instance, opening and closing of ionic channels. In [25] [26], it was shown that the Poisson-Nernst-Planck model with three or more ionic species of different charge may admit multiple steady state solutions, and the existence of multiple steady state solutions is important to study transitions between such states which may be related to the gating (switching between open and closed states of ionic channels) and selectivity of ion channels.
Acknowledgements
The author thank the reviewers for their invaluable comments and suggestions which help improve the manuscript. The author is supported by the NNSFC 11971477.