[gmx-developers] MPI crash with ghost molecules

[Jason de Joannis]

>> simulations it is). It would probably be better to set the masses to 
>> zero, but I was afraid of creating singularities in integration of 
>> the bonded forces. i.e. dv/dt = f/m.
>
> IIRC, atoms with zero mass are supported; they do not move at all.
> But, this may not be what you want, either.
I tried zeroing the masses and invmass. After about 3000 steps I 
started getting LINCS warnings. This happens even if I use 1e-30 
instead of zero for invmass.

-- 
\ Jason de Joannis  \
\ Emory University  \
  \ jdejoan at emory.edu \
   \ 404-402-1332      \