Dear Marcella,<div>Thanks for your response, I have one more question about calculating homo-lumo gap<br></div><div>for the attached results, I calculated this way: -0.067428+0.075159=0.0077</div><div>Am I calculating it the right way?</div><div><br></div><div>Best regards,</div><div><img alt="cp2k g.PNG" data-iml="439513.0999999965" width="280px" height="179px" src="cid:d0ce1ede-429a-4283-9293-04b8d92c6d27"><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">في الجمعة، 12 فبراير 2021 في تمام الساعة 4:14:35 م UTC+3، كتب Marcella Iannuzzi رسالة نصها:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-right: 1px solid rgb(204, 204, 204); padding-right: 1ex;">Dear Malin Khalil<div><br></div><div>An easy way is by calculating the PDOS</div><div>Regards</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, February 12, 2021 at 1:40:47 PM UTC+1 <a href data-email-masked rel="nofollow">mal...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear CP2K users<br><div><div>I am working on studying the homo lumo gap of carbon nanotubes, but the output file does not show it.</div><div>I have already perform calculations without smearing and the output was showing the homo lumo gap</div><div>could you please tell me where is the wrong?</div><div><br></div><div><br></div><div>best regards, </div><div>Malin Khalil</div></div><div><br></div><div>my input:</div><div><div>&GLOBAL</div><div>  PROJECT CNT</div><div>  RUN_TYPE ENERGY</div><div>  PRINT_LEVEL LOW</div><div><span style="white-space:pre">  </span>PREFERRED_FFT_LIBRARY FFTW</div><div>&END GLOBAL</div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &SUBSYS</div><div>    &CELL </div><div>ABC 34.000 34.000 34.000</div><div>PERIODIC NONE</div><div><br></div><div>    &END CELL</div><div>    &TOPOLOGY</div><div><span style="white-space:pre"> </span>&CENTER_COORDINATES</div><div><span style="white-space:pre">       </span>&END CENTER_COORDINATES</div><div>      COORD_FILE_NAME <a href="http://24.xyz" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=ar&q=http://24.xyz&source=gmail&ust=1613326046249000&usg=AFQjCNHyorPPf3GtWFizkCsdXC5dsUiH2A">24.xyz</a></div><div>      COORD_FILE_FORMAT XYZ</div><div>    &END</div><div> <span style="white-space:pre">         </span>&KIND C</div><div>      BASIS_SET DZVP-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div><span style="white-space:pre">     </span>  <span style="white-space:pre">  </span>&END KIND</div><div><span style="white-space:pre">         </span>&KIND H</div><div>      BASIS_SET DZVP-GTH</div><div>      POTENTIAL GTH-PBE-q1</div><div><span style="white-space:pre">     </span>  <span style="white-space:pre">  </span>&END KIND</div><div><br></div><div>  &END SUBSYS</div><div>    &DFT<span style="white-space:pre"> </span></div><div>   BASIS_SET_FILE_NAME  HFX_BASIS</div><div>   POTENTIAL_FILE_NAME  POTENTIAL</div><div>        </div><div><span style="white-space:pre">            </span>MULTIPLICITY 1</div><div><span style="white-space:pre">                </span>CHARGE 0</div><div><span style="white-space:pre">              </span>        </div><div><span style="white-space:pre"> </span>&MGRID</div><div>      NGRIDS 4</div><div>      CUTOFF 300</div><div>      REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &QS</div><div><br></div><div>       METHOD GPW </div><div>     EPS_DEFAULT 1.0E-14</div><div><span style="white-space:pre">     </span>EPS_FILTER_MATRIX 0</div><div>      EPS_PGF_ORB 1E-32</div><div>    &END</div><div>    &PRINT</div><div>       &E_DENSITY_CUBE ON</div><div>       &END E_DENSITY_CUBE</div><div><br></div><div>      &MO_CUBES</div><div>         NHOMO  5</div><div>         NLUMO  5</div><div>WRITE_CUBE .True.</div><div><br></div><div>      &END MO_CUBES</div><div>    &END</div><div><br></div><div>    &SCF</div><div>     SCF_GUESS ATOMIC</div><div>     EPS_SCF 8.0E-5</div><div>     MAX_SCF 100</div><div>     ADDED_MOS  50</div><div>     CHOLESKY INVERSE</div><div>     &SMEAR  ON</div><div>         METHOD FERMI_DIRAC</div><div>         ELECTRONIC_TEMPERATURE [K] 800</div><div>     &END SMEAR</div><div>     &DIAGONALIZATION</div><div>          ALGORITHM STANDARD</div><div>     &END DIAGONALIZATION</div><div>     &MIXING</div><div>          METHOD BROYDEN_MIXING</div><div>            ALPHA 0.2</div><div>            BETA 1.5</div><div>            NBROYDEN 8</div><div>      &END MIXING</div><div>    &END SCF</div><div>    &XC</div><div>    &XC_FUNCTIONAL</div><div>       &LYP</div><div>         SCALE_C 0.81</div><div>       &END</div><div>       &BECKE88</div><div>         SCALE_X 0.72</div><div>       &END</div><div>       &VWN</div><div>         FUNCTIONAL_TYPE VWN5</div><div>         SCALE_C 0.19</div><div>       &END</div><div>       &XALPHA</div><div>         SCALE_X 0.08</div><div>       &END</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div><br></div><div><span style="white-space:pre">    </span>&POISSON</div><div><span style="white-space:pre">  </span>PERIODIC NONE</div><div><span style="white-space:pre"> </span>POISSON_SOLVER WAVELET</div><div><span style="white-space:pre">        </span>&END POISSON</div><div>  &END DFT</div><div>  &END FORCE_EVAL</div></div><div><br></div></blockquote></div></blockquote></div>