Dear all:
I am using the latest lammps version: 16 Mar 2018 to run a Dissipative Particle Dynamics. But when I followed the manual and used a command: fix mvv/dpd, I got a error: Fix mvv/dpd requires atom attribute vest. Then I check the manual, but I didn’t see any solution. And I also search the mailing list archives and also didn’t find a way to slove it. Here is my input file:
units lj
atom_style angle
comm_modify vel yes
boundary p p p
pair_style dpd 1.0 1.0 14915642
bond_style harmonic
special_bonds lj 1.0 1.0 1.0
read_data polymer.data
pair_coeff 1 1 25.000000 4.500000
pair_coeff 1 2 40.000000 4.500000
pair_coeff 2 2 25.000000 4.500000
bond_coeff 1 4.000000 0.000000
bond_coeff 2 4.000000 0.000000
bond_coeff 3 4.000000 0.000000
neighbor 2.0 bin
neigh_modify delay 2
group mobile union all
timestep 0.06
fix 1 mobile nve
fix 2 mobile mvv/dpd 0.65
velocity mobile create 1.0 14915642 mom yes rot yes dist gaussian
dump 1 all atom 10000 polymer.lammpstrj
dump_modify 1 scale no
thermo_style custom step cpu etotal temp
thermo 1000
thermo_modify flush yes
run 1000000
If exclude the fix mvv/dpd, this could be run properly. But when including it, it gave a error.
I will really appreciate for your reply,
best wishes,
Fang.