<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>