Hi, there are at least a couple of problems in your input files. <div><br /></div><div>The first one is the one mentioned by marcella: The termination is a 100, that is going to reconstruct to increase the lateral density (to a triangular close packed layer). </div><div>The other one is that your cell has a lateral dimension <br /><br /> &CELL<br /> ABC 22.9836444855 40.0 22.9836444855 <br /> PERIODIC XZ<br /> &END CELL</div><div><br /></div><div>And since the default (see manual<br /><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/CELL.html">https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/CELL.html</a></div><div><br />) is Angstrom, your system has a lot of empty space laterally, and this is seen in the optimization trajectory where the not-fixed atoms are behaving like a cluster and not as a slab.</div><div><br /></div><div>HTH<br /><br />Daniele</div><div><br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, July 9, 2025 at 11:50:20 AM UTC+2 Marcella Iannuzzi wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi Hamid, <div><br></div><div>Are you sure that the Au coordinates correspond to a Au(111) slab?</div><div>I just had a quick look, but the staking looks weird.</div><div><br></div><div>Regards</div><div>Marcella<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, July 9, 2025 at 6:47:51 AM UTC+2 <a href data-email-masked rel="nofollow">hamidzabihi...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<p>Dear CP2K community,</p>
<p>I am performing a geometry optimization of a <strong>C₆₀ molecule adsorbed on an Au(111) slab</strong>, using <strong>PBE+D3</strong> and <strong>DZVP-MOLOPT-SR-GTH</strong> basis sets.</p>
<p>The initial configuration had C₆₀ placed ~3 Å above the gold surface. However, after optimization, I noticed that <strong>several top-layer Au atoms were pulled significantly upward</strong>—toward the C₆₀. This distortion seems nonphysical and may indicate issues with my setup or parameters.</p>
<p>I would be grateful for any advice on:</p>
<ul><li>
<p>Best practices for modeling <strong>molecule-surface interactions</strong> in CP2K</p>
</li><li>
<p>Ensuring <strong>slab stability</strong> and correct <strong>adsorption behavior</strong></p>
</li><li>
<p>Whether this is a known issue with this CP2K version or functional</p>
</li></ul>
<p><strong>System details:</strong></p>
<ul><li>
<p>CP2K version: 8.0 (development)</p>
</li><li>
<p>Functional: PBE with D3 dispersion correction</p>
</li><li>
<p>Basis sets: DZVP-MOLOPT-SR-GTH</p>
</li><li>
<p>Slab: 4-layer Au(111), ~15 Å vacuum in y-direction, bottom layers unconstrained</p>
</li></ul>
<p>I’ve attached the following files:</p>
<ul><li>
<p><span>cp2k_c60_au_geoopt.inp</span>: Input file</p>
</li><li>
<p><span>out-C60_Au_GEO_OPT.log</span>: Optimization output</p>
</li><li>
<p><span>C60_Au_opt-trajectory.xyz-pos-1</span>: Final trajectory</p>
</li></ul>
<p>Thank you in advance for your help!</p>
<p>Best regards,<br>
<strong>Hamid Zabihi</strong></p>
<br></blockquote></div></blockquote></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:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/361844b1-fa13-4309-8cc9-40fd6d48122bn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/361844b1-fa13-4309-8cc9-40fd6d48122bn%40googlegroups.com</a>.<br />