[Pw_forum] Symmetry in PWSCF 1.3.0 - were there any bugs?

Paolo Giannozzi giannozz at nest.sns.it
Fri Jan 23 10:27:33 CET 2004


On Thursday 22 January 2004 21:41, Konstantin Kudin wrote:

> I seem to recall some discussion on the list about
> some symmetry issues

there is a rather lengthy discussion about symmetry issues
in the manual. I am quite confident that there is nothing
wrong with the symmetry-searching algorithm. If nothing of 
the following explains your problem, please post your data.

Paolo
---
\item {\em pw.x does not find all the symmetries you expected.}
      {\tt pw.x} determines first the symmetry operations (rotations)
      of the Bravais lattice; then checks which of these are
      symmetry operations of the system (including if needed
      fractional translations). This is done by rotating
      (and translating if needed) the atoms in the unit cell
      and verifying if the rotated unit cell coincides with
      the original one.

      You may miss symmetry operations because:
      \begin{enumerate}
      \item
      The number of significant figures in the atomic positions
         is not big enough. In {\tt PW/eqvect.f90}, the variable
         {\tt accep} is used to decide whether a rotation is a
         symmetry operation. Its current value ($10^{-5}$) is quite
         strict: a rotated atom must coincide with another atom to 5
         significant digits. You may change the value of {\tt accep}
         and recompile.
      \item
       They are not acceptable symmetry operations of the Bravais
         lattice. This is the case for C$_{60}$, for instance:
         the $I_h$ icosahedral group of C$_{60}$ contains 5-fold
         rotations that are incompatible with translation symmetry
      \item
       The system is rotated with respect to symmetry axis.
         For instance: a C$_{60}$ molecule in the fcc lattice will
         have 24 symmetry operations ($T_h$ group) only if the double
         bond is aligned  along one of the crystal axis; if C$_{60}$
         is rotated in some arbitrary way,  pw.x may not find any
         symmetry, apart from inversion.
      \item
       They contain a fractional translation that is incompatible
         with the FFT grid. Typical fractional translations are 1/2 or 1/3
         of a lattice vector. If the FFT grid is not divisible
         respectively by 2 or by 3, the symmetry operation will not
         transform the FFT grid into itself. Such symmetry operations
         are disabled to prevent problems with symmetrization.
         Note that if you change cutoff or unit cell volume, the
         automatically computed FFT grid changes. and this may explain
         changes in symmetry (and number of k-points as a consequence)
         for no apparent good reason (only if you have fractional
         translations in the system, though).
      \item
       A fractional translation, without rotation, is a symmetry
         operation of the system (meaning that the cell is actually
         a supercell). In this case, all symmetry operations containing
         fractional translations are disabled. The reason is that in this
         rather exotic case there is no (simple) way to select those
         symmetry operations forming a group (in the strict mathematical
         sense).
-- 
Paolo Giannozzi             e-mail:  giannozz at nest.sns.it
Scuola Normale Superiore    Phone:   +39/050-509876, Fax:-563513 
Piazza dei Cavalieri 7      I-56126 Pisa, Italy




More information about the Pw_forum mailing list