Overall, the optimized force field parameters successfully improved the accuracy of the binding energy profiles at the equilibrium distances, demonstrating the effectiveness of our approach in capturing the cation-π interactions. Given that our model performs well for small molecules, we aim to extend its application to larger molecules such as "Free TBT monomer" and "PTBT constituent monomer".