Dear all,
I am doing a coarse-grained simulation on the self-assembly of specific peptides in explicit solvent.
I have used an existing model for coarse-graining amino acids and to find the Lennard Jones parameters for the CG beads. In this model, the cut-off radius for the LJ potential is defined to be 15 angstroms (four times the smallest bead).
Considering that in the simulations, the cut-off radius for the LJ is larger than 2^(1/6)*sigma, is it still suggested to use r_cut = 2^(1/6)*sigma for DPD?