<div dir="ltr"><div>Hi,<br>please ignore it directly.<br></div>Zhendong<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 12, 2016 at 11:08 AM,  <span dir="ltr"><<a href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<br>
<br>
this warning is related to a problem described in detail in<br>
<br>
M. Guidon et al J. Chem. Theory Comput. 2009, 5, 3010–3021<br>
<br>
Your options are:<br>
<br>
1) ignore it and hope your numerics is OK, very well possible as<br>
   you are using a small ADMM basis<br>
<br>
2) go to a larger simulation box<br>
<br>
3) Switch to the truncated PBE0 functional as described in the paper.<br>
<br>
Also have a look at the different integral screening options in<br>
order to make the calculation more stable.<br>
<br>
regards<br>
<br>
Juerg<br>
<br>
--------------------------------------------------------------<br>
Juerg Hutter                         Phone : <a href="tel:%2B%2B41%2044%20635%204491" value="+41446354491">++41 44 635 4491</a><br>
Institut für Chemie C                FAX   : <a href="tel:%2B%2B41%2044%20635%206838" value="+41446356838">++41 44 635 6838</a><br>
Universität Zürich                   E-mail: <a href="mailto:hut...@chem.uzh.ch">hut...@chem.uzh.ch</a><br>
Winterthurerstrasse 190<br>
CH-8057 Zürich, Switzerland<br>
---------------------------------------------------------------<br>
<br>
-----<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a> wrote: -----To: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a><br>
From: Apostolidou Christina<br>
Sent by: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a><br>
Date: 02/12/2016 10:03AM<br>
Subject: [CP2K:7424] I need help for HSE-functional calculation<br>
<div class="HOEnZb"><div class="h5"><br>
Dear cp2k developers,<br>
<br>
Please help me.<br>
<br>
I am trying to perform an AIMD-calculation with the HSE-functional.<br>
<br>
I get the following warning:<br>
<br>
 *** 15:27:23 WARNING in hfx_types:hfx_create_neighbor_cells :: Periodic ***<br>
 *** Hartree Fock calculation requested with use of a truncated or       ***<br>
 *** shortrange potential. The cutoff radius is larger than half the     ***<br>
 *** minimal cell dimension. This may lead to unphysical total energies. ***<br>
 *** Reduce the cutoff radius in order to avoid possible problems.       ***<br>
 *** hfx_types.F line 1448  <br>
<br>
This is my input:<br>
<br>
&GLOBAL<br>
 PROJECT auxpothse_funco19<br>
 RUN_TYPE MD<br>
 PRINT_LEVEL LOW<br>
 FFTW_PLAN_TYPE PATIENT<br>
&END GLOBAL<br>
<br>
&FORCE_EVAL<br>
 &DFT<br>
  BASIS_SET_FILE_NAME BASIS_ALL<br>
  POTENTIAL_FILE_NAME POTENTIAL<br>
 UKS<br>
 multiplicity 2<br>
  &MGRID<br>
   CUTOFF 280<br>
   REL_CUTOFF 40<br>
   NGRIDS 5<br>
  &END<br>
 &QS<br>
   EPS_FILTER_MATRIX 1.0E-12<br>
   EPS_PGF_ORB 1.0E-12<br>
  &END<br>
 &AUXILIARY_DENSITY_MATRIX_METHOD<br>
   METHOD BASIS_PROJECTION<br>
   ADMM_PURIFICATION_METHOD MO_DIAG<br>
  &END<br>
  &SCF<br>
   SCF_GUESS ATOMIC<br>
#  SCF_GUESS RESTART<br>
   MAX_SCF 1000<br>
   EPS_SCF 1.0E-5<br>
   &OT<br>
    MINIMIZER CG<br>
    PRECONDITIONER FULL_ALL<br>
 ENERGY_GAP 0.002<br>
   &END<br>
   &OUTER_SCF<br>
    EPS_SCF 1.0E-5<br>
    MAX_SCF 20<br>
   &END<br>
  &END<br>
#  &LOCALIZE<br>
#   METHOD CRAZY<br>
#   MAX_ITER 2000<br>
#   USE_HISTORY<br>
#   &PRINT<br>
#    &WANNIER_CENTERS<br>
#     &EACH<br>
#      MD 8<br>
#     &END<br>
#     IONS+CENTERS<br>
#     FILENAME =<a href="http://ohradikal-wannier.xyz" rel="noreferrer" target="_blank">ohradikal-wannier.xyz</a><br>
#    &END<br>
#   &END<br>
#  &END<br>
 &XC<br>
 &XC_FUNCTIONAL<br>
  &XWPBE<br>
     SCALE_X -0.25<br>
     SCALE_X0 1.0<br>
     OMEGA 0.2<br>
   &END<br>
   &PBE<br>
      SCALE_X 0.0<br>
      SCALE_C 1.0<br>
   &END PBE<br>
 &END XC_FUNCTIONAL<br>
 &XC_GRID<br>
    XC_DERIV NN10_SMOOTH<br>
    XC_SMOOTH_RHO NN10<br>
   &END<br>
  &HF<br>
   &INTERACTION_POTENTIAL<br>
    POTENTIAL_TYPE SHORTRANGE<br>
    OMEGA 0.19<br>
    &END<br>
       &SCREENING<br>
          EPS_SCHWARZ 1.0E-10<br>
        &END<br>
       FRACTION 0.25<br>
      &END<br>
   &VDW_POTENTIAL<br>
    DISPERSION_FUNCTIONAL PAIR_POTENTIAL<br>
    &PAIR_POTENTIAL<br>
     TYPE DFTD3<br>
     PARAMETER_FILE_NAME dftd3.dat<br>
     REFERENCE_FUNCTIONAL PBE<br>
    &END<br>
   &END<br>
  &END<br>
#  &PRINT<br>
#   &E_DENSITY_CUBE<br>
#    &EACH<br>
#     MD 8<br>
#    &END<br>
#    FILENAME =Phenol-density.cube<br>
#    APPEND<br>
#    STRIDE 2 2 2<br>
#   &END<br>
#  &END<br>
 &END<br>
<br>
 &SUBSYS<br>
  &CELL<br>
   ABC 9.988226 9.988226 9.988226<br>
  &END<br>
  &COORD             <br>
 &END<br>
  &KIND H<br>
   BASIS_SET DZVP-MOLOPT-SR-GTH<br>
  AUX_FIT_BASIS_SET cFIT3<br>
   POTENTIAL GTH-PBE-q1<br>
  &END<br>
#  &KIND C<br>
#   BASIS_SET DZVP-MOLOPT-SR-GTH<br>
#   POTENTIAL GTH-B3LYP-q4<br>
#  &END<br>
  &KIND O<br>
   BASIS_SET DZVP-MOLOPT-SR-GTH<br>
   AUX_FIT_BASIS_SET cFIT3<br>
   POTENTIAL GTH-PBE-q6<br>
  &END<br>
 &END<br>
&END<br>
<br>
&MOTION<br>
 &MD<br>
  ENSEMBLE NVT<br>
  STEPS 10000<br>
  TIMESTEP 0.5<br>
  &THERMOSTAT<br>
   TYPE NOSE<br>
   REGION MASSIVE<br>
   &NOSE<br>
    TIMECON 10<br>
   &END<br>
  &END<br>
  TEMPERATURE 370.15<br>
 &END<br>
 &PRINT<br>
  &RESTART<br>
   &EACH<br>
    MD 1<br>
   &END<br>
  &END<br>
 &END<br>
&END<br>
<br>
#&EXT_RESTART<br>
# RESTART_FILE_NAME -1.restart<br>
# RESTART_THERMOSTAT F<br>
#&END<br>
<br>
<br>
I<br>
 tried also to vary the cutoff-radii in the SE-Section: There are the<br>
following subsections in which you can define the cutoff-radius:<br>
COULOMB, EXCHANGE, SCREENING, LR_CORRECTION. (I tried it also with very small cutoff-radii)<br>
<br>
Nothing helps!<br>
Does anyone has another idea or knows this problem?<br>
Maybe I forgot one keyword which defines the cutoff-radius?<br>
<br>
I<br>
 need to take the HSE06 functional because I need the screened exchange<br>
in combination with the MIC in order to calculate my system as a<br>
condensed phase.<br>
<br>
Thank you!<br>
<br>
Best regards<br>
<br>
Christina Apostolidou<br>
<br>
<br>
<br>
--<br>
<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="noreferrer" target="_blank">https://groups.google.com/group/cp2k</a>.<br>
<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="noreferrer" target="_blank">https://groups.google.com/d/optout</a>.<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="noreferrer" target="_blank">https://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="noreferrer" target="_blank">https://groups.google.com/d/optout</a>.<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr">Zhendong Guo <div><div>Phone: +41786943316</div><div>Email:  <a href="mailto:zhendon...@gmail.com" target="_blank">zhendon...@gmail.com</a></div><div>            <a href="mailto:zhendo...@epfl.ch" target="_blank">zhendo...@epfl.ch</a>-------------------------------------------------------------------------------------------------</div></div></div></div>
</div>