Hello,<br><br> I'm working with a polarizable model of Water (TTM2.1F) and I try to transcribe it to Gromacs' "grammar" using a dummy shell particle to describe polarizability. Both the model with shell particle and without it worked properly in the NVT ensemble. In the NVE ensemble, nevertheless, the model without shell particle worked well conserving energy but the model with the shell particle didn't. (with the same set of parameters specially those of Ewald and Shake tolerance).<br> So is it possible that the iterative process to arrive a minimum of the electric field in the shell polarization model really maps the whole system from t to t+dt in a non conservative way?<br> Could I ask this question to the programmer? (did He/She implement the model only for the NVT ensemble?)<br> Do someone know some paper of a NVE Ensemble employing shell polarizability of any chemical nature --not only water--?<br><br> I would be very thankful with the
information<br><br> Juan S. Medina Álvarez, Ph. D.<br>IUSIANI, ULPGC, SPAIN<br>tlazcala@yahoo.com<br><p> 
<hr size=1><br><font face="Verdana" size="-2">Enviado desde <a href="http://us.rd.yahoo.com/mailuk/taglines/isp/control/*http://us.rd.yahoo.com/evt=52431/*http://es.docs.yahoo.com/mail/overview/index.html">Correo Yahoo!</a><br>La bandeja de entrada más inteligente.<br></font>