Dear all, I am trying to equilibrate a solvent of pure ionic liquid. The system keeps exploding (after 2-5 ns) and I am not sure why, though I believe coulombic interactions are to blame. This is because the Coul-SR term is negative, but the Coul. recip term is very positive throughout the entire run (giving the entire system a positive potential energy). I think this means that the short-range electrostatics are okay, but the long range electrostatics (calculated with PME) are not. Does anybody have any suggestions as to why this would happen? I have used the exact same PME input parameters for another ionic liquid that works just fine. They are listed below.
rcoulomb = 1.2 fourierspacing = 0.12 pme_order = 4 ewald_rtol = 1e-5 Thanks! Denzil Frost
-- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists