<div dir="ltr">Dear Sir,<br><br>Thank you very much for your detailed reply. It is very informative and I got the results what we expected.<br>A comparison for :FER and :GAP (Ry) is below:<br><br>                                    :FER                               :GAP<br>option 15           0.5740619683                         0.0478      >>> chemical potential (mue): 0.5501619683<br>option 25 :         0.5820501038                           0.0464<br>option 26 :        0.5820700726                            0.0464<br>EX/VX_EV93 EC/VC_NONE: 0.6497801651       0.0554<br><br>The above comparison confirms that option 15 gives the results as you told (it includes EC_xx and VC_XX also) and that's why I have a clear difference between FER obtained from option 15 and EX_EV93 EC_NONE VX_EV93 VC.<br><br>Now we can do a better comparison with the provided test case:<br><br>The FER reported by Madsen is:0.55475  in case.intrance file which is nearly equal is we replace it with the value of mue obtained from option 15 that is:0.5501619683.<br><br>I ran Boltztrap using both values (FER and mue) but my case.trace/case.condtens are very different from the provided with Madsen.<div><br></div><div>I shall mention that I did not optimize the structure. I used the same structure as provided by Madsen with test case.</div><div><br></div><div>One obtained with FER is attached here : <a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/data.trace">case.trace</a>.</div><div><br></div><div>Still I am doing something wrong.</div><div><br></div><div>Warm regards</div><div>Fatima</div><div><br><br><br>On Tue, May 16, 2017 at 12:01 PM, Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>> wrote:<br>><br>> My earlier remark (probably it appears more than once or in different ways in the mailing list)<br>> means that you cannot use the Engel-Vosko (EV) based functional for optimitzation of the crystal structure.<br>><br>> If you read the paper of Engel and Vosco, you find that they optimized the exchange part of the functional.<br>> The method they used for that is the optimized potential model that is sometimes also called exact exchange method.<br>> If you check the contributions to the energy you will find that the EV contributes a much different amount of exchange energy Ex to the exchange-correlation energy Exc<br>> compared to other GGA's say Perdew-Burke-Ernzerhof (PBE).<br>><br>> Note that (nearly) all exchange-correlation functionals have an exchange and a correlation part.<br>> (There might be methods where one or the other is not needed to be calculated, but that's a different story)<br>> Exception appear if the two parts cannot be divided from each other. Examples may be the Slater X_alpha fuctional or the Kohn-Sham functional that is sometimes called an exchange-only functional.<br>> If you read  the paper of Engel and Vosco carefully, you will find that the latter, the Kohn-Sham xc functional, still enters<br>> the exchange potential and energy (v_x and E_x) of the Engel Vosko GGA as v^LDA, E^LDA.<br>> Indeed, one can use the Engel-Vosko functional as an exchange-only functional, however, in the old option 15 it was used together with the Perdew-Wang correlation part<br>> and thus the results will be different if you compare the "exchange-only" calculation to an "exchange-correlation" calculation.<br>><br>> Note: The EV93 + PW91 functional may lead to strange magnetic states in spin polarized calculations.<br>> I use it only for semiconductors  to "correct" the size of the band gap, as it is as a simple GGA much faster than other methods.<br>><br>> For more you should always read and understand the original papers on the different functionals before using them.<br>><br>> Ciao<br>> Gerhard<br>><br>> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>> "I think the problem, to be quite honest with you,<br>> is that you have never actually known what the question is."<br>><br>> ====================================<br>> Dr. Gerhard H. Fecher<br>> Institut of Inorganic and Analytical Chemistry<br>> Johannes Gutenberg - University<br>> 55099 Mainz<br>> and<br>> Max Planck Institute for Chemical Physics of Solids<br>> 01187 Dresden<br>> ________________________________________<br>> Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>] im Auftrag von fatima DFT [<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>]<br>> Gesendet: Montag, 15. Mai 2017 14:49<br>> An: A Mailing list for WIEN2k users<br>> Betreff: Re: [Wien] Coult not reproduced results of CoSb3<br>><br>> 1. For numeric options:<br>><br>> Let me do a more careful calculation and will update.<br>><br>> 2. For EX_EV93 EC_NONE VX_EV93 VC_NONE:<br>><br>> I read one of your thread<<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg03993.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg03993.html</a>> where you mentioned that:  Note that the Engel-Vosko GGA is not designed for total energies but optimized for exact exchange only.<br>><br>> Then I looked to UG and in table 7.4 and 7.6 where I got the information for EX_xx and VX_xx only. So I considered EC_xx and VC_xx should be EC_NONE and VC_NONE, respectively.<br>><br>> Regards<br>> Fatima<br>><br>> On Mon, May 15, 2017 at 5:22 PM, Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>>> wrote:<br>> I did a fast test with Si, and I have clearly different band gaps when switching between PBE (option 13) and  EV exchange parts of the GGA with different Correlation parts (15, 25, 26, or the corresponding EX, EC, VX, VC switches)<br>><br>> It seems to me that you do not realy know what you are doing when you are switching to different exchange-correlation functionals.<br>><br>> Who told you that EX_EV93 EC_NONE VX_EV93 VC_NONE  is supposed for EV approach ?<br>><br>> Georg Madsen pointed in his paper only on (E. Engel, S.H. Vosko, Phys. Rev. B 47 (1993) 13164–13174.) but there is no hint on the correlation part,<br>> from the cited Engel Vosko paper one expects that this was taken from the Perdew Wang 1991 GGA, however, there are also hints on GEA.<br>> Option 15 corresponds to the Engel Vosko exchange part combined with Perdew Wang 1991 correlation part<br>>           25 corresponds to the Engel Vosko exchange part combined with the LSDA correlation term<br>>           26 to a combination of 15 and 25 (roughly speaking).<br>><br>> PS.: I would wonder if for CoSb3 all different functionals would give the same result as PBE.<br>><br>> Ciao<br>> Gerhard<br>><br>> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>> "I think the problem, to be quite honest with you,<br>> is that you have never actually known what the question is."<br>><br>> ====================================<br>> Dr. Gerhard H. Fecher<br>> Institut of Inorganic and Analytical Chemistry<br>> Johannes Gutenberg - University<br>> 55099 Mainz<br>> and<br>> Max Planck Institute for Chemical Physics of Solids<br>> 01187 Dresden<br>> ________________________________________<br>> Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>>] im Auftrag von fatima DFT [<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>>]<br>> Gesendet: Montag, 15. Mai 2017 13:04<br>> An: A Mailing list for WIEN2k users<br>> Betreff: Re: [Wien] Coult not reproduced results of CoSb3<br>><br>> Sir, See my pointwise answer:<br>><br>> Did you save case.in0 after you changed to option 15 ?<br>><br>> Yes, I saved after case.in0 after replacing PBE with 15.<br>><br>> and did you NOT run again the initialisation (for PBE) after you changed to option 15 in case.in0 ?<br>><br>> I did not run init_lapw again after replacing PBE with 15.<br>><br>> If you use the new switches EX_... etc, and the results differ from earlier work, may be you have to set to the same correlation energy and potential  (EC_xxx, VC_xxx) , that was used by Madsen<br>><br>> Sir, I do not know which switched Pof. Madsen used for EC_xx and VC_xx. It is not mentioned in the paper so simpley I used NONE. Could you please tell me what switched I should use with EC/VC_xxx?<br>><br>> Regards<br>> Fatima<br>><br>><br>> Ciao<br>> Gerhard<br>><br>> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>> "I think the problem, to be quite honest with you,<br>> is that you have never actually known what the question is."<br>><br>> ====================================<br>> Dr. Gerhard H. Fecher<br>> Institut of Inorganic and Analytical Chemistry<br>> Johannes Gutenberg - University<br>> 55099 Mainz<br>> and<br>> Max Planck Institute for Chemical Physics of Solids<br>> 01187 Dresden<br>> ________________________________________<br>> Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>>>] im Auftrag von fatima DFT [<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>>>]<br>> Gesendet: Montag, 15. Mai 2017 11:56<br>> An: A Mailing list for WIEN2k users<br>> Betreff: Re: [Wien] Coult not reproduced results of CoSb3<br>><br>> Sir, I did not test the options 25 and 26 but I tested option 15 which is exactly same as PBE. So, the option 15 is not the correct option.<br>>  I initialised the case with option 15. when I saw case.in0 file PBE was printed there. I replaced that PBE with "15" and calculated further properties.<br>> In another directory, I ran the case with PBE option and compared results obtained from both and both are exactly same.<br>><br>> Option EX_EV93 EC_NONE VX_EV93 VC_NONE gave me quite different results but then mu results are not in agreement with the reported in the test case of Boltztrap.<br>><br>> Please correct me where I have mistaken? My thermoelectric properties are quite different that reported in test case.<br>><br>> Warm regards<br>> Fatima<br>><br>><br>> On Mon, May 15, 2017 at 12:24 PM, Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>>><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>>>>> wrote:<br>> If I read the code correctly then options 15, 25, and 26 should still be related to the functional of Engel and Vosco<br>><br>> Ciao<br>> Gerhard<br>><br>> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>> "I think the problem, to be quite honest with you,<br>> is that you have never actually known what the question is."<br>><br>> ====================================<br>> Dr. Gerhard H. Fecher<br>> Institut of Inorganic and Analytical Chemistry<br>> Johannes Gutenberg - University<br>> 55099 Mainz<br>> and<br>> Max Planck Institute for Chemical Physics of Solids<br>> 01187 Dresden<br>> ________________________________________<br>> Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>>><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>>>>] im Auftrag von fatima DFT [<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>>><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a><mailto:<a href="mailto:fatimadft5@gmail.com">fatimadft5@gmail.com</a>>>>]<br>> Gesendet: Samstag, 13. Mai 2017 23:18<br>> An: A Mailing list for WIEN2k users<br>> Betreff: [Wien] Coult not reproduced results of CoSb3<br>><br>> Dear Sir<br>><br>> I tried to reproduce the case.trace and case.condtens provide with Boltztrap test case.<br>> I took CoSb3.struct<<a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CoSb3_ex.struct">https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CoSb3_ex.struct</a>> file and run it with  EX_EV93 EC_NONE VX_EV93 VC_NONE (which is supposed for EV approach), RKMax=5.5, numk 300 shifted, non-sp,  -ecut -6 Gmax=12,   and then after scf I took  24000 non-shifted mesh after as mentioned in the calculation part of the original paper of Boltztrap: Madsen et al. Commun. 2006 under section: 3.1.2. Test case: CoSb3.<br>><br>> scf was run using:  run_lapw -p -i 300 -ec 0.00001 -cc 0.0001<br>><br>> I see in the case of CoSb3 the authors reported FER: 0.55475(Ry) for a thermoelectric part but here I am getting 0.6497801651 with a gap 0.0554 Ry.<br>> .<br>> A complete scf run is here: CoSb3_15.tar.gz<<a href="https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CoSb3_15.tar.gz">https://sites.google.com/site/dftfatima5/letter/mailing-list-queries/CoSb3_15.tar.gz</a>>,<br>><br>><br>> I also tested scf with option 15 (EV-GGA as suggested by Prof Fecher Fecher<<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg04252.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg04252.html</a>> and Dr. <<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg03990.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg03990.html</a>> Tran<<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg03990.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg03990.html</a>>  got FER:  0.4971964326,  GAP 0.0443(Ry). If I run init_lapw for option 15 then it runs for PBE-GGA while prof Fecher and Dr. Tran suggested option 15 for EV-GGA.<br>><br>> Could you please correct me where I did a mistake?<br>><br>> Warm regards<br>><br>> Fatima<br>> _______________________________________________<br>> Wien mailing list<br>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>>>><br>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>>><br>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a></div></div>