Donor−acceptor heterojunctions composed of thiophene oligomers and C 60 fullerene were investigated with computational methods. Benchmark calculations were performed with time-dependent density functional theory. The effects of varying the density functional, the number of oligomers, the intermolecular distance, the medium polarization, and the chemical functionalization of the monomers were analyzed. The results are presented in terms of diagrams where the electronic states are classified as locally excited states, charge-transfer states, and delocalized states. The effects of each option for computational simulations of realistic heterojunctions employed in photovoltaic devices are evaluated and discussed.