<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-9"
      http-equiv="Content-Type">
  </head>
  <body text="#000000" bgcolor="#ffffff">
    On 5/04/2011 3:57 PM, ahmet yżldżrżm wrote:
    <blockquote
      cite="mid:BANLkTik9Z8NaDUuW9L5yH2x2BsQR4WRO0w@mail.gmail.com"
      type="cite">Dear Tsjerk,<br>
      <br>
      <blockquote style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px
        solid rgb(204, 204, 204); padding-left: 1ex;"
        class="gmail_quote">Hi Ahmet,<br>
        <br>
        As suggested, it's better to break up your molecule into smaller<br>
        charge groups. Note that charge groups don't need to have zero
        charge,<br>
        nor integer charge. In your case, I'd suggest two COH groups for
        EDO,<br>
        which will have zero net charge each, and for TRS I'd take the
        COH<br>
        groups as separate charge groups. I also note that the COH
        groups,<br>
        although chemically identical - H3NC(COH)3, right?-, have
        different<br>
        charges. That doesn't seem proper.<br>
        <br>
        Hope it helps,<br>
        <br>
        Tsjerk<br>
      </blockquote>
      <br>
      nonrevised .itp file:<br>
      EDO      3<br>
      <br>
      [ atoms ]<br>
      ;   nr      type  resnr resid  atom  cgnr   charge     mass<br>
           1        OA     1  EDO     OAB     1   -0.111  15.9994   <br>
           2         H     1  EDO     HAE     1    0.031   1.0080   <br>
           3       CH2     1  EDO     CAA     1    0.080  14.0270   <br>
           4       CH2     1  EDO     CAC     1    0.080  14.0270   <br>
           5        OA     1  EDO     OAD     1   -0.111  15.9994   <br>
           6         H     1  EDO     HAF     1    0.031   1.0080  <br>
      <br>
      nonrevised .itp file:<br>
      EDO      3<br>
      <br>
      [ atoms ]<br>
      ;   nr      type  resnr resid  atom  cgnr   charge     mass<br>
           1        OA     1  EDO     OAB     1   -0.111  15.9994   <br>
           2         H     1  EDO     HAE     1    0.031   1.0080   <br>
           3       CH2     1  EDO     CAA     1    <b>0.000</b> 
      14.0270   <br>
           4       CH2     1  EDO     CAC     1    <b>0.000</b> 
      14.0270   <br>
           5        OA     1  EDO     OAD     1   -0.111  15.9994   <br>
           6         H     1  EDO     HAF     1    0.031   1.0080  <br>
      <br>
      can you show me on the itp file? how do I seperate two COH groups?
      Please help me<br>
    </blockquote>
    <br>
    We didn't say to use zero *charges*. We said you need small charge
    groups. Since EDO is neutral, then by symmetry each half must also
    be neutral.<br>
    <br>
    Have a look at the way charge groups are used in the .rtp file for
    your forcefield. You need to do similar things. If you don't
    understand the file formats yet, you need to read up on chapter 5 of
    the manual.<br>
    <br>
    Mark<br>
  </body>
</html>