Theoretical investigation for the geometric and energetic properties, rotational constants, harmonic vibrational frequencies, and binding energies of nucleic acid base pair, cytosine dimer, are carried out by using the density functional theory method. The dimer structures resulting from both the keto and the enol (cis/trans) tautomers are investigated in the present study. Various isomers are considered to find the stable structures of the cytosine dimer. The planar cytosine dimer, K-K-3 with C-2h, symmetry, resulting from nonplanar keto tautomers, is found to be thermodynamically most stable out of the four different stable isomers and having the highest binding energy value, 19.51 kcal/mol (including basis set superposition error correction). The vibrational frequency analysis also suggests a red shift of 367.97 cm(-1) for the hydrogen-bonding K-K3 symmetric dimer with two hydrogen bond lengths, each of length 1.913 angstrom. Moreover, charge distribution (ChelpG charges), Laplacian electronic density distribution, and the dimerization equilibrium for the most stable dimer, K-K3, have also been investigated using the same method and the basis set.