Global conservation constraints obtained from the nonlinear Vlasov-Maxwell equations are used to derive a three-dimensional kinetic stability theorem for an intense non-neutral ion beam (or charge bunch) propagating in the z direction with average axial velocity upsilon(b) = const and characteristic kinetic energy (gamma(b)-1)mc(2) in the laboratory frame. Here, gamma(b) = (1 - upsilon(b)(2)/c(2))(-1/2) is the relativistic mass factor, and a perfectly conducting cylindrical wall is located at radius r=r2(w), where r=(x(2) +y(2))(1/2) is the radial distance from the beam axis. The particle motion in the beam frame (''primed'' coordinates) is assumed to be nonrelativistic, and the beam is assumed to have sufficiently high directed axial velocity that upsilon(b)much greater than\v'\. Space-charge effects and transverse electromagnetic effects are incorporated into the analysis in a fully self-consistent manner. The nonlinear Vlasov-Maxwell equations are Lorentz-transformed to the beam frame, and the applied focusing potential is assumed to have the (time-stationary) form psi(sf)'(x') = (gamma(b)m/2)[omega(beta perpendicular to)(2)(x'(2) + y'(2)) + omega(beta z)(2)z'(2)], where omega(beta perpendicular to) and up, are constant focusing frequencies. It is shown that a sufficient condition for linear and nonlinear stability for perturbations with arbitrary polarization about a beam equilibrium distribution f(eq)(x',p') is that f(eq) be a monotonically decreasing function of the single-particle energy, i.e., partial derivative f(eq)(H')/partial derivative H' less than or equal to 0. Here, H' = p'(2)/2m + psi(sf)'(x') + phi(eq)(x'), where phi(eq)(x') is the space-charge potential. (C) 1998 American Institute of Physics. [S1070-664X(98)01909-0]