<div dir="ltr">Dear Katarina,<div><br></div><div>I believe you are doing a good job in how you treat the exchange and number of quadrature points. A few points:</div><div>1) I would not use D3 together with RPA, as it is double counting: RPA is responsible for intermolecular interactions. I see that you don't actually do it, which is good.</div><div>2) You may actually reduce CUTOFF      900  in the WFC_CORRELATION section to speed-up you calculation and without much effect on the accuracy.</div><div>3) If 40 quadrature points the largest value? With 50-60 it should converge.</div><div>4) SCREEN_ON_INITIAL_P affects efficiency, not the result.</div><div>5) EPS_DEFAULT are sensitive parameters EPS_SCHWARZ (but you already checked them). They should be sufficiently small.</div><div>6) What you have done is actually RPA with exchange <i style="box-sizing: border-box; outline: none; color: rgb(0, 0, 0); font-family: Roboto, arial, sans-serif; font-size: 17px;">E</i><span class="stack" style="box-sizing: border-box; outline: none; font-size: 12.75px; line-height: 0; position: relative; vertical-align: baseline; bottom: -0.25em; color: rgb(0, 0, 0); font-family: Roboto, arial, sans-serif;">tot</span><span class="stack" style="box-sizing: border-box; outline: none; font-size: 12.75px; line-height: 0; position: relative; vertical-align: baseline; top: -0.5em; color: rgb(0, 0, 0); font-family: Roboto, arial, sans-serif;">EXX/RPA</span> (thus HF section in WFC section). This is an accurate method, but I am not sure whether you computed the same in VASP. </div><div>7) You .xyz file contains only water molecule, so I can't say much about the system. You may try to increase the z dimension of the cell to make more empty space between you periodic images.</div><div>8) Did VASP employ k-point sampling? CP2K is Gamma-point only.</div><div>9) Yes, it may still be basis set incompleteness responsible for the difference. But there are options to check (see above).</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir</div><div><br>пятница, 9 августа 2019 г., 12:03:19 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">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></blockquote></div></div>