A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator

1.   Introduction

In this note we consider the following magnetic, quantum Euler equations:

in the unknowns $(\rho, J):{{\mathbb R}}_t\times{{\mathbb R}}^3\to {{\mathbb R}}_{\geqslant 0}\times{{\mathbb R}}^3$, 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\wedge B$ is the magnetic force associated to the field $B$ (we are neglecting here the self-induced electromagnetic field). We assume $B = \nabla\wedge A$ for some real-valued, divergence-free potential $A$. Finally, $P(\rho)$ is a barotropic pressure term, which assumes the form $P(\rho) = \frac{\gamma-1}{\gamma}\rho^{\gamma}$, with $\gamma > 1$. From now on, we focus without loss of generality only on forward-in-time solutions (i.e., defined for $t\geqslant 0$) and we denote the initial data at time $t = 0$ by $(\rho_0, J_0)$.

1.1.   Physical background and motivations

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\wedge 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.

1.2.   Rigorous setup and the notion of finite-energy, weak solution

System (1.1) admits two quantities that are formally preserved along the flow: the total charge $\|\rho\|_{L^1} = \|\sqrt{\rho}\|_{L^2}^2$ and the total energy

where $\Lambda: = \frac{J}{\sqrt{\rho}}$ and $f(\rho): = \frac{\rho^\gamma}{\gamma}$ is the internal energy density. Observe that, for $\gamma\leqslant 3$, we have $f(\rho)\in L^1({{\mathbb R}}^3)$ whenever $\sqrt{\rho}\in H^1({{\mathbb R}}^3)$. We restrict to the subcritical regime $\gamma < 3$, while the critical case $\gamma = 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

We work in a quasi-irrotational framework. More precisely, for almost every time we require the generalized irrotationality condition

to be satisfied as an identity in $\mathcal{D}'({{\mathbb R}}^3)$. In order to clarify the meaning of the above condition, consider a smooth solution $(\rho, J)$ to System (1.1), for which we can define the velocity field $v$ through the relation $J = \rho v$. A direct computation then shows that System (1.3) is equivalent to the identity

which implies in particular that the relative velocity field $v_A: = v+A$ is irrotational outside vacuum, namely $\nabla\wedge v_A = 0$ on $\{\rho > 0\}$. On the other side, condition (1.3) is compatible with the occurrence of quantum vortices, namely solutions such that $v_A$ 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\wedge 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,

which guarantees that $J\wedge B$ is locally integrable only when $B\in L^3_{\mathrm{loc}}({{\mathbb R}}^3)$. 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\in L^{2+}({{\mathbb R}}^3)$.

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

while the quantum term can be written in divergence form as

Definition 1.1. Let $0 < T\leqslant\infty$, $\rho_0, J_0\in L^1_{loc}({{\mathbb R}}^3)$. We say that $(\rho, J)$ is a finite-energy weak solution to System (1.1) on the space-time slab $[0, T)\times{{\mathbb R}}^3$, with initial data $(\rho_0, J_0)$, if the following conditions are satisfied:

(i) (finite energy hydrodynamic state) there exists a pair $(\sqrt{\rho}, \Lambda)$, with $\sqrt{\rho}\in L_{\mathrm{loc}}^{\infty}(0, T; H^1({{\mathbb R}}^3))$ and $\Lambda\in L_{\mathrm{loc}}^{\infty}(0, T;L^2({{\mathbb R}}^3))$, such that $\rho = (\sqrt{\rho})^2$, $J = \sqrt{\rho}\Lambda$;

(ii) (continuity and momentum equations - weak formulations) for every $\eta\in\mathcal{C}^{\infty}_c([0, T)\times{{\mathbb R}}^3)$ and $\zeta\in\mathcal{C}^{\infty}_c([0, T)\times{{\mathbb R}}^3;{{\mathbb R}}^3)$,

(iv) (well-defined magnetic force) $J\wedge B\in L^1_{\mathrm{loc}}((0, T)\times{{\mathbb R}}^3)$;

(v) (generalized irrotationality condition) for a. e. $t\in(0, T)$, the generalized irrotationality condition (1.3) holds true in $\mathcal{D}'({{\mathbb R}}^3)$.

When $T = \infty$ we say that the solution is global.

2.   Connection with magnetic nonlinear Schrödinger (NLS) and main result

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 $\psi$ satisfies the magnetic NLS equation

where $\Delta_A: = \nabla_A\cdot\nabla_A: = (\nabla- \mathrm{i} A)\cdot(\nabla- \mathrm{i} A) = \Delta-2 \mathrm{i} A\cdot\nabla -|A|^2$ is the magnetic Laplacian, then its Madelung transform

formally solves System (1.1). Suppose indeed that we can write $\psi = \sqrt{\rho}e^{ \mathrm{i}\theta}$ so that $J = \rho(\nabla\theta-A)$ for smooth functions $\rho, \theta$, with $\rho > 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\wedge 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

for which a satisfactory well-posedness theory at the finite energy level (i.e., $H^1$-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:

As anticipated, in terms of the magnetic field $B = \nabla\wedge A$ we are just above the finite-energy level $B\in L^2({{\mathbb R}}^3)$. 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

coincides (with equivalence of norms) with the classical Sobolev space $H^1({{\mathbb R}}^3)$. Since the energy of the magnetic NLS (2.1) is given by

we then deduce (in the sub-critical regime $\gamma < 3$) that $H^1({{\mathbb R}}^3)$ is the energy space for Eq (2.1).

Concerning the initial data $(\rho_0, J_0)$, we impose that they are the Madelung transform of a finite energy wave function:

The above condition is natural from the quantum mechanics perspective, as it is equivalent ^[30] to assume that the (relative) initial velocity field $v_{0, A}$ has quantized circulations around vacuum, namely

for any closed, smooth path in ${\rho > 0}$.

We are now able to state our main result.

Theorem 2.1. Assume that the magnetic potential $A$ and the initial data $(\rho_0, J_0)$ 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 $H^1({{\mathbb R}}^3)$ 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^{\infty}$-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 $2 \mathrm{i} A\cdot \nabla$ 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\lesssim B$ 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})$.

3.   Proof of the main result

Our starting point is the existence of a global, finite energy weak solution

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

of polar factors of $\psi$. Given $\varphi\in P(\psi)$, let us define

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

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

Using assumptions (2.4) on $A$ we get

Moreover, since $\gamma < 3$, a direct application of Sobolev embedding, Gagliardo-Niremberg and Hölder inequalities yields

for some $\varepsilon = \varepsilon(\sigma)\in(0, 1)$. Since $\|\psi\|_{L_T^{\infty}H^1}\lesssim_T 1$, we then obtain

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

(ⅱ) 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$,

which implies in particular $J\wedge B\in L^1_{\mathrm{loc}}({{\mathbb R}}\times{{\mathbb R}}^3)$.

(ⅳ) We start by observing that

which yields

On the other hand, given a polar factor $\varphi$ of $\psi$, we have

Moreover, we have the chain of identities

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.

4.   Conclusions

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

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

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.

Use of AI tools declaration

The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

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".

Conflict of interest

The author declares there is no conflict of interest.