<div dir="ltr">Hi all,<div><br></div><div>When I use the example regtest-gapw-1/c8_broy_gapw_all.inp to print out the MO_CUBES. It works fine. But if I change the basis set from 6-31Gxx to 6-311Gxx it gave me the cube files with NaNs in it. Have you ever seen this problem or do you have an idea of where is the problem?</div><div><br></div><div>Yi</div><div><br></div><div>I use cp2k 3.0</div><div><br></div><div>the modified input</div><div>====</div><div><div>&GLOBAL</div><div>  PROJECT  c8_broy_gapw_all</div><div>  PRINT_LEVEL MEDIUM</div><div>  RUN_TYPE ENERGY</div><div>  FLUSH_SHOULD_FLUSH</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  &PRINT</div><div>    &FORCES ON</div><div>    &END</div><div>  &END</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME EMSL_BASIS_SETS</div><div>    POTENTIAL_FILE_NAME POTENTIAL</div><div>    LSD</div><div>    &MGRID</div><div>      NGRIDS 4</div><div>      CUTOFF 100</div><div>    &END MGRID</div><div>    &QS</div><div>      METHOD GAPW</div><div>      MAP_CONSISTENT</div><div>      EXTRAPOLATION PS</div><div>      EXTRAPOLATION_ORDER 2</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>      MAX_SCF   50</div><div>      EPS_SCF 5.0e-3</div><div>      &SMEAR</div><div>        METHOD FERMI_DIRAC</div><div>        ELECTRONIC_TEMPERATURE   500.</div><div>        FIXED_MAGNETIC_MOMENT  0.0</div><div>      &END</div><div>      &MIXING</div><div>          METHOD BROYDEN_MIXING</div><div>          ALPHA   0.1</div><div>          BETA  1.0</div><div>          NBUFFER 8</div><div>      &END</div><div>      ADDED_MOS   20</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div>    &PRINT</div><div>      #&MO</div><div>      #  &EACH</div><div>      #    QS_SCF 100</div><div>      #  &END</div><div>      #  ADD_LAST NUMERIC</div><div>      #  EIGENVALUES</div><div>      #  OCCUPATION_NUMBERS</div><div>      #&END</div><div>      &MO_CUBES</div><div>        NHOMO  64</div><div>        NLUMO 20</div><div>        WRITE_CUBE T</div><div>      &END</div><div>    &END</div><div>  &END DFT</div><div><br></div><div>  &SUBSYS</div><div>    &CELL</div><div>       ABC   5.42858871335 5.42858871335 5.42858871335</div><div>    &END CELL</div><div><br></div><div>    &KIND  C</div><div>      BASIS_SET  <span style="background-color: rgb(255, 255, 0);">6-311Gxx</span></div><div>      POTENTIAL  ALL</div><div>    &END</div><div><br></div><div>    &TOPOLOGY</div><div>      COORD_FILE_NAME ../sample_xyz/C_8.xyz</div><div>      COORDINATE XYZ</div><div>      CONNECTIVITY OFF</div><div>    &END TOPOLOGY</div><div><br></div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div></div><div>====</div></div>