Hi Juerg,<div><br></div><div>  Thanks for your advice. I will try to cast my potential into GTH analytical form to save coding work.</div><div><br></div><div>Hanning</div><div><br></div><div><br><div class="gmail_quote">On Sun, Sep 13, 2009 at 10:54 AM, Juerg Hutter <span dir="ltr"><<a href="mailto:hut...@pci.uzh.ch">hut...@pci.uzh.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><br>
<br>
Hi<br>
<br>
if you can cast you potential into the form of the GTH<br>
pseudopotentials (see litterature for analytic form)<br>
your set. However, if this is not the case or<br>
if the interaction is also state<br>
specific then you have to do some programming.<br>
<br>
regards<br>
<br>
Juerg Hutter<br>
<br>
----------------------------------------------------------<br>
Juerg Hutter                   Phone : ++41 44 635 4491<br>
Physical Chemistry Institute   FAX   : ++41 44 635 6838<br>
University of Zurich           E-mail: <a href="mailto:hut...@pci.uzh.ch">hut...@pci.uzh.ch</a><br>
Winterthurerstrasse 190<br>
CH-8057 Zurich, Switzerland<br>
----------------------------------------------------------<br>
<div><div></div><div class="h5"><br>
<br>
On Thu, 10 Sep 2009, Hanning Chen wrote:<br>
<br>
> Dear CP2K community,<br>
>  I need to define an unusual electron - nuclei interaction, instead of the<br>
> usual Coulomb's law. Specifically, the nuclei-electron energy of a given<br>
> electron orbital does not follow the equation:   <e|1/(R-e)|e>, where |e> is<br>
> the electron orbital and (R-e) is the vector connecting nuclei and electron.<br>
><br>
>  I am wondering what subroutines of CP2K handle the calculation of<br>
> nuclei-electron (or core-electron) energy ? Is it possible to define a new<br>
> type of such interaction in the input file without hacking the source code?<br>
><br>
>  Thanks for any advice.<br>
><br>
> Hanning<br>
><br>
> Department of Chemistry<br>
> Northwestern University<br>
> Evanston, IL 60201<br>
><br>
> ><br>
><br>
<br>
</div></div><br>
</blockquote></div><br></div>