Dear David,<br>
<br>
I followed your instructions and calculated Heat of vaporization of my 
alkane once with one molecule in gas phase (no cutoff) and once with 
equivalent number of molecules as in liquid phase as Justin suggested. 
Results are as follows:<br>
<br>
<b>one single molecule in gas phase</b><br>
<br>
Energy                      Average   Err.Est.       RMSD  Tot-Drift<br>
-------------------------------------------------------------------------------<br>
LJ (SR)                    -2.24473      0.073      1.292   0.342696  (kJ/mol)<br>



Coulomb (SR)                11.5723       0.55    2.17577   -2.33224  (kJ/mol)<br>
Potential                  <b>  59.244 </b>      0.94    10.9756    6.35631  (kJ/mol)<br>
Total Energy                106.647          1    15.4828    6.78792  (kJ/mol)<br>
<br>
<b>equivalent number of molecules as in liquid</b> ( large box 20 nm)<br>
<br>
Statistics over 1000001 steps [ 0.0000 through 2000.0000 ps ], 4 data sets<br>



All statistics are over 100001 points<br>
<br>
Energy                      Average   Err.Est.       RMSD  Tot-Drift<br>
-------------------------------------------------------------------------------<br>
LJ (SR)                    -2.16367      0.053   0.171542   0.374027  (kJ/mol)<br>



Coulomb (SR)                11.2894       0.23    0.49105   -1.44437  (kJ/mol)<br>
Potential                 <b>  63.2369    </b>    1.1    2.47211    7.69756  (kJ/mol)<br>
Total Energy                114.337        1.1    2.65547    7.72258  (kJ/mol)<br>
<br>
<br>
 Since pbc is set to NO molecules leave the box and I dont know if this all right. I hope the difference is acceptable...!<br>
<br>
0- I am going to do the same calculation but for some polymers solvated 
in the alkane. For binary system do I need to look at nonboded terms? 
and then run a simulation for a single polymer in vacuum? <br>
<br>
Can you please provide me with a recipe for Delta Hvap of the solute in a solvent?<br>
<br>
<br>
1- If I want to look at nonboded interactions only, do I have to add  Coul. recip.  to [ LJ (SR)  + Coulomb (SR) ] ?<br>
<br>
<br>
coulombtype         =  PME                      <br>
vdw-type            =  Shift               <br>
rcoulomb-switch     =  0                   <br>
rvdw-switch         =  0.9 <br>
rlist               =  1.2                <br>
rcoulomb            =  1.2             <br>
rvdw                =  1.0        <br>
<br>
<br>
ff: OPLSAA<br>
<br>
<br>
Appreciate your help,<br>
Best,<br>
<br>