Check the Forcefield of the electrolyte system EC DEC DMC

GROMACS version: 2018
GROMACS modification: No
I need help checking my electrolyte system’s Forcefield (OPLS-AA). The simulation results don’t match expected values:

EC = Lower coordination number, solvation shell at a slightly different distance. (Green line)
RDF =Deviates from literature.

The simulation box looks stable after the run (no broken bonds).
My electrolyte system is EC:DEC:DMC
How can I pinpoint the problem in the Forcefield?

gmx.top (288.1 KB)


Does this paper, https://pubs.acs.org/doi/10.1021/acs.jpcb.3c05591, help? It seems like lithium parameters vary quite a lot between different force fields. OPLS might not be the most suitable.