[Pw_forum] CP simulation with ensemble-dft algorithm

pascal boulet pascal.boulet at univ-provence.fr
Sun Jun 19 10:21:47 CEST 2011


Dear Giovanni,

Humm... hard to say because the data are not printed in the output file.
Here is what I get, for instance:

>>>>
  nfi    ekinc  temph  tempp        etot      enthal       econs
econt    vnhh   xnhh0    vnhp   xnhp0
        nfi   tempp            E -T.S-mu.nbsp         +K_p #Iter
Step    20       0   -502.49200   -502.62143   -502.62143     6
 cg_sub: missed minimum, case 1, iteration           1
 cg_sub: missed minimum, case 1, iteration           2
 cg_sub: missed minimum, case 1, iteration           3
 cg_sub: missed minimum, case 1, iteration           4
 cg_sub: missed minimum, case 1, iteration           5
NOTE: eigenvalues are not computed without ortho
Occupations  :
2.000000 2.000000 2.000000 2.000000 2.000000 2.000000 1.999999 1.999999
1.999999 1.999998
 1.999997 1.999995 1.999994 1.999989 1.999983 1.999977 1.999958 1.999933
1.999884 1.999793 ... 0.001473 0.000546 0.000502 0.000278 0.000134
0.000079 0.000033 0.000019 0.000000
        nfi   tempp            E -T.S-mu.nbsp         +K_p #Iter
Step    21       0   -502.49200   -502.62143   -502.62143     6

<<<<<
etc.


So, there is nothing about ekinc, temph, tempp, ...


However, the data "nfi    ekinc  temph  tempp        etot     enthal
   econs       econt    vnhh   xnhh0    vnhp   xnhp0" look
like those given in the file xxx.evp (except for one missing column).
If this is right, then in this file ekinc = 0. during the whole
simulation. That's strange!
So, I guess there is something wrong somewhere in my input.

Pascal


> Date: Fri, 17 Jun 2011 10:27:09 +0200 (CEST)
> From: Giovanni La Penna <glapenna at iccom.cnr.it>
> Subject: Re: [Pw_forum] CP simulation with ensemble-dft algorithm
> To: PWSCF Forum <pw_forum at pwscf.org>
> Message-ID: <alpine.LNX.2.00.1106171015110.8707 at biurca.iccom.cnr.it>
> Content-Type: TEXT/PLAIN; format=flowed; charset=US-ASCII
>
>
> Dear Pascal,
>
> The input shows that you are trying to
> carry the system to 627 K starting with
> random velocities and a time-step of 20 au (about 0.5 fs).
> I expect a huge Kel after a few time-steps, due
> to huge current atomic velocities, not manageable
> with the given time-step.
> Please, tell me (us) if this won't happen.
>
>               Giovanni La Penna
>
>

    *
    * Anglais
    * Français

    * Anglais
    * Français

 <javascript:void(0);>



More information about the users mailing list