<div dir="ltr">Dear Frederick,<div><div>Thank you for your early reply and detailed explanation. My calculations run smoothly under this setting.</div><div>Sincerely,</div><div>Tianhua Wang</div><div><br></div><br>在 2019年10月16日星期三 UTC+8下午4:18:51，Frederick Stein写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Tianhua,<div><br></div><div>Either use FRACTION 0.54 with SCALE_COULOMB 1.0 (default value for SCALE_COULOMB) or FRACTION 1.0 with SCALE_COULOMB 0.54. Both should work. The second way is more convenient when you mix interaction potentials. Basically you can use every combination of SCALE_COULOMB and FRACTION as long as the product of their input values is 0.54.</div><div><br></div><div>Best,</div><div>Frederick<br><br>Am Mittwoch, 16. Oktober 2019 03:44:37 UTC+2 schrieb Tianhua Wang:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Frederick,<div>I also want to use the M06-2X functional in CP2K through the first option you provided above, because of the version (5.1) I'm using, and the &XC_FUNCTIONAL section in the input file is written as</div><div>&XC_FUNCTIONAL</div><div><div>  &LIBXC</div><div>    FUNCTIONAL XC_MGGA_X_M06_2X XC_MGGA_C_M06_2X</div><div>  &END LIBXC</div><div>&END XC_FUNCTIONAL</div><div><br></div><div>However, I'm hesitating at the keyword FRACTION in the &HF section. The M06-2X functional has 54% HF, and I wonder if I should write "FRACTION 0.54" in this section as</div><div>&HF</div><div><div>  FRACTION 0.54</div><div>  &SCREENING</div><div>    EPS_SCHWARZ 1.0E-6</div><div>    SCREEN_ON_INITIAL_P TRUE</div><div>  &END SCREENING</div><div>  &INTERACTION_POTENTIAL</div><div>    POTENTIAL_TYPE TRUNCATED</div><div>    CUTOFF_RADIUS 6.0</div><div>    T_C_G_DATA t_c_g.dat</div><div>  &END INTERACTION_POTENTIAL</div><div>  &MEMORY</div><div>    MAX_MEMORY 5000</div><div>    EPS_STORAGE_SCALING 0.1</div><div>  &END MEMORY</div><div>&END HF</div></div><div><br></div><div>Thank you for your attention to my consultation and looking forward to your reply.</div><div>Sincerely,</div><div>Tianhua Wang</div><div><br></div><br>在 2019年9月3日星期二 UTC+8下午8:15:53，Frederick Stein写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Subhasish,<div><br></div><div>The DFT part of M06-2X should be XC_HYB_MGGA_XC_M06_2X. It is basically the combination of XC_MGGA_X_M06_2X and XC_MGGA_C_M06_2X. When you want to use M06-2X within CP2K, you have the following options:</div><div>1) Within 6.1.0 use XC_MGGA_X_M06_2X and XC_MGGA_C_M06_2X.</div><div>2) Within 6.1.0 adjust the code by either using the patches in <a href="https://groups.google.com/forum/#!searchin/cp2k/wb97x%7Csort:date/cp2k/ARya5uM-pA4/2wt5PgtkAQAJ" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/forum/#!searchin/cp2k/wb97x%7Csort:date/cp2k/ARya5uM-pA4/2wt5PgtkAQAJ';return true;" onclick="this.href='https://groups.google.com/forum/#!searchin/cp2k/wb97x%7Csort:date/cp2k/ARya5uM-pA4/2wt5PgtkAQAJ';return true;">https://groups.google.com/<wbr>forum/#!searchin/cp2k/wb97x%<wbr>7Csort:date/cp2k/ARya5uM-pA4/<wbr>2wt5PgtkAQAJ</a> by Stanislav or applying the patches from PRs 304 and 422 (<a href="https://github.com/cp2k/cp2k/pull/304" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgithub.com%2Fcp2k%2Fcp2k%2Fpull%2F304\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEen-bT38tekm6G-oAqtDqjxcaTvQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgithub.com%2Fcp2k%2Fcp2k%2Fpull%2F304\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEen-bT38tekm6G-oAqtDqjxcaTvQ';return true;">https://github.com/cp2k/cp2k/<wbr>pull/304</a>, <a href="https://github.com/cp2k/cp2k/pull/422" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgithub.com%2Fcp2k%2Fcp2k%2Fpull%2F422\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHm-vZwSy9xb5671ebqM30zi2sUcA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgithub.com%2Fcp2k%2Fcp2k%2Fpull%2F422\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHm-vZwSy9xb5671ebqM30zi2sUcA';return true;">https://github.com/<wbr>cp2k/cp2k/pull/422</a>) and use XC_HYB_MGGA_XC_M06_2X.</div><div>3) Use the latest trunk from CP2K with XC_HYB_MGGA_XC_M06_2X. (The PRs from 2 have already been included.)</div><div>Please keep in mind that you need an additional HF section in all of these three cases.</div><div><br></div><div>Best regards,</div><div>Frederick<br><br>Am Dienstag, 3. September 2019 12:59:46 UTC+2 schrieb SUBHASISH MALLICK:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>I want to use M06-2X DFT functional in conjunction with 6-31G* basis sets in cp2k version 6.1.0. In this purpose I used libxc library, where, only XC_MGGA_C_M06_2X functional is supported in cp2k version 6.1.0 and the  XC_MGGA_X_M06_2X and <br></div><div>XC_HYB_MGGA_XC_M06_2X are not supported in the current version of cp2k. May I ask you among the above three M06-2X functionals mentioned in libxc which is equivallent to the gaussian M06-2X and why other two functionals are not supported in this version? </div><div>Is there any way to implement the HYB_MGGA functionals in cp2k version 6.1.0.</div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div>