Use direct method of lyapunov to determine stability of a system I am working a problem to determine the stability of the given system:
$\ddot{q} + \alpha\dot{q} + \beta q^{3}$
$\alpha \neq 0$
$\beta \neq 0$
With the given candidate lyapunov function:
$V(q,\dot{q}) = \frac{\dot{q^{2}}}{2} + \frac{\beta q^{4}}{4} $
I have worked out the problem but I am having difficulty understanding my results. 
From the state space representation
$
f(x) = 
\begin{bmatrix}
    x_{2} \\
    -\alpha x_2 - \beta x_1^{3} \\
\end{bmatrix}
= 0
$
There is an equilibrium point at $x_e =[0\ 0]^T$
I checked the suitability of the candidate lyapunov function by checking 
$[\partial xV(x)]\cdot f(x) \leq  0 $ for all $ x \in D $
Which reduced down to:
$-\alpha x_{2}^{2} \leq 0 $
Which holds for all x, provided $\alpha > 0$, so the function is fine as a candidate. 
I then checked if $x_{e}$ was locally stable using the method below:
$V(x) > V(x_{e}) $ for all $x \in D \setminus \{x_e\}$  
$\frac{\dot{q^{2}}}{2} + \frac{\beta q^{4}}{4} > 0$
Which the lyapunov function satisfies this condition.
Then for local asymptotic stability I looked at the following condition
$\partial{x}V(x)f(x) < 0$ for all $x \in D \setminus \{x_e\}$
$-\alpha x_{2}^{2} < 0$
Which led me to believe that the system was locally asymptotically stable, but when I used the indirect method of lyapunov the system linearized about the equilibrium was not Hurwitz.
Clearly I am missing some understanding here.
 A: Given your system as
$$
\begin{split}
\dot{x}_1 &= x_2 \\
\dot{x}_2 &= -\alpha x_2 - \beta x_1^3 \,.
\end{split}
$$
Using the indirect method of linearization around the equilibrium $\textbf{x}_e = \textbf{0}$, we get
$$
\dot{\textbf{x}} = A(\textbf{x}) |_{\textbf{x} = \textbf{0}} \textbf{x} = \begin{bmatrix}
0 & 1 \\
-3\beta x_1^2 & -\alpha 
\end{bmatrix}\Bigg\rvert_{\textbf{x} = \textbf{0}} \textbf{x} = \begin{bmatrix}
0 & 1 \\
0 & -\alpha \\
\end{bmatrix} \textbf{x}\,.
$$
Notice that the eigenvalues of $A(\textbf{0})$ are located at $\lambda_1 = 0$, $\lambda_2 = -\alpha$.
This means we have the special case here that one eigenvalue lies on the imaginary axis. This means that no conclusion about stability can be drawn  using the indirect method. The system could be either stable or unstable. 
This is the classical case where investigating stability by looking at the linearization of the system fails.
To check for stability, you can use the direct method as you already did. Notice however that the function you used only shows stability, not asymptotic stability since $\dot{V}(\textbf{x})$ is only negative semi definite.
However, you can use LaSalle's theorem to show the system is actually asymptotically stable (see the comment by @Kwin van der Veen). Let
$$
S=\{\mathbf{x}\,| \, \dot{V}(\mathbf{x}) = 0 \} \,.
$$
Since for any $\mathbf{x} \in S$ with $\mathbf{x} \neq \mathbf{0}$ we get $\dot{x}_2 \neq 0$, no trajectory other than the trivial one can stay in $S$. Thus, by LaSalle's theorem, the origin is asymptotically stable.
