<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <meta name="Generator" content="Microsoft Word 14 (filtered medium)"> <style></style> </head> <body lang="EN-GB" link="blue" vlink="purple"> <div class="WordSection1"> Dear Sudhir<o:p></o:p> <o:p> </o:p> You enforce a smearing of the projected d orbital occupations until an SCF convergence threshold of 1.0E-8 is reached (and for at least 10 SCF iteration steps). This will imply that the occupation constraint is never released during the SCF iteration loop unless your SCF convergence limit is even smaller than 1.0E-8 (well, you possibly never get there). That is also most likely the reason for the persistence of the warning. I suggest to comment/remove the EPS_SCF 1.0E-8 line. After 10 SCF iteration steps the SCF cycle will then evolve without occupation constraint. You may even set “&ENFORCE_OCCUPATION off” and drop the occupation constraint completely, since you do not want to enforce any specific d orbital occupation pattern obviously as you are smearing over all d orbitals anyway. You may also drop the U_RAMPING unless that is crucially needed to achieve convergence.<o:p></o:p> <o:p> </o:p> Best regards<o:p></o:p> <o:p> </o:p> Matthias<o:p></o:p> <o:p> </o:p> From: cp...@googlegroups.com [mailto:cp...@googlegroups.com] On Behalf Of Sudhir K. Sahoo Sent: 18 September 2017 16:40 To: cp...@googlegroups.com Subject: Re: [CP2K:9424] DFT+U Calculation<o:p></o:p> <o:p> </o:p> <div> <h3>Dear Matthias, <o:p></o:p></h3> Many thanks for your reply and explanation. I have done one calculation with a lower value of U (2 eV), but got the same warning. <o:p></o:p> Below is the input for DFT+U. <o:p></o:p> &DFT_PLUS_U T L 2 U_MINUS_J [eV] 2.0 INIT_U_RAMPING_EACH_SCF U_RAMPING [eV] 0.10 &ENFORCE_OCCUPATION EPS_SCF 1.0E-8 MAX_SCF 10 ORBITALS -2 -1 0 1 2 SMEAR &END ENFORCE_OCCUPATION &END DFT_PLUS_U<o:p></o:p> I will give a try to print the detail orbital occupation. <o:p></o:p> Thanks, <o:p></o:p> Sudhir <o:p></o:p> <div> <o:p> </o:p> <div> <o:p> </o:p> <div> On Mon, Sep 18, 2017 at 4:03 PM, Krack Matthias (PSI) <<a href="mailto:matthia...@psi.ch" target="_blank">matthia...@psi.ch</a>> wrote:<o:p></o:p> <div> <div> Dear Sudhir<o:p></o:p> <o:p></o:p> This warning pops up when the Mulliken population analysis (MPA) returns an orbital occupation larger than the maximum number of electrons allowed for the actual l value specified in the +U section, e.g. more than 10 electrons are assigned by the MPA to the Cu 3d (l=2) orbitals. This can be an indicator for a too large U value enforcing a too high on-site localisation. Please, note that the U(eff) values successfully employed by PW codes are not transferable to CP2K. The corresponding U values for CP2K are usually smaller, roughly by a factor of two, e.g. with respect to their impact on the band gap. You can specifically print the detailed orbital occupations of the +U atoms using this <a href="https://manual.cp2k.org/cp2k-4_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/PLUS_U.html" target="_blank"> CP2K print key</a> with “QS_SCF 1” in its &EACH section. The warning appears easily for systems like ZnO with a 3d10 configuration when the MPA assigns for instance 10.04 electrons to the Zn 3d orbitals of some Zn atoms. The warning can then be ignored possibly for such spurious violations of the maximum orbital occupation. It is up to the user eventually to decide, if the obtained results are still physically meaningful. There is always the option to employ another method when DFT+U troubles, e.g. DFT with hybrid functionals.<o:p></o:p> <o:p></o:p> Best regards<o:p></o:p> <o:p></o:p> Matthias<o:p></o:p> <o:p></o:p> From: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> [mailto:<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>] On Behalf Of sudhir kumar Sahoo Sent: 18 September 2017 15:03 To: cp2k Subject: [CP2K:9422] DFT+U Calculation<o:p></o:p> <o:p></o:p> <div> Dear cp2k users, I am using DFT+U method for modelling CuInSe2 (U=5.0 eV for 3d electrons of Cu), particularly, using Mulliken method as implemented in cp2k program. However, a warning is repeatedly written in the output file. *** WARNING in dft_plus_u.F:1455 :: DFT+U energy contibution is negative *** *** possibly due to unphysical Mulliken charges. Check your input, if *** *** this warning persists or try a different method! *** After reading the cp2k google group discussion, I have tried to perform DFT+U calculations as explained below, but each time I got the same warning in the output file. 1) Restarting the calculation U=5 eV, from an optimized wavefunction obtained from U=0 eV. 2) Point 1 + u_ramping option 3) Point 2 + &enforce_occupation Can someone please help me in this. Thanking You, Sudhir<o:p></o:p> </div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsu...@googlegroups.com" target="_blank">cp2k+unsu...@googlegroups.com</a>. To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank"> cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank"> https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout" target="_blank"> https://groups.google.com/d/optout</a>.<o:p></o:p> </div> </div> <div> <div> -- You received this message because you are subscribed to a topic in the Google Groups "cp2k" group. To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/t0tJh9KSCBs/unsubscribe" target="_blank"> https://groups.google.com/d/topic/cp2k/t0tJh9KSCBs/unsubscribe</a>. To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+unsu...@googlegroups.com" target="_blank"> cp2k+unsu...@googlegroups.com</a>. To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank"> cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank"> https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout" target="_blank"> https://groups.google.com/d/optout</a>.<o:p></o:p> </div> </div> </div> -- <o:p></o:p> <div> <div> <div> Sudhir Kumar Sahoo C/O - Prof. Thomas D. Kühne Post-doc Researcher <o:p></o:p> <div> Department of Chemistry<o:p></o:p> </div> <div> University of Paderborn <o:p></o:p> </div> <div> <o:p> </o:p> </div> </div> </div> </div> </div> </div> </div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsu...@googlegroups.com">cp2k+unsu...@googlegroups.com</a>. To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k">https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout">https://groups.google.com/d/optout</a>.<o:p></o:p> </div> </body> </html>