<div dir="ltr">Dear all, <div><br></div><div><br></div><div>I have a question regarding CP2K MD simulations using PLUMED. This might become lengthy to explain, but I hope that someone can perhaps give some advice.</div><div>I will make corresponding files available through google drive.</div><div><br></div><div>The aim is to simulate the rearrangement of a substrate coordinated to a metal centre. Specifically, this involves a "flip" of the coordination site of a double bond. </div><div>From previous static calculations I have found out a possible pathway for this, which serves as a guidance for setting up the MD simulations. I have defined 3 structures (start, transitions state, final) in the file <b>path_flip.pdb</b><b>. </b>This should give you an idea what process I am interested in, essentially the substrate flips around. As a note aside, here I am interested in a molecular system, i.e. these simulations will use a cubic supercell, and perhaps this will introduce some problems (see below). Later on, I wish to return to the actual system, and embed this molecular portion back into it's solid state environment (for various reasons it is good to look at the molecular system first).</div><div><br></div><div>Now, since this process can be viewed as a conformational change, an approach that appears reasonable here is the sampling method proposed by Parinello (J. Chem. Phys., 2007, 126, 054103), which uses as the collective variable an expression that measures the root-mean-square-deviation (RMSD) of the sampled structures along a defined reaction path (as defined in <b>path_flip.pdb</b>). This method appears to work quite nicely with PLUMED. Indeed, a quick MD simulation (500 steps) using the movingconstraints keyword in PLUMED and forcing the trajectory along the defined path shows that the transition is smooth and follows the defined path (see https://drive.google.com/open?id=0B3k7Eu2qqsBBNUl2dFlneWROajA). </div><div><br></div><div>The reference point is always provided in the <b>path_flip.pdb</b> file. Note also that only the atoms of the substrate are taken into account, in order to remove "noise" from the movement of the other atoms (this is specified in the .pdb file). Having established that this is a suitable collective variable describing the rearrangement under study, I then moved on to sample the trajectory using metadynamics (see https://drive.google.com/open?id=0B3k7Eu2qqsBBQ21hczVpbUt1azQ).</div><div><br></div><div>I have done 5 MD runs so far (4 restarts), and this is where things get confusing. I have uploaded the files corresponding to the last restart. The file <b>movie.xyz </b>file contains the complete trajectory, which I have concatenated from the individual simulations. In each case I would follow the following protocol:</div><div><br></div><div>(i) copy the CP2K .restart file from the previous md run to new folder</div><div>(ii) copy path_flip.pdb, plumed.dat, HILLS and COLVAR files to this folder as well</div><div>(iii) resubmit job.</div><div><br></div><div><b>plumed.dat</b> contains the RESTART keyword, and I hope this is sufficient information to continue with the simulation. The PLUMED.OUT file reports that 192 Gaussians were read, which implies that the restart has worked. Previous runs restarted from fewer numbers of hills, naturally. </div><div><br></div><div>What are the problems:</div><div><br></div><div>(i) From the COLVAR file it appears that the collective variable as defined above (=path) is changing and the substrate is undergoing some conformational changes. In the start the value is around 1.0, fluctuates around this value for a while, then changes rapidly at around t=13-14 ps, appears to reach the transition phase at 20 ps and is a t the final structure thereafter (although there is some strong fluctuations going on there). However, none of this is reflected in the actual trajectory, where the substrate moves around but never undergoes the transition. There appears to be a dichotomy between the values of the CV and the actual structure. </div><div><br></div><div>(ii) It also appears that there are jumps in the sampling, as seen from the COLVAR and HILL file. For example after 8 ps the next hills are only spawned at 9.44 ps. How does this happen? Are my settings not right?</div><div><br></div><div>(iii) Another point of concern relates to the fact that this simulation is molecular, and uses a cubic box (25 Ang). You can see that the molecule rotates a lot, and I wonder if this is affects the estimation of the RMSD value (i.e. the CV). I believe that this method can deal with translations, but not so sure about rotations. HAs anyone experience with this, when the molecule moves a lot during the simulation. I am also using the center coordinates keyword (WAVELET!), how does this affect the results? Is it a problem that the reference structures are in fact sitting on one corner of the box, while the molecule for simulations is at the center (I think not, since this should be accounted for through transnational symmetry). Anyways, I would appreciate any reassurance here. </div><div><br></div><div>It appears that this is a neat way of doing these simulations, but it seems to be also quite tricky to get it right. I am fairly new to this kind of simulations, so any help would be very much appreciated. The initial simulations hold promise, now I am trying to get the actual mtd right to get some reliable results. Ultimately I would like to sample the system in the solid state, i.e. using periodicity, in order to model this rearrangement in dependence of the crystal structure. It would be nice to work out the reason for variations in the barriers for this conformational change in different systems, which we also study experimentally. </div><div><br></div><div>I am also using umbrella sampling as an alternative to the mtd part. I am also aware that other collective variables could be defined to sample this process (even without resorting to PLUMED), but I am simply interested in making this work. </div><div><br></div><div>Thanks for your help, and apologies for this lengthy post. I hope that there are some experts out there that can help.</div><div><br></div><div><br></div><div>Tobias</div></div>