  We consider the nonlinear Dirac equation in one dimension, also known as the
Soler model in (1+1) dimensions, or the massive Gross-Neveu model:
$i\partial_t\psi=-i\alpha\partial_x\psi+m\beta\psi-f(\psi^\ast\beta\psi)\beta\psi$,
$\psi(x,t)\in\C^2$, $x\in\R$, $f\in C^\infty(\R)$, $m>0$, where $\alpha$,
$\beta$ are $2\times 2$ hermitian matrices which satisfy $\alpha^2=\beta^2=1$,
$\alpha\beta+\beta\alpha=0$.
  We study the spectral stability of solitary wave solutions
$\phi_\omega(x)e^{-i\omega t}$. More precisely, we study the presence of point
eigenvalues in the spectra of linearizations at solitary waves of arbitrarily
small amplitude, in the limit $\omega\to m$. We prove that if
$f(s)=s^k+O(s^{k+1})$, $k\in\N$, with $k\ge 3$, then one positive and one
negative eigenvalue are present in the spectrum of linearizations at all
solitary waves with $\omega$ sufficiently close to $m$. This shows that all
solitary waves of sufficiently small amplitude are linearly unstable. The
approach is based on applying the Rayleigh-Schr\"odinger perturbation theory to
the nonrelativistic limit of the equation. The results are in formal agreement
with the Vakhitov-Kolokolov stability criterion. Let us mention a similar
independent result [Guan-Gustafson] on linear instability for the nonlinear
Dirac equation in three dimensions, with cubic nonlinearity (this result is
also in formal agreement with the Vakhitov-Kolokolov stability criterion).
