<div dir="ltr">Dear Matthias,<br><br>By using LBFGS, the simulation stops in 27 steps and the result is better than I got after 1000 steps. But I still can not be sure whether it is finally optimized. However, thank you very much! I think it works very well.<br><br>Best regards,<br><br>在 2014年3月3日星期一UTC+1下午12时14分13秒，Matthias Krack写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Enlong,<br><br>either you provide a full setup (i.e. including the atomic coordinates) for the actual case, or you provide the input and the corresponding output file (at least part of it) showing the problem.<br>Did you try LBFGS?<br><br>Best,<br><br>Matthias<br><br>On Monday, March 3, 2014 11:57:27 AM UTC+1, Enlong Liu wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Matthias,<br><br>Thanks!<br><br>But I find out that the input I upload is a wrong one. So maybe you can have look at the right one and give me some suggestion. <br><br>Best regards,<br>Enlong<br><br>在 2014年3月3日星期一UTC+1上午11时12分52秒，<wbr>Matthias Krack写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Enlong,,<br><br>I would suggest to use BFGS (or LBFGS for large systems) instead of CG. Moreover, you may change the PRINT_LEVEL to medium to get more info. However, even if the run will converge somehow, the results won't be reasonable given the cell size of only 2.5A in z direction in a gamma point calculation.<br><br>Matthias<br><br>On Monday, March 3, 2014 10:41:11 AM UTC+1, Enlong Liu wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Fellows,<br><br>I am practicing CP2K by doing cell optimization and geometry optimization to CNT. In the cell_opt, I choose 'DIRECT_CELL_OPT<i>: Performs a geometry and cell optimization at the same time. The stress tensor is computed at every step</i>
                        '. And the max SCF iteration is 1000. However, every time I run the program, the program always used all 1000 iteration to get the result, which in my opinion, means the result does not converge. It seems like for every step, the threshold for energy is satisfied but between steps the threshold for geometry can not be satisfied. Could you help me to solve this issue?<br><br>PS: the input file is in the attachment.<br><br>Thanks a lot for your help!<br><br>Best regards,<br>Enlong<br></div></blockquote></div></blockquote></div></blockquote></div></blockquote></div>