<div><br></div><div><br></div>Dear Matthew<div><br></div><div>I am not a purist of basis set and as long as relevant properties as the representation of the electronic density, the resulting density of states, binding energies are converged (within the level of theory of choice) I am not that concerned in taking specifically optimised basis sets. </div><div><br></div><div>1) The massive thermostat is a Local thermostat, i.e. it enforces the correct distribution of <span style="font-style: italic;">each component </span>of the momentum</div><div>This means that the velocities of the particles and thus their dynamics are modified, with respect to the microcanonical ensemble, which could potentially lead to thermostat-dependent dynamical artifacts. In general NVE is better for IR spectra</div><div>https://pubs.acs.org/doi/10.1021/ct400109a</div><div><br></div><div>Best</div><div>Marcella<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, May 13, 2022 at 6:44:06 AM UTC+2 mhvg...@gmail.com 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>Hi Marcella</div><div><br></div>Thank you so much, I can't believe I missed the UZH Ar basis set! To be honest, I didn't consider it because I thought that UZH was some theoretical method I'd never heard of (like HFX or ADMM), and I only just found out it stood for University of Zurich...<div><br></div><div>Unfortunately, there's no HFX basis set for argon (even on the git repository), so I can't use that for my system. Would I be right in assuming that, if I use BASIS_ADMM-UZH basis sets, I should also use POTENTIAL_UZH and the BASIS_MOLOPT_UZH? Originally, I was going to use the DZVP-GTH basis sets because (I believe) they are smaller. I also don't know whether the fact that the UZH sets were optimised against the PBE0 will cause issues with me using B3LYP as my functional. What do you think?<br><div><br></div><div>1) In that case, why would someone equilibrate their system with a MASSIVE thermostat first, and then again with a GLOBAL thermostat? This is what is done in the IR spec tutorial...</div></div><div>2) I thought that might be the case. Thank you!</div><div>3) Of course, that would make life too easy ;)</div><div><br></div><div>Thank you so much for your help! I have no experience with any type of bulk calculation, so your advice is very much appreciated</div><div><br></div><div>Regards,</div><div><br></div><div>Matthew</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, May 13, 2022 at 1:37:18 AM UTC+8 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">Dear Matthew<div><br></div><div>There are Ar ADMM basis sets in BASIS_ADMM_UZH in cp2k/data (at least present git repository)<br>The ADMM basis set can be tested against full HFX calculations.</div><div><br></div><div>1) Using MASSIVE there is one  thermostat attached to each degree of freedom, with GLOBAL only one thermostat controls all degrees of freedom</div><div>2) Larger cutoff are needed to collocate on the real space grid functions with larger exponents. In this respect, yes there is a basis set dependency</div><div>3) No rule to thumb, it depends what the simulation goals are</div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, May 11, 2022 at 9:15:36 AM UTC+2 <a rel="nofollow">mhvg...@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">Hi there,<div><br></div><div>I was wondering if I could get some advice/help on running an AIMD calculation to run? I've <i>no</i> experience with AIMD, MD, or any plane wave calculations (my only experience is with ab initio methods in ORCA and Gaussian), and no one else in my department has any experience with AIMD calculations, either... I've finally gotten to the stage where I can start an AIMD calculation (it's taken weeks to get to this stage!), but I have come across a persistent issue which I do not have the expertise to solve.</div><div><br></div><div>Essentially, I want to simulate my molecular system in a solid argon matrix, and calculate IR spectra to compare with experiment. So far, I've optimised my cutoff/rel_cutoff values, gotten my input files sorted, gotten CP2K installed installed on our local HPC, and tested that it works for a few of the test systems provided for the IR/RAMAN and HFX tutorials.</div><div><br></div><div>Based on results from a recent study, I want to use B3LYP with a double zeta basis set, which I am very well aware will take a long time to complete... According to the CP2K website (and my own testing), it looks like using either the ADMM or RI approximations might speed up this calculation to a rate which will be acceptable for our HPC. The problem is that I cannot find any <i>matching</i> primary, potential and auxiliary basis sets for argon. In fact, I cannot find any auxiliary basis sets for argon at all.</div><div><br></div><div>Is this because argon doesn't require an aux basis set for RI or ADMM? Does anyone know of any basis set families which I could use for my calculations? They'd need to be of at least double zeta quality, and include basis sets for H, C, N, S and Ar.</div><div><br></div><div>Also, three unrelated questions (but things which I've been confused/unsure about when trying to learn all this stuff):</div><div><ol><li>what's the difference between the 'global' and 'massive' thermostats? The academic (yes, there's only one) at my university who has experience with MD had never seen those terms before, so I'm not sure how they differ.<br></li><li>Are the cutoff/rel_cutoff values basis set dependent? As in, if I swap the basis sets I ran my testing on based on advice I get here, would I need to rerun my testing for the new basis set, or could I use the values I determined using the original basis sets?<br></li><li>Are there any 'rules of thumb' as to how big a cell should be used for a system like mine? Essentially, I've embedded my molecule in a crystal of argon. I'm <i>not</i> trying to extrapolate this to a larger crystal structure (I think that means it's 'non periodic'?)</li></ol><div>I'd really appreciate any help I can get with this; I'm very much out of my depth here...</div><div><br></div><div>Thanks in advance!</div><div><br></div><div>Matthew~</div></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/8285cd89-82bb-4950-a344-f1d2dfbf8191n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/8285cd89-82bb-4950-a344-f1d2dfbf8191n%40googlegroups.com</a>.<br />