We analyze a system of self-gravitating identical bosons by means of a semirelativistic Hamiltonian comprising the relativistic kinetic energies of the involved particles and added (instantaneous) Newtonian gravitational pair potentials. With the help of an improved lower bound to the bottom of the spectrum of this Hamiltonian, we are able to enlarge the known region for relativistic stability for such boson systems against gravitational collapse and to sharpen the predictions for their maximum stable mass.

For two bound-state equations derived as simplified forms of the Bethe-Salpeter equation with confining interaction, stability of all solutions is rigorously shown.

Analytic energy bounds for N-boson systems governed by semirelativistic Hamiltonians of the form H=\sum_{i=1}^N(p_i^2 + m^2)^{1/2} - sum_{1=i<j}^N v/r_{ij}, with v>0, are derived by use of Jacobi relative coordinates. For gravity v=c/N, these bounds are substantially tighter than earlier bounds and they are shown to coincide with known results in the nonrelativistic limit.