<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">

<html xmlns="http://www.w3.org/1999/xhtml">
  <head>
    <meta content="text/html; charset=UTF-8" http-equiv="Content-Type" />
    <title></title>
  </head>

  <body>
    <p style="margin: 0px;"><span>&#160;</span></p>

    <div style="margin: 5px 0px 5px 0px;">
      <br />
       On February 14, 2011 at 3:28 PM Moeed &lt;lecielll@googlemail.com&gt; wrote:<br />
      <br />

      <blockquote type="cite" style="margin-left: 0px; padding-left: 10px; border-left: solid 1px blue;">
        <br />
         Dear experts,<br />
        <br />
         I am going to do solvation FE of polymer (polyethylene) in a hydrocarbon solvent. I have prepared a system consisting of 4 polymer chains and 480 hexane molecules with the actual density of polymer solution (~ 0.5 g/cm3).<br />
        <br />
         1- For such a study I dont know how many polymers I need to have in my system. If FE can be done with only one chain, am I making system bigger in vain? Does this matter affect the accuracy of results?<br />
      </blockquote>
    </div>

    <p style="margin: 0px;">How close are the polymers?</p>

    <p style="margin: 0px;">&#160;</p>

    <div style="margin: 5px 0px 5px 0px;">
      <blockquote type="cite" style="margin-left: 0px; padding-left: 10px; border-left: solid 1px blue;">
        2- I have switched off electrostatics so I am using<br />
        <br />
         free_energy&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; yes<br />
         init_lambda&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; 0&#160;<br />
         delta_lambda&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; 0<br />
         sc_alpha&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; 0.5<br />
         sc-power&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; 1<br />
         sc_sigma&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; 0.3<br />
         couple-lambda0&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; vdw<br />
         couple-lambda1&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160; none<br />
         couple-intramol&#160;&#160;&#160;&#160;&#160; =&#160;&#160; no<br />
        <br />
         In David Mobley&#39;s turorial the last three lines are not included. I wanted to know if I am to run say 10 simulations for different lambda, what purpose does the last three lines serve in 4.0.7&#160; ? I got very close values in that tutorial without these settings. ( I know what these lines mean, just curious how these three lines affect the results in 4 X +).<br />
      </blockquote>
    </div>

    <p style="margin: 0px;">&#160;</p>

    <p style="margin: 0px;">Have you thought of using Bennett&#39;s Acceptance Ratio? This make for a little more work but knowledge of your error is nice.</p>

    <p style="margin: 0px;">&#160;</p>

    <div style="margin: 5px 0px 5px 0px;">
      <blockquote type="cite" style="margin-left: 0px; padding-left: 10px; border-left: solid 1px blue;">
        Please let me know your comments/point of view about the system and setting I am using.<br />
        <br />
         Thanks<br />
         Moeed<br />
      </blockquote>
    </div>

    <p style="margin: 0px;">&#160;</p>

    <p style="font-family: monospace; white-space: nowrap; margin: 5px 0px 5px 0px;">TJ Mustard<br />
    Email: mustardt@onid.orst.edu</p>
  </body>
</html>