<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=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:"Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Balloon Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.BalloonTextChar
{mso-style-name:"Balloon Text Char";
mso-style-priority:99;
mso-style-link:"Balloon Text";
font-family:"Tahoma","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:85.05pt 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>I am trying to apply Hubbard U (DFT+U) to Ce in cerium based oxides including CeO2 and Ce2O3. However, as you see from header of Ce psp (see [I] below), our Ce psp (RRKJ, generated by OPIUM) is designed to have zero 4f occupation. When I applied +U to Ce using QE V4.3.1, I got following error message:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>from offset_atom_wfc : error # 9<o:p></o:p></p><p class=MsoNormal>wrong offset<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>So, I searched previous threads and I figured out that I needed to modify QE code (http://www.democritos.it/pipermail/pw_forum/2010-July/017622.html). <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I did modification and recompiled QE code. It looked fine but I found weird energetics for reduction of CeO2 to Ce2O3 in comparison with data from PRB 75, 035109 (2007). So, I looked into output files and found out that occupation matrix for both CeO2 and Ce2O3 are somewhat unphysical (at least for me):<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>*U_eff=0.01 eV <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Tr[ns(na)] <o:p></o:p></p><p class=MsoNormal>CeO2 = 1.1690424<o:p></o:p></p><p class=MsoNormal>Ce2O3 = 1.3908759<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>*U_eff=3.0 eV <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Tr[ns(na)] <o:p></o:p></p><p class=MsoNormal>CeO2 = 0.9963736<o:p></o:p></p><p class=MsoNormal>Ce2O3 = 1.2650403 <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>*U_eff=9.0 eV <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Tr[ns(na)] <o:p></o:p></p><p class=MsoNormal>CeO2 = 0.5972449 <o:p></o:p></p><p class=MsoNormal>Ce2O3 = 1.1056369<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I thought that even with non-orthogonal terms, at least Tr[ns(na)] for CeO2 should be very close to 0.0 as Ce exists as Ce4+ in CeO2 (with ionic picture) and should have empty 4f state. With analyzing Ce psp itself, I wanted to know whether the modification of QE discussed above is safe or not. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Answer of Dr. Sclauzero to my question was as follows:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>"Briefly, I think that your problem stems exactly from that fix which I proposed and is now included in the latest PWscf version. That fix will actually introduce a problem if the atomic wave function which you want to use for making projections have zero occupations. I thought that this would never be the case, but probably one should not exclude it a priori. A further problem arises if you generate such a pseudopotential (Ce with zero 4f occupation) and want to use it for making the LDA+U projections on the empty states, since the LD1 code usually saves in the PP file only the wavefunctions with nonzero occupations. First, you should check if you PP has the 4f wavefunctions. If it does, you can remove the check upf(nt)%oc(n) > 0.D0 as a quick fix. If not, then you need to regenerate the PP".<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Based on his comment, I did some test calculations after modifying our Ce psp--simply made zero 4f occupation as nonzero, 0.01 (In (I), 430 0.00 --> 430 0.01). However, there is no total E change between the calculation with Ce psp of zero 4f occupation and that with nonzero 4f occupation (E difference <0.001/CeO2). Also, change in Tr[ns(na)] is minute (e.g. 0.9938046 --> 0.9963736)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>My questions to QE users are:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>- Maybe I need to fill 4f more? <o:p></o:p></p><p class=MsoNormal>- Do you have any explanations about why is occupation matrix sum close to 1.0 not 0.0? Again, while there are some non-orthogonal terms, I think that for CeO2 Tr[ns(na)] for Ce should be close to zero not one. Did I miss something?<o:p></o:p></p><p class=MsoNormal>- How about using starting_ns_eigenvalue? By the way, while I tried to use this, I couldn't understand instruction. For instance, in example below ([II]), if I want to make initial occupations as zero, how can I do? In starting_ns_eigenvalue(m,ispin,I), I can see that I should have starting_ns_eigenvalue(m,1,1) and starting_ns_eigenvalue(m,2,1) if U is applied to spin-up and spin-down for atomic species 1. Then how about m? m=1,2,3,4,5,6, and 7? If not, what will be desirable values for m?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Thank you for your advice and comments in advance.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>P.S. I post this while there is existing post regarding the questions above because:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>1) The title of the existing post is not suitable and will never be easy for QE users to search. <o:p></o:p></p><p class=MsoNormal>2) I am not comfortable with adding threads under the post which contains my personal emails (not for the forum) which violate my privacy. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>++++++++++++++++++++++++++++++<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>[I]<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>[Atom]<o:p></o:p></p><p class=MsoNormal>Ce<o:p></o:p></p><p class=MsoNormal>15 # norb: number of orbitals<o:p></o:p></p><p class=MsoNormal>100 2.00 -643.0 # nlm occ eigen(- means auto-generate)<o:p></o:p></p><p class=MsoNormal>200 2.00 -77.4<o:p></o:p></p><p class=MsoNormal>210 6.00 -68.4<o:p></o:p></p><p class=MsoNormal>300 2.00 -10.0<o:p></o:p></p><p class=MsoNormal>310 6.00 -6.0<o:p></o:p></p><p class=MsoNormal>320 10.0 -6.4<o:p></o:p></p><p class=MsoNormal>400 2.00 -6.4<o:p></o:p></p><p class=MsoNormal>410 6.00 -6.0<o:p></o:p></p><p class=MsoNormal>420 10.0 -6.0<o:p></o:p></p><p class=MsoNormal>500 2.00 -2.4<o:p></o:p></p><p class=MsoNormal>510 6.00 -2.0<o:p></o:p></p><p class=MsoNormal>600 0.00 -1.0<o:p></o:p></p><p class=MsoNormal>610 0.00 -1.0<o:p></o:p></p><p class=MsoNormal>520 0.00 -1.0<o:p></o:p></p><p class=MsoNormal>430 0.00 -1.9<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>+++++++++++++++++++++++++++++++<o:p></o:p></p><p class=MsoNormal>[II]<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>U( 1) = 3.0000 U( 2) = 0.0000<o:p></o:p></p><p class=MsoNormal>alpha( 1) = 0.0000 alpha( 2) = 0.0000<o:p></o:p></p><p class=MsoNormal>atom 1 Tr[ns(na)]= 2.0000000<o:p></o:p></p><p class=MsoNormal>atom 1 spin 1<o:p></o:p></p><p class=MsoNormal>eigenvalues: 0.1428571 0.1428571 0.1428571 0.1428571 0.1428571 0.1428571 0.1428571<o:p></o:p></p><p class=MsoNormal> eigenvectors<o:p></o:p></p><p class=MsoNormal> 1 1.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000<o:p></o:p></p><p class=MsoNormal> 2 0.0000000 1.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000<o:p></o:p></p><p class=MsoNormal> 3 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 0.0000000 0.0000000<o:p></o:p></p><p class=MsoNormal> 4 0.0000000 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 0.0000000<o:p></o:p></p><p class=MsoNormal> 5 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000<o:p></o:p></p><p class=MsoNormal> 6 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 0.0000000<o:p></o:p></p><p class=MsoNormal> 7 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000<o:p></o:p></p><p class=MsoNormal> occupations<o:p></o:p></p><p class=MsoNormal> 0.143 0.000 0.000 0.000 0.000 0.000 0.000<o:p></o:p></p><p class=MsoNormal> 0.000 0.143 0.000 0.000 0.000 0.000 0.000<o:p></o:p></p><p class=MsoNormal> 0.000 0.000 0.143 0.000 0.000 0.000 0.000<o:p></o:p></p><p class=MsoNormal> 0.000 0.000 0.000 0.143 0.000 0.000 0.000<o:p></o:p></p><p class=MsoNormal> 0.000 0.000 0.000 0.000 0.143 0.000 0.000<o:p></o:p></p><p class=MsoNormal> 0.000 0.000 0.000 0.000 0.000 0.143 0.000<o:p></o:p></p><p class=MsoNormal> 0.000 0.000 0.000 0.000 0.000 0.000 0.143<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><i><span style='font-size:10.0pt;font-family:"Arial","sans-serif"'>-------<br>Chan-Woo Lee, Ph.D.<o:p></o:p></span></i></p><p class=MsoNormal><i><span style='font-size:10.0pt;font-family:"Arial","sans-serif"'><o:p> </o:p></span></i></p><p class=MsoNormal><i><span style='font-size:10.0pt;font-family:"Arial","sans-serif"'>Department of Chemistry<br>University of Pennsylvania<br>231 South 34th Street<br>Philadelphia, PA 19104-6323 <br>Phone: 1-215-898-3564 (Office)<br></span></i><span style='font-size:10.0pt;font-family:"Arial","sans-serif"'>Email: <a href="mailto:leechanw@sas.upenn.edu"><span style='color:blue'>leechanw@sas.upenn.edu</span></a> / <a href="mailto:cwandtj@gmail.com" target="_blank"><span style='color:blue'>cwandtj@gmail.com</span></a><o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>