<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    On 30/03/2012 12:40 AM, rama david wrote:
    <blockquote
cite="mid:CAD=-SYEN0DuaPoPDCev1e1d7Ltckpp7ES_V7vaYE9a31uezmMw@mail.gmail.com"
      type="cite">Hi mark ..
      <div>thank you for your&nbsp;suggestion..</div>
      <div>My command line is g_mindist &nbsp;-f &nbsp;..xtc -s ..tpr -od
        minimal.xvg &nbsp;-pi&nbsp;</div>
      <div>&nbsp; protein group contain information about the 284 atom (i.e.
        All protein atoms)</div>
    </blockquote>
    <br>
    So clearly some part of some molecule is protruding across periodic
    boundaries, since the nearest distance between any periodic image of
    any atom to any original atom is so small. Since your initial
    structure isn't in the square arrangement you think it is in, I
    suggest you go back and set up your system better.<br>
    <br>
    Mark<br>
    <br>
    <blockquote
cite="mid:CAD=-SYEN0DuaPoPDCev1e1d7Ltckpp7ES_V7vaYE9a31uezmMw@mail.gmail.com"
      type="cite">
      <div>..</div>
      <div>So please suggest &nbsp;me the right way ..</div>
      <div><br>
        <div class="gmail_quote">On Thu, Mar 29, 2012 at 2:02 PM, 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>
          wrote:<br>
          <blockquote class="gmail_quote" style="margin:0 0 0
            .8ex;border-left:1px #ccc solid;padding-left:1ex">
            <div class="im">On 29/03/2012 7:20 PM, rama david wrote:<br>
              <blockquote class="gmail_quote" style="margin:0 0 0
                .8ex;border-left:1px #ccc solid;padding-left:1ex">
                Hi everybody ,<br>
                &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; I run simulation of 4 same &nbsp;molecule keep
                apart in box<br>
                of 4 4 4 dimension ..( 71 atom in one molecule = 71 * 4
                = total atom are 284 )<br>
                force field = gromacs96 &nbsp; 53a6<br>
                <br>
                COM (center of mass) infirmation of molecules<br>
                system size : &nbsp;1.255 &nbsp;1.577 &nbsp;1.883<br>
                box vectors : &nbsp;4.000 &nbsp; 4.000 &nbsp; 4.000 (nm)<br>
                &nbsp;mol1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: 2.057 &nbsp; &nbsp;0.844 &nbsp;1.941<br>
                &nbsp;mol 2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp;2.057 &nbsp;0.844 &nbsp;3.141<br>
                &nbsp;mol 3 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: 2.057 &nbsp; 3.244 &nbsp; 0.744<br>
                &nbsp;mol 4 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: 2.057 &nbsp; 3.244 &nbsp; 3.141<br>
                <br>
                (Four molecule are kept at the four corner of square<br>
                &nbsp;of each side 2.4 nm<br>
                &nbsp; &nbsp;four molecule are catenated in same pdb &nbsp; &nbsp;)<br>
              </blockquote>
              <br>
            </div>
            No they're not starting on a square. Look at those z
            coordinates above.
            <div class="im"><br>
              <br>
              <blockquote class="gmail_quote" style="margin:0 0 0
                .8ex;border-left:1px #ccc solid;padding-left:1ex">
                <br>
                &nbsp;my md.mdp input is like the ..<br>
                <br>
                ;Neighborsearching<br>
                ns_type &nbsp; &nbsp; &nbsp; &nbsp;= grid &nbsp; &nbsp; &nbsp; &nbsp;; search neighboring grid
                cells<br>
                nstlist &nbsp; &nbsp; &nbsp; &nbsp;= 5 &nbsp; &nbsp; &nbsp; &nbsp;; 10 fs<br>
                rlist &nbsp; &nbsp; &nbsp; &nbsp;= 1.0 &nbsp; &nbsp; &nbsp; &nbsp;; short-range neighborlist
                cutoff (in nm)<br>
                rcoulomb &nbsp; &nbsp;= 1.0 &nbsp; &nbsp; &nbsp; &nbsp;; short-range electrostatic
                cutoff (in nm)<br>
                rvdw &nbsp; &nbsp; &nbsp; &nbsp;= 1.0 &nbsp; &nbsp; &nbsp; &nbsp;; short-range van der Waals
                cutoff (in nm)<br>
                ; Electrostatics<br>
                coulombtype &nbsp; &nbsp;= PME &nbsp; &nbsp; &nbsp; &nbsp;; Particle Mesh Ewald for
                long-range electrostatics<br>
                pme_order &nbsp; &nbsp;= 4 &nbsp; &nbsp; &nbsp; &nbsp;; cubic interpolation<br>
                fourierspacing &nbsp; &nbsp;= 0.16 &nbsp; &nbsp; &nbsp; &nbsp;; grid spacing for FFT<br>
                <br>
                With &nbsp;command &nbsp;g_mindist , select option - 1 (Protein )<br>
                I got the following result..<br>
                <br>
                The shortest periodic distance is 0.141718 (nm) at time
                7692 (ps),<br>
                between atoms 26 and 111<br>
                <br>
                Is the simulation is behaving abnormal(I.e &nbsp;simulation
                is wrong ) or &nbsp;I<br>
                &nbsp;have to select the option system on prompting ??? I am
                very new to these simulation field..<br>
                so all suggestion are appreciable ...<br>
              </blockquote>
              <br>
            </div>
            Shrug. You've measured the inter-group distance between a
            group and itself. Curiously enough, that's the length of a
            C-C bond, or similar. g_mindist -h is required reading.
            Also, we don't even know what's in group Protein, or your
            g_mindist command line.<span class="HOEnZb"><font
                color="#888888"><br>
                <br>
                Mark<br>
                -- <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></span></blockquote>
        </div>
        <br>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
    </blockquote>
    <br>
  </body>
</html>