Dear Vega,<br><br>The CP program prints different energy in the output file *.evp. In this file the eighth column is the constant of motion which should remain a constant. Plot this eighth column vs the timesteps and see what fluctuations are you getting in the plot. <br>
<br>This is I guess what Stefano asked you. <br><br>Regards,<br><br>Somesh<br><br><div class="gmail_quote">On Fri, Nov 27, 2009 at 12:42 PM, vega lew <span dir="ltr"><<a href="mailto:quantumdft@gmail.com">quantumdft@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear sir,<br><br>Did you ask whether the (etot+ekinc) conserve during the dynamic process?<br>
<br>I sum the corresponding columns. I found it is not changing very much.<br><br>For example, the (ekinc+ etot) is -4369.40879 (-4369.44669+0.0379) at 2200 step.<br>
Then it changes to -4369.41514 (-4369.45479+0.03965) at 2300 step.<br>And it sometimes goes a little up and sometimes a little down. But the value is in range of<br>-4369.41 and -4369.45.<br><br>thank you for reading.<br>
<font color="#888888">
<br>vega</font><div><div></div><div class="h5"><br><br><br><div class="gmail_quote">On Fri, Nov 27, 2009 at 2:02 AM, Stefano de Gironcoli <span dir="ltr"><<a href="mailto:degironc@sissa.it" target="_blank">degironc@sissa.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
is the total (kinetic + potential) energy conserved ?<br>
stefano<br>
<div><div></div><div><br>
vega lew wrote:<br>
> dear all,<br>
><br>
> I want to calculate certain heterogeneous reaction on a solid surface.<br>
> First input the initial structure from previous pw.x relaxations.<br>
> Second I fixed ions and let electron damping calculation run;<br>
> Then un-fixed the ions and let both electron and ions damping<br>
> calculations run.<br>
> Third I gave a random motion to the atoms in the system with<br>
> amprp(1)~(3) = 0.01<br>
> and then repeated second step to bring the system to the ground state.<br>
> Fourth I added temperature (300K) to the system, using verlet schedule.<br>
><br>
> But the temperature continuously going down, from 1000K to 150K during<br>
> the fisrt 3000 steps.<br>
> I think there might be something wrong with my cp.x calcuation.<br>
> Do you think the amprp(1)~(3) = 0.01 of random motion is too small?<br>
><br>
> Here is the input of my dynamics process,<br>
> &CONTROL<br>
><br>
> title = 'CP'<br>
> ,<br>
> calculation = 'cp'<br>
> ,<br>
> restart_mode = 'restart'<br>
> ,<br>
> outdir = '/disk5/nanjing/vega/cp/300/' ,<br>
> wfcdir = '/tmp/'<br>
> ,<br>
> pseudo_dir =<br>
> '/home/software/nanjing/vega/espresso-4.0.1/pseudo/' ,<br>
> prefix = 'cp_molecular' ,<br>
> nstep = 99999 ,<br>
> dt = 10,<br>
> iprint = 5,<br>
> isave =<br>
> 150,<br>
> /<br>
><br>
> &SYSTEM<br>
><br>
> ibrav =<br>
> 8,<br>
> celldm(1) =24.8624,<br>
> celldm(2) = 0.8520,<br>
> celldm(3) = 1.6964,<br>
> nat =<br>
> 147,<br>
> ntyp = 3,<br>
> nr1b = 20,<br>
> nr2b = 20,<br>
> nr3b = 20,<br>
> nbnd =<br>
> 630,<br>
> ecutwfc = 30,<br>
> ecutrho = 300,<br>
> /<br>
><br>
><br>
> &ELECTRONS<br>
> electron_dynamics = 'verlet',<br>
> /<br>
><br>
> &IONS<br>
><br>
> ion_dynamics = 'verlet'<br>
> ion_temperature = 'nose',<br>
> tempw = 300,<br>
> /<br>
> ..................<br>
> ..................<br>
><br>
> best wishes,<br>
><br>
> vega<br>
><br>
> --<br>
> ==================================================================================<br>
> Vega Lew ( weijia liu)<br>
> Graduate student<br>
> State Key Laboratory of Materials-oriented Chemical Engineering<br>
> College of Chemistry and Chemical Engineering<br>
> Nanjing University of Technology, 210009, Nanjing, Jiangsu, China<br>
> ******************************************************************************************************************<br>
</div></div>> Email: <a href="mailto:vegalew@gmail.com" target="_blank">vegalew@gmail.com</a> <mailto:<a href="mailto:vegalew@gmail.com" target="_blank">vegalew@gmail.com</a>><br>
<div>> Office: Room A705, Technical Innovation Building, Xinmofan Road 5#,<br>
> Nanjing, Jiangsu, China<br>
> ******************************************************************************************************************<br>
><br>
</div>> ------------------------------------------------------------------------<br>
><br>
> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
> <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
><br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>==================================================================================<br>Vega Lew ( weijia liu)<br>Graduate student<br>State Key Laboratory of Materials-oriented Chemical Engineering<br>
College of Chemistry and Chemical Engineering<br>Nanjing University of Technology, 210009, Nanjing, Jiangsu, China<br>******************************************************************************************************************<br>
Email: <a href="mailto:vegalew@gmail.com" target="_blank">vegalew@gmail.com</a><br></div></div><div><div></div><div class="h5">Office: Room A705, Technical Innovation Building, Xinmofan Road 5#, Nanjing, Jiangsu, China <br>
****************************************************************************************************************** <br>
</div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Somesh Kr. Bhattacharya<br>Post Doctoral Fellow<br>Room No. 263,<br>Leonardo Building,<br>The Abdus Salam International Centre for Theoretical Physics<br>Strada Costiera, 11<br>
I-34014 Trieste<br>Italy<br>Phone: +39-040-2240399<br><a href="http://portal.ictp.it/cmsp/members/postdoctoral-fellows/somesh-kumar-bhattacharya/">http://portal.ictp.it/cmsp/members/postdoctoral-fellows/somesh-kumar-bhattacharya/</a><br>