<div dir="ltr">Dear Vladimir,<br>thank you for your reply. I applied BSSE correction and for my system I got BSSE error around 4kJ/mol.  The Eads I got from CP2K is around -88 kJ/mol while the VASP data that I have available predicts Eads for a water = -74 kJ/mol. As I do not have much experience with RPA, I use the VASP data as a reference to "tune" my CP2K setup. Please, do you have some suggestion what can influence the accuracy (or to be the source of the discrepancy)? As an input I use PBE+d3 (I tried also PBE only) wavefunction optimized with the same basis set and I've already tested the following criteria:<br>EPS_DEFAULT;<br>EPS_SCHWARZ;<br>RPA_NUM_QUAD_POINTS;<br>SCREEN_ON_INITIAL_P;<br>CUTOFF RADIUS for the Coulomb potential (I use max. value, i.e. 1/2 of the unit cell)<br><div><br></div><div>Attached are my input files for both PBE-D3 pre-optimization as well for RPA.</div><div><br></div>Any help is really apreciated!<br>Thank you<br>Kind regards,<br>Katarina<br><br><br>Dňa utorok, 6. augusta 2019 14:51:39 UTC+2 Vladimir Rybkin napísal(-a):<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Katarina,<div><br></div><div>cc-TZ are correlation-consistent triple-zera quality basis sets. Normally, they are OK, however, sometimes one needs to go for the complete basis set extrapolation, e.g. for cohesion energies. In data/<span style="font-family:monospace"><span style="color:rgb(0,0,0)">BASIS_RI_cc-TZ there are cc-TZ basis sets for O, H, Si and Al. So, I would try these first. You may also (and actually should) do BSSE correction. Bear in mind that 4- and 5-zeta basis sets will make the computation much more expensive, so it may happen that you will not be able to run those calculations. In any case, if you really need the 4- and 5-zeta basis, please write again.</span></span></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)"><br></span></span></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)">Yours,</span></span></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)"><br></span></span></div><div><span style="font-family:monospace"><span style="color:rgb(0,0,0)">Vladimir </span><br></span><br>понедельник, 29 июля 2019 г., 17:54:59 UTC+2 пользователь katarína stančiaková написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear all,</div><div>I would like to perform RPA and MP2 calculations on a periodic system of cca 40 atoms to get reference adsorption energy for small molecules in alomunosillicates. In my system I have H, O, Al and Si atoms. Which basis set would you recommend to use? I tried cc-TZ in combination with RI-TZ (with GTH_POTENTIALS), however, I was wondering whether there is something "better". I found a tutorial how to make better RI basis set, however, I could not find bigger primary basis set. I don't have much experience with these techniques, so any help would be very much appreciated!</div><div><br></div><div>Thank you</div><div>Kind regards,</div><div>Katarina<br></div></div></blockquote></div></div></blockquote></div>