State-of-the-art models used for drift-diffusion simulations of organic bulk heterojunction solar cells based on band transport are not capable of reproducing the voltage dependence of dark current density and carrier concentration of such devices, as determined by current-voltage and charge-extraction measurements. Here, we show how to correctly reproduce this experimental data by including an exponential tail of localized states into the density of states for both electrons and holes, and allowing recombination to occur between free charge carriers and charge carriers trapped in these states. When this recombination via tail states is included, the dependence of charge-carrier concentration on voltage is distinctly different from the case of band-to-band recombination and the dependence of recombination current on carrier concentration to a power higher than 2 can be explained.