Re: [gmx-users] x2top hangs

2009-07-18 Thread Justin A. Lemkul



Scott Milner wrote:

Folks -- I cannot manage to get x2top to work on even the simplest example.
Here is "ethane.pdb", a PDB file for ethane (built using Avogadro, with 
the CONECT records stripped out):


COMPNDUNNAMED
AUTHORGENERATED BY OPEN BABEL 2.2.2b1
ATOM  1  C   LIG 1  -3.466  -1.641  -0.001  1.00  
0.00   C
ATOM  2  C   LIG 1  -1.954  -1.689   0.025  1.00  
0.00   C
ATOM  3  H   LIG 1  -3.817  -0.605   0.026  1.00  
0.00   H
ATOM  4  H   LIG 1  -3.882  -2.168   0.863  1.00  
0.00   H
ATOM  5  H   LIG 1  -3.850  -2.112  -0.910  1.00  
0.00   H
ATOM  6  H   LIG 1  -1.570  -1.218   0.934  1.00  
0.00   H
ATOM  7  H   LIG 1  -1.538  -1.163  -0.839  1.00  
0.00   H
ATOM  8  H   LIG 1  -1.603  -2.726  -0.002  1.00  
0.00   H

END


I try to run the following command under GROMACS 4.0,

x2top -f ethane.pdb -o ethane.top

with the most recent version of the file ffoplsaa.n2t in its proper 
place (/usr/local/gromacs/share/gromacs/top/),

and I get the following output:

...
Opening library file /usr/local/gromacs/share/gromacs/top/aminoacids.dat
WARNING: masses will be determined based on residue and atom names,
this can deviate from the real mass of the atom type
Opening library file /usr/local/gromacs/share/gromacs/top/atommass.dat
Entries in atommass.dat: 178
WARNING: vdwradii will be determined based on residue and atom names,
this can deviate from the real mass of the atom type
Opening library file /usr/local/gromacs/share/gromacs/top/vdwradii.dat
Entries in vdwradii.dat: 28
Opening library file /usr/local/gromacs/share/gromacs/top/dgsolv.dat
Entries in dgsolv.dat: 7
Opening library file /usr/local/gromacs/share/gromacs/top/electroneg.dat
Entries in electroneg.dat: 71
Opening library file /usr/local/gromacs/share/gromacs/top/elements.dat
Entries in elements.dat: 218
Looking whether force field files exist
Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaa.rtp
Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaa.n2t
Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaa.n2t
There are 23 name to type translations
Generating bonds from distances...
atom 0

... and then the program hangs.

Can someone confirm that x2top actually works, and send a simple working 
example to me?




The 4.0.x versions of x2top do not work, as far as I know.  The code is 
undergoing major changes.  Version 3.3.3 should be operable, if memory serves.


-Justin



Scott Milner
William H. Joyce Professor
The Pennsylvania State University
Department of Chemical Engineering
120 Fenske Laboratory
University Park, PA 16802
(814) 863-9355
smil...@engr.psu.edu

___
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use thewww 
interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/mailing_lists/users.php



--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


___
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] x2top hangs

2009-07-18 Thread Scott Milner
Folks -- I cannot manage to get x2top to work on even the simplest  
example.
Here is "ethane.pdb", a PDB file for ethane (built using Avogadro,  
with the CONECT records stripped out):


COMPNDUNNAMED
AUTHORGENERATED BY OPEN BABEL 2.2.2b1
ATOM  1  C   LIG 1  -3.466  -1.641  -0.001  1.00   
0.00   C
ATOM  2  C   LIG 1  -1.954  -1.689   0.025  1.00   
0.00   C
ATOM  3  H   LIG 1  -3.817  -0.605   0.026  1.00   
0.00   H
ATOM  4  H   LIG 1  -3.882  -2.168   0.863  1.00   
0.00   H
ATOM  5  H   LIG 1  -3.850  -2.112  -0.910  1.00   
0.00   H
ATOM  6  H   LIG 1  -1.570  -1.218   0.934  1.00   
0.00   H
ATOM  7  H   LIG 1  -1.538  -1.163  -0.839  1.00   
0.00   H
ATOM  8  H   LIG 1  -1.603  -2.726  -0.002  1.00   
0.00   H

END


I try to run the following command under GROMACS 4.0,

x2top -f ethane.pdb -o ethane.top

with the most recent version of the file ffoplsaa.n2t in its proper  
place (/usr/local/gromacs/share/gromacs/top/),

and I get the following output:

...
Opening library file /usr/local/gromacs/share/gromacs/top/aminoacids.dat
WARNING: masses will be determined based on residue and atom names,
this can deviate from the real mass of the atom type
Opening library file /usr/local/gromacs/share/gromacs/top/atommass.dat
Entries in atommass.dat: 178
WARNING: vdwradii will be determined based on residue and atom names,
this can deviate from the real mass of the atom type
Opening library file /usr/local/gromacs/share/gromacs/top/vdwradii.dat
Entries in vdwradii.dat: 28
Opening library file /usr/local/gromacs/share/gromacs/top/dgsolv.dat
Entries in dgsolv.dat: 7
Opening library file /usr/local/gromacs/share/gromacs/top/electroneg.dat
Entries in electroneg.dat: 71
Opening library file /usr/local/gromacs/share/gromacs/top/elements.dat
Entries in elements.dat: 218
Looking whether force field files exist
Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaa.rtp
Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaa.n2t
Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaa.n2t
There are 23 name to type translations
Generating bonds from distances...
atom 0

... and then the program hangs.

Can someone confirm that x2top actually works, and send a simple  
working example to me?



Scott Milner
William H. Joyce Professor
The Pennsylvania State University
Department of Chemical Engineering
120 Fenske Laboratory
University Park, PA 16802
(814) 863-9355
smil...@engr.psu.edu

___
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] x2top hangs up

2008-09-14 Thread Vitaly Chaban
Hello,

Trying to generate a topology from PDB file my x2top (gromacs 4.0beta) hangs up.
I have the coordinates of carbon nanotube in the PDB file (generated
with MOLEKEL). I also have

ffencadv.n2t with:
C C 1 C ; CNT Carbon with one bond
C C 2 C C ; CNT double bonded Carbon

===
ffgmxbon.itp with:
[ bondtypes ]
; i j func b0 kb
C C 1 0.14210 478900.
[ angletypes ]
; i j k func th0 cth
C C C 1 120.000 397.480
[ dihedraltypes ]
; i l func q0 cq
C C 1 0.000 167.360

x2top -f cnt.pdb -o cnt.top -ff select
Then I select  "8: Encad all-atom force field, using scaled-down vacuum charges"

Opening library file /usr/local/gromacs/share/gromacs/top/ffencadv.rtp
Opening library file /usr/local/gromacs/share/gromacs/top/ffencadv.n2t
Opening library file /usr/local/gromacs/share/gromacs/top/ffencadv.n2t
There are 0 name to type translations
Generating bonds from distances...
atom 0

And the programme hangs up. The processor load is 99% all the time
before I terminate the job.

The same data input works with GROMACS 3.3.1 as described by C. Stiles
but not with 3.3.3 and 4.0b however.
If it is not a bug then what should be modified to make a topology?
Has anyone tried to make SWCNT topology in the versions newer that
3.3.1?

What essential things have been changed in x2top since 3.3.1?

Thanks a lot for suggestions!

-- 
Vitaly V. Chaban
School of Chemistry
National University of Kharkiv
Svoboda sq., 4, Kharkiv 61077, Ukraine
email: [EMAIL PROTECTED]
skype: vvchaban
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php