<div dir="ltr">Hi all,<div><br></div><div>I'm using the code to print the stresses on a system that I am running. Currently, the system has the input file I show below. For some reason, if I take off the evaluation of the stress tensor, it proceeds fine and carries out one optimization step after the initial SCF Wavefunction optimization. However, when I desire to evaluate the stress tensor, it just keeps doing SCF Wavefunction optimizations over and over. Could someone shed some insight into this?</div><div><br></div><div>Could it have to do with me using a numerical form of the stress tensor?</div><div><br></div><div>Thanks,</div><div>Rafael</div><div>=================</div><div>INPUT CODE</div><div><br></div><div><div>&FORCE_EVAL</div><div>STRESS_TENSOR NUMERICAL</div><div> &DFT</div><div>  &QS</div><div>   METHOD DFTB</div><div>   &DFTB</div><div>    SELF_CONSISTENT T</div><div>    DO_EWALD <span class="Apple-tab-span" style="white-space:pre">   </span>    F</div><div>    DISPERSION      F</div><div>    &PARAMETER</div><div>     PARAM_FILE_PATH /home/ras536/dftb_parameter_sets/cp2k/scc/</div><div>     PARAM_FILE_NAME parameter_table</div><div>     UFF_FORCE_FIELD uff_table</div><div>    &END PARAMETER</div><div>   &END DFTB</div><div>  &END QS</div><div>  &SCF</div><div>   SCF_GUESS ATOMIC</div><div>   MAX_SCF   200</div><div>   &OT ON</div><div>    LINESEARCH <span class="Apple-tab-span" style="white-space:pre">           </span>3PNT</div><div>    PRECONDITIONER <span class="Apple-tab-span" style="white-space:pre">    </span>FULL_SINGLE</div><div>   &END OT</div><div>  &END SCF</div><div> &END DFT</div><div> &SUBSYS</div><div>  &CELL</div><div># Change with the system coordinates</div><div>   ABC      21.50 22.14 100 </div><div>   PERIODIC XYZ</div><div># End of changes</div><div>  &END CELL</div><div>  &COORD</div><div># Paste here system coordinates</div><div># End of system coordinates</div><div>  &END COORD</div><div> &END SUBSYS</div><div> &PRINT</div><div>  &STRESS_TENSOR HIGH</div><div>   FILENAME system_stress.log </div><div>    &EACH</div><div>     GEO_OPT 1</div><div>    &END EACH</div><div>  &END STRESS_TENSOR          </div><div> &END PRINT</div><div>&END FORCE_EVAL</div><div><br></div><div><br></div><div>&GLOBAL</div><div>  PROJECT     geo_opt_syst</div><div>  RUN_TYPE    GEO_OPT</div><div>  PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div><br></div><div><br></div><div>&MOTION</div><div> &GEO_OPT</div><div>  MAX_ITER 1000</div><div> &END GEO_OPT</div><div># Uncomment to allow fixed atoms</div><div> &CONSTRAINT</div><div>  &FIXED_ATOMS</div><div>   COMPONENTS_TO_FIX <span class="Apple-tab-span" style="white-space:pre">    </span>X</div><div>   LIST <span class="Apple-tab-span" style="white-space:pre">          </span>91 100 101 110 111 120 121 130 131 140 141 150 151 160 161 170 171 180</div><div>  &END FIXED_ATOMS</div><div>  &FIXED_ATOMS</div><div>   COMPONENTS_TO_FIX <span class="Apple-tab-span" style="white-space:pre">        </span>Y</div><div>   LIST<span class="Apple-tab-span" style="white-space:pre">                   </span>1 2 3 4 5 6 7 8 9 10 171 172 173 174 175 176 177 178 179 180</div><div>  &END FIXED_ATOMS</div><div> &END CONSTRAINT</div><div># End of uncomment</div><div>&END MOTION</div></div></div>