160 likes | 418 Views
Magnetic Materials in PWSCF. Tutorial session, Thursday September 29 2005. Ralph Gebauer. Spin polarized systems: The collinear magnetic case Exercise: Electronic structure and DOS of Ni Non-collinear magnetism: extension of the formalism Exercise: Magnetism in bcc-iron
E N D
Magnetic Materials in PWSCF Tutorial session, Thursday September 29 2005 Ralph Gebauer • Spin polarized systems: The collinear magnetic caseExercise: Electronic structure and DOS of Ni • Non-collinear magnetism: extension of the formalismExercise: Magnetism in bcc-iron • Non-collinear magnetism: Constraints using a penalty functionalExercise: Angle-dependence of magnetism in bcc-iron
Spin polarized systems So far in this course: we “neglected” spin of the electrons: Implicitly, we have assumed that spin-up and spin-down wavefunctions are the same: Advantage: we need only Nel /2 orbitals
Spin polarized systems For the case of spin-polarized systems, we drop this constraint and Use N orbitals:
Calculate the functional derivative in order to obtain the Hamiltonian: Spin polarized systems Look at how the total energy depends on the two charge densities:
Spin polarized systems We can also write the exchange-correlation energy as . Since we have: The Hamiltonian looks like:
Spin polarized systems Exercise: Use the spin-polarized formalism to calculate The electronic structure of Ni
Non-collinear magnetism The exchange-correlation energy depends only on the amountof magnetization, not on its direction in space. The axis for the spin-up and spin-down projection is completelyarbitrary and not related to the system’s geometry. This axis is UNIQUE all over space (collinear magnetism). It is possible to generalize the concept of magnetization to complex,non-collinear magnetic geometries.
With these spinors, we define the vector magnetization as: The Pauli matrices are defined as: Non-collinear magnetism We introduce two-component spinors:
Non-collinear magnetism As before, the exchange-correlation energy depends onlyon the modulus of m. The generalized Schrödinger equation reads now:
Non-collinear magnetism Exercise: Use the non-collinear program to calculate the magnetic structure of bcc-iron.
Non-collinear magnetism • It is possible to constrain the magnetic moments of some atoms • to some fixed value • to a fixed direction. This can be done using a penalty functional: Which results in an additional term in the Hamiltonian:
Non-collinear magnetism For λ→∞, the magnetic moment takes the value mfixed.In that case, Bpen is finite. Therefore the penalty energy E pen →0.
Non-collinear magnetism Exercise: Use the non-collinear program to constrainthe magnetic moments in bcc-iron, and plot a curve of the total energy as a function of the angle.