Numerical stress ternsor requires many single point energy calculations to estimate the stress tensor using finite differences. Should pretty much only be used for debugging.<br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, March 29, 2021 at 7:20:17 AM UTC+1 Wei Chen wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hi Rafael,</div><div>I got the same problem, but when I changed "NUMERICAL" to "<span>ANALYTICAL</span>", everything is fine. While, <span><span><span>I don't know the reason.</span></span></span> <br></div><div><br></div><div>Wei<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, May 2, 2014 at 12:54:13 AM UTC+8 <a href data-email-masked rel="nofollow">ras...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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 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 style="white-space:pre">                </span>3PNT</div><div>    PRECONDITIONER <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 style="white-space:pre"> </span>X</div><div>   LIST <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 style="white-space:pre"> </span>Y</div><div>   LIST<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></blockquote></div></blockquote></div>