<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;"><br>
Dear Gromacs Developers,<br>
I would like to allow Gromacs to do Markovian Milestoning because it seems a very<br>
effective algorithm for the estimation of both free energy and kinetics. To achieve<br>
this goal I have to confine the simulation inside a Voronoi cell.<br>
To do this confinment, in litterature it is suggested (see Venturoli and Vanden Eijden "Markovian Milestoning") that one can, upon Voronoi cell violation, invert velocities and reset the position to the previous integration step such that the system can remain
inside the Voronoi cell.<br>
<br>
I decided to interface to Gromacs via plumed2 such that collective variables can be<br>
used to compute the metric that rules the Voronoi cell.<br>
<br>
Following plumed2 philosophy I got velocities just before the do_force call in md.c.<br>
Assuming a velocity verlet integrator which are the velocities that I get in this point of the code?<br>
I mean, which time step?<br>
<br>
Do you suggest a specific strategy to invert velocities and set the previous positions?<br>
<br>
Gromacs expects an atom in a domain decomposition cell; after the reset of the position to the previous step, this cell may have changed thus leading to an inconsitency. I got this kind of problems.<br>
<br>
How can I set the positions to the previous integration step and at the same time assuring
<br>
that Gromacs is aware of this change of positions from the domain decomposition stand point?<br>
<br>
This problem can be avoided by using particle decomposition, but it is probably suboptimal.<br>
<br>
I have done many trials and I have other issues. <br>
Often lincs algorithm fails; when I compute<br>
the 'reflection' on the Voronoi facet I use subset of atoms because the metric is defined in a<br>
subset of atoms. To get the maximal stability should I change velocities and reset positions<br>
of the whole system? Still lincs should be enabled or disabled at all?<br>
<br>
Thanks,<br>
Cheers,<br>
Sergio Decherchi<br>
</div>
</body>
</html>