Given the molecule, I would go with the HF method, where defining the exact exchange interaction would be relatively less difficult than finding the exchange correlation function. Though it might be a computationally heavy task, it might still be suitable to solve the Slater determinant than the ground state density functional.
Though I might be wrong with this approach.