![Author ORCID: We display the ORCID iD icon alongside authors names on our website to acknowledge that the ORCiD has been authenticated when entered by the user. To view the users ORCiD record click the icon. [opens in a new tab]](https://chemrxiv.org/engage/assets/public/chemrxiv/images/logos/orcid.png)
Abstract
Lithium-sulfur (Li/S) batteries are emerging as a next-generation energy storage technology due to their high theoretical energy density and cost-effectiveness. π-conjugated organosulfur polymers, such as poly(4-(thiophene-3-yl)benzenethiol) (PTBT), have shown promise in overcoming challenges such as the polysulfide shuttle effect by providing a conductive framework and enabling sulfur copolymerization. In these cathodes, cation-π interactions significantly influence Li+ diffusion and storage properties in π-conjugated cathodes, but classical OPLS-AA force fields fail to capture these effects. This study employs a bottom-up approach based on density functional theory (DFT) to optimize the nonbonded interaction parameters (OPLS-AA/corr.), particularly for the Li+-π interactions with the PTBT polymer. Following prior work, we used an ion-induced dipole potential to model the cation-π interactions. The impact of the solvent on the PTBT monomers was examined by computing the potential of mean force (PMF) between PTBT monomers and Li+ ions in both explicit and implicit solvents using the Boltzmann inversion of probability distributions close to room temperature. In the implicit solvent case, the magnitude of the binding free energy decreased with increasing dielectric constant, as the dominant electrostatics scaled with the dielectric constant. In contrast, in the explicit solvent case, considering the mixtures of organic solvent DME and DOL, the binding free energy shows minimal dependence on solvent composition due to the competing interaction of TBT and Li+ with the solvent molecules. However, increasing salt concentration decreases the binding free energy due to Debye-Hückel screening effects. In general, this work suggests that the optimized parameters can be widely used in the simulation of polymers in electrolytes for the Li/S battery to enhance the representation of cation-π interactions for a fixed charge force field.
Supplementary materials
Title
Electronic Supporting Information:\\Optimizing Cation-π Force-fields for Molecular Dynamics Studies of Competitive Solvation in Conjugated Organosulfur Polymers for Lithium-Sulfur Batteries
Description
It includes sensitivity analysis and optimized parameters for benzene, thiophene, benzenethiol, TBT, and TriTBT, presented in Tables S1 to S5, respectively. Additionally, Fig. S1 shows the atom types for benzene, benzenethiol, and thiophene. The 2D representations of binding energy for thiophene–Li+ and benzenethiol–Li+, calculated using DFT and MS with and without correction, are provided in Figs. S2 and S3. Fig. S4 depicts the clustering behavior of LiTFSI salt in implicit solvent for varying salt concentrations. Free energy profiles for TBT and TriTBT at fixed salt concentrations and various solvent compositions for explicit solvent are shown in Figs. S5 and S6. Furthermore, the free energy profiles for TBT and a two-TBT setup at three different dielectric constants are presented in Fig. S7. Finally, Fig. S8 illustrates the free energy profiles for TBT and TriTBT at a fixed salt concentration of 1.5 M and a fixed solvent composition of 0.55.
Actions
