[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