From: Jim Phillips (jim_at_ks.uiuc.edu)
Date: Mon Sep 25 2006 - 14:05:27 CDT
NAMD takes bonds seriously, since they are used to derive exclusions.
Your best option is to implement that H-O restraint in TclForces.
On Mon, 25 Sep 2006, Luca Bellucci wrote:
> Hi Jerome,
> Thank you.
> I use a not standard topology. I created my psf-file with psfgen than I
> modified it.
> I want fix distance atoms like this:
> A_Molecule-O---H---O-Molecule_B where O--H and H--O are 1.5 Ang.
> I do not want fix atoms but only two distance with stretching coordinate.
> At this moment I see that there is a "sanity check" in hydrogen bond in
> routine Molecule_*.C of namd source.
> This not permit extra bond for only hydrogen atom...
>> Most likely, there is a way to generate this PSF file with the extra bond,
>> and it will have a higher chance of being consistent. Was the file created
>> using psfgen? Or CHARMM? Are you using a custom topology? In most cases,
>> you should be able to either add the bond directly to the residue topology,
>> or use a patch.
>> On Monday 25 September 2006 09:18, Luca Bellucci wrote:
>>> Hi all,
>>> How add bond in psf_file?
>>> When I try to add bond in psf_file I obtain
>>> this ERROR:
>>> FATAL ERROR: Atom 1300 has bad group size. Check for duplicate bonds.
>>> ------------- Processor 0 Exiting: Called CmiAbort ------------
>>> Reason: FATAL ERROR: Atom 1300 has bad group size. Check for duplicate
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:42:37 CST