<div dir="ltr">You will need to post your input and output files. Do you do GGA or hybrid calculations, how do you account for the vdW forces?<br><br>On Friday, December 5, 2014 4:17:45 PM UTC+1, bharat wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hello,<div>This is friendly reminder. Any suggestions?</div><div><br></div><div>Bharat<br><br>On Wednesday, December 3, 2014 12:56:10 PM UTC-5, bharat wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Experts,<div style="line-height:21.299999237060547px"><br style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px"><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">I am trying to reproduce Enthalpy of vaporization using PM3.</div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">Here are my calculation: </div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">(I took Total energy value from CP2K output file as an Internal energy, average is calculated over the MD. Experimental density is used for the constant volume for both liquid and vapor calculation). Am I taking correct energy value for internal energy? </div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px"><br></div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">U_vapor(Avg) = -20797.3 eV </div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">U_liquid(Avg) = -20802.9 eV</div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px"><br></div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">Delta_U = 5.6 eV <br></div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">I divided with 64 (# of water molecules) and converted to kcal/mol</div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">= 2.02 kcal/mol</div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px"><br></div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">I got the half value (reported value is 4.00 kcal/mol in G. Murdachaew et al. J. Phys. Chem. A 115, 6046 (2011)). Because of this value, I got only ~2.60 kcal/mol (RT= ~0.60 kcal/mol) of enthalpy of vaporization (Delta_H = delta_U + RT). </div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px"><br></div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">Can anyone please tell me where I am doing wrong? Where is the factor 2 missing?</div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px"><br></div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">Thank you.</div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px"><br></div><div style="color:rgb(68,68,68);font-family:Calibri,sans-serif;font-size:15px">Bharat </div><div><br></div></div></div></blockquote></div></div></blockquote></div>