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

Strange stress and force variations in Fe16H O-site calculations

Open ieiue opened this issue 1 year ago • 2 comments

Describe the bug

I want to make cell-relax calculation of Fe16H O-site. I tried 2 options: 1. cell-relax after relax (to get the results faster); 2. only cell-relax. But it can't converge. And the change trends of TOTAL PRESSURE and LARGEST GRAD are so strange. ERROR: image

  1. cell-relax after relax: image image

  2. only cell-relax: image image

Expected behavior

No response

To Reproduce

Here are calculation files:

  1. cell-relax after relax: There are 2 STRU in this files: STRU_relax for the first only relax calculation, and STRU_cell-relax is the result STRU_ION_D of the first step after converged. Fe16H O-site relax-cell_relax.zip
  2. only cell-relax: Fe16H O-site cell_relax.zip

Environment

ABACUS v.3.6.0

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 avatar Apr 16 '24 10:04 ieiue

@pxlxingliang, could you have a look?

WHUweiqingzhou avatar Apr 18 '24 04:04 WHUweiqingzhou

I will test it with v3.9.0

dyzheng avatar Dec 31 '24 09:12 dyzheng