This thread is locked.Only browsing is available.
Top Page > Browsing
Berry phase calculation failure
Date: 2017/02/06 17:06
Name: John Sham   <>

Dear Prof.Ozaki

I just tried to perform berry-phase calculation towards NaCl system according to the mannual.

The initial SCF step is normally finished, but unfortunately when I use "polB" and file"

nacl.scfout"to calculate the polarization of the system, it failed and comes to an end like

"found an overcomplete basis set".

I also check the source code of "polB" and the related line of this error:

diagonalize Bloch matrix at k-points, k1 and k2

for (ik=0; ik<2; ik++){

if (recalc[ik]){



for (l=1; l<=fsize; l++){
if (ko[l]<OLP_eigen_cut){
printf("found an overcomplete basis set\n");

Maythis indicate that the eigenvector k cut of overlap matrix is supposed to be modified?

Would you please help me?

Many Thanks!


John Sham
Page: [1]

Re: Berry phase calculation failure ( No.1 )
Date: 2017/03/02 22:39
Name: T. Ozaki


I found that the problem happens when the optimization level of compiler is set
to higher level. Since 'Eigen_HH()' is the routine causing the problem, I modified
polB by directly applying #pragma optimization_level 1 to 'Eigen_HH()'.
I've confirmed that the problem seems to disappear without depending on the
optimization level of compiler users specify.

A patch is now available at



Page: [1]