<div dir="ltr">Dear Nikhil,<div><br></div><div>I have a couple things to suggest:</div><div>1) Your CUTOFF is definitely too small. Please, make sure your model is converged with respect to your CUTOFF by running a series of single-point (ENERGY) calculations, as described in: <a href="https://www.cp2k.org/howto:converging_cutoff">https://www.cp2k.org/howto:converging_cutoff</a>. In general, the values for CUTOFF should be in the range of (higher) triple digits, 400 - 800 Ry, or even more. Whenever possible, make sure you are converged.</div><div>2) You are not using any pseudopotential for your Pb and Cs atoms. Look into the GTH_POTENTIALS file and choose the appropriate one for your choice of basis sets.</div><div>3) Are you sure about your cell parameters? If your plan is to study a periodic system, you must specify periodic symmetry and choose appropriate cell parameters. If you plan to simulate an isolated nanoparticle, consider reducing the size of your cell. 50 angstrom in each direction  might be an overkill and needlessly slow down your jobs.</div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Sep 26, 2019 at 9:01 PM Nikhil Maroli <<a href="mailto:scin...@gmail.com">scin...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">The bonds are breaking during the optimization, Initially, i have tried to optimize the structure with gaussian but it started to break bonds even at correct distances. Anyone can give any suggestion to avoid bond breaking in CP2K and to achieve optimized structure.<div>Inputs are given below for your information.</div><div><div>&GLOBAL</div><div> PROJECT job2</div><div> RUN_TYPE GEO_OPT</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div><br></div><div> METHOD QS</div><div> &DFT</div><div># Please copy the GTH_BASIS_SETS and POTENTIAL files</div><div># which may have been in cp2k-x.x.x/data or cp2k-x.x.x/tests/QS directory</div><div># to the CP2K executable directory if you use these default files.</div><div>   CHARGE 0</div><div>   MULTIPLICITY 2</div><div>  BASIS_SET_FILE_NAME  GTH_BASIS_SETS</div><div>  POTENTIAL_FILE_NAME  POTENTIAL</div><div>  &MGRID</div><div>   NGRIDS 4</div><div>   CUTOFF 30</div><div>  &END MGRID</div><div>  &XC</div><div>   &XC_FUNCTIONAL PBE</div><div>   &END XC_FUNCTIONAL</div><div>  &END XC</div><div> &END DFT</div><div> &SUBSYS</div><div>  &COORD</div><div>    Pb         1.94602       -2.12632        1.25157</div><div>    Br         3.01219       -0.06256        2.35356</div><div>    Pb         2.21556        2.03217        1.09683</div><div>    Br        -0.12169        3.23840        1.92449</div><div>    Pb        -2.34714        2.39328        0.99984</div><div>    Br        -3.57990        0.20869        1.95980</div><div>    Pb        -2.65498       -2.14276        1.00547</div><div>    Br        -0.34094       -3.29325        1.98965</div><div>    Pb         2.06448       -1.33415       -1.21263</div><div>    Pb         2.13521        1.12177       -1.15878</div><div>    Br         0.06365        1.19694       -0.08612</div><div>    Pb        -1.79865        1.47428       -1.68114</div><div>    Br        -0.57583       -0.01056       -3.59712</div><div>    Pb        -1.98293       -1.29020       -1.66395</div><div>    Br         0.01062       -1.18331       -0.10788</div><div>    Br        -3.82309       -3.16020       -1.27765</div><div>    Br         3.19291       -3.45179       -0.50642</div><div>    Br        -3.41601        3.49721       -1.28741</div><div>    Cs        -0.28609        0.02294        3.08499</div><div>  &END COORD</div><div><br></div><div>   &KIND Br</div><div>    BASIS_SET DZVP-GTH</div><div>    POTENTIAL GTH-PBE-q7</div><div>   &END KIND</div><div><br></div><div>   &KIND Cs</div><div>    BASIS_SET DZVP-GTH</div><div>    POTENTIAL # type pseudopotential functional name !</div><div>   &END KIND</div><div><br></div><div>   &KIND Pb</div><div>    BASIS_SET DZVP-GTH</div><div>    POTENTIAL # type pseudopotential functional name !</div><div>   &END KIND</div><div><br></div><div>  &CELL</div><div>   ABC 50 50 50 # Provisional Cell Param. Please "Add Unit Cell" !</div><div>  &END CELL</div><div> &END SUBSYS</div><div>&END FORCE_EVAL</div></div><div><br></div></div>

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/96a565a1-0f85-4692-98fe-68f8d43c59e9%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/96a565a1-0f85-4692-98fe-68f8d43c59e9%40googlegroups.com</a>.<br>
</blockquote></div>