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 Abowrote: > 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. > > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: >
Re: [Wien] Si under strain
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. 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<mailto: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 t
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
what is the space group of the diamond structure ? 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 12:51 An: A Mailing list for WIEN2k users Betreff: 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<mailto: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.496943
Re: [Wien] Si under strain
It is 216_F-43m On Thu, Sep 10, 2015 at 1:51 PM, Muhammad Sajjadwrote: > 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 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
Re: [Wien] Si under strain
http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg02749.html On 9/10/2015 7:33 AM, Muhammad Sajjad wrote: 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? ___ 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
Gavin you are right, I still use it as an example on how to make Si a metal. Seems nowadys symmetry is not teached anymore :-(( 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 Gavin Abo [gs...@crimson.ua.edu] Gesendet: Donnerstag, 10. September 2015 16:21 An: A Mailing list for WIEN2k users Betreff: Re: [Wien] Si under strain http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg02749.html On 9/10/2015 7:33 AM, Muhammad Sajjad wrote: > 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? ___ 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] Si under strain
and your strained lattice with (a =5.49 , b= 5.43 , c = 5.38) is tetragonal and also not F -43m which is cubic, isn't it ? I guess you should find something in space group I 4_1/AMD (141), but it's your work to figure out how and why. 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 15:33 An: A Mailing list for WIEN2k users Betreff: 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<mailto: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<mailto:wien-boun...@zeus.theochem.tuwien.ac.at> [wien-boun...@zeus.theochem.tuwien.ac.at<mailto:wien-boun...@zeus.theochem.tuwien.ac.at>] im Auftrag von Muhammad Sajjad [sajja...@gmail.com<mailto: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><mailto: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
[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] Si under strain
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.4969432792132343 0.4969432792132343 0.49694327921323430. 0.4969432792132343 0.49694327921323430.4969432792132343 0. 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. ___ 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