[Pw_forum] a question about MnSe (was Pw_forum Digest, Vol 43, Issue 37)

Stefano de Gironcoli degironc at sissa.it
Fri Jan 21 10:12:08 CET 2011


Dear Bongani Ngwenya,

Please DO NOT reply to a Digest message .. this messes up thread search.

If one is interested in a AFM solution it is correct to define 2 Mn 
species (Mn1 and Mn2) in the ATOMIC_SPECIES card and use the two labels 
in ATOMIC_POSITIONS to distinguish the two types.

stefano



On 01/21/2011 10:04 AM, Bongani Ngwenya wrote:
> I can see that you need to put Mn once in your ATOMIC SPECIES and use Mn
> twice to indicate ATOMIC POSITIONS than using Mn1 and Mn2.
>
> On Fri, Jan 21, 2011 at 10:34 AM,<pw_forum-request at pwscf.org>  wrote:
>
>> Send Pw_forum mailing list submissions to
>>         pw_forum at pwscf.org
>>
>> To subscribe or unsubscribe via the World Wide Web, visit
>>         http://www.democritos.it/mailman/listinfo/pw_forum
>> or, via email, send a message with subject or body 'help' to
>>         pw_forum-request at pwscf.org
>>
>> You can reach the person managing the list at
>>         pw_forum-owner at pwscf.org
>>
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of Pw_forum digest..."
>>
>>
>> Today's Topics:
>>
>>    1. Re: a question about MnSe (jia chen)
>>    2. Re: a question about MnSe (Lorenzo Paulatto)
>>    3. Re: a question about MnSe (mohnish pandey)
>>    4. Re: a question about MnSe (Prasenjit Ghosh)
>>
>>
>> ----------------------------------------------------------------------
>>
>> Message: 1
>> Date: Thu, 20 Jan 2011 13:37:16 -0500
>> From: jia chen<jiachen at princeton.edu>
>> Subject: Re: [Pw_forum] a question about MnSe
>> To: PWSCF Forum<pw_forum at pwscf.org>
>> Message-ID:
>>         <AANLkTinL+iV_PozddAO0KjEx1CqTALsDnYqLFU0MoLg1 at mail.gmail.com<AANLkTinL%2BiV_PozddAO0KjEx1CqTALsDnYqLFU0MoLg1 at mail.gmail.com>
>> Content-Type: text/plain; charset="iso-8859-1"
>>
>> Dear All,
>>
>> I also have a question about vc-relax. How do you know if your cutoff is
>> large enough to overcome Pulay stress, or you can do something else to get
>> rid of it in you calculations. Thanks.
>>
>> Best Wishes
>> jia
>>
>> On Wed, Jan 19, 2011 at 11:57 AM, mohnish pandey<mohnish.iitk at gmail.com
>>> wrote:
>>> Dear QE users,
>>>                             I am trying to do bulk calculations for MnSe
>> in
>>> rocksalt and wurtzite structure. Its has AFM ground state in rocksalt
>>> structure but after checking convergence I am finding "rocksalt"
>> structure
>>> to be energetically higher than wurtzite structure. I am pasting my
>> segments
>>> of input and output files. Can anybody please help me in figuring out the
>>> problem.
>>>
>>> *INPUT FOR ROCKSALT*
>>>
>>> &control
>>>      calculation = 'vc-relax'
>>>      restart_mode='restart',
>>>      wf_collect = .true.
>>>      prefix='mnse',
>>>      pseudo_dir = '/home/mohnish/mnse-rocksalt-bulk',
>>>      outdir='/home/mohnish/mnse-rocksalt-bulk',
>>>      tprnfor = .true.,
>>>      tstress=.true.
>>>   /
>>>   &system
>>>      ibrav=  0, celldm(1)= 10.3026, nat=  4, ntyp= 3,
>>>      ecutwfc = 40.0, ecutrho = 320.0,
>>>      occupations='smearing', smearing='gauss', degauss=0.01,
>>>      nspin=2,
>>>      starting_magnetization(1)= 0.0,
>>>      starting_magnetization(2)= 0.5,
>>>      starting_magnetization(3)=-0.5,
>>>   /
>>>   &electrons
>>>     diagonalization='david'
>>>     mixing_mode = 'plain'
>>>     mixing_beta = 0.4
>>>     conv_thr = 1.0d-8
>>>     startingpot = 'file'
>>>     startingwfc = 'file'
>>>   /
>>> &IONS
>>>    ion_dynamics='bfgs'
>>>    trust_radius_max = 0.40
>>>    trust_radius_ini = 0.20
>>> /
>>> &CELL
>>>    cell_dynamics='bfgs',
>>> /
>>> ATOMIC_SPECIES
>>>    Se  78.960 Se.pbe-van.UPF
>>>    Mn1 54.938 Mn.pbe-sp-van_mit.UPF
>>>    Mn2 54.938 Mn.pbe-sp-van_mit.UPF
>>> CELL_PARAMETERS
>>> 0.50 0.50 1.00
>>> 0.50 1.00 0.50
>>> 1.00 0.50 0.50
>>> ATOMIC_POSITIONS {crystal}
>>>   Se  0.25 0.25 0.25
>>>   Se  0.75 0.75 0.75
>>>   Mn1 0.0  0.0  0.0
>>>   Mn2 0.5  0.5  0.5
>>> K_POINTS {automatic}
>>> 6 6 6 0 0 0
>>>
>>> *OUTPUT FOR ROCKSALT*
>>>
>>>   bfgs converged in   8 scf cycles and   4 bfgs steps
>>>       (criteria: energy<  0.10E-03, force<  0.10E-02, cell<  0.50E+00)
>>>
>>>       End of BFGS Geometry Optimization
>>>
>>>       Final enthalpy =    -464.4688263706 Ry
>>> Begin final coordinates
>>>       new unit-cell volume =    533.55466 a.u.^3 (    79.06464 Ang^3 )
>>>
>>> CELL_PARAMETERS (alat= 10.30260000)
>>>     0.494008772   0.494008772   0.990589056
>>>     0.494008772   0.990589056   0.494008772
>>>     0.990589056   0.494008772   0.494008772
>>>
>>> ATOMIC_POSITIONS (crystal)
>>> Se       0.249999854   0.249999854   0.249999854
>>> Se       0.750000146   0.750000146   0.750000146
>>> Mn1      0.000000000   0.000000000   0.000000000
>>> Mn2      0.500000000   0.500000000   0.500000000
>>> End final coordinates
>>>
>>> *INPUT FOR WURTZITE*
>>>
>>> &control
>>>      calculation = 'vc-relax',
>>>      restart_mode='restart',
>>>      verbosity = 'high'
>>>      wf_collect = .true.
>>>      outdir='/home/mohnish/mnse-wurtzite-bulk',
>>>      pseudo_dir='/home/mohnish/mnse-wurtzite-bulk'
>>>      prefix='mnse',
>>>      tstress = .true.,
>>>      tprnfor = .true.,
>>> /
>>>   &system
>>>      ibrav=  4, a = 4.178, c = 6.783,nat= 4, ntyp= 3,
>>>      ecutwfc =
>>> 40,ecutrho=320,occupations='smearing',degauss=0.01,smearing='gaussian',
>>>      nspin =2,starting_magnetization(1)=0.5,starting_magnetization(2)=
>>> -0.5,
>>> /
>>>   &electrons
>>>     diagonalization='david'
>>>     mixing_mode = 'plain'
>>>     mixing_beta = 0.4
>>>     conv_thr = 1.0d-8
>>>     startingpot = 'file'
>>>     startingwfc = 'file'
>>>   /
>>> &IONS
>>>    ion_dynamics='bfgs'
>>>    trust_radius_max = 0.40
>>>    trust_radius_ini = 0.20
>>> /
>>> &CELL
>>>    cell_dynamics='bfgs',
>>> /
>>> ATOMIC_SPECIES
>>>    Mn1 54.938 Mn.pbe-sp-van_mit.UPF
>>>    Mn2 54.938 Mn.pbe-sp-van_mit.UPF
>>>    Se  78.960 Se.pbe-van.UPF
>>> ATOMIC_POSITIONS (crystal)
>>>    Mn1  0.000000000  0.000000000000  0.000000000
>>>    Mn2  0.333333333  0.666666666667  0.500000000
>>>    Se   0.000000000  0.000000000000  0.345000000
>>>    Se   0.333333333  0.666666666667  0.845000000
>>>   K_POINTS (automatic)
>>>   6 6 4 0 0 0
>>>
>>> *OUTPUT FOR WURTZITE*
>>>
>>>   End of BFGS Geometry Optimization
>>>
>>>       Final enthalpy =    -464.4799667612 Ry
>>> Begin final coordinates
>>>       new unit-cell volume =    681.36082 a.u.^3 (   100.96726 Ang^3 )
>>>
>>> CELL_PARAMETERS (alat=  7.89527578)
>>>     1.001414763   0.000000000   0.000000000
>>>    -0.500707382   0.867250625   0.000000000
>>>     0.000000000   0.000000000   1.594103933
>>>
>>> ATOMIC_POSITIONS (crystal)
>>> Mn1      0.000000000   0.000000000  -0.014999731
>>> Mn2      0.333333333   0.666666667   0.485592670
>>> Se       0.000000000   0.000000000   0.359333481
>>> Se       0.333333333   0.666666667   0.860073579
>>> End final coordinates
>>>
>>>
>>>
>>>
>>>
>>> --
>>> Regards,
>>> MOHNISH,
>>> -----------------------------------------------------------------
>>> Mohnish Pandey
>>> Y6927262,5th Year dual degree student,
>>> Department of Chemical Engineering,
>>> IIT KANPUR, UP, INDIA
>>> -----------------------------------------------------------------
>>>
>>> _______________________________________________
>>> Pw_forum mailing list
>>> Pw_forum at pwscf.org
>>> http://www.democritos.it/mailman/listinfo/pw_forum
>>>
>>>
>>
>> --
>> Jia Chen
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL:
>> http://www.democritos.it/pipermail/pw_forum/attachments/20110120/28ba57da/attachment.html
>>
>> ------------------------------
>>
>> Message: 2
>> Date: Fri, 21 Jan 2011 08:10:41 +0100
>> From: "Lorenzo Paulatto"<lorenzo.paulatto at impmc.upmc.fr>
>> Subject: Re: [Pw_forum] a question about MnSe
>> To: "PWSCF Forum"<pw_forum at pwscf.org>
>> Message-ID:<op.vpncv1mm5jfbqb at paulax>
>> Content-Type: text/plain; charset=utf-8; format=flowed; delsp=yes
>>
>> In data 20 gennaio 2011 alle ore 19:37:16, jia chen
>> <jiachen at princeton.edu>  ha scritto:
>>> I also have a question about vc-relax. How do you know if your cutoff is
>>> large enough to overcome Pulay stress, or you can do something else to
>>> get
>>> rid of it in you calculations. Thanks.
>> The vc-relax algorithm keeps number of plane waves constant (i.e. it
>> changes ecutwfc as ~ V^3/2) in order to avoid Pulay stress.
>>
>> best regards
>>
>>
>> --
>> Lorenzo Paulatto
>> post-doc @ IMPMC/UPMC - Universit? Paris 6
>> phone: +33 (0)1 44 27 74 89
>> www:   http://www-int.impmc.upmc.fr/~paulatto/
>>
>> previously (take note of the change!):
>> phd student @ SISSA&   DEMOCRITOS (Trieste)
>> phone: +39 040 3787 511
>> www:   http://people.sissa.it/~paulatto/
>>
>>
>> ------------------------------
>>
>> Message: 3
>> Date: Fri, 21 Jan 2011 13:00:16 +0530
>> From: mohnish pandey<mohnish.iitk at gmail.com>
>> Subject: Re: [Pw_forum] a question about MnSe
>> To: PWSCF Forum<pw_forum at pwscf.org>
>> Message-ID:
>>         <AANLkTikazpo9j53oUmWw_CFoLh=_KYzaVWLbPP72CECB at mail.gmail.com>
>> Content-Type: text/plain; charset="iso-8859-1"
>>
>> Dear Lorenzo !
>>                         Can you give me a direction for my question about
>> MnSe. I am really struck in it...
>>
>> On Fri, Jan 21, 2011 at 12:40 PM, Lorenzo Paulatto<
>> lorenzo.paulatto at impmc.upmc.fr>  wrote:
>>
>>> In data 20 gennaio 2011 alle ore 19:37:16, jia chen
>>> <jiachen at princeton.edu>  ha scritto:
>>>> I also have a question about vc-relax. How do you know if your cutoff
>> is
>>>> large enough to overcome Pulay stress, or you can do something else to
>>>> get
>>>> rid of it in you calculations. Thanks.
>>> The vc-relax algorithm keeps number of plane waves constant (i.e. it
>>> changes ecutwfc as ~ V^3/2) in order to avoid Pulay stress.
>>>
>>> best regards
>>>
>>>
>>> --
>>> Lorenzo Paulatto
>>> post-doc @ IMPMC/UPMC - Universit? Paris 6
>>> phone: +33 (0)1 44 27 74 89
>>> www:   http://www-int.impmc.upmc.fr/~paulatto/
>>>
>>> previously (take note of the change!):
>>> phd student @ SISSA&   DEMOCRITOS (Trieste)
>>> phone: +39 040 3787 511
>>> www:   http://people.sissa.it/~paulatto/
>>> _______________________________________________
>>> Pw_forum mailing list
>>> Pw_forum at pwscf.org
>>> http://www.democritos.it/mailman/listinfo/pw_forum
>>>
>>
>>
>> --
>> Regards,
>> MOHNISH,
>> -----------------------------------------------------------------
>> Mohnish Pandey
>> Y6927262,5th Year dual degree student,
>> Department of Chemical Engineering,
>> IIT KANPUR, UP, INDIA
>> -----------------------------------------------------------------
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL:
>> http://www.democritos.it/pipermail/pw_forum/attachments/20110121/b6eb6dcc/attachment-0001.htm
>>
>> ------------------------------
>>
>> Message: 4
>> Date: Fri, 21 Jan 2011 13:44:05 +0530
>> From: Prasenjit Ghosh<prasenjit.jnc at gmail.com>
>> Subject: Re: [Pw_forum] a question about MnSe
>> To: PWSCF Forum<pw_forum at pwscf.org>
>> Message-ID:
>>         <AANLkTinESV_Kk7=FUEDkM41J=ywHbi7wo_+ykj=-Pra3 at mail.gmail.com>
>> Content-Type: text/plain; charset=ISO-8859-1
>>
>> Dear Mohnish,
>>
>> One thing you may try to do is to change the value of starting
>> magnetization and redo the calculation.
>> Sometimes in case of magnetic systems (at least those which are
>> ferromagnetic), the code gets stuck into some local minima. Starting
>> from a different point might help you to reach the correct ground
>> state.
>>
>> Prasenjit
>>
>> On 21 January 2011 13:00, mohnish pandey<mohnish.iitk at gmail.com>  wrote:
>>> Dear Lorenzo !
>>> ?? ? ? ? ? ? ? ? ? ? ? ?Can you give me a direction for my question about
>>> MnSe. I am really struck in it...
>>>
>>> On Fri, Jan 21, 2011 at 12:40 PM, Lorenzo Paulatto
>>> <lorenzo.paulatto at impmc.upmc.fr>  wrote:
>>>> In data 20 gennaio 2011 alle ore 19:37:16, jia chen
>>>> <jiachen at princeton.edu>  ha scritto:
>>>>> I also have a question about vc-relax. How do you know if your cutoff
>> is
>>>>> large enough to overcome Pulay stress, or you can do something else to
>>>>> get
>>>>> rid of it in you calculations. Thanks.
>>>> The vc-relax algorithm keeps number of plane waves constant (i.e. it
>>>> changes ecutwfc as ~ V^3/2) in order to avoid Pulay stress.
>>>>
>>>> best regards
>>>>
>>>>
>>>> --
>>>> Lorenzo Paulatto
>>>> post-doc @ IMPMC/UPMC - Universit? Paris 6
>>>> phone: +33 (0)1 44 27 74 89
>>>> www: ? http://www-int.impmc.upmc.fr/~paulatto/
>>>>
>>>> previously (take note of the change!):
>>>> phd student @ SISSA ?&  ?DEMOCRITOS (Trieste)
>>>> phone: +39 040 3787 511
>>>> www: ? http://people.sissa.it/~paulatto/
>>>> _______________________________________________
>>>> Pw_forum mailing list
>>>> Pw_forum at pwscf.org
>>>> http://www.democritos.it/mailman/listinfo/pw_forum
>>>
>>>
>>> --
>>> Regards,
>>> MOHNISH,
>>> -----------------------------------------------------------------
>>> Mohnish Pandey
>>> Y6927262,5th Year dual degree student,
>>> Department of Chemical Engineering,
>>> IIT KANPUR, UP, INDIA
>>> -----------------------------------------------------------------
>>>
>>> _______________________________________________
>>> Pw_forum mailing list
>>> Pw_forum at pwscf.org
>>> http://www.democritos.it/mailman/listinfo/pw_forum
>>>
>>>
>>
>>
>> --
>> PRASENJIT GHOSH,
>> Assistant Professor,
>> IISER Pune,
>> First floor, Central Tower, Sai Trinity Building
>> Garware Circle, Sutarwadi, Pashan
>> Pune, Maharashtra 411021, India
>>
>> Phone: +91 (20) 2590 8203
>> Fax: +91 (20) 2589 9790
>>
>>
>> ------------------------------
>>
>> _______________________________________________
>> Pw_forum mailing list
>> Pw_forum at pwscf.org
>> http://www.democritos.it/mailman/listinfo/pw_forum
>>
>>
>> End of Pw_forum Digest, Vol 43, Issue 37
>> ****************************************
>>
>
>
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110121/f4b3b5b9/attachment-0001.htm 


More information about the Pw_forum mailing list