Re: [Wien] Changing magnetization direction

2018-06-22 Thread Fecher, Gerhard
Have a look into the optimization script, you just change the case.inso files instead of the case.struct files indeed, you have to create them by hand and you have to include some so specific things that may look like x symmetso # -c cp case.struct_socase.struct cp

Re: [Wien] ELF calculation

2018-06-06 Thread Fecher, Gerhard
d H. Fecher Institut of Inorganic and Analytical Chemistry Johannes Gutenberg - University 55099 Mainz and Max Planck Institute for Chemical Physics of Solids 01187 Dresden Von: Wien [wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Fecher, Gerhard [fec...@un

Re: [Wien] ELF calculation

2018-06-06 Thread Fecher, Gerhard
Dear Victor, just for curiosity, is there any example how to calculate the ELF with Critic2 from Wien2k data (which ? and how to produce them ?) (I never tried and thus never searched for it up to todays discussion) Seems this was a part of the question. BTW.: I have to reduce very often the

Re: [Wien] spagh2rho Output Files

2018-05-23 Thread Fecher, Gerhard
Maybe a look into the code helps if(invert.eq.0) then write(iunit,*) nx,ny,x0,y0 write(iunit,11) ((a2(i,j),j=1,ny),i=1,nx) else write(iunit,*) 2*nx-1,2*ny-1,x0,y0 write(iunit,11)((a2(abs(i)+1,abs(j)+1),j=-ny+1,ny-1),i=-nx+1,nx-1) endif 11

Re: [Wien] (no subject)

2018-05-23 Thread Fecher, Gerhard
Dear Sudipta, then you will need to do calculations with spin orbit interaction take care that the quantiisation axis in case.inso is the same as the magnetisation direction in case.inorb if you like to have different angles then you may start to rotate the direction of the quantisation axis

Re: [Wien] ferrimagnetic

2018-05-21 Thread Fecher, Gerhard
Increase the number of k points and see what happens Btw, an ferrimagnet may also exist if the atoms are of the same kind, it is a matter of symmetry, In antiferromagnets the sites are equivalent, in (compensated) ferrimagnets they are not. Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers

Re: [Wien] Problem in volume optimization

2018-05-19 Thread Fecher, Gerhard
NN error ==> RMT to big See recent post few days ago Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question is." Dr. Gerhard H. Fecher

Re: [Wien] Problems when trying to plot E vs c/a

2018-05-19 Thread Fecher, Gerhard
] Gesendet: Freitag, 18. Mai 2018 22:57 An: A Mailing list for WIEN2k users Betreff: Re: [Wien] Problems when trying to plot E vs c/a I can plot E vs volume but not E vs c/a with the previous optimize.pl<http://optimize.pl>. On Sat, 19 May 2018, 2:09 am Fecher, Gerhard, <fec...@uni-

Re: [Wien] Problems when trying to plot E vs c/a

2018-05-18 Thread Fecher, Gerhard
; your attached one the I get an error Failed to exec /home/dipraj/wien2k/SRC_w2web/htdocs/exec/optimize.pl<http://optimize.pl> : Permission denied On 17 May 2018 at 17:00, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: Hallo Peter, thanks for the fil

Re: [Wien] Problems when trying to plot E vs c/a

2018-05-17 Thread Fecher, Gerhard
Re: [Wien] Problems when trying to plot E vs c/a Thanks for the report. Modified eplot_lapw and SRC_w2web/htdocs/exec/optimize.pl attached. On 05/16/2018 04:20 PM, Fecher, Gerhard wrote: > Dear c/a fitters, > This concerns the latest Wien2k version > I recei

[Wien] Problems when trying to plot E vs c/a

2018-05-16 Thread Fecher, Gerhard
Dear c/a fitters, This concerns the latest Wien2k version I receive only the content of test_opt.analysis when I try with w2web to plot E vs c/a but neither the result of the fit nor the plot are shown, this seems to be a problem with the present version of the eplot script when I

Re: [Wien] error in symmetry step (not seen by sgroup)

2018-05-15 Thread Fecher, Gerhard
Strange, just as a diagnostic it seems the error appears only when running batch initialisation I do not see it when using step by step initialisation in w2web or running x symmetry directly from the command line also strange, when converting the structure to P1 then sgroup finds space group 146

Re: [Wien] allowing arbitrary name of the directory containing case files

2018-05-14 Thread Fecher, Gerhard
use the save to a folder feature that is, save each of your calculation to an appropriate folder inside Fe, say Fe_5x5x5 or rkm7_5x5x5 or whatever you like to have Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that

Re: [Wien] Fe3O4 run scf ROTDEF error

2018-05-10 Thread Fecher, Gerhard
of index 0.7448000 0.5052000 0.5052000 I am posting the cif file I have used. On 10 May 2018 at 18:13, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: Then you are doing something wrong, that I do not know (It works on my computer without problems) Ciao G

Re: [Wien] Fe3O4 run scf ROTDEF error

2018-05-10 Thread Fecher, Gerhard
same error. There was no warning regarding space group also. On 9 May 2018 at 20:30, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: maybe it is the setup (choice of origin) of space group 227 use sgroup and accept (!) the new structure, then initialize SO as usual be s

Re: [Wien] Fe3O4 run scf ROTDEF error

2018-05-09 Thread Fecher, Gerhard
maybe it is the setup (choice of origin) of space group 227 use sgroup and accept (!) the new structure, then initialize SO as usual be sure that you have only a case.struct and not any case.struct_st files in your directory (principally you can start from the struct you did send, maybe it

Re: [Wien] Structural optimization, volume optimization

2018-05-03 Thread Fecher, Gerhard
87 Dresden Von: Wien [wien-boun...@zeus.theochem.tuwien.ac.at] im Auftrag von Lyudmila Dobysheva [lyuk...@mail.ru] Gesendet: Donnerstag, 3. Mai 2018 11:59 An: wien@zeus.theochem.tuwien.ac.at Betreff: Re: [Wien] Structural optimization, volume optimization 02.05.2018 11:20, Fecher, Ger

Re: [Wien] Installation with MPI and GNU compilers

2018-05-02 Thread Fecher, Gerhard
y Johannes Gutenberg - University 55099 Mainz and Max Planck Institute for Chemical Physics of Solids 01187 Dresden Von: Pavel Ondračka [pavel.ondra...@email.cz] Gesendet: Mittwoch, 2. Mai 2018 12:05 An: Fecher, Gerhard Betreff: Re: [Wien] Installation with M

Re: [Wien] Installation with MPI and GNU compilers

2018-05-02 Thread Fecher, Gerhard
I never checked that: does the -DINTEL_VML switch correspond to the VML library routines of MKL or to the SVML library routines of the compiler this makes a difference, the svml routines are automatically invoked by the INTEL compiler if one uses -O2 optimization or higher. (check also the usage

Re: [Wien] Structural optimization, volume optimization

2018-05-02 Thread Fecher, Gerhard
The volume optimization is part of an optimization of the structure, but only one part of it. To optimize the structure you need - optimization of the lattice parameters (a,b,c, alpha, beta, gamma, or a set out of them depending on the space group) - optimization of the atomic positions (if the

Re: [Wien] optimization of pseudocubic structure

2018-03-29 Thread Fecher, Gerhard
If it is not cubic but orthorhombic (or hexagonal, or tetragonal) then you have to optimize as orthorhombic (etc.) or for whatever structure you mean with pseudocubic Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with

Re: [Wien] Band structure and Dos contradiction of an alloy

2018-03-19 Thread Fecher, Gerhard
his should not be the problem since spaghetti gives no gap. On Monday 2018-03-19 10:38, Fecher, Gerhard wrote: >Date: Mon, 19 Mar 2018 10:38:31 >From: "Fecher, Gerhard" <fec...@uni-mainz.de> >Reply-To: A Mailing list for WIEN2k users <wien@zeus.theochem.tuwien.ac.at> &

Re: [Wien] Band structure and Dos contradiction of an alloy

2018-03-19 Thread Fecher, Gerhard
Did you check the spaghettis for all relevant directions ? If not, plot the band structure for other directions. Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question

Re: [Wien] Reg: File to get the number of electrons in majority and minority spin cases for d-orbital

2018-03-04 Thread Fecher, Gerhard
you find the charges and magnetic moments inside the spheres as well as the interstitial in the file case.scfm (indeed also in case.scf) Your remark " As these compounds are metallic in nature the arrived magnetism is due to the itinerant electrons." might not be the truth, in case of Mn you

Re: [Wien] query on optimization

2018-02-08 Thread Fecher, Gerhard
what does it mean " the program do not produce the exact relaxed structure on the optimized lattice parameter" ? what is wrong ? how did you decide that it is wrong ? where do you have your initial parameters from ? what material do you calculate ? read the last question first ! Ciao

Re: [Wien] a note on lattice parameter optimization

2018-02-08 Thread Fecher, Gerhard
As someone who is doing both, experiment and calculations, I do not agree for the following reasons a) For example, in the Pearson database you will find many compounds where (in case of "free" parameters) the atom positions of the prototype are given but not the "real" ones of the actual

Re: [Wien] zigzag potential interpretation

2018-01-03 Thread Fecher, Gerhard
and we are testing this new version with Peter and Pavel at this moment. Happy new year to you Gerhard Xavier Le 03/01/2018 à 14:24, Fecher, Gerhard a écrit : > Dear Stefaan, > I am not realy sure what difference you expect, > I do not see why at two seemingly same surfaces the size of the

Re: [Wien] zigzag potential interpretation

2018-01-03 Thread Fecher, Gerhard
Dear Stefaan, I am not realy sure what difference you expect, I do not see why at two seemingly same surfaces the size of the magnetic moment (orbital or spin) should depend on their orientation in the sense that it is parallel or antiparallel to the surface normal. I wonder about the

Re: [Wien] stability issue

2017-12-22 Thread Fecher, Gerhard
Isn't it simply possible that the binaries AO + BO2 are more stable than the ternary ABO3 ? BTW.: you find tons of papers where quaternaries or ternaries are predicted to exist even so the binaries (possibly in combination with elements) are more stable and a phase separation will occur. Ciao

Re: [Wien] different MLD for bcc structure for magnetic equivalent directions M001, M010 and M100

2017-11-26 Thread Fecher, Gerhard
this would be very helpful Thank you and with my best regards Jaroslav On 25/11/17 14:13, Fecher, Gerhard wrote: > Hi Jaroslav, > > with SO, 001 is not equivalent to 001 or 010, if the magnetisation is along > 001 > this you see easily from the changed symmetry after

Re: [Wien] different MLD for bcc structure for magnetic equivalent directions M001, M010 and M100

2017-11-25 Thread Fecher, Gerhard
Hi Jaroslav, with SO, 001 is not equivalent to 001 or 010, if the magnetisation is along 001 this you see easily from the changed symmetry after initializing SO (symmetso) regards from Dresden Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be

Re: [Wien] Calculations of Slater parameters (or Slater integrals) for open-shell (3d or 4f) ions

2017-11-05 Thread Fecher, Gerhard
You may like to read http://susi.theochem.tuwien.ac.at/reg_user/textbooks/Constraint_U.pdf that contains some notes by Georg K. H. Madsen and Pavel Novac on constrained LDA for those parameters Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the

Re: [Wien] negative eps-1 and very large value of eps-2

2017-10-28 Thread Fecher, Gerhard
But you did not read it completely ==> runsp -so # spin-polarized + SO includes lapw1 -dn, isn't it ? what are youre values for n and k (see case.refraction) what happens in the case k > n and what is its meaning ? Why and when does it appear in metals ? how are the optical

Re: [Wien] why we need to divide plasma frequency by root(2)?

2017-10-27 Thread Fecher, Gerhard
The optics programm tells about the SP Plasma frequencies: w_pl = sqrt( w_pl^2(up-spin) + w_pl^2(dn-spin) ) where does this come from remember the classical approach where the plasma frequency is proportional to root(n) with n being the density of free electrons (number of electrons per

Re: [Wien] negative eps-1 and very large value of eps-2

2017-10-27 Thread Fecher, Gerhard
If your case is spinpolarised why don't you run lapw1 -dn before running lapwso ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question is."

Re: [Wien] optimization of an orthorhombic cell

2017-10-12 Thread Fecher, Gerhard
you have three quantities V=a*b*c, b=a*(b/a) and c=a*(c/a) so what is the problem ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question is."

Re: [Wien] Optical properties with SO coupling

2017-10-06 Thread Fecher, Gerhard
nt program. Then, x kram -up is just fine. With my regards Jaroslav On 05/10/17 15:14, Fecher, Gerhard wrote: > Hi Jaroslav, > if you check only case.jointup it has possibly only half the value because > the other half is supposed to be in case.jointdn > (with SO they should be the

Re: [Wien] Optical properties with SO coupling

2017-10-05 Thread Fecher, Gerhard
g16054.html On 10/5/2017 7:14 AM, Fecher, Gerhard wrote: Hi Jaroslav, if you check only case.jointup it has possibly only half the value because the other half is supposed to be in case.jointdn (with SO they should be the same) Did you try to copy case.jointup to case.jointdn (or run in additio

Re: [Wien] Optical properties with SO coupling

2017-10-05 Thread Fecher, Gerhard
Hi Jaroslav, if you check only case.jointup it has possibly only half the value because the other half is supposed to be in case.jointdn (with SO they should be the same) Did you try to copy case.jointup to case.jointdn (or run in addition everything for dn) and then addjoint then the factor 2

Re: [Wien] Optical properties with SO coupling

2017-10-04 Thread Fecher, Gerhard
I can not reproduce your problem, I always receive ALL columns that I ask for in case.inop (indeed the number of columns in case.injoint should match) did you change your case.inop between the steps you performed and forgot to change the number of choices ? Ciao Gerhard DEEP THOUGHT in D.

Re: [Wien] Fwd: how to simulate the energy of Oxygen molecule

2017-09-22 Thread Fecher, Gerhard
I think 4.4 Angström between two oxygen atoms would be a rather large bond length for O2 shouldn't you give the z parameter in multiples (fractions) of the lattice parameters ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with

[Wien] Problem with spaghetti (bz_lin.f)

2017-09-18 Thread Fecher, Gerhard
I experienced a problem with the plot of the band structure in version 17: Additional (needless) vertical lines appear in the plots for hexagonal (and maybe similar in other) band structures between the M and K points A comparison to the older version of bz_lin.f reveals that the tolerance was

Re: [Wien] About the magnetic moment of vanadium in vanadium sulphide

2017-09-12 Thread Fecher, Gerhard
Did you read and understand Rationale for mixing exact exchange with density functional approximations John P. Perdew, Matthias Ernzerhof, and Kieron Burke The Journal of Chemical Physics 105 , 9982 (1996) or similar articles and references there ? Usually, the parameters of the

Re: [Wien] Total magnetic moment non-integer ??!

2017-09-12 Thread Fecher, Gerhard
In addition, what size of the gap do you have, see e.g. in case.scf2 Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question is." Dr.

Re: [Wien] About the magnetic moment of vanadium in vanadium sulphide

2017-09-08 Thread Fecher, Gerhard
First, as Gavin told, your convergence criteria are too bad to decide about your magnetic moments just try to improve them in steps say up to ec 10^-6 cc 10^-3 and see what happens with the moments (that was why I asked fror it) however be sure you have enough k-points I assume you used the

Re: [Wien] About the magnetic moment of vanadium in vanadium sulphide

2017-09-08 Thread Fecher, Gerhard
Now I am a little confused do you think that it was the convergence criteria that I found a magnetic moment of 2.5 mu_B instead of 0.05 mu_B at the V ? or did I make another mistake ;-) Indeed, you are right, the 0.05 mu_B might be a result of to bad convergence however, to check if such a

Re: [Wien] About the magnetic moment of vanadium in vanadium sulphide

2017-09-07 Thread Fecher, Gerhard
On Behalf Of Fecher, Gerhard Sent: 07 September 2017 08:12 To: A Mailing list for WIEN2k users Subject: Re: [Wien] About the magnetic moment of vanadium in vanadium sulphide 0.05 muB does not mean that it is antiferromagnetic ! what was your charge convergence criterion ? You did never answer m

Re: [Wien] About the magnetic moment of vanadium in vanadium sulphide

2017-09-07 Thread Fecher, Gerhard
0.05 muB does not mean that it is antiferromagnetic ! what was your charge convergence criterion ? You did never answer my question whether you started the EECE calculation from a converged GGA calculation. Why do you like to have an afm state when the experiment tells it is not ? Ciao

Re: [Wien] About the magnetic moment of vanadium in vanadium sulphide

2017-09-04 Thread Fecher, Gerhard
About what moment are you talking, the total magnetic moment or the magnetic moment of the Cr atoms ? Did you start your EECE calculation from a regular GGA calculation that had no magnetic moments at the Cr ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think

Re: [Wien] How to know if the M+2 has a high spin from calculation?

2017-08-28 Thread Fecher, Gerhard
Hallo Abderrahmane, you find the occupation of the states (INSIDE THE SPERES !) in the case.scf2up/dn files. See the lines with e.g.: :CHA002: TOTAL VALENCE CHARGE INSIDE SPHERE :PCS002: PARTIAL CHARGES SPHERE = 2 . :QTL002: the content should be selfexplaining. Take care, you

Re: [Wien] How to replace the availiable session list

2017-08-24 Thread Fecher, Gerhard
You may continue your work in an existing directory if there were no troubles before If you create a session for an existing case (in an old folder), then you need to check that the session info is ok (see session mangement in w2web) in particular complex or not, spin polarized or not, etc.

Re: [Wien] How to simulate the ionic state of a compound?

2017-08-23 Thread Fecher, Gerhard
Funny discussion, I wonder wether iron (Z=26) knows that it is Nickel just if one gives it a new name ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question is."

Re: [Wien] Problem with dstart in new version 17.1

2017-08-14 Thread Fecher, Gerhard
aN's being output. I blame this on Intel, although I may be wrong. On Mon, Aug 14, 2017 at 7:49 AM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: I experienced a problem with dstart of version 17.1, an error appears as given below when running a calcul

[Wien] Problem with dstart in new version 17.1

2017-08-14 Thread Fecher, Gerhard
I experienced a problem with dstart of version 17.1, an error appears as given below when running a calculation for the attached structure (initialisation with defaults). A closer look, however, shows that this error does not result from lapw2 but from dstart. (I found NANs in some files used

Re: [Wien] True half métal and false half métal

2017-08-14 Thread Fecher, Gerhard
Please read and understand C. Felser and G. H. Fecher (Eds.) Spintronics : From Materials to Devices in particular the chapters on Heusler compounds ISBN 978-90-481-3831-9 ISBN 978-90-481-3832-6 (eBook) DOI 10.1007/978-90-481-3832-6 Springer Dordrecht Heidelberg New York London Ciao Gerhard

Re: [Wien] experimental vs relaxed parameters using wien2k

2017-07-31 Thread Fecher, Gerhard
If you do a calculation then the lattice parameter(s) for the optimized structure is (are) correct for the used functional as it results in the lowest energy, indeed only for this functional (supposed you realy searched and found the absolute minimum for the correct magnetic order). That's it.

Re: [Wien] PBE vs mBJ

2017-07-31 Thread Fecher, Gerhard
This question is just not scientific sound. You choose the functional and this functional gives you some answer on the electronic structure, this answer is not the universal answer, it is just that what you asked the computer to do. Remember that the A in LDA, etc. means approximation ! It is

Re: [Wien] energy limit problems

2017-07-28 Thread Fecher, Gerhard
I wonder whether such an optimization makes sense at all, this would optimize the volume of the vacuum, is it that what you want ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known

Re: [Wien] experimental vs relaxed parameters using wien2k

2017-07-27 Thread Fecher, Gerhard
Dear Tristan, don't you think that this is a rather strange advice, it might suggest that the optimized lattice parameters are bad, whereas any arbitrary choosen lattice parameters are good. I do not think it is a good idea to vary the lattice parameters until an arbitrarily choosen property

Re: [Wien] Questions about imposing external magnetic field on no-magnetic system

2017-07-27 Thread Fecher, Gerhard
do you search for something like SKEAF by P.M.C. Rourke found in the software goodies ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question is."

Re: [Wien] experimental vs relaxed parameters using wien2k

2017-07-25 Thread Fecher, Gerhard
On Mon, Jul 24, 2017 at 2:09 PM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: I would first check if the structure was realy completely solved in previous work. I found that e.g. in the Pearson's crystal data base the structural parameters (atomic positions

Re: [Wien] experimental vs relaxed parameters using wien2k

2017-07-24 Thread Fecher, Gerhard
I would first check if the structure was realy completely solved in previous work. I found that e.g. in the Pearson's crystal data base the structural parameters (atomic positions) for the prototype are reported for many compounds. If one uses these parameters (and they do not accidentially fit

[Wien] IRREP and SYMMETSO with ferromagnetic materials

2017-07-20 Thread Fecher, Gerhard
Hi Peter, Maybe this is a little longer but I hope it helps to correct or improve SYMMETSO and IRREP. It concerns for example https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg14673.html I did not realise that I was sticking with that problem not that long ago

Re: [Wien] Questions about imposing external magnetic field on no-magnetic system

2017-07-19 Thread Fecher, Gerhard
An additional thought about the B field effect. The answer to what a B field is doing to the electronic structure might have also some more subtle aspects. One aspect is -- as usual -- the role of symmetry. Similar to the case of spin-orbit interaction in a ferromagnet, the application of a

Re: [Wien] magnetic susceptibility for a ferromagnetic metal or for a ferromagnetic semiconductor

2017-07-18 Thread Fecher, Gerhard
How do you define the magnetic susceptibility ? Did you hear about the field dependence of the magnetisation and hysteresis in ferromagnets, isn't the permeability of an ferromagnet depending on the applied magnetic field ? I suggest to read some tectbook on magnetism, e.g. by M. Coeye or by J.

Re: [Wien] symmetry problem for spontaneous polarization

2017-07-17 Thread Fecher, Gerhard
you need first to initialize the "low symmetry" structure with the displaced atom(s), and afterwards run the "high" symmetry structure with the atoms placed back (do not accept symmetry changes during initialisation). Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I

[Wien] irrep has a problem with SO of ferromagnetic materials: X not equal for all elements in the class

2017-07-11 Thread Fecher, Gerhard
Dear all, irrep cannot find the irrep of a state in a ferromagnet with SO but produces (here for fcc Ni): bnd ndg eigval E 2C4 C2 2C2`2C2" I2IC4 IC22IC2` 2IC2" 1 1 -4.267251 1.00+0.00i 0.71-0.71i STOP: X=(

Re: [Wien] Questions in "initso_lapw" when using "runsp_c_lapw" to run my task

2017-07-11 Thread Fecher, Gerhard
Another "misuse" of U might be to shift the binding energies of the Lu 4f electrons in "nonmagnetic" compounds without U the Lu 4f doublet appears, roughly speeking, centered at about -5eV, with a moderate U of 0.5 Ry one can shift it to the value found in photoelectron spectroscopy at about

Re: [Wien] error in vorb continues...

2017-07-06 Thread Fecher, Gerhard
at 7:30 PM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: I guess youre rmt's are too large, did you reduce them during initialisation (when saving the structure) to account for the 10% volume reduction ? Did you receive other error files that were not empty

[Wien] Problem with clmextrapol_lapw / clmaddsub in spin polarized SO cases

2017-07-06 Thread Fecher, Gerhard
This issue probably appeared already in older versions, here I found it with the latest 17.1 When trying to use clmextrapol_lapw in a optimize job for a spin polarized case with SO, it fails. (note using dstart up/dn instead of clmextrapol works without problems) I used Fe as simple test case

Re: [Wien] error in vorb continues...

2017-07-06 Thread Fecher, Gerhard
alculation?? Looking forward to your suggestions. with regards, On Thu, Jul 6, 2017 at 2:18 PM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: whoever you are, I think you are answering the wrong topic Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide t

Re: [Wien] 答复: Error information during installation of WIEN2k_16.1 on CentOs7

2017-07-06 Thread Fecher, Gerhard
tp://aka.ms/weboutlook> 发件人: Wien <wien-boun...@zeus.theochem.tuwien.ac.at> 代表 Fecher, Gerhard <fec...@uni-mainz.de> 发送时间: 2017年7月6日 8:56 收件人: A Mailing list for WIEN2k users 主题: Re: [Wien] Error information during installation of WIEN2k_16.1 on Cen

Re: [Wien] Error information during installation of WIEN2k_16.1 on CentOs7

2017-07-06 Thread Fecher, Gerhard
件人: Wien <wien-boun...@zeus.theochem.tuwien.ac.at> 代表 Fecher, Gerhard <fec...@uni-mainz.de> 发送时间: 2017年7月6日 6:06 收件人: A Mailing list for WIEN2k users 主题: Re: [Wien] error in vorb continues... Does the error appear only during optimization or already in a regular scf cycle ? Is your structure correct ? (

Re: [Wien] Error information during installation of WIEN2k_16.1 on CentOs7

2017-07-06 Thread Fecher, Gerhard
there is nothing wrong do not mix up compilation errors (there is none) from ifort, cc, or linking with messages from make. Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the

Re: [Wien] 答复: 答复: error in vorb continues...

2017-07-06 Thread Fecher, Gerhard
ate folder and whether or not it converges'? Yang 发送自 Outlook<http://aka.ms/weboutlook> 发件人: Wien <wien-boun...@zeus.theochem.tuwien.ac.at> 代表 Fecher, Gerhard <fec...@uni-mainz.de> 发送时间: 2017年7月6日 7:09 收件人: A Mailing list for WIEN2k users 主题:

Re: [Wien] 答复: error in vorb continues...

2017-07-06 Thread Fecher, Gerhard
chem.tuwien.ac.at<mailto:wien-boun...@zeus.theochem.tuwien.ac.at>> 代表 Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> 发送时间: 2017年7月6日 6:06 收件人: A Mailing list for WIEN2k users 主题: Re: [Wien] error in vorb continues... Does the error appear only during optimization or alre

Re: [Wien] error in vorb continues...

2017-07-06 Thread Fecher, Gerhard
Does the error appear only during optimization or already in a regular scf cycle ? Is your structure correct ? (in most cases of questions here in the forum it is not) Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is

Re: [Wien] Convergence problem in mbj potential calculation

2017-06-14 Thread Fecher, Gerhard
mis.northwestern.edu Partner of the CFW 100% gender equity project, www.cfw.org/100-percent<http://www.cfw.org/100-percent> Co-Editor, Acta Cryst A On Jun 14, 2017 04:24, "Fecher, Gerhard" <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: Indeed, your compound

Re: [Wien] Convergence problem in mbj potential calculation

2017-06-14 Thread Fecher, Gerhard
eached the convergence >> criteria cc=0.0001 & ec=0.0001). >> >> My structure is correct as it is giving the same xrd as the experimental >> xrd. >> >> Should I just stop the simulation & starts with MSR1 mixing scheme by >> removing broyden files. >> &

Re: [Wien] Convergence problem in mbj potential calculation

2017-06-13 Thread Fecher, Gerhard
eV ( metallic ) :GAP :0.0126 Ry = 0.172 eV (provided you have a proper k-mesh) :GAP :0.0273 Ry = 0.371 eV (provided you have a proper k-mesh) :GAP : -9. Ry = -.eV ( metallic ) :GAP : -9. Ry = -.eV ( metallic ) with regards, On Tue, Jun 13, 2017 at 1

Re: [Wien] Convergence problem in mbj potential calculation

2017-06-13 Thread Fecher, Gerhard
potential calculation Dear Gerhard, Convergence criteria is : CC=0.0001, ec=0.0001 8 atoms , K-points = 125, rk max=7 with regards, On Tue, Jun 13, 2017 at 11:35 AM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: No answer possible, without kn

Re: [Wien] Ferromagnetic compound has non-magnetic moments

2017-06-13 Thread Fecher, Gerhard
what is the magnetic moment per Ni atom and the total magnetic moment per cell from the experiment ? (any reference available ?) what type of "magnetic" does the experiment find: ferro-, ferri-, antiferro-magnetic, spiral spin structure or any other noncollinear spin structure ? Was your

Re: [Wien] Convergence problem in mbj potential calculation

2017-06-13 Thread Fecher, Gerhard
No answer possible, without knowing your convergence criteria are they just too tight ? Number of k-points too small ? rkmax ? number of atoms ? Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never

Re: [Wien] Exchange energy

2017-06-06 Thread Fecher, Gerhard
As by itself, without any context, it seems to be rather useless Hmm, the lower the negative energy the higher its absolute value, seems to be mathematically correct |-5| > |-4|, isn't it. Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be

Re: [Wien] Energy ordering of the Ni-3d

2017-05-27 Thread Fecher, Gerhard
i, May 26, 2017 at 7:59 AM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: The point group of fcc Ni in the magnetized case (along the easy axis [111]) is D3d(C3i) (Shubnikov group 166.101), read for example A. P. Cracknell "Magnetism in Crystalline Ma

Re: [Wien] Energy ordering of the Ni-3d

2017-05-26 Thread Fecher, Gerhard
The point group of fcc Ni in the magnetized case (along the easy axis [111]) is D3d(C3i) (Shubnikov group 166.101), read for example A. P. Cracknell "Magnetism in Crystalline Materials", Pergamon Press, 1975 Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the

Re: [Wien] XC-functional Option 15 in version 16

2017-05-21 Thread Fecher, Gerhard
What just repeats last weeks the discussion with Fatima. Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the question is." Dr. Gerhard H.

Re: [Wien] installation error in lapwso and lapw1

2017-05-18 Thread Fecher, Gerhard
just a note, in the old thread it was mentioned: "Seems like the new mkl does not installs cluster libs under a non-commertial license." If I remember correctly, it (Clusterlibs, scalapack, and others) was just not installed automatically in some versions but I had to install it by hand.

Re: [Wien] Coult not reproduced results of CoSb3

2017-05-16 Thread Fecher, Gerhard
for exact exchange only. Then I looked to UG and in table 7.4 and 7.6 where I got the information for EX_xx and VX_xx only. So I considered EC_xx and VC_xx should be EC_NONE and VC_NONE, respectively. Regards Fatima On Mon, May 15, 2017 at 5:22 PM, Fecher, Gerhard <fec...@uni-mainz.de<m

Re: [Wien] Coult not reproduced results of CoSb3

2017-05-15 Thread Fecher, Gerhard
e where I have mistaken? My thermoelectric properties are quite different that reported in test case. Warm regards Fatima On Mon, May 15, 2017 at 12:24 PM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de><mailto:fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>

Re: [Wien] Coult not reproduced results of CoSb3

2017-05-15 Thread Fecher, Gerhard
ferent that reported in test case. Warm regards Fatima On Mon, May 15, 2017 at 12:24 PM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: If I read the code correctly then options 15, 25, and 26 should still be related to the functional of Engel and Vosco Ciao G

Re: [Wien] Error with NMR calculation

2017-05-15 Thread Fecher, Gerhard
just a question, why do you expect a metal to be a diamagnetic material ? please check your input files whether or not they are correct Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually

Re: [Wien] Coult not reproduced results of CoSb3

2017-05-15 Thread Fecher, Gerhard
If I read the code correctly then options 15, 25, and 26 should still be related to the functional of Engel and Vosco Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the

Re: [Wien] paramagnetic or diamagnetic

2017-05-09 Thread Fecher, Gerhard
spin-polarized calc. and an external field. >> Usually this part is paramagnetic. And then you have to see, which >> part dominates >> >> See also our NMR package. >> >> Am 08.05.2017 um 16:28 schrieb Fecher, Gerhard: >>> I am afraid that this questio

Re: [Wien] paramagnetic or diamagnetic

2017-05-08 Thread Fecher, Gerhard
nd that the non-magnetic phase is the most stable. so how can I know if the studied material is a paramagnetic or diamagnetic material? Thank you in advance 2017-05-08 8:06 GMT+02:00 Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>>: What distinguishes a paramagneti

Re: [Wien] could not reproduced Elastic properties for CsCl

2017-05-08 Thread Fecher, Gerhard
is ~4.1Ang. Let me do an another calculation and see how results are changing. Thank you Sir, Regards Fatima On Mon, May 8, 2017 at 11:13 AM, Fecher, Gerhard <fec...@uni-mainz.de<mailto:fec...@uni-mainz.de>> wrote: are you sure that this is CsCl ? I always thought that CsCl is simp

Re: [Wien] could not reproduced Elastic properties for CsCl

2017-05-07 Thread Fecher, Gerhard
are you sure that this is CsCl ? I always thought that CsCl is simple cubic with a=4.1 AA (and not fcc like NaCl), isn't it. Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest with you, is that you have never actually known what the

Re: [Wien] How to know the user id & password for graphical interface of wien2k

2017-04-01 Thread Fecher, Gerhard
You installed Wien2k You started W2web You had to give a user id and a password and typed it in And now you ask the forum about the user id and the password ? Thanks, its April 1st Ciao Gerhard DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy: "I think the problem, to be quite honest

Re: [Wien] Help about mBJ

2017-03-29 Thread Fecher, Gerhard
NO you cannot make such a conclusion, it would imply that always the method is correct that gives the largest band gap, this is nonsense. What is the band gap from experiment ? I guess it is not known because the compounds do not exist, do they ? Ciao Gerhard DEEP THOUGHT in D. Adams;

Re: [Wien] effective mass calculation from electronic band structure

2017-03-28 Thread Fecher, Gerhard
for this purpose I made my own k-lists by hand and used the energies from the band structure (case.spaghetti_ene); principally one can also use the case.energy files if one has a proper k-list. What to use may depend on the format you like to have for the method you use to analyse the E(k)

<    1   2   3   4   5   >