[lammps-users] Use colloid in DPD

Hi,

I am trying to create colloidal particles that will group of DPD particles sticking together and act as a rigid body. Most of the paper uses a stiff harmonic force to create a bigger particle from a cluster of DPD bead. In lammps , what will be the best command to do that for DPD simulation?

LAMMPS has fix rigid to directly simulate rigid objects without workarounds, which can be mixed with non-rigid objects.
Of course, you can also add bonds, either by preparing a suitable data file or by using the corresponding command.

Thank you Axel for your answer. Does the bond_style harmonic with a high stiffness value work?

“does it work?” is a very unscientific question (about as bad as saying that something “does not work”, which you should have seen me complain about repeatedly, if you are subscribed to this mailing list for a while).

LAMMPS will compute what you ask it to compute and will simulate the model you feed it.
Whether that model is giving you meaningful results is a different question and not exactly a LAMMPS problem, but a problem of the choice of model and correctly applying it.

If so then is there any difference between this two commands?

again, that is not LAMMPS’ concern, but something that you can easily figure out: just set up a representative simulation with each approach and compare the results. or read up on people’s work in the published literature who have done such calculations and discuss the differences.

axel.

Thank you Axel for your answer. Does the bond_style harmonic with a high stiffness value work? If so then is there any difference between this two commands?

Regards,

Tan