Error: xhegvd failed when running cell-relax of Fe2 with USPP

[ieiue]

Describe the bug

Running cell-relax of bcc Fe2 with USPP, there was an error: xhegvd failed. After the first step of relaxation, the calculation killed. I changed the ks_solver from cg to dav, it was still killed.

Expected behavior

How to complete this calculation.

To Reproduce

Fe2_uspp.zip

Environment

ABACUS v.3.6.4

Additional Context

No response

Task list for Issue attackers (only for developers)

• Verify the issue is not a duplicate.
• Describe the bug.
• Steps to reproduce.
• Expected behavior.
• Error message.
• Environment details.
• Additional context.
• Assign a priority level (low, medium, high, urgent).
• Assign the issue to a team member.
• Label the issue with relevant tags.
• Identify possible related issues.
• Create a unit test or automated test to reproduce the bug (if applicable).
• Fix the bug.
• Test the fix.
• Update documentation (if necessary).
• Close the issue and inform the reporter (if applicable).

[ieiue]

Additionally, if 'basis_type' is set to 'lcao', the calculation continues until it produces an incorrect result, rather than stopping and reporting an error.

[WHUweiqingzhou]

@YuLiu98 could you share any update?

[YuLiu98]

Compared to QE, it seems to arise from the unstable process of init_after_vc for uspp. I will fix it later.