Hello to all <br> I have tried to reproduce the hydration free
energy (TI) of the ethanol from Hess and van der Vegt (JPCB, 110, 17616). <br>The value I have
obtained is around 20kJ/mol while the reference value is -20.1kJ/mol
(if not the sign ...). <br> If someone can help me find the mistake I would be very grateful. <br> Below are the simulation details. <br> <br> I followed the protocol of the paper, Berk: <br> I used 47 lambda values (because of hydrogen bonds between solute-solvent) (dense near lambda=0 and between 0.46 and 0.72). <br>
I am turning off the LJ and Coulomb terms separately.<br>Softcore (alpha = 0.5, power = 1), OPLS-AA<br> Timestep = 2fs, <br>sd = integrator, <br>PME<br>constrained = none<br>trajectories = 40ps (NVT) (I know that is small, but I looking for qualitative results) after 20ps of equilibration (NPT).<br>
<br>I have found a value of -4.8kJ/mol for the DeltaG(vacuum) (relative to mutate the ethanol to dummy in vacuum).<br>My dv/dl curve (for DeltaG(water) is below. Should I expect this form?<br>I think there are many &quot;high&quot; positive values. The numerical integration is 15.9kJ.<br>
So, DG(hyd) = DG(wat) - DG(vac) = 15.9 - (-4.8) = 20.7kJ<br><br>I believe that this protocol is OK (but I want to confirm that). <br><br>@&nbsp;&nbsp;&nbsp; title &quot;dG/d\8l\4&quot;<br>@&nbsp;&nbsp;&nbsp; xaxis&nbsp; label &quot;Time (ps)&quot;<br>@&nbsp;&nbsp;&nbsp; yaxis&nbsp; label &quot;dG/d\8l\4 (kJ mol\S-1\N [\8l\4]\S-1\N)&quot;<br>
@TYPE xy<br>0.000&nbsp;&nbsp; 4.431073e+02<br>0.005&nbsp;&nbsp; 3.660083e+02<br>0.010&nbsp;&nbsp; 2.987275e+02<br>0.015&nbsp;&nbsp; 2.023206e+02<br>0.000&nbsp;&nbsp; 4.431073e+02<br>0.005&nbsp;&nbsp; 3.660083e+02<br>0.010&nbsp;&nbsp; 2.987275e+02<br>0.015&nbsp;&nbsp; 2.023206e+02<br>0.020&nbsp;&nbsp; 1.905198e+02<br>
0.030&nbsp;&nbsp; 1.086573e+02<br>0.040&nbsp;&nbsp; 4.217960e+01<br>0.050&nbsp;&nbsp; 4.101856e+01<br>0.060&nbsp;&nbsp; 2.314358e+01<br>0.070&nbsp;&nbsp; 2.223774e+01<br>0.080&nbsp;&nbsp; 2.669019e+01<br>0.090&nbsp;&nbsp; 1.696131e+01<br>0.100&nbsp;&nbsp; 6.089735e+00<br>0.110&nbsp;&nbsp; 9.569030e+00<br>0.120&nbsp;&nbsp; 1.505562e+01<br>
0.130&nbsp;&nbsp; 2.785974e+00<br>0.140&nbsp;&nbsp; 2.440906e+00<br>0.150&nbsp;&nbsp; 3.050463e+00<br>0.160&nbsp;&nbsp; 8.578429e-01<br>0.200&nbsp;&nbsp; 3.577259e+00<br>0.240&nbsp; -3.532969e+00<br>0.280&nbsp;&nbsp; 5.842623e+00<br>0.320&nbsp;&nbsp; 7.912565e+00<br>0.360&nbsp;&nbsp; 9.322502e+00<br>0.400&nbsp; -2.930754e+00<br>
0.440&nbsp;&nbsp; 2.797944e+00<br>0.460&nbsp;&nbsp; 5.113900e+00<br>0.480&nbsp; -3.919665e+00<br>0.500&nbsp; -9.000282e+00<br>0.520&nbsp; -7.636168e+00<br>0.540&nbsp; -7.574299e+00<br>0.560&nbsp; -2.144124e+01<br>0.580&nbsp; -6.618514e+00<br>0.600&nbsp; -2.755925e+01<br>0.620&nbsp; -1.491331e+01<br>
0.640&nbsp; -2.541398e+01<br>0.660&nbsp; -5.213154e+01<br>0.680&nbsp; -2.364612e+01<br>0.700&nbsp; -2.535817e+01<br>0.720&nbsp; -7.767502e+00<br>0.760&nbsp; -1.096362e+01<br>0.800&nbsp;&nbsp; 1.136031e+01<br>0.840&nbsp;&nbsp; 2.883933e+01<br>0.880&nbsp;&nbsp; 4.465743e+01<br>0.920&nbsp;&nbsp; 5.461566e+01<br>
0.960&nbsp;&nbsp; 6.882000e+01<br>1.000&nbsp;&nbsp; 7.673736e+01<br clear="all"><br>-- <br>_______________________________________<br>Eudes Eterno Fileti<br>Centro de Ciência Naturais e Humanas<br>Universidade Federal do ABC<br>Rua Catequese, 242 - 3º Andar<br>
09090-400 Santo André - SP Brasil<br>Tel: +55 11 4437-1600 ramal 408<br>skype: eefileti<br><a href="http://cromo.ufabc.edu.br/~fileti/">http://cromo.ufabc.edu.br/~fileti/</a>