  We discuss the stability properties of the solutions of the general nonlinear
Schroedinger equation (NLSE) in 1+1 dimensions in an external potential
derivable from a parity-time (PT) symmetric superpotential $W(x)$ that we
considered earlier [Kevrekedis et al Phys. Rev. E 92, 042901 (2015)]. In
particular we consider the nonlinear partial differential equation $\{ i
\partial_t + \partial_x^2 - V^{-}(x) +| \psi(x,t) |^{2\kappa} \} \, \psi(x,t) =
0$, for arbitrary nonlinearity parameter $\kappa$. We study the bound state
solutions when $V^{-}(x) = (1/4- b^2)$ sech$^2(x)$, which can be derived from
two different superpotentials $W(x)$, one of which is complex and $PT$
symmetric. Using Derrick's theorem, as well as a time dependent variational
approximation, we derive exact analytic results for the domain of stability of
the trapped solution as a function of the depth $b^2$ of the external
potential. We compare the regime of stability found from these analytic
approaches with a numerical linear stability analysis using a variant of the
Vakhitov-Kolokolov (V-K) stability criterion. The numerical results of applying
the V-K condition give the same answer for the domain of stability as the
analytic result obtained from applying Derrick's theorem. Our main result is
that for $\kappa>2$ a new regime of stability for the exact solutions appears
as long as $b > b_{crit}$, where $b_{crit}$ is a function of the nonlinearity
parameter $\kappa$. In the absence of the potential the related solitary wave
solutions of the NLSE are unstable for $\kappa>2$.
