<div dir="ltr"><div><div>Dear Prof. peter Blaha<br></div>Thanks for your reply. I want to calculate U from constrained DFT and for that I am going through the following example:<br><br><a href="http://www.wien2k.at/reg_user/textbooks/Constraint_U.pdf">http://www.wien2k.at/reg_user/textbooks/Constraint_U.pdf</a><br>
<br></div>First I have made the structure of Nio.struct which is following:<br>NiO<br>F LATTICE,NONEQUIV.ATOMS: 2225_Fm-3m<br>MODE OF CALC=RELA unit=bohr<br> 7.927000 7.927000 7.927000 90.000000 90.000000 90.000000<br>
ATOM -1: X=0.00000000 Y=0.00000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>ATOM 2: X=0.50000000 Y=0.50000000 Z=0.50000000<br> MULT= 1 ISPLIT= 8<br>O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0 file is given below:<br>
LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br> 0 NUMBER OF SYMMETRY OPERATIONS<br><br><div><div><br>
</div><div>Next I have run the supercell for this strcuture and generated Nio_super.struct file with 2*2*2 and 0 shift in fractional coordianate.In the NiO_super.struc example file, there are 5 non-equivalant atom but here I am getting 16 non equivalent atom. I have deleted all the atom except first Ni atom and rename it with Ni1 in the super.struc file. I have run the init_lapw with this super.struc file and getting the error. My Nio_super.struct file is given below:<br>
NiO<br>F LATTICE,NONEQUIV. ATOMS 16<br>MODE OF CALC=RELA unit=bohr<br> 15.854000 15.854000 15.854000 90.000000 90.000000 90.000000<br>ATOM 1: X=0.00000000 Y=0.00000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>
Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>
ATOM 2: X=0.50000000 Y=0.00000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>ATOM 3: X=0.25000000 Y=0.25000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>ATOM 4: X=0.75000000 Y=0.25000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>ATOM 5: X=0.25000000 Y=0.00000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>
Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>
ATOM 6: X=0.75000000 Y=0.00000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>ATOM 7: X=0.00000000 Y=0.25000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>ATOM 8: X=0.50000000 Y=0.25000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>Ni NPT= 781 R0=0.00005000 RMT= 2.1200 Z: 28.0<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>ATOM 9: X=0.25000000 Y=0.25000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>
O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br>
ATOM 10: X=0.75000000 Y=0.25000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br>
0.0000000 0.0000000 0.0000000<br>ATOM 11: X=0.00000000 Y=0.00000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br>
0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br>ATOM 12: X=0.50000000 Y=0.00000000 Z=0.25000000<br> MULT= 1 ISPLIT= 8<br>O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>
LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br>ATOM 13: X=0.00000000 Y=0.25000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>
O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br>
ATOM 14: X=0.50000000 Y=0.25000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br>
0.0000000 0.0000000 0.0000000<br>ATOM 15: X=0.25000000 Y=0.00000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br>
0.0000000 0.0000000 0.0000000<br> 0.0000000 0.0000000 0.0000000<br>ATOM 16: X=0.75000000 Y=0.00000000 Z=0.00000000<br> MULT= 1 ISPLIT= 8<br>O NPT= 781 R0=0.00010000 RMT= 1.8200 Z: 8.0<br>
<br>0 NUMBER OF SYMMETRY OPERATIONS<br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Jul 15, 2013 at 3:11 PM, Peter Blaha <span dir="ltr"><<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Why did you "shift" all atoms by 2 ???<br>
Don't shift at all. (0 0 0)<br>
<br>
What did you do in the vi ... step ???<br>
In principle you would have to assign to each Ni whether it is spin-up or dn ....<br>
<br>
I don't know what you want to do.<br>
<br>
If you want to simulate the AFM2 magnetic structure of NiO, use the CoO example in the structfile_examples in $WIENROOT.<br>
Change to Ni, adjust the lattice parameters to NiO (they are a/sqrt(2) and a*sqrt(3) and set a smaller R0 value (one more 0; use replace, not insert !!)<div class="im"><br>
<br>
<br>
Am 15.07.2013 08:55, schrieb wasim raja Mondal:<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear wien2k experts<div><div class="h5"><br>
<br>
I am doing constrained DFT calculation of NiO following discussion by<br>
laurence marks in the wien2k mailing list. The following steps I have done:<br>
<br>
(1) I have created usual NiO.struct file using crystal editor which is<br>
following:<br>
<br>
NiO<br>
F LATTICE,NONEQUIV.ATOMS: 2225_Fm-3m<br>
MODE OF CALC=RELA unit=bohr<br>
7.927000 7.927000 7.927000 90.000000 90.000000 90.000000<br>
ATOM -1: X=0.00000000 Y=0.00000000 Z=0.00000000<br>
MULT= 1 ISPLIT= 8<br>
Ni1 NPT= 781 R0=0.00005000 RMT= 2.3000 Z: 28.0<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>
ATOM 2: X=0.50000000 Y=0.50000000 Z=0.50000000<br>
MULT= 1 ISPLIT= 8<br>
O 2 NPT= 781 R0=0.00010000 RMT= 1.6500 Z: 8.0<br>
LOCAL ROT MATRIX: 0.0000000 0.0000000 0.0000000<br>
0.0000000 0.0000000 0.0000000<br>
0.0000000 0.0000000 0.0000000<br>
0 NUMBER OF SYMMETRY OPERATIONS<br>
<br>
<br>
(2) Next I have generated supercell program and getting following errors:<br>
<br>
<br>
Program generates supercell from a WIEN struct file.<br>
<br>
Filename of struct file:<br>
NiO.struct<br>
<br>
Number of cells in x direction:<br>
2<br>
Number of cells in y direction:<br>
2<br>
Number of cells in z direction:<br>
2<br>
Optional shift all atoms by the same amount (fractional coordinates).<br>
Please enter x shift:<br>
2<br>
Please enter y shift:<br>
2<br>
Please enter z shift:<br>
2<br>
<br>
Current structure has lattice type F<br>
Enter your target lattice type: (P,B,F)<br>
F<br>
Target lattice type will be F<br>
<br>
Supercell generated sucessfully.<br>
Stored in struct file: NiO_super.struct<br>
<br>
You may need to replace an atom by an impurity or distort the<br>
positions, ....<br>
raja@viper5:~/BaRuO3_wien2k/<u></u>constr_DFT/NiO$ ls<br>
fort.21 :log NiO.struct NiO.struct_ii NiO_super.struct nn.def<br>
setrmt.bva setrmt.nnshells setrmt.outputnn setrmt.struct<br>
setrmt.struct_nn setrmt.struct_setrmt<br>
raja@viper5:~/BaRuO3_wien2k/<u></u>constr_DFT/NiO$ vi NiO.struct<br>
raja@viper5:~/BaRuO3_wien2k/<u></u>constr_DFT/NiO$ init_lapw<br>
next is setrmt<br>
Automatic determination of RMTs. Please specify the desired RMT reduction<br>
compared to almost touching spheres.<br>
Typically, for a single calculation just hit enter, for force minimization<br>
use 1-5; for volume effects you may need even larger reductions.<br>
<br>
Enter reduction in %<br>
0<br>
Use old or new scheme (o/N)<br>
N<br>
specify nn-bondlength factor: (usually=2) [and optionally dlimit,<br>
dstmax (about<br>
1.d-5, 20)]<br>
DSTMAX: 20.0000000000000<br>
iix,iiy,iiz 5 5 5 39.6350000000000<br>
39.6350000000000 39.6350000000000<br>
NAMED ATOM: Ni1 Z changed to IATNR+999 to determine equivalency<br>
forrtl: severe (24): end-of-file during read, unit 20, file<br>
/home/raja/BaRuO3_wien2k/<u></u>constr_DFT/NiO/NiO.struct<br>
Image PC Routine Line Source<br>
nn 00000000004955EE Unknown Unknown Unknown<br>
nn 0000000000494086 Unknown Unknown Unknown<br>
nn 000000000044B462 Unknown Unknown Unknown<br>
nn 000000000040F1BC Unknown Unknown Unknown<br>
nn 000000000040E6DC Unknown Unknown Unknown<br>
nn 000000000042F7B5 Unknown Unknown Unknown<br>
nn 000000000042D648 Unknown Unknown Unknown<br>
nn 000000000040597A MAIN__ 218 nn.f<br>
nn 000000000040349C Unknown Unknown Unknown<br>
libc.so.6 00002B233B4AA76D Unknown Unknown Unknown<br>
nn 0000000000403399 Unknown Unknown Unknown<br>
0.0u 0.0s 0:00.00 0.0% 0+0k 0+16io 0pf+0w<br>
error: command /home/raja/wien2k_<u></u>installation/nn nn.def failed<br>
atom Z RMT-max RMT<br>
Do you want to accept these radii; discard them; or rerun setRmt (a/d/r):<br>
a<br>
> nn (17:45:12) specify nn-bondlength factor: (usually=2) [and<br>
optionally dlimit, dstmax (about<br>
1.d-5, 20)]<br>
2<br>
DSTMAX: 20.0000000000000<br>
iix,iiy,iiz 5 5 5 39.6350000000000<br>
39.6350000000000 39.6350000000000<br>
NAMED ATOM: Ni1 Z changed to IATNR+999 to determine equivalency<br>
forrtl: severe (24): end-of-file during read, unit 20, file<br>
/home/raja/BaRuO3_wien2k/<u></u>constr_DFT/NiO/NiO.struct<br>
Image PC Routine Line Source<br>
nn 00000000004955EE Unknown Unknown Unknown<br>
nn 0000000000494086 Unknown Unknown Unknown<br>
nn 000000000044B462 Unknown Unknown Unknown<br>
nn 000000000040F1BC Unknown Unknown Unknown<br>
nn 000000000040E6DC Unknown Unknown Unknown<br>
nn 000000000042F7B5 Unknown Unknown Unknown<br>
nn 000000000042D648 Unknown Unknown Unknown<br>
nn 000000000040597A MAIN__ 218 nn.f<br>
nn 000000000040349C Unknown Unknown Unknown<br>
libc.so.6 00002B8A5151076D Unknown Unknown Unknown<br>
nn 0000000000403399 Unknown Unknown Unknown<br>
0.0u 0.0s 0:04.47 0.0% 0+0k 0+0io 0pf+0w<br>
error: command /home/raja/wien2k_<u></u>installation/nn nn.def failed<br>
<br>
stop error<br>
<br>
<br>
<br>
<br>
Regards<br>
wasim<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br></div></div><div class="im">
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
<br>
</div></blockquote><span class="HOEnZb"><font color="#888888">
<br>
-- <br>
Peter Blaha<br>
Inst.Materials Chemistry<br>
TU Vienna<br>
Getreidemarkt 9<br>
A-1060 Vienna<br>
Austria<br>
+43-1-5880115671</font></span><div class="HOEnZb"><div class="h5"><br>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
</div></div></blockquote></div><br></div>