1 Introduction

In the recent years, the following Cauchy problem for the wave equation with scale invariant damping spreads a new line of research on variable coefficient-type equations. More precisely, we are dealing with

$$\begin{aligned} {\left\{ \begin{array}{ll} v_{tt}(t,x)-\varDelta v(t,x)+\frac{\mu }{1+t}v_t(t,x)+\frac{\nu }{(1+t)^2}v(t,x)={\left| v(t,x)\right| }^p, \quad &{}t\ge 0,\quad x\in {\mathbb {R}}^n,\\ v(0,x)=0, \\ v_t(0,x)=\varepsilon g(x), \end{array}\right. } \end{aligned}$$
(1)

with \(n\ge 2\), \(\mu ,\,\nu \in {\mathbb {R}}\), \(p>1\) and g a radial smooth function. In [2, 3, 10, 11], some results on the global existence of a solution for (1) with non-compactly initial data appeared assuming a suitable decay behavior for g. Many other results concern blowup and global existence for this equation, see [12] for a summary of this problem. The main point is to find a critical exponent, fixed a suitable space of data. More precisely, a level \({\bar{p}}\) is critical if for \(p >\bar{p}\), one can prove that for \(\epsilon >0\) sufficiently small and for any g chosen in the fixed space, there exists a unique global (in time) solution of the problem, and conversely, if \(p\in (1, {\bar{p}})\), for any \(\epsilon >0\) there exist some g in this space such that the local solution cannot be prolonged over a finite time. Coming back to (1), in dependence on \(\mu ,\,\nu \) and n,  a competition between two critical exponents appeared. In some cases, the Strauss exponent is dominant; it is given by the wave equation theory, it will be denoted by \(p_S(d)\), and it is the positive root of the quadratic equation

$$\begin{aligned} (d-1)p^2-(d+1)p-2=0\,. \end{aligned}$$

For other assumptions, the equation goes to an heat equation and a Fujita-type exponent \(p_F(h):=1+\frac{2}{h}\) appears. In all known results, the quantities \(d>1, h> 0\) depend on \(\nu ,\,\mu \) and n. Changing the space of data, a change of critical exponent may appear. The novelty of our result consists in showing that if one takes into account the decay rate of the initial data, then the Fujita-type exponent depends also on such decay rate. In addition, we give an upper bound estimate for the lifespan of the solution, in terms of the power of the nonlinear term, the size and the growth of the initial data. Let us recall that the lifespan of the solution is a function of \(\varepsilon \) which gives the maximal existence time:

$$\begin{aligned} T(\varepsilon ):=\sup \,\{ T>0 \text { such that the local solution }\, u \,\text { to } (1) \text { is defined on } [0,T)\times {\mathbb {R}}^{n}\,\}. \end{aligned}$$

Finally, we will prove the following.

Theorem 1

Let \(n\ge 2\). Let \(\epsilon >0\) and g be a radial smooth function satisfying

$$\begin{aligned} g(|x|)\ge \frac{M}{(1+|x|)^{{\bar{k}}+1},} \text { with } \bar{k}>-1, \end{aligned}$$
(2)

for some \(M>0\) and for any \(x \in {\mathbb {R}}^n\). Assuming in addition that

$$\begin{aligned} \bar{k}+\frac{\mu }{2}>0,\quad \frac{\mu }{2}\left( \frac{\mu }{2}-1\right) \ge \nu \end{aligned}$$

and

$$\begin{aligned} 1<p<p_F\left( {\bar{k}}+\frac{\mu }{2}\right) , \end{aligned}$$

then the classical solution of (1) blows up. More precisely, the lifespan of the solution \(T(\varepsilon )>0\) is finite and satisfies

$$\begin{aligned} T(\varepsilon )\le C \epsilon ^{-\frac{2(p-1)}{4-(\mu +2\bar{k})(p-1)}}, \end{aligned}$$
(3)

with \(C>0\), independent of \(\varepsilon \).

Remark 1

Recently, Ikeda, Tanaka, Wakasa in [6] consider a similar question for cubic convolution nonlinearity and a critical decay appears.

Remark 2

In [4], we will also consider a variant of problem (1), in which the nonlinearity depends on \(v,t,v_t\) combined in a suitable way.

Remark 3

The lifespan estimate for the same equation with compactly supported data and \(\nu \ne \frac{\mu }{2}\left( \frac{\mu }{2}-1\right) \) has been considered in [13]. If \(\nu \le \frac{\mu }{2}\left( \frac{\mu }{2}-1\right) \), the lifespan estimate is different from (3) due to the compactness of the support of the initial data.

The paper is organized as follows: In Sect. 2, we give an overview of the known results and we state an auxiliary theorem; in Sect. 3, we prove the main results.

2 Motivations

2.1 The case \(\mu =2\), \(\nu =0\)

Let us start with a quite simple case

$$\begin{aligned} {\left\{ \begin{array}{ll} v_{tt}-\varDelta v+\frac{2}{1+t}v_t={\left| v\right| }^p, \quad &{}t\ge 0,\quad x\in {\mathbb {R}}^n,\\ v(0,x)=0, \\ v_t(0,x)=\varepsilon g(x). \end{array}\right. } \end{aligned}$$
(4)

The global existence of small data solutions for this problem was first solved in [1] for a suitable range of n and p. Some nonexistence results were also established for \(p<p_F(n):=1+\frac{2}{n}\). Except for the one-dimensional case, a gap between this value and the admissible exponents in [1] appeared. In [3] for dimension \(n=2,3,\) this gap was covered with an unexpected result. Indeed, in that paper the Strauss exponent came into play. Afterward, the global existence of small data solutions to (4) has been proved for any \(p>p_S(n+2)\) also in odd dimension \(n\ge 5\) in [2] and in even dimension \(n\ge 4\) in [10].

Moreover, we know that the exponent \(p_2(n):=\max \{p_S(n+2),p_F(n)\}\) is optimal; in fact, in [3], the authors prove the blowup of solutions of (4) for each \(1<p\le p_2(n)\) in each dimension \(n\in {\mathbb {N}}\). In [2, 3, 10], the authors prove a global existence result not necessarily when the initial datum \(g=g(x)\) has compact support. More precisely, let \(n\ge 3\), given a radial initial datum \(g(x)=g(|x|)\) with \(g\in C^1({\mathbb {R}})\), for any \(p>p_S(n+2)\) it is possible to choose \({\bar{k}}>0\) and \(\epsilon _0>0\) such that (4) admits a radial global solution \(u\in C([0,\infty )\times {\mathbb {R}}^n)\cap C^2([0,\infty )\times {\mathbb {R}}^n{\setminus } \{0\})\) provided

$$\begin{aligned} \varepsilon \le \varepsilon \,\, 0,\quad \qquad |g^{(h)}(r)|\le \epsilon \langle r \rangle ^{-(\bar{k}+1+h)} \text { for } h=0,\,1. \end{aligned}$$
(5)

In the present paper, we discuss the dependence of \({\bar{k}}\) from n and p. In (5), the exponent \(\bar{k}\) has to belong to a suitable interval \( [k_1(n,p),k_2(n,p)]\). It is interesting to investigate the case of \({\bar{k}}\not \in [k_1(n,p), k_2(n,p)]\). In the sequel, we will see that the bound \(k_2(n,p)\) can be easily improved (see Remark 4). On the contrary, if \(k<k_1(n,p)\), then a new result appears. The known situation is the following:

  • \(k_1(3,p)=\max \big \{\frac{3-p}{p-1}, \frac{1}{p-1}\big \}\) and \(k_2(3,p)=2(p-1)\), see [3].

  • \(k_1(n,p)= \max \big \{\frac{3-p}{p-1}, \frac{n-1}{2}\big \}\) and \(k_2(n,p)= \min \big \{\frac{(n+1)p}{2}-2,\frac{n^2-2n+13}{2(n-3)}\big \}\) if \(n\ge 5\) odd, see [2].

  • \(k_1(n,p)= \max \big \{\frac{3-p}{p-1},\frac{n-1}{2}\big \}\) and \(k_2(n,p)= \min \left\{ \frac{(n+1)p}{2}-2, n-1\right\} \) if \(n\ge 4\), see [10].

We can write in a different way the previous conditions. Firstly, we concentrate on the case \(n=3\). For \(p\in (1,2)\), we have \({\bar{k}}\ge \frac{3-p}{p-1}\) that is equivalent to

$$\begin{aligned} p\ge 1+\frac{2}{{\bar{k}}+1}=p_F({\bar{k}}+1). \end{aligned}$$

From above, we have \({\bar{k}}\le 2(p-1)\) that is

$$\begin{aligned} p\ge \frac{{\bar{k}}}{2}+1. \end{aligned}$$

The intersection of \(p=p_F({\bar{k}}+1)\) and \(p=1+{\bar{k}}/2\) is exactly in \({\bar{k}}= \frac{-1+\sqrt{17}}{2}\) and \(p=p_S(5)\). We summarize the situation in Fig. 1. In the following graphs, we denote in blue the zone of the known global existence results, in red the zone of the known blow-up results. In this paper, we want to cover the white zones.

Fig. 1
figure 1

\(n=3,\; \mu =2,\; \nu =0\)

Reading [2], we see that the same situation appears for any odd \(n\ge 5\). The critical curve

$$\begin{aligned} p=p_F({\bar{k}}+1) \end{aligned}$$

intersects the line

$$\begin{aligned} p=\frac{2({\bar{k}}+2)}{n+1} \end{aligned}$$

in the Strauss couple

$$\begin{aligned} \left( {\bar{k}}_0,\frac{2({\bar{k}}_0+2)}{n+1} \right) = \left( \frac{n-5+\sqrt{n^2+14n+17}}{4}, p_S(n+2) \right) \,. \end{aligned}$$

The only difference with the case \(n=3\) is that, in the global existence zone, a bound from above appears for p and this has some influence on \(k_2(n,p)\). More precisely, one can take

$$\begin{aligned} p\le \frac{n+1}{n-3},\quad {\bar{k}}\le \frac{n^2-2n+13}{2(n-3)}\,\quad \text { if } n\ge 7, \end{aligned}$$

and \(p\le 2\), \({\bar{k}}\le 3\) if \(n=5\). Hence, the result of such paper can be represented as in Fig. 2. Our aim is to prove blowup in the white zone below the Fujita curve.

Fig. 2
figure 2

\(n\ge 5 \text { odd, }\; \mu =2,\; \nu =0\)

Even dimension is more delicate. In [10], the global existence result is established in the blue zone below the line \(p=\frac{n+5}{n+1}\) except on the curve \(p=p_F({\bar{k}}+1)\). For convenience of the reader, we precise that in the notation of [10], the role of \({\bar{k}}\) is taken by the quantity \(k+\frac{n+1}{2}\).

2.2 The case \(\mu >2\) and \(\nu =\frac{\mu }{2}(\frac{\mu }{2}-1)\)

In [10, 11], the author considers the Cauchy problem (1) for the semilinear wave equation with scale invariant damping and mass terms, that is \(\nu =\frac{\mu }{2}(\frac{\mu }{2}-1)\ge 0\). We see that for \(\mu =2,\) it reduces to (4). Global existence of solutions to (1) holds under the conditions

$$\begin{aligned} \mu \in [2,M(n)], \qquad M(n)=\frac{n-1}{2}\left( 1+\sqrt{\frac{n+7}{n-1}}\right) \,. \end{aligned}$$

In the even case [10], the initial data satisfy (5) for \({\bar{k}}\in (k_1(n,p,\mu ),k_2(n,p,\mu )]\) such that

$$\begin{aligned}&k_1(n,p,\mu ) = \max \Big \{\frac{n-1}{2}, \frac{2}{p-1}-\frac{\mu }{2}\Big \}\,; \end{aligned}$$
(6)
$$\begin{aligned}&k_2(n,p,\mu ) = \min \Big \{n-1, \frac{n+\mu -1}{2}p-\frac{\mu +2}{2}\Big \}. \end{aligned}$$
(7)

Rewriting these conditions in terms of p, we find that

$$\begin{aligned} p>p_F\left( {\bar{k}}+\frac{\mu }{2}\right) , \qquad p\ge \frac{2{\bar{k}}+\mu +2}{n+\mu -1}. \end{aligned}$$
Fig. 3
figure 3

\(n\ge 4 \text { even, }\; \mu =2,\; \nu =0\)

The intersection of the curves those define the global existence zone gives \(p=p_S(n+\mu )\). Hence, the condition \(p>p_S(n+\mu )\) appears. Moreover, another bound from above appears:

$$\begin{aligned} p<{\bar{p}}:=\min \left\{ p_F(\mu ), p_F\left( \frac{n+\mu -1}{2}\right) \right\} . \end{aligned}$$

This means that different results for large \(\mu \) and small \(\mu \) hold. This influences the positions of \(k_1\) and \(k_2\). For our purpose, it is sufficient to say that for \(\mu \not =2\) and even n, the situation is similar to Fig. 3. More precisely, in Fig. 4\(p_S(n+\mu )\) appears. The blow-up result is indeed given in [9]. The zone between \(p=p_S(n+\mu )\) and \(p=p_F\left( {\bar{k}} +\frac{\mu }{2}\right) \) is not covered by any known result.

Fig. 4
figure 4

\(n\ge 4 \text { even, }\;\nu =\frac{\mu }{2}(\frac{\mu }{2}-1)\ge 0\)

The corresponding global existence result for the Cauchy problem (1) in odd space dimension \(n\ge 1\) is studied in [11] for radial and small data, assuming condition (5) with \({\bar{k}} \in [k_1(n,p,\mu ),k_2(n,p,\mu )]\) where \(k_2\) satisfies (7) and it holds:

$$\begin{aligned}&k_1(3,p,\mu )=\max \Big \{1, \frac{2}{p-1}-\frac{\mu }{2}, \frac{1}{p-1}\Big \}\,;\\&k_1(n,p,\mu )=\max \Big \{\frac{n-1}{2}, \frac{2}{p-1}-\frac{\mu }{2}\Big \}\,, \quad n\ge 5\, \quad \mu \in [2,n-1];\\&k_1(n,p,\mu )=\max \Big \{\frac{n-1}{2}, \frac{2}{p-1}-\frac{\mu }{2}, \frac{1}{p-1}, \Big \}\,,\, n\ge 5\quad \mu \in (n-1,M(n)]. \end{aligned}$$

In any case, the condition \(p>p_F({\bar{k}}+\frac{\mu }{2})\) appears. Hence, in odd space dimension \(n\ge 5\) the situation is not different from Fig. 4.

Reading Theorem 1 in the case \(\nu =\frac{\mu }{2}(\frac{\mu }{2}-1)\), it is clear that the aim of this paper is to find blowing-up solutions to (1) even for \(p>p_S(n+\mu )\) by considering initial data with slow decay. More precisely, let us consider

$$\begin{aligned} g(x)\simeq \frac{M}{(1+{\left| x\right| })^{ \bar{k}+1}}\,,\quad \text { for } \quad \frac{n-1}{2}<{\bar{k}}<{\bar{k}}_0\,, \end{aligned}$$
(8)

where \({\bar{k}}_0\) is such that

$$\begin{aligned} p_F\left( {\bar{k}}_0+\frac{\mu }{2}\right) =p_S(n+\mu ). \end{aligned}$$

We will prove the blow-up result in the left white side zones in Figs. 1, 2, 3 and 4 where \({\bar{k}}< {\bar{k}}_0\), \(p>p_S(n+\mu )\) and \(p<p_F({\bar{k}}_0+\frac{\mu }{2})\). Under the same assumption on g, the quoted results assure that for \(p\ge p_F({\bar{k}}+\frac{\mu }{2})\) and \(p>p_S(n+\mu )\), there is global existence. Hence, \(p=p_F({\bar{k}}+\frac{\mu }{2})\) is a critical curve for the Cauchy problem (1), provided \(\nu =\frac{\mu }{2}(\frac{\mu }{2}-1)\ge 0\).

Remark 4

Still fixing \(\nu =\frac{\nu }{2}(\frac{\nu }{2}-1)\ge 0\), let us consider \(\bar{k}>\bar{k}_0\) and \(p>p_S(n+\mu )\). As discussed, the global existence results in the previous literature require p above a line which depends on \({\bar{k}}\), because of a restriction of type \(\bar{k}\le k_2(n,p,\mu )\) which every time appears. Actually, this restriction can be avoided; indeed, if the initial datum satisfies (8) with \({\bar{k}}> k_2(n,p,\mu )\), then we can say that the initial datum also satisfies (5) with \({\bar{k}}= k_2(n,p,\mu )\). Hence, the global existence of a solution to (1) follows from the known results.

Remark 5

For \(\nu =\frac{\mu }{2}(\frac{\mu }{2}-1)\ge 0,\) Theorem 1 provides some new information about the solution of (1) also when p belongs to the red zone of Figs. 1, 2, 3, 4 and 5. In fact, for

$$\begin{aligned} p<\min \left\{ p_S(n+\mu ), p_F\left( {\bar{k}}+\frac{\mu }{2}\right) \right\} \end{aligned}$$

by the previous literature we know that the solution blows up in finite time, whereas Theorem 1 gives a lifespan estimate in the case of radial initial data with non-compact support, relating this estimate with the decay rate of the data.

2.3 The case \(\mu =0\) and \(\nu =0\)

In Fig. 5, we summarize the wave equation case \(\mu =\nu =0\). The red blow-up zone was covered by many authors, see [14] and the reference therein for the whole list of blow-up results. For \(\mu =\nu =0\), the global existence result has been completely solved in [5], where the interested reader can find a long bibliography of previous contributes. In particular, the blue zone for radial solution without compact support assumption for the initial data has been exploited by Kubo, see, for example, [7, 8]. Before these papers, Takamura obtained a blow-up result in the green zone. In [15], the point is to find a critical decay level \(k_0=\frac{2}{p-1},\) equivalently \(p\le 1+\frac{2}{k_0}.\) We underline that this is a Fujita-type exponent.

Fig. 5
figure 5

\(\mu =\nu =0\)

In Theorem 1, we generalize Takamura’s result when \(\mu \ne 0\) and \(\nu \le \frac{\mu }{2}(\frac{\mu }{2}-1)\). To this aim, it is sufficient to consider a peculiar wave equation with nonlinear term having a decaying time-dependent variable coefficient. This means that we will deduce Theorem 1 from the following result.

Theorem 2

Let \(n\ge 2\). Given a smooth function \(g=g(|x|)\) with \(x\in {\mathbb {R}}^n\), we set \(r=|x|\) and we consider \(g=g(r)\) satisfying

$$\begin{aligned} g(r)\ge \frac{M}{(1+r)^{{\bar{k}}+1}}, \text { with } {\bar{k}}>-1, \end{aligned}$$
(9)

for some \(M>0\). Let \(u=u(t,r)\) be the radial local solution to

$$\begin{aligned} {\left\{ \begin{array}{ll} u_{tt}-u_{rr}-\frac{n-1}{r}u_r=(1+t)^{-\frac{\mu }{2}(p-1)}{\left| u\right| }^p, \quad &{}r>0,\\ u(0,r)=0, \\ u_t(0,r)=\varepsilon g(r). \end{array}\right. } \end{aligned}$$
(10)

with \(p>1\) and \(p<p_{F}(\frac{\mu }{2}-1)\) if \(\mu >2\). Assume in addition that

$$\begin{aligned} -1<{\bar{k}} < \frac{2}{p-1}-\frac{\mu }{2}. \end{aligned}$$
(11)

Then, given \(\epsilon >0\), the lifespan \(T(\varepsilon )>0\) of classical solutions to (10) satisfies

$$\begin{aligned} T(\varepsilon )\le C \epsilon ^{-\frac{2(p-1)}{4-(\mu +2\bar{k})(p-1)}}, \end{aligned}$$
(12)

with \(C>0\), independent of \(\varepsilon \).

Remark 6

The assumption \(p<p_{F}(\frac{\mu }{2}-1)\) if \(\mu >2\) guarantees that the range of admissible \({\bar{k}}\) in (11) is not empty.

In the case \(\mu =0\), Theorem 2 coincides with Takamura’s result in [15]. In the proof of Theorem 2, we will follow the same approach of that paper.

3 Proof of the main results

3.1 Proof of Theorem 2

We recall the crucial lemma of [15].

Lemma 1

Let \(n\ge 2\) and \(m=[n/2]\). Given a smooth function \(g=g(|x|)\) with \(x\in {\mathbb {R}}^n\), we set \(r=|x|\) and we consider \(g=g(r)\). Let us denote by \(u^{0}(t,r)\) the solution of the free wave problem

$$\begin{aligned} {\left\{ \begin{array}{ll} \square u^0=0 \quad &{} (t,r)\in [0,\infty )\times [0,\infty )\\ u^0(0,r)=0,~u^0_t(0,r)=g(r)\,. \end{array}\right. } \end{aligned}$$

Let \(u=u(t,r)\) be a solution to

$$\begin{aligned} u_{tt}-u_{rr}-\frac{n-1}{r}u_r=F(t,u) \end{aligned}$$
(13)

with the initial condition

$$\begin{aligned} u(0,r)=0, ~ u_t(0,r)=\varepsilon g(r), \quad r\in [0,\infty ). \end{aligned}$$
(14)

If F is nonnegative, there exists a constant \(\delta _m> 0\) such that

$$\begin{aligned} u(t,r)&\ge \varepsilon u^0(t,r)+\frac{1}{8r^m}\int _0^t d\tau \int _{r-t+\tau }^{r+t+\tau }\lambda ^m F(t,u(t,\lambda ))d \lambda , \end{aligned}$$
(15)
$$\begin{aligned} u^0(t,r)&\ge \frac{1}{8r^m} \int _{r-t}^{r+t} \lambda ^m g(\lambda ) d\lambda , \end{aligned}$$
(16)

provided

$$\begin{aligned} r-t\ge \frac{2}{\delta _m}t>0. \end{aligned}$$

The constant \(\delta _m\) in the previous lemma is described in [15, Lemma 2.5]; it depends on the space dimension, and in particular, it changes accordingly with the different representations of the free wave solution in odd and even dimension.

We are ready to prove that if (9) holds, then the solution of (10) blows up in finite time even for small \(\varepsilon .\)

Let us fix \(\delta >0\); we define a blow-up set,

$$\begin{aligned} \Sigma _\delta =\left\{ (t,r)\in (0,\infty )^2: r-t\ge \max \left\{ \frac{2}{\delta _m}t,\delta \right\} \right\} , \end{aligned}$$
(17)

where \(\delta _m>0\) is the constant given in Lemma 1. Combining the assumption (10) with the formulas (15) and (16), for any \((t,r)\in \Sigma _\delta ,\) it holds

$$\begin{aligned} u(t,r)\ge \varepsilon u^0(t,r)\ge \frac{\varepsilon }{8r^m}\int _{r-t}^{r+t}\lambda ^m g(\lambda )d\lambda \ge \frac{M\varepsilon }{8r^m}\int _{r-t}^{r+t}\lambda ^m (1+\lambda )^{-({\bar{k}} +1)}d\lambda \,. \end{aligned}$$

Then, (17) implies that

$$\begin{aligned} u(t,r)&\ge \frac{M\varepsilon }{8r^m}\left( \frac{1+\delta }{\delta }\right) ^{-({\bar{k}}+1)}\int _{r-t}^{r+t}\lambda ^{m-({\bar{k}}+1)} d\lambda \\&\ge \frac{M\varepsilon }{8r^m}\left( \frac{1+\delta }{\delta }\right) ^{-({\bar{k}}+1)}(r+t)^{-({\bar{k}}+1)}\int _{r-t}^{r+t}\lambda ^{m} d\lambda \\&\ge \frac{M\varepsilon }{8}\left( \frac{1+\delta }{\delta }\right) ^{-({\bar{k}}+1)} \frac{(r-t)^{m}2t}{r^m(r+t)^{{\bar{k}}+1}}\,. \end{aligned}$$

Since \((t,r)\in \Sigma _\delta \), we have

$$\begin{aligned} u(t,r)\ge \frac{C_0t^{m+1}}{r^m(r+t)^{{\bar{k}}+1}}\,, \end{aligned}$$

where we set

$$\begin{aligned} C_0=\varepsilon \frac{2^{m-2}M}{\delta _m^{m}}\left( \frac{\delta }{1+\delta }\right) ^{{\bar{k}}+1}>0. \end{aligned}$$
(18)

Now, we assume an estimate of the form

$$\begin{aligned} u(t,r)\ge \frac{Ct^a}{r^m(r+t)^b} \text { for } (t,r)\in \Sigma _\delta , \end{aligned}$$
(19)

where a, b and C are positive constant. In particular, (19) holds true for \(a=m+1\), \(b={\bar{k}}+1\) and \(C=C_0\).

Being \(g\ge 0\), from (16) we deduce \(u^0\ge 0\). Combining (15) and (19), for \((t,r)\in \Sigma _\delta \), we get

$$\begin{aligned} u(t,r)&\ge \frac{1}{8r^m}\int _0^t d\tau \int _{r-t+\tau }^{r+t-\tau }\frac{\lambda ^m}{(1+\tau )^{\frac{\mu }{2}(p-1)}}|u(\tau , \lambda )|^p d\lambda \nonumber \\&\ge \frac{C^p}{8r^m}\int _0^t \frac{\tau ^{pa}}{(1+\tau )^{\frac{\mu }{2}(p-1)}}d\tau \int _{r-t+\tau }^{r+t-\tau }\lambda ^{m(1-p)}(\lambda +\tau )^{-pb} d\lambda \nonumber \\&\ge \frac{C^p}{8r^m(r+t)^{pb+m(p-1)}}\int _0^t \frac{\tau ^{pa}}{(1+\tau )^{\frac{\mu }{2}(p-1)}}d\tau \int _{r-t+\tau }^{r+t-\tau }d\lambda \nonumber \\&\ge \frac{C^p}{4r^m(r+t)^{pb+m(p-1)}}\int _0^t \frac{(t-\tau )}{(1+\tau )^{\frac{\mu }{2}(p-1)}}\tau ^{pa}d\tau . \end{aligned}$$
(20)

By means of integration by parts, we obtain

$$\begin{aligned}&\quad \int _0^t \frac{(t-\tau )\tau ^{pa}}{(1+\tau )^{\frac{\mu }{2}(p-1)}}d\tau \ge \frac{1}{(1+t)^{\frac{\mu }{2}(p-1)}}\int _0^t (t-\tau )\tau ^{pa}d\tau \\&\quad \ge \frac{1}{(1+t)^{\frac{\mu }{2}(p-1)}} \frac{t^{pa+2}}{(pa+1)(pa+2)}. \end{aligned}$$

While searching a finite lifespan of a solution, it is not restrictive to assume \(t>1\). We have

$$\begin{aligned} \int _0^t \frac{(t-\tau )\tau ^{pa}}{(1+\tau )^{p-1}}d\tau \ge \frac{t^{p(a-\frac{\mu }{2})+2+\frac{\mu }{2}}}{2^{p-1}(pa+1)(pa+2)}. \end{aligned}$$
(21)

Let \((t,r)\in \Sigma _\delta \), from (19)–(21), we can conclude

$$\begin{aligned} u(t,r)\ge \frac{C^*t^{a^*}}{r^m(r+t)^{b^*}} \text { for }(t,r)\in \Sigma _\delta , \end{aligned}$$
(22)

with

$$\begin{aligned} a^*=p\Big (a-\frac{\mu }{2}\Big )+2+\frac{\mu }{2}, b^*=pb+m(p-1), C^*=\frac{(C/2)^p}{2(pa+2)^2}. \end{aligned}$$

Let us define the sequences \(\{a_k\},\,\{b_k\},\,\{C_k\}\) for \(k\in {\mathbb {N}}\) by

$$\begin{aligned}&a_{k+1}=p\Big (a_k-\frac{\mu }{2}\Big )+2+\frac{\mu }{2}, a_1=m+1, \end{aligned}$$
(23)
$$\begin{aligned}&b_{k+1}=pb_k+m(p-1), b_1={\bar{k}}+1, \end{aligned}$$
(24)
$$\begin{aligned}&C_{k+1}=\frac{(C_k/2)^p}{2(pa_k+2)^2}, C_1=C_0, \end{aligned}$$
(25)

where \(C_0\) is defined by (18). Hence, we have

$$\begin{aligned}&a_{k+1}=p^k\left( m+1-\frac{\mu }{2}+\frac{2}{p-1}\right) +\frac{\mu }{2}-\frac{2}{p-1}, \end{aligned}$$
(26)
$$\begin{aligned}&b_{k+1}=p^k({\bar{k}}+1+m)-m, \end{aligned}$$
(27)
$$\begin{aligned}&C_{k+1}\ge K\frac{C_k^p}{p^{2k}} \end{aligned}$$
(28)

for some constant \(K=K(p,\mu , m)>0\) independent of k. The relation (28) implies that for any \(k\ge 1\), it holds

$$\begin{aligned}&C_{k+1}\ge \exp \left( p^{k}\left( \log (C_0)-S_p(k)\right) \right) , \end{aligned}$$
(29)
$$\begin{aligned}&S_p(k)=\Sigma _{j=0}^k d_j, \end{aligned}$$
(30)
$$\begin{aligned}&d_0=0 \text { and } d_j=\frac{j\log (p^2)-\log K}{p^{j}} \text { for } j\ge 1. \end{aligned}$$
(31)

We note that \(d_j>0\) for sufficiently large j. Since \(\lim _{j\rightarrow \infty }d_{j+1}/d_j=1/p,\) the sequence \( S_p(k)\) converges for \(p>1\) by using the ratio criterion for series with positive terms. Hence, there is a positive constant \(S_{p,K}\ge S_p(k)\) for any \(k\in {\mathbb {N}}\), so that

$$\begin{aligned} C_{k+1}\ge \exp (p^{k}(\log (C_0)-S_{p,K})). \end{aligned}$$
(32)

Therefore, by (22), (26)–(29), we obtain

$$\begin{aligned} u(r,t)\ge \frac{(r+t)^m}{r^m t^{-\frac{\mu }{2}+\frac{2}{p-1}}}\exp (p^k J(t,r)), \end{aligned}$$
(33)

where

$$\begin{aligned} J(t,r):=\log (C_0)-S_{p,K}+\Big (m+1-\frac{\mu }{2}+\frac{2}{p-1}\Big )\log t-({\bar{k}}+1+m)\log (r+t). \end{aligned}$$

Thus, if we prove that there exists \((t_0,r_0)\in \Sigma _\delta \) such that \(J(t_0,r_0)>0\), then we can conclude that the solution to (10) blows up in finite time, in fact

$$\begin{aligned} u(t_0,r_0)\rightarrow \infty \text { for } k\rightarrow \infty . \end{aligned}$$

By the definition of \(J=J(t,r)\), we find that \(J(t,r)>0\) if

$$\begin{aligned} \Big (\frac{2}{p-1}-\frac{\mu }{2}-{\bar{k}}\Big )\log t&> \log \Big (\frac{e^{S_{p,K}}}{C_0}\Big (2+\frac{r-t}{t}\Big )^{{\bar{k}}+1+m}\Big ). \end{aligned}$$

In particular, we can take \((t, r)=(t, t+\max \{\frac{2t}{\delta _m},\delta \})\in \Sigma _\delta \); then, it is enough to prove that

$$\begin{aligned} \Big (\frac{2}{p-1}-\frac{\mu }{2}-{\bar{k}}\Big )\log t >\log \Big (\frac{e^{S_{p,K}}}{C_0}\Big (2+\frac{2}{\delta _m}\Big )^{{\bar{k}}+1+m}\Big ). \end{aligned}$$

Now, the crucial assumption (9) comes into play. The coefficient in the left side is positive, and by using (18), we find that \(J(t,r)>0\) provided

$$\begin{aligned} t> C \epsilon ^{-\left( \frac{2}{p-1}-\frac{\mu }{2}-{\bar{k}}\right) ^{-1}}, \end{aligned}$$
(34)

where

$$\begin{aligned} C=\Big (\frac{e^{S_{p,K}}\delta _m^{m}}{2^{m-2}M}\Big (\frac{1+\delta }{\delta }\Big )^{{\bar{k}}+1}\Big (2+\frac{2}{\delta _m}\Big )^ {1+{\bar{k}}+m}\Big )^{\frac{1}{\frac{2}{p-1}-\frac{\mu }{2}-{\bar{k}}}}, \end{aligned}$$

which is positive. As by-product, the inequality (34) gives the lifespan estimate (312) and concludes the proof of Theorem 2.

3.2 Proof of Theorem 1

We start rewriting the Cauchy problem (1) as a nonlinear wave equation with a time-dependent potential. Let \(v=v(t,x)\) be a solution of (1); we define

$$\begin{aligned} u(t,x):=(1+t)^{\frac{\mu }{2}}v(t,x). \end{aligned}$$

Then, the function \(u=u(t,x)\) is a solution of the Cauchy problem

$$\begin{aligned} {\left\{ \begin{array}{ll} u_{tt}-\varDelta u=(1+t)^{-\frac{\mu }{2}(p-1)}{\left| u\right| }^p+(\frac{\mu }{2}(\frac{\mu }{2}-1)-\nu )\frac{u}{(1+t)^2}, \quad &{}t\ge 0,\quad x\in {\mathbb {R}}^n,\\ u(0,x)=0, \\ u_t(0,x)=\varepsilon g(x). \end{array}\right. } \end{aligned}$$
(35)

If g is radial, then u is radial and it satisfies Eqs. (13) and (14) with

$$\begin{aligned} F(t,u)=(1+t)^{-\frac{\mu }{2}(p-1)}{\left| u\right| }^p+\left( \frac{\mu }{2} \left( \frac{\mu }{2}-1\right) -\nu \right) \frac{u}{(1+t)^2}. \end{aligned}$$

Let us fix \(\delta >0\) and use the same notation of the proof of Theorem 2. Since we are assuming \(\frac{\mu }{2}(\frac{\mu }{2}-1)-\nu \ge 0\), by comparison lemma, see [15, Lemma 2.9], we deduce \(u>0\) in \(\Sigma _\delta .\) Then, it holds

$$\begin{aligned} F(t,u)\ge (1+t)^{-\frac{\mu }{2}(p-1)}{\left| u\right| }^p; \end{aligned}$$

hence, by formula (15) in Lemma 1 we still derive the estimate (20). Thus, the proof of Theorem 2 guarantees the result of Theorem 1.