<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
OK, including a mass colum, as&nbsp;Justin A. Lemkul suggested, worked well. Now grompp does not give errors.<BR>
Why should I not use this model for urea? I have found also another model in the contributions section, for a&nbsp;urea/water box 10M, but why should I trust more in a user contribution than in a *itp included in Gromacs?<BR>
I don´t have much experience with Gromacs, so any coment will be wellcome. <BR>
Thank you very much for your help and suggestions.<BR>
Best wishes,<BR>
&nbsp;<BR>
Rebeca García <BR>
Parc Cientific de Barcelona<BR>
<A href="mailto:regafan@hotmail.com">regafan@hotmail.com</A><BR><BR>&gt; Date: Thu, 18 Dec 2008 19:50:03 +0100<BR>&gt; From: spoel@xray.bmc.uu.se<BR>&gt; To: gmx-users@gromacs.org<BR>&gt; Subject: Re: [gmx-users] Urea topology problem<BR>&gt; <BR>&gt; Rebeca García Fandiño wrote:<BR>&gt; &gt; Hello,<BR>&gt; &gt; I don´t understand the correction I should do from itp.<BR>&gt; &gt; I have removed from the urea original urea.itp these lines,<BR>&gt; <BR>&gt; don't use this model. do a proper literature search.<BR>&gt; <BR>&gt; &gt; <BR>&gt; &gt; #ifdef Boek<BR>&gt; &gt; 1 C 1 UREA C1 1 0.38<BR>&gt; &gt; 2 O 1 UREA O2 1 -0.38<BR>&gt; &gt; 3 NT 1 UREA N3 2 -0.83<BR>&gt; &gt; 4 H 1 UREA H4 2 0.415<BR>&gt; &gt; 5 H 1 UREA H5 2 0.415<BR>&gt; &gt; 6 NT 1 UREA N6 3 -0.83<BR>&gt; &gt; 7 H 1 UREA H7 3 0.415<BR>&gt; &gt; 8 H 1 UREA H8 3 0.415<BR>&gt; &gt; #else<BR>&gt; &gt; #endif[ moleculetype ]<BR>&gt; &gt; <BR>&gt; &gt; so now, my urea.itp file is:<BR>&gt; &gt; <BR>&gt; &gt; ; name nrexcl<BR>&gt; &gt; Urea 3<BR>&gt; &gt; [ atoms ]<BR>&gt; &gt; ; nr type resnr residu atom cgnr charge<BR>&gt; &gt; 1 C 1 UREA C1 1 0.683<BR>&gt; &gt; 2 O 1 UREA O2 1 -0.683<BR>&gt; &gt; 3 NT 1 UREA N3 2 -0.622<BR>&gt; &gt; 4 H 1 UREA H4 2 0.346<BR>&gt; &gt; 5 H 1 UREA H5 2 0.276<BR>&gt; &gt; 6 NT 1 UREA N6 3 -0.622<BR>&gt; &gt; 7 H 1 UREA H7 3 0.346<BR>&gt; &gt; 8 H 1 UREA H8 3 0.276<BR>&gt; &gt; [ bonds ]<BR>&gt; &gt; ; ai aj funct b0 kb<BR>&gt; &gt; 3 4 1 1.000000e-01 3.744680e+05<BR>&gt; &gt; 3 5 1 1.000000e-01 3.744680e+05<BR>&gt; &gt; 6 7 1 1.000000e-01 3.744680e+05<BR>&gt; &gt; 6 8 1 1.000000e-01 3.744680e+05<BR>&gt; &gt; 1 2 1 1.230000e-01 5.020800e+05<BR>&gt; &gt; 1 3 1 1.330000e-01 3.765600e+05<BR>&gt; &gt; 1 6 1 1.330000e-01 3.765600e+05<BR>&gt; &gt; [ pairs ]<BR>&gt; &gt; ; ai aj funct c6 c12<BR>&gt; &gt; 2 4 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; 2 5 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; 2 7 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; 2 8 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; 3 7 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; 3 8 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; 4 6 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; 5 6 1 0.000000e+00 0.000000e+00<BR>&gt; &gt; [ angles ]<BR>&gt; &gt; ; ai aj ak funct th0 cth<BR>&gt; &gt; 1 3 4 1 1.200000e+02 2.928800e+02<BR>&gt; &gt; 1 3 5 1 1.200000e+02 2.928800e+022 1 6 1 1.215000e+02 5.020800e+02<BR>&gt; &gt; 3 1 6 1 1.170000e+02 5.020800e+02<BR>&gt; &gt; [ dihedrals ]<BR>&gt; &gt; ; ai aj ak al funct phi cp mult<BR>&gt; &gt; 2 1 3 4 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; 6 1 3 4 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; 2 1 3 5 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; 6 1 3 5 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; 2 1 6 7 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; 3 1 6 7 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; 2 1 6 8 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; 3 1 6 8 1 1.800000e+02 3.347200e+01 2.000000e+00<BR>&gt; &gt; [ dihedrals ]<BR>&gt; &gt; ; ai aj ak al funct q0 cq<BR>&gt; &gt; 3 4 5 1 2 0.000000e+00 1.673600e+02<BR>&gt; &gt; 6 7 8 1 2 0.000000e+00 1.673600e+02<BR>&gt; &gt; 1 3 6 2 2 0.000000e+00 1.673600e+02<BR>&gt; &gt; <BR>&gt; &gt; 4 3 5 1 1.200000e+02 3.347200e+02<BR>&gt; &gt; 1 6 7 1 1.200000e+02 2.928800e+02<BR>&gt; &gt; 1 6 8 1 1.200000e+02 2.928800e+02<BR>&gt; &gt; 7 6 8 1 1.200000e+02 3.347200e+02<BR>&gt; &gt; 2 1 3 1 1.215000e+02 5.020800e+02<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; Following Chapter 5 in the manual (page 102), the file described for <BR>&gt; &gt; urea.itp is the same as mine. For the topology (page 103)<BR>&gt; &gt; #include "ffgmx.itp" is used, but I don´t see any more different.<BR>&gt; &gt; With this modification in the urea.itp I get the same error. Any idea of <BR>&gt; &gt; what could be the problem?<BR>&gt; &gt; Thank you very much for your help,<BR>&gt; &gt; <BR>&gt; &gt; Rebeca Garcia <BR>&gt; &gt; Parc Cientific de Barcelona<BR>&gt; &gt; regafan@hotmail.com &lt;mailto:regafan@hotmail.com&gt;<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; &gt; Date: Thu, 18 Dec 2008 09:21:41 -0500<BR>&gt; &gt; &gt; From: jalemkul@vt.edu<BR>&gt; &gt; &gt; To: gmx-users@gromacs.org<BR>&gt; &gt; &gt; Subject: Re: [gmx-users] Urea topology problem<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Rebeca García Fandiño wrote:<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; &lt;snip&gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; ERROR 1 [file solvated.top, line 39]:<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; atom C1 (Res UREA-1) has mass 0<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; In urea.itp, no masses are defined. If you correct the format of this <BR>&gt; &gt; file (see<BR>&gt; &gt; &gt; Chapter 5 of the manual), then this issue should be resolved.<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Which force field does this urea.itp correspond to? Where should it <BR>&gt; &gt; look<BR>&gt; &gt; &gt; &gt; for the atomtypes of urea?<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; The .atp file corresponding to the force field you are using <BR>&gt; &gt; (ffG43a2). The<BR>&gt; &gt; &gt; atomtypes for urea appear to be generic for use with the Gromos96 <BR>&gt; &gt; force fields.<BR>&gt; &gt; &gt; There are specific atomtypes within ffG53a6 for urea, if you want to <BR>&gt; &gt; use the<BR>&gt; &gt; &gt; newer force field (check the ffG53a6.rtp file for the urea parameters).<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; -Justin<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Thank you very much for your help,<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Rebeca Garcia<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Parc Cientific de Barcelona<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; regafan@hotmail.com<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; <BR>&gt; &gt; ------------------------------------------------------------------------<BR>&gt; &gt; &gt; &gt; ¿Aún no tienes Internet Explorer 7? Bájatelo y consigue un regalo <BR>&gt; &gt; gratis<BR>&gt; &gt; &gt; &gt; &lt;http://vivelive.com/ieak7/&gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; <BR>&gt; &gt; ------------------------------------------------------------------------<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; _______________________________________________<BR>&gt; &gt; &gt; &gt; gmx-users mailing list gmx-users@gromacs.org<BR>&gt; &gt; &gt; &gt; http://www.gromacs.org/mailman/listinfo/gmx-users<BR>&gt; &gt; &gt; &gt; Please search the archive at http://www.gromacs.org/search before <BR>&gt; &gt; posting!<BR>&gt; &gt; &gt; &gt; Please don't post (un)subscribe requests to the list. Use the<BR>&gt; &gt; &gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<BR>&gt; &gt; &gt; &gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; --<BR>&gt; &gt; &gt; ========================================<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Justin A. Lemkul<BR>&gt; &gt; &gt; Graduate Research Assistant<BR>&gt; &gt; &gt; Department of Biochemistry<BR>&gt; &gt; &gt; Virginia Tech<BR>&gt; &gt; &gt; Blacksburg, VA<BR>&gt; &gt; &gt; jalemkul[at]vt.edu | (540) 231-9080<BR>&gt; &gt; &gt; http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; ========================================<BR>&gt; &gt; &gt; _______________________________________________<BR>&gt; &gt; &gt; gmx-users mailing list gmx-users@gromacs.org<BR>&gt; &gt; &gt; http://www.gromacs.org/mailman/listinfo/gmx-users<BR>&gt; &gt; &gt; Please search the archive at http://www.gromacs.org/search before <BR>&gt; &gt; posting!<BR>&gt; &gt; &gt; Please don't post (un)subscribe requests to the list. Use the<BR>&gt; &gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<BR>&gt; &gt; &gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; ------------------------------------------------------------------------<BR>&gt; &gt; Descárgate gratis el nuevo Windows Live Messenger <BR>&gt; &gt; &lt;http://download.live.com/&gt;<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; ------------------------------------------------------------------------<BR>&gt; &gt; <BR>&gt; &gt; _______________________________________________<BR>&gt; &gt; gmx-users mailing list gmx-users@gromacs.org<BR>&gt; &gt; http://www.gromacs.org/mailman/listinfo/gmx-users<BR>&gt; &gt; Please search the archive at http://www.gromacs.org/search before posting!<BR>&gt; &gt; Please don't post (un)subscribe requests to the list. Use the <BR>&gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<BR>&gt; &gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>&gt; <BR>&gt; <BR>&gt; -- <BR>&gt; David van der Spoel, Ph.D., Professor of Biology<BR>&gt; Molec. Biophys. group, Dept. of Cell &amp; Molec. Biol., Uppsala University.<BR>&gt; Box 596, 75124 Uppsala, Sweden. Phone: +46184714205. Fax: +4618511755.<BR>&gt; spoel@xray.bmc.uu.se spoel@gromacs.org http://folding.bmc.uu.se<BR>&gt; _______________________________________________<BR>&gt; gmx-users mailing list gmx-users@gromacs.org<BR>&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<BR>&gt; Please search the archive at http://www.gromacs.org/search before posting!<BR>&gt; Please don't post (un)subscribe requests to the list. Use the <BR>&gt; www interface or send it to gmx-users-request@gromacs.org.<BR>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR><BR><br /><hr />Comparte hasta 500 fotos en un solo email con  <a href='http://download.live.com/' target='_new'>Windows Live</a></body>
</html>