Fixed position structure breaking apart when water is near

GROMACS version:
GROMACS modification: No
Here post your question

Hello all

I am trying to fille a fixed structure, a tube, with water.

I have employed these sigma and epsilon parameters:
0.115817 0.011000

Whenever water reaches the inner surface of the tube the simulation crashes. My hope was to create parameters that repel the water, but also don’t crash the system.

I am unsure why water near the inner surface will make my simulation crash.

So the question is why does my simulation crash when water reaches the inner surface, and how do you fix that problem?


If the simulation is unstable, then the physical interactions are poorly parametrized. What is the source of your LJ parameters?

I have used multiple LJ parameters, some have been imitating other atomtypes(carbon alpha) looking for just a stable combination. One LJ set drawn from a paper had a sigma of 0.04 and an epsilon as 233 that did not produce even a stable initial structure. I have tried just arbitrarily changing each parameter by orders of magnitude and see how that changes the behaviour but I have yet to find a realistic working scenario.

Now you see why it is so difficult to obtain well-balanced force field parameters. You can’t necessarily expect arbitrary values to produce physically meaningful results. There are tons of published carbon nanotube simulations that you should be able to use to draw some inspiration (i.e. try their parameters, because they work in those cases) but if your system is geometrically different, there’s no reason to necessarily expect those parameters will be valid.

1 Like

my next step was carbon nanotubes, thanks for your confirmation it is most helpful!