<DIV>
<DIV>Dear Sax developers:</DIV>
<DIV> I am now doing GW calculation with Sax-1.0.0 , When I test the code with examples (example3)by calculating the band gap of Silicon ,It was running well, I got a reasonable bandgap of 0.98eV with 8*8*8 k-mesh .</DIV>
<DIV> Problem came out during I calculting the bandgap of P3HT <<FONT color="#000000">poly(</FONT>3-hexylthiophene)> which has 1-D long polymer chain structure.As I increase the number of k-points ,The band gap decreases relatively . A 0.26eV band gap of which experimental value is 1.9eV was obtained with 12*1*1 k-mesh. <BR> Similar problem happens when calculating the band gap of polyacetylene. <BR> I am not sure if I wrongly selected the parameters of GW calcultion or SAX code does not fit calculating band gap of oranic semiconductor(or some bugs exist) . </DIV>
<DIV> Hope you can give me some recommendation.</DIV>
<DIV> My input files are as follows:</DIV>
<DIV> </DIV>
<DIV> <FONT color="#f00000">PW1.inp</FONT></DIV>
<DIV>&CONTROL<BR>title="P3HT"<BR>calculation='scf',<BR>restart_mode = 'from_scratch',<BR>tprnfor=.t.,<BR>PSEUDO_DIR='/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/',<BR>prefix='P3HT',<BR>/<BR> &SYSTEM<BR> ibrav = 0,<BR> celldm(1) = 15.0,<BR> nat = 50,<BR> ntyp = 3,<BR> ecutwfc = 30.0 , <BR> ecutrho = 180.0,<BR> nosym=.true.,<BR> nbnd = 120,<BR> /<BR> &ELECTRONS<BR> /<BR>CELL_PARAMETERS cubic<BR> 1.000000000 0.000000000 0.000000000<BR> 0.000000000 3.00000000 0.000000000<BR> 0.000000000 0.000000000 1.100000000<BR>ATOMIC_SPECIES<BR> C 12.0107 C.GG.UPF<BR> S 32.0660 S.pz-bhs.UPF<BR> H 1.00079 H.vbc.new.UPF<BR>ATOMIC_POSITIONS angstrom<BR>C 1.217895241 0.002248390 0.001143191<BR>C 1.765269676 1.286141466 0.000076174<BR>C 3.177952059 1.363894578 0.000132124<BR>C 3.787482281 0.100166278 0.001199236<BR>C 3.973235161 2.646483093 -0.001256180<BR>C 3.185858486 3.954894242 0.000112415<BR>C 4.090657661 5.191548144 -0.001076051<BR>C 3.282609183 6.487323099 0.000637591<BR>C 4.093987701 7.782833685 -0.001508102<BR>C 3.182660520 9.009024821 0.000878701<BR>S 2.551626394 -1.144435997 0.002672455<BR>H 4.650002196 2.646646813 0.871941232<BR>H 4.647483056 2.646548858 -0.876429889<BR>H 2.527397896 3.995646753 0.883065789<BR>H 2.524171981 3.996054016 -0.880430210<BR>H 4.757007903 5.165517889 0.878170567<BR>H 4.754501263 5.165922255 -0.882226937<BR>H 2.617837552 6.490529999 0.881640276<BR>H 2.613018269 6.490269877 -0.876737817<BR>H 4.759799962 7.808734653 0.876661741<BR>H 4.755523982 7.808731675 -0.882894271<BR>H 2.531398653 9.009126194 0.887301738<BR>H 2.526293733 9.009243840 -0.881799573<BR>H 3.747595410 9.950782334 -0.000595651<BR>C 5.185008991 -0.284008523 0.000832983<BR>C 5.722371253 -1.568956760 0.000218186<BR>C 7.133923165 -1.645651506 0.000380522<BR>C 7.753277987 -0.383102429 0.000989843<BR>C 7.916880779 -2.934145369 -0.000021308<BR>C 7.122472324 -4.240039327 -0.000077630<BR>C 8.025195543 -5.479704033 -0.000153593<BR>C 7.220865197 -6.778318100 -0.000041689<BR>C 8.042955689 -8.066952378 -0.001012927<BR>C 7.149026332 -9.307014617 0.000136136<BR>S 6.517046608 0.861883019 0.001245678<BR>H 8.589929155 -2.935288364 0.874944728<BR>H 8.589925918 -2.934984712 -0.875010814<BR>H 6.462100374 -4.279709479 0.881715428<BR>H 6.462008825 -4.279716881 -0.881803967<BR>H 8.690588297 -5.452773545 0.880062815<BR>H 8.690429673 -5.452740713 -0.880481541<BR>H 6.554715753 -6.784987086 0.879680222<BR>H 6.553368869 -6.784354545 -0.878739903<BR>H 8.708859075 -8.085706520 0.877670127<BR>H 8.706861590 -8.085785792 -0.881197479<BR>H 6.495783018 -9.318176822 0.885110623<BR>H 6.496059032 -9.320576794 -0.885035711<BR>H 7.728684857 -10.239295533 0.001631764<BR>H 1.135198120 2.173715970 -0.001216108<BR>H 5.082477354 -2.450566116 -0.000504931<BR>K_POINTS crystal<BR> 2<BR> 0.250000000 0.500000000 0.500000000 0.500000000<BR> 0.750000000 0.500000000 0.500000000 0.500000000</DIV>
<DIV> </DIV>
<DIV><FONT color="#f00000">GW1.inp</FONT></DIV>
<DIV><FONT color="#000000"><input><BR><sax_options<BR> calculation_kind = "GW"<BR> start_from = "DFT"<BR> gw_integration_method = "SSHF"<BR> lorentzian_broadening = "0.01"<BR> nbandmin = "1"<BR> nbandmax = "120"<BR> nelec = "120.0"<BR> convert_from_pwscf="T"<BR> convert_from_pwscf_file="P3HT.export/index.xml"<BR> system_kind="1D"<BR> calc_polarizability="T"<BR> cutoff_polarizability="6.0"<BR> calc_w="T"<BR> cutoff_fock="30.0"<BR> cutoff_density="6.0" <BR> cutoff_vloc= "6.0"<BR> cutoff_vxc ="6.0"<BR>/></FONT></DIV>
<DIV><FONT color="#000000"><structure><BR> <direct<BR> a1 = " 15.0 , 0.0 , 0.0"<BR> a2 = " 0.0 , 45.0 , 0.0"<BR> a3 = " 0.0 , 0.0 , 16.5" /><BR></structure></FONT></DIV>
<DIV><FONT color="#000000"><atoms><BR><types ntypes="3"><BR>C "/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/C.GG.UPF" 30.0<BR>S "/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/S.pz-bhs.UPF" 30.0<BR>H "/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/H.vbc.new.UPF" 30.0<BR></types><BR><positions natoms="50" units="angstrom" alat="1.0"><BR>C 1.217895241 0.002248390 0.001143191<BR>C 1.765269676 1.286141466 0.000076174<BR>C 3.177952059 1.363894578 0.000132124<BR>C 3.787482281 0.100166278 0.001199236<BR>C 3.973235161 2.646483093 -0.001256180<BR>C 3.185858486 3.954894242 0.000112415<BR>C 4.090657661 5.191548144 -0.001076051<BR>C 3.282609183 6.487323099 0.000637591<BR>C 4.093987701 7.782833685 -0.001508102<BR>C 3.182660520 9.009024821 0.000878701<BR>S 2.551626394 -1.144435997 0.002672455<BR>H 4.650002196 2.646646813 0.871941232<BR>H 4.647483056 2.646548858 -0.876429889<BR>H 2.527397896 3.995646753 0.883065789<BR>H 2.524171981 3.996054016 -0.880430210<BR>H 4.757007903 5.165517889 0.878170567<BR>H 4.754501263 5.165922255 -0.882226937<BR>H 2.617837552 6.490529999 0.881640276<BR>H 2.613018269 6.490269877 -0.876737817<BR>H 4.759799962 7.808734653 0.876661741<BR>H 4.755523982 7.808731675 -0.882894271<BR>H 2.531398653 9.009126194 0.887301738<BR>H 2.526293733 9.009243840 -0.881799573<BR>H 3.747595410 9.950782334 -0.000595651<BR>C 5.185008991 -0.284008523 0.000832983<BR>C 5.722371253 -1.568956760 0.000218186<BR>C 7.133923165 -1.645651506 0.000380522<BR>C 7.753277987 -0.383102429 0.000989843<BR>C 7.916880779 -2.934145369 -0.000021308<BR>C 7.122472324 -4.240039327 -0.000077630<BR>C 8.025195543 -5.479704033 -0.000153593<BR>C 7.220865197 -6.778318100 -0.000041689<BR>C 8.042955689 -8.066952378 -0.001012927<BR>C 7.149026332 -9.307014617 0.000136136<BR>S 6.517046608 0.861883019 0.001245678<BR>H 8.589929155 -2.935288364 0.874944728<BR>H 8.589925918 -2.934984712 -0.875010814<BR>H 6.462100374 -4.279709479 0.881715428<BR>H 6.462008825 -4.279716881 -0.881803967<BR>H 8.690588297 -5.452773545 0.880062815<BR>H 8.690429673 -5.452740713 -0.880481541<BR>H 6.554715753 -6.784987086 0.879680222<BR>H 6.553368869 -6.784354545 -0.878739903<BR>H 8.708859075 -8.085706520 0.877670127<BR>H 8.706861590 -8.085785792 -0.881197479<BR>H 6.495783018 -9.318176822 0.885110623<BR>H 6.496059032 -9.320576794 -0.885035711<BR>H 7.728684857 -10.239295533 0.001631764<BR>H 1.135198120 2.173715970 -0.001216108<BR>H 5.082477354 -2.450566116 -0.000504931<BR></positions><BR></atoms></FONT></DIV>
<DIV><FONT color="#000000"><kmesh><BR><mesh nk="2 1 1" shift="0.5 0.5 0.5"/><BR></kmesh></FONT></DIV>
<DIV><FONT color="#000000"></input></FONT></DIV>
<DIV><FONT color="#f00000">PW2.inp</FONT></DIV>
<DIV>&CONTROL<BR>title="P3HT"<BR>calculation='scf',<BR>restart_mode = 'from_scratch',<BR>tprnfor=.t.,<BR>PSEUDO_DIR='/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/',<BR>prefix='P3HT',<BR>/<BR> &SYSTEM<BR> ibrav = 0,<BR> celldm(1) = 15.0,<BR> nat = 50,<BR> ntyp = 3,<BR> ecutwfc = 30.0 , <BR> ecutrho = 180.0,<BR> nosym=.true.,<BR> nbnd = 120,<BR> /<BR> &ELECTRONS<BR> /<BR>CELL_PARAMETERS cubic<BR> 1.000000000 0.000000000 0.000000000<BR> 0.000000000 3.00000000 0.000000000<BR> 0.000000000 0.000000000 1.100000000<BR>ATOMIC_SPECIES<BR> C 12.0107 C.GG.UPF<BR> S 32.0660 S.pz-bhs.UPF<BR> H 1.00079 H.vbc.new.UPF<BR>ATOMIC_POSITIONS angstrom<BR>C 1.217895241 0.002248390 0.001143191<BR>C 1.765269676 1.286141466 0.000076174<BR>C 3.177952059 1.363894578 0.000132124<BR>C 3.787482281 0.100166278 0.001199236<BR>C 3.973235161 2.646483093 -0.001256180<BR>C 3.185858486 3.954894242 0.000112415<BR>C 4.090657661 5.191548144 -0.001076051<BR>C 3.282609183 6.487323099 0.000637591<BR>C 4.093987701 7.782833685 -0.001508102<BR>C 3.182660520 9.009024821 0.000878701<BR>S 2.551626394 -1.144435997 0.002672455<BR>H 4.650002196 2.646646813 0.871941232<BR>H 4.647483056 2.646548858 -0.876429889<BR>H 2.527397896 3.995646753 0.883065789<BR>H 2.524171981 3.996054016 -0.880430210<BR>H 4.757007903 5.165517889 0.878170567<BR>H 4.754501263 5.165922255 -0.882226937<BR>H 2.617837552 6.490529999 0.881640276<BR>H 2.613018269 6.490269877 -0.876737817<BR>H 4.759799962 7.808734653 0.876661741<BR>H 4.755523982 7.808731675 -0.882894271<BR>H 2.531398653 9.009126194 0.887301738<BR>H 2.526293733 9.009243840 -0.881799573<BR>H 3.747595410 9.950782334 -0.000595651<BR>C 5.185008991 -0.284008523 0.000832983<BR>C 5.722371253 -1.568956760 0.000218186<BR>C 7.133923165 -1.645651506 0.000380522<BR>C 7.753277987 -0.383102429 0.000989843<BR>C 7.916880779 -2.934145369 -0.000021308<BR>C 7.122472324 -4.240039327 -0.000077630<BR>C 8.025195543 -5.479704033 -0.000153593<BR>C 7.220865197 -6.778318100 -0.000041689<BR>C 8.042955689 -8.066952378 -0.001012927<BR>C 7.149026332 -9.307014617 0.000136136<BR>S 6.517046608 0.861883019 0.001245678<BR>H 8.589929155 -2.935288364 0.874944728<BR>H 8.589925918 -2.934984712 -0.875010814<BR>H 6.462100374 -4.279709479 0.881715428<BR>H 6.462008825 -4.279716881 -0.881803967<BR>H 8.690588297 -5.452773545 0.880062815<BR>H 8.690429673 -5.452740713 -0.880481541<BR>H 6.554715753 -6.784987086 0.879680222<BR>H 6.553368869 -6.784354545 -0.878739903<BR>H 8.708859075 -8.085706520 0.877670127<BR>H 8.706861590 -8.085785792 -0.881197479<BR>H 6.495783018 -9.318176822 0.885110623<BR>H 6.496059032 -9.320576794 -0.885035711<BR>H 7.728684857 -10.239295533 0.001631764<BR>H 1.135198120 2.173715970 -0.001216108<BR>H 5.082477354 -2.450566116 -0.000504931<BR>K_POINTS crystal<BR> 2<BR> 0.000000000 0.000000000 0.000000000 0.500000000<BR> 0.500000000 0.000000000 0.000000000 0.500000000</DIV>
<DIV><FONT color="#f00000"></FONT> </DIV>
<DIV><FONT color="#f00000">GW2.inp</FONT></DIV>
<DIV><input><BR><sax_options<BR> calculation_kind = "GW"<BR> start_from = "DFT"<BR> gw_integration_method = "SSHF"<BR> lorentzian_broadening = "0.01"<BR> nbandmin = "1"<BR> nbandmax = "120"<BR> nelec = "120.0"<BR> convert_from_pwscf="T"<BR> convert_from_pwscf_file="P3HT.export/index.xml"<BR> system_kind="1D"<BR> calc_polarizability="F"<BR> cutoff_polarizability="6.0"<BR> calc_w="F"<BR> cutoff_fock="30.0" <BR> cutoff_density="6.0" <BR> cutoff_vloc= "6.0"<BR> cutoff_vxc ="6.0"<BR> calc_sigma_x="T"<BR> calc_sigma_c="T"<BR> calc_energies = "T"<BR> calc_sp_hmatrix = "T"<BR> sigma_nbmin = "1"<BR> sigma_nbmax = "80"<BR> diagonal = "T"<BR> sigma_first_order = "T"<BR> calc_QP_states = "T"<BR>/></DIV>
<DIV><structure><BR> <direct<BR> a1 = " 15.0 , 0.0 , 0.0"<BR> a2 = " 0.0 , 45.0 , 0.0"<BR> a3 = " 0.0 , 0.0 , 16.5" /><BR></structure></DIV>
<DIV><atoms><BR><types ntypes="3"><BR>C "/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/C.GG.UPF" 30.0<BR>S "/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/S.pz-bhs.UPF" 30.0<BR>H "/home/2/xiexiaoyin/sax-1.0.0/sax/examples/pseudos/H.vbc.new.UPF" 30.0<BR></types><BR><positions natoms="50" units="angstrom" alat="1.0"><BR>C 1.217895241 0.002248390 0.001143191<BR>C 1.765269676 1.286141466 0.000076174<BR>C 3.177952059 1.363894578 0.000132124<BR>C 3.787482281 0.100166278 0.001199236<BR>C 3.973235161 2.646483093 -0.001256180<BR>C 3.185858486 3.954894242 0.000112415<BR>C 4.090657661 5.191548144 -0.001076051<BR>C 3.282609183 6.487323099 0.000637591<BR>C 4.093987701 7.782833685 -0.001508102<BR>C 3.182660520 9.009024821 0.000878701<BR>S 2.551626394 -1.144435997 0.002672455<BR>H 4.650002196 2.646646813 0.871941232<BR>H 4.647483056 2.646548858 -0.876429889<BR>H 2.527397896 3.995646753 0.883065789<BR>H 2.524171981 3.996054016 -0.880430210<BR>H 4.757007903 5.165517889 0.878170567<BR>H 4.754501263 5.165922255 -0.882226937<BR>H 2.617837552 6.490529999 0.881640276<BR>H 2.613018269 6.490269877 -0.876737817<BR>H 4.759799962 7.808734653 0.876661741<BR>H 4.755523982 7.808731675 -0.882894271<BR>H 2.531398653 9.009126194 0.887301738<BR>H 2.526293733 9.009243840 -0.881799573<BR>H 3.747595410 9.950782334 -0.000595651<BR>C 5.185008991 -0.284008523 0.000832983<BR>C 5.722371253 -1.568956760 0.000218186<BR>C 7.133923165 -1.645651506 0.000380522<BR>C 7.753277987 -0.383102429 0.000989843<BR>C 7.916880779 -2.934145369 -0.000021308<BR>C 7.122472324 -4.240039327 -0.000077630<BR>C 8.025195543 -5.479704033 -0.000153593<BR>C 7.220865197 -6.778318100 -0.000041689<BR>C 8.042955689 -8.066952378 -0.001012927<BR>C 7.149026332 -9.307014617 0.000136136<BR>S 6.517046608 0.861883019 0.001245678<BR>H 8.589929155 -2.935288364 0.874944728<BR>H 8.589925918 -2.934984712 -0.875010814<BR>H 6.462100374 -4.279709479 0.881715428<BR>H 6.462008825 -4.279716881 -0.881803967<BR>H 8.690588297 -5.452773545 0.880062815<BR>H 8.690429673 -5.452740713 -0.880481541<BR>H 6.554715753 -6.784987086 0.879680222<BR>H 6.553368869 -6.784354545 -0.878739903<BR>H 8.708859075 -8.085706520 0.877670127<BR>H 8.706861590 -8.085785792 -0.881197479<BR>H 6.495783018 -9.318176822 0.885110623<BR>H 6.496059032 -9.320576794 -0.885035711<BR>H 7.728684857 -10.239295533 0.001631764<BR>H 1.135198120 2.173715970 -0.001216108<BR>H 5.082477354 -2.450566116 -0.000504931<BR></positions><BR></atoms></DIV>
<DIV><kmesh><BR><mesh nk="2 1 1" shift="0.0 0.0 0.0"/><BR></kmesh></DIV>
<DIV></input></DIV></DIV><br><!-- footer --><br><span title="neteasefooter"/><hr/>
<a href="http://email.163.com">网易邮箱,中国第一大电子邮件服务商</a>
</span>