PM: (x_2=1, x_3=0). The corresponding eigenvalues are -2, -1, and thus this point is a local sink, i.e. stable point

P¯M: (x_2=-1, x_3=0). The corresponding eigenvalues are 2, 1, and thus this point is a local source, i.e. unstable point

PM: (x_2=1, x_3=0, y_1=0). The corresponding eigenvalues are 2, -1, -1, and therefore this point is saddle

P¯M: (x_2=-1, x_3=0, y_1=0). The corresponding eigenvalues are -2, 1, 1, and thus this point is saddle

PM: (x_2=1, x_3=0, y_1=0, \lambda=\lambda_{\rm c}). The corresponding eigenvalues are 2, -1, -1, 0

P¯M: (x_2=-1, x_3=0, y_1=0, \lambda=\lambda_{\rm c}). The corresponding eigenvalues are -2, 1, 1, 0

Q(λ): (x_2, x_3, y_1, \lambda)= \left(\frac{{\lambda^\star}^2}{\sqrt{{\lambda^\star}^4-4}}, \frac{2}{\sqrt{{\lambda^\star}^4-4}}, \frac{\sqrt{2} \lambda ^\star}{\sqrt{{\lambda^\star}^4-4}}, \lambda ^\star\right)

Q¯(λ): (x_2, x_3, y_1, \lambda)= \left(-\frac{{\lambda^\star}^2}{\sqrt{{\lambda^\star}^4-4}}, -\frac{2}{\sqrt{{\lambda^\star}^4-4}}, -\frac{\sqrt{2} \lambda ^\star}{\sqrt{{\lambda^\star}^4-4}}, \lambda ^\star\right)

PH: \left(0, 1, 0, \lambda_{\rm c}\right) with \lambda_{\rm c}\in\mathbb{R}. The eigenvalues are 2, 2, 0, 0, thus the point is nonhyperbolic. The center manifold of PH is the 2-dimensional surface defined by the boundary of M

X_2+X_3-\sqrt{1-X_2^2}=0, X_2^2+X_3^2+Y_1^2=1, \lambda\in\mathbb{R} and the critical point is unstable (see the full analysis in appendix, A.5: #cqgaa70fcapp1-5). This point corresponds to the horizon as analyzed previously

PS: \left(\frac{2}{\sqrt{5}}, -\frac{1}{\sqrt{5}}, 0, \lambda_{\rm c}\right) with \lambda_{\rm c}\in\mathbb{R}. The eigenvalues are \frac{6}{\sqrt{5}}, \frac{2}{\sqrt{5}}, 0, 0, thus the point is nonhyperbolic. The center manifold of PS is the same as for PH and it is also unstable (see the full analysis in appendix, A.6: #cqgaa70fcapp1-6). As seen, in the massless case, this point corresponds to a singularity where A(r)=B(r)\simeq r^{-1}

P(λ): \left(\sqrt{\frac{2}{{\lambda ^\star}^2+4}}, \sqrt{\frac{2}{{\lambda ^\star}^2+4}}, \frac{\lambda ^\star}{\sqrt{{\lambda ^\star}^2+4}}, \lambda ^\star\right). It exists for f(\lambda ^\star)=0, {\lambda ^\star}^2\geqslant 6. The eigenvalues are \frac{{\lambda ^\star}^2-6}{\sqrt{2({\lambda ^\star}^2+4)}}, \frac{{\lambda ^\star}^2-6}{\sqrt{2({\lambda ^\star}^2+4)}}, \frac{\sqrt{2} \left({\lambda ^\star}^2-2\right)}{\sqrt{{\lambda ^\star}^2+4}}, -\frac{\sqrt{2} \lambda ^\star f^{\prime}(\lambda ^\star)}{\sqrt{{\lambda ^\star}^2+4}}. It is unstable for f^{\prime}(\lambda^\star)<0, \lambda^\star>\sqrt{6}, or f^{\prime}(\lambda^\star)>0, \lambda^\star<-\sqrt{6}, and saddle otherwise. This point has already been studied in the case of an exponential potential

C(λ): is given by the conditions f(\lambda ^\star)=0, 2X_2^2+2X_2X_3+X_3^2=1, X_2^2+X_3^2+Y_1^2=1. This line has already been studied in the exponential potential case where it was represented by the red curve in figure, 3: #cqgaa70fcf03. It is the boundary of M and exists only if f(\lambda) has roots. Notice that for a given root of f(\lambda), this line corresponds to the center manifold of PH and PS. By using the same analysis than the exponential case, we found that the line is unstable for f^{\prime}(\lambda ^\star)<0, \lambda ^\star<\frac{2\cos{\theta}+2\sin\theta}{\sqrt{-1-\cos(2\theta)+2\sin(2\theta)}} and Y1 > 0 or f^{\prime}(\lambda ^\star)>0, \lambda ^\star>-\frac{2\cos{\theta}+2\sin\theta}{\sqrt{-1-\cos(2\theta)+2\sin(2\theta)}} and Y1 < 0 while \arccos(2/\sqrt{5}) < \theta < \pi/2

PAdS: \left(\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0, 0\right). This point belongs to the phase space (, 60: #cqgaa70fceqn092), that is, it exists only for negative potential V(\psi) < 0. The eigenvalues are -\frac{3}{\sqrt{2}}, -\sqrt{2}, -\frac{\sqrt{9-12 f(0)}+3}{2 \sqrt{2}}, -\frac{3-\sqrt{9-12 f(0)}}{2 \sqrt{2}}. It is a sink for f(0) > 0 and corresponds to AdS spacetime. Thus, the solutions are asymptotically AdS if we impose the condition V(\psi) < 0 and f(0) > 0 (see appendix, A.7: #cqgaa70fcapp1-7 for complete analysis)

M,M¯: 2X_2^2+2X_2X_3+X_3^2=1, v=0. These surfaces are the two boundaries encountered in the exponential case and they correspond to V(\psi)=0. The eigenvalues are 0, 0, -\sqrt{2} Y_1, \sqrt{2} Y_1 (3 \beta_{1}-\beta_{2}), where \beta_1=\lim_{v\rightarrow 0} v^2 f\left(\frac{1}{v}\right), \beta_2=\lim_{v\rightarrow 0} v f^{\prime}\left(\frac{1}{v}\right)\!. These surfaces are saddle for 3 \beta_{1}-\beta_{2}>0. For 3 \beta_{1}-\beta_{2}<0 and depending on the sign of Y1, the surface will be a 2D unstable manifold (respectively, stable manifold) for Y1 < 0 (respectively, Y1 > 0)

P(∞): X_2=0, X_3=0, Y_1=1, v=0. It is the endpoint of the line P(\lambda):=\left(\frac{\sqrt{2}}{\sqrt{\lambda ^2+4}}, \frac{\sqrt{2}}{\sqrt{\lambda ^2+4}}, \frac{\lambda }{\sqrt{\lambda ^2+4}}, \frac{1}{\lambda }\right) as \lambda\rightarrow \infty. The eigenvalues are \sqrt{2}, \frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, \sqrt{2} (3 \beta_{1}-\beta_{2}). This point is unstable if 3 \beta_{1}-\beta_{2}>0, or saddle otherwise

P¯(∞): X_2=0, X_3=0, Y_1=-1, v=0. It is the endpoint of the line \bar P(\lambda):=\left(-\frac{\sqrt{2}}{\sqrt{\lambda ^2+4}}, -\frac{\sqrt{2}}{\sqrt{\lambda ^2+4}}, -\frac{\lambda }{\sqrt{\lambda ^2+4}}, \frac{1}{\lambda }\right) as \lambda\rightarrow \infty. The eigenvalues are -\sqrt{2}, -\frac{1}{\sqrt{2}},-\frac{1}{\sqrt{2}}, -\sqrt{2} (3 \beta_{1}-\beta_{2}). This point is stable if 3 \beta_{1}-\beta_{2}>0, or saddle otherwise