Hi LAMMPS users, I am performing particle dynamics simulations to model biological cell and tissue mechanics using the LAMMPS 29 October 2020 version. The cytoplasmic environment of the cell is highly viscous, hence the motion of particles is overdamped. I am using the following command lines in the input script: fix 1 all langevin 300 300 6.20343e-10 236509 fix 2 all bd 6.2034…

I am trying to use a conditional statement to add forces in my simulation box. I am having difficulty implementing the if statements with the fix addforce command. For example, I would like my particle when found between -6.5 and -5 z coordinate to have a force of -20 in the z coordinate Using this if statement is not working for me if “-6.5<z<-5” then " fix 2 addforce 0 0 …

In the Sturat paper the Vtors has half angle of the cos but looking at the source code of airebo.cpp line number 1035 it is given as: cw2 = (.5*(1.0-cw)). But in order to find the half angle square of the cos shouldnt there be positive sign? Or I am missing something. 1 post - 1 participant Read full topic