<div dir="ltr"><div>Hello everyone, </div><div><br></div><div>I am doing geometry optimization using cp2k for isopropanol alcohol (gas phase).  However, I don't know how to fix the error below: Can anyone help me with it? (I centered my molecule using         &CENTER_COORDINATES)</div><div>*** WARNING in pw/ps_wavelet_methods.F:236 :: Density non-zero on the edges of the unit cell: wrong results in WAVELET solver </div><div><br></div><div><div>Will the cell size affect the geometry optimization results for single-molecule in the gas phase?  and do I have to use a cubic cell?</div><div>       &CELL                         </div><div>       PERIODIC NONE                </div><div>       ABC 14.83 14.83 14.83 </div><div>       ALPHA_BETA_GAMMA 60.0 60.0 60.0</div><div>      &END CELL</div></div><div><br></div><div>Here is my input file: (my output file and xyz file are attached) </div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT IPA</div><div>  RUN_TYPE GEO_OPT</div><div>  PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL                           </div><div> </div><div>   METHOD Quickstep                   </div><div>   &DFT</div><div>    BASIS_SET_FILE_NAME /ihome/crc/build/cp2k/cp2k-6.1-intel-mkl/data/BASIS_MOLOPT </div><div>    &POISSON                      </div><div>      PERIODIC NONE</div><div>      PSOLVER  WAVELET             </div><div>    &END POISSON</div><div>    &QS                            </div><div>      METHOD GAPW              </div><div>    &END QS</div><div>      </div><div>      </div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 1.0E-05</div><div>      MAX_SCF 200</div><div>      &DIAGONALIZATION T</div><div>        ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &MIXING T</div><div>        ALPHA 0.5</div><div>        METHOD PULAY_MIXING</div><div>        NPULAY 5</div><div>      &END MIXING</div><div>      &PRINT</div><div>        &RESTART OFF</div><div>        &END RESTART</div><div>      &END PRINT</div><div>    &END SCF </div><div>    </div><div>    </div><div>    &XC</div><div>      &XC_FUNCTIONAL</div><div>        &PBE</div><div>        &END</div><div>      &END XC_FUNCTIONAL</div><div>      ! add Grimme D3 correction (without C9 terms, default)</div><div>      &VDW_POTENTIAL</div><div>        POTENTIAL_TYPE PAIR_POTENTIAL</div><div>        &PAIR_POTENTIAL</div><div>          PARAMETER_FILE_NAME dftd3.dat</div><div>          TYPE DFTD3</div><div>          REFERENCE_FUNCTIONAL PBE</div><div>          R_CUTOFF [angstrom] 16</div><div>        &END</div><div>      &END VDW_POTENTIAL</div><div>    &END XC</div><div> </div><div>   &END DFT</div><div> </div><div> </div><div>   &SUBSYS                            </div><div>     &CELL                         </div><div>       PERIODIC NONE                </div><div>       ABC 14.83 14.83 14.83 </div><div>       ALPHA_BETA_GAMMA 60.0 60.0 60.0</div><div>      &END CELL</div><div><br></div><div>      &TOPOLOGY</div><div>        &CENTER_COORDINATES</div><div>        &END CENTER_COORDINATES</div><div>        COORD_FILE_NAME IPA.xyz</div><div>        COORD_FILE_FORMAT XYZ</div><div>      &END</div><div><br></div><div>      &KIND C</div><div>        BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>        POTENTIAL GTH-PBE-q4</div><div>      &END KIND</div><div>      &KIND H</div><div>        BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>        POTENTIAL GTH-PBE-q1</div><div>      &END KIND</div><div>      &KIND O</div><div>        BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>        POTENTIAL GTH-PBE-q6</div><div>      &END KIND</div><div>      &KIND Zr</div><div>        BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>        POTENTIAL GTH-PBE-q12</div><div>      &END KIND</div><div>    &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div>&MOTION</div><div>  &GEO_OPT</div><div>    OPTIMIZER BFGS</div><div>    MAX_ITER 1500</div><div>    MAX_DR  [bohr] 0.003</div><div>    &BFGS</div><div>    &END</div><div>  &END</div><div>  &PRINT</div><div>    &RESTART</div><div>      BACKUP_COPIES 0</div><div>    &END RESTART</div><div>  &END PRINT</div><div>&END</div></div><div><br></div><div>Thanks!!!</div><div>Xin</div><div><br></div></div>