<div dir="ltr"><div>Dear Conrad,</div><div><br></div><div>I appreciate your quick reply. But I still have some uncertainty.</div><div><br></div><div>The paper [1] introduced the GGA_XC_OPBE_D, refitting procedure applied to PBE to optimize the (meta-)GGA functional. The prefix ‘o’ stands for optimized.</div><div>While the OPBE functional is Handy-Cohen’s 2001 OPTX correction for exchange plus Perdew-Burke-Ernzerhof’s 1996 correction for correlation. [2]</div><div>OPBE is already built into the Gaussian[3] and ADF[4], which shows some accuracy in capturing the features of metal complex. One paper [5] reports that employing OPBE calculate the Ru complex in CP2K, whereas the OPBE haven’t built in CP2K, I am really want to know the setting about OPBE is right or not. </div><div><br></div><div><div><br></div><div>[1] Goerigk, Lars, and Stefan Grimme. “A General Database for Main Group Thermochemistry, Kinetics, and Noncovalent Interactions − Assessment of Common and Reparameterized (meta-)GGA Density Functionals.” Journal of Chemical Theory and Computation 6, no. 1 (January 12, 2010): 107–26. doi:10.1021/ct900489g.</div><div>[2] Liao, Meng-Sheng, John D. Watts, and Ming-Ju Huang. “Electronic Structure of Some Substituted Iron(II) Porphyrins. Are They Intermediate or High Spin?” The Journal of Physical Chemistry A 111, no. 26 (July 1, 2007): 5927–35. doi:10.1021/jp070734z. </div><div>[3] Andrikopoulos, Prokopis C., Carine Michel, Sandra Chouzier, and Philippe Sautet. “In Silico Screening of Iron-Oxo Catalysts for CH Bond Cleavage.” ACS Catalysis 5, no. 4 (April 3, 2015): 2490–99. doi:10.1021/cs500996k.</div><div>[4] Swart, Marcel, Andreas W. Ehlers, and * Koop Lammertsma. “Performance of the OPBE Exchange-Correlation Functional.” Molecular Physics 102, no. 23–24 (December 10, 2004): 2467–74. doi:10.1080/0026897042000275017.</div><div>[5] Vallés-Pardo, José Luis, Marieke C. Guijt, Marcella Iannuzzi, Khurram S. Joya, Huub J. M. de Groot, and Francesco Buda. “Ab Initio Molecular Dynamics Study of Water Oxidation Reaction Pathways in Mono-Ru Catalysts.” ChemPhysChem 13, no. 1 (January 16, 2012): 140–46. doi:10.1002/cphc.201100546.</div><div><br></div></div><div><br></div>On Monday, July 18, 2016 at 1:35:50 PM UTC-4, Conrad wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Jianmin,<div><br></div><div><br></div><div>When you want to use an exotic XC functional, it is better to use the LIBXC library.</div><div>All available functionals are shown here: </div><div><a href="http://www.tddft.org/programs/octopus/wiki/index.php/Libxc_3.0_functionals" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.tddft.org%2Fprograms%2Foctopus%2Fwiki%2Findex.php%2FLibxc_3.0_functionals\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEPerN-b5ltOM-sv8g30cam1sF3Gw';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.tddft.org%2Fprograms%2Foctopus%2Fwiki%2Findex.php%2FLibxc_3.0_functionals\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEPerN-b5ltOM-sv8g30cam1sF3Gw';return true;">http://www.tddft.org/programs/<wbr>octopus/wiki/index.php/Libxc_<wbr>3.0_functionals</a><br></div><div>(Assuming you built CP2K with LIBXC version 3.0. Older versions are elsewhere on the website.)</div><div><br></div><div>In your input you include the functional like this:</div><div><div><br></div><div>            &XC_FUNCTIONAL</div><div>                &LIBXC</div><div>                    FUNCTIONAL GGA_XC_OPBE_D</div><div>                &END LIBXC</div><div>            &END XC_FUNCTIONAL</div></div><div><br></div><div>In this case, as it is a GGA functional there is nothing else to be done. If it was a hybrid functional, for example, you would need also to make a HF section. </div><div><br></div><div>For your pair potential section, set REFERENCE FUNCTIONAL to OPBE. </div><div><br></div><div>I'm not familiar with this functional, but this should be sufficient. Try this and see if it works for you.</div><div>It runs on my machine. Do some benchmarking before doing any production runs, of course. </div><div><br></div><div>Best wishes,</div><div>Conrad</div><div><br></div><div><br></div></div></blockquote></div>