For vc-md, only ion_temperature = 'rescaling' is implemented... ,<br> and with  tempw =  use tolp= and nraise= to control temprature tolerance ,  <br><br><div class="gmail_quote">On Fri, Oct 12, 2012 at 4:26 AM, Riping Wang <span dir="ltr"><<a href="mailto:wang.riping.81@gmail.com" target="_blank">wang.riping.81@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Forum,<br><br>It try to heat vc-md calculations.<br>But I can not reach the target temperature.<br>When I set tempw = 100 ,  it will stable at about 38K<br>
When I set tempw=1000, it will stable at about 400K.<br>I pasted part of my input bellow.<br>
Any one have some idea?<br><br>Thank you very much.<br><br>'''<br>&control<br>    calculation = 'vc-md' ,<br>    prefix = 'alpo4-10' ,<br>    restart_mode = 'from_scratch' ,<br>    nstep = 30000 , <br>

    /iprint = 1 , <br>    dt = 20 , <br>/<br>&ions <br>    pot_extrapolation='second-order'<br>    wfc_extrapolation='second-order' <br>    ion_temperature = 'berendsen' ,<br>    tempw = 100 ,  <br>

/<br>&cell <br>    cell_dynamics = 'pr' ,  <br>    press = 0 , <br>'''<span class="HOEnZb"><font color="#888888"><br clear="all"><br>WANG Riping<br>2012.10.12<br><br>-- <br><div>******************************************************************************<br>

WANG Riping<br>Ph.D student, </div>
<div>Institute for Study of the Earth's Interior,Okayama University,<br>827 Yamada, Misasa, Tottori-ken 682-0193, Japan<br>Tel: <a href="tel:%2B81-858-43-3739" value="+81858433739" target="_blank">+81-858-43-3739</a>(Office), <a href="tel:%2B81-858-43-1215" value="+81858431215" target="_blank">+81-858-43-1215</a>(Inst)</div>

<div>E-mail: <a href="mailto:wang.riping.81@gmail.com" target="_blank">wang.riping.81@gmail.com</a><br>******************************************************************************<br><br><br></div><br>
</font></span><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br>