In Japan, major and minor bimodal seasonal patterns of varicella have been observed. To investigate the underlying mechanisms of seasonality, we evaluated the effects of the school term and temperature on the incidence of varicella in Japan. We analyzed epidemiological, demographic and climate datasets of seven prefectures in Japan. We fitted a generalized linear model to the number of varicella notifications from 2000 to 2009 and quantified the transmission rates as well as the force of infection, by prefecture. To evaluate the effect of annual variation in temperature on the rate of transmission, we assumed a threshold temperature value. In northern Japan, which has large annual temperature variations, a bimodal pattern in the epidemic curve was observed, reflecting the large deviation in average weekly temperature from the threshold value. This bimodal pattern was diminished with southward prefectures, gradually shifting to a unimodal pattern in the epidemic curve, with little temperature deviation from the threshold. The transmission rate and force of infection, considering the school term and temperature deviation from the threshold, exhibited similar seasonal patterns, with a bimodal pattern in the north and a unimodal pattern in the south. Our findings suggest the existence of preferable temperatures for varicella transmission and an interactive effect of the school term and temperature. Investigating the potential impact of temperature elevation that could reshape the epidemic pattern of varicella to become unimodal, even in the northern part of Japan, is required.
Citation: Ayako Suzuki, Hiroshi Nishiura. Seasonal transmission dynamics of varicella in Japan: The role of temperature and school holidays[J]. Mathematical Biosciences and Engineering, 2023, 20(2): 4069-4081. doi: 10.3934/mbe.2023190
[1] | Mustafa Mudhesh, Hasanen A. Hammad, Eskandar Ameer, Muhammad Arshad, Fahd Jarad . Novel results on fixed-point methodologies for hybrid contraction mappings in Mb-metric spaces with an application. AIMS Mathematics, 2023, 8(1): 1530-1549. doi: 10.3934/math.2023077 |
[2] | Muhammad Tariq, Muhammad Arshad, Mujahid Abbas, Eskandar Ameer, Saber Mansour, Hassen Aydi . A relation theoretic m-metric fixed point algorithm and related applications. AIMS Mathematics, 2023, 8(8): 19504-19525. doi: 10.3934/math.2023995 |
[3] | Huaping Huang, Bessem Samet . Two fixed point theorems in complete metric spaces. AIMS Mathematics, 2024, 9(11): 30612-30637. doi: 10.3934/math.20241478 |
[4] | Hanadi Zahed, Zhenhua Ma, Jamshaid Ahmad . On fixed point results in F-metric spaces with applications. AIMS Mathematics, 2023, 8(7): 16887-16905. doi: 10.3934/math.2023863 |
[5] | Umar Ishtiaq, Fahad Jahangeer, Doha A. Kattan, Manuel De la Sen . Generalized common best proximity point results in fuzzy multiplicative metric spaces. AIMS Mathematics, 2023, 8(11): 25454-25476. doi: 10.3934/math.20231299 |
[6] | Leyla Sağ Dönmez, Abdurrahman Büyükkaya, Mahpeyker Öztürk . Fixed-point results via αji-(DC(PˆE))-contractions in partial ♭-metric spaces. AIMS Mathematics, 2023, 8(10): 23674-23706. doi: 10.3934/math.20231204 |
[7] | Jamshaid Ahmad, Abdullah Shoaib, Irshad Ayoob, Nabil Mlaiki . Common fixed points for (κGm)-contractions with applications. AIMS Mathematics, 2024, 9(6): 15949-15965. doi: 10.3934/math.2024772 |
[8] | Muhammad Waseem Asghar, Mujahid Abbas, Cyril Dennis Enyi, McSylvester Ejighikeme Omaba . Iterative approximation of fixed points of generalized αm-nonexpansive mappings in modular spaces. AIMS Mathematics, 2023, 8(11): 26922-26944. doi: 10.3934/math.20231378 |
[9] | Tahair Rasham, Muhammad Nazam, Hassen Aydi, Abdullah Shoaib, Choonkil Park, Jung Rye Lee . Hybrid pair of multivalued mappings in modular-like metric spaces and applications. AIMS Mathematics, 2022, 7(6): 10582-10595. doi: 10.3934/math.2022590 |
[10] | Amer Hassan Albargi, Jamshaid Ahmad . Fixed point results of fuzzy mappings with applications. AIMS Mathematics, 2023, 8(5): 11572-11588. doi: 10.3934/math.2023586 |
In Japan, major and minor bimodal seasonal patterns of varicella have been observed. To investigate the underlying mechanisms of seasonality, we evaluated the effects of the school term and temperature on the incidence of varicella in Japan. We analyzed epidemiological, demographic and climate datasets of seven prefectures in Japan. We fitted a generalized linear model to the number of varicella notifications from 2000 to 2009 and quantified the transmission rates as well as the force of infection, by prefecture. To evaluate the effect of annual variation in temperature on the rate of transmission, we assumed a threshold temperature value. In northern Japan, which has large annual temperature variations, a bimodal pattern in the epidemic curve was observed, reflecting the large deviation in average weekly temperature from the threshold value. This bimodal pattern was diminished with southward prefectures, gradually shifting to a unimodal pattern in the epidemic curve, with little temperature deviation from the threshold. The transmission rate and force of infection, considering the school term and temperature deviation from the threshold, exhibited similar seasonal patterns, with a bimodal pattern in the north and a unimodal pattern in the south. Our findings suggest the existence of preferable temperatures for varicella transmission and an interactive effect of the school term and temperature. Investigating the potential impact of temperature elevation that could reshape the epidemic pattern of varicella to become unimodal, even in the northern part of Japan, is required.
In this paper, we are concerned with the existence of positive solutions for the following Kirchhoff equation in RN:
−(a+b∫RN|∇u|2dx)Δu+λu=μ|u|q−2u+|u|p−2u,(K) |
having prescribed mass
∫RN|u|2dx=c2, | (1.1) |
where c>0 is prescribed, N⩾4, μ∈R, a>0 and b>0 are positive parameters, 2<q<p⩽2∗, 2∗ is the Sobolev critical exponent, 2∗=2NN−2 if N⩾3, 2∗=+∞ if N=1,2.
Problem (K) is analogous to the stationary case of the following equation:
utt−(a+b∫Ω|∇u|2dx)Δu=f(x,u),(x,t)∈Ω×R+, |
where Ω⊂RN is a smooth domain, u stands for the displacement, f∈C(Ω×R,R) is the external force, a is the initial tension and b is related to the intrinsic properties of the string. Such a hyperbolic equation is a general version of the following equations
ρ∂2u∂t2−(P0h+E2L∫L0|∂u∂x|2dx)∂2u∂x2=0, |
proposed by Kirchhoff [1] as an extension of the classical D'Alembert's wave equation for transversal oscillations of a stretched string, particularly, taking into account the subsequent change in string length during the oscillations. In addition, nonlocal problems similar to (K) also arise in biological systems that can be applied to describe the growth and movement of a specific species[2]. In this situation, u may describe a process that depends on the average of itself, for example, the population density.
After the pioneering work of [3], problem (K) with fixed λ, or even with an additional external and fixed potential V(x) has been extensively studied in the last ten years, see e.g., [2,4,5,6,7] and the references therein. Compared with Schrödinger type problem, the main difficulty for Kirchhoff type problem is how to obtain the compactness of the Palais-Smale sequence. Tang-Cheng [5] proposed a new approach to recover the compactness for the Palais-Smale sequences. Tang-Chen [6] proposed a new approach to recover the compactness for the minimizing sequences.
In this situation, it is well known that solutions of (K) are the critical points of the associated energy functional
Φλ(u):=12∫RN(a|∇u|2+λ|u|2)dx+b4(∫RN|∇u|2dx)2−μq∫RN|u|qdx−1p∫RN|u|pdx. |
In the present paper, motivated by the fact that physicists are more interested in solutions satisfying the prescribed mass (1.1) which are called normalized solutions, we search for solutions of (K) having a prescribed L2-norm. This type of solutions can be obtained as a constrained critical point of the following functional
I(u)=a2∫RN|∇u|2dx+b4(∫RN|∇u|2dx)2−μq∫RN|u|qdx−1p∫RN|u|pdx, | (1.2) |
on the L2-sphere
S(c)={u∈H1(RN):‖u‖L2(RN)=c}, |
where the Hilbert space H1(RN) is defined as
H1(RN):={u∈L2(RN):∫RN|∇u|2dx<∞}, |
with the inner product and norm given by
(u,v)H1(RN)=∫RN(∇u∇v+uv)dx,‖u‖H1(RN)=(‖∇u‖2L2(RN)+‖u‖2L2(RN))12. |
When a=1,b=μ=0, problem (K) is reduced to the classical Schrödinger equation with pure power nonlinearity:
−Δu+λu=|u|p−2u,x∈RN, | (1.3) |
where N⩾1, p∈(2,2∗). It is well known that the L2-critical exponent
ˉp=2+4N, |
plays an important role in dealing with the Schrödinger equation. For the pure power nonlinearity |u|p−2u satisfying L2-subcritical growth, i.e., 2<p<ˉp, then the energy functional associated to (1.3) is bounded from below. Cazenave-Lions [8] and Shibata [9] proved the existence of normalized solution via minimizing methods. We also refer the reader to the references [10,11]. While for L2-supercritical and Sobolev subcritical case ˉp<p<2∗, the energy functional associated to (1.3) cannot be minimized on the L2-sphere S(c). Jeanjean [12] proved for the first time the existence of a normalized solution in this case. When a=1,b=0 and μ≠0, problem (K) is reduced to the Schrödinger equation with combined power nonlinearities:
−Δu+λu=μ|u|q−2u+|u|p−2u,x∈RN, | (1.4) |
where N⩾1, μ∈R, the nonlinearity satisfies the mixed growth, i.e., 2<q⩽ˉp⩽p⩽2∗. Recently, a very complete analysis of the various cases that may occur for (1.4), depending on the values of p and q, has been studied exhaustively in [13,14]. The paper [13] deals with the Sobolev subcritical case and [14] with Sobolev critical case.
It is well known that two exponents play a special role for existence and properties of the solutions for the Kirchhoff equation, in addition to the energy (Sobolev) critical exponent
2∗=2NN−2, |
we have the mass (L2-) critical one
2+8N. |
The mass and energy critical exponents define the thresholds for the existence of global minima and H1-solution, respectively. It follows from
2+8N<2NN−2 |
that N⩽3. For this reason possibly, there are many results on normalized solution of Kirchhoff problem in low dimensions N⩽3, we refer the readers to [15,16,17,18,19] and references therein. Howover, to the best of our knowledge, the Kirchhoff type problems with prescribed L2-norm in high dimensions N⩾4 remains unsolved in the existing literature. This motivates us to study this kind of problem in the present paper for the case of N=4,6,8. When N⩾4, since 2∗⩽4, it is difficult to study the geometric structure of the energy functional. Also, it is not easy to prove the boundedness and convergence of the Palais-Smale sequence.
Particularly, in dimension N=4, there is a special phenomenon for Kirchhoff problems that will never happen for the Schrödinger equation, i.e., the mass critical exponent 2+8N is the same as the energy critical exponent 2∗=2NN−2. We call this special phenomenon the new doubly critical case, which is first considered to happen for the Kirchhoff equation in this paper. For other types of double critical exponents, we refer the readers to [20,21] and references therein. Moreover, it is worth mentioning that when N⩾4, if we consider the mass (L2-) supercritical case, the problem will no longer have a variational structure owing to the fact that the mass critical exponent is larger than the energy critical exponent. Thus, there will be no mass supercritical and mixed cases.
Motivated by the works mentioned above, our contribution in this paper is to find appropriate transformation to get the following equivalent system of (K) in different dimension spaces with respect to (v,β)∈RN×R+:
{−Δv+λv=μ|v|q−2v+|v|p−2v,β−a−bβN−22∫RN|∇v|2dx=0,∫RN|v|2dx=˜β−N2Nc2,(PN) |
where
˜βN={a1−d,N=4,1±√1−4ad2d,N=6,2√3dcos(θ+2jπ3)(j=0,2),N=8, | (1.5) |
θ=13arccos(−32√3da)∈(π6,π3), and
d=b‖∇v‖2L2(RN). | (1.6) |
We take two steps to get the existence of solutions for system (PN). First, we solve the equation
{−Δv+λv=μ|v|q−2v+|v|p−2v,∫RN|v|2dx=˜β−N2Nc2, | (1.7) |
for the special given ˜βN. Then, when v is known, coincidentally, ˜βN can solve the algebra equation
β−a−bβN−22∫RN|∇v|2dx=0. | (1.8) |
Therefore, Some known difficulties can be overcome by establishing equivalent results, which is the key ingredient to prove the nonexistence, existence and multiplicity of normalized solutions for the Kirchhoff equation (K).
Theorem 1.1. Let N=4,6,8. Then equation (K) has at least one positive, radially symmetric normalized solution u∈S(c) if and only if the system (PN) has at least one solution (v,˜βN)∈S(˜β−N4Nc)×R+ such that u is positive and radially symmetric, where ˜βN is defined by (1.5).
To state our another main results, we introduce some frequently used constants. For p∈(2,2∗), we use the notation
γp:=N(p−2)2p, |
and we denote by CN,p the best constant in the Gagliardo-Nirenberg inequality[22]
‖u‖Lp(RN)⩽CN,p‖∇u‖γpL2(RN)‖u‖1−γpL2(RN),∀u∈H1(RN). |
In the special case p=ˉp, we denote by ˉCN(β)=(ˉpβ2CˉpN,ˉp)N4 the critical mass. For 2<q⩽ˉp<p<2∗ and μ>0, we denote:
μ∗(β)=(p(2−γqq)2CpN,p(γpp−γqq))2−γqqγpp−2(q(γpp−2)2CqN,q(γpp−γqq))((β−N4c)(1−γq)q)((β−N4c)(1−γp)p)2−γqqγpp−2,μ∗∗(β)=ˉpβ2CˉpN,ˉpc4N. |
For 2<q⩽ˉp<p<2∗ and μ<0, we denote:
μ∗(β)=(1−γpCqN,q(γp−γq))(1γpCpN,p)2−γqqγpp−2((β−N4c)(1−γq)q)((β−N4c)(1−γp)p)2−γqqγpp−2. |
For 2<q<ˉp<p=2∗ and μ>0, we denote:
ˆμ(β)=(2∗S2∗2(2−γqq)2(2∗−γqq))2−γqq2∗−2(q(2∗−2)2CqN,q(2∗−γqq)(β−N4c)(1−γq)q),ˇμ(β)=22∗NγqCqN,q(2∗−γqq)(β−N4c)(1−γq)q(γqqSN22−γqq)2−γqq2. |
For q=ˉp<p=2∗ and μ>0, we denote:
ˉμ(β)=ˉp2CˉpN,ˉpc4N(β−N4c)(1−γq)q. |
For ˉp<q<p=2∗ and μ>0, we denote:
μ∗∗(β)=+∞ if N=4, and μ∗∗(β)=SN4(1−γq)qγq(β−N4c)(1−γq)q if N⩾5, |
where the optimal Sobolev embedding constant S[23] is given by
S=infu∈D1,2(RN)∖{0}‖∇u‖2L2(RN)‖u‖2L2∗(RN), |
where D1,2(RN) denotes the completion of C∞0(RN) with respect to the semi-norm ‖u‖D1,2(RN):=‖∇u‖L2(RN).
Applying Theorem 1.1, we can prove the following results. Firstly, in the Sobolev subcritical case 2<q<p=ˉp or 2<q⩽ˉp<p<2∗, we have:
Theorem 1.2. Suppose that a,b>0 and 2<q<p<2∗ hold. Let ˜βN and d be defined by (1.5) and (1.6), respectively. If one of the following conditions holds:
(1)2<q<p=ˉp, 0<c<ˉCN(˜βN) and μ>0;
(2)2<q⩽ˉp<p<2∗, c>0 and 0<μ<min{μ∗(˜βN),μ∗∗(~βN)};
(3)2<q⩽ˉp<p<2∗, c>0 and −μ∗(˜βN)<μ<0;
then we have the following conclusions:
(i) for N=4, equation (K) has a positive, radially symmetric solution u∈S(c) if d∈(0,1), and has no nontrivial solution if d∈[1,+∞);
(ii) for N=6,8, equation (K) has two positive, radially symmetric solutions u1, u2∈S(c) if d∈(0,dN), has one positive, radially symmetric solution u∈S(c) if d=dN, and has no nontrivial solution if d∈(dN,+∞), where dN=14a when N=6 and dN=427a2 when N=8.
Next, in the Sobolev critical case 2<q<p=2∗, we have the following result:
Theorem 1.3. Suppose that a,b,c>0 and 2<q<p=2∗ hold. Let ˜βN and d be defined by (1.5) and (1.6), respectively. If one of the following conditions holds:
(1)2<q<ˉp<p=2∗ and 0<μ<min{ˇμ(˜βN),ˆμ(˜βN)};
(2)ˉp=q<p=2∗ and 0<μ<ˉμ(˜βN);
(3)ˉp<q<p=2∗ and 0<μ<μ∗∗(˜βN);
then we have the following conclusions:
(i) for N=4, equation (K) has a positive, radially symmetric solution u∈S(c) if d∈(0,1), and has no nontrivial solution if d∈[1,+∞);
(ii) for N=6,8, equation (K) has two positive, radially symmetric solutions u1, u2∈S(c) if d∈(0,dN), has one positive, radially symmetric solution u∈S(c) if d=dN, and has no nontrivial solution if d∈(dN,+∞), where dN=14a when N=6 and dN=427a2 when N=8.
Remark 1.1. Theorem 1.3 covers the doubly critical case that the mass critical exponent is the same as the energy critical exponent i.e., in dimension N=4,
p=2∗=2+8N. |
Remark 1.2. As far as we know, Theorems 1.2 and 1.3 seem to be the first results on the existence and multiplicity of normalized solutions for equation (K) in high dimensions N=4,6,8. We also remark that it is interesting to know whether (K) has a normalized solution for aN−42b small in more general dimensions N⩾4. However, we can not solve this question now. So far, we just know that (K) has at least a positive normalized solution for aN−42b small, and has no nontrivial solution for aN−42b large in the special dimensions N=4,6,8. Moreover, it is reasonable to regard a and b as parameters in (K) since, as already noted before, a is the initial tension and b is related to the intrinsic properties of the string.
Finally, in the doubly Sobolev critical case q=p=2∗, based on the Pohožaev identity, we have the following nonexistence result, which is not related to the equivalent results but interesting in itself.
Theorem 1.4. Suppose that p=q=2∗ and N=4 holds, then equation (K) has no positive solutions u∈S(c) for any μ>0.
This paper is organized as follows: In Section 2, we present some preliminaries. In Section 3, we show some results of (PN). In Section 4, we give the proofs of Theorems 1.1–1.4.
In this section, we introduce some preliminaries that will be used to prove Theorems 1.2–1.4.
Lemma 2.1. Let μ>0, N⩽4 and u∈S(c) be a solution of (K). If u⩾0 and u≢, then \lambda > 0 .
Proof. Arguing indirectly, suppose that \lambda\leqslant 0 . It follows from u\geqslant 0 that
-\left(a+b \int_{\mathbb{R}^{N}}|\nabla u|^2\mathrm{d}x\right) \Delta u = -\lambda u+\mu |u|^{q-2} u+|u|^{p-2} u\geqslant0, |
Thus -\Delta u\geqslant 0 . Moreover, from standard regularity results we have that u\in C^2(\mathbb{R}^N) . By the Liouville's result [24,Lemma A.2], we obtain that u\equiv0 . This contradicts with u\not\equiv0 . Thus \lambda > 0 .
Lemma 2.2. Let a, b > 0 , p, q\in(2, 2^*] and \mu, \lambda\in\mathbb{R} . If u\in H^1(\mathbb{R}^N) is a weak solution of the Kirchhoff equation (\mathcal{K}), then the following Pohožaev identity holds:
\begin{equation} P(u) = a\|\nabla u\|_{L^2(\mathbb{R}^N)}^2+b\|\nabla u\|_{L^2(\mathbb{R}^N)}^4-\mu \gamma_{q}\|u\|^q_{L^q(\mathbb{R}^N)}-\gamma_{p}\|u\|^p_{L^p(\mathbb{R}^N)} = 0. \end{equation} | (2.1) |
Proof. Let M = \|\nabla u\|_{L^2(\mathbb{R}^N)}^2 , then u\in H^1(\mathbb{R}^N) is a weak solution of the following problem
-\Delta u = \frac{1}{a+bM}(\mu |u|^{q-2} u+|u|^{p-2} u-\lambda u). |
Then the elliptic regularity theory implies that u\in C^2(\mathbb{R}^N) . Combining
a\|\nabla u\|_{L^2(\mathbb{R}^N)}^2+b\|\nabla u\|_{L^2(\mathbb{R}^N)}^4 = \mu \|u\|^q_{L^q(\mathbb{R}^N)}+\|u\|^p_{L^p(\mathbb{R}^N)}-\lambda \|u\|_{L^2(\mathbb{R}^N)}^2 |
with
a\|\nabla u\|_{L^2(\mathbb{R}^N)}^2+b\|\nabla u\|_{L^2(\mathbb{R}^N)}^4 = \frac{2N}{N-2}\left[\frac{\mu}{q}\|u\|^q_{L^q(\mathbb{R}^N)}+\frac{1}{p}\|u\|^p_{L^p(\mathbb{R}^N)}-\frac{\lambda}{2} \|u\|_{L^2(\mathbb{R}^N)}^2\right], |
we conclude that (2.1) holds.
Lemma 2.3. For m, n\in \mathbb{R} and (\frac{n}{2})^2+(\frac{m}{3})^3 < 0 , then the following cubic
\begin{equation} x^3+mx+n = 0, \end{equation} | (2.2) |
has three distinct real roots, which can be written in the form
\begin{equation} x_j = 2\sqrt[3]{r}\cos\left(\theta+\frac{2j\pi}{3}\right)(j = 0,1,2), \end{equation} | (2.3) |
where r = \sqrt{-(\frac{m}{3})^3} and \theta = \frac{1}{3}\arccos(-\frac{n}{2r}).
Proof. This lemma is a special case of Cardano's formulas, and for the convenience of the reader, we give a simple proof. Set
\begin{equation} m = -3yz\text{ and } n = -y^3-z^3, \end{equation} | (2.4) |
it follows from the Binomial Theorem that x = y+z always satisfies (2.2). From (2.4), we can obtain that y^3 and z^3 are the two roots of the quadratic
X^2+nX-\frac{m^3}{27} = 0. |
Furthermore, by (\frac{n}{2})^2+(\frac{m}{3})^3 < 0 , without loss of generality, we may assume that
y^3 = -\frac{n}{2}+i\sqrt{-(\frac{n}{2})^2-(\frac{m}{3})^3} = r(\cos(3\theta)+i\sin(3\theta)), |
and
z^3 = -\frac{n}{2}-i\sqrt{-(\frac{n}{2})^2-(\frac{m}{3})^3} = r(\cos(3\theta)-i\sin(3\theta)), |
where r = \sqrt{-(\frac{m}{3})^3} and \theta = \frac{1}{3}\arccos(-\frac{n}{2r})\in (0, \frac{\pi}{3}). It follows that
y = \sqrt[3]{r}\left(\cos(\theta+\frac{2j\pi}{3})+i\sin(\theta+\frac{2j\pi}{3})\right)(j = 0,1,2), |
and
z = \sqrt[3]{r}\left(\cos(\theta+\frac{2j\pi}{3})-i\sin(\theta+\frac{2j\pi}{3})\right)(j = 0,1,2). |
Note that x = y+z , we then conclude that (2.2) has three distinct real roots, which can be expressed as in (2.3).
We devote to solve system \left(\mathcal{P}_{N}\right) in this section. We will prove the following results which will be used to conclude Theorems 1.2 and 1.3.
Firstly, in the Sobolev subcritical case 2 < q < p = \bar{p} or 2 < q\leqslant\bar{p} < p < 2^* , we have:
Proposition 3.1. Suppose that a > 0 and 2 < q < p < 2^* hold. Let \tilde{\beta}_N and d be defined by (1.5) and (1.6), respectively. If one of the following conditions holds:
(1) \; 2 < q < p = \bar{p} , 0 < c < \bar{C}_N(\tilde{\beta}_N) and \mu > 0 ;
(2) \; 2 < q\le\bar{p} < p < 2^* , c > 0 and 0 < \mu < \min\{\mu^*(\tilde{\beta}_N), \mu^{**}(\tilde{\beta}_N)\} ;
(3) \; 2 < q\le \bar{p} < p < 2^* , c > 0 and -\mu_{*}(\tilde{\beta}_N) < \mu < 0 ;
then we can obtain the following results:
(i) for N = 4 , system \left(\mathcal{P}_{N}\right) has a positive, radially symmetric solution (v, \frac{a}{1-d}) if d\in(0, 1) , and has no nontrivial solution if d\in[1, +\infty) ;
(ii) for N = 6, 8 , system \left(\mathcal{P}_{N}\right) has two positive, radially symmetric solutions (v, \tilde{\beta}_{N}^{1}) and (v, \tilde{\beta}_{N}^{2}) if d\in(0, d_N) , has one positive, radially symmetric solution (v, \frac{1}{2d}) if d = d_N , and has no nontrivial solution if d\in (d_N, +\infty) , where d_N = \frac{1}{4a} when N = 6 and d_N = \frac{4}{27a^2} when N = 8 .
Next, in the Sobolev critical case 2 < q < p = 2^* , we have the following result.
Proposition 3.2. Suppose that c > 0 and 2 < q < p = 2^* hold. Let \tilde{\beta}_N and d be defined by (1.5) and (1.6), respectively. If one of the following conditions holds:
(1) \; 2 < q < \bar{p} < p = 2^* and 0 < \mu < \min\{\check{\mu}(\tilde{\beta}_N), \hat{\mu}(\tilde{\beta}_N)\} ;
(2) \; \bar{p} = q < p = 2^* and 0 < \mu < \bar{\mu}(\tilde{\beta}_N) ;
(3) \; \bar{p} < q < p = 2^* and 0 < \mu < \mu_{**}(\tilde{\beta}_N) ;
then we have the following results:
(i) for N = 4 , system \left(\mathcal{P}_{N}\right) has a positive, radially symmetric solution (v, \frac{a}{1-d}) if d\in(0, 1) , and has no nontrivial solution if d\in[1, +\infty) .
(ii) for N = 6, 8 , system \left(\mathcal{P}_{N}\right) has two positive, radially symmetric solutions (v, \tilde{\beta}_{N}^{1}) and (v, \tilde{\beta}_{N}^{2}) if d\in(0, d_N) , has one positive, radially symmetric solution (v, \frac{1}{2d}) if d = d_N , and has no nontrivial solution if d\in (d_N, +\infty) , where d_N = \frac{1}{4a} when N = 6 and d_N = \frac{4}{27a^2} when N = 8 .
In order to obtain conclusions of Proposition 3.1 and 2, we first recall some results about existence of normalized solutions for the equation (1.7).
Theorem 3.1. [13] Suppose that c, \tilde{\beta}_N > 0 , N\geqslant 4 and 2 < q < p < 2^* hold. If one of the following conditions holds:
(1) \; 2 < q < p = \bar{p} , 0 < c < \bar{C}_N(\tilde{\beta}_N) and \mu > 0 ;
(2) \; 2 < q\leqslant\bar{p} < p < 2^* , c > 0 and 0 < \mu < \min\{\mu^*(\tilde{\beta}_N), \mu^{**}(\tilde{\beta}_N)\} ;
(3) \; 2 < q\leqslant \bar{p} < p < 2^* , c > 0 and -\mu_{*}(\tilde{\beta}_N) < \mu < 0 ;
then Eq (1.7) has at least one positive, radially symmetric solution.
Theorem 3.2. [14,Theorem 1.1] Suppose that c, \tilde{\beta}_N > 0 , N\geqslant 4 and 2 < q < p = 2^* hold. If one of the following conditions holds:
(1) \; 2 < q < \bar{p} < p = 2^* and 0 < \mu < \min\{\check{\mu}(\tilde{\beta}_N), \hat{\mu}(\tilde{\beta}_N)\} ;
(2) \; \bar{p} = q < p = 2^* and 0 < \mu < \bar{\mu}(\tilde{\beta}_N) ;
(3) \; \bar{p} < q < p = 2^* and 0 < \mu < \mu_{**}(\tilde{\beta}_N) ;
then Eq (1.7) has at least one positive, radially symmetric solution.
Remark 3.1.
(i) Theorem 3.1 can be obtained from [13] with some minor modifications. Specifically, (i), (ii) and (iii) of Theorem 3.1 can be derived from Theorems 1.1, 1.3 and 1.6, 1.9 of [13], respectively.
(ii) In dimension N = 4 , the right hand side of the condition (1.6) in [14] is positive infinity. In dimension N\geqslant 5 , instead, the value is finit. For more detail, we refer to the formula (6.1) in [14].
Proof of Proposition 3.1. By Theorem 3.1, equation (1.7) has at least one positive, radially symmetric solution v when the assumptions in any case are satisfied. We define the function f_v(\beta) as
\begin{equation} f_v(\beta) = \beta-a-b \beta^{\frac{N-2}{2}} \|\nabla v\|_{L^2( \mathbb{R}^N)}^{2}. \end{equation} | (3.1) |
Case 1: N = 4 . In this case, by (1.6), (3.1), we have
f_v(\beta) = (1-d)\beta-a. |
It follows that if d\in(0, 1) , then
\tilde{\beta}_4 = \frac{a}{1-d}, |
defined as in (1.5), is a solution of f_v(\beta) = 0 , and (v, \tilde{\beta}_4) is a solution of \left(\mathcal{P}_{N}\right). If d\in[1, +\infty) , \left(\mathcal{P}_{N}\right) has no nontrivial solution. This implies that Proposition 3.1 (i) holds.
Case 2: N = 6 . In this case, by (1.6), (3.1), we have
f_v(\beta) = -d\beta^2+\beta-a. |
Then it is easy to know that \lim\limits_{\beta\to+\infty}f_v(\beta) = -\infty for any a, b > 0 and f_v(\beta) < 0 for any \beta\in(0, a] . If d\in(0, \frac{1}{4a}) , then there exist
\tilde{\beta}_{6,1} = \frac{1-\sqrt{1-4ad}}{2d} \in(a,\frac{1}{2d}), |
and
\tilde{\beta}_{6,2} = \frac{1+\sqrt{1-4ad}}{2d}\in(\frac{1}{2d},+\infty), |
such that (v, \tilde{\beta}_{6, 1}) , (v, \tilde{\beta}_{6, 2}) solve \left(\mathcal{P}_{N}\right). If d = \frac{1}{4a} , then (v, \frac{1}{2d}) solves \left(\mathcal{P}_{N}\right). If d\in(\frac{1}{4d}, +\infty) , then \left(\mathcal{P}_{N}\right) has no nontrivial solution. Thus, we can obtain that Proposition 3.1 (ii) holds.
Case 3: N = 8 . In this case, by (1.6), (3.1), we have
f_v(\beta) = -d\beta^3+\beta-a. |
Then it is easy to know that \lim\limits_{\beta\to+\infty}f_v(\beta) = -\infty for any a, b > 0 and f_v(\beta) < 0 for any \beta\in(0, a] . For each v\neq0 , we have
\frac{\mathrm{d}}{\mathrm{d} \beta} f_{v}(\beta) = 1-3d\beta^{2}. |
It is easy to see that f_{v}(\beta) has a unique maximum point \beta_{v} = \frac{1}{\sqrt{3d}} > 0, and
\max\limits_{\beta \in \mathbb{R}^{+}} f_{v}(\beta) = f_{v}\left(\beta_{v}\right) = \frac{2}{3\sqrt{3d}}-a. |
If d\in(0, \frac{4}{27a^2}) , it follows from Lemma 2.3 that the equation f_v(\beta) = 0 has three distinct real solutions
\tilde{\beta}_{8,j} = \frac{2}{\sqrt{3d}}\cos(\theta+\frac{2j\pi}{3})(j = 0,1,2), |
where
\begin{equation} \theta = \frac{1}{3}\arccos\left(-\frac{3a}{2}\sqrt{3d}\right). \end{equation} | (3.2) |
It is easy to see that one of which is negative and two of which are positive. Due to \beta > 0 , we next wish to exclude the negative solution. Note that (3.2) and d\in(0, \frac{4}{27a^2}), it follows that
-\frac{3a}{2}\sqrt{3d}\in(-1,0), |
and then
\theta\in\left(\frac{\pi}{6},\frac{\pi}{3}\right). |
This leads to \tilde{\beta}_{8, 0}, \tilde{\beta}_{8, 2} > 0 \text{ and } \tilde{\beta}_{8, 1} < 0. Moreover, it is easy to see that \tilde{\beta}_{8, 0} > \tilde{\beta}_{8, 2} . Therefore, we conclude that if d\in(0, \frac{4}{27a^2}) , then there exist
\tilde{\beta}_{8,0} = \frac{2}{\sqrt{3d}}\cos\theta\in(\frac{1}{\sqrt{3d}},\frac{1}{\sqrt{d}}), |
and
\tilde{\beta}_{8,2} = \frac{2}{\sqrt{3d}}\cos(\theta+\frac{4\pi}{3})\in(a,\frac{1}{\sqrt{3d}}), |
such that (v, \tilde{\beta}_{8, 0}) , (v, \tilde{\beta}_{8, 2}) solve \left(\mathcal{P}_{N}\right). If d = \frac{4}{27a^2} , the equation f_v(\beta) = 0 has a unique positive real solution
\tilde{\beta}_{8} = \frac{1}{\sqrt{3d}}, |
then (v, \tilde{\beta}_{8}) solves \left(\mathcal{P}_{N}\right). If d\in(\frac{4}{27a^2}, +\infty) , then \left(\mathcal{P}_{N}\right) has no nontrivial solution. Thus, we can obtain that Proposition 3.1 (iii) hold.
Proof of Proposition 3.2. Similar to the proof of Proposition 3.1, and applying Theorem 3.2, we can prove Proposition 3.2.
Proof of Theorem 1.1. Firstly, it is clear that u and v have the same sign and radial symmetry. We next prove the equivalence.
Necessary condition. Suppose the equation (\mathcal{K}) has a normalized solution u\in S(c) , we then have
\begin{aligned} &\left(a+b\|\nabla u\|_{L^2(\mathbb{R}^N)}^2\right)\int_{\mathbb{R}^{N}}\nabla u \nabla \varphi \mathrm{d}x +\lambda \int_{\mathbb{R}^{N}}u\varphi \mathrm{d}x\\ = &\mu \int_{\mathbb{R}^{N}}|u|^{q-2}u\varphi \mathrm{d}x+\int_{\mathbb{R}^{N}}|u|^{p-2}u\varphi \mathrm{d}x, \;\forall \varphi\in H^1( \mathbb{R}^N). \end{aligned} |
Let \beta = a+b\|\nabla u\|_{L^2(\mathbb{R}^N)}^2 and v(x) = u(\beta^{\frac{1}{2}}x) = u(y) , we conclude that
\begin{equation} \beta = a+b \beta^{\frac{N-2}{2}}\|\nabla v\|_{L^2(\mathbb{R}^N)}^2. \end{equation} | (4.1) |
Let \tilde{\beta}_N be defined by (1.5), i.e.,
\tilde{\beta}_N = \begin{cases} \frac{a}{1-d},& N = 4,\\ \frac{1\pm \sqrt{1-4ad}}{2d},& N = 6,\\ \frac{2}{\sqrt{3d}}\cos\left(\theta+\frac{2j\pi}{3}\right)(j = 0,2),&N = 8,\\ \end{cases} |
where d = b\|\nabla v\|_{L^2(\mathbb{R}^N)}^{2} , by a similar argument as used in the proof of Proposition 3.1, we obtain that \beta = \tilde{\beta}_N > 0 can solve (4.1) in dimensions N = 4, 6, 8 . Furthermore, let \phi(x) = \varphi(\beta^{\frac{1}{2}}x) = \varphi(y) , we have
\begin{aligned} &\int_{\mathbb{R}^{N}}\nabla v(x)\nabla\phi(x)\mathrm{d}x+\lambda \int_{\mathbb{R}^{N}}v(x)\phi(x)\mathrm{d}x\\ = &\tilde{\beta}_N^{-\frac{N}{2}} \left(\tilde{\beta}_N \int_{\mathbb{R}^{N}}\nabla u(y)\nabla\varphi(y)\mathrm{d}y+\lambda \int_{\mathbb{R}^{N}}u(y)\varphi(y)\mathrm{d}y\right)\\ = &\tilde{\beta}_N^{-\frac{N}{2}}\left(a+b \int_{\mathbb{R}^{N}}|\nabla u(y)|^2\mathrm{d}y\right) \int_{\mathbb{R}^{N}}\nabla u(y)\nabla\varphi(y)\mathrm{d}y+\lambda\tilde{\beta}_N^{-\frac{N}{2}} \int_{\mathbb{R}^{N}}u(y)\varphi(y)\mathrm{d}y\\ = &\tilde{\beta}_N^{-\frac{N}{2}}\mu \int_{\mathbb{R}^{N}}|u(y)|^{q-2}u(y)\varphi(y) \mathrm{d}y+\tilde{\beta}_N^{-\frac{N}{2}}\int_{\mathbb{R}^{N}}|u(y)|^{p-2}u(y)\varphi(y) \mathrm{d}y\\ = &\mu \int_{\mathbb{R}^{N}}|v(x)|^{q-2}v(x)\phi(x) \mathrm{d}x+\int_{\mathbb{R}^{N}}|v(x)|^{p-2}v(x)\phi(x) \mathrm{d}x. \end{aligned} |
Moreover, \int_{ \mathbb{R}^N}|v|^2\mathrm{d}x = \tilde{\beta}_N^{-\frac{N}{2}}\int_{ \mathbb{R}^N}|u(y)|^2\mathrm{d}y = \tilde{\beta}_N^{-\frac{N}{2}}c^2, which implies that (v, \tilde{\beta}_N)\in S(\tilde{\beta}_N^{-\frac{N}{4}}c)\times \mathbb{R}^+ is a solution of \left(\mathcal{P}_{N}\right).
Sufficient condition. Suppose the system \left(\mathcal{P}_{N}\right) has a solution (v, \tilde{\beta}_N)\in S(\tilde{\beta}_N ^{-\frac{N}{4}}c)\times \mathbb{R}^+ . At least in a weak sense, we then have, in dimensions N = 4, 6, 8,
\int_{\mathbb{R}^{N}}\nabla v\nabla\phi\mathrm{d}x+\lambda \int_{\mathbb{R}^{N}}v\phi\mathrm{d}x = \mu \int_{\mathbb{R}^{N}}|v|^{q-2}v\phi \mathrm{d}x+\int_{\mathbb{R}^{N}}|v|^{p-2}v\phi \mathrm{d}x, \;\forall \phi\in H^1( \mathbb{R}^N), |
and
\tilde{\beta}_N = a+b \tilde{\beta}_N^{\frac{N-2}{2}}\int_{\mathbb{R}^{N}}|\nabla v|^2\mathrm{d}x. |
Let u(x) = v(\tilde{\beta}_N^{-\frac{1}{2}}x) = v(y) and \varphi(x) = \phi(\tilde{\beta}_N^{-\frac{1}{2}}x) = \phi(y) we can obtain that
\begin{aligned} &\left(a+b \int_{\mathbb{R}^{N}}|\nabla u(x)|^2\mathrm{d}x\right) \int_{\mathbb{R}^{N}}\nabla u(x)\nabla\varphi(x)\mathrm{d}x+\lambda\int_{\mathbb{R}^{N}}u(x)\varphi(x)\mathrm{d}x\\ = &\tilde{\beta}_N^{\frac{N}{2}} \left(\tilde{\beta}_N^{-1}\left(a+b\tilde{\beta}_N^{\frac{N-2}{2}} \int_{\mathbb{R}^{N}}|\nabla v(y)|^2\mathrm{d}y\right) \int_{\mathbb{R}^{N}}\nabla v(y)\nabla\phi(y)\mathrm{d}y+\lambda\int_{\mathbb{R}^{N}}v(y)\phi(y)\mathrm{d}y\right)\\ = &\tilde{\beta}_N^{\frac{N}{2}} \left( \int_{\mathbb{R}^{N}}\nabla v(y)\nabla\phi(y)\mathrm{d}y+\lambda\int_{\mathbb{R}^{N}}v(y)\phi(y)\mathrm{d}y\right)\\ = &\tilde{\beta}_N^{\frac{N}{2}} \left( \mu \int_{\mathbb{R}^{N}}|v(y)|^{q-2}v(y)\phi(y) \mathrm{d}y+\int_{\mathbb{R}^{N}}|v(y)|^{p-2}v(y)\phi(y) \mathrm{d}y\right)\\ = &\mu \int_{\mathbb{R}^{N}}|u(x)|^{q-2}u(x)\varphi(x) \mathrm{d}x+\int_{\mathbb{R}^{N}}|u(x)|^{p-2}u(x)\varphi(x) \mathrm{d}x, \end{aligned} |
and
\int_{ \mathbb{R}^N}|u|^2\mathrm{d}x = \tilde{\beta}_N^{\frac{N}{2}}\int_{ \mathbb{R}^N}|v(y)|^2\mathrm{d}y = c^2, |
which implies that u\in S(c) is a solution of (\mathcal{K}).
Proof of Theorem 1.2. The conclusions follow directly from Theorem 1.1 and Proposition 3.1.
Proof of Theorem 1.3. The conclusions follow directly from Theorem 1.1 and Proposition 3.2.
Proof of Theorem 1.4. Suppose that u\in S(c) is a positive solution of (\mathcal{K}) with p = q = 2^* , then
a\|\nabla u\|_{L^2(\mathbb{R}^N)}^2+b\|\nabla u\|_{L^2(\mathbb{R}^N)}^4+\lambda \|u\|_{L^2(\mathbb{R}^N)}^2 = \mu \|u\|^q_{L^q(\mathbb{R}^N)}+\|u\|^p_{L^p(\mathbb{R}^N)}. |
It follows from Lemma 2.2 and u\in S(c) that
\lambda c^2 = \mu(1-\gamma_q)\|u\|^q_{L^q(\mathbb{R}^N)}+(1-\gamma_p)\|u\|^p_{L^p(\mathbb{R}^N)}. |
Since p = q = 2^* , we have \gamma_p = \gamma_q = 1 , it follows that \lambda c^2 = 0 . However, by Lemma 2.1, we have \lambda > 0 , which is a contradiction. Hence, equation (\mathcal{K}) has no positive solutions u\in S(c) for any \mu > 0 .
This work is supported by National Natural Science Foundation of China (No. 12071486).
The authors declare there is no conflicts of interest.
[1] |
A. M. Arvin, Varicella-zoster virus, Clin. Microbiol. Rev., 9 (1996), 361–381. https://doi.org/10.1128/CMR.9.3.361 doi: 10.1128/CMR.9.3.361
![]() |
[2] |
A. A. Gershon, J. Breuer, J. I. Cohen, R. J. Cohrs, M. D. Gershon, D. Gilden, et al., Varicella zoster virus infection, Nat. Rev. Dis. Primers, 1 (2015), 15016. https://doi.org/10.1038/nrdp.2015.16 doi: 10.1038/nrdp.2015.16
![]() |
[3] |
M. Marin, M. Marti, A. Kambhampati, S. M. Jeram, J. F. Seward, Global varicella vaccine effectiveness: A meta-analysis, Pediatrics, 137 (2016), e20153741. https://doi.org/10.1542/peds.2015-3741 doi: 10.1542/peds.2015-3741
![]() |
[4] |
A. Suzuki, H. Nishiura, Reconstructing the transmission dynamics of varicella in Japan: an elevation of age at infection, PeerJ, 10 (2022), e12767. https://doi.org/10.7717/peerj.12767 doi: 10.7717/peerj.12767
![]() |
[5] |
A. Suzuki, H. Nishiura, Transmission dynamics of varicella before, during and after the COVID-19 pandemic in Japan: A modelling study, Math. Biosci. Eng., 19 (2022), 5998–6012. https://doi.org/10.3934/mbe.2022280 doi: 10.3934/mbe.2022280
![]() |
[6] |
H. F. Gidding, M. Brisson, C. R. Macintyre, M. A. Burgess, Modelling the impact of vaccination on the epidemiology of varicella zoster virus in Australia, Aust. N. Z. J. Public Health, 29 (2005), 544–551. https://doi.org/10.1111/j.1467-842x.2005.tb00248.x doi: 10.1111/j.1467-842x.2005.tb00248.x
![]() |
[7] |
M. Karhunen, T. Leino, H. Salo, I. Davidkin, T. Kilpi, K. Auranen, Modelling the impact of varicella vaccination on varicella and zoster, Epidemiol. Infect., 138 (2010), 469–481. https://doi.org/10.1017/S0950268809990768 doi: 10.1017/S0950268809990768
![]() |
[8] |
F. Lienert, O. Weiss, K. Schmitt, U. Heininger, P. Guggisberg, Acceptance of universal varicella vaccination among Swiss pediatricians and general practitioners who treat pediatric patients, BMC Infect. Dis., 21 (2021), 12. https://doi.org/10.1186/s12879-020-05586-3 doi: 10.1186/s12879-020-05586-3
![]() |
[9] |
W. P. London, J. A. Yorke, Recurrent outbreaks of measles, chickenpox and mumps. I. Seasonal variation in contact rates, Am. J. Epidemiol., 98 (1973), 453–468. https://doi.org/10.1093/oxfordjournals.aje.a121575 doi: 10.1093/oxfordjournals.aje.a121575
![]() |
[10] |
N. C. Grassly, C. Fraser, Seasonal infectious disease epidemiology, Proc. Biol. Sci., 273 (2006), 2541–2550. https://doi.org/10.1098/rspb.2006.3604 doi: 10.1098/rspb.2006.3604
![]() |
[11] |
P. E. Fine, J. A. Clarkson, Measles in England and Wales--Ⅰ: An analysis of factors underlying seasonal patterns, Int. J. Epidemiol., 11 (1982), 5–14. https://doi.org/10.1093/ije/11.1.5 doi: 10.1093/ije/11.1.5
![]() |
[12] |
B. Finkenstädt, B. Grenfell, 2000 Time series modelling of childhood diseases: A dynamical systems approach, Appl. Statist., 49 (2000), 187–205. https://doi.org/10.1111/1467-9876.00187 doi: 10.1111/1467-9876.00187
![]() |
[13] |
D. He, D. J. Earn, The cohort effect in childhood disease dynamics, J. R. Soc. Interface, 13 (2016), 20160156. https://doi.org/10.1098/rsif.2016.0156 doi: 10.1098/rsif.2016.0156
![]() |
[14] |
C. Jackson, P. Mangtani, P. Fine, E. Vynnycky, The effects of school holidays on transmission of varicella zoster virus, England and Wales, 1967-2008, PLoS One, 9 (2014), e99762. https://doi.org/10.1371/journal.pone.0099762 doi: 10.1371/journal.pone.0099762
![]() |
[15] | D. L. Heymann, American Public Health Association, in Control of Communicable Diseases Manual (20th edition), American Public Health Association Publications, Washington, 2015. |
[16] |
R. E. Baker, A. S. Mahmud, C. J. E. Metcalf, Dynamic response of airborne infections to climate change: predictions for varicella, Clim. Change, 148 (2018), 547–560. https://doi.org/10.1007/s10584-018-2204-4 doi: 10.1007/s10584-018-2204-4
![]() |
[17] | Statistical Handbook of Japan, Chapter 1 Land and Climate, 2021. Available from: https://www.stat.go.jp/english/data/handbook/c0117.html. |
[18] | Japan Meteorological Agency, General Information on Climate of Japan. Available from: https://www.data.jma.go.jp/gmd/cpd/longfcst/en/tourist.html. |
[19] | Ministry of Health, Labour and Welfare, Available from: https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/kenkou_iryou/kenkou/kekkaku-kansenshou/kekkaku-kansenshou11/01.html#list05 (in Japanese). |
[20] | Japan Meteorological Agency. Available from: https://www.data.jma.go.jp/obd/stats/etrn/ (in Japanese). |
[21] |
K. Harigane, A. Sumi, K. Mise, N. Kobayashi, The role of temperature in reported chickenpox cases from 2000 to 2011 in Japan, Epidemiol. Infect., 143 (2015), 2666–2678. https://doi.org/10.1017/S095026881400363X doi: 10.1017/S095026881400363X
![]() |
[22] | National Institute of Population and Social Security Research, Population Projections for Japan. Available from: http://www.ipss.go.jp/pp-zenkoku/e/zenkoku_e2017/pp29_summary.pdf. |
[23] |
T. Ozaki, Long-term clinical studies of varicella vaccine at a regional hospital in Japan and proposal for a varicella vaccination program, Vaccine, 31 (2013), 6155–6160. https://doi.org/10.1016/j.vaccine.2013.10.060 doi: 10.1016/j.vaccine.2013.10.060
![]() |
[24] |
O. N. Bjørnstad, B. F. Finkenstädt, B. T. Grenfell, Dynamics of measles epidemics: estimating scaling of transmission rates using a time series SIR model, Ecol. Monogr., 72 (2002), 169–184. https://doi.org/10.2307/3100023 doi: 10.2307/3100023
![]() |
[25] |
J. Shaman, V. E. Pitzer, C. Viboud, B. T. Grenfell, M. Lipsitch, Absolute humidity and the seasonal onset of influenza in the continental United States, PLoS Biol., 8 (2010), e1000316. https://doi.org/10.1371/journal.pbio.1000316 doi: 10.1371/journal.pbio.1000316
![]() |
[26] |
J. D. Tamerius, J. Shaman, W. J. Alonso, K. Bloom-Feshbach, C. K. Uejio, A. Comrie, et al., Environmental predictors of seasonal influenza epidemics across temperate and tropical climates, PLoS Pathog., 9 (2013), e1003194. https://doi.org/10.1371/journal.ppat.1003194 doi: 10.1371/journal.ppat.1003194
![]() |
[27] |
J. B. Axelsen, R. Yaari, B. T. Grenfell, L. Stone, Multiannual forecasting of seasonal influenza dynamics reveals climatic and evolutionary drivers, Proc. Natl. Acad. Sci.U. S. A., 111 (2014), 9538–9542. https://doi.org/10.1073/pnas.1321656111 doi: 10.1073/pnas.1321656111
![]() |
[28] |
A. C. Lowen, J. Steel, Roles of humidity and temperature in shaping influenza seasonality, J. Virol., 88 (2014), 7692–7695. https://doi.org/10.1128/JVI.03544-13 doi: 10.1128/JVI.03544-13
![]() |
[29] |
D. E. te Beest, M. van Boven, M. Hooiveld, C. van den Dool, J. Wallinga, Driving factors of influenza transmission in the Netherlands, Am. J. Epidemiol., 178 (2013), 1469–1477. https://doi.org/10.1093/aje/kwt132 doi: 10.1093/aje/kwt132
![]() |
[30] |
H. Yuan, S. C. Kramer, E. H. Y. Lau, B. J. Cowling, W. Yang, Modeling influenza seasonality in the tropics and subtropics, PLoS Comput. Biol., 17 (2021), e1009050. https://doi.org/10.1371/journal.pcbi.1009050 doi: 10.1371/journal.pcbi.1009050
![]() |
[31] |
S. T. Ali, B. J. Cowling, J. Y. Wong, D. Chen, S. Shan, E. H. Y. Lau, et al., Influenza seasonality and its environmental driving factors in mainland China and Hong Kong, Sci. Total Environ., 818 (2022), 151724. https://doi.org/10.1016/j.scitotenv.2021.151724 doi: 10.1016/j.scitotenv.2021.151724
![]() |
[32] |
N. Toyama, K. Shiraki, Epidemiology of herpes zoster and its relationship to varicella in Japan: A 10-year survey of 48,388 herpes zoster cases in Miyazaki prefecture, J. Med. Virol., 12 (2009), 2053–2058. https://doi.org/10.1002/jmv.21599 doi: 10.1002/jmv.21599
![]() |
1. | Amjad Ali, Muhammad Arshad, Eskandar Ameer, Asim Asiri, Certain new iteration of hybrid operators with contractive M -dynamic relations, 2023, 8, 2473-6988, 20576, 10.3934/math.20231049 | |
2. | Muhammad Tariq, Muhammad Arshad, Mujahid Abbas, Eskandar Ameer, Saber Mansour, Hassen Aydi, A relation theoretic m-metric fixed point algorithm and related applications, 2023, 8, 2473-6988, 19504, 10.3934/math.2023995 | |
3. | Imo Kalu Agwu, Naeem Saleem, Umar Isthiaq, A new modified mixed-type Ishikawa iteration scheme with error for common fixed points of enriched strictly pseudocontractive self mappings and ΦΓ-enriched Lipschitzian self mappings in uniformly convex Banach spaces, 2025, 26, 1989-4147, 1, 10.4995/agt.2025.17595 | |
4. | Muhammad Tariq, Sabeur Mansour, Mujahid Abbas, Abdullah Assiry, A Solution to the Non-Cooperative Equilibrium Problem for Two and Three Players Using the Fixed-Point Technique, 2025, 17, 2073-8994, 544, 10.3390/sym17040544 |