  We present a systematic Density Functional Theory (DFT) study of geometries
and energies of the nucleic acid DNA bases (guanine, adenine, cytosine and
thymine) and 30 different DNA base-pairs. We use a recently developed
linear-scaling DFT scheme, which is specially suited for systems with large
numbers of atoms. As a first step towards the study of large DNA systems, in
this work: (i) We establish the reliability of the approximations of our method
(including pseudopotentials and basis sets) for the description of the
hydrogen-bonded base pairs, by comparing our results with those of former
calculations. We show that the interaction energies at Hartree-Fock geometries
are in very good agreement with those of second order M{\o}ller-Plesset (MP2)
perturbation theory (the most accurate technique that can be applied at present
for system of the sizes of the base-pairs). (ii) We perform DFT structural
optimizations for the 30 different DNA base-pairs, only three of which had been
previously studied with DFT. Our results provide information on the effect of
correlation on the structure of the other 27 base pairs, for which only
Hartree-Fock geometries were formerly available.
