<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body text="#000000" bgcolor="#ffffff">
    On 19/02/2011 3:02 PM, C.Y. Chang wrote:
    <blockquote
      cite="mid:AANLkTimB3XMPd9wcg1V+jwKZsSk88U4yP+_0EzArU5je@mail.gmail.com"
      type="cite">Hi,<br>
      <br>
      I have tried to add<br>
      <br>
      [ position restraints ]<br>
    </blockquote>
    <br>
    This is misspelled. Surely grompp warned about this?<br>
    <br>
    <blockquote
      cite="mid:AANLkTimB3XMPd9wcg1V+jwKZsSk88U4yP+_0EzArU5je@mail.gmail.com"
      type="cite">2 1 1000 1000 1000 ; Restrain to a point<br>
      1 1 1000 0 1000 ; Restrain to a line (Y-axis)<br>
      3 1 1000 0 0 ; Restrain to a plane (Y-Z-plane)<br>
      <br>
      in the end of the topology file.<br>
    </blockquote>
    <br>
    This does not even approach the solution to your problem. Position
    restraints inhibit diffusion and structural changes. You need a
    vsite, like you knew earlier.<br>
    <br>
    <blockquote
      cite="mid:AANLkTimB3XMPd9wcg1V+jwKZsSk88U4yP+_0EzArU5je@mail.gmail.com"
      type="cite">In the em. process, I get the eroor msg.<br>
    </blockquote>
    <br>
    No, this error happened in grompp.<br>
    <br>
    <blockquote
      cite="mid:AANLkTimB3XMPd9wcg1V+jwKZsSk88U4yP+_0EzArU5je@mail.gmail.com"
      type="cite">Fatal error:<br>
      Invalid dihedral type 1000<br>
    </blockquote>
    <br>
    GROMACS didn't recognise the mis-spelled directive, and so it's
    trying to make sense of your position restraint lines as dihedrals.<br>
    <br>
    Mark<br>
    <br>
    <blockquote
      cite="mid:AANLkTimB3XMPd9wcg1V+jwKZsSk88U4yP+_0EzArU5je@mail.gmail.com"
      type="cite">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Chia-yun<br>
      <br>
      <br>
      <div class="gmail_quote">2011/2/18 Mark Abraham <span dir="ltr">&lt;<a
            moz-do-not-send="true" href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt;</span><br>
        <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
          0.8ex; border-left: 1px solid rgb(204, 204, 204);
          padding-left: 1ex;">
          <div class="im">On 18/02/2011 1:13 PM, C.Y. Chang wrote:<br>
            <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
              0.8ex; border-left: 1px solid rgb(204, 204, 204);
              padding-left: 1ex;">
              Hi,<br>
              <br>
              I am dealing with the lipid bilayer permeation simulation.<br>
              Most compounds can be finished, but the compounds with CN
              can't be performed simulation.<br>
              I have searched the discussion in the gmx-users
              discussion.<br>
              The "vsite" has been mentioned, and I have refered to the
              gromacs manual.<br>
              But I don't understand that how I can use the "vsite".<br>
            </blockquote>
            <br>
          </div>
          The theory is discussed in chapter 4 and there's a brief
          example in 5.2.2. What have you tried and what went wrong?<br>
          <br>
          Mark
          <div>
            <div class="h5"><br>
              <br>
              <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
                0.8ex; border-left: 1px solid rgb(204, 204, 204);
                padding-left: 1ex;">
                (add the toplogy file or use the command line?)<br>
                I attach the PDB and toplogy file of the molecule.<br>
                First, I performed the molecular dynamic simulation for
                the pure lipid bilayer, and the step was been finished.<br>
                After I insert the molecule, these command lines are
                performed.<br>
                <br>
                grompp -f minim.mdp -c cmplx.pdb -p topol_dppc.top -o
                em.tpr<br>
                mdrun -v -deffnm em<br>
                grompp -f npt_cmplx.mdp -c em.gro -p topol_dppc.top -o
                npt_cmplx.tpr<br>
                nohup mdrun -v -deffnm npt_cmplx &amp;<br>
                <br>
                Thanks for your help.<br>
                Best,<br>
                <br>
                &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
                &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;Chia-yun<br>
              </blockquote>
              <br>
            </div>
          </div>
          <font color="#888888">
            -- <br>
            gmx-users mailing list &nbsp; &nbsp;<a moz-do-not-send="true"
              href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
            <a moz-do-not-send="true"
              href="http://lists.gromacs.org/mailman/listinfo/gmx-users"
              target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
            Please search the archive at <a moz-do-not-send="true"
              href="http://www.gromacs.org/Support/Mailing_Lists/Search"
              target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a>
            before posting!<br>
            Please don't post (un)subscribe requests to the list. Use
            the www interface or send it to <a moz-do-not-send="true"
              href="mailto:gmx-users-request@gromacs.org"
              target="_blank">gmx-users-request@gromacs.org</a>.<br>
            Can't post? Read <a moz-do-not-send="true"
              href="http://www.gromacs.org/Support/Mailing_Lists"
              target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
          </font></blockquote>
      </div>
      <br>
    </blockquote>
    <br>
  </body>
</html>