<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: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;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:0;margin-left:0.8ex;border-left:1px #ccc solid;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>