
Thanks for the reply,<br>buT I am still thinking that should be a bug in the optimizer, since  I also tried to run my system like a metal and nothing gonna change. <br>If I change the optimization method I can reach easly the convergence.<br>Cheers, Piero  <br><br>----- Original Message -----<br>From: Giuseppe Mattioli <giuseppe.mattioli@mlib.ism.cnr.it><br>Date: Friday, March 13, 2009 11:21 am<br>Subject: Re: [Pw_forum] roblem with bfgs algorithm optimizer for spin  polarized system<br>To: PWSCF Forum <pw_forum@pwscf.org><br><br>> <br>> Dear Piero<br>> <br>> >  &ELECTRONS<br>> >  /<br>> The default beta mixing parameter is often too high to fit for <br>> many purposes.<br>> Try to use the keyword mixing_beta = 0.3 (e.g.) within the <br>> &electrons group.<br>> <br>>  >                  occupations = 'fixed' ,<br>> This is not a good choice (at least in my experience) when you <br>> deal with <br>> system in which the magnetization may change. Try to add a <br>> gaussian or a cold <br>> smearing.<br>> <br>> Yours<br>> <br>> Giuseppe<br>>  <br>> On Friday 13 March 2009 01:23:46 pc229@kent.ac.uk wrote:<br>> <br>> > Dear PWscf developers, <br>> > I am trying to  optimize, on the latest version of the <br>> espresso code<br>> > (v4.0.4 ),  the cell parameter and also the ionic coordinate <br>> of hematite. <br>> > Initially I treat hematite like an insulator using  the following<br>> > pseudo-potential :Fe.pbe-sp-van_ak.UPF for  Iron and O.pbe-<br>> van_ak.UPF for<br>> > oxygen (I attached the input at the bottom). To optimize  both <br>> atomic> coordinates and cell parameters I used bfgs algorithm . <br>> As you know<br>> > hematite is an open-shell system, and you need to define the spin<br>> > orientation in order to describe it like a ferro-magnetic or<br>> > anti-ferromagnetic system. I choose the anti-ferromagnetic <br>> one, that means <br>> > to have  the total magnetization equal to 0 every SCF cycle <br>> and absolute<br>> > magnetization ranging from 11-16 (this highly pseudo potential <br>> dependent).> Anyway, for the first SCF cycle the magnetization <br>> carried out is correct,<br>> > while after the first geometry optimization step hematite <br>> starts to lose<br>> > the spin and the energy  also dancing between two value not <br>> allowing the<br>> > SCF to converge. So attempting to managed this problem I try <br>> to run it like<br>> > a metal introducing  a methfessel-paxton  smearing. Nothing <br>> seems to<br>> > change, the SCF of the second optimization step doesn't <br>> converge as<br>> > previously. Then I tried to restart my calculation from the <br>> first point<br>> > coordinates and the first SCF step converged, whilst the next <br>> no as in the<br>> > previous cases . This is very strange and get me frustrated ! <br>> Then I run<br>> > the same calculation changing both optimization algorithm to <br>> damp for the<br>> > ionic coordinate and damp-w for the cell. In this way the SCFs <br>> after the<br>> > first optimization cycle seem to converge with the right <br>> magnetization and<br>> > the next as well. I did also a further trial, because I need <br>> to use the<br>> > bfgs algorithm to optimize my structure. I tried to run the <br>> same system on<br>> > a previous version of espresso (v.2.1.2), and treating it like <br>> a metal <br>> > (with methfessel-paxton  smering) and using bfgs for ionic <br>> coordinate and<br>> > damp for cell optimization (since bfgs wasn't implemented to <br>> optimize the<br>> > cell parameter in this version) the first SCF  reached the <br>> convergence and<br>> > the following as well. How can I solve it ?  I also read in <br>> the documents<br>> > related to PWscf that bfgs  for optimizing cell parameters (vc-<br>> relax) is<br>> > still in a experimental version. Would it be possible that It <br>> has some<br>> > problem with spin polarized systems?? How can I manage this <br>> problem? I need<br>> > desperately to use bfgs to optimize the hematite.<br>> ><br>> >   <br>> > Best Regards, Piero<br>> ><br>> ><br>> >  &CONTROL<br>> >                        title = 'opt_bulk_PBE_van3.inp' ,<br>> >                  calculation = 'vc-relax' ,<br>> >                 restart_mode = 'from_scratch' ,<br>> >                   wf_collect = .true. ,<br>> >                       outdir =<br>> > <br>> '/home/pc229/backup/Counts/PWscf/Fe2O3_bulk/PWSCF/Fe2O3_bulk/opt_bulk_PBE_v>an/' , pseudo_dir =<br>> > <br>> '/home/pc229/backup/Counts/PWscf/Fe2O3_bulk/PWSCF/Fe2O3_bulk/pseudo/' ,<br>> > prefix = 'opt_bulk_PBE_van3.inp' ,<br>> >  /<br>> >  &SYSTEM<br>> >                        ibrav = 5,<br>> >                    celldm(1) = 11.091428,<br>> >                    celldm(4) = 0.587237106,<br>> >                          nat = 10,<br>> >                         ntyp = 3,<br>> >                      ecutwfc = 50 ,<br>> >                        nosym = .true. ,<br>> >                         nbnd = 70,<br>> >                        nelec = 100,<br>> >                   tot_charge = 0.000000,<br>> >                  occupations = 'fixed' ,<br>> >                        nspin = 2 ,<br>> >    starting_magnetization(1) = 1.,<br>> >    starting_magnetization(2) = -1.,<br>> >    starting_magnetization(3) = 0.,<br>> >            tot_magnetization = 0 ,<br>> >  /<br>> >  &ELECTRONS<br>> >  /<br>> >  &IONS<br>> >                 ion_dynamics = 'bfgs' ,<br>> >  /<br>> >  &CELL<br>> >                cell_dynamics = 'bfgs' ,<br>> >                        wmass = 2001.0060 ,<br>> >                  cell_dofree = 'all' ,<br>> >  /<br>> > ATOMIC_SPECIES<br>> >    Fe   55.62000  Fe.pbe-sp-van_ak.UPF <br>> >   Fe1   55.62000  Fe.pbe-sp-van_ak.UPF <br>> >     O   15.99000  O.pbe-van_ak.UPF <br>> > ATOMIC_POSITIONS crystal <br>> >   Fe       0.144358141   0.144342247   0.144393006<br>> >   Fe1      0.355605889   0.355642696   0.355656118<br>> >   Fe1      0.644394003   0.644359012   0.644344959<br>> >   Fe       0.855641399   0.855655859   0.855605070<br>> >   O        0.556470421   0.943529173   0.250007053<br>> >   O        0.250006034   0.556470770   0.943529639<br>> >   O        0.943529367   0.250006711   0.556470944<br>> >   O        0.443532949   0.056467967   0.749992942<br>> >   O        0.749994148   0.443532454   0.056467520<br>> >   O        0.056467649   0.749993110   0.443532750    <br>> > K_POINTS automatic <br>> >   4 4 4   0 0 0<br>> ><br>> ><br>> > -- <br>> > Pieremanuele Canepa<br>> > Room 230<br>> > School of Physical Science <br>> > Ingram Building<br>> > University of Kent,<br>> >  Canterbury,<br>> >  Kent, CT2 7NH <br>> > United Kingdom<br>> > Cell: +44-7768685833<br>> > ID Skype: pieremanuele<br>> > ID msn: pieremanuele83@hotmail.com<br>> > --------------------------------------------------------<br>> <br>> -- <br>> ********************************************************<br>> - Article premier - Les hommes naissent et demeurent<br>> libres et �gaux en droits. Les distinctions sociales<br>> ne peuvent �tre fond�es que sur l'utilit� commune<br>> - Article 2 - Le but de toute association politique<br>> est la conservation des droits naturels et <br>> imprescriptibles de l'homme. Ces droits sont la libert�,<br>> la propri�t�, la s�ret� et la r�sistance � l'oppression.<br>> ********************************************************<br>> <br>>    Giuseppe <br>> Mattioli                            <br>>    CNR - ISTITUTO DI STRUTTURA DELLA <br>> MATERIA   <br>>    v. Salaria Km 29,300 - C.P. <br>> 10                <br>>    I 00016 - Monterotondo Stazione <br>> (RM)          <br>>    Tel + 39 06 90672836 - Fax +39 06 <br>> 90672316    <br>>    E-mail: <giuseppe.mattioli@ism.cnr.it><br>>