The initial and reported magnetic m_y moment is wrong in some subroutines

Bug #1623415 reported by Roberto Robles
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
Siesta
Status tracked in Trunk
4.1
Fix Released
High
Ramon Cuadrado
Trunk
Fix Released
High
Ramon Cuadrado

Bug Description

In noncollinear and spin-orbit calculations the convention for the m_y spin moment is not consistent through the code. Subroutines fill_dscf_from_atom_info (where the magnetic moment is initialized) and spnvec (used to calculate the magnetic moment) are not consistent with the convention used in the rest of the code. The subroutine moments.F calculates the correct magnetic moment.

In addition the actual initialization of the magnetization is not the one intended in the input file.

The bug will be present in any noncollinear or spin-orbit calculation where m_y=/0 (phi=/0). An example is attached where the behaviour is shown. This example proves which is the right convention.

Version: up to trunk_564 and 4.1

Related branches

Revision history for this message
Roberto Robles (roberto-robles) wrote :
description: updated
Changed in siesta:
assignee: nobody → Ramon Cuadrado (ramon-cuadrado)
Revision history for this message
Roberto Robles (roberto-robles) wrote :

The bug has been addressed in revision 567 of Ramon Cuadrado's branch.

Nick Papior (nickpapior)
Changed in siesta:
status: New → Confirmed
importance: Undecided → High
Nick Papior (nickpapior)
Changed in siesta:
milestone: none → 4.1-b2
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Remote bug watches

Bug watches keep track of this bug in other bug trackers.