Density functional calculations of planar DNA base-pairs

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.