<div>Dear All and Prof. Peter;</div> <div> Thank you very much for your kindly and rapidly help!</div> <div> I want to simulate a-Al2O3 Al-L2,3 edge, The supercell is Al11O18Al*. I use the first way you mentioned to solve this problem, I set Al Rmt=2.1, O Rmt=1.39 (R0=0.00005 for all of them). The case.outputst does not give any warnning. Is it ok?</div> <div> As to the second way, I don't understand how to operate it since I am still not very familar with this code. In the following I have attached the structer, in1c,in2c files. Could you show me a demo on the operation? And where you have reduced an electron in this second way? How to merge the 'two windown' calculation results into together at last?</div> <div> Another question: I test a-Al2O3 with kpoint=500 can convergence, how much k-point I should use in the supercell (Al11O18Al*) to keep convergence?</div> <div> Sorry for so many questions.</div>
<div>Yours;</div> <div> </div> <div> </div> <div><EM>2 ways to solve this:<BR><BR>Use as large as possible spheres for Al and choose a separation energy<BR>such that Al-2p is in the core. Check the core leakage, which should not<BR>be too large (but maybe a small miss is tolerable)<BR><BR>Perform a "2-window" calculation. This is a bit "tricky" but should be<BR>the optimal solution.<BR>After the "normal init_lapw" do:<BR>cp case.in1 case.in1s (or in1c -> in1cs)<BR>cp case.in2 case.in2s<BR><BR>Edit these files and put proper EMIN values in case.in1 (so that the Al<BR>2p states are NOT included; this should lead to a message like X<BR>EIGENVALUES below yy) , proper NE values in in2/in2s (The in2<BR>calculation should calculate all valence states except the Al 2p; the<BR>in2s should only calculate the Al 2p states).<BR><BR>I haven't done this for a long time, but it should still work.<BR><BR>For the XSPEC/TELNES calculations you would have to
put the 2p state<BR>into case.inc anyway; or you could also use "optic" to calculate this.<BR><BR>> Dear All;<BR>> I want to simulate Al atom L2,3 edge ELNES. Since L2,3 edge <BR>> corresponde to 2p->3d transition. To consider core hole effect, I have <BR>> to reduce one electron from Al 2p state (one Al in supercell). Am I right?<BR>> However, when we use the recommed energy -6 Ry to seperate core <BR>> state, 2p state does not include in the core state ( I did not find 2p <BR>> state in case.inc ). So how can I resolve this problem?<BR>> Thank you in advance !<BR><BR><BR>
P.Blaha<BR>-----------------------------------------------------------------------------------------------------</EM></div> <div><EM>case.struct</EM></div> <div>Al11O18Al <BR>H LATTICE,NONEQUIV.ATOMS: 18143_P3 <BR>MODE OF CALC=RELA
unit=bohr <BR> 8.950773 8.950773 24.398055 90.000000 90.000000120.000000 <BR>ATOM -1: X=0.66666667 Y=0.33333333 Z=0.48145870<BR> MULT= 1 ISPLIT= 4<BR>Al1 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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.33333333 Y=0.66666667 Z=0.01854130<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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.00000000 Y=0.00000000 Z=0.14812537<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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.66666667 Y=0.33333333 Z=0.18520797<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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.00000000 Y=0.00000000 Z=0.35187463<BR> MULT= 1 ISPLIT=
4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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.33333333 Y=0.66666667 Z=0.31479203<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z:
13.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.00000000 Z=0.64812537<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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.33333333 Y=0.66666667 Z=0.51854130<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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.66666667 Y=0.33333333 Z=0.68520797<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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 -10:
X=0.00000000 Y=0.00000000 Z=0.85187463<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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 -11: X=0.33333333 Y=0.66666667 Z=0.81479203<BR> MULT= 1 ISPLIT=
4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z: 13.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 -12: X=0.66666667 Y=0.33333333 Z=0.98145870<BR> MULT= 1 ISPLIT= 4<BR>Al2 NPT= 781 R0=0.00005000 RMT= 2.1000 Z:
13.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 -13: X=0.36108030 Y=0.33333333 Z=0.08333333<BR> MULT= 3 ISPLIT= 8<BR> -13: X=0.66666667 Y=0.02774697 Z=0.08333333<BR> -13: X=0.97225303 Y=0.63891970 Z=0.08333333<BR>O NPT= 781 R0=0.00005000 RMT= 1.3900 Z:
8.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 -14: X=0.00000000 Y=0.30558637 Z=0.25000000<BR> MULT= 3 ISPLIT= 8<BR> -14: X=0.69441363 Y=0.69441363 Z=0.25000000<BR> -14: X=0.30558637 Y=0.00000000 Z=0.25000000<BR>O NPT= 781 R0=0.00005000 RMT= 1.3900 Z:
8.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 -15: X=0.02774697 Y=0.66666667 Z=0.41666667<BR> MULT= 3 ISPLIT= 8<BR> -15: X=0.33333333 Y=0.36108030 Z=0.41666667<BR> -15: X=0.63891970 Y=0.97225303 Z=0.41666667<BR>O NPT= 781 R0=0.00005000 RMT= 1.3900 Z:
8.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 -16: X=0.36108030 Y=0.02774697 Z=0.58333333<BR> MULT= 3 ISPLIT= 8<BR> -16: X=0.97225303 Y=0.33333333 Z=0.58333333<BR> -16: X=0.66666667 Y=0.63891970 Z=0.58333333<BR>O NPT= 781 R0=0.00005000 RMT= 1.3900 Z:
8.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 -17: X=0.00000000 Y=0.69441363 Z=0.75000000<BR> MULT= 3 ISPLIT= 8<BR> -17: X=0.30558637 Y=0.30558637 Z=0.75000000<BR> -17: X=0.69441363 Y=0.00000000 Z=0.75000000<BR>O NPT= 781 R0=0.00005000 RMT= 1.3900 Z:
8.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 -18: X=0.02774697 Y=0.36108030 Z=0.91666667<BR> MULT= 3 ISPLIT= 8<BR> -18: X=0.63891970 Y=0.66666667 Z=0.91666667<BR> -18: X=0.33333333 Y=0.97225303 Z=0.91666667<BR>O NPT= 781 R0=0.00005000 RMT= 1.3900 Z:
8.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> 3 NUMBER OF SYMMETRY OPERATIONS<BR>-1 1 0 0.0000000<BR>-1 0 0 0.0000000<BR> 0 0 1 0.0000000<BR> 1<BR> 0-1 0 0.0000000<BR> 1-1 0 0.0000000<BR> 0 0 1 0.0000000<BR> 2<BR> 1 0 0 0.0000000<BR> 0 1 0 0.0000000<BR> 0 0 1 0.0000000<BR>
3<BR>--------------------------------------------------------------------------------------------------</div> <div>case.in1c</div> <div> </div> <div>WFFIL (WFPRI, SUPWF)<BR> 8.00 10 4 (R-MT*K-MAX; MAX L IN WF, V-NMT<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL
E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT
1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2
0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 2 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0
-1.53 0.010 CONT 1<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 -1.53 0.010 CONT 1<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 -1.53 0.010 CONT 1<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 3
0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 -1.53 0.010 CONT 1<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 -1.53 0.010 CONT 1<BR> 0 0.30 0.000 CONT 1<BR> 1 0.30 0.000 CONT 1<BR> 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR> 0 -1.53 0.010 CONT 1<BR> 0 0.30 0.000 CONT 1<BR> 1
0.30 0.000 CONT 1<BR>K-VECTORS FROM UNIT:4 -6.0 1.5 emin/emax window<BR>---------------------------------------------------------------------------------------------------------</div> <div>case.in2c</div> <div>TOT (TOT,FOR,QTL,EFG,FERMI)<BR> -9.0 144.0 0.50 0.05 EMIN, NE, ESEPERMIN, ESEPER0<BR>TETRA 0.000 (GAUSS,ROOT,TEMP,TETRA,ALL eval)<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4
0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0
0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 2 0 3 0 3 3 -3 3 4 0 4 3 -4 3 5 0 5 3 -5 3 6 0 6 3 -6 3 6 6 -6 6<BR> 0 0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5 1 -5 1 5 2 -5 2 5 3 -5 3 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6 -6 6<BR> 0
0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5 1 -5 1 5 2 -5 2 5 3 -5 3 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6 -6 6<BR> 0 0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5 1 -5 1 5 2 -5 2 5 3 -5 3 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6 -6 6<BR> 0 0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5
1 -5 1 5 2 -5 2 5 3 -5 3 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6 -6 6<BR> 0 0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5 1 -5 1 5 2 -5 2 5 3 -5 3 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6 -6 6<BR> 0 0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5 1 -5 1 5 2 -5 2 5 3 -5 3 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6 -6
6<BR> 14. GMAX<BR>NOFILE FILE/NOFILE write recprlist<BR></div><p> 
                <hr size=1><a href="http://cn.mail.yahoo.com" target=blank>ÇÀ×¢ÑÅ»¢Ãâ·ÑÓÊÏä-3.5GÈÝÁ¿£¬20M¸½¼þ£¡</a>