<div dir="ltr">Hi<div><br></div><div>I have attached my input.</div><div><br></div><div>SL<br><br>On Saturday, May 9, 2015 at 3:44:24 PM UTC+1, Mengjia He wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Thank you SL.<div><br></div><div>I've no idea that A combination of Cartesian coordinates and HEXAGONAL cell could be a mistake, because I always used it for MD by CP2K. I have some confusions for transforming the Cartesian coordinates input file to fractional coordinates input file.</div><div><br></div><div>Would you please attach the input file for me?</div><div><br></div><div>Setting the mass of carbon atoms to 2.016 is my mistake while I copied it from my another input file.<br><br>Cheers.</div><div>Mengjia<br>On Saturday, May 9, 2015 at 3:31:30 AM UTC+8, S Ling wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I have tried to run your calculation with fractional coordinates instead of Cartesian coordinates, and it looks like the calculation is running fine. A combination of Cartesian coordinates and HEXAGONAL cell could be the reason of your job failure, and your initial error message "ERRORL2 in cp_fm_cholesky:cp_fm_cholesky_<wbr>decompose processor" might indicate that you have very bad geometry at that GEO_OPT step.<div><br></div><div>In addition, can you explain why you set the mass of carbon atoms to 2.016 in your input?</div><div><br></div><div>SL</div><div><br><br>On Friday, May 8, 2015 at 2:36:09 PM UTC+1, Mengjia He wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Thank you Matt.</div><div><br></div><div>YES, the input started running OK, but stopped at 17 steps with the error mentioned above.</div><div><br></div><div>I have carried out the modified inp file removing the STRESS_TENSOR NUMERICAL and increasing the cutoff to 500 Ry. The result seems to be more far away from convergence, though it isn't ending now. </div><div>( I've add the xyz to attanchment.)</div><div><br></div><div>Actually I wanna a CELL_OPT for calcite calculation at first. But the result turned out badly. The Convergence check always showed NO of all and the structures of xyz file became stranger , likely the whole crystal structure dispersing ( almost like the attanchment). So I turned  GEO_OPT to check first. </div><div><br></div><div>Could you please suggest me what else can be done for the GEO_OPT?<br></div><div><br></div><div>PS: As a new user of CP2K, would you like to tell me how to set the cutoff correctly and fast? (I learned the tutorial <a href="http://www.cp2k.org/howto:converging_cutoff" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\46sa\75D\46sntz\0751\46usg\75AFQjCNFnuxmJU7f0gdukVHPM7DfJudfcAA';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\46sa\75D\46sntz\0751\46usg\75AFQjCNFnuxmJU7f0gdukVHPM7DfJudfcAA';return true;">http://www.cp2k.org/howto:<wbr>converging_cutoff</a>. And I wonder if there's any file to summarize the cutoff set?)</div><div><br></div><div>Cheers.<br></div><div>Mengjia</div><br>在 2015年5月7日星期四 UTC+8下午9:16:43，Matt W写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>your input started running OK for me. Couple of things that might help:</div><div><br></div><div>Remove the STRESS_TENSOR NUMERICAL - you never want that, maybe STRESS_TENSOR ANALYTICAL if doing a cell optimization.</div><div> </div><div>Your cutoff is much too low for the BLYP functional and the elements in the system - this might be the cause of the problem later on. Increase the cutoff to at least 500 Ry, probably more. Check that the "Electronic density on regular grids:       -767.9999808982        0.0000191018" is correct to at least 6 or 7 figures as a rough idea.</div><div><br></div><div>Matt </div></div></blockquote></div></div></div></div></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div>