On the estimation of sequestered infected erythrocytes in Plasmodium falciparum malaria patients

  • Received: 01 January 2013 Accepted: 29 June 2018 Published: 01 March 2014
  • MSC : Primary: 93C15, 92D30, 93B07, 93B30; Secondary: 93C41, 92C60.

  • The aim of this paper is to give a method for the estimation of total parasite burden of the patient and the rate of infection in a malaria's intra-host model by using control theory tools. More precisely, we use an auxiliary system, called observer or estimator, whose solutions tend exponentially to those of the original model. This observer uses only the available measurable data, namely, the values of peripheral infected erythrocytes. It provides estimates of the sequestered infected erythrocytes, that cannot be measured by clinical methods. Therefore this method allows to estimate the total parasite burden within a malaria patient. Moreover, our constructed observer does not use the uncertain infection rate parameter $\beta$. In fact, we derive a simple method to estimate this parameter $\beta$.We apply this estimation method using real data that have been collected when malaria was used as therapy for neurosyphilis by the US Public Health Service.

    Citation: Derdei Bichara, Nathalie Cozic, Abderrahman Iggidr. On the estimation of sequestered infected erythrocytes in Plasmodium falciparum malaria patients[J]. Mathematical Biosciences and Engineering, 2014, 11(4): 741-759. doi: 10.3934/mbe.2014.11.741

    Related Papers:

    [1] Ayman Elsharkawy, Clemente Cesarano, Abdelrhman Tawfiq, Abdul Aziz Ismail . The non-linear Schrödinger equation associated with the soliton surfaces in Minkowski 3-space. AIMS Mathematics, 2022, 7(10): 17879-17893. doi: 10.3934/math.2022985
    [2] Zuhal Küçükarslan Yüzbașı, Dae Won Yoon . On geometry of isophote curves in Galilean space. AIMS Mathematics, 2021, 6(1): 66-76. doi: 10.3934/math.2021005
    [3] Nural Yüksel, Burçin Saltık . On inextensible ruled surfaces generated via a curve derived from a curve with constant torsion. AIMS Mathematics, 2023, 8(5): 11312-11324. doi: 10.3934/math.2023573
    [4] Xiaoling Zhang, Cuiling Ma, Lili Zhao . On some $ m $-th root metrics. AIMS Mathematics, 2024, 9(9): 23971-23978. doi: 10.3934/math.20241165
    [5] Hassan Al-Zoubi, Bendehiba Senoussi, Mutaz Al-Sabbagh, Mehmet Ozdemir . The Chen type of Hasimoto surfaces in the Euclidean 3-space. AIMS Mathematics, 2023, 8(7): 16062-16072. doi: 10.3934/math.2023819
    [6] Melek Erdoğdu, Ayșe Yavuz . On differential analysis of spacelike flows on normal congruence of surfaces. AIMS Mathematics, 2022, 7(8): 13664-13680. doi: 10.3934/math.2022753
    [7] Atakan Tuğkan Yakut, Alperen Kızılay . On the curve evolution with a new modified orthogonal Saban frame. AIMS Mathematics, 2024, 9(10): 29573-29586. doi: 10.3934/math.20241432
    [8] Abdelkader Laiadi, Abdelkrim Merzougui . Free surface flows over a successive obstacles with surface tension and gravity effects. AIMS Mathematics, 2019, 4(2): 316-326. doi: 10.3934/math.2019.2.316
    [9] B. B. Chaturvedi, Kunj Bihari Kaushik, Prabhawati Bhagat, Mohammad Nazrul Islam Khan . Characterization of solitons in a pseudo-quasi-conformally flat and pseudo- $ W_8 $ flat Lorentzian Kähler space-time manifolds. AIMS Mathematics, 2024, 9(7): 19515-19528. doi: 10.3934/math.2024951
    [10] Özgür Boyacıoğlu Kalkan, Süleyman Şenyurt, Mustafa Bilici, Davut Canlı . Sweeping surfaces generated by involutes of a spacelike curve with a timelike binormal in Minkowski 3-space. AIMS Mathematics, 2025, 10(1): 988-1007. doi: 10.3934/math.2025047
  • The aim of this paper is to give a method for the estimation of total parasite burden of the patient and the rate of infection in a malaria's intra-host model by using control theory tools. More precisely, we use an auxiliary system, called observer or estimator, whose solutions tend exponentially to those of the original model. This observer uses only the available measurable data, namely, the values of peripheral infected erythrocytes. It provides estimates of the sequestered infected erythrocytes, that cannot be measured by clinical methods. Therefore this method allows to estimate the total parasite burden within a malaria patient. Moreover, our constructed observer does not use the uncertain infection rate parameter $\beta$. In fact, we derive a simple method to estimate this parameter $\beta$.We apply this estimation method using real data that have been collected when malaria was used as therapy for neurosyphilis by the US Public Health Service.


    It is well known that many nonlinear phenomena in physics, chemistry and biology are described by dynamics of shapes, such as curves and surfaces, and the time evolution of a curve and a surface has significations in computer vision and image processing. The time evolution of a curve and a surface is described by flows, in particular inextensible flows of a curve and a surface. Physically, inextensible flows give rise to motion which no strain energy is induced. The swinging motion of a cord of fixed length or of a piece of paper carried by the wind, can be described by inextensible flows of a curve and a surface. Also, the flows arise in the context of many problems in computer visions and computer animations [2,7]

    In [5], Hasimoto studied the relation between integrable systems and geometric curve flows and he showed that the non linear Schrödinger equation is equivalent to the binormal notion flow of space curves by using a transformation relating the wave function of the Schrödinger equation to the curvature and torsion of curves (so-called Hasimoto transformation). In particular, Lamb [8] proved that the mKdV equation and the sine-Gordon equation arise from the invariant curve flows by using the Hasimoto transformation. After, Mohamed [9] investigated the general description of the binormal motion of a spacelike and a timelike curve in a 3-dimensional de-Sitter space and gave some explicit examples of a binormal motion of the curves. Schief and Rogers [11] investigated binormal motions of inextensible curves with a constant curvature or a constant torsion, and introduced two new examples of integrable equations which are derived from the binormal motion of curves and established the Bäcklund transformation and matrix Darboux transformations of the extended Dym equation and the mKdV equation. In [1] the authors studied curve motions by the binormal flow with the curvature and the torsion depending velocity and sweeping out immersed surfaces and obtained filaments evolving with a constant torsion which arise from extremal curves of curvature energy functionals. Curve flows have been studied by many experts and geometers [4,6,13].

    The outline of the paper is organized as follows: In Section 2, we give some geometric concepts of curves and surfaces in an isotropic 3-space. In Section 3, we study inextensible flows of a space curve and give time evolutions of the Frenet frame, the curvature and the torsion of the curve. In the last section, we construct the Bäcklund transformations of the Schrödinger flows and the extended Harry-Dym flows as the binormal flows and give a nonexistence of bi-harmonic Hasimoto surfaces in an isotropic 3-space.

    The three dimensional isotropic space $ \Bbb I^3 $ has been developed by Strubecker in the 1940s, and it is based on the following group $ G_6 $ of an affine transformations $ (x, y, z) \rightarrow(\bar x, \bar y, \bar z) $ in $ \Bbb R^3 $,

    $ ˉx=a+xcosϕysinϕ,ˉy=b+xsinϕ+ycosϕ,ˉz=c+c1x+c2y+z,
    $

    where $ a, b, c, c_1, c_2, \phi \in \Bbb R $. Such affine transformations are called isotropic congruence transformations or isotropic motions of $ \Bbb I^3 $, in this case $ \mathcal{B}_6 $ is denoted by the group of isotropic motions (cf. [12]). Observe that on the $ xy $ plane this geometry looks exactly like the plane Euclidean geometry. The projection of a vector $ {\bf x} = (x_1, y_1, z_ 1) $ in $ \Bbb I^3 $ on the $ xy $ plane is called the top view of $ \bf x $ and we shall denote it by $ \widehat{ {\bf x} } = (x_1, y_1, 0) $. The top view concept plays a fundamental role in the isotropic space $ \Bbb I^3 $, since the $ z $ direction is preserved under action of $ \mathcal {B}_6 $. A line with this direction is called an isotropic line and a plane that contains an isotropic line is said to be an isotropic plane.

    In the sequel, many of metric properties in isotropic geometry (invariants under $ \mathcal{B}_6 $) are Euclidean invariants in the top view such as the isotropic distance, so call $ i $-distance. The isotropic distance of two points $ {P} = (x_1, y_1, z_1) $ and $ {Q} = (x_2, y_2, z_2) $ is defined as the Euclidean distance of the their top views, i.e.,

    $ d(P,Q):=(x1x2)2+(y1y2)2.
    $
    (2.1)

    As a fact, two points $ (x, y, z_i) $ $ (i = 1, 2) $ with the same top views have isotropic distance zero, they called parallel points.

    Let $ \mathbf{x} = (x_1, y_1, z_1) $ and $ \mathbf{y} = (x_2, y_2, z_2) $ be vectors in $ \Bbb I^3 $. The isotropic inner product of $ \mathbf{x} $ and $ \mathbf{y} $ is defined by

    $ x,y={z1z2,ifxi=0andyi=0x1x2+y1y2,if otherwise.
    $
    (2.2)

    We call vector of the form $ {\bf x} = (0, 0, z) $ in $ \Bbb I^3 $ isotropic vector, and non-isotropic vector otherwise.

    Now we introduce some terminology related to curves. A regular curve $ C : I \rightarrow \Bbb I^3 $, i.e., $ C' \neq 0 $, is parametrized by an arc-length $ s $ if $ || C'|| = 1. $ In the following we assume that all curves are parametrized by an arc-length $ s $. In addition, a point $ \alpha(s_0) $ in which $ \{ C'(s_0), C''(s_0)\} $ is linearly dependent is an inflection point and a regular unit speed curve $ C(s) = (x(s), y(s), z(s)) $ with no inflection point is called an admissible curve if $ x' y'' - x'' y' \neq 0 $.

    On the other hand, the (isotropic) unit tangent, principal normal, and curvature function of the curve $ C $ are defined as usual

    $ t(s)=C(s),n(s)=t(s)κ(s),andκ(s)=||t(s)||=||^t(s)||,
    $
    (2.3)

    respectively. As usually happens in isotropic geometry, the curvature $ \kappa $ is just the curvature function of its top view $ \hat{ C} $ and then we may write $ \kappa(s) = (x' y''- x'' y')(s). $ To compute the moving trihedron, we define the binormal vector as the vector $ {\bf b} = (0, 0, 1) $ in the isotropic direction. The Frenet equation corresponding to the isotropic Frenet frame $ \{ {\bf t}, {\bf n}, {\bf b} \} $ can be written as ([12])

    $ dds(tnb)=(0κ0κ0τ000)(tnb),
    $
    (2.4)

    where $ \tau $ is the (isotropic) torsion, that is,

    $ \tau = \frac {\det (C', C'', C''')}{\kappa^2}. $

    Consider a $ C^r $-surface $ M $, $ r\geq 1 $, in $ \Bbb I^3 $ parameterized by

    $ {X}(u, v) = ( x(u, v), y(u, v), z(u, v)). $

    A surface $ {M} $ immersed in $ \Bbb I^3 $ is called admissible if it has no isotropic tangent planes. We restrict our framework to admissible regular surfaces.

    For such a surface, the coefficients $ g_{ij} $ $ (i, j = 1, 2) $ of its first fundamental form are given by

    $ g_{11} = \langle X_u, X_u \rangle, \quad g_{12} = \langle X_u, X_v \rangle, \quad g_{22} = \langle X_v, X_v \rangle, $

    where $ X_u = \frac {\partial X}{\partial u} $ and $ X_v = \frac {\partial X}{\partial v}. $ The coefficients $ h_{ij} $ $ (i, j = 1, 2) $ of the second fundamental form of $ M $ are calculated with respect to the normal vector of $ M $ and they are given by

    $ h_{11} = \frac {\det (X_{uu} X_u X_v) }{\sqrt{\det(g_{ij})}}, \quad h_{12} = \frac {\det (X_{uv} X_u X_v) }{\sqrt{\det(g_{ij})}}, \quad h_{22} = \frac {\det (X_{vv} X_u X_v) }{\sqrt{\det(g_{ij})}}. $

    The isotropic Gaussian curvature $ K $ and the isotropic mean curvature $ H $ are defined by

    $ K=det(hij)det(gij),H=g11h222g12h12+g22h112det(gij).
    $
    (2.5)

    The surface $ {M} $ is said to be isotropic flat (resp. isotropic minimal) if $ {K} $ (resp.$ {H} $) vanishes (cf. [12,14]).

    We assume that $ C : [0, l] \times [0, w] \rightarrow M \subset \Bbb I_3 $ is a one parameter family of a space curve in $ \Bbb I_3 $, where $ l $ is the arc-length of a initial curve. Let $ u $ be the curve parametrization variable, $ 0 \leq u \leq l. $ We put $ v = || \frac {\partial C}{\partial u}|| $, from which the arc-length of $ C $ is defined by $ s(u) = \int_0^u v d u. $ Also, the operator $ \frac {\partial }{\partial s} $ is given in terms of $ u $ by $ \frac {\partial }{\partial s} = \frac 1 v \frac {\partial }{\partial u} $ and the arc-length parameter is given by $ ds = v du. $

    On the Frenet frame $ \{ \bf t, \bf n, \bf b \} $ of the curve $ C $ in $ \Bbb I_3 $, any flow of $ C $ can be given by

    $ Ct=αt+βn+γb,
    $
    (3.1)

    where $ \alpha, \beta, \gamma $ are scalar speeds of the space curve $ C $, respectively. We put $ s(u, t) = \int ^u _0 v du $, it is called the arc-length variation of $ C $. From this, the requirement that the curve not be subject to any elongation or compression can be expressed by the condition

    $ ts(u,t)=u0vtdu=0
    $
    (3.2)

    for all $ u \in [0, l]. $

    Definition 3.1. A curve evolution $ C(u, t) $ and its flow $ \frac {\partial C}{\partial t} $ in $ \Bbb I_3 $ are said to be inextensible if

    $ \frac {\partial }{\partial t}\left|\left|\frac {\partial C}{\partial u}\right|\right| = 0. $

    Theorem 3.2. If $ \frac {\partial C}{\partial t} = \alpha {\bf t} + \beta {\bf n} + \gamma {\bf b} $ is a flow of $ C $ in an isotropic 3-space $ \Bbb I_3 $, then we have the following equation:

    $ vt=αuvκβ.
    $
    (3.3)

    Proof. Since $ v^2 = \left \langle \frac {\partial C}{\partial u}, \frac {\partial C}{\partial u} \right \rangle $ and $ u $, $ t $ are independent coordinates, $ \frac {\partial}{\partial u} $ and $ \frac {\partial}{\partial t} $ commute. So by differentiating of $ v^2 $ with respect to $ t $ and using (2.4) and (3.1) we can easily obtain (3.3).

    Now, we give necessary and sufficient condition for a inextensible flow in an isotropic space and it is useful to get our results.

    Theorem 3.3. Let $ \frac {\partial C}{\partial t} = \alpha {\bf t} + \beta {\bf n} + \gamma {\bf b} $ be a flow of a space curve $ C $ in $ \Bbb I_3 $. Then the flow is inextensible if and only if

    $ αs=κβ.
    $
    (3.4)

    Proof. Suppose that the curve flow of a space curve $ C $ is inextensible. From (3.2) and (3.3) we have

    $ \frac {\partial }{\partial t} s(u, t) = \int^u_0 \frac {\partial v}{\partial t} du = \int^u_0 \left( \frac {\partial \alpha}{\partial u} - v \kappa \beta\right) du = 0, \quad {\text for \; all} \; u \in [0, l]. $

    It follows that

    $ \frac {\partial \alpha}{\partial u} = v \kappa \beta. $

    Since $ \frac {\partial }{\partial s} = \frac 1 v \frac {\partial }{\partial u} $, we can obtain (3.4).

    Conversely, by following a similar way as above, the proof is completed.

    Theorem 3.4. Let $ \frac {\partial C}{\partial t} = \alpha {\bf t} + \beta {\bf n} + \gamma {\bf b} $ be a flow of a space curve $ C $ in $ \Bbb I_3 $. If the flow is inextensible, then a time evolution of the Frenet frame $ \{ \bf t, \bf n, \bf b \} $ along a curve $ C $ is given by

    $ ddt(tnb)=(0φ1φ2φ100000)(tnb),
    $
    (3.5)

    where

    $ φ1=βs+ακ,φ2=γs+βτ.
    $
    (3.6)

    Proof. Noting that

    $ tt=t(Cs)=s(αt+βn+γb)=(βs+ακ)n+(γs+βτ)b.
    $
    (3.7)

    Also,

    $ 0=tt,n=tt,n+t,nt=(βs+ακ)+t,nt,0=tn,b=nt,b,
    $

    which imply that a time evolution of the principal normal vector $ \bf n $ can be expressed as

    $ \frac {\partial \bf n}{\partial t} = - (\frac {\partial \beta}{\partial s} + \alpha \kappa) {\bf t}. $

    This completes the proof.

    Now, we give time evolution equations of the curvature $ \kappa $ and the torsion $ \tau $ of the inextensible space curve $ C $ in $ \Bbb I_3 $.

    Theorem 3.5. Let $ \frac {\partial C}{\partial t} = \alpha {\bf t} + \beta {\bf n} + \gamma {\bf b} $ be a flow of a space curve $ C $ in $ \Bbb I_3 $. Then, the time evolution equations of the functions $ \kappa $ and $ \tau $ for the inextensible space curve $ C $ are given by

    $ κt=φ1s,τt=κφ2.
    $
    (3.8)

    Proof. It is well known that the arc length and time derivatives commute. That is,

    $ s(tt)=s(φ1n+φ2b)=(κφ1)t+(φ1s)n+(φ2s+φ1τ)b,
    $

    and

    $ t(ts)=t(κn)=(κφ1)t+(κt)n.
    $

    Comparing two equations, we find

    $ \frac {\partial \kappa}{\partial t} = \frac {\partial \varphi_1}{\partial s}. $

    Also by using $ \frac {\partial }{\partial s } \left(\frac {\partial {\bf n}}{\partial t} \right) = \frac {\partial }{\partial t } \left(\frac {\partial {\bf n}}{\partial s} \right) $ and following a similar way as above, we can obtain the second equation of (3.8). The proof is completed.

    Remark 3.6. Taking $ \beta = -\kappa_s $, by (3.4) one find $ \alpha = - \frac 1 2 \kappa^2 $. From the time evolution of the curvature (3.8) we get

    $ \kappa_t = - \kappa_{sss} - \frac 3 2 \kappa^2 \kappa_s, $

    it follows that the curvature $ \kappa $ evolves according to the mKdV equation

    $ \kappa_t + \kappa_{sss} + \frac 3 2 \kappa^2 \kappa_s = 0, $

    where $ \kappa_s = \frac {\partial \kappa }{\partial s} $ and $ \kappa_t = \frac {\partial \kappa }{\partial t} $. The corresponding flow of a curve is

    $ C_t = - \frac 1 2 \kappa^2 {\bf t} - \kappa_s {\bf n}, $

    which is the so-called modified KdV flow [3].

    We give a example of inextensible flow of a curve with constant torsion as follows:

    Example 3.7. Taking $ \tau = \tau_0 $ = constant. If we consider

    $ \alpha = \alpha_0, \quad \beta = \beta_0, \quad \gamma = - \beta_0 \tau_0 s + \gamma_0 t, $

    where $ \alpha_0, \beta_0, \gamma_0 $ are non zero constants, then the PDE system (3.8) takes the form

    $ \kappa_t = \alpha_0 \kappa_s. $

    it follows that one solution of the last equation is

    $ \kappa(s, t) = e^{s + \alpha_0 t}. $

    Thus, (2.4) and (3.5) imply $ \alpha_0 {\bf t}_s = {\bf t}_t $ and $ {\bf n} = e^{-(s + \alpha_0 t)} {\;\bf t}_s $. If we take

    $ {\bf t} = ( \cosh( \alpha_0 s + t ), \sinh(\alpha_0 s + t), \alpha_0 s + t), $

    the vector $ {\bf n} $ is given by

    $ {\bf n}(s, t) = \alpha_0 e^{(s + \alpha_0)} ( \sinh( \alpha_0 s + t ), \cosh(\alpha_0 s + t), 1). $

    Thus, we can get the family of curves $ C_t $, so we can determine the surface that is generated by this family of curves.

    Qu et al. [10] studied the Bäcklund transformations of geometric curve flows in a Euclidean 3-space $ \Bbb R^3 $. In this section, we aim to give the Bäcklund transformations of integrable geometric curve flows by using time evolutions of a curve in an isotropic 3-space $ \Bbb I^3 $.

    It is well-known that the Schrödinger flow in a 3-space is given by [5]

    $ Ct=Cs×Css=κb.
    $
    (4.1)

    In this case we take $ (\alpha, \beta, \gamma) = (0, 0, \kappa) $ in (3.1), it follows that (3.5) and (3.8) imply the following theorem:

    Theorem 4.1. The Schrödinger flow (4.1) implies the time evolutions of frame fields, the curvature and the torsion of a space curve $ C $ in an isotropic 3-space $ \Bbb I_3 $ as follows:

    $ ddt(tnb)=(00κs000000)(tnb),
    $
    (4.2)
    $ ddt(κτ)=(00κs0)(κτ).
    $
    (4.3)

    We now construct the Bäcklund transformation of the Schrödinger flow (4.1). Considering another curve in $ \Bbb I^3 $ related $ C $ by

    $ ˜C(s,t)=C(s,t)+μ(s,t)t+ν(s,t)n+ξ(s,t)b,
    $
    (4.4)

    where $ \mu, \nu $ and $ \xi $ are the smooth functions of $ s $ and $ t $. Using (2.4), (4.1) and (4.2), a direct computation leads to

    $ ˜Cs=(1+μsκν)t+(νs+κμ)n+(ξs+τν)b,˜Ct=μtt+νtn+(κ+ξt+κsμ)b.
    $
    (4.5)

    Let $ \tilde{s} $ be the arclength parameter of the curve $ \tilde{C} $. Then

    $ d\tilde{s} = || \tilde{C_s}|| ds = \sqrt{(1 + \mu_s - \kappa \nu )^2 + ( \nu_s + \kappa \mu)^2} ds : = \Omega ds, $

    where $ \Omega $ is a non zero smooth function. It follows that the unit tangent vector of the curve $ \tilde{C} $ is determined by

    $ ˜t=Φ1t+Φ2n+Φ3b,
    $
    (4.6)

    where $ \Phi_1 = \Omega^{-1} (1 + \mu_s - \kappa \nu), \Phi_2 = \Omega^{-1} (\nu_s + \kappa \mu) $ and $ \Phi_3 = \Omega^{-1} (\xi_s + \tau \nu) $. Differentiating (4.5) with respect to $ \tilde{s} $, we get

    $ \tilde{ {\bf t}}_{\tilde s} = \frac {\Phi_{1s} - \kappa \Phi_2}{\Omega} {\bf t} + \frac {\Phi_{2s}+ \kappa \Phi_1}{\Omega} {\bf n} + \frac {\Phi_{3s} + \tau \Phi_2}{\Omega} {\bf b} $

    which gives the curvature of the curve $ \tilde{C} $:

    $ ˜κ=||^˜Cs||=(Φ1sκΦ2)2+(Φ2s+κΦ1)2Ω:=ΘΩ.
    $
    (4.7)

    Thus form (2.3) the Frenet frames of the curve $ \tilde{C} $ are given by

    $ ˜n=Φ1sκΦ2Θt+Φ2s+κΦ1Θn+Φ3s+τΦ2Θb,˜b=(0,0,1).
    $
    (4.8)

    Assume that the curve $ \tilde{C} $ also fulfills the Schrödinger flow, that is,

    $ ˜Ct=˜κ˜b.
    $
    (4.9)

    The Bäcklund transformation of the Schrödinger flow with the help of (4.5), (4.8) and (4.9) turns out to be the following result (cf. [10]):

    Theorem 4.2. The Schrödinger flow (4.1) is invariant with respect to the Bäcklund transformation (4.4) in an isotropic 3-space $ \Bbb I_3 $ if $ \mu, \nu $ and $ \xi $ satisfy the system

    $ μt=0,νt=0,κ+ξt+κsμ=ΘΩ.
    $

    The extended Harry-Dym flow [11]

    $ Ct=τ12b
    $
    (4.10)

    is obtained by setting $ \alpha = 0 $, $ \beta = 0 $ and $ \gamma = \tau^{-\frac 12} $ in a space curve flow (3.1).

    Theorem 4.3. The extended Harry-Dym flow (4.10) implies the time evolutions of frame fields, the curvature and the torsion of a curve $ C $ in an isotropic 3-space $ \Bbb I_3 $ as follows:

    $ ddt(tnb)=(00(τ12)s000000)(tnb),
    $
    (4.11)
    $ ddt(κτ)=(00(τ12)s0)(κτ).
    $
    (4.12)

    We consider the Bäcklund transformation of the extended Harry-Dym flow (4.10)

    $ ˜C(s,t)=C(s,t)+μ(s,t)t+ν(s,t)n+ξ(s,t)b,
    $
    (4.13)

    where $ \mu, \nu $ and $ \xi $ are the smooth functions of $ s $ and $ t $. Differentiating (4.13) with respect to $ s $ and $ t $ and using (4.11) we get

    $ ˜Cs=(1+μsκν)t+(νs+κμ)n+(ξs+τν)b,˜Ct=μtt+νtn+(τ12+ξt+μ(τ12)s)b.
    $
    (4.14)

    Let $ \tilde{s} $ be the arclength parameter of the curve $ \tilde{C} $. Then the unit tangent vector of the curve $ \tilde{C} $ is given by

    $ ˜t=Φ1t+Φ2n+Φ3b,
    $
    (4.15)

    where we put

    $ Φ1=Ω1(1+μsκν),Φ2=Ω1(νs+κμ),Φ3=Ω1(ξs+τν),Ω=(1+μsκν)2+(νs+κμ)2.
    $

    Differentiating (4.15) with respect to $ \tilde{s} $, we get

    $ \tilde{ {\bf t}}_{\tilde s} = \frac {\Phi_{1s} - \kappa \Phi_2}{\Omega} {\bf t} + \frac {\Phi_{2s}+ \kappa \Phi_1}{\Omega} {\bf n} + \frac {\Phi_{3s} + \tau \Phi_2}{\Omega} {\bf b}, $

    it follows that the curvature of the curve $ \tilde{C} $ leads to

    $ ˜κ=||^˜Cs||=(Φ1sκΦ2)2+(Φ2s+κΦ1)2Ω:=ΘΩ.
    $
    (4.16)

    Thus form (2.3) the Frenet frames of the curve $ \tilde{C} $ are given by

    $ ˜n=Φ1sκΦ2Θt+Φ2s+κΦ1Θn+Φ3s+τΦ2Θb,˜b=(0,0,1).
    $
    (4.17)

    Assume that the curve $ \tilde{C} $ also fulfills the extended Harry-Dym flow, that is,

    $ ˜Ct=˜τ12˜b.
    $
    (4.18)

    The Bäcklund transformation of the extended Harry-Dym flow with the help of (4.14) and (4.17) turns out to be the following result:

    Theorem 4.4. The extended Harry-Dym flow (4.10) is invariant with respect to the Bäcklund transformation (4.13) in an isotropic 3-space $ \Bbb I_3 $ if $ \mu, \nu $ and $ \xi $ satisfy the system

    $ μt=0,νt=0,τ12+ξt+μ(τ12)s=ΘΩ.
    $

    Definition 4.5. The surface $ C(s, t) $ wiped out by the Schrödinger flow (4.1) is called a Hasimoto surface.

    Sometimes, Eq. (4.1) is called the vortex filament or smoke ring equation, and can be viewed as a dynamical system on the space curves. Since $ {C}_s = {\bf t} $ and $ {C}_t = \kappa {\bf b}, $ the first and second fundamental forms of the Hasimoto surface are given by

    $ I=ds2+κ2dt2,II=κds2.
    $

    In this case $ \kappa $ is non-vanishing everywhere. Thus, we have

    Theorem 4.6. Let $ {C}(s, t) $ be a Hasimoto surface such that $ {C}(s, t) $ is a unit speed curve for all $ t $ in an isotropic 3-space $ \Bbb I_3 $. Then the surface has zero Gaussian curvature and the mean curvature $ H = -\frac 1 2 \kappa $.

    Corollary 4.7. There is no minimal Hasimoto surface in an isotropic 3-space.

    Theorem 4.8. Let $ {C}(s, t) $ be a Hasimoto surface in an isotropic 3-space $ \Bbb I_3 $. Then the followings are satisfy:

    (1). $ s $-parameter curves of the surface are non-asymptotic curves.

    (2). $ s $-parameter curves of the surface are geodesic curves.

    Proof. Suppose that $ {C}(s, t) $ is a Hasimoto surface such that $ {C}(s, t) $ is a unit speed curve for all $ t $. It is well known that the normal curvature $ \kappa_n $ of a $ s $-parameter curve on the surface $ {C}(s, t) $ is given by

    $ \kappa_n = \langle U, {C}_{ss} \rangle = \langle -{\bf n}, \kappa {\bf n} \rangle = -\kappa, $

    where $ U = \frac {{C}_s \times {C}_t} {||{C}_s \times {C}_t||} $ is a unit normal vector of the surface. Since $ \kappa \neq 0 $, a $ s $-parameter curve is non-asymptotic.

    Also, the geodesic curvature $ \kappa_g $ of a $ s $-parameter curve on the surface $ {C}(s, t) $ is

    $ \kappa_g = \langle U \times {C}_s, {C}_{ss} \rangle = -\langle {\bf n}\times {\bf t}, \kappa {\bf n} \rangle = 0, $

    it follows that a $ s $-parameter curve is geodesic. The proof is completed.

    Theorem 4.9. Let $ {C}(s, t) $ be a Hasimoto surface in an isotropic 3-space $ \Bbb I_3 $. Then the followings are satisfy:

    (1). $ t $-parameter curves of the surface are asymptotic curves.

    (2). $ t $-parameter curves of the surface are geodesic curves.

    Proof. The normal curvature $ \kappa_n $ and the geodesic curvature $ \kappa_g $ of a $ t $-parameter curve on the surface $ {C}(s, t) $ are given by

    $ κn=U,Ctt=n,κtb=0,κg=U×Ct,Ctt=n×κb,κtb=0.
    $

    Thus, a $ t $-parameter curve is asymptotic and geodesic. The proof is completed.

    On the other hand, the mean curvature vector $ \bf H $ is given by

    $ H=HU=12κn,
    $
    (4.19)

    and Laplacian of the mean curvature vector is expressed as

    $ ΔH=(2κκs)t+12κ(κ4κ2sκκss)n12(3κsτ+κτs)b.
    $
    (4.20)

    Theorem 4.10. Let $ {C}(s, t) $ be a Hasimoto surface such that $ {C}(s, t) $ is a unit speed curve for all $ t $ in an isotropic 3-space $ \Bbb I_3 $. If the surface satisfies the condition $ \Delta {\bf H} = \lambda {\bf H} $ for some constant $ \lambda $, then the curve $ {C}(s, t) $ for all $ t $ has a constant curvature and a constant torsion. Furthermore, $ \lambda = \kappa^2 $.

    Proof. Suppose that a Hasimoto surface satisfies the condition $ \Delta {\bf H} = \lambda {\bf H} $. With the help of (4.19) and (4.20), we obtain the following equations:

    $ κκs=0,κ4κ2sκκssλκ2=0,3κsτ+κτs=0,
    $
    (4.21)

    which imply $ \kappa_s = 0 $ and $ \tau_s = 0 $. From (4.3) we also get $ \kappa_t = 0 $ and $ \tau_t = 0 $, that is, $ \kappa $ and $ \tau $ are constant. Furthermore, from (4.21) $ \lambda = \kappa^2 $. Thus, the proof is completed.

    Corollary 4.11. There is no bi-harmonic Hasimoto surface in an isotropic 3-space.

    Example 4.12. Let us consider the surface in $ \Bbb I^3 $ parametrized by

    $ {C}(s, t) = (\cos s, \sin s, \sin s-\cos s +s +t) $

    with $ 0 \leq s \leq 2 \pi $ and $ -2 \leq t \leq 2 $. Then, we have

    $ t(s,t)=(sins,coss,coss+sins+1),n(s,t)=(coss,sins,sins+coss),b(s,t)=(0,0,1),κ(s,t)=1,τ(s,t)=1.
    $

    On the other hand, we can check the following:

    $ Ct=(0,0,1)=κb,κt=0,τt=0.
    $

    Thus the surface is a Hasimoto surface (Figure 1) and it satisfies $ \Delta {\bf H} = {\bf H}. $

    Figure 1.  Hasimoto surface with constant curvature and torsion.

    The authors would like to thank the anonymous referee for useful comments and suggestions that helped for a significant improvement of the exposition.

    This paper was supported by Basic Science Research Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Education (NRF-2018R1D1A1B07046979).

    The authors declare no conflicts of interest.

    [1] Parasitology, 99 (1989), S59-S79.
    [2] Am. J. Trop. Med. Hyg., 61 (1999), 36-43.
    [3] T. Roy. Soc. Trop. Med. H., 95 (2001), 497-501.
    [4] Parasitology, 117 (1998), 97-105.
    [5] Parasitology, 117 (1998), 409-410.
    [6] Parasitology, 110 (1995), 115-122.
    [7] Proc. Natl. Acad. Sci. USA., 95 (1998), 7620-7624.
    [8] J. Theor. Biol., 217 (2002), 137-148.
    [9] Parasitology, 113 (1996), 25-38.
    [10] IEEE Trans. Automat. Control, 37 (1992), 871-875.
    [11] IEEE Trans. Automat. Control, 39 (1994), 1338-1341.
    [12] SIAM J. Appl. Math., 67 (2006), 260-278.
    [13] Proc. Natl. Acad. Sci. USA., 88 (1991), 5111-5113.
    [14] Am. J. Trop. Med. Hyg., 61 (1999), 367-374.
    [15] J. Theor. Biol., 198 (1999), 549-566.
    [16] Proc. Natl. Acad. Sci. USA., 101 (2004), 9161-9166.
    [17] Parassitologia, 41 (2000), 221-231.
    [18] T. Roy. Soc. Trop. Med. H., 96 (2002), 205-209.
    [19] Parasitology, 131 (2005), 449-458.
    [20] Proc. Natl. Acad. Sci. USA., 100 (2003), 3473-3478.
    [21] Parasitol. Today., 12 (1996), 74-79.
  • This article has been cited by:

    1. Yi Zheng, Xiaoqun Wu, Ziye Fan, Kebin Chen, Jinhu Lü, An Improved Topology Identification Method of Complex Dynamical Networks, 2025, 55, 2168-2267, 2165, 10.1109/TCYB.2025.3547772
  • Reader Comments
  • © 2014 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3004) PDF downloads(481) Cited by(15)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog