<div dir="ltr"><div>I suggest you take your time to read this if you have not yet, I am also learning like you. I might try to experiment with you system to see what I get with same pseudopotentials as you used.<br><br><a href="http://urania.sissa.it/xmlui/handle/1963/2507">http://urania.sissa.it/xmlui/handle/1963/2507</a><br><br></div>and then look at this tutorial<br><br><a href="http://www.iiserpune.ac.in/~smr2626/hands_on/week1/july4/">http://www.iiserpune.ac.in/~smr2626/hands_on/week1/july4/</a><br></div><div class="gmail_extra"><br><div class="gmail_quote">2016-10-16 11:54 GMT+02:00 Phanikumar Pentyala <span dir="ltr"><<a href="mailto:phani12.chem@gmail.com" target="_blank">phani12.chem@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:times new roman,serif">Hello QE users<br><br></div><div class="gmail_default" style="font-family:times new roman,serif">I want to find U value for CeO2 which is having 81 atoms. I already calculated this value by using linear response (<a href="http://hjklol.mit.edu/content/calculating-hubbard-u" target="_blank">http://hjklol.mit.edu/<wbr>content/calculating-hubbard-u</a>)<wbr>. I got the U value initially in negative numbers in this case <b>diago_thr_init</b> value was default, so that I changed <b>diago_thr_init</b> to 10^-11 then I got U value as 0.00881. From literature I found this value should be around 4.5 - 7eV  (Phy. Rev. B 75, 035109, 2007͒). I have some doubts regarding this calculation<br><br>Is that correct, occupations means I took value in the name of <b>"N of occupied +U levels" </b>(i.e sum of Tr[ns(na)] values on all 27 Ce atoms)?<br><br></div><div class="gmail_default" style="font-family:times new roman,serif">What is the unit for value I got (0.00881)? (Ry or eV or any other )<br><br></div><span style="font-family:times new roman,serif">Is U parameter value depend on number of atoms in unit cell?</span><br><div class="gmail_default" style="font-family:times new roman,serif"><br><div class="gmail_default" style="font-family:times new roman,serif">Is there any other mistake I done in my input file ? (input file attached)<br></div><br><b>Conditions I applied in my calculations:</b><br></div><div class="gmail_default" style="font-family:times new roman,serif">Surface calculation with vacuum 15A<br></div><div class="gmail_default" style="font-family:times new roman,serif">First three layers relaxed (two oxygen and one ceria layers) remaining are frozen<br></div><div class="gmail_default" style="font-family:times new roman,serif"><br><br>Thanks in advance and I will be greatly appreciated for your help <br></div><div class="gmail_default" style="font-family:times new roman,serif"><br></div><div class="gmail_default" style="font-family:times new roman,serif">(input file and results of my calculation attached)<br></div><div class="gmail_default" style="font-family:times new roman,serif"><br><br></div><div class="gmail_default" style="font-family:times new roman,serif">Phanikumar<br></div><div class="gmail_default" style="font-family:times new roman,serif">Research Scholar <br></div><div class="gmail_default" style="font-family:times new roman,serif">Indian Institute of Technology Kharagpur<br></div><div class="gmail_default" style="font-family:times new roman,serif">India<br></div></div>
<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div>