<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Thanks a lot GS. <br></span></div><div><span>I was waiting for some suggestions. As u have advised I need to study more about collinear and non-collinear formalism. About the system I am working on -- SiC doped with Cr. I am taking various supercells and Cr as a dilute impurity, doing relaxed and unrelaxed calcualtions for magnetic properties. Hope this is enough and would expect some suggestions for the magnetic moment problem.</span></div><div><br><span></span></div><div><span>Thanks and regards,</span></div><div>Padmaja Patnaik<br>Research Scholar<br>Dept of Physics<br>IIT Bombay<br>Mumbai, India<br></div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font face="Arial" size="2"><hr
size="1"><b><span style="font-weight: bold;"></span></b></font>Message: 2<br>Date: Tue, 20 Sep 2011 12:28:58 +0200<br>From: Gabriele Sclauzero <<a ymailto="mailto:gabriele.sclauzero@epfl.ch" href="mailto:gabriele.sclauzero@epfl.ch">gabriele.sclauzero@epfl.ch</a>><br>Subject: Re: [Pw_forum] Doubt in spin-orbit coupling results<br>To: Padmaja Patnaik <<a ymailto="mailto:padmaja_patnaik@yahoo.co.uk" href="mailto:padmaja_patnaik@yahoo.co.uk">padmaja_patnaik@yahoo.co.uk</a>>, PWSCF Forum<br> <<a ymailto="mailto:pw_forum@pwscf.org" href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br>Message-ID: <<a ymailto="mailto:BE7B21AF-2135-4CB9-9A66-6128CBB6A4FE@epfl.ch" href="mailto:BE7B21AF-2135-4CB9-9A66-6128CBB6A4FE@epfl.ch">BE7B21AF-2135-4CB9-9A66-6128CBB6A4FE@epfl.ch</a>><br>Content-Type: text/plain; charset="iso-8859-1"<br><br>Dear Padmaja Patnaik,<br><br> in order to include
spin-orbit interaction in the KS equations (lspinorb=.true.) one needs to use a formalism for noncollinear magnetization (noncolin=.true.). This means that there is no preferential axis for the magnetization, which is a vector field, (m_x(r),m_y(r),m_z(r)) in this case.<br>Therefore, in general one cannot define "spin up" and "spin down" (unless the ground state of the system is collinear, e.g. m_x=m_y=0 for each r). That's why the code does not give two distinct components of the DOS in the noncollinear formalism.<br><br>What you should do depends on what you are looking for... first I would suggest that you try to understand the difference between the collinear and noncollinear formalisms, in order to be able to interpret the results.<br><br>About the total magnetic moment, I can't say anything because you didn't say a word about your system.<br><br><br>GS<br><br>Il giorno 12/set/2011, alle ore 07.06, Padmaja Patnaik ha scritto:<br><br>> Hi
All<br>> <br>> I am doing spin-orbit coupling calculations. To get spin-polarized calculations with spin-orbit coupling I mentioned in the input for scf,<br>> lspinorb=.true.,<br>> noncolin=.true.,<br>> starting_magnetization(1)=0.5,<br>> starting_magnetization(2)=0.5,<br>> After this scf calculations I am doing the nscf calculation and finally the density of states and partial density of states. What I expect is it should give me both up and down spin DOS separately but I didnt get that. The input, nspin=2 in the input and get spin polarized results with both up and down spin density of states given separately which can be plotted. Can anybody please suggest what to do? Another doubt is the total magnetic moment is found to be zero (with spin-orbit interaction) which was earlier (without spin-orbit interaction) found to be 2 BohrMagneton./cell. <br>> <br><br><br></div></div></div></body></html>