<div dir="ltr"><div>Dear Chen,</div><div><br></div><div>A CELL_OPT calculation automatically finds the optimal atomic configuration as well. In general, the DFT lattice parameter will not be identical to the experimental one. Therefore, it makes sense to run a CELL_OPT calculation to determine the correct unit cell parameters and the atomic structure. In practice, such calculations do not take much more time than a GEO_OPT simulation. Note that in the case of bulk Pt, you know that the cell is cubic, and thus can constrain the angles during CELL_OPT by switching the option KEEP_ANGLES in the &CELL_OPT subsection of the &MOTION section to true.</div><div><br></div><div>But regarding your original issue: CP2K, by default, only uses a single k-point, the Gamma point. Whereas in a plane-wave code you would converge any results with respect to the k-point mesh density, in CP2K you must converge the results with respect to the supercell size. I would suggest running a CELL_OPT calculation with at least a 2x2x2, preferably a 3x3x3 supercell, and compare the lattice parameters with literature (still, the number will be off due to the functional used, so compare with PBE results).</div><div><br></div><div>Note that you must specify the MULTIPLE_UNIT_CELL keyword in both the &TOPOLOGY and the &CELL subsections!</div><div><br></div><div>Of course, get rid of that extra charge, you most probably want to study a neutral system. And if you are dealing with bulk systems, keeping fixed the position of an atom is not necessary, as the system will not float around. This might only be relevant in the case of e.g. molecules or 2D sheets in vacuum.</div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 28 Aug 2019 at 12:54, Chn <<a href="mailto:chen...@gmail.com">chen...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Daniele,<div><br></div><div>Thanks for your answer! I annotated the keyword "CHARGE" and the constraint part as you mentioned, but it still doesn't work. Could you please tell me some reasons that why does this situation happen in general? I am wandering that whether it is caused by the input file.</div><div>Does the variation of cell volume you mentioned above mean that I need to do the CELL_OPT before GEO_OPT?</div><div><br></div><div>thanks again,</div><div>chn<br><br>在 2019年8月28日星期三 UTC+8下午5:26:58，Daniele Passerone写道：<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Chn<div><br></div><div>Apart from the size of the cell that could give you some trouble, the problem is the keyword</div><div><br></div><div>CHARGE 1</div><div><br></div><div>that you added in the &DFT section, meaning that you are simulating a charged system. Removing that would work.</div><div><br></div><div>Indeed in the output we see:</div><div><br></div><div>
<p><span><span> </span>Hirshfeld Charges</span></p>
<p><span></span><br></p>
<p><span><span> </span>#Atom<span> </span>Element<span> </span>Kind<span> </span>Ref Charge <span> </span>Population<span> </span>Net charge</span></p>
<p><span><span> </span>1 <span> </span>Pt <span> </span>1<span> </span>18.000 <span> </span>17.750 <span> </span>0.250</span></p>
<p><span><span> </span>2 <span> </span>Pt <span> </span>1<span> </span>18.000 <span> </span>17.750 <span> </span>0.250</span></p>
<p><span><span> </span>3 <span> </span>Pt <span> </span>1<span> </span>18.000 <span> </span>17.750 <span> </span>0.250</span></p>
<p><span><span> </span>4 <span> </span>Pt <span> </span>1<span> </span>18.000 <span> </span>17.750 <span> </span>0.250</span></p>
<p><span></span><br></p>
<p><span><span> </span>Total Charge<span> </span>1.000</span></p><p><span><br></span></p><p><span><br></span></p><p><span><br></span></p></div><div><br></div><div><br></div><div>Also not necessary is the constraint on the position of one atom that you have in the optimization section.</div><div>In general you should also vary the cell volume and find what is the equilibrium volume for the PBE functional.</div><div>Hope this helps</div><div><br></div><div>Daniele</div><div><br></div><div><br>On Wednesday, August 28, 2019 at 10:43:05 AM UTC+2, Chn wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear experts in CP2K,<br><br>I am a learner of CP2K package and I am really thankful if my naive question can be explained.<div>I tried to do geometry optimization toward Pt bulk or Pt(111) by set the keyword "RUN_TYPE" as GEO_OPT and the *.inp file was referred to the exercises "Adsorption of acetylene on an intermetallic surface" about PdGa at CP2K website. But when the calculation done after several hours, I found that the atoms in cell distorted seriously. I think it should be obvious incorrect but I have no idea about how to fix it. Many thanks in advance.<div>The imported .xyz was prepared by using materials studio and VESTA. I have tried to use larger supercell (3*3*3) or to set the keyword "MULTIPLE_UNIT_CELL" as 3 3 3 and tried to optimize other structures like TiO2 but no use..is it that the cell I used still too small?</div><div>The attachments were the input files I used.<br></div><div><br></div><div>Regards,</div><div>chn</div><p style="text-align:center;clear:both"><a style="clear:left;margin-bottom:1em;float:left;margin-right:1em"><img src="https://groups.google.com/group/cp2k/attach/3b8e602657d89/o.png?part=0.5&view=1&authuser=0" alt="o.png" width="200" height="174" border="0"></a><a style="margin-left:1em;margin-right:1em"><img src="https://groups.google.com/group/cp2k/attach/3b8e602657d89/d.png?part=0.6&view=1&authuser=0" alt="d.png" width="200" height="198" border="0"></a></p><div><br></div></div></div></blockquote></div></div></blockquote></div></div>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/7092cb9a-be30-4f63-8246-b988dea87c67%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/7092cb9a-be30-4f63-8246-b988dea87c67%40googlegroups.com</a>.<br>
</blockquote></div>