<div dir="ltr"><div><div>Hi,<br><br></div>I have been seeking an optimized cell for my system lately. But I'm not entirely sure whether my calculation goes right or not. I have used crystal_sg to describe the atoms position. However, as it is known, my output file shows me only de cell parameters and atomic positions and crystal "units". My crystal is a "ferrite", whose symmetry is 227, which should correpond to fcc. Why don't I have a fcc as my output on scf calculations, for instace?<br></div><div><br></div><div>When do I know that my vc-relax calculation is right? I'm sending my input file (on vc-relax) for any mistakes, in any case.</div><div><br></div><div>Best Regards.</div><div><br></div><div><br></div><div>Input file:</div><div><br></div><div> &CONTROL<br>                 calculation = 'vc-relax' ,<br>                restart_mode = 'from_scratch' ,<br>                      outdir = '/home/drumsris/Bulk-227/' ,<br>                  pseudo_dir = '/home/drumsris/Bulk-227/' ,<br>                      prefix = 'pwscf' ,<br>                     tstress = .true. ,<br>                     tprnfor = .true. ,<br>                       lorbm = .false.,<br>          wf_collect = .false.,<br>                       nstep = 20 ,<br> /<br> &SYSTEM<br>                       ibrav = 2,<br>                   celldm(1) = 15.8370388,<br>          space_group = 227,<br>                         nat = 3,<br>                        ntyp = 3,<br>                     ecutwfc = 71 ,<br>                     ecutrho = 564 ,<br>               nosym = .true. ,<br>                   input_dft = 'pbe' ,<br>                 occupations = 'smearing' ,<br>                     degauss = 0.0002 ,<br>                    smearing = 'methfessel-paxton' ,<br>                       nspin = 2 ,<br>   starting_magnetization(1) = -1,<br>   starting_magnetization(2) = 1,<br>   starting_magnetization(3) = 0,<br>          !tot_magnetization = -1,<br> /<br> &ELECTRONS<br>            electron_maxstep = 100,<br>             diagonalization = 'david' ,<br>              diago_thr_init = 1.0e-5 ,<br>            conv_thr = 1.D-5 ,<br>         mixing_beta = 0.7,<br>         mixing_mode = 'plain'<br> /<br> &IONS<br>        ion_dynamics = 'bfgs' ,<br>       pot_extrapolation = 'atomic' ,<br>                 upscale = 10D0 ,<br>               bfgs_ndim = 1 ,<br>        trust_radius_max = 0.8D0 ,<br>        trust_radius_min = 1.D-3 ,<br>        trust_radius_ini = 0.5D0 ,<br>                     w_1 = 0.01D0 ,<br>                     w_2 = 0.5D0 ,<br> /<br> &CELL<br>           cell_dynamics = 'bfgs' ,<br>               press = 0.D0 ,<br>             cell_factor = 2.0 ,<br>          press_conv_thr = 0.5D0 ,<br>         cell_dofree = 'z' ,            <br>/<br>ATOMIC_SPECIES<br>   Ni    58.69000  Ni.pbe-n-kjpaw_psl.1.0.0.UPF <br>   Fe    55.85000  Fe.pbe-spn-kjpaw_psl.1.0.0.UPF <br>    O    16.00000  O.pbe-n-kjpaw_psl.1.0.0.UPF <br>ATOMIC_POSITIONS crystal_sg<br>Fe    16d      0.62500    0.62500    0.62500 0 0 0<br>Ni    8a       0.00000    0.00000    0.00000 0 0 0<br>O     32e      0.37980    0.37980    0.37980 <br>K_POINTS {automatic}<br> 5 5 5 0 0 0<br></div><div><div><div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div><div>Ricardo Afonso<br></div>Student of Magnetism and Superconductivity Group<br></div><div>Federal University of Sao Carlos<br></div></div></div></div></div></div>
</div></div></div></div>