[Wien] Initialization with spin-orbit

2016-01-26 Thread Islam, Md F
Dear Wien2k users, I am having trouble with initializing a job with initso. I am doing a substitutional magnetic impurity calculation in a surface supercell (only 2 symmetry operations). After I relaxed the surface, I did an SCF calculation (more accurate) without spin-orbit which is conve

Re: [Wien] (no subject)

2016-01-26 Thread said chibani
I did "nt find "" "write(22" statements by "write(21"" in mixer.F please give a solution 2016-01-21 18:26 GMT+01:00 Peter Blaha : > If the error is really caused by gfortran and not by an "user-error", > you can fix it by: > > cd $WIENROOT/SRC_mixer/ > edit mixer.F and replace two "writ

[Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Hu, Wenhao
Hi, all: I have a basic question on the DOS calculation of onsite-hybrid functional (ONF) calculation. When I check the UG, I can’t find any notes about the DOS plotting of ONF, so I assume the program flow is the same in ONF as normal (x lapw2 -p -c -up/dn -qtl, then x tetra -p -up/dn). But w

Re: [Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Laurence Marks
Did you fully converge the calculation or just do one iteration? You have to fully converge. Also, what does "grep :MV case.scf | tail" give? On Tue, Jan 26, 2016 at 11:08 AM, Hu, Wenhao wrote: > Hi, all: > > I have a basic question on the DOS calculation of onsite-hybrid functional > (ONF) calc

Re: [Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Hu, Wenhao
Hi, Prof. Laurence: This result is converged to large extent with respect to force, charge and energy. The forces on individual atoms are all below 0.5 and the total force in z is 3.5, which is the direction without symmetry. The :MV flag gives the following information: :MVORD NDM 100 L1

Re: [Wien] wien2k workshop 2016

2016-01-26 Thread Gunakar Das
Respected Blaha, I am staying in the extreme north-part of India keenly interested to join your progrmme. Since dollar (for registration etc,) takes a little bit time for conversion from Rupees ( Indian currency) of some days due to some official procedure, I am just requesting you to give me som

Re: [Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Laurence Marks
That looks converged. Are you using -so, or is the orbital potential being applied in lapw1? (Look for lapw1 -up -orb lines in :log). On Tue, Jan 26, 2016 at 11:43 AM, Hu, Wenhao wrote: > Hi, Prof. Laurence: > > This result is converged to large extent with respect to force, charge and > energy.

Re: [Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Hu, Wenhao
Hi, Prof. Laurence: No, I don’t include the SO effect. The orbital potential is applied in lapw1. Basically, I’m using the -eece flag of runs_lapw. The program flow is as follows: > (runeece_lapw) options: -p Sun Jan 24 14:24:18 CST 2016> (x) lapwdm -up -p -c Sun Jan 24 14:24:40 CST 2016> (x

Re: [Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Laurence Marks
Hard to say. Thoughts: Check case.scf2up/dn and look at the Cr d-occupancy grep :EORB *scforb* If everything is zero then no orbital terms were being applied. Check case.vorbup/dn -- non zero terms? Check case.scforbup/dn Last line should be PRATT mixing scheme with 1.000 head case.in2

Re: [Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Hu, Wenhao
Hi, Prof. Laurence: By following your suggestion, It seems that everything is normal except for the case.in2eece file. My case.in2eece is as follows: TOT EECE (TOT,FOR,QTL,EFG,FERMI) -12.0 52.00 0.50 0.05 1 EMIN, NE, ESEPERMIN, ESEPER0, iqtlsave TETRA0.000 (GAUSS,ROOT,TEMP,TETRA,ALL

Re: [Wien] Basic question about DOS of onsite-hybrid functional calculation

2016-01-26 Thread Laurence Marks
Caveat: you have to modify case.in2, not case.in2eece. Add a few spaces (2-3) before the 452, and then do runeece -p as this just does the part where case.in2eece is created. On Tue, Jan 26, 2016 at 2:51 PM, Hu, Wenhao wrote: > Hi, Prof. Laurence: > > By following your suggestion, It seems that

Re: [Wien] Initialization with spin-orbit

2016-01-26 Thread Peter Blaha
Not exactly sure what can cause the problem. Please search the mailing list, there was some time ago a patch for symmetso. From your description of the problem (optimized/non-optimized structure) it could be a symmetry problem, i.e. during optimization (hexagonal case ???) the symmetry of th