Message of " Sub-Space-Matrix is not hermitian "

Problems running VASP: crashes, internal errors, "wrong" results.


Moderators: Global Moderator, Moderator

Message
Author
TING QIN
Newbie
Newbie
Posts: 5
Joined: Tue Jan 10, 2006 11:13 am

Message of " Sub-Space-Matrix is not hermitian "

#1 Post by TING QIN » Tue May 23, 2006 12:59 pm

Hi, I do calculations for transition metals with respect to primitive monoclinic structure containing 16 atoms. VASP runs in the cluster and uses 4 processors. The INCAR is listed below:
SYSTEM = A12
ENCUT = 330
EDIFF=0.000001
ISMEAR =1
SIGMA = 0.1
IBRION=2
ISIF=3
NSW=1000
EDIFFG = -0.01
AMIX = 0.02

During the electronic structure relaxation, it always display lots of warning messages “WARNING: Sub-Space-Matrix is not hermitian in DAV 4� and eventually the calculation fails with such error message “Error EDDDAV: Call to ZHEGV failed. Returncode = 47 315�

As seen, I set AMIX from the default value of 0.4 to 0.02, and the problem still exists.
What is the reason and how can I avoid such problem?
Many thanks.
Last edited by TING QIN on Tue May 23, 2006 12:59 pm, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Message of " Sub-Space-Matrix is not hermitian "

#2 Post by admin » Tue May 23, 2006 1:45 pm

there may be several reasons for this error, the most common are
1a) the input geometry was not reasonable (error occurs at the very first ionic step) or
1b) the last ionic relaxation step lead to an unreasonable geometry (compare the input and output geometries of the
last ionic relaxation steps). In that case it can be helpful to
--> switch to a different relaxation algorithm (IBRION-tag)
--> reduce the step size of the first step by setting POTIM explicitely
2) If this error shows up at all your calculations: The installation of the LAPACK on your machine was not done properly: use a different LAPACK, e.g. the LAPACK which is delivered with the code
(vasp.4.lib/lapack_double.o)
3) on some architectures (especially SGI) some Lapack
routines are not working properly. However, it is
possible to avoid the usage of the ZHEGV subroutine
by commenting the line #define USE_ZHEEVX in
davidson.F, subrot.F, and wavpre_noio.F and recompiling
VASP.
Last edited by admin on Tue May 23, 2006 1:45 pm, edited 1 time in total.

sridevi
Newbie
Newbie
Posts: 5
Joined: Tue Jan 24, 2006 3:13 pm

Message of " Sub-Space-Matrix is not hermitian "

#3 Post by sridevi » Tue May 30, 2006 2:47 pm

I got the same problem.
If you are using vasp 4.6. Try using "IALGO =48" in the INCAR file.
It solved the problem for me.
Last edited by sridevi on Tue May 30, 2006 2:47 pm, edited 1 time in total.

jcconesa
Newbie
Newbie
Posts: 34
Joined: Mon Jul 03, 2006 10:48 am
Location: Madrid

Message of " Sub-Space-Matrix is not hermitian "

#4 Post by jcconesa » Sun Jan 07, 2007 6:06 pm

Hi,
I got the same problem on a sgi machine when running the Cu benchmark. Changing to IALGO=48 just changed the problem; then I get the messages:
.............
lib-4201: UNRECOVERABLE library error.

Encountered during a direct access unformatted READ from unit 21.
Fortran unit 21 is connected to a direct unformatted unblocked file: "TMPCAR"
IOT trap

core dumped
.....................

These messages appeared after 20 cycles of RMM: and one line beginnig with
1 T= 2080. E= -.90209009E+04 E0=...

Please help to solve. Concerning previous admin comments I should say that in the davidson, subrot and wavpre_noio files the variable "USE_ZHEEVX" does not control whether ZHEGV is used or not; this seems to depend rather on the definition of the variable "gammareal".
Last edited by jcconesa on Sun Jan 07, 2007 6:06 pm, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Message of " Sub-Space-Matrix is not hermitian "

#5 Post by admin » Wed Jan 24, 2007 2:05 pm

1) if you set USE_ZHEEVX
DSYEVX instead of DSYEV is used if gamma_real
ZHEEVX ZHEEV all other cases.
so it does not only affect the gamma-real version.


2) please have a look to the FAQs of the manual:
(http://cms.mpi.univie.ac.at/vasp/vasp/node265.html,
"I am running VASP on a SGI Origin, and the ...)
set IWAVPR=10
Last edited by admin on Wed Jan 24, 2007 2:05 pm, edited 1 time in total.

jcconesa
Newbie
Newbie
Posts: 34
Joined: Mon Jul 03, 2006 10:48 am
Location: Madrid

Message of " Sub-Space-Matrix is not hermitian "

#6 Post by jcconesa » Mon Feb 19, 2007 5:27 pm

Hi,
I continue to have the previously reported problem found when trying to run the Cu benchmark on a sgi machine, i.e. the program aborts and the output ends with a series of lines such as

entering main loop
N E dE d eps ncg rms
rms(c)
WARNING: Sub-Space-Matrix is not hermitian in DAV 1, -18.497193968206293
WARNING: Sub-Space-Matrix is not hermitian in DAV 2, -106.6910638174717
WARNING: Sub-Space-Matrix is not hermitian in DAV 3, -3.4046873909742339
WARNING: Sub-Space-Matrix is not hermitian in DAV 4, -37.403094929979197
WARNING: Sub-Space-Matrix is not hermitian in DAV 5, 0.14502948098981785

followed by

Error EDDDAV: Call to ZHEGV failed. Returncode = 13 1 8

The earlier solution suggested by admin (suppressing the line #define USE_ZHEEVX in davidson.F, subrot.F, and wavpre_noio.F and recompiling VASP) does not work, i.e. the same error messages, and the same indication of ZHEGV failure, still appear. I may add now that the problem appears both with the lapack which comes with VASP and with a system-native lapack library.

The warnings given suggest that the problem actually appears at an earlier stage, in which a matrix is generated with inadequate values which make it nonhermitian, and consequently ZHEGV fails even if working correctly; the solution thus would not be to avoid using ZHEGV, but to avoid an incorrect generation of the said matrix.

Can someone give an idea to really solve the problem?
Last edited by jcconesa on Mon Feb 19, 2007 5:27 pm, edited 1 time in total.

konglt

Message of " Sub-Space-Matrix is not hermitian "

#7 Post by konglt » Tue Feb 20, 2007 3:01 pm

Please try if it works by adding "LSCALAPACK = .FALSE." in your INCAR.
Last edited by konglt on Tue Feb 20, 2007 3:01 pm, edited 1 time in total.

konglt

Message of " Sub-Space-Matrix is not hermitian "

#8 Post by konglt » Tue Feb 20, 2007 3:03 pm

Please try if it works by adding "LSCALAPACK = .FALSE." in your INCAR.
Last edited by konglt on Tue Feb 20, 2007 3:03 pm, edited 1 time in total.

jcconesa
Newbie
Newbie
Posts: 34
Joined: Mon Jul 03, 2006 10:48 am
Location: Madrid

Message of " Sub-Space-Matrix is not hermitian "

#9 Post by jcconesa » Tue Feb 20, 2007 11:04 pm

No; adding "LSCALAPACK = .FALSE." in INCAR makes no dfference, the problem continues the same.
Last edited by jcconesa on Tue Feb 20, 2007 11:04 pm, edited 1 time in total.

ntq1982

Message of " Sub-Space-Matrix is not hermitian "

#10 Post by ntq1982 » Sat Feb 21, 2009 6:11 am

I was successful to fix this problem by using IALGO=48 instead of IALGO=Default
Last edited by ntq1982 on Sat Feb 21, 2009 6:11 am, edited 1 time in total.

xiaomatao
Newbie
Newbie
Posts: 4
Joined: Fri Oct 13, 2006 12:04 pm

Message of " Sub-Space-Matrix is not hermitian "

#11 Post by xiaomatao » Fri Mar 06, 2009 9:15 am

unfortunately, when i set IALGO=48, the new warning is
WARNING in EDDRMM: call to ZHEGV failed, returncode = 6 314
how to solve this problem?
Last edited by xiaomatao on Fri Mar 06, 2009 9:15 am, edited 1 time in total.

User avatar
VASP001
Jr. Member
Jr. Member
Posts: 72
Joined: Tue Oct 27, 2009 2:56 pm
License Nr.: 176, upgraded to paw
Location: Hefei National Laboratory for Physical Sciences at the Microscale

Message of " Sub-Space-Matrix is not hermitian "

#12 Post by VASP001 » Fri Dec 24, 2010 2:03 am

I face the same problem ,but it doesn't work .Any suggestions ?
Last edited by VASP001 on Fri Dec 24, 2010 2:03 am, edited 1 time in total.

rpf

Message of " Sub-Space-Matrix is not hermitian "

#13 Post by rpf » Mon Apr 18, 2011 12:36 pm

I got the same problem and solution is easy:

I was using 32 cpu when I took the error and the error disappeared for 16 cpu.

cheers :)
Last edited by rpf on Mon Apr 18, 2011 12:36 pm, edited 1 time in total.

jxjaghl

Message of " Sub-Space-Matrix is not hermitian "

#14 Post by jxjaghl » Sat Aug 11, 2012 8:38 am

I have the same problem during my calculation, and set IALGO=48 makes no difference either.
Last edited by jxjaghl on Sat Aug 11, 2012 8:38 am, edited 1 time in total.

nmwkim

Message of " Sub-Space-Matrix is not hermitian "

#15 Post by nmwkim » Mon Oct 08, 2012 11:13 am

I got the same warning msgs a lot..
Fortunately after setting "ALGO = Fast" it disappeared.
It seems the same with IALGO = 48 .
<span class='smallblacktext'>[ Edited Mon Oct 08 2012, 12:06PM ]</span>
Last edited by nmwkim on Mon Oct 08, 2012 11:13 am, edited 1 time in total.

Locked