Unfortunately, I was never able to figure what was wrong with the inputs. I ended up using a pure DFT functional for the whole system which was just under 1000 atoms. The simulation was costly but gave me very good data. Depending on the size of your system you may consider this as well. <br /><div><br /></div><div>I also wish I still had access to these files so I could send them to you but here is an explanation of what I remember: Each subsystem file contained coordinates for the corresponding layer of the calculation. For example, "dft_subsys.inc" had coordinates for the solute I was treating with DFT (along with basis sets, pseudopotentials, and cell section). I believe the only difference between this file and "solute_subsys.inc" was that I removed anything pertaining to DFT. The file "combined_subsys.inc" contained coordinates for my solute and all the solvent molecules which I was treating with PM6. The idea was to manually implement ONIOM by adding the DFT energy and PM6 energy of the whole system and then removing the PM6 energy of the solute. I attempted to define this in the main input under the first &FORCE_EVAL section.</div><div><br /></div><div>I'm still not sure if this would have worked if I had resolved the error, but I hope it serves as a decent starting point for you. Feel free to respond with any questions.<br /><br /> -Josh</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, October 18, 2023 at 4:41:24 PM UTC+2 Vittorio Bariosco 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;">Dear Josh,<div><br></div><div>Did you manage to setup the ONIOM input? I would like to set an oniom too, may I ask you what do you defined inside the files: <font color="#000000" style="background-color:white"><span style="font-family:Menlo">dft_subsys.inc,</span><font face="Menlo">combined_subsys.inc <span>and</span> </font><span style="font-family:Menlo">solute_subsys.inc.</span></font></div><div><font color="#000000" style="background-color:white"><font face="Menlo"><br></font></font></div><div><font color="#000000" style="background-color:white"><font face="Menlo">Thank you in advance!</font></font></div><div><font color="#000000" style="background-color:white"><font face="Menlo"><br></font></font></div><div><font color="#000000" style="background-color:white"><font face="Menlo">Best,</font></font></div><div><font color="#000000" style="background-color:white"><font face="Menlo"><br></font></font></div><div><font color="#000000" style="background-color:white"><font face="Menlo">Vittorio<br></font></font><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno giovedì 30 marzo 2023 alle 12:43:45 UTC+2 Marcella Iannuzzi ha scritto:<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 ..<div><br></div><div>The input seems OK<br>From the error message it is possible that there is an inconsistency between the number of atoms effectively found in subsystem 1 </div><div>with respect to what defined as FRAGMENT </div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, March 23, 2023 at 3:39:18 AM UTC+1 Josh W 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>Hello,</div><div><br></div><div>I am trying to perform an ONIOM type calculation. I am getting errors related to how I am mapping my fragments and I'm not sure what's wrong.</div><div><br></div><div>I have two questions:</div><div><br></div><div>1) How do I solve the issue with mapping my fragments?</div><div><br></div><div>2) Is this a reasonable setup for an ONIOM calculation in CP2K?</div><div><br></div><div>Thanks!<img alt=""><img alt=""><img alt=""><img alt=""><img alt=""><img alt=""><img alt=""></div></blockquote></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/40926c47-afe9-4315-8bd7-0daefe788314n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/40926c47-afe9-4315-8bd7-0daefe788314n%40googlegroups.com</a>.<br />