<div>Dear  CP2K users and developers,<br></div><div><br></div><div>I am trying to optimal tow HSi(SiH3)3 molecules, but the result is unbelievable,as can be seen figure 1. Then I only optimal one HSi(SiH3)3 molecule using same input, and the result is good as shown figure 2. I have no idea about this problem, and hope somebody help me to solve the problem. Any help would be very appreciated!</div><div><img alt="figure2.png" data-iml="1124354.2000000477" width="273px" height="104.803px" src="cid:d7c31560-f1eb-4026-8f1b-7c3d726db4c7"><br></div><div>                        figure 1</div><div><img alt="figure1.png" data-iml="1140620.2000000477" width="234px" height="142.854px" src="cid:00064aca-e0f0-4244-a5d0-44d577fdb172"><br></div><div>                     figure 2</div><div><br></div><div>This is my input file :</div><div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &SUBSYS</div><div>    &KIND Si</div><div>      BASIS_SET DZVP-GTH-PADE</div><div>      POTENTIAL GTH-PADE-q4</div><div>    &END KIND</div><div>    &KIND H</div><div>      BASIS_SET DZVP-GTH-PADE</div><div>      POTENTIAL GTH-PADE-q1</div><div>    &END KIND</div><div>    &CELL</div><div>      ABC 12.0000 12.0000 12.0000</div><div>      PERIODIC NONE</div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>      COORD_FILE_NAME 1_2.xyz</div><div>      COORD_FILE_FORMAT XYZ</div><div>    &END TOPOLOGY</div><div>  &END SUBSYS</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME  BASIS_SET</div><div>    POTENTIAL_FILE_NAME  GTH_POTENTIALS</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-12</div><div>    &END QS</div><div>    &POISSON</div><div>       PERIODIC NONE</div><div>       POISSON_SOLVER MT</div><div>     &MT</div><div>      ALPHA 7.0</div><div>      REL_CUTOFF 1.2</div><div>     &END MT</div><div>    &END POISSON</div><div>    &MGRID</div><div>      NGRIDS 4</div><div>      CUTOFF 300</div><div>      REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</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 16</div><div>       &END</div><div>      &END VDW_POTENTIAL</div><div>    &END XC</div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 1.0E-5</div><div>      MAX_SCF 500</div><div>      &OT  ON</div><div>        PRECONDITIONER FULL_ALL</div><div>        MINIMIZER CG</div><div>        LINESEARCH 3PNT</div><div>      &END OT</div><div>      &MIXING  T</div><div>        METHOD BROYDEN_MIXING</div><div>        ALPHA 0.4</div><div>        NBROYDEN 8</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>  &END DFT</div><div>&END FORCE_EVAL</div><div>&MOTION</div><div> &GEO_OPT</div><div>   TYPE MINIMIZATION</div><div>   MAX_DR 1.0E-03</div><div>   MAX_FORCE 1.0E-03</div><div>   RMS_DR 1.0E-03</div><div>   RMS_FORCE 1.0E-03</div><div>   MAX_ITER 200</div><div>   OPTIMIZER BFGS</div><div> &END GEO_OPT</div><div>&END MOTION</div></div><div><br></div><div><br></div><div>Regards,</div><div>Guoying</div>