It appears reasonable to assume that the broad polaronic absorption observed in experiment is a combination of several (multi-)protonated structures, giving rise to a whole range of different polaronic absorption features. The doping induced absorption feature close to the main peak of PCPDT-BT also observed in the experiment appears to have its origin in protonation at a position for which polaron formation is suppressed (e.g., our case Proton-2). Here we simulate a scenario with different contributions as a superposition of the absorption spectra found in the calculations for different structures as shown in Fig. . With increasing concentration of BCF in the experiment, the proportion of double protonated species over single protonated ones will increase, while the proportion of the neutral PCPDT-BT will decrease, which is represented by the different ratios considered for the calculated spectra shown in Fig. . As can be seen in Fig. (dark blue), a ratio of 20% of each single protonated structure, corresponding to a low BCF concentration, mostly only leads to a reduction of the main peak intensity, with a slight broadening. If the amount of protonation is increased by introducing an additional 10% of each doubly protonated structure (Fig. , dark green), the amplitude of the main peak decreases further and the broad polaronic absorption increases. Note that the amplitude of the polaronic absorption appears to be small, but compared to the reduction of the main peak it is in reasobale agreement with the experimental data. Finally, we address a 20% ratio of each single and double protonated structure (Fig. , dashed dark red), corresponding to a high BCF concentration in the experiment.