<div dir="ltr">Dear all,<div>                  I have been working on bismuth ferrite BiFeO3, from past few months.</div><div>I use NC pseudopotential from pseudodojo and made all the optimization with respect to cutoff energy, K-points before performing the calculations.</div><div>The band calculation shows no bandgap (it has a bandgap of about 2.67 eV), hence calculations like Berry phase, dielectric and Turbo Lanczos and others cannot be performed as it is predicted to be a metal.</div><div>I solved the bandgap issue with the implementation of Hubbard potential (up to U= 5 eV for Fe), but this does not solve the issue completely as few of the calculations do not have the implementation of Hubbard potential in them. I wanted to know whether I am making any mistake or is there any workaround so that I can perform the dielectric property calculations for my material. I am pasting the scf input file if it is of any use. Thanks in advance </div><div><br></div><div>/* Input file*/</div><div>&CONTROL<br>    calculation = "scf"<br>restart_mode = 'from _scratch"<br>    max_seconds =  8.64000e+04<br>    pseudo_dir  = "."<br>    tprnfor     = .TRUE.<br>prefix      = "BFO"<br>    outdir      = "./outdir"<br>    tstress     = .TRUE.<br>/<br><br>&SYSTEM<br>    a                         =  5.645<br>    cosab                     =  0.50000<br>    degauss                   =  1.00000e-02<br>    ecutrho                   =  9.60000e+02<br>    ecutwfc                   =  8.00000e+01<br>    ibrav                     = 5<br>    nat                       = 10<br>    nspin                     = 2<br>    ntyp                      = 3<br> lda_plus_u                =.true<br>    Hubbard_U(1)              =2.0<br>    occupations               = "smearing"<br>    smearing                  = "gaussian"<br>    starting_magnetization(1) =  2.00000e-01<br>    starting_magnetization(2) =  0.00000e+00<br>    starting_magnetization(3) =  0.00000e+00<br>/<br><br>&ELECTRONS<br>    conv_thr         =  1.00000e-09<br>    electron_maxstep = 200<br>    mixing_beta      =  4.00000e-01<br>    startingpot      = "atomic"<br>    startingwfc      = "atomic+random"<br>/<br><br>K_POINTS {automatic}<br> 8  8  8  0 0 0<br><br>ATOMIC_SPECIES<br>Fe     55.84500  Fe.UPF<br>Bi    208.98038  Bi.UPF<br>O      15.99940  O.UPF<br><br>ATOMIC_POSITIONS (crystal)<br>Fe       0.777130455   0.777130455   0.777130455<br>Fe       0.277129455   0.277129455   0.277129455<br>Bi       0.509340660   0.509340660   0.509340660<br>Bi       0.009340660   0.009340660   0.009340660<br>O        0.563857130   0.956239313   0.116400442<br>O        0.956239313   0.116400442   0.563857130<br>O        0.116400442   0.563857130   0.956239313<br>O        0.616401442   0.456238313   0.063856130<br>O        0.456238313   0.063856130   0.616401442<br>O        0.063856130   0.616401442   0.456238313<br></div><div>/* end of input file */</div><div><br></div><div><br></div><div>Regards,</div><div>Srihari </div><div>Research scholar</div><div>Department of Physics</div><div>Manipal Institute of Technology</div><div>Manipal</div></div>