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