I'm trying to simulate the behavior of chromophore molecules incorporated into a polymer under an external electric field.

I've tried to use classical molecular dynamics (GROMACS) but it's unable to reproduce some phenomena such as molecule aggregation due to its dispersion nature.

So I'm looking for a dynamic method/software that is able to reproduce dispersion interactions on the scale of thousands of atoms. I'm personally thinking about CP2K or CPMD but I'm not sure.

  • $\begingroup$ First of all, are you sure the molecules do aggregate in reality? Force fields are coarse methods indeed, but they are not very likely to be qualitatively wrong for dispersion. $\endgroup$
    – wzkchem5
    Jul 3, 2021 at 14:42
  • $\begingroup$ I believe yes, they do aggregate. At least I've seen a couple of papers where some phenomena were attributed to aggregation. $\endgroup$
    – Roman
    Jul 3, 2021 at 15:29
  • $\begingroup$ roma, @wzkchem5: this may be a naive question, but dispersion interactions are typically pairwise (e.g. between 2 atoms), so what is meant by "on the scale of thousands of atoms"? $\endgroup$ Jul 3, 2021 at 16:08
  • 1
    $\begingroup$ What forcefields are you using? I would blame the brand of FF, not the use of a FF. $\endgroup$
    – B. Kelly
    Jul 3, 2021 at 22:47
  • 2
    $\begingroup$ I believe gromacs has an option to turn on long range dispersion corrections (DispCorr=EnerPress). However, according the manual it should only be used for force fields which were parameterized with dispersion corrections (e.g. Amber99). For CHARMM the dispersion should not be used. You can also try the software xTB which has fast semi-empirical methods and a polarizable force-field. $\endgroup$
    – S R Maiti
    Jul 4, 2021 at 3:57