Dear CP2K Developers,<div><br /></div><div>I am running Gibbs Ensemble Monte Carlo example from</div><div>https://www.cp2k.org/howto:gemc<br /></div><div><br /></div><div>I have a question of why the PRESSURE in the MC section needs to be specified for the GEMC_NVT ensemble. </div><div><br /></div><div>I have a single component H2O system, and as long as temperature is defined in the Monte Carlo simulations, all the other intensive variables, such as  pressure should come out as the result of the simulation, not the input parameter.</div><div><br /></div><div>Could you please explain the meaning of the PRESSURE in this case?</div><div>If I do not specify it, I get the error as shown below.</div><div><div>The MC block for one of the boxes is also shown below.</div><div><br /></div></div><div>Thank you,</div><div>Slava</div><div>Vyacheslav Bryantsev</div><div><br /></div><div>*******************************************************************************<br /> *   ___                                                                       *<br /> *  /   \                                                                      *<br /> * [ABORT]                                                                     *<br /> *  \___/      Value requested, but no value set getting value from keyword    *<br /> *    |                           PRESSURE of section MC                       *<br /> *  O/|                                                                        *<br /> * /| |                                                                        *<br /> * / \                                        input/input_section_types.F:1097 *<br /> *******************************************************************************<br /></div><div><br /></div><div>The MC block for one of the boxes is shown below.</div><div><br /></div><div>&MOTION<br />  &MC<br />    ENSEMBLE GEMC_NVT<br />    TEMPERATURE 398.0<br />    IPRINT 1<br />    LBIAS no<br />    LSTOP yes<br />    NMOVES 8<br />    NSWAPMOVES 640<br />    NSTEP 500<br />    PRESSURE 1.013<br />    RESTART no<br />    BOX2_FILE_NAME GEMC_NVT_box2.inp<br />    RESTART_FILE_NAME mc_restart_1<br />    ETA 0.0<br />!    VIRIAL_TEMPS 300.0<br />    &MOVE_PROBABILITIES<br />      PMSWAP 0.20<br />      PMTRAION 0.30<br />      PMTRANS 0.60<br />      PMVOLUME 0.05<br />      PMHMC 0.0<br />      PMAVBMC 0.0<br />      &MOL_PROBABILITIES<br />        PMAVBMC_MOL 1.0<br />        PMSWAP_MOL 1.0<br />        PMTRAION_MOL 1.0<br />        PMTRANS_MOL 1.0<br />        PMROT_MOL 1.0<br />      &END MOL_PROBABILITIES<br />      &BOX_PROBABILITIES<br />        PMVOL_BOX 1.0<br />        PMHMC_BOX 1.0<br />      &END BOX_PROBABILITIES<br />    &END MOVE_PROBABILITIES<br />    &MOVE_UPDATES<br />      IUPTRANS 100<br />      IUPVOLUME 100<br />    &END MOVE_UPDATES<br />    &MAX_DISPLACEMENTS<br />      &MOL_DISPLACEMENTS<br />        RMDIHEDRAL 3.0<br />        RMANGLE 3.0<br />        RMBOND 0.074<br />        RMROT 26.0<br />        RMTRANS 0.38<br />      &END MOL_DISPLACEMENTS<br />      &BOX_DISPLACEMENTS<br />        RMVOLUME 100.5<br />      &END BOX_DISPLACEMENTS<br />    &END MAX_DISPLACEMENTS<br />    &AVBMC<br />      AVBMC_ATOM 1<br />      AVBMC_RMIN 1.0<br />      AVBMC_RMAX 5.0<br />      PBIAS 0.5<br />    &END AVBMC<br />  &END MC<br />&END MOTION<br /></div><div><br /></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+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/f70d82a8-3c28-406c-90c2-71b4ba58d700n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/f70d82a8-3c28-406c-90c2-71b4ba58d700n%40googlegroups.com</a>.<br />