[Pw_forum] constrained magnetization calculation to converge at the ferromagnetic state

Arles V. Gil Rebaza arvifis at gmail.com
Thu Aug 25 14:08:57 CEST 2011


Dear Sreekar Guddeti, if you want to calculate a constrained magnetization
with antiparallel direction between the atom 2 and atom 8, you must to put
the "starting_magnetization" tag with opposite value, for example:

  starting_magnetization(2)=  0.3091,
  starting_magnetization(8)= -0.3091,

best wishes

PhD stud. Arles V. Gil Rebaza
Instituto de Física de La Plata
Argentina.

2011/8/25 sreekar guddeti <colonel.sreekar at gmail.com>

> Respected Sir/Madam,
> 1. I would like to perform a constrained magnetization calculation for a
> system, whose unconstrained spin polarized calculation gives an
> antiferromagnetic state, to converge at the ferromagnetic state.
> 2. The system is a zigzag graphene nanoribbon(*
> http://en.wikipedia.org/wiki/Graphene***) and I would like to flip the
> atomic moment at one of the edge site and constrain both the edge site
> moments.
> 3. I calculated the antiferromagnetic edge site moment  using PDOS
> calculation. I used 'lambda' parameter to constrain the moments.
>
> 4.Questions: Is 3) the right way to do it? If yes, what should be value of
> lambda? I used lambda =1 and the scf output (given below) is disturbing... I
> would like to ask what does "constrained moment :   -20.000000" mean and why
> is the constrain energy so high.  If it is not the right way, I would be
> grateful if you could suggest an alternative.
>
> Thanking you,
> Yours sincerely,
> Sreekar Guddeti
> Undergraduate (5th year)
> Engineering Physics
> IIT Bombay
> India
>
> PS: The system namelist for scf calculation and the sample output are given
> below and the complete input file is attached.
>
>
> --------------------------------------
> ...
>
>  ==============================================================================
>      atom number   15 relative position :   -3.5000   6.0622   0.0000
>      charge :     0.983302
>      magnetization :          0.000256
>      magnetization/charge:    0.000260
>      constrained moment :     0.309100
>
>
>  ==============================================================================
>
>
>  ==============================================================================
>      atom number   16 relative position :   -3.0000   5.7735   0.0000
>      charge :     1.783674
>      magnetization :         -0.385046
>      magnetization/charge:   -0.215873
>      constrained moment :   -20.000000
>
>
>  ==============================================================================
> ...
> ...
>
> *constraint energy* (Ryd) =   6095.83663339
> ---------------------------------
>
>
>
> --------------------------------
> &SYSTEM
>                        ibrav = 12,
>                    celldm(1) = 4.608737,
>                    celldm(2) = 12.0,
>                    celldm(3) = 10.0,
>                    celldm(4) = -0.50,
>                          nat = 18,
>                         ntyp =10,
>     starting_magnetization(2)= 0.3091,
>     starting_magnetization(8)= 0.3091,
>                   occupations='smearing',
>                      smearing='gauss',
>                       degauss=0.01,
>                         nspin=2,
>     constrained_magnetization='atomic'
>                       lambda = 1
>                      ecutwfc = 50.D0 ,
>                       report = 2
>
>  /
>
>
> --
> Sreekar Guddeti
>
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>


-- 
###--------->   Arles V.   <---------###
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20110825/01a1aa1e/attachment.html>


More information about the users mailing list