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

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions