To analyze the effects of the DFT parameterization in the solution phase, we calculated the free energy profiles (as a function of radial distance) between groups of the oligomers and the Li + ion in classical MD at T = 300 K. Specifically, we examined TBT and 2 TBTs in the implicit solvent model, and TBT and TriTBT in the explicit solvent model.