[Pw_forum] time-reversal symmetry broken
Timo Thonhauser
thonhauser at MIT.EDU
Wed Nov 21 15:05:42 CET 2007
Thanks Stefano, Nicola,
I already use the no_sym flag, which prevents the charge density from
being symmetrized with respect to all symmetry operations. However,
even with the no_sym flag, only half the k-points are generated (e.g.
for a 4x4x4 mesh the code generates 32 points and all the ones with
negative k_x are missing). Here Stefano's idea of supplying a k-point
list myself, would help. But, it is not that simple, is it? If inversion
symmetry and time reversal are present, the wave functions can be
chosen real (by real I mean that the wave functions only have a constant
arbitrary phase factor). So, if I break time reversal and give the
code my own, complete k-point list, does the code anywhere
still assume that the wave functions are "real"? Maybe breaking
inversion symmetry by hand could help?
Timo
Nicola Marzari wrote:
> Stefano Baroni wrote:
>
>> Timo: I do not have a solution handy, and others may have deeper ideas.
>> I think the problem is that when simmetrizing the charge density the
>> code always assumes e(k)=e(-k) and psi(-k)=psi(k)*, and so it never uses
>> more than 1/2 the BZ. One thing you may try is to pass explicitly all
>> the k-points, and wish that the code is not so smart as to remove 1/2 of
>> them ... Otherwise, I do not see any other options than hacking the code
>> (but, again, others may have smarter ideas). SB
>>
>
> Timo - what about the no_sym flag ? I thought that that would allow
> to force pwscf not to try to symmetrize the points (e.g. if
> one wanted to use a Baldereschi point, such as 1/4 1/4 1/4) - since
> it does not attempt to symmetrize anything, it might not use
> time reversal (rambling hypothesis - I'm sure SdG or PG know better,
> as always the case :-) ).
>
> As a side note - once k and -k are treated separately, even
> if the Hamiltonian is invariant for time reversal,
> in principle it's not psi(-k)=psi(k)* , but you might
> find more generally that psi(-k) = exp (i alpha) psi(k)* ,
> where alpha is random.
>
>
> nicola
>
>
> ---------------------------------------------------------------------
> Prof Nicola Marzari Department of Materials Science and Engineering
> 13-5066 MIT 77 Massachusetts Avenue Cambridge MA 02139-4307 USA
> tel 617.4522758 fax 2586534 marzari at mit.edu http://quasiamore.mit.edu
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
--
______________________________________________________________________
Dr. TIMO THONHAUSER
Department of Materials Science and Engineering (617) 253 6026
MIT, 13-4084, 77 Massachusetts Avenue thonhauser at mit.edu
Cambridge, MA 02139-4307, USA www.mit.edu/~thonhaus
______________________________________________________________________
More information about the Pw_forum
mailing list