  We investigate the charging energy level statistics of disordered interacting
electrons in quantum dots by numerical calculations using the Hartree
approximation. The aim is to obtain a global picture of the statistics as a
function of disorder and interaction strengths. We find Poisson statistics at
very strong disorder, Wigner- Dyson statistics for weak disorder and
interactions, and a Gaussian intermediate regime. These regimes are as expected
from previous studies and fundamental considerations, but we also find
interesting and rather broad crossover regimes. In particular, intermediate
between the Gaussian and Poisson regimes we find a two-sided exponential
distribution for the energy level spacings. In comparing with experiment, we
find that this distribution may be realized in some quantum dots.
