Re: [Wien] HDLO and LVNS

2017-11-13 Thread pieper
High Pablo, I think there is a reasonable amount of information in the UG, including a citation for further reading (which I never followed): UG Version 17.1, ch. 7.5, first par.: When adding such LOs, make sure their E-parameters are far away from each other (typically at least 1 Ry, often

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread Peter Blaha
You must modify your O2 struct file. Make the 2 oxygens "equivalent", i.e. you should have only ONE non-equivalent atomic site (NAT=1), but MULT=2 for the 2 oxygens. This will produce inversion symmetry and all calculations will take less time and memory !!! And yes:

[Wien] qtlb erro of Ba

2017-11-13 Thread chin Sabsu
Dear Sir,I faced QTL-B error and I changed Ef in case.in1 to -0.5 but it automatically restored to the original value after a couple of cycles. Please suggest me what should I change so that the qtlb error do not occur. I also used in1new 2 and the scf went more than 70 cycles and  I again saw

[Wien] XCPOT3 - Error

2017-11-13 Thread chin Sabsu
Hello, What I asked about IFFT parameters in that if they should be same in all cases then I have a problem now: I run ZrO2 (monoclinic) for optimization and getting "XCPOT3 - Error" at the beginning of the one middle scf step while I have used Gmax=24. I am sure if I use IFFt as  0 0 0 2 then

[Wien] QTL WARN differ between lapw2 and lapw2 -qtl

2017-11-13 Thread Jichao Chang
Dear WIEN2k users, My Question: 1) Why "x lapw2" and "x lapw2 -qtl" show a different warning? Description: I caliculate fcc Ni( struct file is attached below) by "run_lapw -p -ec 0.0001 -cc 0.0001"and I got QTL-B warn from scf file: ``` :WARN : QTL-B value eq. 2.82 in Band of energy 0.64306

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread chin Sabsu
Thanks Sir for the nice explanation and introducing one more factor  "IFFT". Sometimes we get XCOPT error and we set it 0 0  0 2. So I was doing wrong. If I was doing three compounds A, B, and C then let for A, and B everything was okay but for C the linear dependency failed then I set IFFT to

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread Peter Blaha
Thanks for the correction. With respect to your other questions: In order to save cpu time, I would do (at least for the more complicated compounds, not necessarily for BaO, which is fast anyway) all calculations with "standard RMTs" (volume optimization, optimization of internal parameters,

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread Laurence Marks
Let me add one important thing to what Peter sent. My experience is that the atomic positions are only weakly sensitive to the exact parameters used, although this may not be true if you are doing phonons. Hence it has always worked well for me to calculate (including MSR1a relaxations) using

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread Laurence Marks
N.N.B., instead of a smaller ecut one can use .lcore. Sometimes this also improves convergence. _ Professor Laurence Marks "Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi www.numis.northwestern.edu On Nov 13, 2017 2:12 PM,

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread Laurence Marks
"besides the statement: be careful with the energy parameter (very fine E-search in case.in1)" N.B., one does not need a fine search in the more recent Wien2k, only in the older versions. I seem to remember that GMax 20 mainly cancels out for LDA/GGA -- I would always use it (or more) for a meta

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread SM Alay-e-Abbas
Good to see more details added to the FAQ on cohesive/formation energies. I have two suggestions for the new FAQ on cohesive energies 1- Computing the atomic energies of atoms like N and P in an FCC cell is ok, however for O atom the high symmetry of the FCC cell results in 1/3 occupancies (for

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread Peter Blaha
1- Computing the atomic energies of atoms like N and P in an FCC cell is ok, however for O atom the high symmetry of the FCC cell results in 1/3 occupancies (for the 4th p electron of O) in the spin down case. Only using a lower symmetry cell (orthorhombic) for O atom eliminates this issue.

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread Peter Blaha
Please check the faq pages at www.wien2k.at/reg_user/faq. I've made a couple of modifications and there is a page for cohesive and formation energies. Any suggestions are welcome. Peter Blaha On 11/13/2017 11:50 AM, Peter Blaha wrote: You must modify your O2 struct file. Make the 2

Re: [Wien] O2 molecule continue (rmt issue)

2017-11-13 Thread chin Sabsu
Thank you very much Sir for a very nice NOTE. I suggest one correction and two more queries: Correction: The RKmax 5.6 in the sentence "Then you should do atom A with RKmax=7/2.5*2.0=5.6" should be replaced by  7/2*2.5=8.75!! Query: 1. I could not understand what is the purpose of another