Error in BZ sampling factor for GW sigma calculation

sponce · Post by **sponce** » Sun Sep 16, 2012 5:36 pm

Dear all,

While performing a Sigma calculation with abinit 6.15.2 I got the following error message:

Subroutine Unknown:0:ERROR
 Not able to give unique norm to order vectors
You should increase factor (currently, factor is = 70)

This error was raised while computing the 27th (out of 30) qpoints:

Code: Select all

     1)     0.00000000E+00  0.00000000E+00  0.00000000E+00       0.00694
     2)     1.66666667E-01  0.00000000E+00  0.00000000E+00       0.04167
     3)     3.33333333E-01  0.00000000E+00  0.00000000E+00       0.04167
     4)     5.00000000E-01  0.00000000E+00  0.00000000E+00       0.02083
     5)     1.66666667E-01  1.66666667E-01  0.00000000E+00       0.04167
     6)     3.33333333E-01  1.66666667E-01  0.00000000E+00       0.04167
     7)     5.00000000E-01  1.66666667E-01  0.00000000E+00       0.04167
     8)     3.33333333E-01  3.33333333E-01  0.00000000E+00       0.01389
     9)     0.00000000E+00  0.00000000E+00  2.50000000E-01       0.01389
    10)     1.66666667E-01  0.00000000E+00  2.50000000E-01       0.04167
    11)     3.33333333E-01  0.00000000E+00  2.50000000E-01       0.04167
    12)     5.00000000E-01  0.00000000E+00  2.50000000E-01       0.04167
    13)    -3.33333333E-01  0.00000000E+00  2.50000000E-01       0.04167
    14)    -1.66666667E-01  0.00000000E+00  2.50000000E-01       0.04167
    15)     1.66666667E-01  1.66666667E-01  2.50000000E-01       0.04167
    16)     3.33333333E-01  1.66666667E-01  2.50000000E-01       0.04167
    17)     5.00000000E-01  1.66666667E-01  2.50000000E-01       0.04167
    18)     3.33333333E-01  3.33333333E-01  2.50000000E-01       0.01389
    19)    -1.66666667E-01  3.33333333E-01  2.50000000E-01       0.04167
    20)    -3.33333333E-01  5.00000000E-01  2.50000000E-01       0.04167
    21)    -1.66666667E-01  5.00000000E-01  2.50000000E-01       0.04167
    22)    -3.33333333E-01 -3.33333333E-01  2.50000000E-01       0.01389
    23)     0.00000000E+00  0.00000000E+00  5.00000000E-01       0.00694
    24)     1.66666667E-01  0.00000000E+00  5.00000000E-01       0.04167
    25)     3.33333333E-01  0.00000000E+00  5.00000000E-01       0.04167
    26)     5.00000000E-01  0.00000000E+00  5.00000000E-01       0.02083
    27)     1.66666667E-01  1.66666667E-01  5.00000000E-01       0.04167

 Subroutine Unknown:0:ERROR
 Not able to give unique norm to order vectors
You should increase factor (currently, factor is = 70)

    28)     3.33333333E-01  1.66666667E-01  5.00000000E-01       0.04167
    29)     5.00000000E-01  1.66666667E-01  5.00000000E-01       0.04167
    30)     3.33333333E-01  3.33333333E-01  5.00000000E-01       0.01389

The input file used for this calculation is join to this post.

I looked at the code and apparently in the routine 56_recipspace/m_bz_mesh.F90:

Code: Select all

 ! WARNING ! This could not be sufficient if tested k1 - k2 has lower
 ! components than smallestlen. The factor 10 is giving us a security margin.
 factor = 10*(int(one/smallestlen)+1)

One should raise this "10" factor to "20" ? Also the code is reporting a value for the variable "factor" of 70 in my case while it should be 90.
The smallestlen beeing 1/8 the factor should be 10*(8+1) = 90

Kind regards,

Samuel

EDIT: I tested with a factor of factor = 15*(int(one/smallestlen)+1) but I still have the exact same problem (except that the error indicate a factor value of 105)