Re: Custom CHARMM forcefield usage problem

From: Jernej Zidar (jernej.zidar_at_gmail.com)
Date: Tue Aug 09 2011 - 23:47:30 CDT

I think I have solved the problem.

The PSFGEN utility in VMD apparently failed to process the
toppar_all27_lipid_cholesterol.str correctly.

After recreating the PSF/PDB files within the Charmm program and
taking care to write an xplor PSF file everything appears to be in
order.

Best regards,
Jernej Zidar

On Wed, Aug 10, 2011 at 12:12, Jernej Zidar <jernej.zidar_at_gmail.com> wrote:
> Hi.
>  Years ago we developed a CHARMM force field for lipids, that works
> flawlessly within the program CHARMM. When I try to use it in NAMD
> though NAMD complains about:
> FATAL ERROR: UNABLE TO FIND BOND PARAMETERS FOR CEL1 CTL1 (ATOMS 8 39)
>
>  In Charmm I would generally set the input like so:
> !read topology
> open read card unit 10 name top_all27_lipid.rtf
> read rtf card unit 10
> close unit 10
> !read parameters
> open read card unit 10 name par_all27_lipid.prm
> read param card unit 10
> close unit 10
>
> stream toppar_all27_lipid_cholesterol.str
>
>  In NAMD I'm using only the par_all27_lipid.prm file for NAMD
> complains if I try to use the stream
> toppar_all27_lipid_cholesterol.str in the conf file:
> parameters          par_all27_lipid.prm
> parameters          toppar_all27_lipid_cholesterol.str
>
>  The error I get is:
> FATAL ERROR: UNABLE TO FIND ANGLE PARAMETERS FOR HT HT OT (ATOMS 44641
> 44643 44642)
>
>  I created the PSF and PDB files for the simulation using the PSFGEN
> tool within VMD and using the aforementioned force field files.
>
>  Any tips on how to solve this will be much appreciated.
>
> Thanks,
> Jernej Zidar
>
>

-- 
Windows: Re-Boot, Linux: Be-Root.

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:57:33 CST