Dear Devyesh Rana,<div><br></div><div>Setting the target  of the constraint imposes that the collective variable assumes that value immediately. </div><div>Setting the target growth, one can determine a rate of change of the target value along the simulation. </div><div>This method is not optimal to estimate free energy profiles reliably. I would rather suggest the thermodynamics integration. </div><div>Anyway CP2K just runs the MD and if required can impose restraints or constraints. How you extract free energy surfaces from the data of the MD simulations is something you need to learn from the literature.  There are books and review articles that can help you. </div><div>Kind regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, November 2, 2020 at 3:30:19 PM UTC+1 Dev Rana 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;"><div dir="ltr">Hi Marcella,<div><br></div><div>Thank you for your suggestions. Is there an example of the 3rd method anywhere? I <b>think </b>this is the technique that I actually need. I have some bond lengths (RDF) from ab initio MD between different species so I'd like to place two atoms within the solvent. Subsequently, constraint 1 atom, and push the other atom toward a target distance (based on RDF). However, I'm not sure how to do this correctly. The closest example I have found is <a href="https://www.cp2k.org/exercises:2020_uzh_acpc2:ex03?s[]=target" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://www.cp2k.org/exercises:2020_uzh_acpc2:ex03?s%5B%5D%3Dtarget&source=gmail&ust=1604424648210000&usg=AFQjCNEczzgIMp_2JEwR8jCrIq7tUrc3Og">https://www.cp2k.org/exercises:2020_uzh_acpc2:ex03?s[]=target</a> and <a href="https://www.cp2k.org/exercises:2017_ethz_mmm:nacl_free_energy?s[]=target" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://www.cp2k.org/exercises:2017_ethz_mmm:nacl_free_energy?s%5B%5D%3Dtarget&source=gmail&ust=1604424648210000&usg=AFQjCNF-1KyPQP2mqsg1YEun2uq7gF1FiQ">https://www.cp2k.org/exercises:2017_ethz_mmm:nacl_free_energy?s[]=target</a>. However, they don't explain what CP2K is doing when you specify a target. The first example seems to do what I think is the proper technique in which initial coordinates differ from the TARGET setpoint. I'll play around with this example. The 2nd example explicitly sets the distance of the starting atom COORD the same as the TARGET, i.e. TARGET = ATOM COORD. So even though this 2nd example is MD, the atom is at the target distance from the beginning of the simulation. </div><div><br></div><div>So my question is for the 1st example, does this method, for example, allow you to place two atoms 10A away and set TARGET to 1.72A? Will MD then forcibly through some sort of steepest descent push these two atoms towards 1.72A? Is it this simple?</div><div><br></div><div>Thank you again!</div><div><br></div><div>Best Regards,</div><div>Devyesh Rana</div></div><div dir="ltr"><div><br>On Monday, November 2, 2020 at 3:23:24 AM UTC-5, Marcella Iannuzzi wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div><br></div>Dear Devyesh Rana<div><br></div><div>You can run a series of MD with restraints, i.e., following an Umbrella Sampling procedure.</div><div>You can run a series of MD with constraints, sampling the corresponding Lagrangian multiplier, i.e., following a Thermodynamic Integration procedure. </div><div>You can run MD with constraints, specifying  a target, a target_growth, and a target_limit to steer the change along the MD.</div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto">On Friday, October 30, 2020 at 9:06:04 PM UTC+1 Dev Rana 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">Is there an easy way of preconfiguring a system for steered molecular dynamics?<div><br></div><div>Specifically, I'd like to understand a reaction in a solvent via AIMD. However, instead of waiting for spontaneous reactions to occur, I'd like to forcibly make particles interact. Using the potential of mean force to obtain transition states etc. Is there an easy way to constrain atoms and automatically have them incrementally move towards each other?</div><div><br></div><div>I think the nudged elastic band method might help, but this may not be the best avenue. </div><div><br></div><div>Any suggestions?</div><div><br></div><div>Best Regards,</div><div>Devyesh Rana</div></blockquote></div></blockquote></div></div></blockquote></div>