From: John Stone (johns_at_ks.uiuc.edu)
Date: Wed Apr 17 2013 - 20:48:35 CDT

Hi guys,
  I think there's some confusion in this discussion.
The "animate" command is a VMD command (not psfgen), and the
only data that it saves is whatever data is in the VMD molecule you're
writing from. This is not the same as the structure that one builds
in psfgen. If you're building a molecule in psfgen, then you have
to use the psfgen "writepsf" and "writepdb" commands, and you
would not use any of the VMD commands to write a structure built inside
of psfgen.

Cheers,
  John Stone
  johns_at_ks.uiuc.edu

On Wed, Apr 17, 2013 at 06:32:56PM -0500, Josh Vermaas wrote:
> Hmm, you appear to have stumbled into a bug with how psfs are written, as
> my tests show that it fails to write them out too if I use guessbonds. If
> you are just trying to delete waters, the psfgen userguide
> (www.ks.uiuc.edu/Research/vmd/plugins/psfgen/ug.pdf) has a handy little
> section devoted to this very topic, and you'd save yourself quite a bit of
> heartache until someone fixes the psf writer.
> -Josh
>
> On 04/17/2013 05:52 PM, Anurag Sharma wrote:
>
> Josh,
>
> I have tried the topo guessbonds command. It gave me the right no. Of bonds which I verified with topo numbonds but when I used animate write psf command and generate psf file. There are no bonds mentioned in that psf file ( all other parameters are present). That means the problem I am facing is in animate command.
> I am deleting water molecules from inside the tube here.
>
> Thanks
>
> Anurag Sharma
>
> Sent from my iPhone
>
> On Apr 17, 2013, at 5:42 PM, "Josh Vermaas" <vermaas2_at_illinois.edu> wrote:
>
>
> If you are satisfied with the number of other terms present in the psf, all you'd need to do is add a line like:
> topo guessbonds
> and then write the psf file with animate. I'd be very careful though, as psfs are easy to screw up unintentionally when building them by hand. As others have pointed out on this list, there is a perfectly serviceable nanotube builder (http://www.ks.uiuc.edu/Research/vmd/plugins/nanotube/) and the usual solvation schemes (http://www.ks.uiuc.edu/Research/vmd/plugins/solvate/) will get you a nanotube in water without needing to rebuild bond tables.
> -Josh
>
> On 04/17/2013 02:54 PM, Anurag Sharma wrote:
>
> Hi,
>
> I am working with a solvated Carbon nanotube model. I have deleted some of the water molecules in the model and try to generate psf file for this new model using topo tools. I used animate write psf script to create this psf file. But in the generated psf file the bonds are missing (0 bonds). Though when I used guess bonds and then numbonds in tcl it shows the correct no. What could be the possible solution for this problem. Please suggest.
>
> Angles, dihedrals and impropers are present in the generated psf file. All other parameters are also present. Only bonds are missing.
>
> Thank you
>
> Anurag Sharma

-- 
NIH Center for Macromolecular Modeling and Bioinformatics
Beckman Institute for Advanced Science and Technology
University of Illinois, 405 N. Mathews Ave, Urbana, IL 61801
http://www.ks.uiuc.edu/~johns/           Phone: 217-244-3349
http://www.ks.uiuc.edu/Research/vmd/