<div dir="ltr">Hi Geoffrey,<br><br>The libxc v2.1 and v2.2 were added to the SVN version of cp2k (https://groups.google.com/forum/#!searchin/cp2k/Marco/cp2k/kyPFdrTV4uc/-skBy_pEMYUJ). The xc_libxc.F 
was updated. Using the latest libxc versions may help. In my previous post regarding meta-GGA functionals, I may have made an incorrect statement. It may not be possible to perform cp2k calculations with meta-GGAs which require 2nd derivatives but single-point and geometry optimization calculations should run. Can you compile CP2K-2.5.1 with libxc version 2.1.0 and test your calculations again. The XC_MGGA_X_M06 and XC_MGGA_C_M06 components should be combined.<br><br>Best regards,<br>Marco<br><br>On Monday, July 21, 2014 9:16:39 AM UTC-4, Geoffrey Wood wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div> </div><div>I was wondering if this problem is resolvable or will be addressed in version 2.6?  I can run Hybrid-GGA functionals (eg B3LYP and PBE0) using the libxc library but hybrid-meta-GGAs all fail at the SCF (i.e not a geometry optimization just a single-point). I've compiled CP2k version 2.5.1 svn:13632 and libxc version 2.01.  I've tried to run M06 and BR89 both of with stop at the SCF initiation. However the TPSS functional, which doesn't have an HF section as given in the CP2K tests does run (i.e  H2O-tpssx_libxc.inp). As a side question the libxc library contains both exchange and correlation parts of these functionals, i.e. XC_MGGA_X_M06 and XC_MGGA_C_M06 are these suppossd to be combined somehow?</div><div><br>Thanks</div><div> </div><div> </div><div><br>On Sunday, April 13, 2014 2:33:23 PM UTC-4, Jingyun Ye wrote:</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><div dir="ltr">Hi, all<div><br></div><div>I am trying to use M06-2X functional in CP2K (PBE potential)  to compare the binding energy of CO2  with Gaussian results( also use M06-2X functional, same basis sets). </div><div><br></div><div>But the calculation stops  at the SCF every time with the following information. But if I change the functional M06-2X to PBE, the job runs well. </div><div><br></div><div>Anyone know what's the problem? Any suggestion will be help. </div><div>Thanks very much.</div><div><br></div><div> </div><div><div> SCF WAVEFUNCTION OPTIMIZATION</div><div><br></div><div>  ------------------------------<wbr>----- OT ------------------------------<wbr>---------</div><div><br></div><div>  Allowing for rotations:  F</div><div>  Optimizing orbital energies:  F</div><div>  Minimizer      : CG                  : conjugate gradient</div><div>  Preconditioner : FULL_ALL            : diagonalization, state selective</div><div>  Precond_solver : DEFAULT</div><div>  Line search    : 2PNT                : 2 energies, one gradient</div><div>  stepsize       :    0.15000000</div><div>  energy_gap     :    0.00100000</div><div><br></div><div>  eps_taylor     :   0.10000E-15</div><div>  max_taylor     :             4</div><div><br></div><div>  mixed_precision    : F</div><div><br></div><div>  ------------------------------<wbr>----- OT ------------------------------<wbr>---------</div><div><br></div><div>  Step     Update method      Time    Convergence         Total energy    Change</div><div>  ------------------------------<wbr>------------------------------<wbr>------------------</div><div><br></div><div>  Trace(PS):                                  209.9999999999</div><div>  Electronic density on regular grids:                   NaN                 NaN</div><div>  Core density on regular grids:              210.0000000000        0.0000000000</div><div>  Total charge density on r-space grids:                 NaN</div><div>  Total charge density g-space grids:                    NaN</div></div><div><br></div></div></blockquote></div></blockquote></div>