<!DOCTYPE html>

<html>

  <head>

    <meta http-equiv="content-type" content="text/html; charset=UTF-8">

  </head>

  <body>

    <p>Dear QE users and developers,</p>

    <p>I am currently trying to reproduce the construction of an

      irreducible k-point set as done by QE. <br>

      For this, I set "verbosity = high" to get the symmetry operations

      printed in the output file. <br>

      I start from a uniform k-point mesh. Then, using the same symmetry

      operations as QE, I transform every k-point and fold it back in

      the first Brillouin zone.<br>

      If the resulting k-point falls on another k-point of the uniform

      grid, it is NOT irreducible. <br>

      In this manner, as also described by BlÃ¶chl et al (Phys. Rev. B <b>49</b>,

      16223, 1994) I find the set of irreducible kpoints.<br>

    </p>

    <p>My code agrees with QE for a simple structure (fcc crystal tested

      and verified) but I have problems with a more complicated case

      (the 2D magnet FGT).<br>

      In this example, 6 symmetry operations are found (see attached

      QE-output file).<br>

      Starting from a 3x3x3 uniform grid, the irreducible set of kpoints

      - according to QE - contains 7 points. However, I find 12

      irreducible k-points.<br>

    </p>

    <p>First, please note, that every point found by QE is also

      contained in my set. But I find additional points which (according

      to QE) should be related by some symmetry operation. By looking at

      the weights, I could figure out which kpoints should belong

      together. <br>

      For instance: According to QE, the kpoints [1/3, 0, 0] and [2/3,

      0, 0] are equivalent, as well as [0, 0, 1/3] and [0, 0, 2/3]

      should be equivalent too. I recognized that all the extra points

      could be transformed into each other by translating the lattice.

      However, applying all the symmetry operations from the QE output

      file (these are exclusively rotations and not translations), I

      cannot transform these points into each other. You might try for

      yourself.</p>

    <p>So the question that I would like to ask is: Are there any

      "hidden" symmetry operations which are not explicitly printed in

      the output file? Could fractional translations be the reason? Is

      it maybe related to differences between point group and space

      group? Any other hints to what I am missing?</p>

    Thank you! Your help would be highly appreciated!<br>

    <br>

    Best,<br>

    Lukas

    <p></p>

  </body>

</html>