<div dir="ltr">Hi Patrick!<div><br></div><div>Thank you so much for your very thorough response. This absolutely answers my question on how to do this. I was trying to avoid having a "globally" charged system i.e. having an electric field, and I believe this does just that. I'm comparing this data directly with a simulation of block of Al and C and the output radial distribution function (RDF), with the aim being to see if there's a bond forming. I was specifically trying to compare the bond length of a simple Al-C system (neutral or charged ions) and see if that bond length that results corresponds to the RDF data. </div><div><br></div><div>Thank you for this! Have a great weekend!</div><div><br></div><div>Best Regards,</div><div>Devyesh Rana</div><div><br>On Saturday, February 22, 2020 at 6:22:06 AM UTC-5, Patrick Gono wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Devyesh,<div><br></div><div>The CHARGE keyword manages the total charge of the simulated system. For example, CHARGE -1 increases the number of electrons by one. Unfortunately, you can't really know where the electron decides to localize, if at all. You can modify your starting geometry to make the structure more conducive to the behavior of the extra charge you want to see. For instance, you could increase or decrease the bond length in your charged simulations.</div><div><br></div><div>You can also modify the core charge of selected ions through the CP2K_INPUT / FORCE_EVAL / SUBSYS / KIND / CORE_CORRECTION keyword. For example, you could try to put extra fractional charge on just the two ions forming the bond. You can then optimize the structure, remove the CORE_CORRECTION and add or subtract an electron by setting CHARGE appropriately. This way, the extra charge should, in principle, move to the bond you wish to study.</div><div><br></div><div>Finally, you can also take a look at constrained DFT, which allows you to enforce constraints on the electronic density: <a href="https://www.cp2k.org/howto:cdft" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Acdft\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHcqenU2A7g8h-hpnRTo2Q8aGqQqQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Acdft\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHcqenU2A7g8h-hpnRTo2Q8aGqQqQ';return true;">https://www.cp2k.org/<wbr>howto:cdft</a></div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono</div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, 21 Feb 2020 at 23:45, Dev Rana <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="TGyZBYW4AwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">d...@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">I was looking to put a charge to a C-Al bond. Specifically if there was a C+ or a C- to the molecule. <div><br></div><div>C+ -- Al</div><div>C -- Al</div><div>C- -- Al</div><div><br></div><div>Does anyone know how I would go about this?</div><div><br></div><div>Is it as simple as using Force_Eval > DFT > Charge 0, -1, +1?</div><div><br></div><div>Best Regards,</div><div>Devyesh Rana</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="javascript:" target="_blank" gdf-obfuscated-mailto="TGyZBYW4AwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/9fc83862-5daa-499a-9ea0-7b97228152f1%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/9fc83862-5daa-499a-9ea0-7b97228152f1%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/9fc83862-5daa-499a-9ea0-7b97228152f1%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/9fc83862-5daa-499a-<wbr>9ea0-7b97228152f1%<wbr>40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div></div>