We study static, spherically symmetric, self-gravitating systems minimally coupled to a scalar field with U(1) gauge symmetry: charged boson stars. We find numerical solutions to the EinsteinMaxwell equations coupled to the relativistic Klein-Gordon equation. It is shown that bound stable configurations exist only for values of the coupling constant less than or equal to a certain critical value. The metric coefficients and the relevant physical quantities such as the total mass and charge, turn out to be in general bound functions of the radial coordinate, reaching their maximum values at a critical value of the scalar field at the origin. We discuss the stability problem both from the quantitative and qualitative point of view. Taking properly into account the electromagnetic contribution to the total mass, the stability issue is faced also following an indication of the binding energy per particle, verifying then the existence of configurations having positive binding energy: objects apparently bound can be unstable against small perturbations, in full analogy with what observed in the mass-radius relation of neutron stars.

The existence of electric fields close to their critical value $E_c=\frac{m_e^2 c^3}{e \hbar}$ has been proved for massive cores of $10^7$ up to $10^{57}$ nucleons using a proton distribution of constant density and a sharp step function at its boundary. We explore the modifications of this effect by considering a smoother density profile with a proton distribution fulfilling a Woods--Saxon dependence. The occurrence of a critical field has been confirmed. We discuss how the location of the maximum of the electric field as well as its magnitude is modified by the smoother distribution.