<div dir="ltr">Hi,<div><br></div><div>there are going to be several competing / contradictory ways of doing these calcs.</div><div><br></div><div>I have only done them for a non-periodic system (slab), then you can apply the external field, as you are doing, and print of the density as a cube file. You can relate the change of density with applied field to get the dielectric 'constant'.</div><div><br></div><div>For a periodic system, you would need to use a periodic electric field and presumably also look at the density response and use electrostatics.</div><div><br></div><div>I haven't used the new implicit psolver method, my feeling is that the dielectric there is more likely to be an applied boundary condition, rather than an output, but it is very flexible, so maybe it can be made to do what you want. It will only be activated if you use the implicit solver, which is probably why you see no output.</div><div><br></div><div>Matt</div><div><br><br>On Monday, April 24, 2017 at 2:35:15 PM UTC+1, Fernan Saiz wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<br>Thanks for the reply. In this case, I guess CP2k should be able to write a CUBE output file with epsilon(x,y,z) once a external potential is applied, shouldn't it?<br><br>I can implement an external potential by writing in the input file: <br><br>&FORCE_EVAL<br>  &DFT<br>    &EXTERNAL_POTENTIAL<br>      FUNCTION (A/B)*X<br>      PARAMETERS A B<br>      VALUES [eV] 5.0 [angstrom] 1000.0<br>    &END <br>  &END <br>&END <br><br>But to print the dielectric constant epsilon(x,y,z) in CUBE files, I am writing the following commands:<br><br>&FORCE_EVAL<br>  &DFT<br>    &PRINT<br><br>      &IMPLICIT_PSOLVER<br>        &DIELECTRIC_CUBE<br>          &EACH<br>            JUST_ENERGY 1<br>          &END EACH<br>          FILENAME diel/data_diel<br>        &END DIELECTRIC_CUBE<br>      &END IMPLICIT_PSOLVER<br><br>    &END <br>  &END <br>&END <br><br>Unfortunately, this second paragraph of commands does not write anything in my working directory. Am I using the wrong commands?<br><br>Regards, <br> - Fernan<br><br>On Tuesday, April 18, 2017 at 3:09:50 PM UTC+1, jgh wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">Hi
<br>
<br>I think all the parts needed to do such a calculation are available,
<br>but there is no such keyword to let the code do it automatically.
<br>
<br>regards
<br>
<br>Juerg
<br>------------------------------<wbr>------------------------------<wbr>--
<br>Juerg Hutter                        <wbr> Phone : ++41 44 635 4491
<br>Institut für Chemie C                FAX   : ++41 44 635 6838
<br>Universität Zürich                   E-<wbr>mail: <a rel="nofollow">hut...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>------------------------------<wbr>------------------------------<wbr>---
<br>
<br>-----<a rel="nofollow">cp...@googlegroups.com</a> wrote: -----To: cp2k <<a rel="nofollow">cp...@googlegroups.com</a>>
<br>From: Fernan Saiz 
<br>Sent by: <a rel="nofollow">cp...@googlegroups.com</a>
<br>Date: 04/15/2017 12:06PM
<br>Subject: [CP2K:8947] dielectric function calculation
<br>
<br>Hello all,
<br>I would like to ask if CP2k is able to calculate the high-frequency dielectric function epsilon(z) as a function of a given direction z as formulated by Fu. et al. PRB 57, 6967 (1998). As a further reference, the all-electron code CRYSTAL14 has it implemented using a finite field perturbation method.
<br>
<br> Thanks,
<br> - Fernan Saiz, PhD
<br>Department of Chemistry
<br>Imperial College London
<br>
<br>
<br>
<br>
<br>-- 
<br>
<br>You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>
<br>To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp2k+...@googlegroups.com</a>.
<br>
<br>To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.
<br>
<br>Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.
<br>
<br>For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>
<br>
<br></blockquote></div></blockquote></div></div>