You need to build the model first and then you can simulate and optimize the parameter using matlab, octave, etc. FYI, I attach my previous simulation about DSSC based on ZnO photoanode. I hope it will help.
I perfectly agreement with Arman and see COMSOL Multiphysics program. But if you can study the electronic properties of organic compounds implying the GAUSSIAN prog.