<div dir="ltr">Dear cp2k user,<div><br><div>          I am attemping to run NPT ensemble to simulate 9-CH4 molecules with <font color="#ff0000">BLYP</font> under the condition of 8atm and 500K, but I failed. I have already read some previous discussions on this forum, most of them are unrelated to my work and some experienced users advise to monitor the stress tensor along the NVT (yet, I don't konw how to do it ). I will very appreciate it if you can give me possible ways, thanks in advance. Below is my input file:</div><div>   </div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  STRESS_TENSOR ANALYTICAL</div><div>  &DFT</div><div>  BASIS_SET_FILE_NAME /share/apps/CP2K/data/BASIS_MOLOPT</div><div>  POTENTIAL_FILE_NAME /share/apps/CP2K/data/POTENTIAL</div><div>  CHARGE 0</div><div>  MULTIPLICITY 1</div><div>  &SCF</div><div>  SCF_GUESS ATOMIC</div><div>  EPS_SCF 1.0E-6</div><div>  MAX_SCF 50</div><div>  &OUTER_SCF</div><div>  MAX_SCF 10</div><div>  &END OUTER_SCF</div><div>  &OT</div><div>  # My scheme</div><div>  PRECONDITIONER FULL_ALL</div><div>  MINIMIZER DIIS</div><div>  &END OT</div><div>  &PRINT</div><div>  &RESTART</div><div>  &EACH</div><div>  MD 20</div><div>  &END EACH</div><div>  &END RESTART</div><div>  &RESTART_HISTORY OFF</div><div>  &END RESTART_HISTORY</div><div>  &END PRINT</div><div>  &END SCF</div><div>  &QS</div><div>      METHOD GPW</div><div>      EPS_DEFAULT 1.0E-11</div><div>    &END QS</div><div>  &MGRID</div><div>  COMMENSURATE</div><div>  CUTOFF 350</div><div>  &END MGRID</div><div>  &XC</div><div>  &XC_FUNCTIONAL BLYP</div><div>  &END XC_FUNCTIONAL</div><div>   &XC_GRID</div><div>       XC_DERIV SPLINE2</div><div>       XC_SMOOTH_RHO NN10</div><div>     &END XC_GRID</div><div>  &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>  &CELL</div><div>    ABC 42.472  42.472 42.472</div><div>    &END CELL</div><div>    &COORD</div><div> C     23.213   2.485  28.226</div><div> H     22.154   2.308  27.944</div><div> H     23.264   3.313  28.963</div><div> H     23.634   1.563  28.677</div><div> H     23.797   2.757  27.323</div><div> C     24.261  26.373  14.721</div><div> H     24.718  27.114  15.409</div><div> H     25.052  25.918  14.091</div><div> H     23.757  25.582  15.313</div><div> H     23.517  26.880  14.073</div><div> C      2.114   8.706  29.027</div><div> H      2.061   9.413  29.880</div><div> H      3.022   8.918  28.426</div><div> H      2.160   7.666  29.411</div><div> H      1.213   8.825  28.391</div><div> C     23.010  13.312   0.946</div><div> H     23.509  13.888   0.139</div><div> H     22.505  12.426   0.510</div><div> H     23.767  12.980   1.685</div><div> H     22.259  13.956   1.448</div><div> C     21.366   1.604  31.128</div><div> H     22.443   1.688  31.383</div><div> H     21.231   0.825  30.350</div><div> H     20.792   1.324  32.036</div><div> H     21.001   2.578  30.744</div><div> C      4.466   0.666  39.061</div><div> H      4.402  -0.206  39.743</div><div> H      4.457   0.314  38.009</div><div> H      3.598   1.336  39.234</div><div> H      5.406   1.221  39.256</div><div> C     38.649  19.173  34.531</div><div> H     37.957  18.359  34.831</div><div> H     38.651  19.962  35.310</div><div> H     39.673  18.763  34.417</div><div> H     38.314  19.606  33.566</div><div> C     22.852  30.055  14.682</div><div> H     22.298  29.280  14.115</div><div> H     23.816  29.637  15.036</div><div> H     23.046  30.927  14.024</div><div> H     22.246  30.377  15.555</div><div> C     21.823  27.509  17.108</div><div> H     21.052  28.051  17.694</div><div> H     22.789  28.050  17.180</div><div> H     21.943  26.483  17.512</div><div> H     21.508  27.454  16.046</div><div>&end COORD</div><div>   &KIND H</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q1</div><div>      POTENTIAL GTH-BLYP-q1</div><div>    &END KIND</div><div>    &KIND C</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q4</div><div>      POTENTIAL GTH-BLYP-q4</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&GLOBAL</div><div>  PROJECT ch4-blyp-8atm-9-5</div><div>  RUN_TYPE MD</div><div>  PRINT_LEVEL MEDIUM</div><div>  EXTENDED_FFT_LENGTHS</div><div>  &TIMINGS</div><div>     THRESHOLD 0.000001</div><div>  &END</div><div>&END GLOBAL</div><div>&MOTION</div><div>  &MD</div><div>    ENSEMBLE NPT_F</div><div>    STEPS 15</div><div>    TIMESTEP 0.5</div><div>    TEMPERATURE 500</div><div>&BAROSTAT</div><div>   PRESSURE 8.0</div><div>   TIMECON [wavenumber_t] 1000</div><div>   &THERMOSTAT</div><div>    TYPE CSVR</div><div>   &END THERMOSTAT</div><div>  &END BAROSTAT</div><div>  &THERMOSTAT</div><div>    TYPE NOSE</div><div>    REGION GLOBAL</div><div>   &NOSE</div><div>    LENGTH                 3</div><div>    YOSHIDA                3</div><div>    TIMECON                1000</div><div>    MTS                         2</div><div>   &END NOSE</div><div>  &END THERMOSTAT</div><div>  &END MD</div><div>&END MOTION</div><div><br></div></div></div>