We present a theoretical approach for simulating single photon double core ionization of small molecules in the dipolar approximation. Our method couple a selected CI computation of wave functions to reduce the number of slater determinants needed to converge the peak relative energies and the use of non-orthogonal molecular orbital basis set for the neutral and the ionised system in order to speed up the convergence of the peak intensities.