A variational calculation of the energy levels of the class of PJ-invariant
quantum mechanical models described by the non-Hermitian Hamiltonian H = p
(2) - (ix)(N) with N positive and x complex is presented. The energy levels
are determined by finding the stationary points of the functional [H](a,b,
c) = (integral(C)dx psi(x) H psi(x))/(integral(C)dx psi(2)(x)), where psi(x
) = (ix)(c)exp(a(ix)(b)) is a three-parameter class of PJ-invariant trial w
ave functions. The integration contour C used to define [H](a,b,c) lies ins
ide a wedge in the complex-x plane in which the wave function falls off exp
onentially at infinity. Rather than having a local minimum the functional h
as a saddle point in the three-parameter (a,b,c)-space. At this saddle poin
t the numerical prediction for the ground-state energy is extremely accurat
e for a wide range of N. The methods of supersymmetric quantum mechanics ar
e used to determine approximate wave functions and energy eigenvalues of th
e excited states of this class of non-Hermitian Hamiltonians. (C) 1999 Publ
ished by Elsevier Science B.V. All rights reserved.