<div dir="ltr"><div dir="ltr">Hi all, <div><br></div><div>Just to clarify, what I refer is that the box updating the progress of the geometry optimization:</div><div><br></div><div><div> --------  Informations at step =     1 ------------                            </div><div>  Optimization Method        =                 BFGS</div><div>  Total Energy               =      -495.6366242369</div><div>  Real energy change         =        -2.2642515093</div><div>  Predicted change in energy =        -2.7904990758</div><div>  Scaling factor             =         0.5267466124</div><div>  Step size                  =         0.4724315332</div><div>  Trust radius               =         0.4724315332</div><div>  Decrease in energy         =                  YES</div><div>  Used time                  =               95.726</div><div><br></div><div>  Convergence check :</div><div>  Max. step size             =         0.4724315332</div><div>  Conv. limit for step size  =         0.0030000000</div><div>  Convergence in step size   =                   NO</div><div>  RMS step size              =         0.1432565779</div><div>  Conv. limit for RMS step   =         0.0015000000</div><div>  Convergence in RMS step    =                   NO</div><div>  Max. gradient              =         0.0773821875</div><div>  Conv. limit for gradients  =         0.0004500000</div><div>  Conv. for gradients        =                   NO</div><div>  RMS gradient               =         0.0232460107</div><div>  Conv. limit for RMS grad.  =         0.0003000000</div><div>  Conv. for gradients        =                   NO</div><div> -----------------------------<wbr>----------------------</div><div><br></div><div>Is gone, and the program just keeps looping over the SCF wavefunctions. It could possibly be that the optimization algorithm changed (since I haven't specified it). Also, I can't find the stresses printed in the log file (or in a separate file if I specify it).</div><div><br></div><div>I checked if the algorithm had changed (say, from BFGS to CG) but it didn't. It's still BFGS for all my trial runs.</div><div><br></div><div>Thanks,</div><div>Rafael</div><br>On Thursday, May 1, 2014 11:54:13 AM UTC-5, Rafael Soler-Crespo wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;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_<wbr>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></div></div>