abacus-develop icon indicating copy to clipboard operation
abacus-develop copied to clipboard

magnetic calculation errors

Open matphy1 opened this issue 8 months ago • 3 comments

Describe the bug

I use abacus LST version to calculate magnetic systems, the wrong is as followes. attachment files:

1-sdw-U-1.5.zip

Image

Expected behavior

No response

To Reproduce

No response

Environment

No response

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).

matphy1 avatar May 05 '25 07:05 matphy1

@dyzheng could you have time to take a look?

mohanchen avatar May 05 '25 09:05 mohanchen

I have repeated your calculation using ABACUS v3.10 LTS with the supplied input files, and the SCF runs smoothly and finishes without errors. The output files are attached.

1-sdw-U-1.5-rerun.zip

Also, using LDA to calculate the magnetic properties confuses me. I recalculated them using PBE instead of LDA from ABACUS with recently released pseudopotentials and orbitals. The SCF converged in the second calculation (directory ending with 02). The calculated magnetic moment of the Ni atoms is essentially the same as that calculated by VASP 5.4.4. The magnetic moment for Ni1-Ni6 in the ABACUS calculation (atom 25-30 in the VASP calculation) is close to 0, and other Ni atoms retain the initial magnetic moment well.

6191-neworb-pbe.zip 6191-vasp.zip

ahxbcn avatar May 08 '25 09:05 ahxbcn

Okay, thank you so much for your help. I'll go check if it's a memory issue for abacus. Once again, thank you for your assistance. @ahxbcn

matphy1 avatar May 09 '25 18:05 matphy1