In this paper, we investigate the global dynamics of a viral infection model with defectively infected cells. The explicit expression of the basic reproduction number of virus is obtained by using the next generation matrix approach, where each term has a clear biological interpretation. We show that the basic reproduction number serves as a threshold parameter. The virus dies out if the basic reproduction number is not greater than unity, otherwise the virus persists and the viral load eventually approaches a positive number. The result is established by Lyapunov's direct method. Our novel arguments for the stability of the infection equilibrium not only simplify the analysis (compared with some traditional ones in the literature) but also demonstrate some correlation between the two Lyapunov functions for the infection-free and infection equilibria.
Citation: Jianquan Li, Xiaoyu Huo, Yuming Chen. Threshold dynamics of a viral infection model with defectively infected cells[J]. Mathematical Biosciences and Engineering, 2022, 19(7): 6489-6503. doi: 10.3934/mbe.2022305
[1] | Kenneth M. Madden, Boris Feldman, Shane Arishenkoff, Graydon S. Meneilly . Bedside tests of muscle mass in older adults with Type 2 diabetes. AIMS Medical Science, 2022, 9(3): 433-446. doi: 10.3934/medsci.2022022 |
[2] | Tyler C. Smith, Besa Smith . Consistency in Physical Activity and Increase in Mental Health in Elderly over a Decade: Are We Achieving Better Population Health?. AIMS Medical Science, 2016, 3(1): 147-161. doi: 10.3934/medsci.2016.1.147 |
[3] | Breanne E. Kunstler, Paul D O’Halloran, Jill L. Cook, Joanne L. Kemp, Caroline F. Finch . “…like you’re pushing the snowball back up hill”—the experiences of Australian physiotherapists promoting non-treatment physical activity: A qualitative study. AIMS Medical Science, 2018, 5(3): 224-237. doi: 10.3934/medsci.2018.3.224 |
[4] | Carol J. Farran, Olimpia Paun, Fawn Cothran, Caryn D. Etkin, Kumar B. Rajan, Amy Eisenstein, and Maryam Navaie . Impact of an Individualized Physical Activity Intervention on Improving Mental Health Outcomes in Family Caregivers of Persons with Dementia: A Randomized Controlled Trial. AIMS Medical Science, 2016, 3(1): 15-31. doi: 10.3934/medsci.2016.1.15 |
[5] | Lena Fleig, Maureen C. Ashe, Jan Keller, Sonia Lippke, Ralf Schwarzer . Putting psychology into telerehabilitation: Coping planning as an example for how to integrate behavior change techniques into clinical practice. AIMS Medical Science, 2019, 6(1): 13-32. doi: 10.3934/medsci.2019.1.13 |
[6] | Christina L. Ekegren, Rachel E. Climie, William G. Veitch, Neville Owen, David W. Dunstan, Lara A. Kimmel, Belinda J. Gabbe . Sedentary behaviour and physical activity patterns in adults with traumatic limb fracture. AIMS Medical Science, 2019, 6(1): 1-12. doi: 10.3934/medsci.2019.1.1 |
[7] | Hidetaka Hamasaki . Dog ownership: is it beneficial for physical activity, cardiovascular disease, and diabetes?. AIMS Medical Science, 2020, 7(4): 311-327. doi: 10.3934/medsci.2020021 |
[8] | Paula Barrett-Brown, Donovan McGrowder, Dalip Ragoobirsingh . Diabetes education—Cornerstone in management of diabetes mellitus in Jamaica. AIMS Medical Science, 2021, 8(3): 189-202. doi: 10.3934/medsci.2021017 |
[9] | Nathacha Garcés, Angel Jara, Felipe Montalva-Valenzuela, Claudio Farías-Valenzuela, Gerson Ferrari, Paloma Ferrero-Hernández, Antonio Castillo-Paredes . Motor performance in children and adolescents with attention deficit and hyperactivity disorder: A systematic review. AIMS Medical Science, 2025, 12(2): 247-267. doi: 10.3934/medsci.2025017 |
[10] | Carol Dillon, Fernando E. Taragano . Special Issue: Activity and Lifestyle Factors in the Elderly: Their Relationship with Degenerative Diseases and Depression. AIMS Medical Science, 2016, 3(2): 213-216. doi: 10.3934/medsci.2016.2.213 |
In this paper, we investigate the global dynamics of a viral infection model with defectively infected cells. The explicit expression of the basic reproduction number of virus is obtained by using the next generation matrix approach, where each term has a clear biological interpretation. We show that the basic reproduction number serves as a threshold parameter. The virus dies out if the basic reproduction number is not greater than unity, otherwise the virus persists and the viral load eventually approaches a positive number. The result is established by Lyapunov's direct method. Our novel arguments for the stability of the infection equilibrium not only simplify the analysis (compared with some traditional ones in the literature) but also demonstrate some correlation between the two Lyapunov functions for the infection-free and infection equilibria.
In this note we consider the following magnetic, quantum Euler equations:
{∂tρ+divJ=0∂tJ+div(J⊗Jρ)+∇P(ρ)=J∧B+12ρ∇(Δ√ρ√ρ), | (1.1) |
in the unknowns (ρ,J):Rt×R3→R⩾0×R3, representing the charge and current densities of a positively charged quantum fluid, interacting with an external, time-independent magnetic field B. The third order nonlinear term in the right hand side of the momentum balance encodes the quantum correction, while J∧B is the magnetic force associated to the field B (we are neglecting here the self-induced electromagnetic field). We assume B=∇∧A for some real-valued, divergence-free potential A. Finally, P(ρ) is a barotropic pressure term, which assumes the form P(ρ)=γ−1γργ, with γ>1. From now on, we focus without loss of generality only on forward-in-time solutions (i.e., defined for t⩾0) and we denote the initial data at time t=0 by (ρ0,J0).
System (1.1), as well as suitable variations, arise in the description of hydrodynamic models in which the quantum statistic becomes relevant at mesoscopic/macroscopic scales, which happens when the thermal de Broglie wavelength is comparable to (i.e., larger than or equal to) the average interatomic distance [1]; important examples include superfluidity [2], quantum plasma [3], and nanoscale semiconductor devices [4]. Since the constituent particles (e.g. fermions, such as electrons, and bosons as the Cooper pairs in superconductivity [5]) of a quantum fluid are typically charged, in many situations it becomes necessary to account for electromagnetic effects in the hydrodynamic description. A relevant model is provided by the quantum Euler-Maxwell equations, which describe the dynamics of a positively charged quantum fluid interacting with its self-induced electromagnetic field, see [5,6,7,8] and references therein for a comprehensive analysis.
On the other hand, as in the classical Ginzburg-Landau theory of superconductivity, the presence of an external magnetic field (especially of large amplitude) can drastically affect the behavior (both dynamical and at equilibrium) of a quantum fluid, even when the self-generated electromagnetic field may be neglected. In this perspective, we mention the work [9], which deals with a quantum fluid model for spinning fermions under the influence of an external electromagnetic field and discusses several applications, including quantum X-ray free-electron lasers and micro-mechanical systems. We also point out that in [9] the quantum fluid equations are not only introduced as effective model, but they are also formally derived from a microscopic, many-body description. In fact, an important issue in hydrodynamic and condensed-matter physics concerns the possibly rigorous derivation of effective equations from first principles. In this direction we mention [10], which gives a rigorous derivation of the incompressible 2D Euler equation from the von Neumann equation with an external magnetic field, and [11] where a 3D quantum fluid model with self-induced electromagnetic field is derived from the bosonic Pauli-Fierz Hamiltonian.
In the search for a rigorous derivation of a quantum hydrodynamic model, a preliminary (yet crucial) step is to have at hand a suitable solution theory, possibly at the whole finite-energy level (i.e., without requiring any extra regularity or smallness assumption on the initial data). Moreover, in order to deal with physically meaningful solutions, all the corresponding observables should be well defined.
Motivated by the above discussion, in this paper we aim to construct finite energy solutions to the magnetic quantum Euler equations (1.1) (without restrictions on the size of the initial data) such that the magnetic force J∧B is well-defined. We will deal only with weak solutions, since this allows us to consider a wide class of (possibly large and non-smooth) magnetic fields, see Remark 2.2.
System (1.1) admits two quantities that are formally preserved along the flow: the total charge ‖ρ‖L1=‖√ρ‖2L2 and the total energy
E(t)=12∫R3|∇√ρ|2+|Λ|2+2f(ρ)dx, | (1.2) |
where Λ:=J√ρ and f(ρ):=ργγ is the internal energy density. Observe that, for γ⩽3, we have f(ρ)∈L1(R3) whenever √ρ∈H1(R3). We restrict to the subcritical regime γ<3, while the critical case γ=3 requires a slightly finer analysis, which for concreteness will not be addressed here. Then, in view of System (1.2), the natural finite-energy condition for solutions to System (1.1) is
√ρ∈H1(R3),Λ∈L2(R3). |
We work in a quasi-irrotational framework. More precisely, for almost every time we require the generalized irrotationality condition
∇∧J+ρB=2∇√ρ∧Λ | (1.3) |
to be satisfied as an identity in D′(R3). In order to clarify the meaning of the above condition, consider a smooth solution (ρ,J) to System (1.1), for which we can define the velocity field v through the relation J=ρv. A direct computation then shows that System (1.3) is equivalent to the identity
ρ(∇∧v+B)=0, |
which implies in particular that the relative velocity field vA:=v+A is irrotational outside vacuum, namely ∇∧vA=0 on {ρ>0}. On the other side, condition (1.3) is compatible with the occurrence of quantum vortices, namely solutions such that vA has non-zero (quantized) circulations around a vacuum region. Formation of vortices plays a fundamental role in the study of quantum fluids [12], hence the importance of developing a theory which embodies solutions carrying a non-trivial vorticity.
Our last requirement is that a solution to System (1.1) must be regular enough to provide a meaning (e.g., local integrability) to the magnetic force J∧B. It is worth pointing out that the finite-energy condition alone does not provide a priori a well-defined magnetic force. Indeed, by Sobolev embedding and Hölder inequality,
(√ρ,Λ)∈H1(R3)×L2(R3)⇒J=√ρΛ∈L32(R3), |
which guarantees that J∧B is locally integrable only when B∈L3loc(R3). We will be able instead, by means of suitable a priori space-time dispersive estimates, to deal with magnetic fields that are almost finite-energy, namely B∈L2+(R3).
As anticipated, the goal of this note is to show the existence of global in time, finite-energy, weak solutions to System (1.1), without assuming any smallness condition on the initial data. Motivated by the above discussion, we provide the rigorous notion of solution. To this aim, observe that the convective term in the momentum equation can be rewritten as
ρ−1(J⊗J)=Λ⊗Λ, | (1.4) |
while the quantum term can be written in divergence form as
12ρ∇(Δ√ρ√ρ)=div(14∇2ρ−∇√ρ⊗∇√ρ). | (1.5) |
Definition 1.1. Let 0<T⩽∞, ρ0,J0∈L1loc(R3). We say that (ρ,J) is a finite-energy weak solution to System (1.1) on the space-time slab [0,T)×R3, with initial data (ρ0,J0), if the following conditions are satisfied:
(i) (finite energy hydrodynamic state) there exists a pair (√ρ,Λ), with √ρ∈L∞loc(0,T;H1(R3)) and Λ∈L∞loc(0,T;L2(R3)), such that ρ=(√ρ)2, J=√ρΛ;
(ii) (continuity and momentum equations - weak formulations) for every η∈C∞c([0,T)×R3) and ζ∈C∞c([0,T)×R3;R3),
∫T0∫R3ρ∂tη+J⋅∇ηdxdt+∫R3ρ0(x)η(0,x)dx=0; | (1.6) |
∫T0∫R3J∂tζ+∇ζ:(∇√ρ⊗∇√ρ+Λ⊗Λ)+(J∧B)⋅ζ−14ρΔdivζ+P(ρ)divζdxdt+∫R3J0(x)⋅ζ(0,x)dx=0; | (1.7) |
(iv) (well-defined magnetic force) J∧B∈L1loc((0,T)×R3);
(v) (generalized irrotationality condition) for a. e. t∈(0,T), the generalized irrotationality condition (1.3) holds true in D′(R3).
When T=∞ we say that the solution is global.
The quantum term in System (1.1) induces a quadratic dispersion relation in the irrotational part of the flow, as can be readily verified by linearizing around a constant state. Nevertheless, this enhanced dispersion (comparing with the linear dispersion of the classical Euler equations) is still not sufficient to directly derive suitable a priori bounds and prove the existence of solutions for large, finite energy initial data. We exploit instead a suitable wave-function dynamics associated to System (1.1). More precisely, if a complex-valued function ψ satisfies the magnetic NLS equation
i∂tψ=−12ΔAψ+f′(|ψ|2)ψ, | (2.1) |
where ΔA:=∇A⋅∇A:=(∇−iA)⋅(∇−iA)=Δ−2iA⋅∇−|A|2 is the magnetic Laplacian, then its Madelung transform
ρ=|ψ|2,J=Im(¯ψ∇Aψ) | (2.2) |
formally solves System (1.1). Suppose indeed that we can write ψ=√ρeiθ so that J=ρ(∇θ−A) for smooth functions ρ,θ, with ρ>0. Substituting this ansatz into the NLS (2.1), and separating real and imaginary part, we obtain a system of equations that reduces after some manipulations to System (1.1).
The connection between quantum fluid models and wave-functions dynamics governed by Schrödinger-type equations, originally observed in the seminal work by Madelung [13], has been widely exploited in the mathematical literature: we mention [14,15] and [16,17,18,19,20], see also [21] and references therein.
In particular, the existence of global, finite-energy solutions for the 3D quantum Euler equations (i.e., System (1.1) without the magnetic component J∧B in the momentum balance) has been proved in [14], for arbitrarily large initial data and allowing for vacuum. In this case, the corresponding wave-function dynamics is governed by the semilinear NLS
i∂tψ=−12Δψ+f′(|ψ|2)ψ, | (2.3) |
for which a satisfactory well-posedness theory at the finite energy level (i.e., H1-regularity) is available (see e.g., the monograph [22] and references therein). Still, the precise construction of the hydro-dynamic variables and the rigorous derivation of the quantum Euler equations require a fine analysis, based upon the polar factorization technique (see the discussion in Section 3). Existence of global weak solutions has been also proved for aforementioned quantum Euler-Maxwell equations [23], see also [24]. The wave-function dynamics is given now by a nonlinear Maxwell-Schrödinger (MS) system, for which well-posedness is available only above the finite energy level [25]. Still, relying on the existence of finite energy weak solutions to MS [26] and suitable a priori estimates, in [23] the quantum Euler-Maxwell equations are derived by directly working with the weak formulations. This approach has been adopted in [18] for the quantum Euler equations, and is in turn inspired by ideas in [27,28]. Here, we will follow an analogous strategy.
In order to state our main result, we first need some assumptions. Concerning the magnetic potential A, we impose the following condition:
A∈L∞(R3)∩Hσ(R3)for some σ>1. | (2.4) |
As anticipated, in terms of the magnetic field B=∇∧A we are just above the finite-energy level B∈L2(R3). The extra integrability assumption on A will be crucial to obtain the appropriate a priori estimates (see Proposition 3.1) for weak solutions to the magnetic NLS (2.1). In addition, Condition (2.4) guarantees [29, Lemma 2.4] that the magnetic Sobolev space
H1A(R3):={ψ∈L2(R3):∇Aψ∈L2(R3)},‖ψ‖H1A=‖ψ‖L2+‖∇Aψ‖L2, |
coincides (with equivalence of norms) with the classical Sobolev space H1(R3). Since the energy of the magnetic NLS (2.1) is given by
E(t)=12∫R3|∇Aψ|2+2f(|ψ|2)dx, | (2.5) |
we then deduce (in the sub-critical regime γ<3) that H1(R3) is the energy space for Eq (2.1).
Concerning the initial data (ρ0,J0), we impose that they are the Madelung transform of a finite energy wave function:
ρ0=|ψ0|2,J0=Im(¯ψ0∇Aψ0) for some ψ0∈H1(R3). | (2.6) |
The above condition is natural from the quantum mechanics perspective, as it is equivalent [30] to assume that the (relative) initial velocity field v0,A has quantized circulations around vacuum, namely
∫γv0,Ads∈2πZ | (2.7) |
for any closed, smooth path in ρ>0.
We are now able to state our main result.
Theorem 2.1. Assume that the magnetic potential A and the initial data (ρ0,J0) satisfy Systems (2.4) and (2.6), respectively. Then there exists a global, finite energy weak solution to System (1.1), in the sense of Definition 1.1.
Remark 2.2. The existence of strong (i.e., continuous in time) solutions to System (1.1) is unknown in general. Indeed, the typical dispersive bounds, such as Strichartz estimates, which would guarantee the well-posedness in H1(R3) of the corresponding NLS (2.1) are available only for some restricted class of magnetic potentials, satisfying either smoothness, smallness, or non-resonant assumptions (see e.g., [31,32,33] and references therein).
Remark 2.3. In the case of the quantum Euler-Maxwell equations, the L∞-bound (suitably averaged over time) on A follows by means of the dispersive estimates for its wave-type evolution, as soon as the initial magnetic potential is slightly more regular than finite energy. Here, dealing with a fixed, external magnetic field, such condition is imposed a priori. An interesting open question that arises is whether the integrability assumption on A could be weakened without imposing any extra smallness, regularity or spectral condition. This would require new technical ideas to deal with the derivative term 2iA⋅∇ in the expansion of the magnetic Laplacian.
The solution constructed in Theorem 2.1 will be obtained as Madelung transform of a time-dependent wave-function satisfying System (2.1). As such, the quantized vorticity condition will be automatically satisfied at (almost every) time.
Finally, we mention that it would be possible to prove a weak-stability result for System (1.1), in the same spirit of [23,Theorem 1.3]. For the sake of concreteness, we do not pursue this direction here.
The next section is devoted to the proof of Theorem 2.1. We preliminary fix some notations. We will write A≲ if A\leqslant CB for some C > 0 , and A\lesssim_k B if the above constant C depends on the parameter k . Given s > 0 , p\geqslant 1 , we denote by H^{s, p}({{\mathbb R}}^3) the classical Sobolev (Bessel-potential) space. Given p\in[1, \infty] , T > 0 and a Banach space \mathcal{X} , we write L_T^p\mathcal{X} to denote the Bochner space L^p((0, T);\mathcal{X}) .
Our starting point is the existence of a global, finite energy weak solution
\begin{equation} \psi\in L_{\mathrm{loc}}^{\infty}({{\mathbb R}},H^1({{\mathbb R}}^3))\cap W^{1,\infty}_{\mathrm{loc}}({{\mathbb R}},H^{-1}({{\mathbb R}}^3)) \end{equation} | (3.1) |
to the magnetic NLS (2.1), with initial data \psi_0 . This has been proved in [29, Theorem 1.2], where a wider class of rough, time-dependent magnetic potentials is actually considered.
We are going to show that the Madelung transform (\rho, J) of \psi , defined by System (2.2), is a global, finite energy, weak solution to System (1.1), in the sense of Definition (1.1). There are two main issues.
First of all, since the velocity field v is not well-defined on the vacuum regions, the hydrodynamic variable \Lambda , formally given by \Lambda = \sqrt{\rho}v , needs to be carefully realized. To this aim, we briefly recall the polar factorization technique (see e.g., [24]). Consider the set
\begin{equation} P(\psi): = \{\varphi\in L^{\infty}({{\mathbb R}}^3)\mbox{ s.t. } \|\varphi\|_{L^{\infty}}\leqslant 1,\,\psi = |\psi|\varphi\mbox{ a.e. in }{{\mathbb R}}^3\} \end{equation} | (3.2) |
of polar factors of \psi . Given \varphi\in P(\psi) , let us define
\Lambda = \text{Im}(\bar{\varphi}\nabla_A\psi)\in L^2({{\mathbb R}}^3). |
Since \nabla_A\psi = 0 almost everywhere on \{\rho = 0\} [24, Theorem 6.19], the above definition tells us that the \Lambda is uniquely defined, with \Lambda = 0 on the vacuum. As a consequence, the identity J = \sqrt{\rho}\Lambda and formulas (1.4) and (1.5) are rigorously justified, and we deduce that the energy (1.2) of the hydrodynamic state (\sqrt{\rho}, \Lambda) coincides with the energy (2.5) of the wave function \psi .
The second issue is related to the rigorous derivation of System (1.1) from the corresponding NLS Equation (2.3), which could be directly justified only for sufficiently regular (above the finite energy level) wave functions. Since we can not appeal to the H^1({{\mathbb R}}^3) -stability (recall indeed that the well-posedness in H^1({{\mathbb R}}^3) for System (2.3) is unknown in general, as discussed in Remark 2.2), we can not rely on any approximation procedure. We derive instead the weak formulations (1.6) and (1.7) of System (1.1) directly from the Duhamel (integral) formula of System (2.3). To this aim, we need a suitable a priori estimate, similar to [23, Proposition 3.1].
Proposition 3.1. For any T\in(0, \infty) it holds the estimate
\begin{equation} \|\psi\|_{L_T^2H^{1/2,6}}\lesssim_T 1. \end{equation} | (3.3) |
Proof. In view of the Koch-Tzvetkov type estimates [25, Lemma 2.4] and the classical inhomogeneous Strichartz estimate [34] for the Schrödinger equation, we deduce
\|\psi\|_{L_T^2H^{\frac12,6}}\lesssim_T \|\psi\|_{L_T^{\infty}H^1}+\|A\nabla\psi+|A|^2\psi\|_{L_T^2L^2}+\big\|f'(|\psi|^2)\psi\big\|_{L_T^2H^{\frac12,\frac65}}. |
Using assumptions (2.4) on A we get
\|A\nabla\psi+|A|^2\psi\|_{L_T^2L^2}\lesssim_T\|A\|_{L^{\infty}}\|\psi\|_{L_T^{\infty}H^1}. |
Moreover, since \gamma < 3 , a direct application of Sobolev embedding, Gagliardo-Niremberg and Hölder inequalities yields
\big\|f'(|\psi|^2)\psi\big\|_{L_T^2H^{\frac12,\frac65}} = \big\||\psi|^{2(\gamma-1)}\psi\big\|_{L_T^2H^{\frac12,\frac65}}\lesssim_T \|\psi\|_{L_T^{\infty}H^1}^{1+\varepsilon}\|\psi\|_{L_T^2H^{\frac12,6}}^{1-\varepsilon} |
for some \varepsilon = \varepsilon(\sigma)\in(0, 1) . Since \|\psi\|_{L_T^{\infty}H^1}\lesssim_T 1 , we then obtain
\|\psi\|_{L_T^2H^{\frac12,6}}\lesssim_T 1+\|\psi\|_{L_T^2H^{\frac12,6}}^{1-\varepsilon}, |
which implies Eq (3.3), as desired.
Remark 3.2. The above argument can be rigorously justified by considering a sufficiently regular \psi , for which the finiteness of \|\psi\|_{L_T^2H^{\frac12, 6}} is guaranteed, and then appealing to a standard density argument. Alternatively, one could exploit a bootstrap procedure as in [23, Proposition 3.1].
We can prove now that the pair of hydrodynamic variables (\rho, J) satisfies conditions (ⅰ)–(ⅳ) of Definition (1.1), with T = \infty .
(ⅰ) Since the map \psi\to (\sqrt{\rho}, \Lambda) is bounded (and actually continuous) [24, Lemma 5.3] from H^1({{\mathbb R}}^3) into H^1({{\mathbb R}}^3)\times L^2({{\mathbb R}}^3) , we then deduce from Eq (3.1) that
(\sqrt{\rho},\Lambda)\in L_{\mathrm{loc}}^{\infty}({{\mathbb R}},H^1({{\mathbb R}}^3)\times L^2({{\mathbb R}}^3)). |
(ⅱ) The weak formulations of the continuity and momentum equations can be obtained following the same lines as in [23, Section 4]: The extra-integrability (almost L^{\infty} ) on \psi provided by Eq (3.3), and the extra-regularity (above finite energy) on B assumed in (2.4) allow in fact to rigorously justify the derivation of Eqs (1.6) and (1.7) from the integral (Duhamel) formulation of the magnetic NLS (2.1).
(ⅲ) In view of Sobolev embedding, H^{\sigma-1}({{\mathbb R}}^3)\hookrightarrow L^{p}({{\mathbb R}}^3) and H^{\frac12, 6}({{\mathbb R}}^3)\hookrightarrow L^{q}({{\mathbb R}}^3) , for suitable p > 2 and q < \infty such that \frac1p+\frac1q = \frac12 . Then, for every T > 0 ,
\begin{equation*} \begin{split} \|J\wedge B\|_{L_T^1L^1}&\lesssim_T \|\psi\|_{L_T^2L^q}\|\nabla\psi\|_{L_T^{\infty}L^2}\|B\|_{L^p}+\|\psi\|_{L_T^{\infty}L^4}^2\|A\|_{L^{\infty}}\|B\|_{L^2}\\ &\lesssim \|\psi\|_{L_T^2H^{\frac12,6}}\|\psi\|_{L_T^{\infty}H^1}\|A\|_{H^{\sigma}}+\|\psi\|_{L_T^{\infty}H^1}^2\|A\|_{L^{\infty}\cap H^1}\lesssim_T 1, \end{split} \end{equation*} |
which implies in particular J\wedge B\in L^1_{\mathrm{loc}}({{\mathbb R}}\times{{\mathbb R}}^3) .
(ⅳ) We start by observing that
\nabla\wedge J = \text{Im}(\nabla\bar{\psi}\wedge\nabla\psi)-\nabla|\psi|^2\wedge A-B|\psi|^2, |
which yields
\begin{equation} \nabla\wedge J+\rho B = \text{Im}(\nabla\bar{\psi}\wedge\nabla\psi)-2\sqrt{\rho}\,\nabla\sqrt{\rho}\wedge A. \end{equation} | (3.4) |
On the other hand, given a polar factor \varphi of \psi , we have
\begin{equation} 2\nabla\sqrt{\rho}\wedge\Lambda = 2\nabla\sqrt{\rho}\wedge\text{Im}(\bar{\varphi}\nabla\psi)-2\nabla\sqrt{\rho}\wedge(A\sqrt{\rho}). \end{equation} | (3.5) |
Moreover, we have the chain of identities
\begin{equation} \begin{split} \text{Im}(\nabla\bar{\psi}\wedge\nabla\psi)& = 2\text{Im}(\bar{\varphi}\nabla\sqrt{\rho}\wedge\sqrt{\rho}\,\nabla\phi) = 2\nabla\sqrt{\rho}\wedge\text{Im}(\bar{\varphi}\nabla\psi). \end{split} \end{equation} | (3.6) |
Combining Eqs (3.4)–(3.6), we deduce that the generalized irrotationality condition (1.3) is satisfied for almost every time as a distributional identity.
The proof of Theorem 2.1 is complete.
We have shown how to construct global, finite energy weak solutions to a quantum fluid system with an external (possibly large and non-smooth) magnetic field, following similar lines as in the case of the quantum Euler-Maxwell system. Analogous arguments can be also used to deal with the case when both an external magnetic field and the self-induced electromagnetic fields are taken into account in the dynamics of a charged quantum fluid. Finally, let us point out that the overall strategy exploited in this work appears to be versatile and could potentially be adapted to deal with more general dispersive hydrodynamic models, whose solution theory is still not well understood. A relevant example is provided by the Euler-Korteweg system
\begin{equation} \begin{cases} \partial_t \rho+\operatorname{div}(\rho v) = 0, \\ \partial_t (\rho v)+\operatorname{div}(\rho v \otimes v)+ \nabla p(\rho) = \rho\nabla \operatorname{div}\Big((k(\rho)\nabla\rho)-{ {\frac{1}{2}}}k'(\rho)|\nabla\rho|^2\Big), \end{cases} \end{equation} | (4.1) |
where k(\rho)\geqslant 0 is the capillarity coefficient. System (4.1) arises, among the others, in the description of capillarity effects in diffuse interfaces [35]. The choice k(\rho) = \frac{1}{4\rho} leads in particular to the quantum Euler equations. Also, the Euler-Korteweg system is connected, through the Madelung transform, to a wave-function dynamics, given in this case by the quasi-linear NLS
\begin{equation} \mathrm{i} \partial_t \psi = { {-\frac{1}{2}}} \Delta \psi-\operatorname{div}\big(\tilde{k}(|\psi|^2) \nabla|\psi|^2\big)\psi+{ {-\frac{1}{2}}}\tilde{k}^{\prime}(|\psi|^2)|\nabla|\psi|^2|^2\psi+f'(|\psi|^2)\psi, \end{equation} | (4.2) |
where \tilde{k}(\rho) = k(\rho)-\frac{1}{4\rho} . Comparing with the magnetic NLS (2.1), where the Laplace operator is perturbed by the derivative term 2 \mathrm{i} A\cdot\nabla , in Eq (4.2) there appears a second-order term with non-constant and \psi -dependent coefficient, which makes it even harder to prove suitable dispersive estimates. In fact, the global well-posedness for the 3D Euler-Korteweg system is only known under the strong ellipticity condition \rho\tilde{k}(\rho)\geqslant c > 0 at sufficiently high regularity ( H^s with s > \frac{7}{2} ), without vacuum, and for small, spatially localized initial data [36]. On the other side, it is conceivable that weak H^1 -solutions to System (4.2) can be constructed by means of compactness arguments and in turn that low regularity, weak solutions to the Euler-Korteweg system can be derived by exploiting the strategy presented here, which only relies upon dispersive estimates with loss of derivatives.
The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.
The author would like to thank the anonymous referees for the useful suggestions. The author acknowledges support to INdAM-GNAMPA through the project "Local and nonlocal equations with lower order terms".
The author declares there is no conflict of interest.
[1] |
M. A. Nowak, S. Bonhoeffer, A. M. Hill, R. Boehme, H. C. Thomas, H. Mcdade, Viral dynamics in hepatitis B virus infection, Proc. Nati. Acad. Sci., 93 (1996), 4398–4402. https://doi.org/10.0000/PMID8633078 doi: 10.0000/PMID8633078
![]() |
[2] |
M. A. Nowak, R. M. Charles, Population dynamics of immune responses to persistent viruses, Science, 272 (1996), 74–74. https://doi.org/10.1126/science.272.5258.74 doi: 10.1126/science.272.5258.74
![]() |
[3] |
S. Bonhoeffer, R. M. May, G. M. Shaw, M. A. Nowak, Virus dynamics and drug therapy, Proc. Natl. Acad. Sci., 94 (1997), 6971–6976. https://doi.org/10.1073/pnas.94.13.6971 doi: 10.1073/pnas.94.13.6971
![]() |
[4] |
S. Pankavich, The effects of latent infection on the dynamics of HIV, Differ. Equ. Dyn. Syst., 24 (2016), 281–303. https://doi.org/10.1007/s12591-014-0234-6 doi: 10.1007/s12591-014-0234-6
![]() |
[5] |
B. Buonomo, C. Vargas-De-León, Global stability for an HIV-1 infection model including an eclipse stage of infected cells, J. Math. Anal. Appl., 385 (2012), 709–720. https://doi.org/10.1016/j.jmaa.2011.07.006 doi: 10.1016/j.jmaa.2011.07.006
![]() |
[6] |
A. M. Elaiw, N. H. Alshamrani, Analysis of a within-host HIV/HTLV-I co-infection model with immunity, Virus Res., 295 (2020), 198204. https://doi.org/10.1016/j.virusres.2020.198204 doi: 10.1016/j.virusres.2020.198204
![]() |
[7] |
J. Lang, M. Li, Stable and transient periodic oscillations in a mathematical model for CTL response to HTLV-I infection, J. Math. Biol., 65 (2012), 181–199. https://doi.org/10.1007/s00285-011-0455-z doi: 10.1007/s00285-011-0455-z
![]() |
[8] |
M. Li, H. Shu, Multiple stable periodic oscillations in a mathematical model of CTL response to HTLV-I infection, Bull. Math. Biol., 73 (2011), 1774–1793. https://doi.org/10.1007/s11538-010-9591-7 doi: 10.1007/s11538-010-9591-7
![]() |
[9] | A. M. Elaiw, A. M. Althiabi, M. A. Alghamdi, N. Bellomo, Dynamical behavior of a general HIV-1 infection model with HAART and cellular reservoirs, J. Comput. Anal. Appl., 24 (2018), 728–743. |
[10] |
A. M. Elaiw, E. K. Elnahary, A. A. Raezah, Effect of cellular reservoirs and delays on the global dynamics of HIV, Adv. Differ. Equ., 2018 (2018), 85–121. https://doi.org/10.1186/s13662-018-1523-0 doi: 10.1186/s13662-018-1523-0
![]() |
[11] |
M. Giovanetti, E. Cella, F. Benedetti, B. R. Magalis, V. Fonseca, S. Fabris, et al., SARS-CoV-2 shifting transmission dynamics and hidden reservoirs potentially limit efficacy of public health interventions in Italy, Commun. Biol., 4 (2021), 489–511. https://doi.org/10.1038/s42003-021-02025-0 doi: 10.1038/s42003-021-02025-0
![]() |
[12] |
A. M. Elaiw, T. O. Alade, S. M. Alsulami, Analysis of latent CHIKV dynamics models with general incidence rate and time delays, J. Biol. Dyn., 12 (2018), 700–730. https://doi.org/10.1080/17513758.2018.1503349 doi: 10.1080/17513758.2018.1503349
![]() |
[13] |
A. M. Elaiw, N. H. Alshamrani, Modeling and stability analysis of HIV/HTLV-I co-infection, Int. J. Biomath., 14 (2021), 2150030. https://doi.org/10.1142/S1793524521500303 doi: 10.1142/S1793524521500303
![]() |
[14] |
J. M. Mutua, F. B. Wang, N. K. Vaidya, Effects of periodic intake of drugs of abuse (morphine) on HIV dynamics: Mathematical model and analysis, Math. Biosci., 326 (2020), 108395. https://doi.org/10.1016/j.mbs.2020.108395 doi: 10.1016/j.mbs.2020.108395
![]() |
[15] |
E. Genoyer, C. B. López, The impact of defective viruses on infection and immunity, Annu. Rev. Virol., 6 (2019), 547–566. https://doi.org/10.1146/annurev-virology-092818-015652 doi: 10.1146/annurev-virology-092818-015652
![]() |
[16] |
S. B. Halstead, E. J. O'Rourke, A. C. Allison, Dengue viruses and mononuclear phagocytes. II. Identity of blood and tissue leukocytes supporting in vitro infection, J. Exp. Med., 146 (1977), 218–229. https://doi.org/10.1084/jem.146.1.218 doi: 10.1084/jem.146.1.218
![]() |
[17] | M. A. Nowak, R. M. May, Virus Dynamics, Oxford University Press, New York, 2000. |
[18] |
D. S. Callaway, A. S. Perelson, HIV-1 infection and low steady state viral loads, Bull. Math. Biol., 64 (2002), 29–64. https://doi.org/10.1006/bulm.2001.0266 doi: 10.1006/bulm.2001.0266
![]() |
[19] |
A. Korobeinikov, Global properties of basic virus dynamics models, Bull. Math. Biol., 66 (2004), 879–883. https://doi.org/10.1016/j.bulm.2004.02.001 doi: 10.1016/j.bulm.2004.02.001
![]() |
[20] |
A. M. Elaiw, Global threshold dynamics in humoral immunity viral infection models including an eclipse stage of infected cells, J. KSIAM, 19 (2015), 137–170. https://doi.org/10.12941/jksiam.2015.19.137 doi: 10.12941/jksiam.2015.19.137
![]() |
[21] |
A. Korobeinikov, Global asymptotic properties of virus dynamics models with dose-dependent parasite reproduction and virulence and non-linear incidence rate, Math. Med. Biol., 26 (2009), 225–239. https://doi.org/10.1093/imammb/dqp006 doi: 10.1093/imammb/dqp006
![]() |
[22] |
A. M. Elaiw, S. A. Azoz, Global properties of a class of HIV infection models with Beddington-DeAngelis functional response, Math. Meth. Appl. Sci., 36 (2013), 383–394. https://doi.org/10.1002/mma.2596 doi: 10.1002/mma.2596
![]() |
[23] |
S. Bates, H. Hutson, J. Rebaza, Global stability of Zika virus dynamics, Differ. Equ. Dyn. Syst., 29 (2021), 657–672. https://doi.org/10.1007/s12591-017-0396-0 doi: 10.1007/s12591-017-0396-0
![]() |
[24] |
Y. Cai, K. Wang, W. Wang, Global transmission dynamics of a Zika virus model, Appl. Math. Lett., 92 (2019), 190–195. https://doi.org/10.1016/j.aml.2019.01.015 doi: 10.1016/j.aml.2019.01.015
![]() |
[25] |
S. Wang, D. Zou, Global stability of in-host viral models with humoral immunity and intracellular delays, Appl. Math. Model., 36 (2012), 1313–1322. https://doi.org/10.1016/j.apm.2011.07.086 doi: 10.1016/j.apm.2011.07.086
![]() |
[26] |
X.Wang, X. Song, Global properties of a model of immune effector responses to viral infections, Adv. Complex Syst., 10 (2007), 495–503. https://doi.org/10.1142/S0219525907001252 doi: 10.1142/S0219525907001252
![]() |
[27] | J. P. LaSalle, The Stability of Dynamical Systems, In: Regional Conference Series in Applied Mathematics, SIAM, Philadelphia, 1976. https://doi.org/10.1137/1.9781611970432 |
[28] |
J. Li, X. Xie, Y. Chen, A new way of constructing Lyapunov functions with application to an SI epidemic model, Appl. Math. Lett., 113 (2021), 106777. https://doi.org/10.1016/j.aml.2020.106777 doi: 10.1016/j.aml.2020.106777
![]() |
[29] |
Y. Chen, J. Li, S. Zou, Global dynamics of an epidemic model with relapse and nonlinear incidence, Math. Meth. Appl. Sci., 42 (2019), 1283–1291. https://doi.org/10.1002/mma.5439 doi: 10.1002/mma.5439
![]() |
[30] |
H. Guo, M. Li, Z. Shuai, A graph-theoretic approach to the method of global Lyapunov functions, Proc. Amer. Math. Soc., 136 (2008), 2793–2802. https://doi.org/10.1090/S0002-9939-08-09341-6 doi: 10.1090/S0002-9939-08-09341-6
![]() |
[31] |
J. Li, Y. Xiao, F. Zhang, Y. Yang, An algebraic approach to proving the global stability of a class of epidemic models, Nonlinear Anal. RWA, 13 (2012), 2006–2016. https://doi.org/10.1016/j.nonrwa.2011.12.022 doi: 10.1016/j.nonrwa.2011.12.022
![]() |
[32] |
C. Mccluskey, Lyapunov functions for disease models with immigration of infected hosts, Discrete Contin. Dyn. Syst. Ser. B, 22 (2021), 4479–4491. https://doi.org/10.3934/dcdsb.2020296 doi: 10.3934/dcdsb.2020296
![]() |
[33] |
A. Korobeinikov, Lyapunov functions and global stability for SIR and SIRS epidemiological models with non-linear transmission, Bull. Math. Biol., 68 (2006), 615–626. https://doi.org/10.1007/s11538-005-9037-9 doi: 10.1007/s11538-005-9037-9
![]() |
[34] |
J. Li, Y. Yang, J. Wu, X. Song, Global stability of vaccine-age/staged-structured epidemic models with nonlinear incidence, Electron. J. Qual. Theory Differ. Equ., 18 (2016), 1–17. https://doi.org/10.14232/ejqtde.2016.1.18 doi: 10.14232/ejqtde.2016.1.18
![]() |
[35] |
S. Ottaviano, M. Sensi, S. Sottile, Global stability of SAIRS epidemic models, Nonlinear Anal. RWA, 65 (2022), 103501. https://doi.org/10.1016/j.nonrwa.2021.103501 doi: 10.1016/j.nonrwa.2021.103501
![]() |
[36] | S. Ottaviano, M. Sensi, S. Sottile, Global stability of multi-group SAIRS epidemic models, preprint, preprint, arXiv: 2202.02993. |
[37] |
Z. Shuai, P. van den Driessche, Global stability of infectious disease models using Lyapunov functions, SIAM J. Appl. Math., 73 (2013), 1513–1532. https://doi.org/10.1137/120876642 doi: 10.1137/120876642
![]() |
[38] |
A. A. Lashari, G. Zaman, Global dynamics of vector-borne diseases with horizontal transmission in host population, Comput. Math. Appl., 61 (2011), 745–754. https://doi.org/10.1016/j.camwa.2010.12.018 doi: 10.1016/j.camwa.2010.12.018
![]() |
[39] |
S. S. Nadim, I. Ghosh, J. Chattopadhyay, Global dynamics of a vector-borne disease model with two transmission routes, Internat. J. Bifur. Chaos Appl. Sci. Engrg., 30 (2020), 2050083, 23 pp. https://doi.org/10.1142/S0218127420500832 doi: 10.1142/S0218127420500832
![]() |
[40] |
M. Ozair, Q. Din, T. Hussain, A. Z. Awan, Qualitative behavior of vector-borne disease model, J. Nonlinear Sci. Appl., 9 (2016), 1382–1395. http://dx.doi.org/10.22436/jnsa.009.03.62 doi: 10.22436/jnsa.009.03.62
![]() |
[41] |
C. Palmer, E. Landguth, E. Stone, T. Johnson, The dynamics of vector-borne relapsing diseases, Math. Biosci., 297 (2018), 32–42. https://doi.org/10.1016/j.mbs.2018.01.001 doi: 10.1016/j.mbs.2018.01.001
![]() |
[42] |
A. Traoré, Analysis of a vector-borne disease model with human and vectors immigration, J. Appl. Math. Comput., 64 (2020), 411–428. https://doi.org/10.1007/s12190-020-01361-4 doi: 10.1007/s12190-020-01361-4
![]() |
[43] |
P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
![]() |
[44] | F. F. Verhulst, Nonlinear Differential Equations and Dynamical Systems, Springer-Verlag, Berlin, Heidelberg, 1996. https://doi.org/10.1007/978-3-642-61453-8 |
[45] |
J. Li, Y. Yang, Y. Zhou, Global stability of an epidemic model with latent stage and vaccination, Nonlinear Anal. RWA, 12 (2011), 2163–2173. https://doi.org/10.1016/j.nonrwa.2010.12.030 doi: 10.1016/j.nonrwa.2010.12.030
![]() |
[46] | H. Guo, M. Y. Li, Z. Shuai, Global stability of the endemic equilibrium of multigroup SIR epidemic models, Can. Appl. Math. Q., 14 (2006), 259–284. |