<span class="Apple-style-span" style="border-collapse: collapse; font-family: arial; font-size: 13px; ">Dear PWscf developers, <br>I am trying to  optimize, on the latest version of the espresso code (v4.0.4 ),  the cell parameter and also the ionic coordinate of hematite. <br>Initially I treat hematite like an insulator using  the following pseudo-potential :Fe.pbe-sp-van_ak.UPF for  Iron and O.pbe-van_ak.UPF for oxygen (I attached the input at the bottom). To optimize  both atomic coordinates and cell parameters I used bfgs algorithm .<br>As you know hematite is an open-shell system, and you need to define the spin orientation in order to describe it like a ferro-magnetic or anti-ferromagnetic system. I choose the anti-ferromagnetic one, that means  to have  the total magnetization equal to 0 every SCF cycle and absolute magnetization ranging from 11-16 (this highly pseudo potential dependent).  <br>Anyway, for the first SCF cycle the magnetization carried out is correct, while after the first geometry optimization step hematite starts to lose the spin and the energy  also dancing between two value not allowing the SCF to converge. <br>So attempting to managed this problem I try to run it like a metal introducing  a methfessel-paxton  smearing. Nothing seems to change, the SCF of the second optimization step doesn't converge as previously. Then I tried to restart my calculation from the first point coordinates and the first SCF step converged, whilst the next no as in the previous cases . This is very strange and get me frustrated !<br>Then I run the same calculation changing both optimization algorithm to damp for the ionic coordinate and damp-w for the cell. In this way the SCFs after the first optimization cycle seem to converge with the right magnetization and the next as well. <br>I did also a further trial, because I need to use the bfgs algorithm to optimize my structure. I tried to run the same system on a previous version of espresso (v.2.1.2), and treating it like a metal  (with methfessel-paxton  smering) and using bfgs for ionic coordinate and damp for cell optimization (since bfgs wasn't implemented to optimize the cell parameter in this version) the first SCF  reached the convergence and the following as well.  <br>How can I solve it ?  I also read in the documents related to PWscf that bfgs  for optimizing cell parameters (vc-relax) is still in a experimental version. Would it be possible that It has some problem with spin polarized systems??<br>How can I manage this problem? I need 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 = '/home/pc229/backup/Counts/<wbr>PWscf/Fe2O3_bulk/PWSCF/Fe2O3_<wbr>bulk/opt_bulk_PBE_van/' ,<br>                  pseudo_dir = '/home/pc229/backup/Counts/<wbr>PWscf/Fe2O3_bulk/PWSCF/Fe2O3_<wbr>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: <a href="mailto:pieremanuele83@hotmail.com" target="_blank" style="color: rgb(6, 88, 181); ">pieremanuele83@hotmail.com</a><br>------------------------------<wbr>--------------------------</span>