Thanks a lot for your reply!<br><br>I've tried for another day, I know it sounds stupid, but I still have something unclear, could you please still help me with that? Sorry for keeping asking such entrance level questions for such a long
time, but I really want to know more, thanks a lot for your help
again!!!<br><br>First, I checked ICSD, yes the bismuth space group is 166 R-3m, however, on American Mineralogist Crystal Structure Database (AMD <a href="http://rruff.geo.arizona.edu/AMS/amcsd.php">http://rruff.geo.arizona.edu/AMS/amcsd.php
</a> ) if I search "bismuth", then I can see some papers show its space group is 166 R-3m, while others show its space group is 160 R3m (h axes). <br><br>Then I downloaded the CIF file from ICSD, it is exactly that one showed in Dr Blaha' last email (6 equivalent atoms, h axes). And I used CIF2STRUCT on this CIF file, it generated a struct file almost the same as in Dr Blaha's last email, but GAMMA angle is 120 instead of 90. I don't know if there is any problem with that.
<br><br>Second, I downloaded CIF file from AMD which says the space group of Bismuth is 160 R3m , I examined that CIF file and find that it is in Hex axes which complies with the required input of CIF2STRUCT, so I run CIF2STRUCT, and a struct like below was generated:
<br><br>==============================<br>blebleble<br>H LATTICE,NONEQUIV.ATOMS 1 160 R3m<br>MODE OF CALC=RELA unit=bohr<br> 8.590695 8.590695 22.415931 90.000000 90.000000120.000000<br>ATOM -1: X=0.00000000 Y=0.00000000
Z=0.23389000<br> MULT= 3 ISPLIT=15<br> -1: X=0.66666667 Y=0.33333333 Z=0.56722333<br> -1: X=0.33333333 Y=0.66666667 Z=0.90055667<br>Bi NPT= 781 R0=.000005000 RMT= 2.00000 Z: 83.00000
<br>LOCAL ROT MATRIX: 1.0000000 0.0000000 0.0000000<br> 0.0000000 1.0000000 0.0000000<br> 0.0000000 0.0000000 1.0000000<br> 0 NUMBER OF SYMMETRY OPERATIONS<br>=================================
<br>Then I continued with Init_lapw, (NN distance is 8.59070 in this case ) when I run sgroup, it told me 166 R-3m symmetry was found and the Bravais lattice has changed. I got a question here: according to the CIF file, the atomic positions in this case is due to the symmetry of 160 R3m space group (in hex settings), but why after "sgroup" it told me 166 R-3m symmetry was found?
<br><br>Last question:<br><br>A CIF file which contains the file like below can still work with CIF2STRUCT directly, although it is in rhombohedral axes, am I right? Thanks a lot!<br>======================<br>Bismuth<br>Wyckoff R W G
<br>Crystal Structures 1 (1963) 7-83<br>Second edition. Interscience Publishers, New York, New York<br>4.7459 4.7459 4.7459 57.237 57.237 57.237 R-3m<br>atom x y z<br>Bi .237 .237 .237<br>==========================
<br><br>Thanks and Best Wishes<br><br>Sincerely yours,<br>Bo Qiu<br><br><div><span class="gmail_quote">On 10/25/07, <b class="gmail_sendername">Peter Blaha</b> <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at
</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">First: You need to learn some basic crystallography.<br>Some spacegroups have more than one setting. Click on the link in the
<br>STRUCTEDITOR.<br><br>Second: Bi is R-3m (number 166) and not 160 R3m !!!<br><br>When you take a cif file which gives the structure in hex<br>cit2struct converts it into a hex lattice:<br><br>blebleble<br><br>H LATTICE,
NONEQUIV.ATOMS 1 166 R-3m<br>MODE OF CALC=RELA unit=bohr<br> 8.569908 8.569908 22.325225 90.000000 90.000000120.000000<br>ATOM -1: X=0.00000000 Y=0.00000000 Z=0.76600000<br> MULT= 6 ISPLIT=15<br>
-1: X=0.00000000 Y=0.00000000 Z=0.23400000<br> -1: X=0.66666667 Y=0.33333333 Z=0.09933333<br> -1: X=0.66666667 Y=0.33333333 Z=0.56733333<br> -1: X=0.33333333 Y=0.66666667 Z=0.43266667<br> -1: X=
0.33333333 Y=0.66666667 Z=0.90066667<br>Bi 0+ NPT= 781 R0=.000005000 RMT= 2.00000 Z: 83.00000<br><br>So there are 6 equivalent positions, not 3 !<br><br>This struct file does not obey the WIEN2k rules for R lattices, so you
<br>MUST NOT save it in STRUCTEDITOR.<br><br>When you run init_lapw with this struct file, sgroup will convert it to<br>a WIEN2k conform structure and create a file like:<br><br>Bismut<br>R LATTICE,NONEQUIV. ATOMS 1<br>
MODE OF CALC=RELA<br> 8.590000 8.590000 22.415000 90.000000 90.000000 90.000000<br>ATOM -1: X=0.23400000 Y=0.23400000 Z=0.23400000<br> MULT= 2 ISPLIT= 4<br>ATOM -1: X=0.76600000 Y=0.76600000 Z=0.76600000
<br>Bi NPT= 781 R0=0.00001000 RMT= 2.8000 Z: 83.0<br><br>You can also find this file in $WIENROOT/example_struct_files.<br><br><br>Bo Qiu schrieb:<br>> Thanks for your instructions!<br>><br>> However, I did have specified the space group to be R3m. Actually what I
<br>> have done is to use hex2rhomb to give this transformation<br>> (0,0,0.23407)-->(0.23407,0.23407,0.23407), once I get this internal<br>> parameters, I go back to w2web and use STRUCTGEN to make the structure
<br>> file, and I specified the space group to be 160_R3m, however I just get<br>> one atom, and the successive initializations failed...<br>><br>> Then I try to download the CIF file from American Mineralogist Crystal
<br>> Structure database, then use cif2struct to generate the file. Then I go<br>> back to w2web and take a look at the structure file, the internal<br>> parameters are:<br>><br>> ATOM -1: X=0.00000000 Y=0.00000000
Z=0.23407000<br>> MULT= 3 ISPLIT=15<br>> -1: X=0.66666667 Y=0.33333333 Z=0.56740333<br>> -1: X=0.33333333 Y=0.66666667 Z=0.90073667<br>><br>> they looked good. However, once I pressed "edit struct file" then "save
<br>> structure" without any change, I found the structure file has been<br>> modified, again, to be:<br>><br>> ATOM -1: X=0.99999999 Y=0.00000000 Z=0.23407000<br>> MULT= 3 ISPLIT= 5
<br>> ATOM -1:X= 0.00000000 Y=0.23407000 Z=0.99999999<br>> ATOM -1:X= 0.23407000 Y=0.99999999 Z=0.00000000<br>><br>> And I run init_lapw, it still failed... I don't really know why.<br>><br>> So later I tried just use the structure file generated by cif2struct and
<br>> run init_lapw, but when it gets to sgroup, it told me that the space<br>> group found was R-3m instead of R3m (while R3m is the specified one in<br>> that cif file). I accepted that and afterwards there is core electron
<br>> leakage when setting separation energy to be -6Ry, then I changed it to<br>> -10Ry. Seemed it worked. And SCF seemed to be working. Then I proceeded<br>> to generate the band structure plotting, and compare it with
<br>> literature... but it seems to have not much similarity...<br>><br>> Could you please help me with that? I really don't know what's wrong<br>> with my setup... could you please tell me how to do that? Thank you very
<br>> very much!!!!<br>><br>> Sincerely yours,<br>><br>> Bo Qiu<br>><br>><br>><br>> On 10/25/07, *Peter Blaha* <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a><br>
> <mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>>> wrote:<br>><br>> But you still need to specify the space group.<br>> Only by knowing the spacegroup, STRUCTGEN can put the second position
<br>> into your struct file.<br>><br>> Bo Qiu schrieb:<br>> > Thank you very much for your reply!!! It really helps me a lot!<br>> ><br>> > However, I still have some problem with that...
<br>> ><br>> > Due to your instructions, I use hex2rhomb to convert the<br>> coordinates:<br>> > (0,0,0.23407)-->(0.23407,0.23407,0.23407)<br>> ><br>> > Then I modified the input:
<br>> ><br>> > lattice parameters:<br>> > a=b=4.5330 A, c=11.797 A, alpha=belta=90 degree, gamma=120 degree.<br>> > Atom 1: Bi , Z=83, RMT= 0.99<br>> > x=0.23407, y=
0.23407, z=0.23407.<br>> ><br>> > After I "save and clean up", I only got one atom... and when I<br>> continued<br>> > with initialization, when it reached the step "sgroup", it reported
<br>> > inconsistency between the space group found and the struct<br>> file... and<br>> > more errors to come...<br>> ><br>> > Could you please tell me what is still wrong in my setup? Thanks
<br>> a lot!<br>> ><br>> > Best Wishes,<br>> ><br>> > Sincerely yours,<br>> > Bo Qiu<br>> ><br>> ><br>> ><br>> > On 10/24/07, *Rocquefelte* <
<a href="mailto:Xavier.Rocquefelte@cnrs-imn.fr">Xavier.Rocquefelte@cnrs-imn.fr</a><br>> <mailto:<a href="mailto:Xavier.Rocquefelte@cnrs-imn.fr">Xavier.Rocquefelte@cnrs-imn.fr</a>><br>> > <mailto:
<a href="mailto:Xavier.Rocquefelte@cnrs-imn.fr">Xavier.Rocquefelte@cnrs-imn.fr</a><br>> <mailto:<a href="mailto:Xavier.Rocquefelte@cnrs-imn.fr">Xavier.Rocquefelte@cnrs-imn.fr</a>>>> wrote:<br>> >
<br>> > Dear Bo Qiu,<br>> ><br>> > In the particular case of rhombohedral structures the<br>> initialization of<br>> > the Wien2k calculation should be done by specifying the
<br>> hexagonal<br>> > lattice parameters and the rhombohedral coordinates in the<br>> structure<br>> > file. An auxiliary program can be used to convert your<br>> coordinates from
<br>> > hexagonal to rhombohedral (hex2rhomb). You will find more<br>> details about<br>> > this question in the userguide (p. 40).<br>> ><br>> > Have a nice day
<br>> ><br>> > Regards<br>> ><br>> > Xavier<br>> ><br>> ><br>> > Bo Qiu a écrit :<br>> > > Dear Wien2k users and developers,
<br>> > ><br>> > > Recently I'm trying to calculate Bismuth, using the data<br>> got from<br>> > > American Mineralogist Crystal Structure database, however,
<br>> I'm not<br>> > > sure about the structrue generated by wien2k, and what's<br>> more, based<br>> > > on that structure generated, calculation failed even at the
<br>> > > initialization step. Could you please help me with that?<br>> Thanks a<br>> > lot!!!<br>> > ><br>> > > Below is the data from American Mineralogist Crystal
<br>> Structure<br>> > database<br>> > > ===============================================<br>> > > Bismuth<br>> > ><br>> > > Schiferl D, Barrett CS
<br>> > > Journal of Applied Crystallography 2 (1969) 30-36<br>> > > The crystal structure of arsenic at 4.2, 78 and 299 K<br>> > > Sample: T = 4.2 K<br>> > >
<br>> > > 4.5330 4.5330 11.797 90 90 120 R3m<br>> > > atom x y z<br>> > > Bi 0 0 .23407<br>> > > ==============================================
<br>> > ><br>> > > Then based on this data, I start a session with 1<br>> inequivalent atom,<br>> > > space group R3m,<br>> > > lattice parameters:
<br>> > > a=b=4.5330 A, c=11.797 A, alpha=belta=90 degree, gamma=120<br>> degree.<br>> > > Atom 1: Bi , Z=83, RMT=0.99<br>> > > x=0, y=0, z=0.23407.<br>> > >
<br>> > > However after I saved the structure, I got this structrue<br>> file:<br>> > ><br>> > > ==================================================<br>> > > bi
<br>> > > R LATTICE,NONEQUIV.ATOMS: 1160_R3m<br>> > > MODE OF CALC=RELA unit=ang<br>> > > 8.566510 8.566510 22.293108 90.000000 90.000000120.000000<br>> > > ATOM -1: X=
0.99999999 Y=0.00000000 Z=0.23407000<br>> > > MULT= 3 ISPLIT= 8<br>> > > ATOM -1:X= 0.00000000 Y=0.23407000 Z=0.99999999<br>> > > ATOM -1:X= 0.23407000
Y=0.99999999 Z=0.00000000<br>> > > Bi NPT= 781 R0= 0.00000500 RMT= 1.1000 Z: 83.0<br>> > > LOCAL ROT MATRIX: 0.0000000 1.0000000 0.0000000<br>> > >
0.0000000 0.0000000 1.0000000<br>> > > 1.0000000 0.0000000 0.0000000<br>> > > 6 NUMBER OF SYMMETRY OPERATIONS<br>> > > ......<br>> > >
<br>> > > ==================================================<br>> > > I'm not quite sure, but I think the coordinates for atom 1,2,3<br>> > might be:<br>> > > (0,0,
0.23407 ),<br>> (2/3,1/3,1/3+0.23407),(1/3,2/3,2/3+0.23407) due to<br>> > > space group symmetry which is different from that one in the<br>> > structure<br>> > > file. So why is the difference? Thanks a lot!
<br>> > ><br>> > > Though I don't understand the structure file, I still<br>> proceeded to<br>> > > initialization. At the first step it seemed to be ok, however,
<br>> > when it<br>> > > reached "x lstart", there is warning about core electron<br>> leakage.<br>> > So I<br>> > > go back to increase RMT, however, after that, NN error
<br>> showed up<br>> > > during initialization... then I tried to keep RMT but to<br>> reduce<br>> > > separation energy, it failed when this energy is set to be<br>> less than
<br>> > > -10 Ry...<br>> > ><br>> > > That's the problem I met. Could you please help me with that?<br>> > Thanks a<br>> > > lot!
<br>> > ><br>> > > The system I'm using is<br>> > > Ifort, mkl 9.0, bash shell, wien2k_07<br>> > ><br>> > > --<br>> > > Thanks and Best Wishes.
<br>> > ><br>> > > Sincerely yours,<br>> > > Bo Qiu<br>> > > <a href="http://web.ics.purdue.edu/~qiub/">http://web.ics.purdue.edu/~qiub/</a><br>> <
<a href="http://web.ics.purdue.edu/%7Eqiub/">http://web.ics.purdue.edu/%7Eqiub/</a>><br>> > ><br>> ><br>> ------------------------------------------------------------------------<br>
><br>> > ><br>> > > _______________________________________________<br>> > > Wien mailing list<br>> > > <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">
Wien@zeus.theochem.tuwien.ac.at</a><br>> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>> > <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">
Wien@zeus.theochem.tuwien.ac.at</a><br>> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>>><br>> > > <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>> <<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>><br>> > >
<br>> ><br>> > _______________________________________________<br>> > Wien mailing list<br>> > <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at
</a><br>> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>> <mailto:<br>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at
</a><br>> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>>><br>> > <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
</a><br>> <<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>><br>> ><br>> ><br>> ><br>> >
<br>> > --<br>> > Thanks and Best Wishes.<br>> ><br>> > Sincerely yours,<br>> > Bo Qiu<br>> > <a href="http://web.ics.purdue.edu/~qiub/">http://web.ics.purdue.edu/~qiub/
</a><br>> ><br>> ><br>> ><br>> ------------------------------------------------------------------------<br>> ><br>> > _______________________________________________
<br>> > Wien mailing list<br>> > <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at
</a>><br>> > <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>> _______________________________________________<br>> Wien mailing list
<br>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>> <<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>><br>><br>>
<br>><br>><br>> --<br>> Thanks and Best Wishes.<br>><br>> Sincerely yours,<br>> Bo Qiu<br>> <a href="http://web.ics.purdue.edu/~qiub/">http://web.ics.purdue.edu/~qiub/</a><br>><br>><br>> ------------------------------------------------------------------------
<br>><br>> _______________________________________________<br>> Wien mailing list<br>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>_______________________________________________<br>Wien mailing list<br><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br></blockquote></div><br><br clear="all"><br>-- <br>Thanks and Best Wishes.<br><br>Sincerely yours,
<br>Bo Qiu<br><a href="http://web.ics.purdue.edu/~qiub/">http://web.ics.purdue.edu/~qiub/</a>