Critical phenomena in gravitational collapse of Husain-Martinez-Nunez scalar field

TL;DR

This paper develops analytical models to study critical phenomena in gravitational collapse of a scalar field, revealing a universal critical exponent and showing that non-self-similar solutions can also exhibit critical collapse features.

Contribution

The authors construct new analytical models for scalar field collapse, demonstrating that critical phenomena occur even in non-self-similar solutions and deriving the critical exponent.

Findings

Critical exponent γ = 0.5 for black hole mass scaling.

Matching solutions at null hypersurfaces ensures metric continuity.

Non-self-similar solutions can exhibit critical collapse features.

Abstract

We construct analytical models to study the critical phenomena in gravitational collapse of the Husain-Martinez-Nunez massless scalar field. We first use the cut-and-paste technique to match the conformally flat solution ($c=0$ ) onto an outgoing Vaidya solution. To guarantee the continuity of the metric and the extrinsic curvature, we prove that the two solutions must be joined at a null hypersurface and the metric function in Vaidya spacetime must satisfy some constraints. We find that the mass of the black hole in the resulting spacetime takes the form $M\propto (p-p^*)^\gamma$, where the critical exponent $\gamma$ is equal to $0.5$. For the case $c\neq 0$, we show that the scalar field must be joined onto two pieces of Vaidya spacetimes to avoid a naked singularity. We also derive the power-law mass formula with $\gamma=0.5$. Compared with previous analytical models constructed from a different scalar field with continuous self-similarity, we obtain the same value of $\gamma$. However, we show that the solution with $c\neq 0$ is not self-similar. Therefore, we provide a rare example that a scalar field without self-similarity also possesses the features of critical collapse.

Equations

$$M\propto (p-p^{\ast }{)}^{\gamma },$$

$$G_{ab}$$

$$T_{ab}$$

$$d{s}^2$$

$$\Phi (r,t)$$

$$\mathcal{R}=\frac{12c{a}^2(r-c)-3a^2{r}^2}{2r^2(at+1{)}^3}{\left(1-\frac{2c}{r}\right)}^{-2-\alpha }+\frac{2c^2(1-{\alpha }^2)}{(at+1)r^4}{\left(1-\frac{2c}{r}\right)}^{-2+\alpha },$$

$${\Theta }_{+}$$

$${\Theta }_{-}$$

$$\frac{\partial }{\partial {\lambda }^{+}}$$

$$\frac{\partial }{\partial {\lambda }^{-}}$$

$${\Theta }_{+}=\left\{a{r}^2{\left(1-\frac{2c}{r}\right)}^{-\alpha }+\left[2r-\frac{2cr(\alpha -1)}{r-2c}\right](at+1)\right\}{r}^{-2}(at+1{)}^{-2},$$

$${\Theta }_{-}=\left\{a{r}^2{\left(1-\frac{2c}{r}\right)}^{-\alpha }-\left[2r-\frac{2cr(\alpha -1)}{r-2c}\right](at+1)\right\}{r}^{-2}(at+1{)}^{-2}.$$

$$\frac{a}{a{t}_{AH}+1}=-\frac{2}{r_{AH}^2}\left[r_{AH}-c(1+\alpha )\right](1-\frac{2c}{r_{AH}}{)}^{\alpha -1}.$$

$$M=\frac{R}{2}\left(1-g^{ab}{\nabla }_{a}R{\nabla }_{b}R\right).$$

$$g^{ab}{\nabla }_{a}R{\nabla }_{b}R=0.$$

$$M$$

$$d{s}_{-}^2=(1+at)\left[-d{t}^2+d{r}^2+r^2\left(d{\theta }^2+\sin {\theta }^{2}d{\varphi }^2\right)\right].$$

$$d{s}_{+}^2=-f_{+}d{U}^2-2dUdR+R^{2}d{\Omega }^2,$$

$$\displaystyle ds_{-}^{2}\Big{|}_{\Sigma}$$

$$\displaystyle ds_{+}^{2}\Big{|}_{\Sigma}$$

$$\displaystyle ds_{-}^{2}\Big{|}_{\Sigma}$$

$$\displaystyle k^{-}_{ab}\Big{|}_{\Sigma}$$

$$k_{ij}=-n_{\mu }\frac{{\partial }^2{x}^{\mu }}{\partial {\xi }^i{\xi }^j}-n_{\nu }{\Gamma }_{\mu \rho }^{\nu }\frac{\partial x^{\mu }}{\partial {\xi }^i}\frac{\partial x^{\rho }}{\partial {\xi }^j}$$

$$k_{\lambda \lambda }^{-}$$

$$k_{\lambda \lambda }^{+}$$

$$k_{\theta \theta }^{-}$$

$$k_{\theta \theta }^{+}$$

$$R$$

$$(at+1)({\dot{t}}^2-{\dot{r}}^2)$$

$$f_{+}\dot{U}+\dot{R}=\frac{ar\dot{r}+2(at+1)\dot{t}}{2\sqrt{at+1}}.$$

$$\dot{R}=\frac{ar\dot{t}}{2\sqrt{(at+1)}}+\sqrt{at+1}\dot{r}.$$

$$\dot{U}$$

$$f_{+}$$

$$m=\frac{a^2{r}^3}{8(at+1{)}^{\frac{3}{2}}}.$$

$$\ddot{U}$$

$$\ddot{R}$$

$$\frac{dm(U)}{dU}$$

Full text

\CJKindent

1 Introduction

Gravitational collapse is the main reason of various galactic structures and it remains one of the most interesting and fundamental problems in general relativity. The end state of gravitational collapse could be a black hole, naked singularity or flat spacetime. In a seminal work by Choptuik [1], some intriguing and universal properties concerning the formation of black holes from massless scalar fields were found. This is called the critical phenomenon. Particularly, near the threshold, the black hole mass can always be expressed in the form of power-law:

[TABLE]

for $p>p^{*}$, where $p$ is a parameter of the initial data to the threshold of black hole formation. Numerical simulations have shown that that the critical exponent $\gamma$ is equal to $0.5$ for solutions with continuous self-similarity (CSS) and $\gamma\approx 0.37$ for solutions with discrete self-similarity (DSS). Details about the critical phenomenon can be learned in [2, 3].

In addition to numerical calculation, analytical models were also built to explore the critical phenomena. Patrick R. Brady[4] studied an exact one parameter family of scalar field solutions which exhibit critical behaviours when black hole forms. J. Soda and K. Hirata[5] analytically studied the collapse of continuous self-similar scalar field in higher dimensional spacetimes and found a general formula for the critical exponents which agrees with the exponent $\gamma=0.5$ for $n=4$. A. Wang et al [6] constructed an analytical model by pasting the BONT model (a massless scalar field) with the Vaidya model. They demonstrated that the black hole mass obeys the power law with $\gamma=0.5$. A. Wang et al [7] also analytically studied the gravitational collapse of a massless scalar field with conformal flatness. They showed that the mass of the black hole without self-similarity turns on with finite nonzero values. Recent developments regarding the critical phenomenon can be found in Refs. [8, 9, 10, 11, 12, 13, 14].

In this paper, we investigate the critical phenomena associated with an exact scalar field solution discovered by Husain-Martinez-Nunez (HMN) [15]. The HMN solution is interesting because it represents a black hole in a FLRW universe[16, 17, 18]. It was pointed out that the HMN solution brings in new phenomenology(S-curve) of apparent horizon[19]. Moreover, the conformally transformed HMN spacetime can be a inhomogeneous vacuum solution in Brans-Dick theory [20, 21, 22]. Following the treatment in [6], we match the HMN solution onto an outgoing Vaidya solution along a null hypersurface. This is to guarantee that the black hole mass is finite. Usually, the hypersurface connecting the two parts of the spacetime is a thin shell, i.e., the extrinsic curvature across the hypersurface is discontinuous. By applying the Darmois-Israel formula [28], one can find the relationship between the jump of the extrinsic curvature and the surface stress-energy tensor of the thin shell. However, by properly choosing the function $m(U)$ in the Vaidya metric, the extrinsic metric can be continuous across the null surface. Therefore, no thin shell forms and the resulting spacetime can be at least $C^{1}$.

After the matching, we calculate the apparent horizon and define the black hole mass as the Komar mass at the intersection of the apparent horizon and the null hypersurface. We first study the case $c=0$ where the HMN solution is conformally flat and has CSS. The mass of the black hole is found in the power-law form $M\propto\sqrt{-a}$ for $a<0$, which means $\gamma=0.5$. When $a=0$, the black hole disappear and the spacetime becomes Minkowski.

The case of $c\neq 0$ is more complicated. This solution has no self-similarity. We still find $M\propto\sqrt{-a}$ for $a<0$. Differing from the case of $c=0$, we show that the limiting spacetime ($a=0$) possesses a naked singularity with non-zero ADM mass. Therefore, there is a mass gap between the black hole spacetime ($a<0$) and the spacetime with naked singularity $a=0$.

The paper is organized as follows. In section 2, we briefly introduce the Husain-Martinez-Nunez (HMN) scalar field solution. In section 3, by using the cut-and-paste method, we match the conformally flat HMN solution ($c=0$ with CSS) onto an outgoing Vaidya spacetime at both timelike and null hypersurfaces to construct an analytical model. It turns out that only the matching along the null hypersurface can guarantee the continuity of the metric and extrinsic curvature across the surface. Then, we use this analytical model to study the critical phenomenon and derive the mass formula. In section 4, we join a general HMN ($a\neq 0,c\neq 0$) with two outgoing Vaidya spacetimes. We show that the mass of the black hole approaches zero for $a<0$. We also find that the critical spacetime ($a=0$) possesses a naked singularity with nonzero ADM mass. Concluding remarks are given in Section 5. In Appendix A, we prove that the HMN solution has CSS only when and $a\neq 0$ and $c=0$.

2 Husain-Martinez-Nunez (HMN) spacetime

The Husain-Martinez-Nunez spacetime [15] satisfies the Einstein-scalar field equations

[TABLE]

The spherically symmetric solution is given by 111 Without loss of generality, we have set $b=1$ in the original metric in [15].

[TABLE]

where $\alpha=\pm\frac{\sqrt{3}}{2}$. From the Ricci scalar

[TABLE]

we see that the curvature singularities are located at $r=2c$ (timelike singularity) and $t=-1/a$(spacelike singularity). Using $\Theta_{\pm}$ to label the expansion of null geodesics, we have

[TABLE]

where $\sqrt{h}=R^{2}\sin\theta=(at+1)r^{2}\left(1-\frac{2c}{r}\right)^{1-\alpha}\sin\theta$, and $\lambda_{\pm}$ is the affine parameter of the null geodesics. The tangent to the null geodesic with affine parameter is

[TABLE]

We can get

[TABLE]

The apparent horizon satisfies $\Theta_{+}=0,\Theta_{-}<0$, which is located at

[TABLE]

The detailed analyses about the HMN sapcetime can be found in [18]. We shall focus on the case $-\infty\leq t\leq-\frac{1}{a}$ and $a<0$ because it corresponds to a black hole solution. When $a\neq 0$ and $\alpha=\frac{\sqrt{3}}{2}$, the apparent horizon in the HMN spacetime has "S-curve" shape. However, this does not make differences in our results. Therefore, in this paper, we just study the case $\alpha=-\frac{\sqrt{3}}{2}$ for simplicity. The results remain true when $\alpha=\frac{\sqrt{3}}{2}$. The MisnerSharp mass[34] is defined by

[TABLE]

where $R$ denotes the areal radius. From [19], we know that in spherically symmetric spacetimes, the apparent horizon satisfies

[TABLE]

Therefore, on the apparent horizon, the Misner-sharp mass becomes

[TABLE]

3 Critical behaviour of HMN scalar field with conformal flatness ($c\neq 0$)

3.1 Matching at a timelike boundary

To study the critical phenomenon of HMN massless scalar field, we start with the simple case $c=0$, where the spacetime is conformally flat[15]. First, we need to join the HMN solution with an outgoing Vaidya solution such that the resulting spacetime is asymptotically flat. In this section, we assume that the boundary connecting the two solutions is a timelike hypersurface. We shall use "$-$" to label the inner HMN spacetime and "$+$" to label the exterior Vaidya spacetime (see Fig. 1).

The interior spacetime is described by Eq. (4) for $c=0$

[TABLE]

The exterior spacetime is described by the outgoing Vaidya metric

[TABLE]

where $f_{+}=1-\frac{2m(U)}{R}$. We choose $\xi^{i}=\{\lambda,\theta,\phi\}$ as the intrinsic coordinates on the hypersurface. $\Sigma$ is determined by functions $t(\lambda),r(\lambda)$ from the interior and $T(\lambda),R(\lambda)$ from the exterior. The induced metric on $\Sigma$ is given by

[TABLE]

Here the do“.” means the derivative with respect to $\lambda$.

We use the Darmois junction conditions to match the solutions across $\Sigma$

[TABLE]

where $k_{ab}$ is the extrinsic curvature of $\Sigma$.

Denote the coordinates of the four-dimensional spacetime by $\{x^{\mu}\}$. $\Sigma$ is determined by the functions $\{x^{\mu}(\xi^{i})\}$. Then the components of $k_{ab}$ can be calculated from

[TABLE]

where $n_{a}$ is the spacelike normal to $\Sigma$. Computing $k_{ij}$ from the interior and exterior, respectively, we obtain the nonvanishing components

[TABLE]

Substituting Eqs. (19) and (20) into Eq. (21), we have

[TABLE]

Substituting Eqs. (26) and (27) into Eq. (22), with the help of Eq. (29), we obtain

[TABLE]

Eq. (28) yields

[TABLE]

Therefore, one can solve Eqs. (29)-(31) and obtains

[TABLE]

Thus,

[TABLE]

To proceed, we calculate the following derivatives:

[TABLE]

where $\hat{t}=t+a^{-1}$. Substitute the above results into Eq. (25) and according to Eq. (22),let the right-hand side of Eq. (24) be equal to the right-hand side of Eq. (25). After a lengthy calculation, we obtain the following result, which is surprisingly simple

[TABLE]

Obviously, the solution is $\dot{r}=\dot{t}$ or $\dot{r}=-\dot{t}$. But this means that the hypersurface is null, inconsistent with our assumption. Therefore, we conclude that the two spacetimes cannot be matched through a timelike hypersurface if the continuity of the extrinsic curvature is required.

3.2 Matching at a null hypersurface

Matching the two solutions at a null hypersurface is more complicated than at a timelike hypersurface. We shall follow the method in [6] and [26]. First we use the coordinate transformation $v=t+r$ to replace the coordinate $r$ in Eq. (4) and obtain the metric in the interior

[TABLE]

where

[TABLE]

Let $\Sigma$ be the null hypersurface $v=v_{0}$. The normal to $\Sigma$ is

[TABLE]

where $s$ is a negative arbitrary function such that $n^{a}_{-}$ is a future directed vector. We can introduce a transverse null vector $N_{a}$ by requiring

[TABLE]

Without loss of generality, we assume that $N_{a}^{-}=N_{v}dv_{a}+N_{r}dr_{a}$. Then it is easy to show that

[TABLE]

Now we choose $s=-\dfrac{1}{a(v-r)+1}\dfrac{\partial R}{\partial r}$. Choose $\xi^{i}=\{R,\theta,\phi\}$ to be the intrinsic coordinates on $\Sigma$. Then $\Sigma$ can be determined by

[TABLE]

where $r(R)$ is determined by Eq. (40) with $v=v_{0}$. We define $e_{(a)}^{-\mu}\equiv\dfrac{\partial x^{\mu}_{-}}{\partial\xi^{a}}$ as given in the[6]. Thus

[TABLE]

Similarly to Eq. (23), the transverse extrinsic curvature for the null surface is given by [26].

[TABLE]

By straightforward calculation, we find

[TABLE]

On the other hand, we need to join the HMN solution to the exterior Vaidya spacetime (18) at $\Sigma$, as shown in Fig. 2. Assume that the null surface $\Sigma$ can be described by $U=U_{0}(R)$ from the Vaidya solution. It follows from (18) that

[TABLE]

The spacetime coordinates $\{x_{+}^{\mu}\}$ can be expressed as functions of $\xi^{i}$:

[TABLE]

Define $e_{(i)}^{+\mu}\equiv\dfrac{\partial x^{\mu}_{+}}{\partial\xi^{i}}$ and it is easy to find

[TABLE]

The normal to $\Sigma$ is given by

[TABLE]

where $\beta$ is a negative function which will be determined later. Then the transverse null vector $N^{a}$ in Eq. (42) is

[TABLE]

The continuity condition on $\Sigma$ requires [26]

[TABLE]

This also guarantees that the normal vectors $n^{a}$ defined on both sides are the the same. According to Eqs. (44),(46)-(46), Eqs. (51), (53), we find that the nontrivial equations in Eq. (54) are

[TABLE]

Hence, we get $\beta=-1$.

Now we can calculate the corresponding transverse extrinsic curvature from Eqs. (47) and obtain

[TABLE]

Since $k_{ab}^{+}=k_{ab}^{-}$, Eqs. (48) and (57) give rise to

[TABLE]

By integration, we obtain $m(r)$ as

[TABLE]

where $c_{1}$ is an integral constant. Using k_{\theta\theta}^{+}\big{|}_{\Sigma}=k_{\theta\theta}^{-}\big{|}_{\Sigma}, we can fix $c_{1}$:

[TABLE]

Consequently,

[TABLE]

Note that

[TABLE]

So, Eq. (61) together with Eq. (49) specifies a unique metric function $m(U)$ in the Vaidya solution. Therefore, we have matched the conformally flat spacetime with the Vaidya spacetime at the null hypersurface.

3.3 Mass of the black hole

From Eq. (14), one can calculate the Misner-Sharp mass for the HMN spacetime described by metric (17) and obtain

[TABLE]

Therefore, $m(r)$ in Eq. (61) is just the Misner-Sharp mass at $v=v_{0}$. The apparent horizon determined by Eq. (13) takes the simple form for $c=0$:

[TABLE]

Note that the null surface is determined by

[TABLE]

Eqs. (64) and (65) immediately gives the coordinates at the intersection of $\Sigma$ and the apparent horizon:

[TABLE]

Since $r>0$ and $a<0$, from Eq. (66), we see that the existence of the intersection requires

[TABLE]

Therefore, the Misner-Sharp mass at the intersection is

[TABLE]

As is known, the event horizon coincides with the apparent horizon in the outgoing Vaidya spacetime as shown in Fig. 3. It is also known that the mass function $m$ in the Vaidya metric is constant alone the event horizon [35]. Thus, it is natural to take the mass in Eq. (69) to be the mass of the black hole.

To investigate the critical behavior as $a\rightarrow 0$, we impose the condition that $r_{i}$ in Eq. (67) does not change with $a$. This means that $v_{0}$ must take the form

[TABLE]

where $V$ is a positive constant independent of $a$. Thus, Eq. (69) gives the mass of black hole:

[TABLE]

Eq. (71) shows that the mass of black hole can be put in the form of Eq. (1) and the scaling exponent is $\gamma=0.5$. Obviously, as $a$ approaches zero, the mass of the black hole vanishes and the spacetime becomes Minkowski.

4 Collapse of the general HMN scalar field

In this section, we shall investigate the gravitational collapse associated with a general HMN scalar field ($a\neq 0,c\neq 0$).

4.1 Matching to an outgoing Vaidya solution at a null hypersurface

Under the coordinate transformation

[TABLE]

Eq. (4) can be rewritten in the form

[TABLE]

Here, the function $h(r)$ satisfies

[TABLE]

The areal radius $R$ takes the form

[TABLE]

Similarly to section 3, we match the solution with an outgoing Vaidya solution at the null hypersurface $v=v_{0}$ (see Fig. 4). Substitution of $v=v_{0}$ into Eq. (75) yields the function $r=r(R)$. By the method in section 3.2, the extrinsic curvature can be calculated as

[TABLE]

where

[TABLE]

It is easy to see that $k_{ab}^{+}$ takes the same form as Eq. (57). Then $k_{ab}^{-}=k_{ab}^{+}$ yields

[TABLE]

Thus, by our construction, the metric of the resulting spacetime is continuous and the extrinsic curvature of the null hypersurface is also continuous. Therefore, we have shown that the general HMN spacetime ( $a\neq 0,c\neq 0$) and Vaidya spacetime can be matched at a null hypersurface as showed in Fig. 4.

4.2 Mass of the black hole

By the argument in section 3.3, one can show that $m(r)$ in Eq. (4.1) is exactly the Misner-sharp mass for the metric in Eq. (73). Therefore, from Eqs. (13),(16) and (75), we can obtain the mass on the apparent horizon $m_{AH}(r)$:

[TABLE]

The ingoing null hypersurface boundary is defined by

[TABLE]

From the Eq. (13) and the Eq. (80), we can get the coordinates $(r_{i},t_{i})$ at the intersection of the apparent horizon and the null hypersurface $v=v_{0}$, which satisfies

[TABLE]

Similarly, we choose the null hypersurface which intersects with the apparent horizon at a fixed radius $r_{i}=r_{0}$, i.e., independent of $a$. Again, we take the Misner-Sharp mass at the intersection as the black hole mass. Then, Eq. (79) gives the mass of the black hole

[TABLE]

where

[TABLE]

Eq. (83) shows clearly that the black hole mass satisfies the power law with $\gamma=0.5$. However, the spacetime for $a<0$ possesses a naked singularity $r=2c$ (see Fig. 4), in violation of the cosmic censorship conjecture. To remove the naked singularity, we join another outgoing Vaidya spacetime at $v=v_{1}(v_{1}<v_{0})$, as shown in Fig. 5. No naked singularity exists in this new spacetime.

When we study the relation between the mass and the parameter $a$, we treat $c$ as a constant. We see that $M_{bh}\rightarrow 0$ as $a\rightarrow 0$. When $a=0$, there is no black hole but a naked singularity as shown in section 4.3.

4.3 "Critical spacetime": $a=0$ and $c\neq 0$

Now we study the critical HMN spacetime. For $a=0$, Eq. (4) becomes

[TABLE]

To calculate the apparent horizon of the spacetime, we first choose two families of radial null vector fields

[TABLE]

where $\lambda^{\pm}$ is the affine parameter of the null geodesic.

According to Eq. (11)

[TABLE]

Therefore, the spacetime has no apparent horizon when $r>2c$. According to Eq. (6), the spacetime possesses a naked singularity at $r=2c$ (see Fig. 6).

Since the spacetime is asymptotically flat, we calculate the ADM mass and find

[TABLE]

When $a\neq 0$, the mass of the black hole is given by Eq. (83), which shows clearly that $M\rightarrow 0$ as $a\rightarrow 0$. However, when $a=0$, as we just discussed, the spacetime is not Minkowski and its ADM mass is nonzero. Therefore, there exists a mass gap between the black hole solution ($a\neq 0$) and its limiting spacetime ($a=0$).

5 Conclusion

In this paper, we have used the "cut and paste" method to construct analytical models and study the critical phenomena of the HMN scalar filed. We have shown that the HMN solution with conformal flatness ($c\neq 0$) can be matched with the Vaidya solution along a null hypersurface, but not a timelike hypersurface. We have derived the differential equation which specifies the metric function in the Vaidya solution. For $c\neq 0$, we have joined the scalar field onto two pieces of Vaidya spacetimes to avoid the naked singularity.

We have studied the gravitational collapse for the HMN scalar field and shown that black hole mass satisfies the power law with $\gamma=0.5$. This is consistent with previous results in the literature. When $c\neq 0$, the HMN spacetime has no CSS and the black hole also turns on at infinitely small mass. The result is different from the model in [7], which shows that the formation of black holes may turn on at finite mass when the gravitational collapse has no self-similarity. On the other hand, the mass gap exists between the black hole and the naked singularity during the gravitational collapse of HMN scalar field when $c\neq 0$ as discussed in section 4. Our work suggests that critical collapse can be studied from analytical models which are constructed by known solutions. More models should be investigated in the future in order to test the universal features in gravitational collapse.

Appendix A Self-similarity of HMN spacetime

In this appendix, we will prove that the HMN spacetime is CSS (continuous self similar) only when $c=0$ and $a\neq 0$. A spacetime is continuous self-similar if there exists a conformal Killing vector field $\xi^{a}$ satisfying

[TABLE]

As a result of the spherical symmetry, we can write

[TABLE]

Here, $x$ and $y$ are functions of $r$ and $t$. Substituting this expression into Eq. (91), we find that

[TABLE]

Here,

[TABLE]

From Eq. (93), we can get

[TABLE]

Substituting Eq. (98) into Eq. (94), we obtain

[TABLE]

where $D(t)$ is a integration function of $t$. Putting Eqs. (98) and (99) into Eq. (96), we find

[TABLE]

Obviously, $C_{0}$ must be a constant independent of $r$ and $t$. So the only solution is

[TABLE]

and consequently

[TABLE]

Now Eqs. (98) and (99) become

[TABLE]

Plugging Eqs. (103) and (104) into Eq. (95), we have

[TABLE]

Hence,

[TABLE]

Thus, we have proven that the HMN spacetime is continuous self-similar only for $c=0$ and $a\neq 0$.