<p>Check the confidence levl, e.g. t-test (i am not an expert with this package). Xrd Reitfeld is as susceptible to GIGO as DFT.</p>
<div class="gmail_quote">On Sep 5, 2011 11:58 AM, &quot;shamik chakrabarti&quot; &lt;<a href="mailto:shamikphy@gmail.com">shamikphy@gmail.com</a>&gt; wrote:<br type="attribution">&gt; Dear Laurence Marks Sir,<br>&gt; <br>
&gt;      Yes this is an oxide and it is an insulator and probably I have made it<br>&gt; a metal by creating *cationic vacancies*. Sir I will let you know the detail<br>&gt; after checking it from case.output2. But Sir, please forgive my curiosity, I<br>
&gt; could not wait to ask you and also going a bit away from wien2k questions:<br>&gt; 1) Sir the goodness of fit (GOF) obtained from Rietveld analysis done<br>&gt; through XPERT Highscore plus : 1.0038<br>&gt;     and except this data we do not have any other proof to show that it<br>
&gt; really have cationic vacancies.<br>&gt;     Then Sir according to you we can not trust this data (considering GOF)<br>&gt; only to consider cationic vacancies?<br>&gt; <br>&gt; With best regards,<br>&gt; <br>&gt; <br>
&gt; <br>&gt; On Mon, Sep 5, 2011 at 9:31 PM, Laurence Marks &lt;<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a>&gt;wrote:<br>&gt; <br>&gt;&gt; Follow the science, not the DFT.<br>&gt;&gt;<br>&gt;&gt; 1) Unless you have a large grain size and synchotron data I really<br>
&gt;&gt; wonder about the logic of removing 1 atom in a 2x1x1 supercell.<br>&gt;&gt; Particularly if this is an oxide, have you made it a metal whereas it<br>&gt;&gt; should be an insulator? Does the fit pass significance tests etc? Do<br>
&gt;&gt; you have other data that implies vacancies?<br>&gt;&gt; 2) In terms of the DFT, I expect that you have made it a metal which<br>&gt;&gt; is harder to converge than an insulator. Check in case.output2 in both<br>&gt;&gt; cases (a quick, pseudo DOOS). If it really does have vacancies, you<br>
&gt;&gt; might need more than 7 kpts, and I suggest using TEMPS.<br>&gt;&gt;<br>&gt;&gt; 2011/9/5 shamik chakrabarti &lt;<a href="mailto:shamikphy@gmail.com">shamikphy@gmail.com</a>&gt;:<br>&gt;&gt; &gt; Dear wien2k users,<br>
&gt;&gt; &gt;         I have performed a 56 atom / unit cell calculation for a<br>&gt;&gt; &gt; ferrimagnetic material. It was well converged both in energy and charge<br>&gt;&gt; by<br>&gt;&gt; &gt; 60 iterations. 14 K points were used for this calculation. Then I have<br>
&gt;&gt; &gt; generated a 2x1x1 supercell having 112 atoms / supercell. This was done<br>&gt;&gt; to<br>&gt;&gt; &gt; achieve better experimental model as obtained through Rietveld analysis.<br>&gt;&gt; &gt; From Rietveld analysis it was also seen that for 112 atoms we can<br>
&gt;&gt; consider a<br>&gt;&gt; &gt; cation (Fe) vacancy and we have generated it by deleting one Fe atom from<br>&gt;&gt; &gt; the supercell. Final structure, hence, contains 111 atoms and all other<br>&gt;&gt; &gt; occupancies at different lattice sites are well matched with Rietveld<br>
&gt;&gt; &gt; analysis of experimental XRD. I am using 2 k points for this 111 atom<br>&gt;&gt; &gt; supercell. But even after 60 iterations the convergence achieved are as<br>&gt;&gt; &gt; follows:<br>&gt;&gt; &gt;    cc: 0.234...    cc was set to 0.001<br>
&gt;&gt; &gt;    ec: 0.063...    ec was set to 0.0001<br>&gt;&gt; &gt; Sir, my question are:<br>&gt;&gt; &gt; 1) As the supercell is merely an extension of the 56 atoms unit cell with<br>&gt;&gt; &gt; more accurate values of atomic occupancies as in the experimental<br>
&gt;&gt; structure,<br>&gt;&gt; &gt; why it should at all taking so much time to converge?<br>&gt;&gt; &gt; 2) We have created one defect by creating a Fe vacancy...is it that the<br>&gt;&gt; &gt; structures having defect take more time to converge?...but as far as our<br>
&gt;&gt; &gt; analysis until now..the structure should have one Fe vacancy per 112<br>&gt;&gt; &gt; atoms...<br>&gt;&gt; &gt; 3) As we have obtained converged solutions for 56 atom cell by using 14 k<br>&gt;&gt; &gt; points..is it that we should use 7K points for 111 atom supercell<br>
&gt;&gt; &gt; calculation?..is 2 K points are actually small in order to achieve<br>&gt;&gt; &gt; convergence of this structure?<br>&gt;&gt; &gt; Any response in this regard will be very helpful for us. Thanks in<br>&gt;&gt; advance.<br>
&gt;&gt; &gt; with regards,<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; --<br>&gt;&gt; &gt; Shamik Chakrabarti<br>&gt;&gt; &gt; Senior Research Fellow<br>&gt;&gt; &gt; Dept. of Physics &amp; Meteorology<br>&gt;&gt; &gt; Material Processing &amp; Solid State Ionics Lab<br>
&gt;&gt; &gt; IIT Kharagpur<br>&gt;&gt; &gt; Kharagpur 721302<br>&gt;&gt; &gt; INDIA<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; _______________________________________________<br>&gt;&gt; &gt; Wien mailing list<br>&gt;&gt; &gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
&gt;&gt; &gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt;&gt; &gt;<br>&gt;&gt; &gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; --<br>
&gt;&gt; Laurence Marks<br>&gt;&gt; Department of Materials Science and Engineering<br>&gt;&gt; MSE Rm 2036 Cook Hall<br>&gt;&gt; 2220 N Campus Drive<br>&gt;&gt; Northwestern University<br>&gt;&gt; Evanston, IL 60208, USA<br>
&gt;&gt; Tel: (847) 491-3996 Fax: (847) 491-7820<br>&gt;&gt; email: L-marks at northwestern dot edu<br>&gt;&gt; Web: <a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a><br>&gt;&gt; Research is to see what everybody else has seen, and to think what<br>
&gt;&gt; nobody else has thought<br>&gt;&gt; Albert Szent-Gyorgi<br>&gt;&gt; _______________________________________________<br>&gt;&gt; Wien mailing list<br>&gt;&gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
&gt;&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt;&gt;<br>&gt; <br>&gt; <br>&gt; <br>&gt; -- <br>&gt; Shamik Chakrabarti<br>&gt; Senior Research Fellow<br>
&gt; Dept. of Physics &amp; Meteorology<br>&gt; Material Processing &amp; Solid State Ionics Lab<br>&gt; IIT Kharagpur<br>&gt; Kharagpur 721302<br>&gt; INDIA<br></div>