  Our recent formulation of the analytic and variational Slater-Roothaan (SR)
method, which uses Gaussian basis sets to variationally express the molecular
orbitals, electron density and the one body effective potential of density
functional theory, is reviewed. Variational fitting can be extended to the
resolution of identity method,where variationality then refers to the error in
each two electron integral and not to the total energy. It is proposed that the
appropriate fitting functions be charge neutral and that all ab initio energies
be evaluated using two-center fits of the two-electron integrals. The SR method
has its root in the Slater's Xalpha method and permits an arbitrary scaling of
the Slater-Gaspar-Kohn-Sham exchange-correlation potential around each atom in
the system. Of several ways of choosing the scaling factors (Slater's exchange
parameters), two most obvious are the Hartree-Fock (HF), alpha_HF, values and
the exact atomic, alpha_EA, values. The performance of this simple analytic
model with both sets for atomization energies of G2 set of 148 molecules is
better than the local density approximation or the HF theory, although the
errors in atomization energy are larger than the target chemical accuracy.
  To improve peformance for atomization energies, the SR method is
reparametrized to give atomization energies of 148 molecules to be comparbale
to those obtained by one of the most widely used generalized gradient
approximations. The mean absolute error in ionization potentials of 49 atoms
and molecules is about 0.5 eV and that in bond distances of 27 molecules is
about 0.02 Angstrom. The overall good performance of the computationally
efficient SR method using any reasonable set of alpha values makes it a
promising method for study of large systems.
