Hi Mike,
Yes, rigid bonds require a feature, called BOND_CONSTRAINT. Make sure you uncomment this before recompiling the code. I'm not very familiar with rigid bonds, but as far as I know they use an algorithm that converges the particle positions in order to satisfy the rigid bond constraint (i.e. how far should the two atoms be). This requires a positional tolerance that will stop the algorithm once the distance between the two atoms is within that tolerance. Rigid bonds also give rise to velocity constraints (the velocity of p1 and the velocity of p2 aren't independent), and require a velocity tolerance.
Best,
Tristan
On Jun 8, 2010, at 11:02 , Mikheil Azatov wrote:
Hi again,
Sorry I've been asking a lot of questions lately, just started using espresso and a lot of questions appear. I tried looking through the mailist archive but I don't think this one was asked.
In the manual there is a rigid_bond bonded interaction but it's without a description. I tried using it but i receive "unknown bonded interaction type "rigid_bond" ". So I was wondering if anyone has used rigid_bond before then what features I might need to enable to make it work? Also what exactly are positional tolerance and velocity tolerance in the definition of rigid_bond.
Thanks again for you help!
Mike