I am studying the files that allow one to perform BD, however, with a single particle, when I plot r^2 vs time, the graph is non-linear. I was wondering whether anyone has seen this problem before and whether I should use fix viscous, even though the langevin fix is supposed to include the damping term (as well as add random collision force)?
Also, I am slighlty confused about the random number seed – if I understand correctly, this is not a gaussian rn generator? I do not have access to the Dunweg paper:
(Dunweg) Dunweg and Paul, Int J of Modern Physics C, 2, 817-27 (1991).
Does anybody have this article?
my input is as follows:
boundary p p p
atom_modify map array sort 1 1
restart 1000000 particle.restart
group all type 1
dump 1 all custom 100000 particle_ts100_delta1m.dump id x y z ix iy iz
dump_modify 1 format “d .5f .5f .5f %d %d %d”
fix 1 all nve
fix 2 all langevin 1.0 1.0 2.0 424242
thermo_style custom step temp pxx pyy pzz vol
thanks in advance.