Teo,<div><br></div><div>  Your answers are really helpful. I visualized the CP2K electrostatic potential cube file, and atomic unit makes sense.</div><div>  </div><div> Thanks.<br><br></div><div>Hanning</div><div><br></div>
<div><br><div class="gmail_quote">On Mon, Apr 6, 2009 at 10:40 AM, Teodoro Laino <span dir="ltr"><<a href="mailto:teodor...@gmail.com">teodor...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
If not specified in the on-line input reference, atomic units have to be<br>
assumed.<br>
Several keywords support different output units. If V_HARTREE_CUBE has<br>
not yet<br>
this feature (you can check again on the reference input on-line) you<br>
need just to wait.<br>
Very probably in the next future it will be available.<br>
<br>
Regards,<br>
Teo<br>
<div><div></div><div class="h5"><br>
NUCP2K wrote:<br>
> Dear CP2K group,<br>
><br>
>   I am a new user of CP2K. I obtained an electrostatic potential cube<br>
> file using directive %V_HARTREE_CUBE% after the wavefunction<br>
> optimization. I am wondering what is the energy unit in the potential<br>
> cube file. Is there any way to change the potential output unit?<br>
> Please see below for my input file.<br>
><br>
>   Thank you very much.<br>
><br>
> Hanning Chen<br>
> Postdoctoral Fellow<br>
> Northwestern University, Evanston Illinois 60208<br>
><br>
> Input File:<br>
><br>
>  &FORCE_EVAL<br>
>   METHOD Quickstep<br>
>   &DFT<br>
>     &MGRID<br>
>       CUTOFF 280<br>
>     &END MGRID<br>
>     &QS<br>
>       EPS_DEFAULT 1.0E-10<br>
>     &END QS<br>
>     &SCF<br>
>       SCF_GUESS RESTART<br>
>       EPS_SCF 1.0E-6<br>
>       MAX_SCF 1000<br>
>       &OT<br>
>        MINIMIZER DIIS<br>
>        PRECONDITIONER FULL_KINETIC<br>
>       &END OT<br>
>     &END SCF<br>
>     &XC<br>
>       &XC_FUNCTIONAL PADE<br>
>       &END XC_FUNCTIONAL<br>
>     &END XC<br>
>     &PRINT<br>
>      &TOT_DENSITY_CUBE ON<br>
>      &END TOT_DENSITY_CUBE<br>
>      &V_HARTREE_CUBE ON<br>
>      &END V_HARTREE_CUBE<br>
>     &END PRINT<br>
>   &END DFT<br>
>   &SUBSYS<br>
>     &CELL<br>
>       ABC 30.0 30.0 30.0<br>
>     &END CELL<br>
>     &COORD<br>
>       ..................<br>
>     &END COORD<br>
>     &KIND C<br>
>       BASIS_SET SZV-GTH-PADE<br>
>       POTENTIAL GTH-BLYP-q4<br>
>     &END KIND<br>
>   &END SUBSYS<br>
> &END FORCE_EVAL<br>
> &GLOBAL<br>
>   PROJECT TESTING<br>
>   RUN_TYPE ENERGY<br>
> &END GLOBAL<br>
><br>
><br>
> ><br>
><br>
<br>
<br>
<br>
</div></div></blockquote></div><br></div>