I try to run a DPD simulation with 4 types of beads. There is an additional FENE potential for the beads of type 1 to model a polymer chain. (See input file)
The first step in the calculation consisted in relaxing the constrains due to the construction stage (I use chain.f given with the LAMMPS distribution). I use a pair_style soft potential and it works well.
Then I change the potential to DPD style (pair_style dpd 1 2.5 34567) and my calculation seems to be stuck. There is no error message or warning but nothing more happens.
I run this simulation in parallel (10 processors) using mpiexec and OpenPBS.
Do you have already met this kind of problem ? If yes, can you help me ?
I thank you ahead for your answers.
input2.dat (2.68 KB)