Re: [Wien] effective mass
By fitting band structure data (vbm and CBM). On Sunday, 30 October 2016, boudiaf khadidja <boudiafkhadidj...@gmail.com> wrote: > how can i do the effective mass "m* e m*h" calculation?, plz. > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] (no subject)
here is optimized SnO2 structure blebleble P LATTICE,NONEQUIV.ATOMS: 2 136_P42/mnm MODE OF CALC=RELA unit=bohr 8.990462 8.990462 6.039906 90.00 90.00 90.00 ATOM -1: X=0. Y=0. Z=0. MULT= 2 ISPLIT= 8 -1: X=0.5000 Y=0.5000 Z=0.5000 Sn4+ NPT= 781 R0=0.1000 RMT= 1.99Z: 50.0 LOCAL ROT MATRIX:0.7071068 0.7071068 0.000 -0.7071068 0.7071068 0.000 0.000 0.000 1.000 ATOM -2: X=0.3071 Y=0.3071 Z=0. MULT= 4 ISPLIT= 8 -2: X=0.6929 Y=0.6929 Z=0. -2: X=0.1929 Y=0.8071 Z=0.5000 -2: X=0.8071 Y=0.1929 Z=0.5000 O 2- NPT= 781 R0=0.0001 RMT= 1.72Z: 8.0 LOCAL ROT MATRIX:0.000-0.7071068 0.7071068 0.000 0.7071068 0.7071068 -1.000 0.000 0.000 16 NUMBER OF SYMMETRY OPERATIONS -1 0 0 0. 0-1 0 0. 0 0-1 0. 1 -1 0 0 0. 0-1 0 0. 0 0 1 0. 2 0-1 0 0. -1 0 0 0. 0 0-1 0. 3 0-1 0 0. -1 0 0 0. 0 0 1 0. 4 0 1 0 0. 1 0 0 0. 0 0-1 0. 5 0 1 0 0. 1 0 0 0. 0 0 1 0. 6 1 0 0 0. 0 1 0 0. 0 0-1 0. 7 1 0 0 0. 0 1 0 0. 0 0 1 0. 8 0 1 0 0.5000 -1 0 0 0.5000 0 0 1 0.5000 9 0-1 0 0.5000 1 0 0 0.5000 0 0-1 0.5000 10 0-1 0 0.5000 1 0 0 0.5000 0 0 1 0.5000 11 -1 0 0 0.5000 0 1 0 0.5000 0 0-1 0.5000 12 -1 0 0 0.5000 0 1 0 0.5000 0 0 1 0.5000 13 1 0 0 0.5000 0-1 0 0.5000 0 0-1 0.5000 14 1 0 0 0.5000 0-1 0 0.5000 0 0 1 0.5000 15 0 1 0 0.5000 -1 0 0 0.5000 0 0-1 0.5000 16 On Fri, Mar 4, 2016 at 1:22 PM, Qasim Mahmood <qasimmahmo...@gmail.com> wrote: > Dear Wien2k users I want to optimize the the tetragonal structure of SnO2 > with space group 136_P42/mnm, the error of OPT_vol_-5.0 > > Stop error will occure at the start, please help me, What should i do to > overcome this error > > Thanks and and regards > qasim > > > > > > > > > > > > > > > > > > *Mr.Qasim Mahmood* > *Ph.D Schollar, PU,Lahore,Pakistan* > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] case.outputeos
but here is the error $ x eosfit ERROR IN OPENING UNIT: 55 FILENAME: opt.vol STATUS: old FORM:formatted OPEN FAILED 0.0u 0.0s 0:00.00 0.0% 0+0k 0+0io 0pf+0w On Sat, Feb 6, 2016 at 1:57 PM, Stefaan Cottenier < stefaan.cotten...@ugent.be> wrote: > on the command line: > > x eosfit > > > Muhammad Sajjad <sajja...@gmail.com> schreef op 6 februari 2016 11:37:14 > CET: >> >> Dear All >> very basic question. after running optimize.job I did not find >> case.outputeos file. Do I need some more command to generate it? Actually >> in w2web no issue but have no idea for remote system. >> >> Kind Regards >> M. Sajjad >> >> -- >> >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> >> > -- > Sent from a tiny keyboard. Please excuse my brevity. > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] case.outputeos
Dear All very basic question. after running optimize.job I did not find case.outputeos file. Do I need some more command to generate it? Actually in w2web no issue but have no idea for remote system. Kind Regards M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] QTLB error in SOC+U calculation
Yes I did save lapw Use initso Wien2k 14.2 On Thursday, 14 January 2016, Gavin Abo <gs...@crimson.ua.edu> wrote: > Did you do a save_lapw between your steps 1 and 2 [ > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg11363.html > ]? > > Step 2, you prepared all your files manually? There is script initso_lapw > that you can use for that [ > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg10625.html > ]. > > Are you using WIEN2k 14.2 with the clmchange patch [ > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg12595.html > ]? > > You checked the error and scf files according to the FAQ and made sure > that you do not have any QTL-B warnings or ghostbands [ > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg06051.html > , > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg13577.html > ]? > > On 1/14/2016 12:43 AM, Muhammad Sajjad wrote: > > What I did is > 1. run simple scf (spin polarized) > 2. prepared inputs for SOC and U. > 3. runsp_lapw -so -orb -cc 0.0001 -ec 0.0001 > The required files are attached herewith. > In structure I have Te3BiPr > > > On Thu, Jan 14, 2016 at 11:38 AM, Lyudmila Dobysheva < > <javascript:_e(%7B%7D,'cvml','lyuk...@mail.ru');>lyuk...@mail.ru > <javascript:_e(%7B%7D,'cvml','lyuk...@mail.ru');>> wrote: > >> The program came to a wrong way. Why? Nobody can guess without >> information. Start from sending files case.inorb, indm, struct; maybe some >> other files also have to be sent... >> Did you used soc in "SImple scf" that runs well? If not - send case.inso. >> >> Best wishes >> Lyudmila Dobysheva >> -- >> Phys.-Techn. Institute of Ural Br. of Russian Ac. of Sci. >> 426001 Izhevsk, ul.Kirova 132 >> RUSSIA >> -- >> Tel.:7(3412) 432045(office), 722529(Fax) >> E-mail: l...@ftiudm.ru <javascript:_e(%7B%7D,'cvml','l...@ftiudm.ru');>, >> lyuk...@mail.ru <javascript:_e(%7B%7D,'cvml','lyuk...@mail.ru');> >> (office) >> lyuk...@gmail.com >> <javascript:_e(%7B%7D,'cvml','lyuk...@gmail.com');> (home) >> Skype: lyuka17 (home), lyuka18 (office) >> http://ftiudm.ru/content/view/25/103/lang,english/ >> -- >> >> 14.01.2016 00:51, Muhammad Sajjad wrote: >> >>> HOw to solve the issue of appearing QTLB erro in soc+U calculation? >>> SImple scf runs well but this appear in soc+U in cycle 2 as: >>> hup: Command not found. >>> LAPW0 END >>> LAPW1 END >>> LAPW1 END >>> LAPWSO END >>> LAPW2 END >>> LAPW2 END >>> LAPWDM END >>> CORE END >>> CORE END >>> MIXER END >>> in cycle 2ETEST: 7.15292635 CTEST: 7.6432934 >>> hup: Command not found. >>> LAPW0 END >>> ORB END >>> ORB END >>> LAPW1 END >>> LAPW1 END >>> LAPWSO END >>> L2main - QTL-B Error >>> >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at >> <javascript:_e(%7B%7D,'cvml','Wien@zeus.theochem.tuwien.ac.at');> >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] QTLB error in SOC+U calculation
Dear Prof. Balah I followed the instructions: 1. run normal scf. It was spin-polarized case and I found different scf.energyup and scf.energydn. 2. initso 3. runsp -so It is done but when I checked the difference of scf.energysoup and scf.energysodn there is no difference. May I know the issue? Kind regards On Thu, Jan 14, 2016 at 1:01 PM, Peter Blaha <pbl...@theochem.tuwien.ac.at> wrote: > >> SImple scf runs well but this appear in soc+U in cycle 2 as: > >> hup: Command not found. > >> LAPW0 END > >> LAPW1 END > >> LAPW1 END > >> LAPWSO END > >> LAPW2 END > >> LAPW2 END > >> LAPWDM END > >> CORE END > >> CORE END > >> MIXER END > >> in cycle 2ETEST: 7.15292635 CTEST: 7.6432934 > > > From these large values of ETEST and CTEST I can see that you either did > not start with the previous scf-solution, or that you still have the bug in > symmetso included (discussed previously in the mailing list), which > produces a bad density when the symmetry is reduced due to SO. > > I'd try a runsp -so (but without -orb) first. If it does not help, you > have to restart from the beginning using: > > instgen (creation of a new case.inst file, because the number of atoms > has changed) > x lstart > x dstart (-up/-dn) > runsp -so > > > > On 01/14/2016 09:42 AM, Muhammad Sajjad wrote: > >> Yes I did save lapw >> Use initso >> Wien2k 14.2 >> >> >> On Thursday, 14 January 2016, Gavin Abo <gs...@crimson.ua.edu >> <mailto:gs...@crimson.ua.edu>> wrote: >> >> Did you do a save_lapw between your steps 1 and 2 [ >> >> http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg11363.html >> ]? >> >> Step 2, you prepared all your files manually? There is script >> initso_lapw that you can use for that [ >> >> http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg10625.html >> ]. >> >> Are you using WIEN2k 14.2 with the clmchange patch [ >> >> http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg12595.html >> ]? >> >> You checked the error and scf files according to the FAQ and made >> sure that you do not have any QTL-B warnings or ghostbands [ >> >> http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg06051.html >> , >> >> http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg13577.html >> ]? >> >> On 1/14/2016 12:43 AM, Muhammad Sajjad wrote: >> >>> What I did is >>> 1. run simple scf (spin polarized) >>> 2. prepared inputs for SOC and U. >>> 3. runsp_lapw -so -orb -cc 0.0001 -ec 0.0001 >>> The required files are attached herewith. >>> In structure I have Te3BiPr >>> >>> >>> On Thu, Jan 14, 2016 at 11:38 AM, Lyudmila Dobysheva >>> <<javascript:_e(%7B%7D,'cvml','lyuk...@mail.ru');>lyuk...@mail.ru >>> <javascript:_e(%7B%7D,'cvml','lyuk...@mail.ru');>> wrote: >>> >>> The program came to a wrong way. Why? Nobody can guess without >>> information. Start from sending files case.inorb, indm, >>> struct; maybe some other files also have to be sent... >>> Did you used soc in "SImple scf" that runs well? If not - send >>> case.inso. >>> >>> Best wishes >>> Lyudmila Dobysheva >>> >>> -- >>> Phys.-Techn. Institute of Ural Br. of Russian Ac. of Sci. >>> 426001 Izhevsk, ul.Kirova 132 >>> RUSSIA >>> >>> -- >>> Tel.:7(3412) 432045(office), 722529(Fax) >>> E-mail: l...@ftiudm.ru >>> <javascript:_e(%7B%7D,'cvml','l...@ftiudm.ru');>, >>> lyuk...@mail.ru >>> <javascript:_e(%7B%7D,'cvml','lyuk...@mail.ru');> (office) >>> lyuk...@gmail.com >>> <javascript:_e(%7B%7D,'cvml','lyuk...@gmail.com');> (home) >>> >>> Skype: lyuka17 (home), lyuka18 (office) >>> http://ftiudm.ru/content/view/25/103/lang,english/ >>> >>> ---
[Wien] QTLB error in SOC+U calculation
Dear All HOw to solve the issue of appearing QTLB erro in soc+U calculation? SImple scf runs well but this appear in soc+U in cycle 2 as: hup: Command not found. LAPW0 END LAPW1 END LAPW1 END LAPWSO END LAPW2 END LAPW2 END LAPWDM END CORE END CORE END MIXER END in cycle 2ETEST: 7.15292635 CTEST: 7.6432934 hup: Command not found. LAPW0 END ORB END ORB END LAPW1 END LAPW1 END LAPWSO END L2main - QTL-B Error > stop error The file uplapw2.error says "'l2main' - QTL-B.GT.15., Ghostbands, check scf files" -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] QTLB error in SOC+U calculation
What I did is 1. run simple scf (spin polarized) 2. prepared inputs for SOC and U. 3. runsp_lapw -so -orb -cc 0.0001 -ec 0.0001 The required files are attached herewith. In structure I have Te3BiPr On Thu, Jan 14, 2016 at 11:38 AM, Lyudmila Dobysheva <lyuk...@mail.ru> wrote: > The program came to a wrong way. Why? Nobody can guess without > information. Start from sending files case.inorb, indm, struct; maybe some > other files also have to be sent... > Did you used soc in "SImple scf" that runs well? If not - send case.inso. > > Best wishes > Lyudmila Dobysheva > -- > Phys.-Techn. Institute of Ural Br. of Russian Ac. of Sci. > 426001 Izhevsk, ul.Kirova 132 > RUSSIA > -- > Tel.:7(3412) 432045(office), 722529(Fax) > E-mail: l...@ftiudm.ru, lyuk...@mail.ru (office) > lyuk...@gmail.com (home) > Skype: lyuka17 (home), lyuka18 (office) > http://ftiudm.ru/content/view/25/103/lang,english/ > ------ > > 14.01.2016 00:51, Muhammad Sajjad wrote: > >> HOw to solve the issue of appearing QTLB erro in soc+U calculation? >> SImple scf runs well but this appear in soc+U in cycle 2 as: >> hup: Command not found. >> LAPW0 END >> LAPW1 END >> LAPW1 END >> LAPWSO END >> LAPW2 END >> LAPW2 END >> LAPWDM END >> CORE END >> CORE END >> MIXER END >> in cycle 2ETEST: 7.15292635 CTEST: 7.6432934 >> hup: Command not found. >> LAPW0 END >> ORB END >> ORB END >> LAPW1 END >> LAPW1 END >> LAPWSO END >> L2main - QTL-B Error >> > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. scf.indmc Description: Binary data scf.inorb Description: Binary data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Tetra-elastt error
Thank you Prof. Blaha I followed the steps: 1. made case.struct 2. ./command_init_lapw 3. ./T_command_run_lapw USUALLY, TETRAGONAL SYMMETRY HAS INTERNAL COORDINATE AND YOU MUST RELAX YOUR COMPOUND. For calculation C11+C12 , C33, and Czz symmetry of tetragonal compound does not change. For finish type q COMMAND RUN : run_lapw -ec 0.0001 -in1new 2 -i 60 -p COMMAND RUN : min_lapw -j"run_lapw -I -fc 1.0 -i 60 -p" COMMAND RUN : q Would you like to use COMMAND RUN that you define for C11+C12,...for C11-C12, C44, and C66? Y 1. ./T_set_elast_lapw 2. cd elast-constant/c11+c12/test/ 3. ./auto_init_lapw 4. /TETRA.job This does not work for c11-c12. Even I found the file new_super.clmsum in the folder. Please suggest the possible solution. On Sun, Jan 3, 2016 at 4:26 PM, Peter Blaha <pbl...@theochem.tuwien.ac.at> wrote: > I am running tetra-elast for elastic constant calculation (for the given >> example PtSi). C11+c12, C33, Czz run quite well but for C11-C12, C44, >> C66 i face the following error. Please suggest the possible solution. >> > > DSTART - Error >> ** dstart crashed! >> cat: No match. >> 0.0u 0.0s 0:00.03 0.0% 0+0k 0+48io 0pf+0w >> error: command /opt/share/WIEN2K/v14.2-ifort-11.1.075.app/dstartpara >> -c dstart.def failed >> forrtl: severe (24): end-of-file during read, unit 14, file >> >> /home/sajjadm/sajjad/codes/tetra-elastic/example/ptsi/test/elast-constant/c66/test/new_super.clmsum >> > > This file should be present when clmextrapol is run. > > It is created by x dstart -super > > In case that dstart will crash because a case.rsp file is missing, you > have to create it using > > x lstart > > (and if lstart crashes, you need case.inst (created by instgen_lapw) > > -- > -- > Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna > Phone: +43-1-58801-165300 FAX: +43-1-58801-165982 > Email: bl...@theochem.tuwien.ac.atWIEN2k: http://www.wien2k.at > WWW: http://www.imc.tuwien.ac.at/staff/tc_group_e.php > -- > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Tetra-elastt error
Dear Users I am running tetra-elast for elastic constant calculation (for the given example PtSi). C11+c12, C33, Czz run quite well but for C11-C12, C44, C66 i face the following error. Please suggest the possible solution. .machine0 : processors running dstart in single mode DSTART - Error ** dstart crashed! cat: No match. 0.0u 0.0s 0:00.03 0.0% 0+0k 0+48io 0pf+0w error: command /opt/share/WIEN2K/v14.2-ifort-11.1.075.app/dstartpara -c dstart.def failed forrtl: severe (24): end-of-file during read, unit 14, file /home/sajjadm/sajjad/codes/tetra-elastic/example/ptsi/test/elast-constant/c66/test/new_super.clmsum Image PCRoutineLineSource clmaddsub 0048628D Unknown Unknown Unknown clmaddsub 00484D95 Unknown Unknown Unknown clmaddsub 0043FCF0 Unknown Unknown Unknown clmaddsub 0040DE6A Unknown Unknown Unknown clmaddsub 0040D660 Unknown Unknown Unknown clmaddsub 00423BAC Unknown Unknown Unknown clmaddsub 004045E7 MAIN__235 clmaddsub.f clmaddsub 00402E7C Unknown Unknown Unknown libc.so.6 2B30599AB76D Unknown Unknown Unknown clmaddsub 00402D79 Unknown Unknown Unknown 0.0u 0.0s 0:00.02 50.0% 0+0k 0+440io 0pf+0w error: command /opt/share/WIEN2K/v14.2-ifort-11.1.075.app/clmaddsub clmaddsub.def failed clmextrapol_lapw has generated a new test.clmsum -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] ELASTIC CONSTANTS OF CUBIC COMPOUND
I have sent you email with step by step guide already got from Gavin Abo. As the file size is big so can not go through this email. On Sun, Dec 20, 2015 at 8:01 AM, Nazia Erum <eru...@hotmail.com> wrote: > Hello, > Prof. Blaha and wien2k user, > > I want to calculate cubic constants of a cubic perovskite compound with > method developed by charpin. I put valid case.struct and case.inst file in > the directory but when I want to init_elast through terminal it prompts > that you don't have valid case.struc file. Then I past this init_elast in > commandline of wien2k that says xnn invalid xnn error.while I am giving > correct bond length factor which is 2. > please guide me in this respect. > > ph.d. student > > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Elastic constant calculation for tetragonal system
Thank you Tomas. I read this paper. My concern is for cubic I used instgen_lapw init_elast elast_setup ./eos.job ./rhomb.job ./tetra.job What should be the steps for tetragonal system? On Thu, Dec 17, 2015 at 1:29 PM, Tomas Kana <k...@seznam.cz> wrote: > Dear Muhammad, try this link: > > http://www.electrochemsci.org/papers/vol8/81112252.pdf > > Best regards > > Tomas > > > ------ Původní zpráva -- > Od: Muhammad Sajjad <sajja...@gmail.com> > Komu: wien <wien@zeus.theochem.tuwien.ac.at> > Datum: 16. 12. 2015 15:49:50 > Předmět: [Wien] Elastic constant calculation for tetragonal system > > Dear all > Can we use same method to calculate elastic constant for tetragonal system > as we do for cubic system? > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Elastic constant calculation for tetragonal system
Please ignore previous email, I found the way. On Thu, Dec 17, 2015 at 2:15 PM, Muhammad Sajjad <sajja...@gmail.com> wrote: > Thank you Tomas. I read this paper. My concern is > for cubic I used > instgen_lapw > init_elast > elast_setup > ./eos.job > ./rhomb.job > ./tetra.job > What should be the steps for tetragonal system? > > On Thu, Dec 17, 2015 at 1:29 PM, Tomas Kana <k...@seznam.cz> wrote: > >> Dear Muhammad, try this link: >> >> http://www.electrochemsci.org/papers/vol8/81112252.pdf >> >> Best regards >> >> Tomas >> >> >> -- Původní zpráva -- >> Od: Muhammad Sajjad <sajja...@gmail.com> >> Komu: wien <wien@zeus.theochem.tuwien.ac.at> >> Datum: 16. 12. 2015 15:49:50 >> Předmět: [Wien] Elastic constant calculation for tetragonal system >> >> Dear all >> Can we use same method to calculate elastic constant for tetragonal >> system as we do for cubic system? >> >> >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> >> >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> >> > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Elastic constant calculation for tetragonal system
Dear all Can we use same method to calculate elastic constant for tetragonal system as we do for cubic system? -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Correct band gap value
It looks same (in shape). I selected lattice type F and applied strain. Only symmetry operations reduce to 16 from 48 (for bulk Si). Also I kept the volume constant for every value of strain. On Thu, Dec 10, 2015 at 11:29 AM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > Did you check how the Brillouin zone looks like in the tetragonal case > after you applied the strain ? > Is your k-path that you used for the xmgrace plot the one that you like to > have to find the correct gap ? > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Donnerstag, 10. Dezember 2015 09:19 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] Correct band gap value > > Manually means I opened the band structure in xmgrace then simply clicked > on the CBM and the went to the text editor to read the exact lowest value. > I am simply using Silicon (with 2% strain) and the plot is along W L G X W. > > On Thu, Dec 10, 2015 at 10:35 AM, <t...@theochem.tuwien.ac.at t...@theochem.tuwien.ac.at>> wrote: > What does "manually" mean? If it means band structure plotting, then > probably the VBM and/or CBM is not along the path in the Brillouin zone > that you chose. > > F. Tran > > > On Thu, 10 Dec 2015, Muhammad Sajjad wrote: > > I selected 2000 k-points. grep GAP *.scf gives 0.923 eV while manually I > found VBM at 0.0 eV and CBM at 1.10023 eV. So which one should be correct? > > On Thu, Dec 10, 2015 at 12:20 AM, <t...@theochem.tuwien.ac.at t...@theochem.tuwien.ac.at>> wrote: > If your k-mesh is not dense enough such that there is no k-point > at the VBM and/or CBM or close to, then :GAP will be larger than > the band gap from a band structure. > > F. Tran > > On Wed, 9 Dec 2015, Muhammad Sajjad wrote: > > Dear AllDoes the command "grep GAP *.scf" gives > correct band gap value? as when I found the bottom > of conduction band in xmgrace it is higher. (bottom > of CB is > read correctly by opening the plot in xmgrace and > then read the date by editing in text editor). > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at Wien@zeus.theochem.tuwien.ac.at> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > > _______ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at<mailto:Wien@zeus.theochem.tuwien.ac.at> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Correct band gap value
Here is the structure with 2% tensile strain. SILICON F LATTICE,NONEQUIV.ATOMS: 1 MODE OF CALC=RELA unit=ang 10.511739 10.204090 10.204090 90.00 90.00 90.00 ATOM -1: X=0.1250 Y=0.1250 Z=0.1250 MULT= 2 ISPLIT=-2 -1: X=0.8750 Y=0.8750 Z=0.8750 Si NPT= 781 R0=0.0001 RMT=2.2200 Z: 14.0 LOCAL ROT MATRIX:0.000 0.000 1.000 1.000 0.000 0.000 0.000 1.000 0.000 16 NUMBER OF SYMMETRY OPERATIONS On Thu, Dec 10, 2015 at 12:05 PM, Muhammad Sajjad <sajja...@gmail.com> wrote: > It looks same (in shape). I selected lattice type F and applied strain. > Only symmetry operations reduce to 16 from 48 (for bulk Si). Also I kept > the volume constant for every value of strain. > > On Thu, Dec 10, 2015 at 11:29 AM, Fecher, Gerhard <fec...@uni-mainz.de> > wrote: > >> Did you check how the Brillouin zone looks like in the tetragonal case >> after you applied the strain ? >> Is your k-path that you used for the xmgrace plot the one that you like >> to have to find the correct gap ? >> >> Ciao >> Gerhard >> >> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: >> "I think the problem, to be quite honest with you, >> is that you have never actually known what the question is." >> >> >> Dr. Gerhard H. Fecher >> Institut of Inorganic and Analytical Chemistry >> Johannes Gutenberg - University >> 55099 Mainz >> and >> Max Planck Institute for Chemical Physics of Solids >> 01187 Dresden >> >> Von: wien-boun...@zeus.theochem.tuwien.ac.at [ >> wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ >> sajja...@gmail.com] >> Gesendet: Donnerstag, 10. Dezember 2015 09:19 >> An: A Mailing list for WIEN2k users >> Betreff: Re: [Wien] Correct band gap value >> >> Manually means I opened the band structure in xmgrace then simply clicked >> on the CBM and the went to the text editor to read the exact lowest value. >> I am simply using Silicon (with 2% strain) and the plot is along W L G X W. >> >> On Thu, Dec 10, 2015 at 10:35 AM, <t...@theochem.tuwien.ac.at> t...@theochem.tuwien.ac.at>> wrote: >> What does "manually" mean? If it means band structure plotting, then >> probably the VBM and/or CBM is not along the path in the Brillouin zone >> that you chose. >> >> F. Tran >> >> >> On Thu, 10 Dec 2015, Muhammad Sajjad wrote: >> >> I selected 2000 k-points. grep GAP *.scf gives 0.923 eV while manually I >> found VBM at 0.0 eV and CBM at 1.10023 eV. So which one should be correct? >> >> On Thu, Dec 10, 2015 at 12:20 AM, <t...@theochem.tuwien.ac.at> t...@theochem.tuwien.ac.at>> wrote: >> If your k-mesh is not dense enough such that there is no k-point >> at the VBM and/or CBM or close to, then :GAP will be larger than >> the band gap from a band structure. >> >> F. Tran >> >> On Wed, 9 Dec 2015, Muhammad Sajjad wrote: >> >> Dear AllDoes the command "grep GAP *.scf" gives >> correct band gap value? as when I found the bottom >> of conduction band in xmgrace it is higher. (bottom >> of CB is >> read correctly by opening the plot in xmgrace and >> then read the date by editing in text editor). >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> >> >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at> Wien@zeus.theochem.tuwien.ac.at> >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> >> >> >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> >> >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at<mailto:Wien@zeus.theochem.tuwien.ac.at> >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> >> >> >> >> -- >>
Re: [Wien] Correct band gap value
Manually means I opened the band structure in xmgrace then simply clicked on the CBM and the went to the text editor to read the exact lowest value. I am simply using Silicon (with 2% strain) and the plot is along W L G X W. On Thu, Dec 10, 2015 at 10:35 AM, <t...@theochem.tuwien.ac.at> wrote: > What does "manually" mean? If it means band structure plotting, then > probably the VBM and/or CBM is not along the path in the Brillouin zone > that you chose. > > F. Tran > > > On Thu, 10 Dec 2015, Muhammad Sajjad wrote: > > I selected 2000 k-points. grep GAP *.scf gives 0.923 eV while manually I >> found VBM at 0.0 eV and CBM at 1.10023 eV. So which one should be correct? >> >> On Thu, Dec 10, 2015 at 12:20 AM, <t...@theochem.tuwien.ac.at> wrote: >> If your k-mesh is not dense enough such that there is no k-point >> at the VBM and/or CBM or close to, then :GAP will be larger than >> the band gap from a band structure. >> >> F. Tran >> >> On Wed, 9 Dec 2015, Muhammad Sajjad wrote: >> >> Dear AllDoes the command "grep GAP *.scf" gives >> correct band gap value? as when I found the bottom >> of conduction band in xmgrace it is higher. (bottom >> of CB is >> read correctly by opening the plot in xmgrace and >> then read the date by editing in text editor). >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> >> >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> >> >> >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> >> > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Stress values
Dear All Can we find stress value in wien2k? -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Correct band gap value
I selected 2000 k-points. grep GAP *.scf gives 0.923 eV while manually I found VBM at 0.0 eV and CBM at 1.10023 eV. So which one should be correct? On Thu, Dec 10, 2015 at 12:20 AM, <t...@theochem.tuwien.ac.at> wrote: > If your k-mesh is not dense enough such that there is no k-point > at the VBM and/or CBM or close to, then :GAP will be larger than > the band gap from a band structure. > > F. Tran > > On Wed, 9 Dec 2015, Muhammad Sajjad wrote: > > Dear AllDoes the command "grep GAP *.scf" gives correct band gap value? as >> when I found the bottom of conduction band in xmgrace it is higher. (bottom >> of CB is >> read correctly by opening the plot in xmgrace and then read the date by >> editing in text editor). >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> >> > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Correct band gap value
Dear All Does the command "grep GAP *.scf" gives correct band gap value? as when I found the bottom of conduction band in xmgrace it is higher. (bottom of CB is read correctly by opening the plot in xmgrace and then read the date by editing in text editor). -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Fwd: Strain
Thank you very much Dr. Martin Pieper. Yes phonon calculations speak about stability of structure but not about breaking point. On Mon, Dec 7, 2015 at 1:22 PM, pieper <pie...@ifp.tuwien.ac.at> wrote: > Hello Muhammad! > > A DFT code will not simply say *crack* at some point of strain. Wien2k > will always distribute the force you apply according to the perfect > translational symmetry of your crystal structure - in contrast to nature > where things like cracking or breaking occure at weak links - points where > this translational symmetry is broken and the bonding is weaker. > > DFT will do its very best to distribute the electrons in any structure you > throw at it and tell you (an upper bound of) the minimum energy of the > configuration. I am no expert but I understand that from a calculated > variation of total energy with lattice parameter you can estimate the > intrinsic strength of a structure within some model for the process. > > In addition Wien2k gives information where the total energy is stored, for > example the forces on atoms in your structure. Therefore, if you look at > the total energy of your strained structure you will see that it is larger > - at 5% distortion perhaps considerably so. Depending on site symmetry some > calculated internal forces may indicate where in the structure the > 'springs' you pulled are located. Furthermore the phonon code might inform > you of instabilities in your structure. > > Good luck, > > Martin Pieper > > > --- > Dr. Martin Pieper > Karl-Franzens University > Institute of Physics > Universitätsplatz 5 > A-8010 Graz > Austria > Tel.: +43-(0)316-380-8564 > > > > Am 04.12.2015 15:25, schrieb Muhammad Sajjad: > >> Dear Tomas and Fecher >> Thank you very much. It is working perfect with lattice type F and >> symmetry operations equal to 16. But one thing is confusing me that I >> applied 5 % strain to bulk Si and still it is semiconductor with gap >> 0.2eV. Is it correct? as 5 % is too much strain and bulk materials >> tolerate very low strain. In my thinking it should break. >> >> On Fri, Dec 4, 2015 at 5:14 PM, Fecher, Gerhard <fec...@uni-mainz.de> >> wrote: >> >> but your structure becomes tetragonal >>> you should reset the symmetry operations to (generate) to be >>> recalculated during initialisation >>> >>> Ciao >>> Gerhard >>> >>> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: >>> "I think the problem, to be quite honest with you, >>> is that you have never actually known what the question is." >>> >>> >>> Dr. Gerhard H. Fecher >>> Institut of Inorganic and Analytical Chemistry >>> Johannes Gutenberg - University >>> 55099 Mainz >>> and >>> Max Planck Institute for Chemical Physics of Solids >>> 01187 Dresden >>> >>> Von: wien-boun...@zeus.theochem.tuwien.ac.at >>> [wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad >>> Sajjad [sajja...@gmail.com] >>> Gesendet: Freitag, 4. Dezember 2015 15:08 >>> An: A Mailing list for WIEN2k users >>> Betreff: Re: [Wien] Fwd: Strain >>> >>> I am sorry having no idea about it. Simply I used optimized >>> structure and increased a by a0*1.001 and decreased b & c by keeping >>> lattice type F. Then run init_lapw -b. >>> >>> On Fri, Dec 4, 2015 at 4:51 PM, Fecher, Gerhard >>> <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: >>> you applied a tetragonal distortion along x >>> but how did you manage that you still have 48 symmetry operations ? >>> >>> Ciao >>> Gerhard >>> >>> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: >>> "I think the problem, to be quite honest with you, >>> is that you have never actually known what the question is." >>> >>> >>> Dr. Gerhard H. Fecher >>> Institut of Inorganic and Analytical Chemistry >>> Johannes Gutenberg - University >>> 55099 Mainz >>> and >>> Max Planck Institute for Chemical Physics of Solids >>> 01187 Dresden >>> >>> Von: >>> >>> wien-boun...@zeus.theochem.tuwien.ac.at> wien-boun...@zeus.theochem.tuwien.ac.at> >> >>> >>> [wien-boun...@zeus.theochem.tuwien.ac.at> wien-boun...@zeus.theochem.tuwien.ac.at>] >> >>> im A
Re: [Wien] Fwd: Strain
Dear Tomas and Fecher Thank you very much. It is working perfect with lattice type F and symmetry operations equal to 16. But one thing is confusing me that I applied 5 % strain to bulk Si and still it is semiconductor with gap 0.2eV. Is it correct? as 5 % is too much strain and bulk materials tolerate very low strain. In my thinking it should break. On Fri, Dec 4, 2015 at 5:14 PM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > but your structure becomes tetragonal > you should reset the symmetry operations to (generate) to be recalculated > during initialisation > > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Freitag, 4. Dezember 2015 15:08 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] Fwd: Strain > > I am sorry having no idea about it. Simply I used optimized structure and > increased a by a0*1.001 and decreased b & c by keeping lattice type F. Then > run init_lapw -b. > > On Fri, Dec 4, 2015 at 4:51 PM, Fecher, Gerhard <fec...@uni-mainz.de > <mailto:fec...@uni-mainz.de>> wrote: > you applied a tetragonal distortion along x > but how did you manage that you still have 48 symmetry operations ? > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > ____ > Von: wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at> [ > wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at>] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com<mailto:sajja...@gmail.com>] > Gesendet: Freitag, 4. Dezember 2015 14:04 > An: wien > Betreff: [Wien] Fwd: Strain > > Dear All > I got Si structure from wien2k examples and after optimization it is > SILICON > F LATTICE,NONEQUIV.ATOMS: 1 > MODE OF CALC=RELA unit=ang > 10.305626 10.305626 10.305626 90.00 90.00 90.00 > ATOM 1: X=0.1250 Y=0.1250 Z=0.1250 > MULT= 2 ISPLIT= 2 >1: X=0.8750 Y=0.8750 Z=0.8750 > Si NPT= 781 R0=0.0001 RMT= 2.22Z: 14.0 > LOCAL ROT MATRIX:1.000 0.000 0.000 > 0.000 1.000 0.000 > 0.000 0.000 1.000 > 48 NUMBER OF SYMMETRY OPERATIONS > > mBJLDA gives band gap of 1.19 eV (perfect). > > Now I apply strain (0.1 %) and structure is > > SILICON > F LATTICE,NONEQUIV.ATOMS: 1 > MODE OF CALC=RELA unit=ang > 10.315932 10.300477 10.300477 90.00 90.00 90.00 > ATOM -1: X=0.1250 Y=0.1250 Z=0.1250 > MULT= 2 ISPLIT=-2 > -1: X=0.8750 Y=0.8750 Z=0.8750 > Si NPT= 781 R0=0.0001 RMT= 2.22Z: 14.0 > LOCAL ROT MATRIX:0.000 0.000 1.000 > 1.000 0.000 0.000 > 0.000 1.000 0.000 > 48 NUMBER OF SYMMETRY OPERATIONS > This time band structure (attached herewith) is very different and Si is > metal. > Kindly guide me whether I am apply strain in right way or not? Band > Structure is > > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at<mailto:Wien@zeus.theochem.tuwien.ac.at> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.a
Re: [Wien] Fwd: Strain
I am sorry having no idea about it. Simply I used optimized structure and increased a by a0*1.001 and decreased b & c by keeping lattice type F. Then run init_lapw -b. On Fri, Dec 4, 2015 at 4:51 PM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > you applied a tetragonal distortion along x > but how did you manage that you still have 48 symmetry operations ? > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Freitag, 4. Dezember 2015 14:04 > An: wien > Betreff: [Wien] Fwd: Strain > > Dear All > I got Si structure from wien2k examples and after optimization it is > SILICON > F LATTICE,NONEQUIV.ATOMS: 1 > MODE OF CALC=RELA unit=ang > 10.305626 10.305626 10.305626 90.00 90.00 90.00 > ATOM 1: X=0.1250 Y=0.1250 Z=0.1250 > MULT= 2 ISPLIT= 2 >1: X=0.8750 Y=0.8750 Z=0.8750 > Si NPT= 781 R0=0.0001 RMT= 2.22Z: 14.0 > LOCAL ROT MATRIX:1.000 0.000 0.000 > 0.000 1.000 0.000 > 0.000 0.000 1.000 > 48 NUMBER OF SYMMETRY OPERATIONS > > mBJLDA gives band gap of 1.19 eV (perfect). > > Now I apply strain (0.1 %) and structure is > > SILICON > F LATTICE,NONEQUIV.ATOMS: 1 > MODE OF CALC=RELA unit=ang > 10.315932 10.300477 10.300477 90.00 90.00 90.00 > ATOM -1: X=0.1250 Y=0.1250 Z=0.1250 > MULT= 2 ISPLIT=-2 > -1: X=0.8750 Y=0.8750 Z=0.8750 > Si NPT= 781 R0=0.0001 RMT= 2.22Z: 14.0 > LOCAL ROT MATRIX:0.000 0.000 1.000 > 1.000 0.000 0.000 > 0.000 1.000 0.000 > 48 NUMBER OF SYMMETRY OPERATIONS > This time band structure (attached herewith) is very different and Si is > metal. > Kindly guide me whether I am apply strain in right way or not? Band > Structure is > > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] About material breaking
Dear Oleg and Marks Thank you for your quick reply. Sorry for late, I was away for two days. @Oleg: Yes I mean ultimate material strength. The sent paper is quite helping. In FIG. 2., stress-strain curve is drawn. We know applied strain but how can we measure stress (F/A). Ok A is known but What is F? Simply how to get stress value? On Thu, Nov 26, 2015 at 6:39 PM, Oleg Rubel <oru...@lakeheadu.ca> wrote: > Do you mean the ultimate material strength when you refer to "breaking > point”? > > Straightforward calculations (like this one > http://arxiv.org/pdf/1303.7259v1.pdf) are possible, but the obtained > result will be orders (!) of magnitude higher than any practical values. > The reason is that it is hard to capture mechanisms of materials plasticity > at such a small scale. > > Oleg > > > On Nov 26, 2015, at 02:35, Muhammad Sajjad <sajja...@gmail.com> wrote: > > > > Dear All > > Can we use WIEN2K fto know about breaking point of material? I know to > calculate phonon spectrum (using phonopy) for the stability of material > but not breaking. Also some information from elastic constant (like C11 = > Young Mod.). > > > > -- > > Kind Regards > > Muhammad Sajjad > > Post Doctoral Fellow > > KAUST, KSA. > > ___ > > Wien mailing list > > Wien@zeus.theochem.tuwien.ac.at > > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] About material breaking
Dear All Can we use WIEN2K fto know about breaking point of material? I know to calculate phonon spectrum (using phonopy) for the stability of material but not breaking. Also some information from elastic constant (like C11 = Young Mod.). -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Bulk vs supercell Plasma frequency
Dear Prof. P Blaha Thank you very much for yous fruitful suggestions. However, I did not change RKmax. For cubic Na only W_p_xx is enough but after making supercell (1 1 6) still I need only w_p_xx and forget about w_p_zz, though the system is tetragonal now. Is it right I understood? Second thing is it principally right that one converges SCF with small k-mesh say 30*30*5 and then use the same scf ao get OPTICAL properties with k-mesh say 90*90*15? I mean no need to converge scf with 90*90*15. Kind Regards On Thu, Oct 29, 2015 at 9:03 PM, Peter Blaha <pbl...@theochem.tuwien.ac.at> wrote: > In principle you are doing thinks now correct. > In reality, however, we use the tetrahedra-method in the integration of > the BZ to calculate a joint-density-of-states. And for a metal !!! here > there is a difference between k-meshes in the small or large cell. (it is > related to the "back-folding" problem and in the supercell there are less > tetrahedra and thus less possibilities to interpolate for the tetrahedron > method than in the small cell). > > In any case, you need to consider two things: > > a) You have to decide to which accurcy a certain number should be > calculated. Typically, a plasma frequency would be quoted with one or two > digits after the comma (i.e. 5.9 or at most 5.94 eV). > > b) Then you have to converge the numbers with respect to k-mesh (but also > Rkmax !!) > Definitely you should use a ratio of 6:1 for the supercell, but eventually > you need to increase the mesh until there are no changes within the desired > accuracy. In your case, 90x15 is almost sufficient, but not quite if you > want two digits accuracy. > > c) in a cubic system, of course it is useless to quote two tensor > components and one is sufficient. > > > > > Am 29.10.2015 um 10:35 schrieb Muhammad Sajjad: > > Dear Prof. Blaha > > I did some calculation for Bulk Na (bcc, space gp #229_Im-3m) wiht k-mesh > 30*30*30 and found the plasma frequency (in case.outputjoint) as > Plasma frequencies: > >w_p_xx w_p_zz[eV] > >5.9446 5.9446 > > Now I constructed a supercell 1*1*6 (definitely symmetry reduced and now > sp.gp is 1_P1), used k-mesh 30*30*5 and found > Plasma frequencies (much deviating from bulk values): > >w_p_xx w_p_zz[eV] > >5.76146 5.3446 > > Then used k-mesh 90*90*15 and found > Plasma frequencies (w_p_zz is deviating from bulk w_p_zz ): > >w_p_xx w_p_zz[eV] > >5.9485 5.8903 > > I have read the previous post ( > https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg09782.html) > but it is not speaking about w_p_zz. So when I plot epsilon data for > w_p_xx then it is matching with that of bulk but when I use w_p_zz then it > is definitely away from bulk. Should I consider w_p_zz or use "number > of choice = 1 in case.inop" and plot epsilon only with w_p_xx ? > > > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > > ___ > Wien mailing > listw...@zeus.theochem.tuwien.ac.athttp://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > > _______ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Bulk vs supercell Plasma frequency
Dear Prof. Blaha I did some calculation for Bulk Na (bcc, space gp #229_Im-3m) wiht k-mesh 30*30*30 and found the plasma frequency (in case.outputjoint) as Plasma frequencies: w_p_xx w_p_zz[eV] 5.9446 5.9446 Now I constructed a supercell 1*1*6 (definitely symmetry reduced and now sp.gp is 1_P1), used k-mesh 30*30*5 and found Plasma frequencies (much deviating from bulk values): w_p_xx w_p_zz[eV] 5.76146 5.3446 Then used k-mesh 90*90*15 and found Plasma frequencies (w_p_zz is deviating from bulk w_p_zz ): w_p_xx w_p_zz[eV] 5.9485 5.8903 I have read the previous post ( https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg09782.html) but it is not speaking about w_p_zz. So when I plot epsilon data for w_p_xx then it is matching with that of bulk but when I use w_p_zz then it is definitely away from bulk. Should I consider w_p_zz or use "number of choice = 1 in case.inop" and plot epsilon only with w_p_xx ? Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Convergence and optical properties
Dear Users Suppose I got a perfect convergence (for bulk) with some k-mesh say 30*30*30 and computed optical properties with k-mesh 90*90*90. Again I got a perfect convergence (for supercell 1*1*6) with k-mesh say 20*20*10 (rather than 30*30*6, just for time saving and luckily convergence is done) and then computed optical properties with k-meshes 90*90*15. My questions are: is it computationally right? Can I used the previous calculation (i.e directory with converged SCF with k-mesh say 20*20*10) to compute optical properties with different k-meshes, like simply pick previous calculation, increase Emax in case.in1, x kgen, lapw1 and so on? -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Fwd: Convergence and optical properties
-- Forwarded message -- From: Muhammad Sajjad <sajja...@gmail.com> Date: Wed, Oct 14, 2015 at 12:05 PM Subject: Convergence and optical properties To: wien <Wien@zeus.theochem.tuwien.ac.at> Dear Users Suppose I got a perfect convergence (for bulk) with some k-mesh say 30*30*30 and computed optical properties with k-mesh 90*90*90. Again I got a perfect convergence (for supercell 1*1*6) with k-mesh say 20*20*10 (rather than 30*30*6, just for time saving and luckily convergence is done) and then computed optical properties with k-meshes 90*90*15. My questions are: is it computationally right? Can I used the previous calculation (i.e directory with converged SCF with k-mesh say 20*20*10) to compute optical properties with different k-meshes, like simply pick previous calculation, increase Emax in case.in1, x kgen, lapw1 and so on? -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Different Plasma Frequency (PF) for bulk and supercell Na
Dear Lyudmila and BHamu Thank you for your suggestions. I have already used 18*18*3 for supercell and 18*18*18, but found difference in plasma frequency for supercell and bulk. Should plasma frequency be same for bulk and supercell? or I should use the new plasma frequency obtained from super cell calculation? @Lyudmila | In terms of accuracy, this is the worst. but why? On Tue, Oct 13, 2015 at 10:55 AM, Dr. K. C. Bhamu <kcbham...@gmail.com> wrote: > > > Dear sajjad > > I am also a beginner but I must suggest you to choose 18*18*3 k points for > 1*1*6 supercell in comparison with bulk. > As volume in real space is inversely proportional to volume of reciprocal > space. Larger the volume in real space lower the volume in reciprocal > space. So while computing properties of 1*1*6 structure you are minimizing > your volume in z direction and hence need less k-points. > > I appreciate any comment on this from DFT experts. > > > Sincerely > Bhamu > > I have computed PF for bulk Na and its value is *5.943 eV* with k-mesh >> *18*18*18*. >> >> Then I constructed a supercell 1*1*6 >> > > ...>>>>>>> use 18*18*3 k points for 1*1*6 supercell. > > and run the calculation with different k-meshes and obtained PF values. >> The details are >> >> *K-meshPF (x/z components in eV)* >> >> *24*24*24 5.7097 / 5.8192* >> >> >> *36*36*24 5.7614 / 5.8098 * >> >> >> *50*50*10 5.9358 / 5.6280 * >> >> *50*50*30 5.9365 / 5.8079* >> >> My question is why PF is very sensitive to k-mesh selection? and which >> value of PF is required for correct dielectric constant calculation? Do I >> need *PF = 5.943 eV* for the supercell and then, definitely, need more >> test calculation for right k-mesh selection? >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html >> >> > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Diference in EFG values
Dear Prof. Stefaan and Blaha Thank you for your reply. Actually I have to reduce separation energy to -8.5 Ry because core leakage issue was appearing with Ge. If I focus on Co, the useful information is: *Magnetic moment * VASPWIEN2K Co1 2.530:MMI001: = 2.35929 Co2 -2.530 :MMI002: = -2.35900 *Core states * VASP WIEN2K VRHFIN = Co: d8 s1 E-up(Ry) E-dn(Ry) Occupancy q/sphere core-state 1S -557.761897 -557.761629 1.00 1.00 1.T 2S -66.048914 -65.9348801.00 1.00 1.T 2P* -57.314278-57.2287991.00 1.00 1.T 2P -56.225866 -56.1379022.00 2.00 1.T 3S -7.605100 -7.363843 1.00 1.00 0.9968 F 3P* -5.003328 -4.767091 1.00 1.00 0.9917 F 3P -4.866118 -4.631354 2.00 2.00 0.9909F 3D* -0.699000 -0.490964 2.00 2.00 0.8948F 3D -0.684578 -0.477338 3.00 0.00 0.9053F 4S -0.423921 -0.381223 1.00 1.00 0.1768F If I use -6 Ry as the separation energy then core states are E-up(Ry) E-dn(Ry) Occupancy q/sphere core-state 1S -557.761897 -557.761629 1.00 1.00 1.T 2S -66.048914 -65.9348801.00 1.00 1.T 2P* -57.314278-57.2287991.00 1.00 1.T 2P -56.225866 -56.1379022.00 2.00 1.T 3S -7.605100 -7.363843 1.00 1.00 0.9968 T 3P* -5.003328 -4.767091 1.00 1.00 0.9917 F 3P -4.866118 -4.631354 2.00 2.00 0.9909F 3D* -0.699000 -0.490964 2.00 2.00 0.8948F 3D -0.684578 -0.477338 3.00 0.00 0.9053F 4S -0.423921 -0.381223 1.00 1.00 0.1768F According to my understanding, core-state with tag F are valence states. AM I right? If yes then why it is taking 3P states as valence states? The valence electrons for Co are 3d7, 4s2. Kind Regards On Mon, Sep 14, 2015 at 5:16 PM, Peter Blaha <pbl...@theochem.tuwien.ac.at> wrote: > And of course, make sure that the same DFT functionals are used (not > comparing DFT+U and DFT alone. This could explain the differences on Co. > > PS: In addition, I'd make sure that the basic "electronic structure" is > identical (magnetic moments, DOS, bands, ...) > > On 09/14/2015 04:07 PM, Stefaan Cottenier wrote: > >> >> First guess (assuming everything is numerically converged): do you >> consider the same electrons as valence electrons both in vasp and wien2k >> ? It could happen that a 'semicore' state is taken as valence in wien2k >> and core in vasp. As long as the EFG contribution of these states are >> small, wien2k and vasp will show no differences. But if you hit a >> crystal structure where these states have a larger EFG contribution, >> then the differences will show up. >> >> You can find out whether or not this is the case by (1) explicitly >> taking the same core/valence assignment in both codes, or (2) examining >> the contribution to the EFG of all different oribitals (or regions of >> the DOS). See http://www.wien2k.at/reg_user/faq/efg2.pdf (top of page 9) >> for the procedure to follow in wien2k. >> >> No warranty, just a guess... >> >> Stefaan >> >> >> Op 14/09/2015 om 15:53 schreef Muhammad Sajjad: >> >>> Dear Users >>> >>> I run some test calculations for EFG values for different compounds >>> (Sc2O3, In, SmCo5). Their EFG values (computed with VASP) are in >>> agreement with that of previous values. Then I computed the EFG v
Re: [Wien] Diference in EFG values
Right I forgot to mention about XC potentials. Identical Potentials were used. I already read the document related to computing EFG by WIEN2K. How to isolate the EFG contribution of Co-3p? On Tue, Sep 15, 2015 at 2:55 PM, Stefaan Cottenier < stefaan.cotten...@ugent.be> wrote: > According to my understanding, core-state with tag F are valence states. >> AM I right? If yes then why it is taking 3P states as valence states? >> The valence electrons for Co are 3d7, 4s2. >> > > You use the words 'valence electrons' in the traditional 'chemical' way as > states outside the noble gas core. In the context of the DFT, the term > 'valence electrons' refers to all states that are not forced to behave as > in a free atom (i.e. they allow to feel the presence of the surrounding > crystal). > > You cannot avoid to use Co-3p as valence states with wien2k, the numerics > would break down if you didn't. > > Your output shows that wien2k has Co-3s and Co-3p as valence states, > whereas vasp considers these as core states (i.e. they will not contribute > to the EFG in vasp). Co-3s will never contribute to the EFG (spherical), > whereas Co-3p could (albeit not very likely). > > First consider the suggestion by Peter Blaha: make sure you have identical > XC-functionals in both codes, and inspect whether all magnetic moments and > the DOS and band structure pictures agree for both codes. > > If that is all right, then use the procedure to which I refered in my > previous post to isolate the EFG contribution of Co-3p. If that is nearly > zero, I don't know what is going on. If it would be equal to the > vasp-wien2k difference, then you have found the origin of your problem. > > Stefaan > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] EFG dependence on magnetic field
Dear Fecher Sorry I can not understand about which energy your are asking about. My intention was to see how B effects on EFG. By the way Without B_ext TOTAL ENERGY IN Ry = -89574.89353469 with B_ext TOTAL ENERGY IN Ry = -89574.89351208 On Wed, Sep 9, 2015 at 1:14 PM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > Did you check how much the energy is changed if you apply an external > field of 8 T > the Bohr magneton is approximately 58 x 10^-6 eV / T > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Montag, 7. September 2015 09:55 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] EFG dependence on magnetic field > > Dear Gavin > Bundle of thanks for such a helping reply. Are the case.inorb and > case.indm files are ok now for B-ext applied at 32 degree angle with x-axis? > > > == case.inorb = > 3 2 0 nmod, natorb, ipr > PRATT 1.0BROYD/PRATT, mixing > 1 1 2 iatom nlorb, lorb > > 2 1 2 iatom nlorb, lorb > > 8. Bext 1. 0.62487 0. direction > == == case.indm > = -9. Emin cutoff energy 2 number of atoms for which > density matrix is calculated 1 1 2 index of 1st atom, number of L's, L1 > > 2 1 2 dtto for 2nd atom, repeat NATOM times > > 0 0 r-index, (l,s)index == > > On Sun, Sep 6, 2015 at 9:14 PM, Gavin Abo <gs...@crimson.ua.edu gs...@crimson.ua.edu>> wrote: > > Many thanks for your guidance. Actually my system has magnetic (2) and > non-magnetic (3) species. As B_ext. means we are apply magnetic field on > the whole system then why do we need to select natorb = 2 ? > > Bext is applied to the iatoms (i.e., in atomic spheres) that you specify > in case.inorb. The program searches for file case.vorbup, if it finds it, > Bext energy is add to Vxc in atomic spheres and in interstitial region [ > http://www.wien2k.at/reg_user/textbooks/orbital_potentials.pdf (section > "4.1 LAPW0 package" on page 6)]. > > Secondly could you please clarify to me about "adjusting the "direction of > Bext in terms of lattice vectors" line in case.inorb. ". Any example please > or guidance that how to make it. > > For example, > > y = x*tan(theta) = 1*tan(32 degrees) = 0.62487 [ > https://en.wikipedia.org/wiki/Trigonometry ] > > Consider a cubic lattice with the "direction of Bext in terms of lattice > vectors" set to: > > 1 0.62487 0 > > Calculation of the angle between vector (1,0,0) and vector (1,0.62487,0) > with octave: > > username@computername:~/wiendata/case$ octave > octave:1> a=[1 0 0] > a = >1 0 0 > octave:2> b=[1 0.62487 0] > b = >1.0 0.62487 0.0 > octave:3> angle_rad=acos(dot(a,b)/(norm(a)*norm(b))) > angle_rad = 0.55851 > octave:4> angle_deg=angle_rad*180/pi > angle_deg = 32.000 > > This gives an angle of 32 degrees with respect to the (100) axis. > > Reference: > http://www.mathworks.com/matlabcentral/newsreader/view_thread/151925 > > [cid:part1.02070908.02000808@crimson.ua.edu] > > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at<mailto:Wien@zeus.theochem.tuwien.ac.at> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Si under strain
Here is the structure file script. Title F LATTICE,NONEQUIV.ATOMS: 2 216_F-43m MODE OF CALC=RELA unit=ang 10.263031 10.263031 10.263031 90.00 90.00 90.00 ATOM 1: X=0. Y=0. Z=0. MULT= 1 ISPLIT= 2 Si NPT= 781 R0=0.0001 RMT= 2.21Z: 14.0 LOCAL ROT MATRIX:1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 ATOM 2: X=0.2500 Y=0.2500 Z=0.2500 MULT= 1 ISPLIT= 2 Si NPT= 781 R0=0.0001 RMT= 2.21Z: 14.0 LOCAL ROT MATRIX:1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 24 NUMBER OF SYMMETRY OPERATIONS 1 0 0 0. 0-1 0 0. 0 0-1 0. 1 1 0 0 0. 0 0-1 0. 0-1 0 0. 2 0 1 0 0. -1 0 0 0. 0 0-1 0. 3 0 0 1 0. -1 0 0 0. 0-1 0 0. 4 0 1 0 0. 0 0-1 0. -1 0 0 0. 5 0 0 1 0. 0-1 0 0. -1 0 0 0. 6 0-1 0 0. 1 0 0 0. 0 0-1 0. 7 0 0-1 0. 1 0 0 0. 0-1 0 0. 8 -1 0 0 0. 0 1 0 0. 0 0-1 0. 9 -1 0 0 0. 0 0 1 0. 0-1 0 0. 10 0-1 0 0. 0 0-1 0. 1 0 0 0. 11 0 0-1 0. 0-1 0 0. 1 0 0 0. 12 0 0 1 0. 0 1 0 0. 1 0 0 0. 13 0 1 0 0. 0 0 1 0. 1 0 0 0. 14 -1 0 0 0. 0 0-1 0. 0 1 0 0. 15 -1 0 0 0. 0-1 0 0. 0 0 1 0. 16 0 0 1 0. 1 0 0 0. 0 1 0 0. 17 0 1 0 0. 1 0 0 0. 0 0 1 0. 18 0 0-1 0. 0 1 0 0. -1 0 0 0. 19 0-1 0 0. 0 0 1 0. -1 0 0 0. 20 0 0-1 0. -1 0 0 0. 0 1 0 0. 21 0-1 0 0. -1 0 0 0. 0 0 1 0. 22 1 0 0 0. 0 0 1 0. 0 1 0 0. 23 1 0 0 0. 0 1 0 0. 0 0 1 0. 24 On Wed, Sep 9, 2015 at 10:52 PM, Gavin Abo <gs...@crimson.ua.edu> wrote: > Did you send the wrong scfmbj.struct? At least for me, "x nn" of WIEN2k > 14.2 gives WARNINGS that the struct file is not okay. > > username@computername:~/wiendata/scfmbj$ ls > scfmbj.struct > username@computername:~/wiendata/scfmbj$ x nn > specify nn-bondlength factor: (usually=2) [and optionally dlimit, dstmax > (about > 1.d-5, 20)] > 2 > DSTMAX: 20.0 > iix,iiy,iiz 4 4 5 41.052124000 >41.05212400051.315155000 > > ATOM 1 Si ATOM 2 Si > RMT( 1)=2.21000 AND RMT( 2)=2.21000 > SUMS TO 4.42000 LT. NN-DIST= 4.44402 > > ATOM 2 Si ATOM 1 Si > RMT( 2)=2.21000 AND RMT( 1)=2.21000 > SUMS TO 4.42000 LT. NN-DIST= 4.44402 > WARNING: Mult not equal. PLEASE CHECK outputnn-file > WARNING: Mult not equal. PLEASE CHECK outputnn-file > WARNING: ityp not equal. PLEASE CHECK outputnn-file > > NN created a new scfmbj.struct_nn file > NN created a new CASE.STRUCT_NN FILE > 0.0u 0.0s 0:01.35 0.0% 0+0k 0+32io 0pf+0w > > > On 9/9/2015 8:48 AM, Muhammad Sajjad wrote: > > Dear user > I computed Si BS for the structure (scfmbj.struct) attached herewith. When > I apply strain (a+a*0.01, b, c-c*0.01) then symmetry changes (space group > 1_P1). WIth this group calculations does not run every times QTLB error > appears and I have done my best to apply each solution for it. > Secondly I got structure from wien2k examples. That works fine but when > strain is applied then it behaves as a metal. Though this work has been > done many times but I can not find my solution. I found VASP study but not > wien2k study. > In VASP POSCAR is > cubic diamond >5.50 > 0.0.49694327921323430.4969432792132343 > 0.49694327921323430.0.4969432792132343 > 0.49694327921323430.49694327921323430. >Si > 2 > Direct > 0.8750 0.8750 0.8750 > 0.1250 0.1250 0.1250 > > I changed this POSCAR to cif and used for wien2k (structure is attached > TEST.struct) , then BS is correct, but after application of strain (1%) > again it becomes metal. > I do not know how to apply strain in POSCAR. > Could you please guide me? > > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > >
Re: [Wien] Si under strain
Actually I downloaded structure from COD base. So you are write its space group is fd3m and structure is Title F LATTICE,NONEQUIV.ATOMS: 1 227_Fd-3m MODE OF CALC=RELA unit=ang 10.263031 10.263031 10.263031 90.00 90.00 90.00 ATOM -1: X=0. Y=0. Z=0. MULT= 4 ISPLIT= 8 ATOM -1:X= 0. Y=0.2500 Z=0.2500 ATOM -1:X= 0.2500 Y=0.2500 Z=0. ATOM -1:X= 0.2500 Y=0. Z=0.2500 Si NPT= 781 R0=0.0001 RMT=1.8000 Z: 14.0 LOCAL ROT MATRIX:1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 0 NUMBER OF SYMMETRY OPERATIONS . But strain application ((a =5.49 , b= 5.43 , c = 5.38)) does not allow fd-3m, definitely symmetry changed. If is make it primitive, it will be ok or not? On Thu, Sep 10, 2015 at 4:24 PM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > Then > https://www.webelements.com/silicon/crystal_structure.html > is wrong ! > > Before you start calculations you should know the crystal structure very > well and you should know about symmetry and space groups. > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Donnerstag, 10. September 2015 14:51 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] Si under strain > > It is 216_F-43m > > On Thu, Sep 10, 2015 at 1:51 PM, Muhammad Sajjad <sajja...@gmail.com > <mailto:sajja...@gmail.com>> wrote: > Here is the structure file script. > > > Title > F LATTICE,NONEQUIV.ATOMS: 2 216_F-43m > MODE OF CALC=RELA unit=ang > 10.263031 10.263031 10.263031 90.00 90.00 90.00 > ATOM 1: X=0. Y=0. Z=0. > MULT= 1 ISPLIT= 2 > Si NPT= 781 R0=0.0001 RMT= 2.21Z: 14.0 > LOCAL ROT MATRIX:1.000 0.000 0.000 > 0.000 1.000 0.000 > 0.000 0.000 1.000 > ATOM 2: X=0.2500 Y=0.2500 Z=0.2500 > MULT= 1 ISPLIT= 2 > Si NPT= 781 R0=0.0001 RMT= 2.21Z: 14.0 > LOCAL ROT MATRIX:1.000 0.000 0.000 > 0.000 1.000 0.000 > 0.000 0.000 1.000 > 24 NUMBER OF SYMMETRY OPERATIONS > 1 0 0 0. > 0-1 0 0. > 0 0-1 0. >1 > 1 0 0 0. > 0 0-1 0. > 0-1 0 0. >2 > 0 1 0 0. > -1 0 0 0. > 0 0-1 0. >3 > 0 0 1 0. > -1 0 0 0. > 0-1 0 0. >4 > 0 1 0 0. > 0 0-1 0. > -1 0 0 0. >5 > 0 0 1 0. > 0-1 0 0. > -1 0 0 0. >6 > 0-1 0 0. > 1 0 0 0. > 0 0-1 0. >7 > 0 0-1 0. > 1 0 0 0. > 0-1 0 0. >8 > -1 0 0 0. > 0 1 0 0. > 0 0-1 0. >9 > -1 0 0 0. > 0 0 1 0. > 0-1 0 0. > 10 > 0-1 0 0. > 0 0-1 0. > 1 0 0 0. > 11 > 0 0-1 0. > 0-1 0 0. > 1 0 0 0. > 12 > 0 0 1 0. > 0 1 0 0. > 1 0 0 0. > 13 > 0 1 0 0. > 0 0 1 0. > 1 0 0 0. > 14 > -1 0 0 0. > 0 0-1 0. > 0 1 0 0. > 15 > -1 0 0 0. > 0-1 0 0. > 0 0 1 0. > 16 > 0 0 1 0. > 1 0 0 0. > 0 1 0 0. > 17 > 0 1 0 0. > 1 0 0 0. > 0 0 1 0. > 18 > 0 0-1 0. > 0 1 0 0. > -1 0 0 0. > 19 > 0-1 0 0. > 0 0 1 0. > -1 0 0 0. > 20 > 0 0-1 0. > -1 0 0 0. > 0 1 0 0. > 21 > 0-1 0 0. > -1 0 0 0. > 0 0 1 0. > 22 > 1 0 0 0. > 0 0 1 0. > 0 1 0 0.000
Re: [Wien] Si under strain
It is 216_F-43m On Thu, Sep 10, 2015 at 1:51 PM, Muhammad Sajjad <sajja...@gmail.com> wrote: > Here is the structure file script. > > > Title > F LATTICE,NONEQUIV.ATOMS: 2 216_F-43m > MODE OF CALC=RELA unit=ang > 10.263031 10.263031 10.263031 90.00 90.00 90.00 > ATOM 1: X=0. Y=0. Z=0. > MULT= 1 ISPLIT= 2 > Si NPT= 781 R0=0.0001 RMT= 2.21Z: 14.0 > LOCAL ROT MATRIX:1.000 0.000 0.000 > 0.000 1.000 0.000 > 0.000 0.000 1.000 > ATOM 2: X=0.2500 Y=0.2500 Z=0.2500 > MULT= 1 ISPLIT= 2 > Si NPT= 781 R0=0.0001 RMT= 2.21Z: 14.0 > LOCAL ROT MATRIX:1.000 0.000 0.000 > 0.000 1.000 0.000 > 0.000 0.000 1.000 > 24 NUMBER OF SYMMETRY OPERATIONS > 1 0 0 0. > 0-1 0 0. > 0 0-1 0. >1 > 1 0 0 0. > 0 0-1 0. > 0-1 0 0. >2 > 0 1 0 0. > -1 0 0 0. > 0 0-1 0. >3 > 0 0 1 0. > -1 0 0 0. > 0-1 0 0. >4 > 0 1 0 0. > 0 0-1 0. > -1 0 0 0. >5 > 0 0 1 0. > 0-1 0 0. > -1 0 0 0. >6 > 0-1 0 0. > 1 0 0 0. > 0 0-1 0. >7 > 0 0-1 0. > 1 0 0 0. > 0-1 0 0. >8 > -1 0 0 0. > 0 1 0 0. > 0 0-1 0. >9 > -1 0 0 0. > 0 0 1 0. > 0-1 0 0. > 10 > 0-1 0 0. > 0 0-1 0. > 1 0 0 0. > 11 > 0 0-1 0. > 0-1 0 0. > 1 0 0 0. > 12 > 0 0 1 0. > 0 1 0 0. > 1 0 0 0. > 13 > 0 1 0 0. > 0 0 1 0. > 1 0 0 0. > 14 > -1 0 0 0. > 0 0-1 0. > 0 1 0 0. > 15 > -1 0 0 0. > 0-1 0 0. > 0 0 1 0. > 16 > 0 0 1 0. > 1 0 0 0. > 0 1 0 0. > 17 > 0 1 0 0. > 1 0 0 0. > 0 0 1 0. > 18 > 0 0-1 0. > 0 1 0 0. > -1 0 0 0. > 19 > 0-1 0 0. > 0 0 1 0. > -1 0 0 0. > 20 > 0 0-1 0. > -1 0 0 0. > 0 1 0 0. > 21 > 0-1 0 0. > -1 0 0 0. > 0 0 1 0. > 22 > 1 0 0 0. > 0 0 1 0. > 0 1 0 0. > 23 > 1 0 0 0. > 0 1 0 0. > 0 0 1 0. > 24 > > > On Wed, Sep 9, 2015 at 10:52 PM, Gavin Abo <gs...@crimson.ua.edu> wrote: > >> Did you send the wrong scfmbj.struct? At least for me, "x nn" of WIEN2k >> 14.2 gives WARNINGS that the struct file is not okay. >> >> username@computername:~/wiendata/scfmbj$ ls >> scfmbj.struct >> username@computername:~/wiendata/scfmbj$ x nn >> specify nn-bondlength factor: (usually=2) [and optionally dlimit, dstmax >> (about >> 1.d-5, 20)] >> 2 >> DSTMAX: 20.0 >> iix,iiy,iiz 4 4 5 41.052124000 >>41.05212400051.315155000 >> >> ATOM 1 Si ATOM 2 Si >> RMT( 1)=2.21000 AND RMT( 2)=2.21000 >> SUMS TO 4.42000 LT. NN-DIST= 4.44402 >> >> ATOM 2 Si ATOM 1 Si >> RMT( 2)=2.21000 AND RMT( 1)=2.21000 >> SUMS TO 4.42000 LT. NN-DIST= 4.44402 >> WARNING: Mult not equal. PLEASE CHECK outputnn-file >> WARNING: Mult not equal. PLEASE CHECK outputnn-file >> WARNING: ityp not equal. PLEASE CHECK outputnn-file >> >> NN created a new scfmbj.struct_nn file >> NN created a new CASE.STRUCT_NN FILE >> 0.0u 0.0s 0:01.35 0.0% 0+0k 0+32io 0pf+0w >> >> >> On 9/9/2015 8:48 AM, Muhammad Sajjad wrote: >> >> Dear user >> I computed Si BS for the structure (scfmbj.struct) attached herewith. >> When I apply strain (a+a*0.01, b, c-c*0.01) then symmetry changes (space >> group 1_P1). WIth this group calculations does not run every times QTLB >> error appears and I have done my best to apply each solution for it. >> Secondly I got structure from wien2k examples. That works fine but when >> strain is applied then it behaves as a metal. Though this work has been >> done many times but I can not find my solution. I found VASP s
Re: [Wien] EFG dependence on magnetic field
Sorry I can not understand well. As Magnetic Energy = (0.2)(meu)(B). With meu = 58 x 10^-6 eV / T and B = 8T, Magnetic Energy added to my system = 0.000232 eV And What I found from ground state energies of my system before B application (-89574.89353469 Ry) and after B application (-89574.89351208 Ry). Their difference E_before - E_after gives −0.000307496 eV. It means supposed energy to be added to my system was 0.000232 eV but computational value is −0.000307496 eV. I am sorry it seems I am total wrong in replying your answer. Please correct me. Many thanks On Thu, Sep 10, 2015 at 2:33 PM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > As Gavin told you, one does not add a magnetic field but a magnetic energy > to the up and down parts of the potential. > I just wanted to know whether you thought about the added energy in > comparison to the magnetic energy of your system. > For example, in Fe you have a "magnetic energy" in the ordre of 1 Ry and I > guess the energy added by a field of 8 T will not change the result of a > calculation. > > Note that the symmetry reduction by the magnetic field is not accounted > for when you do not use an SO calculation. > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Donnerstag, 10. September 2015 09:34 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] EFG dependence on magnetic field > > Dear Fecher > Sorry I can not understand about which energy your are asking about. My > intention was to see how B effects on EFG. > > By the way > > Without B_ext > TOTAL ENERGY IN Ry = -89574.89353469 > > with B_ext > TOTAL ENERGY IN Ry = -89574.89351208 > > On Wed, Sep 9, 2015 at 1:14 PM, Fecher, Gerhard <fec...@uni-mainz.de > <mailto:fec...@uni-mainz.de>> wrote: > Did you check how much the energy is changed if you apply an external > field of 8 T > the Bohr magneton is approximately 58 x 10^-6 eV / T > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at> [ > wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at>] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com<mailto:sajja...@gmail.com>] > Gesendet: Montag, 7. September 2015 09:55 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] EFG dependence on magnetic field > > Dear Gavin > Bundle of thanks for such a helping reply. Are the case.inorb and > case.indm files are ok now for B-ext applied at 32 degree angle with x-axis? > > > == case.inorb = > 3 2 0 nmod, natorb, ipr > PRATT 1.0BROYD/PRATT, mixing > 1 1 2 iatom nlorb, lorb > > 2 1 2 iatom nlorb, lorb > > 8. Bext 1. 0.62487 0. direction > == == case.indm > = -9. Emin cutoff energy 2 number of atoms for which > density matrix is calculated 1 1 2 index of 1st atom, number of L's, L1 > > 2 1 2 dtto for 2nd atom, repeat NATOM times > > 0 0 r-index, (l,s)index == > > On Sun, Sep 6, 2015 at 9:14 PM, Gavin Abo <gs...@crimson.ua.edu gs...@crimson.ua.edu><mailto:gs...@crimson.ua.edu gs...@crimson.ua.edu>>> wrote: > > Many thanks for your guidance. Actually my system has magnetic (2) and > non-magnetic (3) species. As B_ext. means we are apply magnetic field on > the whole system then why do we need to select natorb = 2 ? > > Bext is applied to the iatoms (i.e., in atomic spheres) that you s
Re: [Wien] EFG dependence on magnetic field
Dear Fecher Sorry I can not understand about which energy your are asking about. My intention was to see B effects on EFG. By the way Without B_ext TOTAL ENERGY IN Ry = -89574.89353469 with B_ext TOTAL ENERGY IN Ry = -89574.89351208 On Wed, Sep 9, 2015 at 1:14 PM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > Did you check how much the energy is changed if you apply an external > field of 8 T > the Bohr magneton is approximately 58 x 10^-6 eV / T > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Montag, 7. September 2015 09:55 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] EFG dependence on magnetic field > > Dear Gavin > Bundle of thanks for such a helping reply. Are the case.inorb and > case.indm files are ok now for B-ext applied at 32 degree angle with x-axis? > > > == case.inorb = > 3 2 0 nmod, natorb, ipr > PRATT 1.0BROYD/PRATT, mixing > 1 1 2 iatom nlorb, lorb > > 2 1 2 iatom nlorb, lorb > > 8. Bext 1. 0.62487 0. direction > == == case.indm > = -9. Emin cutoff energy 2 number of atoms for which > density matrix is calculated 1 1 2 index of 1st atom, number of L's, L1 > > 2 1 2 dtto for 2nd atom, repeat NATOM times > > 0 0 r-index, (l,s)index == > > On Sun, Sep 6, 2015 at 9:14 PM, Gavin Abo <gs...@crimson.ua.edu gs...@crimson.ua.edu>> wrote: > > Many thanks for your guidance. Actually my system has magnetic (2) and > non-magnetic (3) species. As B_ext. means we are apply magnetic field on > the whole system then why do we need to select natorb = 2 ? > > Bext is applied to the iatoms (i.e., in atomic spheres) that you specify > in case.inorb. The program searches for file case.vorbup, if it finds it, > Bext energy is add to Vxc in atomic spheres and in interstitial region [ > http://www.wien2k.at/reg_user/textbooks/orbital_potentials.pdf (section > "4.1 LAPW0 package" on page 6)]. > > Secondly could you please clarify to me about "adjusting the "direction of > Bext in terms of lattice vectors" line in case.inorb. ". Any example please > or guidance that how to make it. > > For example, > > y = x*tan(theta) = 1*tan(32 degrees) = 0.62487 [ > https://en.wikipedia.org/wiki/Trigonometry ] > > Consider a cubic lattice with the "direction of Bext in terms of lattice > vectors" set to: > > 1 0.62487 0 > > Calculation of the angle between vector (1,0,0) and vector (1,0.62487,0) > with octave: > > username@computername:~/wiendata/case$ octave > octave:1> a=[1 0 0] > a = >1 0 0 > octave:2> b=[1 0.62487 0] > b = >1.0 0.62487 0.0 > octave:3> angle_rad=acos(dot(a,b)/(norm(a)*norm(b))) > angle_rad = 0.55851 > octave:4> angle_deg=angle_rad*180/pi > angle_deg = 32.000 > > This gives an angle of 32 degrees with respect to the (100) axis. > > Reference: > http://www.mathworks.com/matlabcentral/newsreader/view_thread/151925 > > [cid:part1.02070908.02000808@crimson.ua.edu] > > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at<mailto:Wien@zeus.theochem.tuwien.ac.at> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] EFG dependence on magnetic field
Dear Fecher Sorry I can not understand about which energy your are asking about. My intention was to see B effects on EFG. Without B_ext By the way with B_ext TOTAL ENERGY IN Ry = -89574.89351208 On Wed, Sep 9, 2015 at 1:14 PM, Fecher, Gerhard <fec...@uni-mainz.de> wrote: > Did you check how much the energy is changed if you apply an external > field of 8 T > the Bohr magneton is approximately 58 x 10^-6 eV / T > > Ciao > Gerhard > > DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: > "I think the problem, to be quite honest with you, > is that you have never actually known what the question is." > > > Dr. Gerhard H. Fecher > Institut of Inorganic and Analytical Chemistry > Johannes Gutenberg - University > 55099 Mainz > and > Max Planck Institute for Chemical Physics of Solids > 01187 Dresden > > Von: wien-boun...@zeus.theochem.tuwien.ac.at [ > wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Muhammad Sajjad [ > sajja...@gmail.com] > Gesendet: Montag, 7. September 2015 09:55 > An: A Mailing list for WIEN2k users > Betreff: Re: [Wien] EFG dependence on magnetic field > > Dear Gavin > Bundle of thanks for such a helping reply. Are the case.inorb and > case.indm files are ok now for B-ext applied at 32 degree angle with x-axis? > > > == case.inorb = > 3 2 0 nmod, natorb, ipr > PRATT 1.0BROYD/PRATT, mixing > 1 1 2 iatom nlorb, lorb > > 2 1 2 iatom nlorb, lorb > > 8. Bext 1. 0.62487 0. direction > == == case.indm > = -9. Emin cutoff energy 2 number of atoms for which > density matrix is calculated 1 1 2 index of 1st atom, number of L's, L1 > > 2 1 2 dtto for 2nd atom, repeat NATOM times > > 0 0 r-index, (l,s)index == > > On Sun, Sep 6, 2015 at 9:14 PM, Gavin Abo <gs...@crimson.ua.edu gs...@crimson.ua.edu>> wrote: > > Many thanks for your guidance. Actually my system has magnetic (2) and > non-magnetic (3) species. As B_ext. means we are apply magnetic field on > the whole system then why do we need to select natorb = 2 ? > > Bext is applied to the iatoms (i.e., in atomic spheres) that you specify > in case.inorb. The program searches for file case.vorbup, if it finds it, > Bext energy is add to Vxc in atomic spheres and in interstitial region [ > http://www.wien2k.at/reg_user/textbooks/orbital_potentials.pdf (section > "4.1 LAPW0 package" on page 6)]. > > Secondly could you please clarify to me about "adjusting the "direction of > Bext in terms of lattice vectors" line in case.inorb. ". Any example please > or guidance that how to make it. > > For example, > > y = x*tan(theta) = 1*tan(32 degrees) = 0.62487 [ > https://en.wikipedia.org/wiki/Trigonometry ] > > Consider a cubic lattice with the "direction of Bext in terms of lattice > vectors" set to: > > 1 0.62487 0 > > Calculation of the angle between vector (1,0,0) and vector (1,0.62487,0) > with octave: > > username@computername:~/wiendata/case$ octave > octave:1> a=[1 0 0] > a = >1 0 0 > octave:2> b=[1 0.62487 0] > b = >1.0 0.62487 0.0 > octave:3> angle_rad=acos(dot(a,b)/(norm(a)*norm(b))) > angle_rad = 0.55851 > octave:4> angle_deg=angle_rad*180/pi > angle_deg = 32.000 > > This gives an angle of 32 degrees with respect to the (100) axis. > > Reference: > http://www.mathworks.com/matlabcentral/newsreader/view_thread/151925 > > [cid:part1.02070908.02000808@crimson.ua.edu] > > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at<mailto:Wien@zeus.theochem.tuwien.ac.at> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > > > > -- > Kind Regards > Muhammad Sajjad > Post Doctoral Fellow > KAUST, KSA. > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Si under strain
Dear user I computed Si BS for the structure (scfmbj.struct) attached herewith. When I apply strain (a+a*0.01, b, c-c*0.01) then symmetry changes (space group 1_P1). WIth this group calculations does not run every times QTLB error appears and I have done my best to apply each solution for it. Secondly I got structure from wien2k examples. That works fine but when strain is applied then it behaves as a metal. Though this work has been done many times but I can not find my solution. I found VASP study but not wien2k study. In VASP POSCAR is cubic diamond 5.50 0.0.49694327921323430.4969432792132343 0.49694327921323430.0.4969432792132343 0.49694327921323430.49694327921323430. Si 2 Direct 0.8750 0.8750 0.8750 0.1250 0.1250 0.1250 I changed this POSCAR to cif and used for wien2k (structure is attached TEST.struct) , then BS is correct, but after application of strain (1%) again it becomes metal. I do not know how to apply strain in POSCAR. Could you please guide me? Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. TEST.struct Description: Binary data scfmbj.struct Description: Binary data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] EFG dependence on magnetic field
Dear Gavin Bundle of thanks for such a helping reply. Are the case.inorb and case.indm files are ok now for B-ext applied at 32 degree angle with x-axis? == case.inorb = 3 2 0 nmod, natorb, ipr PRATT 1.0BROYD/PRATT, mixing 1 1 2 iatom nlorb, lorb 2 1 2 iatom nlorb, lorb 8. Bext 1. 0.62487 0. direction == == case.indm = -9. Emin cutoff energy 2 number of atoms for which density matrix is calculated 1 1 2 index of 1st atom, number of L's, L1 2 1 2 dtto for 2nd atom, repeat NATOM times 0 0 r-index, (l,s)index == On Sun, Sep 6, 2015 at 9:14 PM, Gavin Abo <gs...@crimson.ua.edu> wrote: > > Many thanks for your guidance. Actually my system has magnetic (2) and > non-magnetic (3) species. As B_ext. means we are apply magnetic field on > the whole system then why do we need to select natorb = 2 ? > > > Bext is applied to the iatoms (i.e., in atomic spheres) that you specify > in case.inorb. The program searches for file case.vorbup, if it finds it, > Bext energy is add to Vxc in atomic spheres and in interstitial region [ > http://www.wien2k.at/reg_user/textbooks/orbital_potentials.pdf (section > "4.1 LAPW0 package" on page 6)]. > > Secondly could you please clarify to me about "adjusting the "direction > of Bext in terms of lattice vectors" line in case.inorb. ". Any example > please or guidance that how to make it. > > > For example, > > y = x*tan(theta) = 1*tan(32 degrees) = 0.62487 [ > https://en.wikipedia.org/wiki/Trigonometry ] > > Consider a cubic lattice with the "direction of Bext in terms of lattice > vectors" set to: > > 1 0.62487 0 > > Calculation of the angle between vector (1,0,0) and vector (1,0.62487,0) > with octave: > > username@computername:~/wiendata/case$ octave > octave:1> a=[1 0 0] > a = >1 0 0 > octave:2> b=[1 0.62487 0] > b = >1.0 0.62487 0.0 > octave:3> angle_rad=acos(dot(a,b)/(norm(a)*norm(b))) > angle_rad = 0.55851 > octave:4> angle_deg=angle_rad*180/pi > angle_deg = 32.000 > > This gives an angle of 32 degrees with respect to the (100) axis. > > Reference: > http://www.mathworks.com/matlabcentral/newsreader/view_thread/151925 > > > > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] EFG dependence on magnetic field
Thank you veru much Dear Gavin On Mon, Sep 7, 2015 at 11:25 AM, Gavin Abo <gs...@crimson.ua.edu> wrote: > Yes > > On 9/7/2015 1:55 AM, Muhammad Sajjad wrote: > > Dear Gavin > Bundle of thanks for such a helping reply. Are the case.inorb and > case.indm files are ok now for B-ext applied at 32 degree angle with x-axis? > > == case.inorb = > 3 2 0 nmod, natorb, ipr > PRATT 1.0BROYD/PRATT, mixing > 1 1 2 iatom nlorb, lorb > > 2 1 2 iatom nlorb, lorb > > 8. Bext 1. 0.62487 0. direction > == == case.indm > = -9. Emin cutoff energy 2 number of atoms for which > density matrix is calculated 1 1 2 index of 1st atom, number of L's, L1 > > 2 1 2 dtto for 2nd atom, repeat NATOM times > > 0 0 r-index, (l,s)index == > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] EFG dependence on magnetic field
Dear Gavin Many thanks for your guidance. Actually my system has magnetic (2) and non-magnetic (3) species. As B_ext. means we are apply magnetic field on the whole system then why do we need to select natorb = 2 ? Secondly could you please clarify to me about "adjusting the "direction of Bext in terms of lattice vectors" line in case.inorb. ". Any example please or guidance that how to make it. On Sat, Sep 5, 2015 at 5:18 PM, Gavin Abo <gs...@crimson.ua.edu> wrote: > See comments below: > > Thank you very much for your prompt response. What I understood is > 1. run init_lapw > 2. add case inorb and case.indm files with format as (I have two > magnetic ions Co) > > == case.inorb = > 3 1 0 nmod, natorb, ipr > PRATT 1.0BROYD/PRATT, mixing > 1 1 2 iatom nlorb, lorb > > 2 1 2 iatom nlorb, lorb > > 8. Bext 0. 0. 1. direction == > == case.indm = -9. Emin cutoff energy > 1 number of atoms for which density matrix is calculated 1 1 2 index of > 1st atom, number of L's, L1 > > 2 1 2 index of 1st atom, number of L's, L1 > > 0 0 r-index, (l,s)index == > > 3. use the command runsp_lapw -p -orb -cc 0.0001 -ec 0.0001 > > Am I right? > > > I think your steps look fine, but your inorb and indm files look > problematic. You have two iatom lines (natorb of 2) in case.inorb, but you > only have natorb set for one iatom. This is likely why you are getting the > error "forrtl: severe (59): list-directed I/O syntax error, unit 5, file > /scratch/sajjadm/work/EFG/mfieldBaCoGeO/001B/001B.inorb". In case.indm, > you have two "index of x atom" lines, but you have the "number of atoms" > line set to 1 instead of 2. > > It is also noted that the effect of 8 T might be very small [ > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg11093.html > ]. > > > Can we apply B_ext at some angle say 32 degree with [100]? > > > You should be able to by adjusting the "direction of Bext in terms of > lattice vectors" line in case.inorb. > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] EFG dependence on magnetic field
Dear Gavin Thank you very much for your prompt response. What I understood is 1. run init_lapw 2. add case inorb and case.indm files with format as (I have two magnetic ions Co) == case.inorb = 3 1 0 nmod, natorb, ipr PRATT 1.0BROYD/PRATT, mixing 1 1 2 iatom nlorb, lorb 2 1 2 iatom nlorb, lorb 8. Bext 0. 0. 1. direction == == case.indm = -9. Emin cutoff energy 1 number of atoms for which density matrix is calculated 1 1 2 index of 1st atom, number of L's, L1 2 1 2 index of 1st atom, number of L's, L1 0 0 r-index, (l,s)index == 3. use the command runsp_lapw -p -orb -cc 0.0001 -ec 0.0001 Am I right? Can we apply B_ext at some angle say 32 degree with [100]? On Fri, Sep 4, 2015 at 12:27 AM, Gavin Abo <gs...@crimson.ua.edu> wrote: > To apply external magnetic field Bext inside atomic sphere [1-3] and > control its direction, see the section "Input for interaction with Bext > (nmod=3)" on page 113 in section "7.3.3 Input" of the WIEN2k 14.2 > usersguide [4]. > > [1] > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg01508.html > [2] > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg09205.html > [3] > http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg09002.html > [4] http://www.wien2k.at/reg_user/textbooks/usersguide.pdf > > P.S. The image in your email that I received was broken (missing). > > On 9/3/2015 2:06 PM, Muhammad Sajjad wrote: > >> Dear user >> I am interested to check how electric field gradient behaves when >> external magnetic field is applied. My question is how can we apply >> magnetic field in wien2k? I remember in soc calculation my choose some >> magnetisation direction in case.inrob file. Is it magnetic field or what? >> Can we apply magnetic field at some angle say H is making an angle of 32 >> degree with (100) axis as showmen below? >> >> Thank you >> M. Sajjad >> KAUST, KSA. >> >> >> -- >> Kind Regards >> Muhammad Sajjad >> Post Doctoral Fellow >> KAUST, KSA. >> > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] I/O syntax error in B_ext calculation
Dear Users I am running scf for Ba2CoSi2O7. simply I run 1. init_lapw 2. applied B_ext of 5T in case.inorb file. 3. run the command line runsp_lapw -p -orb -cc 0.001 -ec 0.001 I came across the following error forrtl: severe (59): list-directed I/O syntax error, unit 5, file /scratch/sajjadm/work/EFG/mfieldBaCoGeO/001B/001B.inorb Image PCRoutineLineSource orb004A99DD Unknown Unknown Unknown orb004A84E5 Unknown Unknown Unknown orb004557D0 Unknown Unknown Unknown orb0041C5CA Unknown Unknown Unknown orb0041BDC0 Unknown Unknown Unknown orb0043C20B Unknown Unknown Unknown orb00439D5A Unknown Unknown Unknown orb0040D7E0 init_ 193 init.f orb004039F9 MAIN__103 main.f orb0040341C Unknown Unknown Unknown libc.so.6 00398E81ECDD Unknown Unknown Unknown orb00403319 Unknown Unknown Unknown The case.inorb file is 3 1 0 nmod, natorb, ipr PRATT 1.0BROYD/PRATT, mixing 1 1 2 iatom nlorb, lorb 2 1 2 iatom nlorb, lorb 5. Bext 0. 0. 1.direction I did not run a regular scf before applying B_ext. Is it necessary ?? -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] EFG dependence on magnetic field
Dear user I am interested to check how electric field gradient behaves when external magnetic field is applied. My question is how can we apply magnetic field in wien2k? I remember in soc calculation my choose some magnetisation direction in case.inrob file. Is it magnetic field or what? Can we apply magnetic field at some angle say H is making an angle of 32 degree with (100) axis as showmen below? Thank you M. Sajjad KAUST, KSA. -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] SOC and U calculation
.inorb and .indm 1f. runsp_lapw -orb 1g. save -d GGAUwithoutSO 1h. Select LDA in case.in0 ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] SOC and U calculation
Dear users I am intended to perform FM+U+SOC calculation for Ca2CoSiO7. The step that I understood to perform for such calculations are: 1. run simple SCF with LDA using spin polarized calculations 2. save and run initso_lapw it asks about magnetization direction, Emax value, add RLO for NONE/ ALL/ CHOOSE elements, spin-polarized case, selection of new structure for SO, Number of kpoints. 3. Import .inorb and .indm from SRC_templates for LDA+U calculation and edit them by requirement. 4. Finally use command runsp_lapw -p -so -orb -cc 0.0001 -ec 0.0001 -NI. My queries are: a. Is it the right approach for the mentioned calculation? b. In second step, can we change magnetization direction? how much Emax should be increased (default value is 5 Ry)? Which option is correct for adding RLO? do we need to increase kpoints (before it was 200 for simple scf)? c. Is there any way (like in mBJ we run init_mbj_lapw that imports inm_vresp files and makes changes) to import .inorb and .indm file by command line? d. Is the final command correct? and if someone wants to use GGA+U then simply it is required to run simple SCF by GGA and then apply U or any other method? -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Fixing of layers in relaxation + and Surface reconsrtuction
Dear User How to fix an layer in wien2k during relaxation? In many cases I read the we fixed bottom but allowed the top layer to be relaxed, although it is not very much recommended when we talk about isolated surface, but I guess we can fix the middle layer and allow the top and bottom layers to be relaxed. Other query is the after Si (100) reconstruction one of the two dangling bonds for each Si is left at the surface. Why we do not need to add H (though we called this a stable surface but why?) ? Will it not add metallic behavior? Last one: Is it right that for reconstruction we need to perturb positions of surface Si a few times and the to required to run relaxation? Suppose unperturbed positions for surface two Si (in p(2*1) phase) are (0.960004, 0.960004, 15.7517) and (4.80002, 0.960004, 15.7517) then how much perturbation is necessary and should it be in all direction (x, y , z)? -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Reconstructed Si 100 surface
Thank you very much Fecher and Michael. On Mon, Jul 27, 2015 at 12:11 PM, Michael Sluydts michael.sluy...@ugent.be wrote: It is recommended to make both your slab surfaces identical where possible. The p21, p22 and p24 refer to the type of supercell used to represent the surface. Since you are using a periodic model you are imposing a strong symmetry on the surface, not all reconstructions can be made within your unit cell, indeed for your dimer you will need a 2x1 expansion, i.e. two cells next to eachother along the a axis. The symbols (and possible angles) represent the shape and rotation of the resulting cell (in this case p=primitive, c=centered). You can find more info in for instance: Introduction to Solid State Physics - Kittel - about the description of cells and Density Functional Theory: A Practical Introduction - Sholl and Steckel - for the basics of how to model surfaces Muhammad Sajjad schreef op 27/07/2015 om 10:51: Dear Delamora and Michael Many thanks for informative suggestions. I am intended to study electronic properties of Si (100) of thickness 1 nm. I am asked to introduce dimers instead of using H. But my question is one surface of the Si 100 must be passivated by hydrogen or not?Like If ones introduce dimer on the top (maximum Z) then bottom (minimum Z) must be passivated or not? Can we introduce dimers on both surfaces (top bottom). I am reading the paper ( http://pubs.acs.org/doi/pdf/10.1021/jp509095t ). What does it means P21, P22, P24 reconstructions? Many thanks On Sun, Jul 26, 2015 at 10:51 PM, Michael Sluydts michael.sluy...@ugent.be wrote: Muhammad, When you move the two Si atoms with the dangling bond towards eachother you will create a dimer on the surface, no extra passivation by H atoms is needed, just ensure you make your slab symmetric where possible. If you google Si recontructions you should find very exact geometrical information since this has been done many times before. Best regards, Michael Sluydts Muhammad Sajjad schreef op 26/07/2015 om 17:59: Thank you Pablo But I am required to connect two Si atoms at the surface to reduce one dangling bond and then put H or may be H is not needed if both of the dangling bonds have possibility of omitting. I am reading the paper but could not understand although Fig. 2 is looking easy to understand. WHat is 2*1 or 2*2 or 2*4 ? are the super-cells like 2*1*1 or 2*2*1 etc. On Sun, Jul 26, 2015 at 6:22 PM, delamora delam...@unam.mx wrote: This is the same that I simplified and symmetized. If you are going to add H then you need set the space group as P1 (#1) so the H can move freely and not in a symmetrical path. Now, if you use my proposed structure then you can put inversion (SG P-1, #2) and add two H, one on each surface, but in a symmetric way. Other thing, maybe my proposed structure is too small and your structure is ok, but I would symmetrize as I explained in an earlier message. Pablo -- *De:* wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at en nombre de Muhammad Sajjad sajja...@gmail.com *Enviado:* domingo, 26 de julio de 2015 04:44 a. m. *Para:* wien *Asunto:* [Wien] Reconstructed Si 100 surface Dear Users Kindly guide me how to reconstruct Si 100 layer? The Si 100 structure is attached here with. It has two dangling bonds on both surfaces (Top and Bottom). My intension is to keep one dangling bond on the top and bottom surfaces that will be saturated by H addition. -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing listw...@zeus.theochem.tuwien.ac.athttp://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- ir. Michael Sluydts Center for Molecular Modeling Ghent University Technologiepark 903 9052 Zwijnaarde, Belgium tel. +32 (0)9 264 66 19https://molmod.ugent.be ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing listw...@zeus.theochem.tuwien.ac.athttp://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail
Re: [Wien] Reconstructed Si 100 surface
Dear Delamora and Michael Many thanks for informative suggestions. I am intended to study electronic properties of Si (100) of thickness 1 nm. I am asked to introduce dimers instead of using H. But my question is one surface of the Si 100 must be passivated by hydrogen or not?Like If ones introduce dimer on the top (maximum Z) then bottom (minimum Z) must be passivated or not? Can we introduce dimers on both surfaces (top bottom). I am reading the paper ( http://pubs.acs.org/doi/pdf/10.1021/jp509095t ). What does it means P21, P22, P24 reconstructions? Many thanks On Sun, Jul 26, 2015 at 10:51 PM, Michael Sluydts michael.sluy...@ugent.be wrote: Muhammad, When you move the two Si atoms with the dangling bond towards eachother you will create a dimer on the surface, no extra passivation by H atoms is needed, just ensure you make your slab symmetric where possible. If you google Si recontructions you should find very exact geometrical information since this has been done many times before. Best regards, Michael Sluydts Muhammad Sajjad schreef op 26/07/2015 om 17:59: Thank you Pablo But I am required to connect two Si atoms at the surface to reduce one dangling bond and then put H or may be H is not needed if both of the dangling bonds have possibility of omitting. I am reading the paper but could not understand although Fig. 2 is looking easy to understand. WHat is 2*1 or 2*2 or 2*4 ? are the super-cells like 2*1*1 or 2*2*1 etc. On Sun, Jul 26, 2015 at 6:22 PM, delamora delam...@unam.mx wrote: This is the same that I simplified and symmetized. If you are going to add H then you need set the space group as P1 (#1) so the H can move freely and not in a symmetrical path. Now, if you use my proposed structure then you can put inversion (SG P-1, #2) and add two H, one on each surface, but in a symmetric way. Other thing, maybe my proposed structure is too small and your structure is ok, but I would symmetrize as I explained in an earlier message. Pablo -- *De:* wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at en nombre de Muhammad Sajjad sajja...@gmail.com *Enviado:* domingo, 26 de julio de 2015 04:44 a. m. *Para:* wien *Asunto:* [Wien] Reconstructed Si 100 surface Dear Users Kindly guide me how to reconstruct Si 100 layer? The Si 100 structure is attached here with. It has two dangling bonds on both surfaces (Top and Bottom). My intension is to keep one dangling bond on the top and bottom surfaces that will be saturated by H addition. -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing listw...@zeus.theochem.tuwien.ac.athttp://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- ir. Michael Sluydts Center for Molecular Modeling Ghent University Technologiepark 903 9052 Zwijnaarde, Belgium tel. +32 (0)9 264 66 19https://molmod.ugent.be ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Reconstructed Si 100 surface
Dear Users Kindly guide me how to reconstruct Si 100 layer? The Si 100 structure is attached here with. It has two dangling bonds on both surfaces (Top and Bottom). My intension is to keep one dangling bond on the top and bottom surfaces that will be saturated by H addition. -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. SIopt.struct Description: Binary data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Metallic Behavior of Si 100 layer
Many thanks Delamora for sending the structure. @ Prof. L. Marks My intention is to study electronic band structure of Si (100), say thickness is 0.384 nm, by varying the strain. I am adding H on the top and bottom to saturate the bonds otherwise free electrons will make Si (100) metallic. I guess the received structures will work fine with my problem. Am I right? True Regards On Mon, Jul 20, 2015 at 5:14 PM, Laurence Marks l-ma...@northwestern.edu wrote: This discussion appears to be useful, and hopefully the structures people are sending will help. That said, there is one deeper issue with this thread. If the intent is to calculate a thin hydrogen terminated silicon layer for which there is experimental data then the model(s) are reasonable. It is also relevant if the idea is to learn about surface calculations. However, if the intent is to model a hydrogen terminated Si (001) surface it is extremely unlikely that any of the models are relevant. --- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University http://www.numis.northwestern.edu Corrosion in 4D http://MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi On Jul 20, 2015 08:54, delamora delam...@unam.mx wrote: Here I send a possible structure, but it needs optimization. Pablo de la Mora -- *De:* wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at en nombre de Muhammad Sajjad sajja...@gmail.com *Enviado:* lunes, 20 de julio de 2015 04:22 a. m. *Para:* A Mailing list for WIEN2k users *Asunto:* Re: [Wien] Metallic Behavior of Si 100 layer Dear Delamore and Oleg Thank you very much for your kind suggestions. @Delamore: The attached structure is in fine symmetry and I am trying this for further calculations and will update once done with the calculations. @Oleg: I understood your point. Dangling bonds and less than required Si-Si distance are prominent reason of sending the Si 100 in metallic category. Could you please guide me how to add second H atom at the surface with Si? If possible please add the structure your made. True Regards M. Sajjad KAUST, KSA. On Mon, Jul 20, 2015 at 12:08 AM, Oleg Rubel oru...@lakeheadu.ca wrote: In addition to that: The surface is truncated in a way that each Si has *two* dangling bonds. Thus two H-atoms per one Si at the surface are needed to passivate dangling bonds. I enclosed a sketch (I hope it will come through) where the proposed changes are shown for one Si atom. Note that the H-atoms are oriented to keep tetrahedral coordination of Si. I hope it will help Oleg On Jul 19, 2015, at 3:00 PM, delamora delam...@unam.mx wrote: Sorry, I wanted to say symmetrize De: wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at en nombre de delamora delam...@unam.mx Enviado: domingo, 19 de julio de 2015 02:56 p. m. Para: A Mailing list for WIEN2k users Asunto: Re: [Wien] Metallic Behavior of Si 100 layer If you systematize your system you will have a much simpler problem; Here in the SiH100.struct I moved the central plane to 0 0 1/2, then I averaged the positions of the atoms above the plane and the equivalent atom below this plane. After this sgroup changed the structure to a simpler and more symmetric one. Now, here you have 4 Si layers, maybe you can grow the thickness of your system. Pablo De: wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at en nombre de Muhammad Sajjad sajja...@gmail.com Enviado: domingo, 19 de julio de 2015 06:16 a. m. Para: A Mailing list for WIEN2k users Asunto: Re: [Wien] Metallic Behavior of Si 100 layer Thank you Oleg, the structure file is attached herewith. On Sun, Jul 19, 2015 at 3:05 PM, Oleg Rubel oru...@lakeheadu.ca wrote: Including the structure file would be helpful. Oleg On Jul 19, 2015 12:24 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear User After relaxation, the Si 100 layer saturated with H is behaving like a metal. I am unable to find its reason or some wrong approach followed in my calculations. Could you please guide me? Many thanks Muhammad Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post
Re: [Wien] Metallic Behavior of Si 100 layer
Dear Delamore and Oleg Thank you very much for your kind suggestions. @Delamore: The attached structure is in fine symmetry and I am trying this for further calculations and will update once done with the calculations. @Oleg: I understood your point. Dangling bonds and less than required Si-Si distance are prominent reason of sending the Si 100 in metallic category. Could you please guide me how to add second H atom at the surface with Si? If possible please add the structure your made. True Regards M. Sajjad KAUST, KSA. On Mon, Jul 20, 2015 at 12:08 AM, Oleg Rubel oru...@lakeheadu.ca wrote: In addition to that: The surface is truncated in a way that each Si has *two* dangling bonds. Thus two H-atoms per one Si at the surface are needed to passivate dangling bonds. I enclosed a sketch (I hope it will come through) where the proposed changes are shown for one Si atom. Note that the H-atoms are oriented to keep tetrahedral coordination of Si. I hope it will help Oleg On Jul 19, 2015, at 3:00 PM, delamora delam...@unam.mx wrote: Sorry, I wanted to say symmetrize De: wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at en nombre de delamora delam...@unam.mx Enviado: domingo, 19 de julio de 2015 02:56 p. m. Para: A Mailing list for WIEN2k users Asunto: Re: [Wien] Metallic Behavior of Si 100 layer If you systematize your system you will have a much simpler problem; Here in the SiH100.struct I moved the central plane to 0 0 1/2, then I averaged the positions of the atoms above the plane and the equivalent atom below this plane. After this sgroup changed the structure to a simpler and more symmetric one. Now, here you have 4 Si layers, maybe you can grow the thickness of your system. Pablo De: wien-boun...@zeus.theochem.tuwien.ac.at wien-boun...@zeus.theochem.tuwien.ac.at en nombre de Muhammad Sajjad sajja...@gmail.com Enviado: domingo, 19 de julio de 2015 06:16 a. m. Para: A Mailing list for WIEN2k users Asunto: Re: [Wien] Metallic Behavior of Si 100 layer Thank you Oleg, the structure file is attached herewith. On Sun, Jul 19, 2015 at 3:05 PM, Oleg Rubel oru...@lakeheadu.ca wrote: Including the structure file would be helpful. Oleg On Jul 19, 2015 12:24 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear User After relaxation, the Si 100 layer saturated with H is behaving like a metal. I am unable to find its reason or some wrong approach followed in my calculations. Could you please guide me? Many thanks Muhammad Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html Si-surface.pdf ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Metallic Behavior of Si 100 layer
Thank you Oleg, the structure file is attached herewith. On Sun, Jul 19, 2015 at 3:05 PM, Oleg Rubel oru...@lakeheadu.ca wrote: Including the structure file would be helpful. Oleg On Jul 19, 2015 12:24 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear User After relaxation, the Si 100 layer saturated with H is behaving like a metal. I am unable to find its reason or some wrong approach followed in my calculations. Could you please guide me? Many thanks Muhammad Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. SiH100.struct Description: Binary data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] core leakage problem
Dear L. Dobysheva Thank you for your suggestions. I have over come the problem. The forces have reached around 1 mRy/**. Now the issue is the Si 100 layer saturated with H but the system is behaving like a metal. I am unable to find its reason or some wrong approach followed in my calculations. On Tue, Jul 14, 2015 at 12:09 PM, Lyudmila Dobysheva lyuk...@mail.ru wrote: On 14.07.2015 10:29, Muhammad Sajjad wrote: I used second method according to my understanding. After structure generation and init_lapw, applied the switch *MSEC1a in case.inm* and The userguide says that after init lapw you should make run lapw -fc 1 [another runXX script or additional options are of course also possible] then run the calculations using *run_lapw -fc 0.5 -ec 0.0001 -cc 0.001 80* I don't know what is 80 here. But the calculations are not proceeding and what I found from*case.dayfile* is Calculating 2relax in /home/sajjadm/sajjad/flexibleSi/silayer/2relax on kw12165 with PID 9457 using WIEN2k_14.2 (Release 15/10/2014) in /opt/share/WIEN2K/v14.2-ifort-11.1.075.app 2relax.inM and .minrestart have been created by pairhess lapw1 -p -c(21:44:23) starting parallel lapw1 at Mon Jul 13 21:44:23 AST 2015 - starting parallel LAPW1 jobs at Mon Jul 13 21:44:23 AST 2015 running LAPW1 in parallel mode (using .machines) 5 number_of_parallel_jobs [1] 14781 [2] 14797 [3] 14813 [4] 14829 [5] 14845 In terminal, top command is not showing 5 processors running lapw1 but before it was going on. I am afraid nobody is able to find the reason with this info. Maybe this is an occasional system crash or something else. Try and find it. Please clarify that how to use* the shell script min lapw, together with the program mini*?I mean which* command* do I need to run? As far as I understand you, this is a first method, while you used a second one? Userguide: When using the second method we recommend you read carefully $WIEN-ROOT/SRC mixer/README 5.2.pdf. Section 5.3.2 Minimization of internal parameters describes the first method. Best wishes Lyudmila Dobysheva -- Phys.-Techn. Institute of Ural Br. of Russian Ac. of Sci. 426001 Izhevsk, ul.Kirova 132 RUSSIA -- Tel.:7(3412) 432045(office), 722529(Fax) E-mail: l...@ftiudm.ru, lyuk...@mail.ru (office) lyuk...@gmail.com (home) Skype: lyuka17 (home), lyuka18 (office) http://ftiudm.ru/content/view/25/103/lang,english/ -- ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Metallic Behavior of Si 100 layer
Dear User After relaxation, the Si 100 layer saturated with H is behaving like a metal. I am unable to find its reason or some wrong approach followed in my calculations. Could you please guide me? Many thanks Muhammad Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] core leakage problem
I used the command run_lapw -p -I -i 60 -fc 1. DO I need to change it ? On Mon, Jul 13, 2015 at 11:33 AM, Laurence Marks l-ma...@northwestern.edu wrote: No. On Jul 13, 2015 4:37 PM, Muhammad Sajjad sajja...@gmail.com wrote: Thank you professor L.Marks I followed your first two points and ran calculation. I did not fix the center of slab but relax the whole structure and the minimized forces are :FOR001: 1.ATOM 60.004776 -0.02 -0.03 60.004776 total forces :FOR002: 2.ATOM 33.431081 -0.03 0.03 33.431081 total forces :FOR003: 3.ATOM 33.431078 0.03 -0.03 33.431078 total forces :FOR004: 4.ATOM 60.004776 0.03 0.03 60.004776 total forces :FOR005: 5.ATOM 33.434049 0.02 0.01 -33.434049 total forces :FOR006: 6.ATOM 60.056705 0.01 -0.01 -60.056705 total forces :FOR007: 7.ATOM 60.056708 -0.01 0.00 -60.056708 total forces :FOR008: 8.ATOM 33.434051 -0.02 -0.01 -33.434051 total forces :FOR009: 9.ATOM 32.291125 0.00 0.00 -32.291125 total forces :FOR010: 10.ATOM 32.291121 0.00 0.00 -32.291121 total forces :FOR011: 11.ATOM 32.312836 0.00 0.00 32.312836 total forces :FOR012: 12.ATOM 32.312836 0.00 0.00 32.312836 total forces Is this calculation reliable? On Sun, Jul 12, 2015 at 4:33 PM, Laurence Marks l-ma...@northwestern.edu wrote: A few points: a) Your Si-H distance is too small, it should be about 1.45 Angstroms and you have 1.0 Angstroms. With a better value most of your problems will probably go away. b) To do a surface calculation you must use the DFT equilibrium lattice parameters, not the bulk ones. If you do not the results will be wrong. c) For a surface calculation you need the distance between the two surfaces across the vacuum gap to be equal to or smaller than the distance between the two surfaces across the crystal. I personally prefer to make the bulk material twice as wide as the vacuum gap, so the atomic relaxations at the center of the crystal are as small as possible. (Some people will fix the center of the slab, but I think that is bad physics.) Your bulk crystal is way too small. On Sun, Jul 12, 2015 at 7:27 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear All I am performing structural relaxation for Si (100) with H at its top and bottom (structure is attached, vacuum is 12 A). I have inspect the mailing list in detail (like http://www.wien2k.at/reg_user/textbooks/DFT_and_LAPW_2nd.pdf AND http://www.mail-archive.com/search?q=core+leakagel=wien%40zeus.theochem.tuwien.ac.at) and have resolved the issue by selecting SE -12, but the calculation hang out after LAPW0. Inspection of lapw1_1.error message is Error in LAPW1. Please suggest me the possible solution. One more thing I found in the list If you have core leakage, inspect in case.outputst whether you can take some high-lying core states as valence states instead Is it possible to make it manually with same SE ? or simple done by lowering down SE as I did? Many thanks -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu Corrosion in 4D: MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] core leakage problem
Thank you Prof. L. Marks As I am not expert in using Wine2k, please suggest me the change, For me I can just change convergence criteria like 0.5 instead of 1. On Mon, Jul 13, 2015 at 1:10 PM, Laurence Marks l-ma...@northwestern.edu wrote: Yes. Please think. --- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University http://www.numis.northwestern.edu Corrosion in 4D http://MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi On Jul 13, 2015 03:50, Muhammad Sajjad sajja...@gmail.com wrote: I used the command run_lapw -p -I -i 60 -fc 1. DO I need to change it ? On Mon, Jul 13, 2015 at 11:33 AM, Laurence Marks l-ma...@northwestern.edu wrote: No. On Jul 13, 2015 4:37 PM, Muhammad Sajjad sajja...@gmail.com wrote: Thank you professor L.Marks I followed your first two points and ran calculation. I did not fix the center of slab but relax the whole structure and the minimized forces are :FOR001: 1.ATOM 60.004776 -0.02 -0.03 60.004776 total forces :FOR002: 2.ATOM 33.431081 -0.03 0.03 33.431081 total forces :FOR003: 3.ATOM 33.431078 0.03 -0.03 33.431078 total forces :FOR004: 4.ATOM 60.004776 0.03 0.03 60.004776 total forces :FOR005: 5.ATOM 33.434049 0.02 0.01 -33.434049 total forces :FOR006: 6.ATOM 60.056705 0.01 -0.01 -60.056705 total forces :FOR007: 7.ATOM 60.056708 -0.01 0.00 -60.056708 total forces :FOR008: 8.ATOM 33.434051 -0.02 -0.01 -33.434051 total forces :FOR009: 9.ATOM 32.291125 0.00 0.00 -32.291125 total forces :FOR010: 10.ATOM 32.291121 0.00 0.00 -32.291121 total forces :FOR011: 11.ATOM 32.312836 0.00 0.00 32.312836 total forces :FOR012: 12.ATOM 32.312836 0.00 0.00 32.312836 total forces Is this calculation reliable? On Sun, Jul 12, 2015 at 4:33 PM, Laurence Marks l-ma...@northwestern.edu wrote: A few points: a) Your Si-H distance is too small, it should be about 1.45 Angstroms and you have 1.0 Angstroms. With a better value most of your problems will probably go away. b) To do a surface calculation you must use the DFT equilibrium lattice parameters, not the bulk ones. If you do not the results will be wrong. c) For a surface calculation you need the distance between the two surfaces across the vacuum gap to be equal to or smaller than the distance between the two surfaces across the crystal. I personally prefer to make the bulk material twice as wide as the vacuum gap, so the atomic relaxations at the center of the crystal are as small as possible. (Some people will fix the center of the slab, but I think that is bad physics.) Your bulk crystal is way too small. On Sun, Jul 12, 2015 at 7:27 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear All I am performing structural relaxation for Si (100) with H at its top and bottom (structure is attached, vacuum is 12 A). I have inspect the mailing list in detail (like http://www.wien2k.at/reg_user/textbooks/DFT_and_LAPW_2nd.pdf AND http://www.mail-archive.com/search?q=core+leakagel=wien%40zeus.theochem.tuwien.ac.at) and have resolved the issue by selecting SE -12, but the calculation hang out after LAPW0. Inspection of lapw1_1.error message is Error in LAPW1. Please suggest me the possible solution. One more thing I found in the list If you have core leakage, inspect in case.outputst whether you can take some high-lying core states as valence states instead Is it possible to make it manually with same SE ? or simple done by lowering down SE as I did? Many thanks -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu Corrosion in 4D: MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA
[Wien] core leakage problem
Dear All I am performing structural relaxation for Si (100) with H at its top and bottom (structure is attached, vacuum is 12 A). I have inspect the mailing list in detail (like http://www.wien2k.at/reg_user/textbooks/DFT_and_LAPW_2nd.pdf AND http://www.mail-archive.com/search?q=core+leakagel=wien%40zeus.theochem.tuwien.ac.at) and have resolved the issue by selecting SE -12, but the calculation hang out after LAPW0. Inspection of lapw1_1.error message is Error in LAPW1. Please suggest me the possible solution. One more thing I found in the list If you have core leakage, inspect in case.outputst whether you can take some high-lying core states as valence states instead Is it possible to make it manually with same SE ? or simple done by lowering down SE as I did? Many thanks Sicopy.struct Description: Binary data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Change in Si structure
Dear Users I uploaded cif file (attached herewith) in Wine2k 14 and after initialization it changes (attached in .struct formate) completely. What will be reason behind it? I have followed the post https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg02749.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. si227.cif Description: CIF chemical test opt12.struct Description: Binary data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Fwd: EFG and MM
Dear Users Are the values of asymmetry parameter (eta) or Vxx, Vyy, Vzz printed in some file? I found VZZ001 is printed in case.scf file like :VZZ001: EFG INSIDE SPHERE 1 = 3.735466 UP TO R = 2.5 I am trying to understand http://www.wien2k.at/reg_user/faq/efg2.pdf but facing some problems to understand it very well. M. Sajjad On Thu, Jul 2, 2015 at 4:40 PM, t...@theochem.tuwien.ac.at wrote: The last value of :EFG001, :EFG002, etc. in case.scf is the one from the last iteration, thus this is this value that you have to consider. The EFG is a quantity assigned to a particular nucleus and there is no total EFG defined for the unit cell. The EFG can be negative or positive. Maybe you should read more literature on EFG, e.g., http://www.wien2k.at/reg_user/faq/efg2.pdf On Thu, 2 Jul 2015, Muhammad Sajjad wrote: Dear Tran I guess the second value (like in :EFG001:Co1 EFG = 11.74814 *10**21 V / m**2 :EFG001:Co1 EFG = 11.74810 *10**21 V / m**2) is the correct value of EFG that is 11.74810 *10**21 V / m**2 in this case. Actually, SCF converged after 25 cycles and for each atom 25 values of EFG printed. Am I right? M. Sajjad On Thu, Jul 2, 2015 at 12:19 PM, Muhammad Sajjad sajja...@gmail.com wrote: Dear TranMany thanks for your prompt reply. I have completely understood how to compute MM. About EFG: I guess we can not write total EFG for the whole system, I mean only one value of EFG ? Am I right, as it is a related to every nucleus behavior against external EF? I used the functional WCGGA. How can we select the EFG value for each nucleus if it is repeating as in my case like :EFG001:Co1EFG =11.74814 *10**21 V / m**2 :EFG001:Co1EFG =11.74810 *10**21 V / m**2 :EFG002:Co2EFG =11.74791 *10**21 V / m**2 :EFG002:Co2EFG =11.74777 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19830 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19828 *10**21 V / m**2 Why are the values of EFG for Ba negative? I have run calculation for EFG of In using LDA and the out put is again repeating as :EFG001:EFG =-0.00432 *10**21 V / m**2 :EFG001:EFG = 0.13141 *10**21 V / m**2 :EFG001:EFG = 0.47071 *10**21 V / m**2 :EFG001:EFG = 2.66825 *10**21 V / m**2 :EFG001:EFG = 3.74518 *10**21 V / m**2 :EFG001:EFG = 3.75462 *10**21 V / m**2 :EFG001:EFG = 3.73242 *10**21 V / m**2 True Regards On Wed, Jul 1, 2015 at 3:52 PM, t...@theochem.tuwien.ac.at wrote: Hi, :EFGn (n=001, 002, etc.) is the value of the EFG at the nucleus with index 1, 2, etc. as listed in case.struct (the integer next to ATOM). I can see that the EFG is printed twice for each atom, which is the case when lapw0 is run twice per iterations like for hybrid functionals. Maybe you are running hybrid functionals? :MMIn is the spin magnetic moment inside the sphere surrounding nucleus with index n. :MMTOT is equal to :MMTOT = :MMI001*mult(001) + :MMI002*mult(002) + ... + :MMINT where mult(n) is the multiplicity of atom n (MULT= in case.struct) F. Tran On Wed, 1 Jul 2015, Muhammad Sajjad wrote: Dear Users I am supposed to compute Electric field gradient and Magnetic moment for Ba2CoSi2O7. The spin quantum numbers for Ba and Co are 3/2 and 7/2 both definitely 1. SO their nuclear quadrupole will interact with EFG. I run the calculations as suggest in UG for In. After using the command grep :EFG *.scf I got many values (given below). It is confusing which value should be selected and for which atom Ba or Co. :EFG001:EFG =11.74814 *10**21 V / m**2 :EFG001:EFG =11.74810 *10**21 V / m**2 :EFG002:EFG =11.74791 *10**21 V / m**2 :EFG002:EFG =11.74777 *10**21 V / m**2 :EFG003:EFG = -15.19830 *10**21 V / m**2 :EFG003:EFG = -15.19828 *10**21 V / m**2 :EFG004:EFG = 7.27133 *10**21 V / m**2 :EFG004:EFG = 7.27130 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG005
Re: [Wien] EFG and MM
Dear Users Are the values of asymmetry parameter (eta) or Vxx, Vyy, Vzz printed in some file? I found VZZ001 is printed in case.scf file like :VZZ001: EFG INSIDE SPHERE 1 = 3.735466 UP TO R = 2.5 I am trying to understand http://www.wien2k.at/reg_user/faq/efg2.pdf but facing some problems to understand it very well. M. Sajjad On Thu, Jul 2, 2015 at 4:40 PM, t...@theochem.tuwien.ac.at wrote: The last value of :EFG001, :EFG002, etc. in case.scf is the one from the last iteration, thus this is this value that you have to consider. The EFG is a quantity assigned to a particular nucleus and there is no total EFG defined for the unit cell. The EFG can be negative or positive. Maybe you should read more literature on EFG, e.g., http://www.wien2k.at/reg_user/faq/efg2.pdf On Thu, 2 Jul 2015, Muhammad Sajjad wrote: Dear Tran I guess the second value (like in :EFG001:Co1 EFG = 11.74814 *10**21 V / m**2 :EFG001:Co1 EFG = 11.74810 *10**21 V / m**2) is the correct value of EFG that is 11.74810 *10**21 V / m**2 in this case. Actually, SCF converged after 25 cycles and for each atom 25 values of EFG printed. Am I right? M. Sajjad On Thu, Jul 2, 2015 at 12:19 PM, Muhammad Sajjad sajja...@gmail.com wrote: Dear TranMany thanks for your prompt reply. I have completely understood how to compute MM. About EFG: I guess we can not write total EFG for the whole system, I mean only one value of EFG ? Am I right, as it is a related to every nucleus behavior against external EF? I used the functional WCGGA. How can we select the EFG value for each nucleus if it is repeating as in my case like :EFG001:Co1EFG =11.74814 *10**21 V / m**2 :EFG001:Co1EFG =11.74810 *10**21 V / m**2 :EFG002:Co2EFG =11.74791 *10**21 V / m**2 :EFG002:Co2EFG =11.74777 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19830 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19828 *10**21 V / m**2 Why are the values of EFG for Ba negative? I have run calculation for EFG of In using LDA and the out put is again repeating as :EFG001:EFG =-0.00432 *10**21 V / m**2 :EFG001:EFG = 0.13141 *10**21 V / m**2 :EFG001:EFG = 0.47071 *10**21 V / m**2 :EFG001:EFG = 2.66825 *10**21 V / m**2 :EFG001:EFG = 3.74518 *10**21 V / m**2 :EFG001:EFG = 3.75462 *10**21 V / m**2 :EFG001:EFG = 3.73242 *10**21 V / m**2 True Regards On Wed, Jul 1, 2015 at 3:52 PM, t...@theochem.tuwien.ac.at wrote: Hi, :EFGn (n=001, 002, etc.) is the value of the EFG at the nucleus with index 1, 2, etc. as listed in case.struct (the integer next to ATOM). I can see that the EFG is printed twice for each atom, which is the case when lapw0 is run twice per iterations like for hybrid functionals. Maybe you are running hybrid functionals? :MMIn is the spin magnetic moment inside the sphere surrounding nucleus with index n. :MMTOT is equal to :MMTOT = :MMI001*mult(001) + :MMI002*mult(002) + ... + :MMINT where mult(n) is the multiplicity of atom n (MULT= in case.struct) F. Tran On Wed, 1 Jul 2015, Muhammad Sajjad wrote: Dear Users I am supposed to compute Electric field gradient and Magnetic moment for Ba2CoSi2O7. The spin quantum numbers for Ba and Co are 3/2 and 7/2 both definitely 1. SO their nuclear quadrupole will interact with EFG. I run the calculations as suggest in UG for In. After using the command grep :EFG *.scf I got many values (given below). It is confusing which value should be selected and for which atom Ba or Co. :EFG001:EFG =11.74814 *10**21 V / m**2 :EFG001:EFG =11.74810 *10**21 V / m**2 :EFG002:EFG =11.74791 *10**21 V / m**2 :EFG002:EFG =11.74777 *10**21 V / m**2 :EFG003:EFG = -15.19830 *10**21 V / m**2 :EFG003:EFG = -15.19828 *10**21 V / m**2 :EFG004:EFG = 7.27133 *10**21 V / m**2 :EFG004:EFG = 7.27130 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG005
Re: [Wien] EFG and MM
Dear Tran Many thanks for your prompt reply. I have completely understood how to compute MM. About EFG: I guess we can not write total EFG for the whole system, I mean only one value of EFG ? Am I right, as it is a related to every nucleus behavior against external EF? I used the functional WCGGA. How can we select the EFG value for each nucleus if it is repeating as in my case like :EFG001:Co1EFG =11.74814 *10**21 V / m**2 :EFG001:Co1EFG =11.74810 *10**21 V / m**2 :EFG002:Co2EFG =11.74791 *10**21 V / m**2 :EFG002:Co2EFG =11.74777 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19830 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19828 *10**21 V / m**2 Why are the values of EFG for Ba negative? I have run calculation for EFG of In using LDA and the out put is again repeating as :EFG001:EFG =-0.00432 *10**21 V / m**2 :EFG001:EFG = 0.13141 *10**21 V / m**2 :EFG001:EFG = 0.47071 *10**21 V / m**2 :EFG001:EFG = 2.66825 *10**21 V / m**2 :EFG001:EFG = 3.74518 *10**21 V / m**2 :EFG001:EFG = 3.75462 *10**21 V / m**2 :EFG001:EFG = 3.73242 *10**21 V / m**2 True Regards On Wed, Jul 1, 2015 at 3:52 PM, t...@theochem.tuwien.ac.at wrote: Hi, :EFGn (n=001, 002, etc.) is the value of the EFG at the nucleus with index 1, 2, etc. as listed in case.struct (the integer next to ATOM). I can see that the EFG is printed twice for each atom, which is the case when lapw0 is run twice per iterations like for hybrid functionals. Maybe you are running hybrid functionals? :MMIn is the spin magnetic moment inside the sphere surrounding nucleus with index n. :MMTOT is equal to :MMTOT = :MMI001*mult(001) + :MMI002*mult(002) + ... + :MMINT where mult(n) is the multiplicity of atom n (MULT= in case.struct) F. Tran On Wed, 1 Jul 2015, Muhammad Sajjad wrote: Dear Users I am supposed to compute Electric field gradient and Magnetic moment for Ba2CoSi2O7. The spin quantum numbers for Ba and Co are 3/2 and 7/2 both definitely 1. SO their nuclear quadrupole will interact with EFG. I run the calculations as suggest in UG for In. After using the command grep :EFG *.scf I got many values (given below). It is confusing which value should be selected and for which atom Ba or Co. :EFG001:EFG =11.74814 *10**21 V / m**2 :EFG001:EFG =11.74810 *10**21 V / m**2 :EFG002:EFG =11.74791 *10**21 V / m**2 :EFG002:EFG =11.74777 *10**21 V / m**2 :EFG003:EFG = -15.19830 *10**21 V / m**2 :EFG003:EFG = -15.19828 *10**21 V / m**2 :EFG004:EFG = 7.27133 *10**21 V / m**2 :EFG004:EFG = 7.27130 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG007:EFG = 9.37876 *10**21 V / m**2 :EFG007:EFG = 9.37872 *10**21 V / m**2 :EFG008:EFG = 9.37871 *10**21 V / m**2 :EFG008:EFG = 9.37872 *10**21 V / m**2 Also, am I right that the total local magnetic moment for Co and O will be simply obtained by adding for each site like for Co it comes to be 2.43 + 2.43 = 4.68 and similarly for O? :MMINT: MAGNETIC MOMENT IN INTERSTITIAL = 0.46707 :MMI001: MAGNETIC MOMENT IN SPHERE 1 (Co) = 2.42639 :MMI002: MAGNETIC MOMENT IN SPHERE 2 (Co) = 2.42640 :MMI005: MAGNETIC MOMENT IN SPHERE 5 (O) = 0.01014 :MMI006: MAGNETIC MOMENT IN SPHERE 6 (O) = -0.00276 :MMI007: MAGNETIC MOMENT IN SPHERE 7 (O) = 0.07897 :MMI008: MAGNETIC MOMENT IN SPHERE 8 (O) = 0.07898 :MMTOT: SPIN MAGNETIC MOMENT IN CELL = 6.6 -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA
Re: [Wien] EFG and MM
Dear Tran I guess the second value (like in :EFG001:Co1 EFG = 11.74814 *10**21 V / m**2 :EFG001:Co1 EFG = 11.74810 *10**21 V / m**2) is the correct value of EFG that is 11.74810 *10**21 V / m**2 in this case. Actually, SCF converged after 25 cycles and for each atom 25 values of EFG printed. Am I right? M. Sajjad On Thu, Jul 2, 2015 at 12:19 PM, Muhammad Sajjad sajja...@gmail.com wrote: Dear Tran Many thanks for your prompt reply. I have completely understood how to compute MM. About EFG: I guess we can not write total EFG for the whole system, I mean only one value of EFG ? Am I right, as it is a related to every nucleus behavior against external EF? I used the functional WCGGA. How can we select the EFG value for each nucleus if it is repeating as in my case like :EFG001:Co1EFG =11.74814 *10**21 V / m**2 :EFG001:Co1EFG =11.74810 *10**21 V / m**2 :EFG002:Co2EFG =11.74791 *10**21 V / m**2 :EFG002:Co2EFG =11.74777 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19830 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19828 *10**21 V / m**2 Why are the values of EFG for Ba negative? I have run calculation for EFG of In using LDA and the out put is again repeating as :EFG001:EFG =-0.00432 *10**21 V / m**2 :EFG001:EFG = 0.13141 *10**21 V / m**2 :EFG001:EFG = 0.47071 *10**21 V / m**2 :EFG001:EFG = 2.66825 *10**21 V / m**2 :EFG001:EFG = 3.74518 *10**21 V / m**2 :EFG001:EFG = 3.75462 *10**21 V / m**2 :EFG001:EFG = 3.73242 *10**21 V / m**2 True Regards On Wed, Jul 1, 2015 at 3:52 PM, t...@theochem.tuwien.ac.at wrote: Hi, :EFGn (n=001, 002, etc.) is the value of the EFG at the nucleus with index 1, 2, etc. as listed in case.struct (the integer next to ATOM). I can see that the EFG is printed twice for each atom, which is the case when lapw0 is run twice per iterations like for hybrid functionals. Maybe you are running hybrid functionals? :MMIn is the spin magnetic moment inside the sphere surrounding nucleus with index n. :MMTOT is equal to :MMTOT = :MMI001*mult(001) + :MMI002*mult(002) + ... + :MMINT where mult(n) is the multiplicity of atom n (MULT= in case.struct) F. Tran On Wed, 1 Jul 2015, Muhammad Sajjad wrote: Dear Users I am supposed to compute Electric field gradient and Magnetic moment for Ba2CoSi2O7. The spin quantum numbers for Ba and Co are 3/2 and 7/2 both definitely 1. SO their nuclear quadrupole will interact with EFG. I run the calculations as suggest in UG for In. After using the command grep :EFG *.scf I got many values (given below). It is confusing which value should be selected and for which atom Ba or Co. :EFG001:EFG =11.74814 *10**21 V / m**2 :EFG001:EFG =11.74810 *10**21 V / m**2 :EFG002:EFG =11.74791 *10**21 V / m**2 :EFG002:EFG =11.74777 *10**21 V / m**2 :EFG003:EFG = -15.19830 *10**21 V / m**2 :EFG003:EFG = -15.19828 *10**21 V / m**2 :EFG004:EFG = 7.27133 *10**21 V / m**2 :EFG004:EFG = 7.27130 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG007:EFG = 9.37876 *10**21 V / m**2 :EFG007:EFG = 9.37872 *10**21 V / m**2 :EFG008:EFG = 9.37871 *10**21 V / m**2 :EFG008:EFG = 9.37872 *10**21 V / m**2 Also, am I right that the total local magnetic moment for Co and O will be simply obtained by adding for each site like for Co it comes to be 2.43 + 2.43 = 4.68 and similarly for O? :MMINT: MAGNETIC MOMENT IN INTERSTITIAL = 0.46707 :MMI001: MAGNETIC MOMENT IN SPHERE 1 (Co) = 2.42639 :MMI002: MAGNETIC MOMENT IN SPHERE 2 (Co) = 2.42640 :MMI005: MAGNETIC MOMENT IN SPHERE 5 (O) = 0.01014 :MMI006: MAGNETIC MOMENT IN SPHERE 6 (O) = -0.00276 :MMI007: MAGNETIC MOMENT IN SPHERE 7 (O) = 0.07897 :MMI008: MAGNETIC MOMENT IN SPHERE 8 (O) = 0.07898 :MMTOT: SPIN MAGNETIC
Re: [Wien] EFG and MM
Yes Tran Thank you for valuable suggestions . On Thu, Jul 2, 2015 at 4:40 PM, t...@theochem.tuwien.ac.at wrote: The last value of :EFG001, :EFG002, etc. in case.scf is the one from the last iteration, thus this is this value that you have to consider. The EFG is a quantity assigned to a particular nucleus and there is no total EFG defined for the unit cell. The EFG can be negative or positive. Maybe you should read more literature on EFG, e.g., http://www.wien2k.at/reg_user/faq/efg2.pdf On Thu, 2 Jul 2015, Muhammad Sajjad wrote: Dear Tran I guess the second value (like in :EFG001:Co1 EFG = 11.74814 *10**21 V / m**2 :EFG001:Co1 EFG = 11.74810 *10**21 V / m**2) is the correct value of EFG that is 11.74810 *10**21 V / m**2 in this case. Actually, SCF converged after 25 cycles and for each atom 25 values of EFG printed. Am I right? M. Sajjad On Thu, Jul 2, 2015 at 12:19 PM, Muhammad Sajjad sajja...@gmail.com wrote: Dear TranMany thanks for your prompt reply. I have completely understood how to compute MM. About EFG: I guess we can not write total EFG for the whole system, I mean only one value of EFG ? Am I right, as it is a related to every nucleus behavior against external EF? I used the functional WCGGA. How can we select the EFG value for each nucleus if it is repeating as in my case like :EFG001:Co1EFG =11.74814 *10**21 V / m**2 :EFG001:Co1EFG =11.74810 *10**21 V / m**2 :EFG002:Co2EFG =11.74791 *10**21 V / m**2 :EFG002:Co2EFG =11.74777 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19830 *10**21 V / m**2 :EFG003:Ba1EFG = -15.19828 *10**21 V / m**2 Why are the values of EFG for Ba negative? I have run calculation for EFG of In using LDA and the out put is again repeating as :EFG001:EFG =-0.00432 *10**21 V / m**2 :EFG001:EFG = 0.13141 *10**21 V / m**2 :EFG001:EFG = 0.47071 *10**21 V / m**2 :EFG001:EFG = 2.66825 *10**21 V / m**2 :EFG001:EFG = 3.74518 *10**21 V / m**2 :EFG001:EFG = 3.75462 *10**21 V / m**2 :EFG001:EFG = 3.73242 *10**21 V / m**2 True Regards On Wed, Jul 1, 2015 at 3:52 PM, t...@theochem.tuwien.ac.at wrote: Hi, :EFGn (n=001, 002, etc.) is the value of the EFG at the nucleus with index 1, 2, etc. as listed in case.struct (the integer next to ATOM). I can see that the EFG is printed twice for each atom, which is the case when lapw0 is run twice per iterations like for hybrid functionals. Maybe you are running hybrid functionals? :MMIn is the spin magnetic moment inside the sphere surrounding nucleus with index n. :MMTOT is equal to :MMTOT = :MMI001*mult(001) + :MMI002*mult(002) + ... + :MMINT where mult(n) is the multiplicity of atom n (MULT= in case.struct) F. Tran On Wed, 1 Jul 2015, Muhammad Sajjad wrote: Dear Users I am supposed to compute Electric field gradient and Magnetic moment for Ba2CoSi2O7. The spin quantum numbers for Ba and Co are 3/2 and 7/2 both definitely 1. SO their nuclear quadrupole will interact with EFG. I run the calculations as suggest in UG for In. After using the command grep :EFG *.scf I got many values (given below). It is confusing which value should be selected and for which atom Ba or Co. :EFG001:EFG =11.74814 *10**21 V / m**2 :EFG001:EFG =11.74810 *10**21 V / m**2 :EFG002:EFG =11.74791 *10**21 V / m**2 :EFG002:EFG =11.74777 *10**21 V / m**2 :EFG003:EFG = -15.19830 *10**21 V / m**2 :EFG003:EFG = -15.19828 *10**21 V / m**2 :EFG004:EFG = 7.27133 *10**21 V / m**2 :EFG004:EFG = 7.27130 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG007:EFG
[Wien] EFG and MM
Dear Users I am supposed to compute Electric field gradient and Magnetic moment for Ba2CoSi2O7. The spin quantum numbers for Ba and Co are 3/2 and 7/2 both definitely 1. SO their nuclear quadrupole will interact with EFG. I run the calculations as suggest in UG for In. After using the command grep :EFG *.scf I got many values (given below). It is confusing which value should be selected and for which atom Ba or Co. :EFG001:EFG =11.74814 *10**21 V / m**2 :EFG001:EFG =11.74810 *10**21 V / m**2 :EFG002:EFG =11.74791 *10**21 V / m**2 :EFG002:EFG =11.74777 *10**21 V / m**2 :EFG003:EFG = -15.19830 *10**21 V / m**2 :EFG003:EFG = -15.19828 *10**21 V / m**2 :EFG004:EFG = 7.27133 *10**21 V / m**2 :EFG004:EFG = 7.27130 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG005:EFG = 9.02430 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG006:EFG = 9.76192 *10**21 V / m**2 :EFG007:EFG = 9.37876 *10**21 V / m**2 :EFG007:EFG = 9.37872 *10**21 V / m**2 :EFG008:EFG = 9.37871 *10**21 V / m**2 :EFG008:EFG = 9.37872 *10**21 V / m**2 Also, am I right that the total local magnetic moment for Co and O will be simply obtained by adding for each site like for Co it comes to be 2.43 + 2.43 = 4.68 and similarly for O? :MMINT: MAGNETIC MOMENT IN INTERSTITIAL = 0.46707 :MMI001: MAGNETIC MOMENT IN SPHERE 1 *(Co)* = 2.42639 :MMI002: MAGNETIC MOMENT IN SPHERE 2 *(Co)* = 2.42640 :MMI005: MAGNETIC MOMENT IN SPHERE 5 *(O) * = 0.01014 :MMI006: MAGNETIC MOMENT IN SPHERE 6 *(O) * = -0.00276 :MMI007: MAGNETIC MOMENT IN SPHERE 7 *(O) * = 0.07897 :MMI008: MAGNETIC MOMENT IN SPHERE 8 *(O) * = 0.07898 :MMTOT: SPIN MAGNETIC MOMENT IN CELL = 6.6 -- Kind Regards Muhammad Sajjad Post Doctoral Fellow KAUST, KSA. ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Deviation of elastic constants from reference values
Dear Users I am using WINE2K 11 version and computing elastic constants for Fe. The computed elastic constants are c11 = 304.760951 c12 = 118.317319 c44 = 367.255487 I want to compare these values with C11 = 297.83 GPa, C12 = 141.89 GPa, C144 = 106.73 GPa Please suggest the possible solution. I have used 3000 as well as 5000 k-points. In all kinds of calculation, the potential WC.GGA is used. Also I used the structure (attached here with) downloaded from COD crystallography site. True Regards Dr. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Convergence problem
Dear Tran and F. Gerhard Many thanks for your help. With me the Cohesive Energies for NI and Fe are 5.53 eV/atom and 5.92eV/atom On Thu, Jan 1, 2015 at 3:36 PM, t...@theochem.tuwien.ac.at wrote: From the total energies it's difficult to say. It's better to compare the cohesive energies. My cohesive energies for Fe and Ni with the WC functional are 5.54 and 5.38 eV/atom, respectively, with an error bar of 0.05 eV/atom maximum. In this paper you can find PBE cohesive energies calculated with the VASP code (my PBE values agree very well with them): http://journals.aps.org/prb/abstract/10.1103/PhysRevB.87.214102 You have to know which error is acceptable for your purposes. On Thu, 1 Jan 2015, Muhammad Sajjad wrote: Dear Tran Thank you for your further help. I am also attaching here with the structure file. Please let know the E0 values you calculated for Ni and Fe. With me these values are -3040.43215615 and -2544.39535597 (in Ry) respectively. On Tue, Dec 30, 2014 at 4:52 AM, t...@theochem.tuwien.ac.at wrote: Hi, The struct file that I used is attached. For my purpose, the size of the unit cell was large enough to avoid spurious interactions between neighbouring cells. Note that a, b and c are different in order to avoid to high symmetry. F. Tran On Mon, 29 Dec 2014, Muhammad Sajjad wrote: Dear Prof. Marks and F. TranThank you so much for your helpful suggestions. I was already doing the spin polarization calculations. I have got the convergence by using mixing factor 0.1, starting calculation with PRATT and then switched to MSR1 after 7 cycles, and the command runsp_lapw -cc 0.1 -in1ef -i 150 The obtained magnetic moment was 2.9. Dear Tran I would really appreciate if you share some more details about case.struct file. It will definitely be helpful for me as well as for others. Kind Regards Dr. Sajjad On Mon, Dec 29, 2014 at 3:50 AM, t...@theochem.tuwien.ac.at wrote: The calculation for an isolated atom with a code which uses periodic boundary conditions (like WIEN2k) is not trivial. This is what I have done recently for most atoms of the periodic table (excluding f-systems), and for some of the transition-metal atoms this was extremely difficult to achieve convergence. A few hints (that I used): 1) It is important (and necessary depending on the atom) to reduce the symmetry from cubic to, e.g., orthorhombic to be able to access the electronic configuration with the lowest energy (this is the procedure followed by several research groups like VASP for instance). If necessary I can give more details about the case.struct that I used. 2) lapw0 requires a lot of memory, while lapw1 requires both memory and time. To reduce computer time for lapw1, I was using iterative digonalization (this was my command for all atoms): runsp_lapw -ec 0.0001 -cc 0.0001 -it -i 1000 -NI 3) For the Ni atom the magnetic moment should be 2: :MMTOT: SPIN MAGNETIC MOMENT IN CELL =2.00238 4) I was using the default setting for mixer. F. Tran On Sun, 28 Dec 2014, Laurence Marks wrote: Also, only 1 k-point (Gamma), an RMT and RKMax similar to what you use for bulk Ni, particularly as I assume you are doing the calculation to get an enthalpy of formation. You may have to use the mpi versions as it is probably too large for a non-mpi run. If you are doing WC+U (or -ineece) the U (or on-site hybrid) removes the phase transition so convergence should be simple. Straight WC is not a simple calculation because the physics for an isolated Ni atom is wrong. With wrong physics there is in fact no guarantee that the calculation will ever converge! ___ Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought
Re: [Wien] Convergence problem
Dear Tran Thank you for your further help. I am also attaching here with the structure file. Please let know the E0 values you calculated for Ni and Fe. With me these values are -3040.43215615 and -2544.39535597 (in Ry) respectively. On Tue, Dec 30, 2014 at 4:52 AM, t...@theochem.tuwien.ac.at wrote: Hi, The struct file that I used is attached. For my purpose, the size of the unit cell was large enough to avoid spurious interactions between neighbouring cells. Note that a, b and c are different in order to avoid to high symmetry. F. Tran On Mon, 29 Dec 2014, Muhammad Sajjad wrote: Dear Prof. Marks and F. TranThank you so much for your helpful suggestions. I was already doing the spin polarization calculations. I have got the convergence by using mixing factor 0.1, starting calculation with PRATT and then switched to MSR1 after 7 cycles, and the command runsp_lapw -cc 0.1 -in1ef -i 150 The obtained magnetic moment was 2.9. Dear Tran I would really appreciate if you share some more details about case.struct file. It will definitely be helpful for me as well as for others. Kind Regards Dr. Sajjad On Mon, Dec 29, 2014 at 3:50 AM, t...@theochem.tuwien.ac.at wrote: The calculation for an isolated atom with a code which uses periodic boundary conditions (like WIEN2k) is not trivial. This is what I have done recently for most atoms of the periodic table (excluding f-systems), and for some of the transition-metal atoms this was extremely difficult to achieve convergence. A few hints (that I used): 1) It is important (and necessary depending on the atom) to reduce the symmetry from cubic to, e.g., orthorhombic to be able to access the electronic configuration with the lowest energy (this is the procedure followed by several research groups like VASP for instance). If necessary I can give more details about the case.struct that I used. 2) lapw0 requires a lot of memory, while lapw1 requires both memory and time. To reduce computer time for lapw1, I was using iterative digonalization (this was my command for all atoms): runsp_lapw -ec 0.0001 -cc 0.0001 -it -i 1000 -NI 3) For the Ni atom the magnetic moment should be 2: :MMTOT: SPIN MAGNETIC MOMENT IN CELL =2.00238 4) I was using the default setting for mixer. F. Tran On Sun, 28 Dec 2014, Laurence Marks wrote: Also, only 1 k-point (Gamma), an RMT and RKMax similar to what you use for bulk Ni, particularly as I assume you are doing the calculation to get an enthalpy of formation. You may have to use the mpi versions as it is probably too large for a non-mpi run. If you are doing WC+U (or -ineece) the U (or on-site hybrid) removes the phase transition so convergence should be simple. Straight WC is not a simple calculation because the physics for an isolated Ni atom is wrong. With wrong physics there is in fact no guarantee that the calculation will ever converge! ___ Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi On Dec 28, 2014 11:47 AM, Laurence Marks l-ma...@northwestern.edu wrote: You almost certainly need to run spin polarized, probably MSR1 (GREED 0.1), TEMPS (room temp). The convergence is complicated for WC ( simple GGAs) due to an electronic phase transition between sp d occupation near the fixed point. ___ Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi On Dec 28, 2014 10:47 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear Users I am running SCF calculation for Ni with lattice constant of 30 Bohr and using WC-GGA. The calculation is not converging even upto 100 iterations and more. To solve the problem I have performed following steps
Re: [Wien] Convergence problem
Dear Prof. Marks and F. Tran Thank you so much for your helpful suggestions. I was already doing the spin polarization calculations. I have got the convergence by using mixing factor 0.1, starting calculation with PRATT and then switched to MSR1 after 7 cycles, and the command runsp_lapw -cc 0.1 -in1ef -i 150 The obtained magnetic moment was 2.9. Dear Tran I would really appreciate if you share some more details about case.struct file. It will definitely be helpful for me as well as for others. Kind Regards Dr. Sajjad On Mon, Dec 29, 2014 at 3:50 AM, t...@theochem.tuwien.ac.at wrote: The calculation for an isolated atom with a code which uses periodic boundary conditions (like WIEN2k) is not trivial. This is what I have done recently for most atoms of the periodic table (excluding f-systems), and for some of the transition-metal atoms this was extremely difficult to achieve convergence. A few hints (that I used): 1) It is important (and necessary depending on the atom) to reduce the symmetry from cubic to, e.g., orthorhombic to be able to access the electronic configuration with the lowest energy (this is the procedure followed by several research groups like VASP for instance). If necessary I can give more details about the case.struct that I used. 2) lapw0 requires a lot of memory, while lapw1 requires both memory and time. To reduce computer time for lapw1, I was using iterative digonalization (this was my command for all atoms): runsp_lapw -ec 0.0001 -cc 0.0001 -it -i 1000 -NI 3) For the Ni atom the magnetic moment should be 2: :MMTOT: SPIN MAGNETIC MOMENT IN CELL =2.00238 4) I was using the default setting for mixer. F. Tran On Sun, 28 Dec 2014, Laurence Marks wrote: Also, only 1 k-point (Gamma), an RMT and RKMax similar to what you use for bulk Ni, particularly as I assume you are doing the calculation to get an enthalpy of formation. You may have to use the mpi versions as it is probably too large for a non-mpi run. If you are doing WC+U (or -ineece) the U (or on-site hybrid) removes the phase transition so convergence should be simple. Straight WC is not a simple calculation because the physics for an isolated Ni atom is wrong. With wrong physics there is in fact no guarantee that the calculation will ever converge! ___ Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi On Dec 28, 2014 11:47 AM, Laurence Marks l-ma...@northwestern.edu wrote: You almost certainly need to run spin polarized, probably MSR1 (GREED 0.1), TEMPS (room temp). The convergence is complicated for WC ( simple GGAs) due to an electronic phase transition between sp d occupation near the fixed point. ___ Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu MURI4D.numis.northwestern.edu Co-Editor, Acta Cryst A Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi On Dec 28, 2014 10:47 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear Users I am running SCF calculation for Ni with lattice constant of 30 Bohr and using WC-GGA. The calculation is not converging even upto 100 iterations and more. To solve the problem I have performed following steps switched to TEMPS = 0.005 from TETRA Changed MSR1 to PRATT changed mixing factor (increase and decrease from 0.2) Thanking in advance. Dr. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Convergence problem
Dear Users I am running SCF calculation for Ni with lattice constant of 30 Bohr and using WC-GGA. The calculation is not converging even upto 100 iterations and more. To solve the problem I have performed following steps switched to TEMPS = 0.005 from TETRA Changed MSR1 to PRATT changed mixing factor (increase and decrease from 0.2) Thanking in advance. Dr. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] SCF Convergence with mBJ
Dear P. Bala Thank you for your suggestions. Dr. Muhammad Sajjad On Sun, Dec 7, 2014 at 3:51 PM, Peter Blaha pbl...@theochem.tuwien.ac.at wrote: Looks as if he is still using PRATT with a small mixing. As I said before, update to wien2k_14.2 Am 07.12.2014 06:43, schrieb Muhammad Sajjad: Dear F. Tran and P. Bala The user Qaasim has some problem with his email and can not see your emails. I am posting on his behalf. The charge criteria is used (runsp_lapw -cc 0.1 -in1new 2 -i 200) with Wirn2k 11 version. Analysis of parameter: The grep :DIS case.scf generates :DIS in SCFNi50.scf (showing last 10 / 1 lines) --- DIS --- :DIS : CHARGE DISTANCE ( 0.5610332 for atom1 spin 2) 0.2378662 :DIS : CHARGE DISTANCE ( 0.5625705 for atom1 spin 2) 0.2377033 :DIS : CHARGE DISTANCE ( 0.5596997 for atom1 spin 2) 0.2372355 :DIS : CHARGE DISTANCE ( 0.5610506 for atom1 spin 2) 0.2370034 :DIS : CHARGE DISTANCE ( 0.5580173 for atom1 spin 2) 0.2364869 :DIS : CHARGE DISTANCE ( 0.5592747 for atom1 spin 2) 0.2362548 :DIS : CHARGE DISTANCE ( 0.5563927 for atom1 spin 2) 0.2357996 :DIS : CHARGE DISTANCE ( 0.5578176 for atom1 spin 2) 0.2356113 :DIS : CHARGE DISTANCE ( 0.5550182 for atom1 spin 2) 0.2351496 :DIS : CHARGE DISTANCE ( 0.5562844 for atom1 spin 2) 0.2349745 Please suggest the possible solution. Kind Regards Muhammad Sajjad On Fri, Dec 5, 2014 at 5:10 PM, Peter Blaha pbl...@theochem.tuwien.ac.at mailto:pbl...@theochem.tuwien.ac.at wrote: Are you using wien2k_14.x ??? Am 05.12.2014 18:20, schrieb Qasim Mahmood: Dear User Could you please let me know what changes we can make to converge our calculations with mBJ ( at50% doping). I have done almost all the steps that i know, like to change mixing factor, switch to MSRI from PRATT within 6 to 10 cycles, use dense k-mesh etc. My system is a magnetic ternary alloy where the magnetic ion is Co. With many thanks and true regards */ Mr.Qasim Mahmood /* */Ph.D Schollar, PU,Lahore,Pakistan/* _ Wien mailing list w...@zeus.theochem.tuwien.ac.__at mailto:Wien@zeus.theochem. tuwien.ac.at http://zeus.theochem.tuwien.__ac.at/mailman/listinfo/wien http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/__ w...@zeus.theochem.tuwien.ac.__at/index.html http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac. at/index.html -- Peter Blaha Inst.Materials Chemistry TU Vienna Getreidemarkt 9 A-1060 Vienna Austria +43-1-5880115671 tel:%2B43-1-5880115671 _ Wien mailing list w...@zeus.theochem.tuwien.ac.__at mailto:Wien@zeus.theochem. tuwien.ac.at http://zeus.theochem.tuwien.__ac.at/mailman/listinfo/wien http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/__ w...@zeus.theochem.tuwien.ac.__at/index.html http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac. at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html -- - Peter Blaha Inst. Materials Chemistry, TU Vienna Getreidemarkt 9, A-1060 Vienna, Austria Tel: +43-1-5880115671 Fax: +43-1-5880115698 email: pbl...@theochem.tuwien.ac.at - ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] SCF Convergence with mBJ
Dear F. Tran and P. Bala The user Qaasim has some problem with his email and can not see your emails. I am posting on his behalf. The charge criteria is used (runsp_lapw -cc 0.1 -in1new 2 -i 200) with Wirn2k 11 version. Analysis of parameter: The grep :DIS case.scf generates :DIS in SCFNi50.scf (showing last 10 / 1 lines) --- DIS --- :DIS : CHARGE DISTANCE ( 0.5610332 for atom1 spin 2) 0.2378662 :DIS : CHARGE DISTANCE ( 0.5625705 for atom1 spin 2) 0.2377033 :DIS : CHARGE DISTANCE ( 0.5596997 for atom1 spin 2) 0.2372355 :DIS : CHARGE DISTANCE ( 0.5610506 for atom1 spin 2) 0.2370034 :DIS : CHARGE DISTANCE ( 0.5580173 for atom1 spin 2) 0.2364869 :DIS : CHARGE DISTANCE ( 0.5592747 for atom1 spin 2) 0.2362548 :DIS : CHARGE DISTANCE ( 0.5563927 for atom1 spin 2) 0.2357996 :DIS : CHARGE DISTANCE ( 0.5578176 for atom1 spin 2) 0.2356113 :DIS : CHARGE DISTANCE ( 0.5550182 for atom1 spin 2) 0.2351496 :DIS : CHARGE DISTANCE ( 0.5562844 for atom1 spin 2) 0.2349745 Please suggest the possible solution. Kind Regards Muhammad Sajjad On Fri, Dec 5, 2014 at 5:10 PM, Peter Blaha pbl...@theochem.tuwien.ac.at wrote: Are you using wien2k_14.x ??? Am 05.12.2014 18:20, schrieb Qasim Mahmood: Dear User Could you please let me know what changes we can make to converge our calculations with mBJ ( at50% doping). I have done almost all the steps that i know, like to change mixing factor, switch to MSRI from PRATT within 6 to 10 cycles, use dense k-mesh etc. My system is a magnetic ternary alloy where the magnetic ion is Co. With many thanks and true regards */ Mr.Qasim Mahmood /* */Ph.D Schollar, PU,Lahore,Pakistan/* ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html -- Peter Blaha Inst.Materials Chemistry TU Vienna Getreidemarkt 9 A-1060 Vienna Austria +43-1-5880115671 ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] How to get exact 25% doping ?
Yep 31 atoms are not possible in 2*2*2 supercell. I think there should be 64 atoms (replace 8 atoms for 25% doping), as in 1*1*1 supercell 8 atoms(replace 1 atom), in 2*1*1 supercell 16 atoms (replace 2 atoms) and in 2*2*1 super cell 32 atoms (replace 4 atoms). On Sun, May 25, 2014 at 4:44 AM, t...@theochem.tuwien.ac.at wrote: It's not possible to get 31 atoms in a 2*2*2 supercell, unless you manually removed one. I guess it should be 32 in your case. On Sun, 25 May 2014, Naseem Hassan wrote: Dear community Fellows I created a 2*2*2 supercell and got 31 Cd atoms. I want to replace Cd atoms by Zn atoms so that to get 25% of doping. 25% of 31 atoms is 7.75 which is not a whole number. How to obtain whole number to get exact 25% of doping ?. Thank You All Naseem ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] calculation of elastic constants + software installation
Dear user. I am using WIEN2K 12 version on Dell Core I5 desktop (OS is Centos 6.4, 64 bit). As suggested in UG that elastic properties can be computed directly from WIEN2K but for cubic structure. My question is can I do this job with the mentioned version of WIEN2K for Magnetic material in Cubic phase? With Thanks M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] question regarding
Dear user I switched to new work place in another lab with different internet settings. when i wrote w2web in terminal, i received the following message (w2web server started, now point your web browser to http://msajjad:7890), same a before. When I opened the address (http://msajjad:7890), it is not running (The server at www.msajjad.com is taking too long to respond). Do I need to delete w2web folder and set it again or any other solution please? With thanks M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] question regarding
Dear Pieper Thank you for your response. I build new w2web after removing the old one, and now it is working. Kind Regards M. Sajjad On Wed, May 7, 2014 at 6:09 PM, pieper pie...@ifp.tuwien.ac.at wrote: It seems, Wien2k is running on your local machine? My guess is, that you changed the name of your machine ('msajjad') when you moved and should setup w2web new. If Wien2k is not running on your local machine there is probably now a firewall between your new place and the machine running Wien2k. Look into the firewall settings of the machine running Wien2k and the network servers in between. Good look! --- Dr. Martin Pieper Karl-Franzens University Institute of Physics Universitätsplatz 5 A-8010 Graz Austria Tel.: +43-(0)316-380-8564 Am 07.05.2014 09:56, schrieb Muhammad Sajjad: Dear user I switched to new work place in another lab with different internet settings. when i wrote w2web in terminal, i received the following message (w2web server started, now point your web browser to http://msajjad:7890 [1]), same a before. When I opened the address (http://msajjad:7890 [1]), it is not running (The server at www.msajjad.com [2] is taking too long to respond). Do I need to delete w2web folder and set it again or any other solution please? With thanks M. Sajjad Links: -- [1] http://msajjad:7890 [2] http://www.msajjad.com ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Fermi + Syntax error
Dear Sir P. Bala Thank you for your suggestion. I am attaching here the case.scf1 (up down) files as I dont have much idea what to do. I have used Emax = 4.5 and what maximum limit I can use for Emax? Other parameters I have used are RMT*Kmax = 8, Gmax. = 16 18, K-points = 300. Kind Regards M. Sajjad On Wed, Apr 16, 2014 at 2:19 PM, Peter Blaha pbl...@theochem.tuwien.ac.atwrote: Check the eigenvalues in case.scf1 if they are meaningful (enough states for the number of electrons). You may compare them to a successful calculation for another volume. Maybe you need to increase EMAX (case.in1), or simply your starting density for a large volume change is so bad that you cannot run with this density. Try x dstart and then rerun the job. On 04/16/2014 08:07 AM, Muhammad Sajjad wrote: Thank you Sir Gavin Abo. Actually, I am using this version for long time and working quite well. I dont think so that this version is creating problem as it is performing right calculations for other alloys (even from -30 to 15 % vol.change) as well as for the same alloy from -20 to 20 % volume change. As Prof. P. Bala has suggested earlier to change k mesh and Emax. I have done these two jobs but the problem persists. On Wed, Apr 16, 2014 at 4:13 AM, Gavin Abo gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu wrote: Your LOWER and UPPER BOUND values are the same and SOS values are 0. I have seen something similar before. One user was able to solve the problem by using the latest Wien2k version and not an old buggy version [http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac. at/msg09737.html]. I think another user had to increase NUME in SRC_lapw1/param.inc using siteconfig (siteconfig recommends that NUME = sqrt(MEMORY/10)/10, so the computer should have enough random access MEMORY to do so) [http://www.mail-archive.com/wien%40zeus.theochem.tuwien. ac.at/msg07933.html]. A -25% volume change seems quite large, so maybe that could also cause the problem. I would suggest you first try and make sure you can successfully calculate some smaller volume change like -5%, before you continue towards such a large change in volume. On 4/15/2014 12:19 PM, Muhammad Sajjad wrote: Dear Sir Gavin Abo Thank you for you kind response. I have checked these links. What I understood is to increase Emax value. Any other option please? The error file (uplapw2.error) has following message Error in LAPW2 'FERMI' - EFERMI OUT OF ENERGY RANGE 'FERMI' - STOP IN EFI 'FERMI' - ENERGY OF LOWER BOUND : 0.94772 'FERMI' - NUMBER OF STATES AT THE LOWER BOUND : 68.88544 'FERMI' - ENERGY OF UPPER BOUND : 0.94772 'FERMI' - NUMBER OF STATES AT THE UPPER BOUND : 68.87432 'FERMI' - ADD 68.20313 'FERMI' - SOS 0..........000 'FERMI' - NOS ** Kind Regards On Wed, Apr 16, 2014 at 1:21 AM, Gavin Abo gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu wrote: There is usually a more detailed message for 'FERMI - Error' in one of the *lapw2.error files: http://www.mail-archive.com/wien%40zeus.theochem.tuwien. ac.at/msg02361.html http://www.mail-archive.com/wien%40zeus.theochem.tuwien. ac.at/msg08805.html You checked all error files in the case directory (cat *.error), and they were empty? On 4/15/2014 9:38 AM, Muhammad Sajjad wrote: Dear Wien2k Users Can you please pull me out of this error? Many times it is appearing while performing Optimization calculations. The wien2k version is 12 and installed on fedora 17 64 bit. Lattice constant I used is 4.397 and it is ternary alloy with 25 % doping of a TM. The error message is ERROR status in optrnew438_vol_-25.0 stop error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error FERMI - Error LAPW1 END LAPW1 END LAPW0 END in cycle 8ETEST: .12554520 CTEST: .0760248 MIXER END CORE END True Regards M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at mailto:Wien@zeus.theochem. tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http
Re: [Wien] Fermi + Syntax error
Dear Pro. P. Bala Thank you so much for your suggestions. I have finished the calculations without any error. Kind Regards M. Sajjad On Wed, Apr 16, 2014 at 2:50 PM, Peter Blaha pbl...@theochem.tuwien.ac.atwrote: Looks ok. Last suggestion: change in case.in2(c) TETRA 0.000 to TEMP 0.004 (Most likely, you can change this back after a couple of iterations). On 04/16/2014 08:46 AM, Muhammad Sajjad wrote: Dear Sir P. Bala Thank you for your suggestion. I am attaching here the case.scf1 (up down) files as I dont have much idea what to do. I have used Emax = 4.5 and what maximum limit I can use for Emax? Other parameters I have used are RMT*Kmax = 8, Gmax. = 16 18, K-points = 300. Kind Regards M. Sajjad On Wed, Apr 16, 2014 at 2:19 PM, Peter Blaha pbl...@theochem.tuwien.ac.at mailto:pbl...@theochem.tuwien.ac.at wrote: Check the eigenvalues in case.scf1 if they are meaningful (enough states for the number of electrons). You may compare them to a successful calculation for another volume. Maybe you need to increase EMAX (case.in1), or simply your starting density for a large volume change is so bad that you cannot run with this density. Try x dstart and then rerun the job. On 04/16/2014 08:07 AM, Muhammad Sajjad wrote: Thank you Sir Gavin Abo. Actually, I am using this version for long time and working quite well. I dont think so that this version is creating problem as it is performing right calculations for other alloys (even from -30 to 15 % vol.change) as well as for the same alloy from -20 to 20 % volume change. As Prof. P. Bala has suggested earlier to change k mesh and Emax. I have done these two jobs but the problem persists. On Wed, Apr 16, 2014 at 4:13 AM, Gavin Abo gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu wrote: Your LOWER and UPPER BOUND values are the same and SOS values are 0. I have seen something similar before. One user was able to solve the problem by using the latest Wien2k version and not an old buggy version [http://www.mail-archive.com/__wien@zeus.theochem.tuwien.ac. __at/msg09737.html http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac. at/msg09737.html]. I think another user had to increase NUME in SRC_lapw1/param.inc using siteconfig (siteconfig recommends that NUME = sqrt(MEMORY/10)/10, so the computer should have enough random access MEMORY to do so) [http://www.mail-archive.com/__wien%40zeus.theochem.tuwien._ _ac.at/msg07933.html http://www.mail-archive.com/wien%40zeus.theochem.tuwien. ac.at/msg07933.html]. A -25% volume change seems quite large, so maybe that could also cause the problem. I would suggest you first try and make sure you can successfully calculate some smaller volume change like -5%, before you continue towards such a large change in volume. On 4/15/2014 12:19 PM, Muhammad Sajjad wrote: Dear Sir Gavin Abo Thank you for you kind response. I have checked these links. What I understood is to increase Emax value. Any other option please? The error file (uplapw2.error) has following message Error in LAPW2 'FERMI' - EFERMI OUT OF ENERGY RANGE 'FERMI' - STOP IN EFI 'FERMI' - ENERGY OF LOWER BOUND : 0.94772 'FERMI' - NUMBER OF STATES AT THE LOWER BOUND : 68.88544 'FERMI' - ENERGY OF UPPER BOUND : 0.94772 'FERMI' - NUMBER OF STATES AT THE UPPER BOUND : 68.87432 'FERMI' - ADD 68.20313 'FERMI' - SOS 0......__....000 'FERMI' - NOS **__ Kind Regards On Wed, Apr 16, 2014 at 1:21 AM, Gavin Abo gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu mailto:gs...@crimson.ua.edu wrote: There is usually a more detailed message for 'FERMI - Error' in one of the *lapw2.error files: http://www.mail-archive.com/__wien%40zeus.theochem.tuwien.__ ac.at/msg02361.html http://www.mail-archive.com/wien%40zeus.theochem.tuwien. ac.at/msg02361.html http://www.mail-archive.com/__wien
[Wien] Fermi + Syntax error
Dear Wien2k Users Can you please pull me out of this error? Many times it is appearing while performing Optimization calculations. The wien2k version is 12 and installed on fedora 17 64 bit. Lattice constant I used is 4.397 and it is ternary alloy with 25 % doping of a TM. The error message is ERROR status in optrnew438_vol_-25.0 stop error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error FERMI - Error LAPW1 END LAPW1 END LAPW0 END in cycle 8ETEST: .12554520 CTEST: .0760248 MIXER END CORE END True Regards M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Fermi + Syntax error
Dear Sir Gavin Abo Thank you for you kind response. I have checked these links. What I understood is to increase Emax value. Any other option please? The error file (uplapw2.error) has following message Error in LAPW2 'FERMI' - EFERMI OUT OF ENERGY RANGE 'FERMI' - STOP IN EFI 'FERMI' - ENERGY OF LOWER BOUND : 0.94772 'FERMI' - NUMBER OF STATES AT THE LOWER BOUND : 68.88544 'FERMI' - ENERGY OF UPPER BOUND : 0.94772 'FERMI' - NUMBER OF STATES AT THE UPPER BOUND : 68.87432 'FERMI' - ADD 68.20313 'FERMI' - SOS 0..........000 'FERMI' - NOS ** Kind Regards On Wed, Apr 16, 2014 at 1:21 AM, Gavin Abo gs...@crimson.ua.edu wrote: There is usually a more detailed message for 'FERMI - Error' in one of the *lapw2.error files: http://www.mail-archive.com/wien%40zeus.theochem.tuwien. ac.at/msg02361.html http://www.mail-archive.com/wien%40zeus.theochem.tuwien. ac.at/msg08805.html You checked all error files in the case directory (cat *.error), and they were empty? On 4/15/2014 9:38 AM, Muhammad Sajjad wrote: Dear Wien2k Users Can you please pull me out of this error? Many times it is appearing while performing Optimization calculations. The wien2k version is 12 and installed on fedora 17 64 bit. Lattice constant I used is 4.397 and it is ternary alloy with 25 % doping of a TM. The error message is ERROR status in optrnew438_vol_-25.0 stop error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error (standard_in) 1: syntax error FERMI - Error LAPW1 END LAPW1 END LAPW0 END in cycle 8ETEST: .12554520 CTEST: .0760248 MIXER END CORE END True Regards M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Persisting LOPW and LAPW1 errors
Dear Peter Bala and other wien2k users I am running SCF calculations for Se (structure is attached) with WCGGA and LSDA using Wien2k 12 version. With WCGGA, an error LOPW (plan waves exhausted) appear and I try to solve it as suggested in M. list (like increasing Rkmax,switching to broading method TEMP=0.002, decreasing kpoints) ut the proplem is still appearing. Then I used LSDA and tis time the error arises in LAPW 1 and its statement is Error in LAPW1 'SELECT' - no energy limits found for atom 1 L= 2 'SELECT' - E-bottom -200.0 E-top -200.0 According to UG and M.list suggestions ( decresing Rkmax, switching to broading method (TEMP=0.02), increasing Gmax values, reducing mixing factor to 0.1 and switching to PRATT mixing scheme, increasing difference between Gmax and Gmin) but the same error persists. Kindly suggest me the possible solution. With thanks and regards. M. Sajjad Te.struct Description: Binary data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Persisting LOPW and LAPW1 errors - bug of nn?
Dear Wien2k Users Thank you so much for your helping suggestions. There was a problem with structure. One of my friends make the structure and right now I found that structure is same as attached by Gerhard. Again many thanks. True Regards M.Sajjad On Thu, Feb 6, 2014 at 12:05 AM, Fecher, Gerhard fec...@uni-mainz.dewrote: Nevertheless, the structure is wrong as was already mentioned by Martin, the correct one is attached. The cif file is from the Pearson database and contains standardized crystallographic data and a reference on the origin. Maybe space group 152 was mixed with its enantiomorphous pair that is 154 (to Martin: it is 1/6 in space group 154) Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: I think the problem, to be quite honest with you, is that you have never actually known what the question is. Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-boun...@zeus.theochem.tuwien.ac.at [ wien-boun...@zeus.theochem.tuwien.ac.at]quot; im Auftrag von quot;Laurence Marks [l-ma...@northwestern.edu] Gesendet: Mittwoch, 5. Februar 2014 19:22 Cc: A Mailing list for WIEN2k users Betreff: Re: [Wien] Persisting LOPW and LAPW1 errors - bug of nn? Line 303 of nn.f can be changed to IF((DIST.LT..001).and.(K.eq.index)) GO TO 110 Note: this will now trap many cases with inappropriate precision with 1/3, 5/6 or similar. On Wed, Feb 5, 2014 at 12:06 PM, Laurence Marks l-ma...@northwestern.edu wrote: Correction, I thought I found the bug but was wrong. On Wed, Feb 5, 2014 at 11:24 AM, Laurence Marks l-ma...@northwestern.edu wrote: There is a bug in SRC_nn, somehow associated with the subroutine ORD2 (probably an array that goes out of bounds). I will send a patch to Peter for him to check. On Wed, Feb 5, 2014 at 10:48 AM, Fecher, Gerhard fec...@uni-mainz.de wrote: The Te (prototype Se,hP3,152) structure has the Pearsonsymbol hP3 but the structure file contains 6 atoms, the Wyckoff position is 3a and that has a multiplicity of 3 and not of 6. I guess the z-parameter or anything else is wrong with the structure that is most probably not a bug in the program Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: I think the problem, to be quite honest with you, is that you have never actually known what the question is. Dr. Gerhard H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz Von: wien-boun...@zeus.theochem.tuwien.ac.at [ wien-boun...@zeus.theochem.tuwien.ac.at]quot; im Auftrag von quot;Laurence Marks [l-ma...@northwestern.edu] Gesendet: Mittwoch, 5. Februar 2014 17:00 An: A Mailing list for WIEN2k users Betreff: Re: [Wien] Persisting LOPW and LAPW1 errors - bug of nn? You are right, there seems to be a nn bug. This may explain everything including the LOPW exhausted error. On Wed, Feb 5, 2014 at 8:48 AM, Lyudmila Dobysheva lyuk...@mail.ru wrote: On 05.02.2014 14:29, Muhammad Sajjad wrote: I am running SCF calculations for Se (structure is attached) with WCGGA and LSDA using Wien2k 12 version. With WCGGA, an error LOPW (plan waves exhausted) appear... Then I used LSDA and tis time the error arises in LAPW 1 and its It is rather strange that program cannot find the error at the stage of initialization, because let's see: atoms 1 and 5, 2 and 4, 3 and 6 are very-very close ATOM -1: X=0.2690 Y=0. Z=0.8333 MULT= 6 ISPLIT= 8 -1: X=0. Y=0.2690 Z=0.1666 -1: X=0.7310 Y=0.7310 Z=0.49996667 -1: X=0. Y=0.2690 Z=0.1667 -1: X=0.2690 Y=0. Z=0.83336667 -1: X=0.7310 Y=0.7310 Z=0.5000 I'd say this is almost bug of nn. Best wishes Lyudmila Dobysheva -- Phys.-Techn. Institute of Ural Br. of Russian Ac. of Sci. 426001 Izhevsk, ul.Kirova 132 RUSSIA -- Tel.:7(3412) 218988(office), 722529(Fax) E-mail: l...@ftiudm.ru, lyuk...@mail.ru (office) lyuk...@gmail.com (home) Skype: lyuka17 (home), lyuka18 (office) http://fti.udm.ru/content/view/25/103/lang,english/ -- ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Professor Laurence Marks Department of Materials Science
[Wien] Problem in Optimization of an AFM structure
Dear Users I am running optimization for a ternary half metallic ferromagnet at 75 doping. But I am unable to get the minimum energy point in the parabolic curve. , Gmax is 12, and Rkmax is 7. 100 k-points are given and generated k-points were 12 like this: symmetry operations without inversion inversion added (non-spinpolarized non-so calculation) NUMBER OF K-POINTS IN WHOLE CELL: (0 allows to specify 3 divisions of G) length of reciprocal lattice vectors: 0.832 0.832 0.294 6.564 6.564 2.321 Shift of k-mesh allowed. Do you want to shift: (0=no, 1=shift) 12 k-points generated, ndiv= 6 6 2 KGEN ENDS 0.000u 0.008s 0:00.00 0.0% 0+0k 0+48io 0pf+0w The opt.graph at -5, 0, 5---40 is. [image: Inline image 2] In another calculation I have used Gmax is 16, and Rkmax is 8 but the results same. Also I found the there is an irragular up and down fluctuation in CTEST and ETEST like this [image: Inline image 3] Please suggest me the possible solution. I have flip the proper spins. True regards and thanks M. Sajjad image.pngimage.png___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Convergence problem
Dear users I am running optimization calculation fro a ternary alloy at 25% doping concentration. The lattice constant is 4.87. for -5, 0, 5 values , convergence is not taking place for these values (-5, 0, 5). I have increased and decreased the mixing factor from 0.20 in case.in0 file, also i have increased the enhancing parameter in case.inm file. One more thing I am using only 100 k-points. Anyone please suggest the possible solution. True regards M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Convergence problem
Dear Oleg and Laurence thank you so much for suggestions. I used larger k-mesh and it is working now quite well.yes there were typo mistakes while writing the file names. True Regards M. Sajjad On Tue, Jan 14, 2014 at 12:54 AM, Laurence Marks l-ma...@northwestern.eduwrote: I am not certain I fully agree with the idea of using 0.003 in case.in2 and TEMP. I normally use TEMPS and 0.018 (which is room temperature), TEMPS is much more rigorous than TEMP. That said, without information about why/how it is not converging it is not possible to say much beyond changing the greed (the term mixing factor is obsolete) in case.inm is not normally useful, and the enhancement in case.in0 is more for accuracy than convergence in most cases (there are exceptions). I assume that .in0 and .inm in the original email were flipped as a typo. On Mon, Jan 13, 2014 at 10:43 AM, Oleg Rubel oru...@lakeheadu.ca wrote: In case your host is an insulator, alloying may turn it into a metal (depending in details). Then the convergence on a coarse k-mesh will be bad. As a work around, you can try to add temperature smearing (TEMP option in case.in2) with the value of 0.003 Ry. The obtained total energy should be taken with caution, especially when comparing to another case without smearing. Oleg On Mon, Jan 13, 2014 at 11:31 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear users I am running optimization calculation fro a ternary alloy at 25% doping concentration. The lattice constant is 4.87. for -5, 0, 5 values , convergence is not taking place for these values (-5, 0, 5). I have increased and decreased the mixing factor from 0.20 in case.in0 file, also i have increased the enhancing parameter in case.inm file. One more thing I am using only 100 k-points. Anyone please suggest the possible solution. True regards M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu 1-847-491-3996 Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] NN error
Dear all I am running optimization with wien2k 12 version. I am facing the following error ERROR status in OPTF123_vol__-5.0 stop error NN - Error LAPW0 Kindlyhelp me to solve this propblem. Kind Regards ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] NN error
Dear Tran, Thank you so much for your kind response. Please let me know in which file i can do the mentioned job? I mean how can i choose smaller spheres readii? True Regards Sajjad On Sat, Dec 21, 2013 at 10:16 PM, t...@theochem.tuwien.ac.at wrote: This is because the atomic spheres overlap. So, you have to choose smaller spheres radii. On Sun, 22 Dec 2013, Muhammad Sajjad wrote: Dear all I am running optimization with wien2k 12 version. I am facing the following error ERROR status in OPTF123_vol__-5.0 stop error NN - Error LAPW0 Kindlyhelp me to solve this propblem. Kind Regards ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] NN error
Dear Tran, Thank you so much. It is working now. On Sat, Dec 21, 2013 at 11:30 PM, t...@theochem.tuwien.ac.at wrote: At the beginning of the initiatlization (init_lapw) you are asked to enter a reduction in %. Choose a larger value than you did before. More informations here: http://www.wien2k.at/reg_user/textbooks/usersguide.pdf On Sat, 21 Dec 2013, Muhammad Sajjad wrote: Dear Tran, Thank you so much for your kind response. Please let me know in which file i can do the mentioned job? I mean how can i choose smaller spheres readii? True Regards Sajjad On Sat, Dec 21, 2013 at 10:16 PM, t...@theochem.tuwien.ac.at wrote: This is because the atomic spheres overlap. So, you have to choose smaller spheres radii. On Sun, 22 Dec 2013, Muhammad Sajjad wrote: Dear all I am running optimization with wien2k 12 version. I am facing the following error ERROR status in OPTF123_vol__-5.0 stop error NN - Error LAPW0 Kindlyhelp me to solve this propblem. Kind Regards ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/ wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] About mBJ convergence
Thank you Laurence, it is working now. I am dealing two cases and both of your suggestions working. Regards On Mon, Sep 9, 2013 at 2:27 AM, Laurence Marks l-ma...@northwestern.eduwrote: With no information nobody can tell you anything useful. On Sun, Sep 8, 2013 at 12:09 PM, Muhammad Sajjad sajja...@gmail.com wrote: I have done the first job with mixing factor 0.1 and PRATT then switched to MSRI, it also takes about 160 iterations. I will try your second suggested job. Can you tell me please that after how many iterations, the convergence is considered to be good? Thanks and True regards On Mon, Sep 9, 2013 at 12:47 AM, Laurence Marks l-ma...@northwestern.edu wrote: With the older mixer, to start mBJ going it was best to run some cycles with 0.1 and PRATT first, then switch to MSR1. With the current mixer this should be better, perhaps sometimes you might also want to start with a few (5-10) PRATT steps. In general MSR1 converges much better. On Sun, Sep 8, 2013 at 11:44 AM, Muhammad Sajjad sajja...@gmail.com wrote: Dear Okoye Thank you for your suggestion, I have tried this already and repeated it again with mixing factor 0.3, but the problem is still there. Any other suggestion Okoye and other Wien2k users please? Regards On Sat, Sep 7, 2013 at 9:36 PM, okoye...@yahoo.com wrote: I think you should increase the mixing factor after some initial scf cycles. This is contained in the UG. Regards C M I Okoye Sent from my BlackBerry wireless device from MTN -Original Message- From: Muhammad Sajjad sajja...@gmail.com Sender: wien-boun...@zeus.theochem.tuwien.ac.at Date: Sat, 7 Sep 2013 21:29:42 To: wienWien@zeus.theochem.tuwien.ac.at Reply-To: A Mailing list for WIEN2k users wien@zeus.theochem.tuwien.ac.at Subject: [Wien] About mBJ convergence ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu 1-847-491-3996 Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu 1-847-491-3996 Research is to see what everybody else has seen, and to think what nobody else has thought Albert Szent-Gyorgi ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] About mBJ convergence
Dear Users I have question about mBJ calculations. When we execute scf calculations with LSDA (or other functionals) then mostly over calculation finishes with in 30 iterations. But when we apply mBJ over it, then it finishes about 250 iterations. Is this going fine? I mean results obtained are correct? I am doing this on Wien2k12.1 version, and running the mBJ calculations with the method given in user guide. In case.inm file I switches MRSI to PRAAT and mixing factor to 0.1? my case is spin-polarized and sometimes with in 60 iterations the charge is converged. please guide me in this regard. With thanks M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Error in mBJGGA
Dear Saeid Jalali Thank you so much, it is really helping. True Regards M Sajjad On Tue, Sep 3, 2013 at 8:16 PM, Saeid Jalali s_jalal...@yahoo.com wrote: Dear Sajjad, I made a mistake in case.in0_grr by selecting indxc value 55 instead of 50, and I am repeating this mistake for two days :( In order to avoid such a mistake, you can use init_mbj_lapw script twice. This script automatically (in its second run) edits case.in0 and sets indxc to 28 and copes case.in0 to case.ino_grr and sets indxc to 50 in case.in0_grr. For more information see 4.5.9 section of the users guide. Sincerely yours, S. Jalali *From:* Muhammad Sajjad sajja...@gmail.com *To:* A Mailing list for WIEN2k users wien@zeus.theochem.tuwien.ac.at *Sent:* Tuesday, September 3, 2013 1:15 PM *Subject:* Re: [Wien] Error in mBJGGA Dear F. Tran Thank you for correction. I made a mistake in case.in0_grr by selecting indxc value 55 instead of 50, and I am repeating this mistake for two days :( True Regards M. Sajjad On Tue, Sep 3, 2013 at 3:57 PM, t...@theochem.tuwien.ac.at wrote: Hi, you have probably not selected the correct value for indxc in case.in0 or case.in0_grr. It should be 28 in case.in0 and 50 in case.in0_grr. I guess that you did it correctly for 0.25 doping. F. Tran On Tue, 3 Sep 2013, Muhammad Sajjad wrote: Dear Wien2k users I am am running mBJGGA calculations for ternary alloy. the super cell is of 8 atoms and doping is 0.75. For 0.25 doping, the mBJGGA calculations (spin is involved) have completed with no error, but for 0.75 doping the following error is appearing. [msajjad@msajjad SCF75]$ runsp_lapw -cc 0.1 -in1new 2 -i 100 -NI forrtl: severe (24): end-of-file during read, unit 60, file /home/msajjad/3rdpaper/MgVTe/**SCF75/SCF75.inhf Image PCRoutineLineSource lapw0 005405AD Unknown Unknown Unknown lapw0 0053F0B5 Unknown Unknown Unknown lapw0 004DF760 Unknown Unknown Unknown lapw0 0049E7AA Unknown Unknown Unknown lapw0 0049DFA0 Unknown Unknown Unknown lapw0 004BDE9C Unknown Unknown Unknown lapw0 00441273 MAIN__255 lapw0.F lapw0 00403BAC Unknown Unknown Unknown libc.so.6 0034BF01ECDD Unknown Unknown Unknown lapw0 00403AA9 Unknown Unknown Unknown stop error Please help me to overcome this problem. With thanks M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Error in mBJGGA
Dear Wien2k users I am am running mBJGGA calculations for ternary alloy. the super cell is of 8 atoms and doping is 0.75. For 0.25 doping, the mBJGGA calculations (spin is involved) have completed with no error, but for 0.75 doping the following error is appearing. [msajjad@msajjad SCF75]$ runsp_lapw -cc 0.1 -in1new 2 -i 100 -NI forrtl: severe (24): end-of-file during read, unit 60, file /home/msajjad/3rdpaper/MgVTe/SCF75/SCF75.inhf Image PCRoutineLineSource lapw0 005405AD Unknown Unknown Unknown lapw0 0053F0B5 Unknown Unknown Unknown lapw0 004DF760 Unknown Unknown Unknown lapw0 0049E7AA Unknown Unknown Unknown lapw0 0049DFA0 Unknown Unknown Unknown lapw0 004BDE9C Unknown Unknown Unknown lapw0 00441273 MAIN__255 lapw0.F lapw0 00403BAC Unknown Unknown Unknown libc.so.6 0034BF01ECDD Unknown Unknown Unknown lapw0 00403AA9 Unknown Unknown Unknown stop error Please help me to overcome this problem. With thanks M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Error in mBJGGA
Dear F. Tran Thank you for correction. I made a mistake in case.in0_grr by selecting indxc value 55 instead of 50, and I am repeating this mistake for two days :( True Regards M. Sajjad On Tue, Sep 3, 2013 at 3:57 PM, t...@theochem.tuwien.ac.at wrote: Hi, you have probably not selected the correct value for indxc in case.in0 or case.in0_grr. It should be 28 in case.in0 and 50 in case.in0_grr. I guess that you did it correctly for 0.25 doping. F. Tran On Tue, 3 Sep 2013, Muhammad Sajjad wrote: Dear Wien2k users I am am running mBJGGA calculations for ternary alloy. the super cell is of 8 atoms and doping is 0.75. For 0.25 doping, the mBJGGA calculations (spin is involved) have completed with no error, but for 0.75 doping the following error is appearing. [msajjad@msajjad SCF75]$ runsp_lapw -cc 0.1 -in1new 2 -i 100 -NI forrtl: severe (24): end-of-file during read, unit 60, file /home/msajjad/3rdpaper/MgVTe/**SCF75/SCF75.inhf Image PCRoutineLine Source lapw0 005405AD Unknown Unknown Unknown lapw0 0053F0B5 Unknown Unknown Unknown lapw0 004DF760 Unknown Unknown Unknown lapw0 0049E7AA Unknown Unknown Unknown lapw0 0049DFA0 Unknown Unknown Unknown lapw0 004BDE9C Unknown Unknown Unknown lapw0 00441273 MAIN__255 lapw0.F lapw0 00403BAC Unknown Unknown Unknown libc.so.6 0034BF01ECDD Unknown Unknown Unknown lapw0 00403AA9 Unknown Unknown Unknown stop error Please help me to overcome this problem. With thanks M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] How to check K-points parallel installation
Dear Jia Yalei Thank you for your help. Regards M. Sajjad Beijing, China On Thu, Aug 1, 2013 at 9:31 AM, 贾亚磊 jia_ya...@163.com wrote: Dear Sajjad, Though I have not tested the parallel installation successfully, in my opinion, the answer is in *http://www.wien2k.at/reg_user/faq/pbs.html*. Change the example shell script there(or in WIENROO/qsub-job1-*_lapw) according to you queuing system. The submit procedure is: 1. Do the init_lapw in your working directory [ or through w2web if you do not have *.struct]. 2. Copy your changed shell script to your working directory. 3. Submit you shell script to your queuing system. Best wishes Jia Yalei ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] How to check K-points parallel installation
Dear Wien2k Group Members I am using Wien2k 12 with k-points parallel installation. Anyone please can let me know that how to confirm k-points parallel installation by giving any command in terminal? Regards M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] how to converge charge in mBJ
Dear Wien2k Users I am running mBJ calculations. Unfortunately every time charge is not converging. I have used 0.5, 0.2, 0.1, 0.05 and 0.02 mixing factor values in case.inm file. Also I used the command *runsp_lapw -cc 0.0001 -in1new 2 -i 100 -NI * for all attempts. What should I do to overcome this proplem? With Thanks and Regard M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] QTL -B error in mBJ
Dear Wien2k Users I am trying to run mBJ for a ternary alloy. everything was right but every time QTL -B error appears. the error statement is *[msajjad@msajjad SCF]$ runsp_lapw -cc 0.0001 -in1new 2 -i 100 -NI* * LAPW0 END* * LAPW0 END* * LAPW1 END* * LAPW1 END* *L2main - QTL-B Error* * * Please someone let me know the solution to this problem? With thanks in advance M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] QTL -B error in mBJ
Dear Users I have found its solution from the mailing list, and up to now the calculations are running error-less. Thank you On Wed, May 29, 2013 at 10:48 PM, Muhammad Sajjad sajja...@gmail.comwrote: Dear Wien2k Users I am trying to run mBJ for a ternary alloy. everything was right but every time QTL -B error appears. the error statement is *[msajjad@msajjad SCF]$ runsp_lapw -cc 0.0001 -in1new 2 -i 100 -NI* * LAPW0 END* * LAPW0 END* * LAPW1 END* * LAPW1 END* *L2main - QTL-B Error* * * Please someone let me know the solution to this problem? With thanks in advance M. Sajjad ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html