Dear all,<div><br /><div>I'm trying to carry out umbrella sampling of a dihedral angle (torsion) using the MOTION / CONSTRAINT / COLLECTIVE / RESTRAINT section in CP2K. I came across a problem with the unit of the force constant K. I understand that users can specify energy units like kJ/mol to the constant, but the actual unit of the constant should contain the unit of the torsion angle as well, say kJ/(mol-rad^2) or kJ/(mol-deg^2). I am wondering what unit is used for the torsion part in the unit of K. Does it use the same user-specified unit for TARGET when running MD simulations? Or does it use the default angle unit for the constant?</div><div><br /></div><div>As an example, below is my input for the CONSTRAINT section:</div><div><br /></div><div>  &CONSTRAINT<br />    &COLLECTIVE<br />      COLVAR 1<br />      TARGET [deg] 10<br />      INTERMOLECULAR .TRUE.<br />      &RESTRAINT<br />         K [kjmol] 100.0<br />      &END RESTRAINT<br />    &END COLLECTIVE<br />  &END CONSTRAINT<br /></div><div><br /></div><div>from which my intention was to use a 100 kJ/(mol-deg^2) force constant for my bias. Yet what I found from my restart file (generated automatically by CP2K using its default unit) is:</div></div><div><br /></div><div>   &CONSTRAINT<br />     &COLLECTIVE<br />       COLVAR  1<br />       INTERMOLECULAR  T<br />       TARGET     1.7453292519943295E-01<br />       &RESTRAINT<br />         K     3.8087988796022196E-02<br />       &END RESTRAINT<br />     &END COLLECTIVE<br />     &COLVAR_RESTART<br />           1.7453292519943295E-01<br />     &END COLVAR_RESTART<br />   &END CONSTRAINT<br /></div><div><br /></div><div>in which the force constant seems to be 100 kJ/(mol-rad^2), which is 0.03 kJ/(mol-deg^2). </div><div>(K = 3.808E-02 hartree = 100kJ/mol, and TARGET = 1.745E-01 rad = 10 degree)</div><div><br /></div><div>I would really appreciate any help in understanding the units here. Thank you!</div><div><br /></div><div>With regards,</div><div>Yuqing Fu</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/59270e08-d26f-47ff-9650-064c226c9e84n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/59270e08-d26f-47ff-9650-064c226c9e84n%40googlegroups.com</a>.<br />