Dear Paolo,<br>  <br>    Thank you for your reply, your comment is very useful to me. I shall have a try to change some part of the code 'bp_c_phase.f90' to include the noncolin=.true. case.<br><br>Best regards<br>
<br>Xun-Wang<br><br><div class="gmail_quote">2009/3/28 Paolo Giannozzi <span dir="ltr"><<a href="mailto:giannozz@democritos.it">giannozz@democritos.it</a>></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im"><br>
On Mar 28, 2009, at 6:25 , Xun-Wang Yan wrote:<br>
<br>
>  nspin=1 and nspin=2 cases is supported,   there is not<br>
> noncolin=.true.  case.<br>
> how to modify the bp_c_phase.f90 to include the noncolin=.true.<br>
> case for electric<br>
> polarization calculation?<br>
<br>
</div>I know nothing about polarization and noncolinear calculations, but<br>
basically<br>
what you need in a polarization calculation is a sequence of scalar<br>
products<br>
between electronic states at different k-points. In the case of<br>
noncolinear<br>
calculations, electronic states are spinors and are expanded into plane<br>
waves with both spin-up and spin-down components. Coefficients 1 to npw<br>
are for spin-up plane waves, coefficients npwx+1 to npwx+npw are for<br>
spin-down plane waves. Scalar products are summed over the two sets.<br>
<br>
Paolo<br>
---<br>
<font color="#888888">Paolo Giannozzi, Democritos and University of Udine, Italy<br>
</font><div><div></div><div class="h5"><br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br>