<html><body><div style="color:#000; background-color:#fff; font-family:tahoma, new york, times, serif;font-size:12pt"><div><span><br></span></div><div><br></div> <div style="font-family: tahoma,new york,times,serif; font-size: 12pt;"> <div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight: bold;">From:</span></b> Hajar Nejati <hajar.nejatipoor@yahoo.com><br> <b><span style="font-weight: bold;">To:</span></b> A Mailing list for WIEN2k users <wien@zeus.theochem.tuwien.ac.at> <br> <b><span style="font-weight: bold;">Sent:</span></b> Tuesday, February 26, 2013 11:17 AM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [Wien] mbj scf error<br> </font> </div> <br><div id="yiv498267432"><div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: tahoma,new york,times,serif; font-size:
12pt;">dear Dr. Blaha<br>I do mbj calculations for other case (GeO2) without problem.<br>there is in the grapphane.output0 NaNs. why? what should I do?<br>summary of graphane.output0:<br> XC-potentials inside spheres (XCPOT1)<br> Gauss-Legendre grid of 14 x 27<br> INDEX V-XC FIT DIFFERENCE<br> 1 NaN NaN NaN<br> 2 NaN
NaN NaN<br> 4 NaN NaN NaN<br> ...<br> 378 NaN NaN NaN<br>ATOM 1 AT RMT: SIGMA OF V-XC FIT: NaN<br>ATOM 1 LARGEST SIGMA AT R( 2)= 0.000 and MAX DIFF: 0.00000E+00
0.00000E+00<br> Gauss-Legendre grid of 14 x 27<br> INDEX V-XC FIT DIFFERENCE<br> 1 NaN NaN NaN<br> 2 NaN NaN NaN<br> ....<br> 378 NaN NaN NaN<br>ATOM 2 AT RMT: SIGMA OF V-XC
FIT: NaN<br>ATOM 2 LARGEST SIGMA AT R( 2)= 0.000 and MAX DIFF: 0.00000E+00 0.00000E+00<br><br> SELECTED FOURIERCOEFF. OF V-XC <br> 0 0 0 NaN 0.00000E+00 NaN 0.00000E+00<br> 0 0
1 NaN NaN NaN NaN<br> 0 0 2 NaN NaN NaN NaN<br> ....<br> 0
0191 NaN NaN NaN NaN<br> v-mean, v0,u0 (NaN,NaN) (NaN,NaN) (552.845591947925,0.000000000000000E+000)<br> (NaN,0.000000000000000E+000)<br>1<br> ------------------------------<br> T O T A L E N E R G
Y<br> 1. DENSITY-POTENTIAL INTEGRALS (H)<br> SPHERE= 1 C <br> L= 0 : NaN SPHERE CHARGE<br> ....<br> L= 6 : NaN SPHERE CHARGE<br> TOTAL= NaN<br> SPHERE=
2 C <br><br> L= 0 : NaN SPHERE CHARGE<br> . ...<br> L= 6 : NaN SPHERE CHARGE<br> TOTAL= NaN<br> SPHERE= 3 H <br><br> L= 0 : NaN SPHERE
CHARGE<br> L= 1
: NaN SPHERE CHARGE<br> ...<br> L= 6
: NaN SPHERE CHARGE<br> TOTAL= NaN<br> SPHERE= 4 H <br> L= 0 : NaN SPHERE CHARGE<br> L= 1 : NaN SPHERE CHARGE<br> ...<br> L= 6 : NaN SPHERE
CHARGE<br> TOTAL= NaN<br> PLANE WAVE
CONTRIBUTION<br><br> TOTAL= NaN<br>2. CONSTANT TERMS (H)<br> SPHERE 1 C = -8.212096<br> SPHERE 2 H = -0.451525<br> --------------------------------------------------------------------------------<br><br> TOTAL VALUE = NaN (H)<br>:DEN : DENSITY INTEGRAL =
NaN (Ry)<br> ELS_POTENTIAL_AT Z=0 and Z=0.5: 0.17793 -1.77017<br> ELS_POTENTIAL_AT Y=0 and Y=0.5: 0.00000 0.00000<br>:VZERO:v0,v0c,v0x NaN 0.17793 NaN v5,v5c,v5x NaN -1.77017 NaN<br>:VZERY:v0,v0c,v0x NaN 0.00000 NaN v5,v5c,v5x NaN 0.00000 NaN<br>:VZERX:v0,v0c,v0x NaN 0.35728 NaN v5,v5c,v5x NaN -0.11638 NaN<br><br><br><div><br></div> <div
style="font-family: tahoma,new york,times,serif; font-size: 12pt;"> <div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight: bold;">From:</span></b> Peter Blaha <pblaha@theochem.tuwien.ac.at><br> <b><span style="font-weight: bold;">To:</span></b> A Mailing list for WIEN2k users <wien@zeus.theochem.tuwien.ac.at> <br> <b><span style="font-weight: bold;">Sent:</span></b> Tuesday, February 26, 2013 11:01 AM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [Wien] mbj scf error<br> </font> </div> <br>
Have you tried another (simple) case, like Si, .... ?<br><br>Do again a PBE calculation until self-consistency. Once there is a <br>problem with mBJ, one cannot easily recover since the mBJ potential <br>depends on the potential of the previous iteration.<br><br>Change mixing in case.inm to PRATT with mixing 0.03 for about 15 <br>iterations. Then switch to MSR1. (I should note, that the newest (net <br>yet released) version of mixer seems to handle also mBJ pretty well in <br>MSR1).<br><br>An if a problem occurs, we need more info. Don't rerun it 5 times. As <br>mentioned before, with mBJ this will not work but make things worse. The <br>failure in lapw1 is usually only a follow up problem with roots <br>somewhere before.<br><br>Check case.output0. Are there NaNs ?<br><br><br>Am 26.02.2013 07:43, schrieb Hajar Nejati:<br>> hi dear ali<br>> I use wien12.<br>> I could not solve my problem up to now.<br>> can you help me,
please?<br>><br>><br>> ------------------------------------------------------------------------<br>> *From:* ali ghafari <<a rel="nofollow" ymailto="mailto:aaghafari@yahoo.com" target="_blank" href="mailto:aaghafari@yahoo.com">aaghafari@yahoo.com</a>><br>> *To:* A Mailing list for WIEN2k users <<a rel="nofollow" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br>> *Sent:* Monday, February 25, 2013 11:10 PM<br>> *Subject:* Re: [Wien] mbj scf error<br>><br>> Which version of Wien2k you are using?<br>><br>> Ali<br>><br>> ------------------------------------------------------------------------<br>> *From:* Hajar Nejati <<a rel="nofollow" ymailto="mailto:hajar.nejatipoor@yahoo.com" target="_blank" href="mailto:hajar.nejatipoor@yahoo.com">hajar.nejatipoor@yahoo.com</a>><br>> *To:* A Mailing list for
WIEN2k users <<a rel="nofollow" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br>> *Sent:* Monday, February 25, 2013 12:46 PM<br>> *Subject:* Re: [Wien] mbj scf error<br>><br>> Thank you very much dear Gavin.<br>><br>> ------------------------------------------------------------------------<br>> *From:* Gavin Abo <<a rel="nofollow" ymailto="mailto:gsabo@crimson.ua.edu" target="_blank" href="mailto:gsabo@crimson.ua.edu">gsabo@crimson.ua.edu</a>><br>> *To:* A Mailing list for WIEN2k users <<a rel="nofollow" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br>> *Sent:* Sunday, February 24, 2013 2:52 PM<br>> *Subject:* Re: [Wien] mbj scf error<br>><br>> These old posts might be helpful:<br>><br>>
http://zeus.theochem.tuwien.ac.at/pipermail/wien/2010-March/012700.html<br>>
http://zeus.theochem.tuwien.ac.at/pipermail/wien/2010-March/012789.html<br>><br>> On 2/24/2013 4:04 AM, Hajar Nejati wrote:<br>>> dear hossein<br>>> I calculate scf without parallel calculations, but the error repeat.<br>>> lapw1.error :<br>>> Error in LAPW1<br>>> Cholesky INFO = 1<br>>> 'SECLR4' - POTRF (Scalapack/LAPACK) failed.<br>>><br>>><br>>> ------------------------------------------------------------------------<br>>> *From:* hossien rahnama <<a rel="nofollow" ymailto="mailto:h_rahnamay@yahoo.com" target="_blank" href="mailto:h_rahnamay@yahoo.com">h_rahnamay@yahoo.com</a>><br>>> <mailto:<a rel="nofollow" ymailto="mailto:h_rahnamay@yahoo.com" target="_blank" href="mailto:h_rahnamay@yahoo.com">h_rahnamay@yahoo.com</a>><br>>> *To:* A Mailing list for WIEN2k users<br>>> <<a rel="nofollow"
ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br>>> <mailto:<a rel="nofollow" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br>>> *Sent:* Sunday, February 24, 2013 1:28 PM<br>>> *Subject:* Re: [Wien] mbj scf error<br>>><br>>> Dear Hajar,<br>>> Please try without parallelcalculations.<br>>><br>>> Dr. H.A.Rahnamaye Aliabad,<br>>> Assistant Professor,<br>>> Department of physics,Hakim Sabzevari University,<br>>> Sabzevar,Iran<br>>> Tel:+98-571-4003155<br>>> Fax:+98-571-4411161<br>>> Academic email:<a rel="nofollow" ymailto="mailto:rahnama@sttu.ac.ir" target="_blank" href="mailto:rahnama@sttu.ac.ir">rahnama@sttu.ac.ir</a> <mailto:email:<a rel="nofollow"
ymailto="mailto:rahnama@sttu.ac.ir" target="_blank" href="mailto:rahnama@sttu.ac.ir">rahnama@sttu.ac.ir</a>><br>>> *From:* Hajar Nejati <<a rel="nofollow" ymailto="mailto:hajar.nejatipoor@yahoo.com" target="_blank" href="mailto:hajar.nejatipoor@yahoo.com">hajar.nejatipoor@yahoo.com</a>><br>>> <mailto:<a rel="nofollow" ymailto="mailto:hajar.nejatipoor@yahoo.com" target="_blank" href="mailto:hajar.nejatipoor@yahoo.com">hajar.nejatipoor@yahoo.com</a>><br>>> *To:* wien2k (for send question) <<a rel="nofollow" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br>>> <mailto:<a rel="nofollow" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br>>> *Sent:* Sunday, February 24, 2013 1:14 PM<br>>> *Subject:*
[Wien] mbj scf error<br>>><br>>> dear all wien users<br>>> I want to calculate band structure of graphane using mbj method for<br>>> XC-potential.<br>>> At first, I copy case.inm_vresp in my directory(to<br>>> graphane.inm_vresp), then i
edit graphane.in0 (change NR2V to R2V) and<br>>> run single cycle.<br>>> I run a scf (run_lapw -p -ec 0.0001), after edition of graphane.in0<br>>> [option xc(default=13) to 28] and graphane.in0_grr [option xc to 50].<br>>> this error appear after one or two cycle and the scf stoped:<br>>> LAPW0 END<br>>> LAPW0 END<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> SECLR4 - Error<br>>> ...<br>>> please help me about this error.<br>>> hajar nejati<br>>><br>><br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank"
href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<mailto:<a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>> <a rel="nofollow" target="_blank" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>><br>><br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a> <mailto:<a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>> <a rel="nofollow" target="_blank"
href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>><br>><br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a> <mailto:<a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>> <a rel="nofollow" target="_blank" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>><br>><br>><br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank"
href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>> <a rel="nofollow" target="_blank" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>><br><br>-- <br>Peter Blaha<br>Inst.Materials Chemistry<br>TU Vienna<br>Getreidemarkt 9<br>A-1060
Vienna<br>Austria<br>+43-1-5880115671<br>_______________________________________________<br>Wien mailing list<br><a rel="nofollow" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br><a rel="nofollow" target="_blank" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br><br><br> </div> </div> </div></div></div><br><br> </div> </div> </div></body></html>