[Wien] Why the antiferromagnetic structure of Cr is still cubic

Dear wien2k users I have read that the antiferromagnetic ordering is accompanied with a little structural distortion as in the case of NiO. This latter is characterized by a distortion from NaCl to rhomboherdral structure; but for the chromide Cr the structure is still cubic when there is a magnet

[Wien] Why the antiferromagnetic structure of Cr is still cubic

Dear wien2k users I have read that the antiferromagnetic ordering is accompanied with a little structural distortion as in the case of NiO. This latter is characterized by a distortion from NaCl to rhomboherdral structure; but for the chromide Cr the structure is still cubic when there is a magnet

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thanks Dr Delamora for your reply Now I inderstand that and according to this raison "geometrical frustration" , is there any structural distortion in the planes (100) and (110) for the NaCl structure. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Al

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thanks again Dr Delamora for your explanation Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wi

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thanks again and I have read that the lifting of the Geometrical Frustration lead to Ferrimagnetism. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.a

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thanks Mr Khan Manuel for your contribution to enrich the subject of geometrical frustration and magnetism . After reading your contribution, I made the following remarks: * you have said:In the case of NiO the magnetic order reduces the symmetry to monoclinic, but what we know all is that the ma

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thank you Pr Stefaan Cottenier for your contribution According to what your are saying , the real magnetic structures for Chromide is tetragonal and orthorhombic and not cubic as it's mentionned in the userguide. I think the cubic magnetic structure fo Cr is given just to simply the calculation f

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thanks again Pr Stephan for the rich information Now , I understand the origins of the tetragonal and the orthorombic distortions for the cubic structure of Chromium Cr. I know how to add the spin-orbit coupling in the wien2k code to get the tertragonal structure , but I don't know how to take t

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thanks again Pr Staphan for the procedure and the article. Unfortunately I have no idea about the Spin Density Wave and I have to learn more about it. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie _

Re: [Wien] Why the antiferromagnetic structure of Cr is still cubic

Thanks again Mr Juan Manuel for the rich information Would you tell me to what extent this information is just: *The character of the unit cell distortion is determined by the orientation of the easy magnetization direction (EMD) of the compounds for binary RCo2 [6]. A close correlation between

[Wien] How to simulate the spin reorientation in the magnetic study

Dear wien2k users we know that the spin reorientation in the magnetic structure has an effect on the physicak properties of materials. How to simulate that using the wien2k code Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie

Re: [Wien] How to simulate the spin reorientation in the magnetic study

Thanks Mr Victor for your reply I have learned from the literature that spin-orbit coupling is involved in the spin orientation . Even for a ferromagnetic or a antiferromagnetic state , the spins have different orientations. How can we determine the orientation of the spins according to the use o

Re: [Wien] How to simulate the spin reorientation in the magnetic study

Thanks Mr Victor for the advise, but right now I need a special workshop that talk about my interesting. thanks again -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem.tuwie

Re: [Wien] How to simulate the spin reorientation in the magnetic study

Thanks Mr Victor for the clarification because I was wondering if I am not able to ask good questions. I am from Algeria and I am working on the antiferromagnetism in some materials and I wanted to know how to simulate magnetic phase transition with spin reorientation. Thank you for the link of

Re: [Wien] Constrained LDA calculations to find effective U

Dear all I have got result similar to that of Mr Gavin with the version 14.2 and using the script but I need to know how to change the imput files for another compound (for example : NiS) Best Regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie

[Wien] A question about the magnetic structure of a hexagonal crystal structure

Dear wien users For the purpose of determing the ground state magnetic structure of a NiAs structure, I have studied the three possibles magnetic structures, but for one of them I have found an othorhombic structure. Is it possible to get a such structure? Best regards -- Mr: A.Reggad Laboratoi

[Wien] Which notation for the spin direction in hexagonal structure

Dear Wien2k users I want to introduce the spin-orbit coupling in my calculation for a hexagonal structure. Would you inform me which notation should we use to introduce the spin direction : that of 3 indicesor that of 4 indices. NB: The spin direction for my case is the z direction Best regards

Re: [Wien] Which notation for the spin direction in hexagonal structure

Thanks Dr Xavier for your quick reply After typing the command : initso_lapw I got this error algerien1970@linux-1h0d:~/Bureau/NiS> initso_lapw The file NiS.in2c has been generated automatically >Please select the direction of the moment ( h k l ) (For R-lattice in R coordinat

Re: [Wien] Which notation for the spin direction in hexagonal structure

Hello again I have resolved the problem. Now I ovserved that as in the case of NiO , the spin-orbit coupling doesn't improve the gap for my compound. My question is: Which quantities are influenced by the spin-orbit coupling? Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Univers

[Wien] A new package for the Convergence Tests

Dear users of Wien2k Our group, because of the tedious task of convergence tests operation, has asked to our colleague Mr Benabdellah to develop a code that can do the job for us. The result is this beta version that needs improvement according to your remarks after trying it. You can download t

[Wien] A new package for the Convergence Tests

This the corrected link for the use of the package [Wien] A new package for the Convergence Tests -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus

Re: [Wien] A new package for the Convergence Tests

Sorry for the bad link This is the correct link Link Good Luck -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wie

[Wien] What's the difference between the spin-polarized and the non spin-polarized calculations

Dear wien2k users I have checked many discussions about the difference between the SP and the non SP calculations and I didn't find a sufficient explanation that removes the ambiguity . So I want to know the difference between the 2 calculations at the level of imput and results for a non magneti

Re: [Wien] What's the difference between the spin-polarized and the non spin-polarized calculations

Thank you Dr Gavin for your reply and also for your interesting for my questions. I have checked the 2 presentations but I didn't find what I look for . It's mentionned that in non spin-polarized calculation the spin-up density = the spin-down density . Which parameters are they monitored to make

Re: [Wien] What's the difference between the spin-polarized and the non spin-polarized calculations

Thank you Dr Gavin and Dr Gerhard for your answers. The were really very fruitful for me. We can consider that the non spin-polarized is a special case of the spin-polarized case. therefore we can use the equations of the general case to get the equations of the special case. So we can make B_ef

[Wien] Is there an effect of SO on the ground state

Dear wien2k users I want to determine the ground state between the non magnetic , ferromagnetic and different antiferromagnetic states and I want to know if the inclusion of spin-orbit coupling affect the ground state or not. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université

Re: [Wien] Is there an effect of SO on the ground state

Hello again I am waiting for an answer to my question . My question is about the effect of the inclusion of the Spin-Orbit Coupling on the ground state energy. I want to know if the SO affect the ground state energy also or It only causes the splitting of the degenerate state energies. Any info

Re: [Wien] Is there an effect of SO on the ground state

Thanks Dr pieper for the rich information within your answer. Perhaps I didn't formulated my question well . I am not interesting exactly to the ground state energy , but to the magnetic ground state. As I mentionned before , I want to determine the magnetic ground state from 5 configurations: nm,

Re: [Wien] Is there an effect of SO on the ground state

Thank you Dr Pieper for the detailed and fruitful information. For my compound with the hexagonal structure, I know beforehand the magnetic ground structure which is still hexagonal with magnetic order AF1, but for the 2 other magnetic structures for AF2 and AF3 , they become orthorhombic. Since

[Wien] What does it mean the commment " Unmatched " during SCF cycle

Dear Wien2k Users During the Volume optimization and after adding another point "-9 per cent " to get a good curve , I got this comment "Unmatched "even the calculation still running without error: The following show part of the SCF cycle with the comment "Unmatched " in cycle 5ETEST: .06445

Re: [Wien] What does it mean the commment " Unmatched " during SCF cycle

Thanks Dr Assmann for your quick reply I repetad the calculation and I same problem and I want to note that I modified the optimization script as follows: runsp_c_lapw -orb -o -ec 0.0001 Perhaps it's related to the use of "runsp_c_lapw " because it's the first time I used it I want to know if i

Re: [Wien] What does it mean the commment " Unmatched " during SCF cycle

Thanks your Dr Xavier for your answer It's just a mistake here the command line is as follows: runsp_c_lapw -orb -p -ec 0.0001 I note that the calculation ends without problem Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie _

[Wien] Difference between case.inst files for non spin-polarized and spin-polarized cases

Dear Wien2k Users I want to know if there is a difference between the case.inst files for the non spin-polarized and spin-polarized cases. We take the example of Ni: Ni Ar 3 3, 2,2.0 N 3, 2,2.0 N 3,-3,3.0 N 3,-3,1.0 N 4,-1,1.0 N 4,-1,1.0 N Best regards -- Mr: A.Reggad Laboratoire de Géni

Re: [Wien] Difference between case.inst files for non spin-polarized and spin-polarized cases

Thank your Prof Blaha for your quick answer I have used the command "instgen -ask " and I have found the difference between them as follows: Non spin-polarized case: Ni Ar 3 3, 2,2.0 N 3, 2,2.0 N 3,-3,2.0 N 3,-3,2.0 N 4,-1,1.0 N 4,-1,1.0 N Thank your Prof Blaha for your quick answer I h

Re: [Wien] Difference between case.inst files for non spin-polarized and spin-polarized cases

But when I have used the spin polarized atoms for a non spin-polarized calculation I have found a different energy . I have done a non spin polarized initialization using the expert mode and when comparing the results with two other calculations fm and afmI to determine the most stable magnetic

Re: [Wien] Difference between case.inst files for non spin-polarized and spin-polarized cases

I think there is a misunderstanding 1/ The problem lies in the expert mode of the graphical interface w2web. I think I have always to use the instgen -ask script before doing the initialization within the expert mode even for a non spin-polarized calculation. 2/ I didn't understand this statement

Re: [Wien] Difference between case.inst files for non spin-polarized and spin-polarized cases

Thank you Prof Blaha for the fruitful information and sorry for the misunderstanding Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.t

[Wien] Minimization before Volume Optimization ?

Dear Wien2k Users I want to do volume optimization for 5 structures (nm (hex);afmI(hex),afmII(orth) and afmIII(orth)) to determine the magnetic ground state. Should I do minimization of the atomic positions for the orhorhombic structures (because of the free parameter) before doing volume optimiz

[Wien] How to use mBJ method to determine the magnetic ground state

Dear Wien Users I have used the PBE+U method (with U=4.5 eV) to determine the magnetic ground structure which is AFM1 for 3d compound. I have succeded to get this magnetic ground state but the problem is that I have got a gap bigger than the expermimental one (0.1 and 1.3 eV). Now, I want to know

Re: [Wien] How to use mBJ method to determine the magnetic ground state

Dear Azam Do you think that the adding of SO coupling will reduce the value of the found gap (1.3 ev ) to the exeperimental value of 0.1 eV. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien m

Re: [Wien] Minimization before Volume Optimization ?

Hello again I am still waiting for any answer for my query . every effort will be helful Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://z

Re: [Wien] Minimization before Volume Optimization ?

Dear Gerhard Yes I did but my question is about the determination of the magnetic ground state . Is it necessary the minimization of the atomic positions for this purpose? Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie __

Re: [Wien] Minimization before Volume Optimization ?

Thank your Gerhard and Dobysheva for your answers * Do you think the ground state is the state with arbitrarily choosen positions of the atoms ? - No,of course bur I want to be sure * When you make a calculation you take _a_model_ instead of an actual sample. The choice of the model usually tak

Re: [Wien] Minimization before Volume Optimization ?

Now I want to know if a low symmetry structure (orthorhombic) derived from its high symmetry structure (hexagonal NiAs )has its atomic positions optimized (because they are created with sgroup program) or need to be optimized . Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Universi

Re: [Wien] Minimization before Volume Optimization ?

Thanks Mr pieper for the rich information . So , I have to minimize the atomic positions of the new low symmetry structure with respect to internal forces, or energy. Now the question is: Is there any difference between using the internal forces or the energy ? or it's the same ? -- Mr: A.Regga

[Wien] No convergence during Volume Optimization

Dear Wien Users Sometimes , It happens that during volume optimization , the calculation doesn't converge to the chosen energy criteria and jumps to the next point . In this case, should we delete this point and choose another point or we could accept this calculation. For example , with energy c

Re: [Wien] No convergence during Volume Optimization

Thank you Dr Pieper for your interesting to my questions. I have optimized the atomic positions before doing calculation. Tha thing that I didn't understand is that the convergence is reached for some points but not for others. For the "optimization notes " , there is no mention on the differenc

[Wien] What's the correct procedure for mBJ calculation

Dear Prof Tran In the page 14 of WIEN2k_lecture-notes_2013 , It's mentionned that the saving of the PBE calculation precedes the first mBJ initialization ( http://www.wien2k.at/reg_user/textbooks/WIEN2k_lecture-notes_2013/Exercises_13.pdf) while in the userguide ( http://www.wien2k.at/reg_user/tex

Re: [Wien] What's the correct procedure for mBJ calculation

Thank you Mr Polak for your answer Yes, you are about the userguide as a unique reference. I have followed the procedure of the first reference and I got results. Another question: I have got a gap value of 1.3 ev for my compound which is bigger than the experimental value (0.1 ev) with the d

Re: [Wien] What's the correct procedure for mBJ calculation

Hi Polack According to the reply of Prof Tran ( check this link http://zeus.theochem.tuwien.ac.at/pipermail/wien/2016-March/024447.html ) , we can modify the c parameter until we get the experimental value of gap. I used the default parameters and it gives wrong value . This value of gap I have

Re: [Wien] No convergence during Volume Optimization

Thank you Dr Pieper for the rich information. Yes I have a antiferromagnetic material and for the convergence problem ; the calculation jumps to the next point before reaching the energy criteria (foe example with energy criteria of 0.0001 , it jumps when reaching energy difference of 0.0008 ). I

Re: [Wien] No convergence during Volume Optimization

Thank you Prof Laurence for your answer As I said in my answer to Dr Pieper I have resolved the problem by using another cell parameters and with other points. This problem happens usually when doing volume optimization and I my procedure to resolve the problem was to delete the points with conv

Re: [Wien] What's the correct procedure for mBJ calculation

I have used the unmodified version of Beck-Johnson and It gives me a little value for the gap 0.047 eV. and I also used the other parametrizations and they give a similar value as the original version of mBJ. I want to know how to reproduce a little value of 0.1 eV. Best regards -- Mr: A.Regga

[Wien] How to get the experimental gap with mBj, PBE+U and EECE methods

Dear Wien Users I have a material with an experimental gap value of 0.1 ev and when using different methods I found values far from this value as follows: * mBJ method: - The original BJ parametrization : Eg=0.06 ev ---> There is no parameter to manipulate - The modified BJ para

Re: [Wien] How to get the experimental gap with mBj, PBE+U and EECE methods

Thank you Prof Tran for your answer You think that this value for gap value is enough . For the other quantity I am interested which is magnetic moment its value is good and doesn't change a lot for all the cases. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldou

Re: [Wien] How to get the experimental gap with mBj, PBE+U and EECE methods

To Prof Tran Since the unmodified version of Beck-Johnson allows me to get 0.06 ev which corresponds to the value of C parameter equal to the unity (01). I think that to increase of this value would improve the gap value to 1 ev. The question is: Which value for C parameter shall I use to get thi

Re: [Wien] How to get the experimental gap with mBj, PBE+U and EECE methods

Dear prof Tran I am sorry for the error of the gap value . The gap value that i want to reproduce is 0.1 ev . Thank you for your intersting -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list W

Re: [Wien] How to get the experimental gap with mBj, PBE+U and EECE methods

Thank you Prof Laurence for your precious notes. In my strudy and at the first step, I am intersting to 2 quantities: energy gap and the magnetic moment. I am using 3 methods which are the GGA+U, mBJ and EECE and I want to know the effect of the exchange energy on these 2 quantities through these

Re: [Wien] How to get the experimental gap with mBj, PBE+U and EECE methods

Thank you again Prof Tran and Prof Laurence for your notes and guidenesses. For the procedure to searching the best method giving the best gap I am looking for . Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie _

[Wien] Discrepancy in the simulation of the paramagnetic state

Dear Wien2k and Prof Blaha According to Prof Blaha ( https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg01133.html ) , in the paramagnetic state, the local magnetic moments still exists bur oriented in arbitrary directions. And according to this approach , there are infinite configura

Re: [Wien] Discrepancy in the simulation of the paramagnetic state

Thank you Prof Blaha for your quick answer. The Ni atom is 3d transition metal . But my question is about the simulation of the paramagnetic state. There are many people that considere that the paramagnetic state is the non-spin polarierd one and the magnetic moment is zero, but you say no and the

Re: [Wien] Discrepancy in the simulation of the paramagnetic state

Thank you All for your contributions to enrich the subject . >From the beginning, I was convinced with the statement of prof Blaha, But when you find in literature some work where they simulate the paramagnetic state by the non-spin-polarized calculation this causes a troube. On this base, I need

[Wien] How to include the localized d orbitals in the atomic spheres?

Dear wien users The orbital potential in DFT+U and EECE for localized orbitals will be included only inside the atomic spheres . To include these localized orbitals, we need a value for energy cut off larger than -6 Ry (maybe -2 Ry ). Is it the only way to include the localized orbitals in the a

Re: [Wien] How to include the localized d orbitals in the atomic spheres?

Sorry for my question I realized that the energy cut off determine the valence and the core states. The question is now as follows: Wen we use the maximum values for the Rmt such a way the spheres become touched. Does that guarantee that the 3d electrons are all inside atomic spheres? Best rega

Re: [Wien] How to include the localized d orbitals in the atomic spheres?

Thank you Prof Cottenier for your answer My question is made according to the following statement: " The DFT+U and EECE are applied only inside atomic spheres " What does it mean that and how to realize it ? Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - T

Re: [Wien] How to include the localized d orbitals in the atomic spheres?

Thank you Dr Pieper for the clarification. How much does the including of the localized d orbitals inside the atomic spheres improve the band gap and the magnetic moment ? And how can we include all the localized d orbitals inside the atomic spheres? Best regards Mr: A.Reggad Laboratoire de Géni

Re: [Wien] How to include the localized d orbitals in the atomic spheres?

Thank you Dr pieper again for the clarification Now, It's clear for me Another question is as follows: How to consider electrons (for example d electrons ) as localized or itinerant or half half or it's considered within the calculation? Best regards -- Mr: A.Reggad Laboratoire de Génie Phys

Re: [Wien] How to include the localized d orbitals in the atomic spheres?

Thank you Dr Assmann for the clarification Now , Can we consider a part of d electrons as localized and the other part as itinerant (for example half and half )? Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie __

Re: [Wien] How to include the localized d orbitals in the atomic spheres?

Thank you a lot Dr Pieper for your contribution I have already given the question to Dr Elias and I am wainting Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem

[Wien] Can we eliminate the correlation-exchange effect

Dear wien users I want to study the effect of correlation and exchange upon some properties separately and I want to know if it's possible to do that . i.e. eliminate the exchange effect and maintin the correlation effect and the inverse . NB: I know that we can't do it completly since the both

Re: [Wien] Can we eliminate the correlation-exchange effect

Hello again I am wainting for an answer to my question Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Algerie ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/

[Wien] Which formalism is adopted in the DFT+U method

Dear wien users We know all that there are two formalism for the DFT+U method : the Dudarev and the Liechtenstein formalisms. I want to know is adopted in the wien2k code for the DFT+U method. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse:

Re: [Wien] Which formalism is adopted in the DFT+U method

Thank you Fecher for your answer I think the both formalisms are not implemented in the wien2k package Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel: +213(0)561861963 Algerie _

[Wien] The self-interaction-correction still exists for U=0

Dear wien users I have some queries about the LDA+U method of Anisimov 1- Is the SIC version of LDA+U method the only one put by Anisimov ? because there are many works of Anisimov where there is no indication for the version of LDA+U method. 2- If we put U=0 , the self-interaction correction

Re: [Wien] The self-interaction-correction still exists for U=0

Hello again Now, to answer my question I need from the wien2k developpers to provide us the Hamiltonian of the LDA+U(SIC) introduced by Anisimov and implemented in the wien code. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTA

Re: [Wien] The self-interaction-correction still exists for U=0

Thank you Prof Blaha for your answer I have found this textbook presented by your person ( textbook

Re: [Wien] The self-interaction-correction still exists for U=0

Hello again I am wainting for an answer for my last question and I think Mr Novak is the most qualified to do that since he's whose implemented the approach LDA+U in the wien2k code. -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEF

Re: [Wien] The self-interaction-correction still exists for U=0

Thank you Dr Fecher for your reply 1- My question "why " is related to the fact normally DFT+U method with U=J=0 leads to the DFT method but for me despite the trying to control the convergence I got a little difference. I wanted to make sure that the fact to make U=0 and J=0 would lead the make t

[Wien] Convergence problem within the method GGA+U

Dear wien2k users To determine to magnetic ground state among three states , I have made three calculations for each state . The problem encontred is to get the convergence for all the case even with the increade of k points number and the energy convergence criteria. I have joined the SCF cycl

Re: [Wien] Convergence problem within the method GGA+U

Hello again I have deleted the following option of charge extrapolation because with it the convergence was worst clmextrapol_lapw if (-e NiS-afmI-optV_U0-ec1.clmup && ! -z NiS-afmI-optV_U0-ec1.clmup) then clmextrapol_lapw -up clmextrapol_lapw -dn endif Best regards -- Mr:

Re: [Wien] Convergence problem within the method GGA+U

Hello again When you choose a convergence criterion for energy equal to 0.1 , do you think it's normal that the scf cycle stops arbitrarily ; onetime stops when dE= 0.59 and another time stops when dE= 0.23 and so on. I didn't see such a situation in the previous versions Best regard

Re: [Wien] Convergence problem within the method GGA+U

Thank you Prof Blaha for your answer 1- "From the file on the web, you cannot determine the final convergence" I haven't been knowing that before and thanks before for this information and I realize that's the source of confusion. 2- "The numbers you quote are for the last but one iteration,

Re: [Wien] Convergence problem within the method GGA+U

Thanks again Dr Fecher for your answers - it can never stop precisely at the number you give That's what happened in my case but according to the statement of Prof Blaha we can't refer to the web screen printing and we have to check the case.daytime during every cycle Best regards -- Mr: A.

[Wien] Different convergences in Non-SP and SP methods

Dear wien users I have done the same calculation with 2 different and equivalent methods and they gave me different results as it's shown in the joined picture. 1- Non Spin Polarized calculation GGA using the command run_lapw 2- Spin Polarized calculation GGA using the command runsp_c_lapw This

Re: [Wien] Different convergences in Non-SP and SP methods

Here is the picture https://3.bp.blogspot.com/-nQPv4kYEP6E/WGTv3rUyu_I/B6s/XBiwlkWG7oIeG73NeliKhvwaTzqYHrlsQCEw/s640/E-V-Plot2.png -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel: +213(0)561861963 Algerie

Re: [Wien] Different convergences in Non-SP and SP methods

Thank you Prof Blaha for your reply I have used the default convergence criterium (-ec 0.0001) and I need to use the values to get the same curve. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel: +213(0)5618619

Re: [Wien] Different convergences in Non-SP and SP methods

Thank you Mr delamora for your reply You have got a 0.003Ry/Cr difference between Cr ferromagnetic and antiferromagnetic structures But my problem is that I have found a difference between 2 equivalent calculations for the same structure: non magnetic Best regards -- Mr: A.Reggad Laboratoire

Re: [Wien] Different convergences in Non-SP and SP methods

Thank you Dr Pieper for you answer According to the results with a difference of 0.001 Ry , can we considere the 2 calculations are the same ? Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel: +213(0)561861963 Alg

Re: [Wien] Different convergences in Non-SP and SP methods

Thank you again Prof Blaha for your answer Sorry for difference value. it's 0.0001 Ry. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel: +213(0)561861963 Algerie ___ Wien

[Wien] The results of the 2Doptimize package

Dear wien users I have optimized an hexagonal structure with the package 2Doptimize and I want to know in which file we find the optimal values of the volume V and the ratio c/a. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF

[Wien] LIBXC 3.0

You may download the compiled version (compiled with the compiler ifort ) from this link : https://drive.google.com/open?id=0B-QmTukdEFfQX29LcHRSbzktYjQ and you may follow this video for the installation of the new version of wien2k_2016.1 https://www.youtube.com/watch?feature=player_embedded&v=

Re: [Wien] LIBXC 3.0

You may download the compiled version (compiled with the compiler ifort ) from this link : https://drive.google.com/open?id=0B-QmTukdEFfQX29LcHRSbzktYjQ and you may follow this video for the installation of the new version of wien2k_2016.1 https://www.youtube.com/watch?feature=player_embedded&v=V

Re: [Wien] The results of the 2Doptimize package

I have found in which file I get the optimal values of the lattice parameters but the question is now what's the best fit order to use. There are big differences between the result of different fit orders # Order of fit 3 a0= 6.60816 , c0= 10.1455 bohr a0= 3.49

Re: [Wien] The results of the 2Doptimize package

Thank you Dr Guechi for your answer Now how to solve this problem Iteration 23 WSSR: 5.95105e-25 delta(WSSR)/WSSR : -0.343082 delta(WSSR) : -2.0417e-25 limit for stopping : 1e-05 lambda : 1.12698e+07 resultant parameter values a1 = -1536.32 a2

Re: [Wien] Tr : Boltztrap with symmetry operation "inversion"!

You can use the boltztrap code for every structure What's the problem of inversion symmetry ? you will find the following links beneficial for you http://wien2k-algerien1970.blogspot.com/2016/09/the-installation-and-execution-of.html http://wien2k-algerien1970.blogspot.com/2016/09/the-installat

Re: [Wien] The results of the 2Doptimize package

Thank you Prof Blaha for these guidelines . They will be of great benefit for me. Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel: +213(0)561861963 Algerie ___ Wien mai

Re: [Wien] The results of the 2Doptimize package

Dear again For the Optimize.job script as for the 2Doptimize.job script it's necessary to uncomment the following line for spin-polarized case : # recommended option: use charge extrapolation *clmextrapol_lapw --> # **clmextrapol_lapw* if (-e NiS-afmI-opt2D_E09.clmup && \ ! -z NiS-afmI

Re: [Wien] The results of the 2Doptimize package

Thanks De Ghechi for your answer But I didn't the statement to uncomment this line only for SP case according the guide of 2Doptimize.job from this link http://susi.theochem.tuwien.ac.at/reg_user/textbooks/2Doptimize.pdf Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn

Re: [Wien] The results of the 2Doptimize package

Thank you very much Dr Guechi for your help Best regards -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel: +213(0)561861963 Algerie ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.

Re: [Wien] The results of the 2Doptimize package

Thank you Prof Blaha This is what I did until I found in the guide of the 2Doptimize package to uncomment this line for the SP case it would be confusing for me. Thank you again -- Mr: A.Reggad Laboratoire de Génie Physique Université Ibn Khaldoun - Tiaret Adresse: BP 144 AL ATTAF AIN DEFLA Tel

  1   2   3   >