<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><div class="h5">
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Thanks David for replying me. I appreciate your attention to comment on the following short inquiries.<br>
<br>
 1- g_energy says:<br>
<br>
you will need to give the number of constraints per molecule –nconstr and</blockquote></div></div></blockquote><div>for water this is 3.  Can you elaborate on how to set this option. When constraints= all-bonds how to set -nconstr for a simlpe alkane for instance? It refers to number of atoms of each molecule on which constraints are applied?<br>
</div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><div class="h5"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">

<br> 
<br>
2- Aside from heat capacities, does g_energy report the following properties correctly?<br>
<br>
Thermal expansion coeff. (NPT): Enthalpy, Vol, Temp<br>
Isothermal compressibility:     Vol, Temp<br>
Adiabatic bulk modulus:         Vol, Temp<br></blockquote></div></div></blockquote><div><br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div><div class="h5"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">3- Since volume appears for calculation of the above properties , does this imply ONLY NPT results in these properties? (because in NPT g_energy reports &quot;volume&quot;).<br>

<br>
4- For thermal expansion: I read one of your earlier messages saying that enthalpy = Etot(total energy) +pv,  [where Etot is potential+kinetic (in g_energy list)]<br></blockquote></div></div></blockquote><div>            and for pv , take ref_p and not measure p (that g_energy reports). If the system is equilibrated to ref_p, then why not taking pv from g_energy? you mean product of pv that        g_energy gives is not correct?<br>
<br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><div class="h5"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">

<br>
Many thanks for you time!
<br>
Best,<br><br></blockquote><font color="#888888"></font></div></div></blockquote></div><br>