<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body dir="auto">
Hi Tommaso
<div><br>
</div>
<div>It seems that you forgot to comment the &END of the commented &THERMOSTAT section above the &ENERGY section which prematurely ends your &MD section.</div>
<div><br>
</div>
<div>Best regards <br>
<br>
<div id="AppleMailSignature">Matthias</div>
<div><br>
Am 07.11.2017 um 19:07 schrieb "<a href="mailto:tfran...@gmail.com">tfran...@gmail.com</a>" <<a href="mailto:tfran...@gmail.com">tfran...@gmail.com</a>>:<br>
<br>
</div>
<blockquote type="cite">
<div>
<div dir="ltr">Dear all,
<div>I'd like to ask help for properly setting the annealing calculation of my molecular crystal. The ultimate goal would be description of the cell dynamic from a high temperature monoclinic phase to low temperature triclinic one. I would like to test both
 heating and cooling procedures, with the intent to identify the collective variables to later use in a metadynamic calculation. </div>
<div><br>
</div>
<div>So, as I set the input (heating mode) following the few examples reported, the calculation crashes saying that in the &PRINT section the &ENERGY block is not recognized. Can you kindly help me fix the problem or tell me I should avoid to put in the input
 file? </div>
<div><br>
</div>
<div>Here the input:</div>
<div><br>
</div>
<div>
<div> &GLOBAL</div>
<div>   PRINT_LEVEL  MEDIUM</div>
<div>   PROJECT_NAME ANN</div>
<div>   RUN_TYPE  MD</div>
<div>   WALLTIME 23:00:00</div>
<div> &END GLOBAL</div>
<div><br>
</div>
<div> &MOTION</div>
<div>    &MD</div>
<div>    ENSEMBLE NVT</div>
<div>    STEPS  2000</div>
<div>    TIMESTEP 1.00</div>
<div>    TEMPERATURE 180.0</div>
<div>    COMVEL_TOL 1.0E-8</div>
<div>   <b> ANNEALING_CELL 0.5</b></div>
<div>#    &THERMOSTAT</div>
<div>#      TYPE CSVR</div>
<div>#      &CSVR</div>
<div>#        TIMECON 50</div>
<div>#      &END CSVR</div>
<div>    &END</div>
<div>    &PRINT</div>
<div>      <font color="#ff0000">&ENERGY</font></div>
<div><font color="#ff0000">        FILENAME =md.ener</font></div>
<div><font color="#ff0000">        &EACH</font></div>
<div><font color="#ff0000">          MD 1</font></div>
<div><font color="#ff0000">        &END EACH</font></div>
<div><font color="#ff0000">      &END ENERGY</font></div>
<div>      &PROGRAM_RUN_INFO</div>
<div>        &EACH</div>
<div>          MD 1</div>
<div>        &END EACH</div>
<div>      &END PROGRAM_RUN_INFO</div>
<div>    &END PRINT</div>
<div>  &END MD</div>
<div>  &PRINT</div>
<div>    &TRAJECTORY</div>
<div>      FILENAME =<a href="http://md.xyz">md.xyz</a></div>
<div>        &EACH</div>
<div>          MD 1</div>
<div>        &END EACH</div>
<div>    &END TRAJECTORY</div>
<div>    &RESTART</div>
<div>      FILENAME =md.restart</div>
<div>      BACKUP_COPIES 0</div>
<div>     &EACH</div>
<div>        MD 10</div>
<div>      &END EACH</div>
<div>    &END RESTART</div>
<div>    &RESTART_HISTORY</div>
<div>      &EACH</div>
<div>        MD 10</div>
<div>      &END EACH</div>
<div>    &END RESTART_HISTORY</div>
<div>  &END PRINT</div>
<div> &END MOTION</div>
<div><br>
</div>
<div> &FORCE_EVAL</div>
<div>   METHOD  QS</div>
<div>   &DFT</div>
<div>     BASIS_SET_FILE_NAME BASIS_SET</div>
<div>     POTENTIAL_FILE_NAME POTENTIAL</div>
<div>     UKS  T</div>
<div>     MULTIPLICITY  1</div>
<div>     &SCF</div>
<div>       MAX_SCF  500</div>
<div>       EPS_SCF     9.9999999999999995E-07</div>
<div>       SCF_GUESS  ATOMIC</div>
<div>       &OT  T</div>
<div>         ALGORITHM  IRAC</div>
<div>         MINIMIZER  DIIS</div>
<div>         PRECONDITIONER  FULL_KINETIC</div>
<div>       &END OT</div>
<div>       &OUTER_SCF  T</div>
<div>         EPS_SCF     9.9999999999999995E-07</div>
<div>         MAX_SCF  10</div>
<div>       &END OUTER_SCF</div>
<div>     &END SCF</div>
<div>     &QS</div>
<div>       EPS_DEFAULT     1.0000000000000000E-10</div>
<div>       EXTRAPOLATION  ASPC</div>
<div>       EXTRAPOLATION_ORDER  4</div>
<div>     &END QS</div>
<div>     &MGRID</div>
<div>       NGRIDS  5</div>
<div>       CUTOFF     6.0000000000000000E+02</div>
<div>       REL_CUTOFF     6.0000000000000000E+01</div>
<div>     &END MGRID</div>
<div>     &XC</div>
<div>       DENSITY_CUTOFF     1.0000000000000000E-10</div>
<div>       GRADIENT_CUTOFF     1.0000000000000000E-10</div>
<div>       TAU_CUTOFF     1.0000000000000000E-10</div>
<div>       &XC_FUNCTIONAL  NO_SHORTCUT</div>
<div>         &PBE  T</div>
<div>         &END PBE</div>
<div>       &END XC_FUNCTIONAL</div>
<div>       &VDW_POTENTIAL</div>
<div>         POTENTIAL_TYPE  PAIR_POTENTIAL</div>
<div>         &PAIR_POTENTIAL</div>
<div>           TYPE  DFTD3(BJ)</div>
<div>           PARAMETER_FILE_NAME ./dftd3.dat</div>
<div>           REFERENCE_FUNCTIONAL PBE</div>
<div>           CALCULATE_C9_TERM  T</div>
<div>           REFERENCE_C9_TERM  T</div>
<div>           VERBOSE_OUTPUT  T</div>
<div>         &END PAIR_POTENTIAL</div>
<div>       &END VDW_POTENTIAL</div>
<div>     &END XC</div>
<div>   &END DFT</div>
<div>   &SUBSYS</div>
<div>     &CELL</div>
<div>       A     2.0731000000000002E+01    0.0000000000000000E+00    0.0000000000000000E+00</div>
<div>       B     1.6600000000000000E-03    1.4442399999999999E+01    0.0000000000000000E+00</div>
<div>       C    -4.6196099999999998E+00   -6.6000000000000000E-04    1.4826380000000004E+01</div>
<div>       PERIODIC  XYZ</div>
<div>       MULTIPLE_UNIT_CELL  1 1 1</div>
<div>     &END CELL</div>
<div>     &COORD</div>
<div> S    7.622173    0.787062    2.602696</div>
<div> S   10.024760   16.523216   11.155081</div>
<div> S    6.768842   -0.348410    9.481586</div>
<div> S   10.878093   17.658689    4.276192</div>
<div> S    9.76797 ....</div>
</div>
<div><br>
</div>
<div>....</div>
<div>
<div>     &KIND N2</div>
<div>       BASIS_SET TZV2P-MOLOPT-GTH</div>
<div>       ELEMENT N</div>
<div>       POTENTIAL GTH-PBE-q5</div>
<div>       &POTENTIAL</div>
<div>2 3</div>
<div>0.2837905100000000E+00 2 -0.1241522559000000E+02  0.1868095920000000E+01</div>
<div>2</div>
<div>0.2554050000000000E+00 1  0.1363026257000000E+02</div>
<div>0.2454945300000000E+00 0</div>
<div>       &END POTENTIAL</div>
<div>       &BS  T</div>
<div>         &ALPHA</div>
<div>           NEL  -2</div>
<div>           L  1</div>
<div>           N  2</div>
<div>         &END ALPHA</div>
<div>         &BETA</div>
<div>           NEL  2</div>
<div>           L  1</div>
<div>           N  2</div>
<div>         &END BETA</div>
<div>       &END BS</div>
<div>     &END KIND</div>
<div>     &TOPOLOGY</div>
<div>       NUMBER_OF_ATOMS  256</div>
<div>       MULTIPLE_UNIT_CELL  1 1 1</div>
<div>     &END TOPOLOGY</div>
<div>   &END SUBSYS</div>
<div> &END FORCE_EVAL</div>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>Thanks in advance for the help!!!</div>
<div>All the Best,</div>
<div>Tommaso</div>
</div>
<p></p>
-- <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+unsu...@googlegroups.com">cp2k+unsu...@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">https://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout">https://groups.google.com/d/optout</a>.<br>
</div>
</blockquote>
</div>
</body>
</html>