[lammps-users] fix rigid torque

LAMMPS LAMMPS Mailing List Mirror
1

Hi all.

I’m trying to use the fix rigid command while turning the torque off in all three components, as documented here:
http://lammps.sandia.gov/doc/fix_rigid.html

It seems that no matter what I do, I get an error. I looked at the source file and found that the command only takes four arguments, rather than eight or more as documented.

INPUT
dimension 3
boundary p p p
units metal
timestep 0.002
atom_style atomic
read_data 3x4x4_neon+tether.struct

group TIP type = 1
group BOB type = 2
fix TIP_rigid TIP rigid single force * 0 0 0
fix TIP_tether TIP spring couple BOB 5 1 1 1 14

velocity BOB set 0.5 0.2 0.8 units box sum no

dump dump_xyz all custom 50 xyz.dmp tag x y z type
thermo 50

run 20000

OUTPUT:
LAMMPS (21 May 2008)
Reading data file …
triclinic box = (-13.0484 -17.3979 -17.3979) to (0 0 0) with tilt (0 0 0)
1 by 1 by 1 processor grid
193 atoms
192 atoms in group TIP
1 atoms in group BOB
ERROR: Illegal fix rigid command

2

OUTPUT:
LAMMPS (21 May 2008)

Your copy of LAMMPS is 8 months old and this fix rigid
feature was added last week. The on-line WWW pages
always reflect the most current, fully-patched version
of the code.

Steve