<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>Il giorno 21/set/2011, alle ore 12.10, Padmaja Patnaik ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; position: static; z-index: auto; "><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></div></blockquote><div><br></div><div>Mmh, not easy to say... if there is a magnetic anisotropy, then you need to choose also the correct orientation of the magnetization (see angle1 and angle2) vector on the impurity in order to be sure that the calculation converges to the ground state. On which atoms do you put the starting magnetization? (what are species 1 and 2?)</div><div>Are you using smearing?</div><div>It would be easier if you post (or send privately) your input file and the relvant part of the output.</div><div><br></div><div>GS</div><div><br></div><br><blockquote type="cite"><div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; position: static; z-index: auto; "><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></div>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline"></span>
</div>
<br></body></html>