Buckingham + 18-6 Lennard-Jones

GROMACS version: 2021
GROMACS modification: No


I would like to try some annealing on amorphos silica; in order to do so I either need a force field with triplets potential (which I rekon would be quite inefficient) or a pair potential that combines Buckingham and 18-6 Lennard-Jones.

Is it possible to do the latter in GROMACS, without recurring to tabulated potentials and/or to a 12-6 LJ approximation? I see I can define Buckingham parameters and combination rules, but I am not sure if (a) I can combine it with Lennard-Jones and (b) if I can define arbitrary powers for the repulsive part of LJ.

The answer is to both question is no. The alternative is to use tabulated potential, but you have to down-grade to version 2018.
The tabulated potentials are not active from 2019 and they will probably come back in the future.

Kind regards


I will try with an approximation then