Re: psfgen readpsf error with no more explanation

From: Vermaas, Joshua (Joshua.Vermaas_at_nrel.gov)
Date: Wed May 17 2017 - 16:12:41 CDT

Hi Sadegh,

There is a non-psfgen way of doing this. Take a look at topotools (http://www.ks.uiuc.edu/Research/vmd/plugins/topotools/), specifically the mergemols command in the utility functions. Also, is this a XPLOR formatted psf? I see numbers in an atomtype column, which makes me think this was made in CHARMM. I'm not sure how NAMD will handle this now, but in the past, this wouldn't simulate regardless.

-Josh

On 05/17/2017 03:03 PM, Sadegh Faramarzi Ganjabad wrote:
Dear NAMD users,

I have two sets of psf/pdb files for a protein and a solvent and want to merge them with VMD with the following simple script

package require psfgen
readpsf protein.psf
readpsf water.psf
coordpdb protein.pdb
coordpdb water.pdb
writepsf all.psf
writepdb all.pdb

both psf files can be read with VMD. However, when I put 'readpsf water.psf' it gives the following error

psfgen) reading structure from psf file water.psf
psfgen) error reading atoms from psf file

MOLECULE DESTROYED BY FATAL ERROR! Use resetpsf to start over.

I am not sure what is wrong with 'water.psf', because it has the same format as 'protein.psf". Here is the beginning of the file

PSF

       1 !NTITLE
 REMARKS VMD-generated NAMD/X-Plor PSF structure file

   23592 !NATOM
       1 WT 1 TIP3 OH2 3 -0.834000 15.9994 0
       2 WT 1 TIP3 H1 1 0.417000 1.0080 0
       3 WT 1 TIP3 H2 1 0.417000 1.0080 0
       4 WT 2 TIP3 OH2 3 -0.834000 15.9994 0
       5 WT 2 TIP3 H1 1 0.417000 1.0080 0
       6 WT 2 TIP3 H2 1 0.417000 1.0080 0
       7 WT 3 TIP3 OH2 3 -0.834000 15.9994 0
******************************************************************************

Am I doing something wrong?

Thanks,
Sadegh

This archive was generated by hypermail 2.1.6 : Sun Dec 31 2017 - 23:21:18 CST