Thank you so much, Frederick!<br /><br /><div>That helped a lot.</div><div><br /></div><div>Michela</div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, September 14, 2024 at 12:53:49 PM UTC-4 Frederick Stein 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;"><div>Dear Michela,</div><div>Given the basis set file BASIS_MOLOPT in your input file, check the q# of the basis set in use. Often, they have two names, a short one and a long one. The latter contains the number of electrons, here 13.</div><div>In case the number of electrons is not provided, check the basis set itself. The manual gives some hints (compare <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND/BASIS.html#CP2K_INPUT.FORCE_EVAL.SUBSYS.KIND.BASIS.DEFAULT_KEYWORD" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND/BASIS.html%23CP2K_INPUT.FORCE_EVAL.SUBSYS.KIND.BASIS.DEFAULT_KEYWORD&source=gmail&ust=1726580005366000&usg=AOvVaw0_v3b3_LJ81PiMJMbWuYpp">https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND/BASIS.html#CP2K_INPUT.FORCE_EVAL.SUBSYS.KIND.BASIS.DEFAULT_KEYWORD</a>). The essence is that each block of functions (in your case one) starts with several numbers in the format</div>
<div><span>n lmin lmax nexp nshell(lmin) nshell(lmin+1) ... nshell(lmax-1) nshell(lmax)</span></div><div><span>The first number is ignored and kept for reasons of compatability. With the name of the basis set, the number of orbitals and the periodic table, you can determine the number of electrons yourself.<br></span></div><div><span>Best,</span></div><div><span>Frederick<br></span></div>

<div class="gmail_quote"><div dir="auto" class="gmail_attr">Michela Benazzi schrieb am Samstag, 14. September 2024 um 18:38:33 UTC+2:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello everyone,<br><br>I am running MD on a system of 500 Ga atoms in a 21x21x21 box (liquid density) at 1 atm and 30 C. I am using GPAW, PBE, single zeta and q3 BLYP potentials.<br><br>I just got this error: "Basis-set and pseudo-potential of atomic kind 'Ga' were optimized for different valence electron numbers.".<br><br><div>I found a thread of someone who had the same issue with a Pt simulation, so I think I should try to increase to q13 BLYP potential, also available in CP2K. However, I want to make a more informed decision: is q13 enough? What is q# for potentials? I tried looking for literature and methods online but could not find an explanation.</div><div><br></div><div>I would love some advice. Thank you kindly!<br><br>Michela<br></div></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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/a746a78c-b20a-4a95-b620-3a51c2a839ban%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/a746a78c-b20a-4a95-b620-3a51c2a839ban%40googlegroups.com</a>.<br />