The MD simulations were carried out on PEDOT-PSS and PEDOT-PSS-sorbitol (10%, 20%, 40%) model systems in a cubic periodic box. PEDOT oligomers with eight repeat units and +2 charge were considered, while the PSS atactic chains consisted of 20 repeat units, with four deprotonated units randomly distributed in the sequence of each chain. To ensure an adequate equilibration of the systems, the following approach was adopted. First, random starting configurations were generated with the PACKMOL program, and each one were allowed to relax with an energy minimization followed by a molecular dynamics simulation of 1 ns within the NVT ensemble at 300K. After that, an NPT simulation of 1ns duration at a temperature of 300K and a pressure of 1 atm was performed to relax the dimensions of the box. Then, the temperature of the system was progressively increased to 700K and decreased again to 300K before the production run of 5ns in the NPT ensemble at 300K and 1 atm. Using that procedure, we generated 8 independent morphologies for each of the sorbitol concentrations using different random starting configurations.