Hi > Is there some bug in solve_e.f90? nothing known. Try to reduce alpha_mix and see if the problem disappears. Are you using mkl libraries ? Try to use lapack instead Paolo