<div dir="ltr">Dear all:<br>    I am working on ZnO structural optimization using the LDA+U method. As in the journal of chemical physics 140, 121105 (2014), the U value on Zn d electrons was set to be 4.5 eV, calculated by vasp code, there was a good agreement with experimaental data in cell parameters by DFT+U than DFT method did.<br>
    the a, c and u in experimental values are 3.250, 5.210, 0.382, respectively. <br>    Calculated by VASP code, <br>    In DFT,the a, c and u values are 3.286, 5.295, 0.3807, respectively.<br>    In DFT+U,the a, c and u values are 3.243, 5.218, 0.3809, respectively.<br>
    It seems that the DFT+U correction works pretty well with the vasp code.<br>    <br>    But when using pwscf in vc-relax mode, I got the fellowing results:<br>    In DFT,the a, c and u values are 3.290, 5.305, 0.379, respectively.<br>
    In DFT+U,the a, c and u values are 3.292, 5.309, 0.379, respectively.<br>    It shows that in pwscf, the DFT+U correction does not much affect the results. also the datas show that there is significant difference in DFT +U results between the two codes.<br>
    <br>    My input file with pwscf attached to the following. I want to know where I went wrong in the pwscf input file.<br>    <br>    <br>    Any help will be appreciated!<br>    <br>    J.J. Gong<br>    Lanzhou University of technology, China<br>
<br><br>the pwscf input file    <br><br> &CONTROL<br>calculation = "vc-relax",<br>restart_mode = "from_scratch",<br>prefix = "ZnO-opt",<br>outdir = "/tmp/",<br>pseudo_dir = "~/pp",<br>
!wf_collect=.TRUE.<br>!etot_conv_thr = 1d-5<br>!forc_conv_thr = 1d-3<br>tstress = .true.<br>tprnfor = .true.<br>/<br>&SYSTEM<br>ibrav=0,<br>celldm(1)=1.88972613,<br>ecutwfc=55, ecutrho=550,<br>nat=4,ntyp=2,<br>lda_plus_U=.true., lda_plus_u_kind = 0, Hubbard_U(1)=4.5, Hubbard_U(2)= 0,<br>
/<br>&ELECTRONS<br>electron_maxstep = 80,<br>mixing_beta=0.5,<br>conv_thr=1D-8,<br>/<br>&ions<br>/<br>&cell<br>/<br><br><br>ATOMIC_SPECIES<br>  Zn  65.3900  Zn.pbe-van.UPF<br>   O  15.9994   o_pbe_v1.2.uspp.F.UPF<br>
  <br>CELL_PARAMETERS (alat=  1.88972613)<br>   3.292518870   0.000394478   0.000000646<br>  -1.640667267   2.850447665  -0.000000288<br>   0.000001154   0.000000426   5.310381081<br><br>ATOMIC_POSITIONS (crystal)<br>O        0.332685749   0.665993276   0.379582752<br>
O        0.667312132   0.334006494   0.879563839<br>Zn       0.333421586   0.666985007   0.000318706<br>Zn       0.666579130   0.333015022   0.500232704<br>K_POINTS (automatic)<br>4 4 3 1 1 0<br><br clear="all"><br>-- <br>
JiJun-Gong<br><br>Department of Applied Physics<br>Lanzhou University of Technology<br>Lanzhou<br>China
</div>