Dear all <div>I am trying to calculate metallic oxide(CoO). However, I have some puzzle in writing vcrelax.in file.</div><div>As the experience shows that CoO is insulator, my vcrelax.in file is showed in the following lines, I calculate it as insulator</div><div><div><b> &control</b></div><div><b>    calculation='vc-relax'</b></div><div><b>    restart_mode='from_scratch',</b></div><div><b>    prefix='CoO'</b></div><div><b>    pseudo_dir = '../../files/pseudo'</b></div><div><b>    outdir='./tmp'</b></div><div><b>    verbosity='high'</b></div><div><b>    etot_conv_thr = 1.0D-4</b></div><div><b>    forc_conv_thr = 1.0D-3</b></div><div><b> /</b></div><div><b> &system</b></div><div><b>    ibrav = 0,</b></div><div><b>    celldm(1) = 8.00,</b></div><div><b>    nat = 4,</b></div><div><b>    ntyp = 3,</b></div><div><b>    ecutwfc = 50.0</b></div><div><b>    ecutrho = 400.0</b></div><div><b>    nspin = 2</b></div><div><b>    tot_magnetization=0.00</b></div><div><b>    occupations = 'fixed',</b></div><div><b>    nbnd = 28</b></div><div><b> /</b></div><div><b> &electrons</b></div><div><b>    conv_thr =  1.d-10</b></div><div><b>    startingwfc ='file'</b></div><div><b>    startingpot ='file'</b></div><div><b> /</b></div><div><b> &ions</b></div><div><b>    ion_dynamics = 'bfgs'</b></div><div><b> /</b></div><div><b> &cell</b></div><div><b>    cell_dynamics = 'bfgs'</b></div><div><b> /</b></div><div><b>ATOMIC_SPECIES</b></div><div><b> Co1  58.933194  co_pbesol_v1.2.uspp.F.UPF</b></div><div><b> Co2  58.933194  co_pbesol_v1.2.uspp.F.UPF</b></div><div><b> O    15.999     O.pbesol-n-kjpaw_psl.0.1.UPF</b></div><div><b>ATOMIC_POSITIONS {crystal}</b></div><div><b>Co1      0.00000000     0.00000000     0.00000000</b></div><div><b>Co2      0.50000000     0.50000000     0.50000000</b></div><div><b>O        0.25000000     0.25000000     0.25000000</b></div><div><b>O        0.75000000     0.75000000     0.75000000</b></div><div><b>CELL_PARAMETERS {alat}</b></div><div><b>   0.570726115   0.570726115   1.031099100</b></div><div><b>   0.570726115   1.031099100   0.570726115</b></div><div><b>   1.031099100   0.570726115   0.570726115</b></div><div><b>K_POINTS {automatic}</b></div><div><b> 6 6 6 0 0 0</b></div><div><b>HUBBARD {ortho-atomic}</b></div><div><b>U Co1-3d  6.7553</b></div><div><b>U Co2-3d  6.7553</b></div></div><div>however, the outfile shows that </div><div><b>highest occupied, lowest unoccupied level (ev):    14.3247   13.7975</b></div><div>I am puzzled, what happened, why is homo higher than lumo? is my calculation correect? Or should I regard the system as matel when calculate varelax?</div><div>Thank you very much</div>