Skip to content

#7664 destroys the y-component of magnetic moment #7831

Description

@maki49

Describe the bug

Before #7664, this symmetry-off Mn3Sn case with in-plane 120° AFM converges correctly to AMAG=19.8. In contrast, it converges to x-colinear state with AMAG=15.2 after #7664, that is, the y-component of magmom collapses.
Mn3Sn-7664-bug.zip

#7664 intended to change the $m_y$ convention consistent with the textbook Pauli matrix. However, DensityMatrix actually stores conj(D)=conj(c)*c instead of D=c*conj(c) to fit the col-major nature of BLAS (and conj(D)=transpose(D) are used as row-major D everywhere), the easiest way to fix the bug is to use the conjugated convention individually at the non-BLAS operators like m_y, i.e. to revert the change of func_xyz_to_updown in #7664.

Expected behavior

/

To Reproduce

/

Environment

/

Additional Context

/

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    BugsBugs that only solvable with sufficient knowledge of DFTcollinear/non-collinear/SOCIssues related to SOC

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions