Hi Hongxia,<div><br></div><div>I guess because you have set DeltaT to 50000K... I'm not surprised the atoms fly apart in this case. </div><div><br></div><div>If this is your first metadynamics run, I'd recommend running an unbiased simulation first (i.e. not WELL_TEMPERED and no DELTA_T set) and see how the CV evolves. You can then set the width of the Gaussians (SCALE) according to variations in the value of the CV with time, for example by taking the standard deviation of the CV and setting this to be SCALE. Once you have set an appropriately sized gaussian, you can then introduce the bias. The bias factor is given by (T+deltaT)/T. Given 298K is your starting temperature and you were using 50000 as deltaT, your bias factor is rather large. A good initial range might be somewhere in the region of 5 or 10 for the bias factor, so deltaT around 1000 - 2000 K. Of course this is system dependent so you need to try several values, the higher the value the more of free-energy space you will explore, though you want to make sure it's comparable to the free energy barriers you want to traverse to avoid getting situations where you explode the simulation as you have done by using 50000K for deltaT.</div><div><br></div><div>Best wishes,</div><div>Matt<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, 20 October 2021 at 00:24:02 UTC+1 hongxi...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear cp2k users and developers,<div><br></div><div>I am trying to do a metadynamics in cp2k. I set 3 collective variables. I want to use CV1 and CV2 as collective variables for metadynamics, while at the same time using CV3 for a restraint to constraint the motion of two molecules not too far away from each other. Could anyone help me set this up appropriately? </div><div><br></div><div>My input is as folllows.</div><div>I set up the three COLVAR in SUBSYS section first. <div><div>&COLVAR 1</div><div> &COORDINATION</div><div> ATOMS_FROM 2 3</div><div> ATOMS_TO 1 4</div><div> R_0 [angstrom] 1.2</div><div> NN 8</div><div> ND 14</div><div> &END COORDINATION</div><div> &END COLVAR</div><div><br></div><div> &COLVAR 2</div><div> &COORDINATION</div><div> ATOMS_FROM 5 6</div><div> ATOMS_TO 1 4</div><div> R_0 [angstrom] 1.2</div><div> NN 8</div><div> ND 14</div><div> &END COORDINATION</div><div> &END COLVAR</div><div><br></div><div> &COLVAR 3</div><div> &DISTANCE</div><div> POINTS 1 2</div><div> &POINT 1</div><div> ATOMS 2 3</div><div> TYPE GEO_CENTER</div><div> &END POINT</div><div><br></div><div> &POINT 2</div><div> ATOMS 5 6</div><div> TYPE GEO_CENTER</div><div> &END POINT</div><div> &END DISTANCE</div><div> &END COLVAR</div></div><div><br></div></div><div>And then I add the constraint section in the MOTION section.</div><div><div>&MOTION</div><div><br></div><div> &CONSTRAINT</div><div> &COLLECTIVE</div><div> COLVAR 3</div><div> INTERMOLECULAR .TRUE.</div><div> TARGET [angstrom] 5.0</div><div> &RESTRAINT</div><div> K 100000</div><div> &END RESTRAINT</div><div> &END COLLECTIVE</div><div> &END CONSTRAINT</div><div><br></div><div> &MD</div><div> ENSEMBLE NVT</div><div> STEPS 400000</div><div> TIMESTEP 0.5</div><div> TEMPERATURE 298.0</div><div> &THERMOSTAT</div><div> TYPE CSVR</div><div> REGION MASSIVE</div><div> &CSVR</div><div> TIMECON [fs] 200.0</div><div> &END</div><div> &END</div><div> &END MD</div><div> &FREE_ENERGY</div><div> &METADYN</div><div><br></div><div> DO_HILLS</div><div> WELL_TEMPERED</div><div> DELTA_T 50000</div><div> NT_HILLS 30</div><div> WW 0.005</div><div><br></div><div> &METAVAR</div><div> SCALE 0.2</div><div> COLVAR 1</div><div> &END METAVAR</div><div><br></div><div> &METAVAR</div><div> SCALE 0.2</div><div> COLVAR 2</div><div> &END METAVAR</div><div><br></div><div> &PRINT</div><div> &COLVAR</div><div> COMMON_ITERATION_LEVELS 2</div><div> &EACH</div><div> MD 1</div><div> &END</div><div> &END</div><div> &HILLS</div><div> COMMON_ITERATION_LEVELS 2</div><div> &EACH</div><div> MD 1</div><div> &END</div><div> &END</div><div> &END</div><div> &END METADYN</div><div> &END FREE_ENERGY</div></div><div><br></div><div>Does anyone know what I did wrong and how to set it up appropriately instead? With this input, I got super huge potential energy and then the geometry fly away.</div><div><br></div><div>Thanks in advance!</div><div><br></div><div>Sincerely</div><div>Hongxia</div><div><br></div><div><br></div></blockquote></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/34f6871e-4d2c-4287-ad2b-e0be107841b1n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/34f6871e-4d2c-4287-ad2b-e0be107841b1n%40googlegroups.com</a>.<br />