<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div>hello<br data-mce-bogus="1"></div><div>i want to apply electric field to bilayer graphene. i am confused why do not have my output fermi energy?<br data-mce-bogus="1"></div><div>input is :<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div> &CONTROL<br> calculation = 'scf' ,<br> restart_mode = 'from_scratch' ,<br> wf_collect = .true. ,<br> outdir = './out' ,<br> wfcdir = './out' ,<br> pseudo_dir = './' ,<br> verbosity = 'high' ,<br> etot_conv_thr = 1.0D-6 ,<br> forc_conv_thr = 1.0D-2 ,<br> tstress = .true. ,<br> tprnfor = .true. ,<br> lelfield=.true., <br> nberrycyc=3 ,<br> gdir=3 , <br>nppstr=1 ,<br> /<br> &SYSTEM<br> ibrav = 4,<br> celldm(1) = 4.511848446,<br> celldm(3) = 9.462870616,<br> nat = 4,<br> ntyp = 1,<br> ecutwfc = 30 ,<br> ecutrho = 180 ,<br> input_dft = 'vdW-DF' ,<br> occupations = 'fixed' ,<br> degauss = 0.02 ,<br> smearing = 'methfessel-paxton' , <br> <br> /<br> &ELECTRONS<br> electron_maxstep = 200,<br> conv_thr = 1.0D-6 ,<br> startingwfc = 'atomic+random' ,<br> mixing_beta = 0.3 ,<br> diagonalization = 'cg' ,<br>efield_cart(1) = 0.d0,<br> efield_cart(2) = 0.d0,<br> efield_cart(3) = 0.019447d0,<br> /<br><br>ATOMIC_SPECIES<br> C 12.00000 C.blyp-mt.UPF <br>ATOMIC_POSITIONS (crystal)<br>C 0.333333330 0.666666660 0.572649573 <br>C 1.000000000 0.000000000 0.572649573 <br>C 0.666666660 0.333333330 0.427350427 <br>C 1.000000000 0.000000000 0.427350427 <br>K_POINTS automatic <br> 12 12 1 0 0 0 <br><br data-mce-bogus="1"></div><div>please help me.<br data-mce-bogus="1"></div></div></body></html>