  We compute energy level correlations in weakly disordered metallic grains
using the fermionic replica method. We use the standard sigma-model approach
and show that non--trivial saddle points, which break replica symmetry, must be
included in the calculation to reproduce the oscillatory behavior of the
correlations. We calculate the correlation functions in all three classical
ensembles GOE, GUE and GSE including the finite-dimensional gradient
corrections. Our results coincides with those obtained by the supersymmetric
sigma-model and the semi-classical trace formula.
