Study of a class of non-polynomial oscillator potentials

We develop a variational method to obtain accurate bounds for the eigenenergies of H = -Delta + V in arbitrary dimensions N>1, where V(r) is the nonpolynomial oscillator potential V(r) = r^2 + lambda r^2/(1+gr^2), lambda in (-infinity,\infinity), g>0. The variational bounds are compared with results previously obtained in the literature. An infinite set of exact solutions is also obtained and used as a source of comparison eigenvalues.