<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="EN-US" link="blue" vlink="purple"> <div class="WordSection1"> BASIS_SET DZVP-GTH<o:p></o:p> POTENTIAL GTH-PADE<o:p></o:p> <o:p> </o:p> The PP GTH-PADE is sufficient to test the code, but obviously not suited to work with the M06-2X functional in a production run.<o:p></o:p> <o:p> </o:p> Best<o:p></o:p> <o:p> </o:p> Matthias<o:p></o:p> <o:p> </o:p> From: cp...@googlegroups.com <...@googlegroups.com> On Behalf Of mdsi...@gmail.com Sent: Freitag, 25. September 2020 15:27 To: cp2k <...@googlegroups.com> Subject: Re: [CP2K:13980] CP2K M06-2X water energy lower vs Other program<o:p></o:p> <o:p> </o:p> Matthias,<o:p></o:p> <div> <o:p> </o:p> </div> <div> What basis set and potential did you use for the M06-2X GPW calculations?<o:p></o:p> </div> <div> <o:p> </o:p> </div> <div> Thanks,<o:p></o:p> </div> <div> Stacey<o:p></o:p> </div> <div> <div> On Friday, September 25, 2020 at 9:00:05 AM UTC-4 Matthias Krack wrote:<o:p></o:p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <div> <div> Hi<o:p></o:p> <o:p></o:p> wrong numbers are obtained for mGGAs (like M06-2X) with GAPW, which is mandatory for all-electron calculations. The mGGA results, including M06-2X, using GPW are correct:<o:p></o:p> <o:p></o:p> DEBUG|======================== BEGIN OF SUMMARY ===============================<o:p></o:p> </div> </div> <div> <div> DEBUG| Atom Coordinate f(numerical) f(analytical) Difference Error [%]<o:p></o:p> </div> </div> <div> <div> DEBUG| 1 x -0.00000000 -0.00000000 -0.00000000 -<o:p></o:p> DEBUG| 1 y -0.00000000 -0.00000003 -0.00000003 -<o:p></o:p> DEBUG| 1 z -0.01487972 -0.01487506 0.00000465 0.03<o:p></o:p> DEBUG| 2 x 0.00000000 0.00000000 0.00000000 -<o:p></o:p> DEBUG| 2 y -0.00408127 -0.00408117 0.00000009 0.00<o:p></o:p> DEBUG| 2 z 0.00661315 0.00661293 -0.00000023 0.00<o:p></o:p> DEBUG| 3 x 0.00000000 0.00000000 0.00000000 -<o:p></o:p> DEBUG| 3 y 0.00408126 0.00408120 -0.00000006 0.00<o:p></o:p> DEBUG| 3 z 0.00661315 0.00661298 -0.00000017 0.00<o:p></o:p> DEBUG| Sum of differences: 0.00000524<o:p></o:p> DEBUG|======================== END OF SUMMARY =================================<o:p></o:p> <o:p></o:p> HTH<o:p></o:p> <o:p></o:p> Matthias<o:p></o:p> </div> </div> <div> <div> <o:p></o:p> From: c...@googlegroups.com <c...@googlegroups.com> On Behalf Of mdsi...@gmail.com Sent: Freitag, 25. September 2020 14:03 To: cp2k <c...@googlegroups.com> Subject: Re: [CP2K:13977] CP2K M06-2X water energy lower vs Other program<o:p></o:p> <o:p></o:p> I'm guessing I probably shouldn't use the M06-2X functional for anything right now or at least until why the energy discrepancy is figured out in the meta functionals?<o:p></o:p> <div> <div> On Thursday, September 24, 2020 at 1:13:45 PM UTC-4 Matthias Krack wrote:<o:p></o:p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <div> <div> Hi Jürg<o:p></o:p> <o:p></o:p> The debug run fails also with<o:p></o:p> <o:p></o:p> &XC_FUNCTIONAL TPSS<o:p></o:p> &END XC_FUNCTIONAL<o:p></o:p> </div> </div> <div> <div> <o:p></o:p> DEBUG| Atom Coordinate f(numerical) f(analytical) Difference Error [%]<o:p></o:p> </div> </div> <div> <div> DEBUG| 1 x -0.00000225 -0.00062077 -0.00061852 -99.64<o:p></o:p> DEBUG| 1 y -0.00000225 -0.00050477 -0.00050253 -99.55<o:p></o:p> DEBUG| 1 z 0.01173212 0.01219262 0.00046050 -3.78<o:p></o:p> </div> </div> <div> <div> <o:p></o:p> Matthias<o:p></o:p> <o:p></o:p> -----Original Message----- From: c...@googlegroups.com <c...@googlegroups.com> On Behalf Of h...@chem.uzh.ch Sent: Donnerstag, 24. September 2020 18:55 To: c...@googlegroups.com Subject: RE: [CP2K:13972] CP2K M06-2X water energy lower vs Other program<o:p></o:p> <o:p></o:p> Hi Matthias<o:p></o:p> <o:p></o:p> thank you. Can I ask you to run the same debug run for M06L and TPSS, two other meta functionals?<o:p></o:p> <o:p></o:p> Juerg<o:p></o:p> --------------------------------------------------------------<o:p></o:p> Juerg Hutter Phone : <a href="tel:+41%2044%20635%2044%2091" target="_blank"> ++41 44 635 4491</a><o:p></o:p> Institut für Chemie C FAX : <a href="tel:+41%2044%20635%2068%2038" target="_blank"> ++41 44 635 6838</a><o:p></o:p> Universität Zürich E-mail: h...@chem.uzh.ch<o:p></o:p> Winterthurerstrasse 190<o:p></o:p> CH-8057 Zürich, Switzerland<o:p></o:p> ---------------------------------------------------------------<o:p></o:p> <o:p></o:p> </div> </div> <div> <div> ----...@<a href="http://googlegroups.com" target="_blank">googlegroups.com</a> wrote: -----<o:p></o:p> </div> </div> <div> <div> To: "'c...@googlegroups.com'" <c...@googlegroups.com><o:p></o:p> From: "Krack Matthias (PSI)" <o:p></o:p> Sent by: c...@googlegroups.com<o:p></o:p> Date: 09/24/2020 06:49PM<o:p></o:p> Subject: RE: [CP2K:13971] CP2K M06-2X water energy lower vs Other program<o:p></o:p> <o:p></o:p> <o:p></o:p> Hi Jürg<o:p></o:p> <o:p></o:p> Unfortunately, there seems to be more than the 5d/6d difference when using the M06-2X functional. A comparison of the H2O structures relaxed with 6-31G** using BLYP, HF, B3LYP, and M06-2X reveals indeed a much larger difference for M06-2X:<o:p></o:p> <o:p></o:p> H2O/6-31G**<o:p></o:p> CCCBDB<o:p></o:p> CP2K<o:p></o:p> BLYP <o:p></o:p> -76.398885<o:p></o:p> -76.396840<o:p></o:p> HF <o:p></o:p> -76.023615<o:p></o:p> -76.023047<o:p></o:p> B3LYP <o:p></o:p> -76.419737<o:p></o:p> -76.418114<o:p></o:p> M06-2X <o:p></o:p> -76.383939<o:p></o:p> -76.378<o:p></o:p> <o:p></o:p> Moreover, the GEO_OPT run with M06-2X shows difference in the OH bonds which indicates wrong forces and indeed a DEBUG run reveals that there exists a problem related to M06-2X using the XC section as suggested:<o:p></o:p> <o:p></o:p> DEBUG| Atom Coordinate f(numerical) f(analytical) Difference Error [%]<o:p></o:p> DEBUG| 1 x 0.00000768 -0.00011877 -0.00012645 -106.47<o:p></o:p> DEBUG| 1 y -0.00006753 0.00169850 0.00176603 -103.98<o:p></o:p> DEBUG| 1 z -0.00096067 -0.00096000 0.00000066 0.07<o:p></o:p> application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0<o:p></o:p> <o:p></o:p> *******************************************************************************<o:p></o:p> * ___ *<o:p></o:p> * / \ *<o:p></o:p> * [ABORT] A mismatch between analytical and numerical forces has been *<o:p></o:p> * \___/ detected. Check the implementation of the analytical force *<o:p></o:p> * | calculation *<o:p></o:p> * O/| *<o:p></o:p> * /| | *<o:p></o:p> * / \ cp2k_debug.F:318 *<o:p></o:p> *******************************************************************************<o:p></o:p> <o:p></o:p> A debug run with B3LYP works fine with errors of only 0.01%.<o:p></o:p> <o:p></o:p> So there is either an input error (&XC section with M06-2X) or a bug in CP2K (or libxc).<o:p></o:p> <o:p></o:p> Best regards<o:p></o:p> <o:p></o:p> Matthias<o:p></o:p> <o:p></o:p> -----Original Message-----<o:p></o:p> From: c...@googlegroups.com <c...@googlegroups.com> On Behalf Of h...@chem.uzh.ch<o:p></o:p> Sent: Donnerstag, 24. September 2020 17:53<o:p></o:p> To: c...@googlegroups.com<o:p></o:p> Subject: Re: [CP2K:13969] CP2K M06-2X water energy lower vs Other program<o:p></o:p> <o:p></o:p> Hi<o:p></o:p> <o:p></o:p> so as I said, the basis sets are different and that makes the results different.<o:p></o:p> <o:p></o:p> JH<o:p></o:p> --------------------------------------------------------------<o:p></o:p> Juerg Hutter Phone : <a href="tel:+41%2044%20635%2044%2091" target="_blank"> ++41 44 635 4491</a><o:p></o:p> Institut für Chemie C FAX : <a href="tel:+41%2044%20635%2068%2038" target="_blank"> ++41 44 635 6838</a><o:p></o:p> Universität Zürich E-mail: h...@chem.uzh.ch<o:p></o:p> Winterthurerstrasse 190<o:p></o:p> CH-8057 Zürich, Switzerland<o:p></o:p> ---------------------------------------------------------------<o:p></o:p> <o:p></o:p> </div> </div> <div> <div> ----...@<a href="http://googlegroups.com" target="_blank">googlegroups.com</a> wrote: -----<o:p></o:p> </div> </div> <div> <div> To: "cp2k" <c...@googlegroups.com><o:p></o:p> From: "mdsi...@gmail.com" <o:p></o:p> Sent by: c...@googlegroups.com<o:p></o:p> Date: 09/24/2020 05:17PM<o:p></o:p> Subject: Re: [CP2K:13969] CP2K M06-2X water energy lower vs Other program<o:p></o:p> <o:p></o:p> Hi Prof Hutter,<o:p></o:p> <o:p></o:p> Here's from CP2K<o:p></o:p> Total number of - Atomic kinds: 2<o:p></o:p> - Atoms: 3<o:p></o:p> - Shell sets: 10<o:p></o:p> - Shells: 12<o:p></o:p> - Primitive Cartesian functions: 21<o:p></o:p> - Cartesian basis functions: 25<o:p></o:p> - Spherical basis functions: 24<o:p></o:p> <o:p></o:p> This is from Gaussian<o:p></o:p> Standard basis: 6-31G(d,p) (6D, 7F)<o:p></o:p> 25 basis functions, 42 primitive gaussians, 25 cartesian basis functions<o:p></o:p> On Thursday, September 24, 2020 at 10:27:08 AM UTC-4 jgh wrote:<o:p></o:p> Hi<o:p></o:p> check the number of basis functions. Part or hopefully all of the energy difference is from the (6d) default in Gaussian. <o:p></o:p> CP2K uses spherical functions only. <o:p></o:p> <o:p></o:p> regards <o:p></o:p> <o:p></o:p> JH<o:p></o:p> -------------------------------------------------------------- <o:p></o:p> Juerg Hutter Phone : <a href="tel:+41%2044%20635%2044%2091" target="_blank"> ++41 44 635 4491</a> <o:p></o:p> Institut für Chemie C FAX : <a href="tel:+41%2044%20635%2068%2038" target="_blank"> ++41 44 635 6838</a> <o:p></o:p> Universität Zürich E-mail: h...@chem.uzh.ch <o:p></o:p> Winterthurerstrasse 190<o:p></o:p> CH-8057 Zürich, Switzerland<o:p></o:p> --------------------------------------------------------------- <o:p></o:p> <o:p></o:p> -----c...@googlegroups.com wrote: -----<o:p></o:p> To: "cp2k" <c...@googlegroups.com><o:p></o:p> From: "mdsi...@gmail.com" <o:p></o:p> Sent by: c...@googlegroups.com<o:p></o:p> Date: 09/24/2020 04:00PM<o:p></o:p> Subject: [CP2K:13967] CP2K M06-2X water energy lower vs Other program <o:p></o:p> <o:p></o:p> Hello, <o:p></o:p> <o:p></o:p> The water energy at the M06-2X/6-31G** level is quite different when compared to Gaussian and the CCCBDB. <o:p></o:p> <o:p></o:p> Gaussian: -76.3831231205 Ha (#n M062X/6-31G** NoSymm SCF=(XQC, maxconventionalcycles=30) Int=UltraFine SP) <o:p></o:p> CCCBDB: -76.383939 Ha (<a href="https://cccbdb.nist.gov/energy3x.asp?method=58&basis=3&charge=0" target="_blank">https://cccbdb.nist.gov/energy3x.asp?method=58&basis=3&charge=0</a>) <o:p></o:p> CP2K: -76.401243254272231Ha <o:p></o:p> <o:p></o:p> Anyone have any ideas on what settings need to be modified or can CP2K match that energy? <o:p></o:p> <o:p></o:p> My CP2K input: <o:p></o:p> ****************************************************************************<o:p></o:p> @SET COORD_FILENAME <a href="http://water.xyz" target="_blank">water.xyz</a> <o:p> </o:p> <o:p></o:p> &GLOBAL <o:p></o:p> PROJECT H2O_Eng <o:p></o:p> RUN_TYPE ENERGY <o:p></o:p> PRINT_LEVEL LOW<o:p></o:p> &END GLOBAL <o:p></o:p> <o:p></o:p> &FORCE_EVAL <o:p></o:p> METHOD Quickstep <o:p></o:p> <o:p></o:p> &DFT <o:p></o:p> BASIS_SET_FILE_NAME EMSL_BASIS_SETS <o:p></o:p> POTENTIAL_FILE_NAME POTENTIAL <o:p></o:p> <o:p></o:p> CHARGE 0 <o:p></o:p> MULTIPLICITY 1 <o:p></o:p> <o:p></o:p> &QS <o:p></o:p> METHOD GAPW <o:p></o:p> EPS_DEFAULT 1.0E-9 <o:p></o:p> EXTRAPOLATION ASPC <o:p></o:p> MAP_CONSISTENT <o:p></o:p> EPSFIT 1.E-4 ! precision to give the extension of a hard gaussian <o:p></o:p> EPSISO 1.0E-12 <o:p></o:p> EPSRHO0 1.E-8 <o:p></o:p> LMAXN0 4 <o:p></o:p> LMAXN1 6 <o:p></o:p> ALPHA0_H 10 ! Exponent for hard compensation charge <o:p></o:p> &END QS <o:p></o:p> <o:p></o:p> &MGRID <o:p></o:p> NGRIDS 5 <o:p></o:p> CUTOFF 600 <o:p></o:p> REL_CUTOFF 100 <o:p></o:p> &END MGRID <o:p></o:p> <o:p></o:p> &SCF <o:p></o:p> MAX_SCF 50 <o:p></o:p> SCF_GUESS ATOMIC <o:p></o:p> EPS_SCF 1E-8 <o:p></o:p> &END SCF <o:p></o:p> <o:p></o:p> &POISSON <o:p></o:p> PERIODIC NONE <o:p></o:p> POISSON_SOLVER WAVELET <o:p></o:p> &END POISSON <o:p></o:p> <o:p></o:p> &XC <o:p></o:p> &XC_GRID <o:p></o:p> XC_DERIV NN10_SMOOTH <o:p></o:p> XC_SMOOTH_RHO NN10 <o:p></o:p> &END XC_GRID <o:p></o:p> &XC_FUNCTIONAL <o:p></o:p> &LIBXC <o:p></o:p> FUNCTIONAL MGGA_C_M06_2X <o:p></o:p> &END LIBXC <o:p></o:p> &LIBXC <o:p></o:p> FUNCTIONAL HYB_MGGA_X_M06_2X <o:p></o:p> &END LIBXC <o:p></o:p> &END XC_FUNCTIONAL <o:p></o:p> &HF <o:p></o:p> FRACTION 0.54 <o:p></o:p> &SCREENING <o:p></o:p> EPS_SCHWARZ 1.0E-10 <o:p></o:p> &END SCREENING <o:p></o:p> &MEMORY <o:p></o:p> MAX_MEMORY 1000 <o:p></o:p> EPS_STORAGE_SCALING 0.1 <o:p></o:p> &END MEMORY <o:p></o:p> &INTERACTION_POTENTIAL <o:p></o:p> POTENTIAL_TYPE COULOMB <o:p></o:p> &END INTERACTION_POTENTIAL <o:p></o:p> &END HF <o:p></o:p> &END XC <o:p></o:p> &END DFT <o:p></o:p> <o:p></o:p> &SUBSYS <o:p></o:p> &CELL <o:p></o:p> ABC 10.0 10.0 10.0 <o:p></o:p> ALPHA_BETA_GAMMA 90.000 90.000 90.000 <o:p></o:p> PERIODIC NONE <o:p></o:p> &END CELL <o:p></o:p> <o:p></o:p> &TOPOLOGY <o:p></o:p> COORD_FILE_FORMAT XYZ <o:p></o:p> COORD_FILE_NAME ${COORD_FILENAME} <o:p></o:p> &CENTER_COORDINATES <o:p></o:p> &END CENTER_COORDINATES <o:p></o:p> &END TOPOLOGY <o:p></o:p> <o:p></o:p> &KIND H <o:p></o:p> ELEMENT H <o:p></o:p> BASIS_SET 6-31Gxx <o:p></o:p> POTENTIAL ALL <o:p></o:p> LEBEDEV_GRID 80 <o:p></o:p> RADIAL_GRID 200 <o:p></o:p> &END KIND <o:p></o:p> <o:p></o:p> &KIND O <o:p></o:p> ELEMENT O <o:p></o:p> BASIS_SET 6-31Gxx <o:p></o:p> POTENTIAL ALL <o:p></o:p> LEBEDEV_GRID 80 <o:p></o:p> RADIAL_GRID 200 <o:p></o:p> &END KIND <o:p></o:p> &END SUBSYS<o:p></o:p> &END FORCE_EVAL <o:p></o:p> <o:p></o:p> <o:p></o:p> --<o:p></o:p> You received this message because you are subscribed to the Google Groups "cp2k" group. <o:p></o:p> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+...@googlegroups.com. <o:p></o:p> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/3d887e97-17e0-4340-835a-1e5db270f0e6n%40googlegroups.com" target="_blank"> https://groups.google.com/d/msgid/cp2k/3d887e97-17e0-4340-835a-1e5db270f0e6n%40googlegroups.com</a>. <o:p></o:p> <o:p></o:p> <o:p></o:p> --<o:p></o:p> You received this message because you are subscribed to the Google Groups "cp2k" group.<o:p></o:p> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+...@googlegroups.com.<o:p></o:p> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/30f65c0c-fe59-4345-bda0-7b142795e86dn%40googlegroups.com" target="_blank"> https://groups.google.com/d/msgid/cp2k/30f65c0c-fe59-4345-bda0-7b142795e86dn%40googlegroups.com</a>.<o:p></o:p> <o:p></o:p> --<o:p></o:p> You received this message because you are subscribed to the Google Groups "cp2k" group.<o:p></o:p> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+...@googlegroups.com.<o:p></o:p> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/OF6FE20CC5.0428E417-ONC12585ED.0057370C-C12585ED.0057370E%40lotus.uzh.ch" target="_blank"> https://groups.google.com/d/msgid/cp2k/OF6FE20CC5.0428E417-ONC12585ED.0057370C-C12585ED.0057370E%40lotus.uzh.ch</a>.<o:p></o:p> <o:p></o:p> --<o:p></o:p> You received this message because you are subscribed to the Google Groups "cp2k" group.<o:p></o:p> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+...@googlegroups.com.<o:p></o:p> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/b66e4f318bca4a2abc719d351e506c5d%40psi.ch" target="_blank"> https://groups.google.com/d/msgid/cp2k/b66e4f318bca4a2abc719d351e506c5d%40psi.ch</a>.<o:p></o:p> <o:p></o:p> -- <o:p></o:p> You received this message because you are subscribed to the Google Groups "cp2k" group.<o:p></o:p> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+...@googlegroups.com.<o:p></o:p> </div> </div> <div> <div> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/OFCAC7426D.F9467A62-ONC12585ED.005CF77A-C12585ED.005CF77D%40lotus.uzh.ch" target="_blank"> https://groups.google.com/d/msgid/cp2k/OFCAC7426D.F9467A62-ONC12585ED.005CF77A-C12585ED.005CF77D%40lotus.uzh.ch</a>.<o:p></o:p> </div> </div> </blockquote> </div> </div> </div> <div> <div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+...@googlegroups.com.<o:p></o:p> </div> </div> <div> <div> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/57263731-b6ee-45bb-bf00-412044a38c38n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/57263731-b6ee-45bb-bf00-412044a38c38n%40googlegroups.com</a>.<o:p></o:p> </div> </div> </blockquote> </div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsu...@googlegroups.com">cp2k+unsu...@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/db190420-3294-40e8-8abd-bb7d3a9efb6en%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/db190420-3294-40e8-8abd-bb7d3a9efb6en%40googlegroups.com</a>.<o:p></o:p> </div> </body> </html>