<P>
Dear GROMACS users,<BR>
<BR>
A method to locate glass transtion temperature of liquids is reported by Stanley et al, PRL 2004, 93, pp 047801: system is cooled initially and then is heated , in the range 0K-400K, in a "modified" NVT ensemble. (Total energies of the system with respect to temp. is plotted, a kink in the plot is refered to the glass transition)<BR>
<BR>
In a standard NVT simulation, To is constant. In this method, the Berendsen thermostat is used but, at everytime step, the bath temperature To is updated such that <BR>
To(t)= To(t=0) + deltaT, with deltaT= cooling/heating * timestep.<BR>
This is done at every time step before rescaling the velocities.<BR>
<BR>
Now, I would like to know how can I 'impliment' this in a GROMACS version; I mean where should I modify the code? Or is there any other way to achieve this without modifying the code?<BR>
<BR>
Is it /usr/local/gromacs/src/mdlib/tgroup.c ?<BR>
<BR>
Thanks for your help,<BR>
Jes<BR>
<BR>
</P>
<br><br>
<Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style='font-family:Verdana;font-size:11px;line-height:15px;'><TR><td><a href='http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-default.htm/1050715198@Middle5/2212064_2204411/2207095/1?PARTNER=3&OAS_QUERY=null' target=new ><img src ='http://imadworks.rediff.com/cgi-bin/AdWorks/adimage.cgi/2212064_2204411/creative_2207095.gif' alt='Sharekhan' border=0></a></td></TR></Table>