<div class="gmail_quote">On Wed, Jun 20, 2012 at 10:49 AM, Virk Naunidh <span dir="ltr"><<a href="mailto:naunidh.virk@epfl.ch" target="_blank">naunidh.virk@epfl.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Dear Users,<br>
<br>
We are trying to project the wavefunction for a series of slab systems with spin orbit using projwfc.x, however for a certain number of those systems the calculation did not work and we were presented with the following error message:<br>

<br>
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
     from d_matrix_so : error #         2<br>
     D_S (j=1/2) for this symmetry operation is not unitary<br>
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
</div></div><div style="direction:ltr;font-size:10pt;font-family:Tahoma"><br></div></blockquote><div><br></div><div>Dear Virk,</div><div>the epsilon used to test for unitarity is very small. Maybe the matrix is indeed unitary but some numerical noise make it looks like it's not. The fact that your system is a slab with plenty of vacuum, and has low symmetry may also contribute. The fact that it only happens for some systems suggest that it is caused by numerical noise. It may also suggest that your calculation is not completely converged wrt the wfc cutoff.</div>
<div><br></div><div>Try to open PP/src/d_matrix_so.f90 and at line 42 change</div><div><div>  REAL(DP), PARAMETER :: eps = 1.0d-9</div></div><div>to something a bit larger. I.e. 1.d-8</div><div><br></div><div>It may also be a good idea to go down to line 282 and following to have the variable "capel" printed on output, which would give you an idea of how much non-unitary the matrix is.</div>
<div><br></div><div>If these suggestions do not solve the problem, please come back with more details, if the do, please confirm.</div><div><br></div><div>bests</div><div> </div></div><div><br></div>-- <br><div>Lorenzo Paulatto IdR @ IMPMC/CNRS & Université Paris 6</div>
<div>phone: +33 (0)1 44275 084 / skype: paulatz</div><div>www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a></div><div>mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05</div>
<br>