<div dir="ltr"><h3 class="gmail-iw">Dear Matthias, </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. Below is the input for DFT+U. &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 I will give a try to print the detail orbital occupation. Thanks, Sudhir <div> <div class="gmail_extra"> <div class="gmail_quote">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: <blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> <div lang="EN-GB"> <div class="gmail-m_-7467694623391386549WordSection1"> Dear Sudhir 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. Best regards Matthias 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 <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 </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+uns...@googlegroups.com" target="_blank">cp2k+unsubscribe@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>. </div> </div><div class="gmail-HOEnZb"><div class="gmail-h5"> -- 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+uns...@googlegroups.com" target="_blank">cp2k+unsubscribe@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>. </div></div></blockquote></div> -- <div class="gmail_signature"><div dir="ltr"><div>Sudhir Kumar Sahoo C/O - Prof. Thomas D. Kühne Post-doc Researcher <div>Department of Chemistry </div><div>University of Paderborn </div><div> </div></div></div></div> </div></div></div>