<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>